[R] need help on nlme()

From: Mingyu Feng <mingyufeng_at_gmail.com>
Date: Mon 20 Feb 2006 - 14:24:35 EST


Hello there,

I am using nlme() to fit a logistic mixed effect model on our data. The outcome variable is binary.
I got the error when I wanted to add a group factor to my model.

My initial model is as below:

model.a <- nlme(response~ 1/(1 + exp( -intercept- u0 - slope*TIME - u1*TIME)),

+ fixed=intercept+slope~1, random= u0+u1~1
|studentID,

+ start=c(slope=.01, intercept=-1), data=log.data,
method='ML')

This works fine on my data. But when i update it by adding a group factor SKILLS,
I got the error message:
"Error in nlme.formula(response ~ 1/(1 + exp( -intercept- u0 - slope*TIME :

        starting values for the fixed component are not the correct length"

The model is as below:
model.a <- nlme(response~ 1/(1 + exp( -intercept- u0 - slope*TIME - u1*TIME)),

+ fixed=intercept+slope ~ SKILLS, random= u0+u1~1
|studentID,

+ start=c(slope=.01, intercept=-1), data=log.data,
method='ML')

Does anybody see anything wrong with the "start" part of this model?

Thanks a lot!!

        [[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 Received on Mon Feb 20 14:29:24 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:36 EST