[R] help with lme function
Deepayan Sarkar
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)
>
>
> Regards,
> Olivier.
__________________________________________________
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
mailing list