Is there an R package that can provide restricted maximum likelihood-empirical Bayes estimates for mixed model? For example, 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)

Yang and Land (2006) [mixed model approach to age-period-cohort analysis of repeated cross-sectional surveys]  used restricted maximum likelihood-empirical Bayes estimates to estimates model statistics and parameters.
Thank you


	[[alternative HTML version deleted]]

