[R] Extracting variance components in lme

Consider the output for the inroductory Rail example in Mixed Effects Models in S and S-PLUS by Pinheiro and Bates: summary(fm1Rail.lme) Linear mixed-effects model fit by REML Data: Rail AIC BIC logLik 128.177 130.6766 -61.0885 Random effects: Formula: ~1 | Rail

[R] Extracting Variance Components fro lme

It seems that what I need to get the within group component is as.numeric(VarCorr(fm1Rail.lme)[2,2]) thanks to Bert Gunter and Peter Alspach. Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton,

Re: [R] Extracting variance components in lme

v - VarCorr(fm1Rail.lme) str(v) # get an idea of how v is structured. This suggests: as.numeric(v[1, 2]) [1] 24.80547 There may be easier and better ways HTH, Simon. At 11:02 AM 3/11/2005, you wrote: Consider the output for the inroductory Rail example in Mixed Effects Models in S and

Re: [R] Extracting variance components in lme

Woops! I should have written: as.numeric(VarCorr(fm1Rail.lme)[1,2]) for the within component. -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: [EMAIL PROTECTED]