[R] help with lme function

[Deepayan Sarkar]

--- 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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._