--- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
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?
Regards,
Olivier.
--
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
Olivier MARTIN
PhD student phone : (33) 04 76 61 53 55
Projet IS2 06 08 67 93 42
INRIA Rhone-Alpes fax : (33) 04 76 61 54 77
655, Av. de l'Europe
Montbonnot e-mail: olivier.martin at inrialpes.fr
38334 Saint Ismier cedex web : http://www.inrialpes.fr/is2/people/omartin
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._