[R] help with lme function
deepayansarkar at yahoo.com
Fri Mar 29 18:17:04 CET 2002
--- Olivier Martin <olivier.martin at inrialpes.fr> wrote:
> Hi all,
> I have some difficulties with the lme function and so this is my problem.
> Supoose i have the following model
> y_(ijk)=beta_j + e_i + epsilon_(ijk)
> where beta_j are fixed effects, e_i is a random effect and
> epsilon_(ijk) is the error.
> If i want to estimate a such model, i execute
> >lme(y~vec.J , random~1 |vec .I )
> where y is the vector of my data, vec.J is a factor object and vec.I
> is the vector for the i indice.
> Now , this is an other model
> y_(ijk)=beta_j + e_(ij) + epsilon_(ijk)
> The parameters estimations can be resolved with
> >lme(y~vec.J , random~1 |vec .IJ )
> My problem concern the following model
> y_(ijk)=beta_j + e_i + e_(ij) + epsilon_(ijk)
> where e_i and e_(ij) are the random effects.
> In this case, what is the command to estimate the parameters?
random = ~1|vec.I/vec.J
(which is interpreted as the grouping being vec.I and vec.J %in% vec.I)
Yahoo! Greetings - send holiday greetings for Easter, Passover
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
More information about the R-help