Extracting variance components

Source: R/utility_functions.R

This function extracts the variance components from a mixed/hierarchical linear model fit using lmer.

varcomp.mer(object)

Arguments

object

a fitted model object of class mer or lmerMod.

Value

A named vector is returned. sigma2 denotes the residual variance. The other variance components are names D** where the trailing digits specify the of that variance component in the covariance matrix of the random effects.

Author

Adam Loy loyad01@gmail.com

Examples

data(sleepstudy, package = "lme4") 
fm1 <- lme4::lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
varcomp.mer(fm1)

#>     sigma2        D11        D21        D22 
#> 654.940008 612.100158   9.604409  35.071714