Is there an R package that can provide empirical Bayes estimates for mixed model? In the example below, the parameters are estimated via ML (REML=FALSE)

library(foreign)
library(lme4)
 
pupcross<-read.dta("http://www.ats.ucla.edu/stat/stata/examples/mlm_ma_hox/pupcross.dta")
 
m1<-lmer(achiev ~ (1|sschool) + (1|pschool), pupcross, REML=FALSE)


Webster Kasongo
PhD Epidemiology program| Vanderbilt University
Nashville, TN
	[[alternative HTML version deleted]]

