Hi All,

I have question about speficifying a constant in gnls() from package nlme.
Here is a testing code:

#############
library(nlme)
 x = exp( rnorm(100))
 y = 1/(1+x) + rnorm(100)/10
 plot( y ~ x)
 fm1 = gnls( y ~ 1/(1+(x/v)^w), start=list( v=1, w=1))

 a =1; b=1;
 fm2 = gnls( y ~ a/(b+(x/v)^w), start=list( v=1, w=1)) #This won't work
because I don't know to set $a$ and $b$ as constant.

############

Thanks,

HQ

	[[alternative HTML version deleted]]

