tests/varcorr.R

In lme4: Linear Mixed-Effects Models using 'Eigen' and S4

library(lme4)

if (.Platform$OS.type != "windows") {
data(Orthodont, package="nlme")
fm1 <- lmer(distance ~ age + (age|Subject), data = Orthodont)
VarCorr(fm1)

fm2ML <- lmer(diameter ~ 1 + (1|plate) + (1|sample), Penicillin, REML=0)
VarCorr(fm2ML)

gm1 <- glmer(cbind(incidence,size-incidence) ~ period + (1|herd),data=cbpp,
             family=binomial)
VarCorr(gm1)

cbpp$obs <- factor(seq(nrow(cbpp)))
gm2 <- update(gm1,.~.+(1|obs))
VarCorr(gm2)

if (FALSE) {
  ## testing lme4/lme4 incompatibility
##  library(lme4)
  VarCorr(fm1)
  lme4:::VarCorr.merMod(fm1) ## OK
}
} ## skip on windows (for speed)