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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.