Skip to content
Prev 66795 / 398525 Next

help im lme()

Hi R people!

I have a doubt in lme().
I use this model:

m2.lme<-lme(log(cmort)~idade+ano,random=~idade|ano,data=dados)

If i use summary I recive this output:

 > summary(m2.lme)
Linear mixed-effects model fit by REML
  Data: dados
        AIC      BIC    logLik
   1139.313 1170.554 -562.6563

Random effects:
  Formula: ~idade | ano
  Structure: General positive-definite, Log-Cholesky parametrization
             StdDev       Corr
(Intercept) 3.131076e-03 (Intr)
idade       3.515018e-05 -0.042
Residual    5.687940e-01

...etc...

I know the value of Residula random effects is m2.lme$sigma but how do I 
find the value of idade or intercept random effects


Thanks in advance

Bernardo Rangel Tura, MD, MSc
National Institute of Cardiology Laranjeiras
Rio de Janeiro Brazil