compare_eb_ls: Visually comparing shrinkage and LS estimates

In HLMdiag: Diagnostic Tools for Hierarchical (Multilevel) Linear Models

Description

This function creates a plot (using qplot()) where the shrinkage estimate appears on the horizontal axis and the LS estimate appears on the vertical axis.

Usage

compare_eb_ls(eb, ols, identify = FALSE, silent = TRUE, ...)

Arguments

eb	a matrix of random effects
ols	a matrix of the OLS estimates found using random_ls_coef
identify	the percentage of points to identify as unusual, FALSE if you do not want the points identified.
silent	logical: should the list of data frames used to make the plots be suppressed.
...	other arguments to be passed to qplot()

Author(s)

Adam Loy [e-mail]

Examples

wages.fm1 <- lme4::lmer(lnw ~ exper + (exper | id), data = wages)
wages.sepLM <- adjust_lmList(lnw ~ exper | id, data = wages)
rancoef.eb <- coef(wages.fm1)$id
rancoef.ols <- coef(wages.sepLM)
compare_eb_ls(eb = rancoef.eb, ols = rancoef.ols, identify = 0.01)

HLMdiag documentation built on May 2, 2021, 9:06 a.m.