[R] help on lme and variance estimation

[Olivier Martin]

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