# Linear mixed model example from lme4
data("sleepstudy", package = "lme4")
mod <- galamm(Reaction ~ Days + (Days | Subject), data = sleepstudy)
# Diagnostic plot
plot(mod)
plot(mod, residuals = "deviance")
Run the code above in your browser using DataLab