[R-sig-ME] Extracting variances from a mer object
Gang Chen
gangchen6 at gmail.com
Wed Jul 7 23:06:00 CEST 2010
This is a simple question. When fitting a model with lmer, how can I
extract all the variances from the output object? For example, with
the following analysis, I just want to extract the column of Variance
under the Random effects table - those 3 variance numbers: 0.63756 for
Subj, 23.03698 for Scan, and 199.50496 for Residual. VarCorr(fm) only
gives me the first two. I couldn't find any slots for this purpose
from str(fm) either.
> (fm<-lmer(Rep~(1|Subj)+(1|Scan), data=Model))
Linear mixed model fit by REML
Formula: Rep ~ (1 | Subj) + (1 | Scan)
Data: Model
AIC BIC logLik deviance REMLdev
166.4 170.3 -79.18 163.2 158.4
Random effects:
Groups Name Variance Std.Dev.
Subj (Intercept) 0.63756 0.79847
Scan (Intercept) 23.03698 4.79969
Residual 199.50496 14.12462
Number of obs: 20, groups: Subj, 10; Scan, 2
Fixed effects:
Estimate Std. Error t value
(Intercept) 27.603 4.642 5.946
Thanks,
Gang
More information about the R-sig-mixed-models
mailing list