[R] Extracting Variance components

[Andrew Robinson]

Murray,

you'll find it in 

VarCorr(fm1OrthF)

Cheers

Andrew

On Mon, Jun 05, 2006 at 04:29:48PM +1200, Murray Jorgensen wrote:
> I can ask my question using and example from Chapter 1 of Pinheiro & Bates.
> 
>  > # 1.4 An Analysis of Covariance Model
>  >
>  > OrthoFem <- Orthodont[ Orthodont$Sex == "Female", ]
>  > fm1OrthF <-
> +   lme( distance ~ age, data = OrthoFem, random = ~ 1 | Subject )
>  > summary( fm1OrthF )
> Linear mixed-effects model fit by REML
>   Data: OrthoFem
>         AIC     BIC    logLik
>    149.2183 156.169 -70.60916
> 
> Random effects:
>   Formula: ~1 | Subject
>          (Intercept)  Residual
> StdDev:     2.06847 0.7800331
> [...etc...]
> 
> I can extract the estimate of the variance component \sigma (0.7800331) via
> 
> sigma <- fm1OrthF$sigma
> 
> How do I extract the other component \sigma_b (2.06847) ?
> 
> Cheers,  Murray Jorgensen
> -- 
> Dr Murray Jorgensen      http://www.stats.waikato.ac.nz/Staff/maj.html
> Department of Statistics, University of Waikato, Hamilton, New Zealand
> Email: maj at waikato.ac.nz                                Fax 7 838 4155
> Phone  +64 7 838 4773 wk    Home +64 7 825 0441    Mobile 021 1395 862
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

-- 
Andrew Robinson  
Department of Mathematics and Statistics            Tel: +61-3-8344-9763
University of Melbourne, VIC 3010 Australia         Fax: +61-3-8344-4599
Email: a.robinson at ms.unimelb.edu.au         http://www.ms.unimelb.edu.au