if (FALSE) {
data(finance, package = "smdata")
lme4fit <- glmer(corr ~ jmeth + (1 | item), data = finance,
family = binomial, nAGQ = 20)
# clusterwise scores
estfun(lme4fit, ranpar = "var")
}
Run the code above in your browser using DataLab