Hi,
It works for me if i encapsulate the formula within a function
a =1; b=1;
f1 - function(x,v,w) {a/(b+(x/v)^w)}
fm2 = gnls( y ~ f1(x,v,w), start=list( v=1, w=1))
Regards,
Jost
hong qin [EMAIL PROTECTED] writes:
Hi All,
I have question about speficifying a constant in gnls() from
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 ~