Hi all,
I have a random effect model that can be written as
y_{ij} = \beta +\alpha_i+ \epslion_{ij}
where \alpha_i ~ N(0,\sigma^2_effect) and \espilon_{ij} ~
N(0,\sigma^2_error)
and i compute
res<-lme(data~1,data=Data,random=~1 | veci)
The estimation of \sigma_error is given by res$sigma but
I don't understand how i can find the estimation of \sigma _effect with
the object res
(I would like to have access to the value of the estimation and not
print the result by using summary(res) or print(res) )
Could you help me ?
Best regards
--
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
