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: [EMAIL PROTECTED] Fax 7 838 4155 > Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 1395 862 > > ______________________________________________ > R-help@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: [EMAIL PROTECTED] http://www.ms.unimelb.edu.au ______________________________________________ R-help@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