[R] LME without convergence

Ralf Sommerfeld sommerfeld at mpil-ploen.mpg.de
Wed Feb 28 12:57:04 CET 2007


Dear R-help list readers,

I am fitting a mixed model using the lme function (R V 2.3.1 for 
Windows). This is an example:

dep<-c(25,40,33.33,60,70.83,72,71.43,50,40,53.33,64,54.17,60,53.57)
yes<-c(0,1,2,3,4,5,6,0,1,2,3,4,5,6)
treat<-c(1,1,1,1,1,1,1,0,0,0,0,0,0,0) #factor

If I now fit a model with random slopes as well as intercepts:

model1<-lme(dep~yes,random=yes|treat)

R encounters the following problem:

Fehler in lme.formula(dep ~ yes, random = ~yes | treat) :
         iteration limit reached without convergence (9)

Following other instructions on that problem, I tried to adjust the 
iteration limits and tolerances using lmeControl(), however this did not 
work.

Can someone help me?

Best regards,
Ralf



More information about the R-help mailing list