R help - May 2009 - Help to extract results from lme.object

Hi list,

Would someone please give me a hint how to extract the between and within
subject StdDev from "lme" fit? Below is the output from the fit, but I
couldn't
find the number 0.3927714 anywhere in the fitted object. Thanks in advance for
your help!


Random effects:
 Formula: ~1 | group
        (Intercept)  Residual
StdDev:   0.3927714 0.7857085



Best,

Auston

Auston wrote:
> Hi list,
>
> Would someone please give me a hint how to extract the between and within
> subject StdDev from "lme" fit? Below is the output from the fit,
but I couldn't
> find the number 0.3927714 anywhere in the fitted object. Thanks in advance
for
> your help!
>
>
> Random effects:
>  Formula: ~1 | group
>         (Intercept)  Residual
> StdDev:   0.3927714 0.7857085
>
>
>
> Best,
>
> Auston
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>   
Hi Auston,

Extracting those standard deviation values isn't so hard.

If you want to extract the variance on your parameter estimates, use

diag(lme.output$varFix) ;

If you want to extract the variance for the random effects, use

as.numeric(VarCorr(lme.output)[,"Variance"]) ;

At last, if you want your within-group error standard deviation, simply try

lme.output$sigma ;

Is there something missing?

Cheers,

-- 
*Luc Villandr?*
/Biostatistician
McGill University Health Center -
Montreal Children's Hospital Research Institute/

Dear Luc,

Thank you so much for your help! The VarCorr() did the trick. Many thanks!

Best,

Auston