Hello,

 

I have a deriv function that I am feeding to nlme. It works, and I can use
it in nls, but when I try to use it in nlme I get Error: subscript out of
bounds.  I can fit the model using SSasympOrig, instead of the deriv
function, but I am trying to reproduce an earlier analysis (done in Splus)
and I get slightly different results with SSasympOrig.

 

Here are my calls:

 

This does fit the nlme model:

L1.lis<-nlsList(SSasympOrig,L1.gd)

L1.nlme<-nlme(L1.lis)                 

 

Here is the derive function -it works for the nls call:

boundedexp<-deriv(~100/(1+exp(-L))*(1-exp(-exp(logR)*x)), 

                        c("L", "logR"), function(x,L,logR){} )

 

L1.nls<-nls(ability ~ boundedexp(age,limit,lograte),

            data=L1.gd, start=c(limit=2,lograte=-3))

summary(L1.nls)

 

. but here I run into problems

L1.nlme2<-nlme(ability~boundedexp(age,limit,lograte),

            data=L1.gd, fixed=limit+lograte~1,start=c(limit=2,lograte=-3))

 

Any suggestions would be well appreciated!


	[[alternative HTML version deleted]]

