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

From: hong qin <alongway_at_gmail.com>
Date: Wed 25 Oct 2006 - 20:10:38 GMT


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]]



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu Oct 26 06:13:48 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 27 Oct 2006 - 02:30:12 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.