1library(lme4) 2 3if (.Platform$OS.type != "windows") { 4data(Orthodont, package="nlme") 5fm1 <- lmer(distance ~ age + (age|Subject), data = Orthodont) 6VarCorr(fm1) 7 8fm2ML <- lmer(diameter ~ 1 + (1|plate) + (1|sample), Penicillin, REML=0) 9VarCorr(fm2ML) 10 11gm1 <- glmer(cbind(incidence,size-incidence) ~ period + (1|herd),data=cbpp, 12 family=binomial) 13VarCorr(gm1) 14 15cbpp$obs <- factor(seq(nrow(cbpp))) 16gm2 <- update(gm1,.~.+(1|obs)) 17VarCorr(gm2) 18 19if (FALSE) { 20 ## testing lme4/lme4 incompatibility 21## library(lme4) 22 VarCorr(fm1) 23 lme4:::VarCorr.merMod(fm1) ## OK 24} 25} ## skip on windows (for speed) 26