library("lme4")
m.p <- lmer(weight ~ week + (1 | id) + (1 | week),
data=Pigs, REML=FALSE,
control=lmerControl(optimizer="bobyqa"))
summary(m.p)
cv(m.p, clusterVariables=c("id", "week"), k=10, seed=8469)
Run the code above in your browser using DataLab