Computation of Maximum Likelihood Estimates for Multiresponse Generalized Linear Mixed Models with Non-nested, Correlated Random Effects
arXiv:1403.7676 · doi:10.1016/j.csda.2013.11.019
Abstract
Estimation of generalized linear mixed models (GLMMs) with non-nested random effects structures requires approximation of high-dimensional integrals. Many existing methods are tailored to the low-dimensional integrals produced by nested designs. We explore the modifications that are required in order to adapt an EM algorithm with first-order and fully exponential Laplace approximations to a non-nested, multiple response model. The equations in the estimation routine are expressed as functions of the first four derivatives of the conditional likelihood of an arbitrary GLMM, providing a template for future applications. We apply the method to a joint Poisson-binary model for ranking sporting teams, and discuss the estimation of a correlated random effects model designed to evaluate the sensitivity of value-added models for teacher evaluation to assumptions about the missing data process. Source code in R is provided in the online supplementary material.
20 pages, 2 figures. Supplementary code in arXiv source package
References in corpus (4)
- Workload forecasting for a call center: Methodology and a case study
- Efficient Maximum Likelihood Estimation of Multiple Membership Linear Mixed Models, with an Application to Educational Value-Added Assessments
- A Correlated Random Effects Model for Nonignorable Missing Data in Value-Added Assessment of Teacher Effects
- The Sensitivity of College Football Rankings to Several Modeling Choices