# NOT RUN {
set.seed(123)
library(dplyr)
# linear mixed effects model
grouped_augment(
data = mutate(MASS::Aids2, interval = death - diag),
grouping.vars = sex,
..f = lme4::lmer,
formula = interval ~ age + (1 | status),
control = lme4::lmerControl(optimizer = "bobyqa")
)
# }
Run the code above in your browser using DataLab