(fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
set.seed(101)
samp1 <- mcmcsamp(fm1, n = 1000)
if (require("coda", quietly = TRUE, character.only = TRUE)) {
densityplot(samp1)
qqmath(samp1)
xyplot(samp1, scales = list(x = list(axs = 'i')))
print(summary(samp1))
print(autocorr.diag(samp1))
}
(eDF <- mean(samp1[,"deviance"]) - deviance(fm1, REML=FALSE)) # potentially useful approximate D.F.
stopifnot(abs(eDF - 7) < 1)
Run the code above in your browser using DataCamp Workspace