compare_eb_ls {HLMdiag}

Visually comparing shrinkage and LS estimates

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)