Re: [R] How to specify a constant in gnls{nlme}

2006-10-26 Thread Jost Burkardt
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

[R] How to specify a constant in gnls{nlme}

2006-10-25 Thread hong qin
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 ~