inst/examples/eg_lmer.R

In ib: Bias Correction via Iterative Bootstrap

## linear mixed-effects regression
library(lme4)
fit_lmm <- lmer(Reaction ~ Days + (Days | Subject), data = sleepstudy, REML = FALSE)
fit_ib <- ib(fit_lmm)
summary(fit_ib)
## correct for variances and correlation
\dontrun{
fit_ib <- ib(fit_lmm, extra_param = TRUE)
summary(fit_ib)
}