inst/testdata/harmel_density.R

hsb <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/hsb.csv')
source("strip_env.R")
library(lme4)
library(lattice)

m31 <- lmer(math ~ ses*meanses + (ses | sch.id), data = hsb)
p <-  strip_profile(profile(m31))
densityplot(p)
saveRDS(p, file="harmel_profile.rds", version=2)
lme4/lme4 documentation built on May 4, 2024, 6:46 p.m.