Hi Dave, Bill,
Thank you. That was it!
Best Regards,
Steven
On Mar 25, 2015, at 1:11 PM, David R Forrest wrote:
> I'm away from a computer that can check syntax, but make sure the return is
> outside of the for loop. To my eye it looks inside.
[[alternative HTML version deleted]]
and the return() then exits the loop and routine with the unmodified result.
>
> Dave
>
> From: R-help [r-help-boun...@r-project.org] on behalf of Steven LeBlanc
> [ores...@gmail.com]
> Sent: Wednesday, March 25, 2015 3:43 PM
> To: r-help@
Greets,
I'm trying to iteratively find solutions to a problem given a range of options.
Code is as follows:
sim.app.wald<-function(B=0.1,min=10,max=50,alpha=0.05){
result<-c(fails=0,n=0)
for(n in min:max){
x<-seq(1,n-1,1)
fhat<-x/n
On Oct 20, 2013, at 9:54 PM, Mark Leeds wrote:
> Bill: I didn't look at the code but I think the OP means that during the
> nlminb algorithm,
> the variance covariance parameters hit values such that the covariance matrix
> estimate becomes negative definite.
Yes, that is what I meant.
>
>
Greets,
I'm trying to use nlminb() to estimate the parameters of a bivariate normal
sample and during one of the iterations it passes a parameter vector to the
likelihood function resulting in an invalid covariance matrix that causes
dmvnorm() to throw an error. Thus, it seems I need to somehow
On Oct 18, 2013, at 1:12 AM, peter dalgaard wrote:
>
> On Oct 18, 2013, at 08:37 , David Winsemius wrote:
>
>>
>> On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
>>
>>> Greets,
>>>
>>> I'm using nlminb() to estimate the paramete
On Oct 17, 2013, at 11:37 PM, David Winsemius wrote:
>
> On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
>
>> Greets,
>>
>> I'm using nlminb() to estimate the parameters of a multivariate normal
>> random sample with missing values and ran into an
Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal random
sample with missing values and ran into an unexpected result from my call to
dmvnorm() within the likelihood function. Particular details are provided
below. It appears that dmvnorm() makes a call to log(eigen
Greets,
I'm trying to learn to use nls and was running the example code for an
exponential model:
> x <- -(1:100)/10
> y <- 100 + 10 * exp(x / 2) + rnorm(x)/10
> nlmod <- nls(y ~ Const + A * exp(B * x))
Error in B * x : non-numeric argument to binary operator
In addition: Warning message:
In nl
Greets,
I'm using a data frame that looks like:
> head(pr2)
X1 X2 X3 X4Y fit res
1 44 33.2 5 30 41.2 39.22201 1.977991
2 43 33.8 4 41 31.7 38.48476 -6.784761
3 48 40.6 3 38 39.4 44.78278 -5.382783
4 52 39.2 7 48 57.5 51.48134 6.018656
5 71 45.5 11 53 74.8 68.25585 6.54415
Greets,
My data looks like:
> p3.18
s xbar subgroup
1 0.84 12.21
2 1.64 11.22
3 2.07 10.63
4 2.49 12.24
5 0.84 11.25
...
Using the command
> qcc(p3.18$xbar,type="xbar",sizes=5,center=mean(p3.18$xbar),std.dev=mean(p3.18$s)/0.94,title="X-bar
> C
11 matches
Mail list logo
