(fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject),
data = sleepstudy))
print(fm2, digits = 10, corr = FALSE) # more precision; no corr.matrix
logLik(fm2)
(V2 <- vcov(fm2))
str(resid(fm2))
## Simulate 'Reaction' according to the fitted model:
dim(ss <- simulate(fm2, nsim = 200, seed = 101)) ## -> 180 x 200
Run the code above in your browser using DataLab