[R] Specify nonlinear equation with nlme
Dieter Menne
dieter.menne at menne-biomed.de
Sat Dec 20 20:35:23 CET 2008
justin bem <justin_bem <at> yahoo.fr> writes:
> I want to estimate a nonlinear mixed effect model using nlme package or lme4.
How can I specify this
>
> y=exp(B'X)/(1+exp(B'X)). the same model can be obtained with logit(y)=B'X but
I have many value of 1 in y logit(1)=Inf.
Using exp in this case is a fairly common approach, and there are some examples
in the Pinheiro/Bates Book, chapter 8.2 (look for orange tree); or, if you do
not have access to it, look at ch08.r in library/nlme/scripts
Dieter
More information about the R-help
mailing list