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