Hi
I suspect that what you want to do is to have a randon intercept for Rails.
i.e.:
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.
Francisco
>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
>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide!
>http://www.R-project.org/posting-guide.html