[R] nls arguments
chenu
chenu at ensam.inra.fr
Mon Dec 15 11:01:40 CET 2003
Hi all,
I've got a problem with the nls function.
I have an adjustment which works when I fix one of the argument of my
function (Xo=150) :
*Xo*=150
f<- function (tt*,Xo*,a,b) ifelse(tt<*Xo*,a*exp(-b**Xo*),a*exp(-b*tt))
ajust<-nls(RER~f(tt,*Xo*,a,b),data=data.frame(tt=Ph2[,2*k],RER=Ph2[,2*k+1]),start=list(a=0.5,b=0.014))
But, when I use it as a "normal" parameter (and even if I do not I use
it in the equation of my "f" function), it does not work :
f<- function (tt,*Xo*,a,b) ifelse(tt<150,a*exp(-b*tt),a*exp(-b*tt))
or : f<- function (tt,*Xo*,a,b)
ifelse(tt<*Xo*,a*exp(-b**Xo*),a*exp(-b*tt)) #it is this
equation I would like to resole
ajust<-nls(RER~f(tt,*Xo*,a,b),data=data.frame(tt=Ph2[,2*k],RER=Ph2[,2*k+1]),start=list(a=0.5,b=0.014,*Xo*=150))
/Error in nlsModel(formula, mf, start) : singular gradient matrix at
initial parameter estimates/
Would you have any idea to help me.
Thanks a lot
Karine Chenu
More information about the R-help
mailing list