[R] about lme

F Z gerifalte28 at hotmail.com
Fri Aug 6 21:44:15 CEST 2004


I suspect that what you want to do is to have a randon intercept for Rails. 

m<-lme(travel ~ Rail, data = Rail, random = ~ 1 | Rail)

You might want to try the documentation under C:\Program 
Files\R\rw1091\library\lme4\doc and also the book "Mixed-effect models in S 
and S-Plus" by Pinheiro and Bates.


>From: "angel hellraiser" <darkjacknife at hotmail.com>
>To: r-help at stat.math.ethz.ch
>Subject: [R] about lme
>Date: Fri, 06 Aug 2004 16:40:43 +0000
>Hi R-users:
>I've got a problem with lme.
>In Rail data,I try to model the next one
>lme(travel ~ Rail, data = Rail, random = ~ Rail | Rail)
>I want travel = Rail(i) + Rail(j) + epsilon(i,j)
>say, an effect fixed for every Rail and other for Rail(j) random, R says
>fewer observations or false convergence.
>Why can't I model an effect fixed and other random?
>Thanks in advance.
>My e-mail is darkjacknife at hotmail.com
>Descárgalo y pruébalo 2 meses gratis.
>R-help at stat.math.ethz.ch mailing list
>PLEASE do read the posting guide! 

More information about the R-help mailing list