Computing Gaussian Likelihoods and Their Derivatives for General Linear Mixed Models
SIAM Journal on Scientific Computing1994Vol. 15(6), pp. 1294–1310
Citations Over TimeTop 10% of 1994 papers
Abstract
Algorithms are described for computing the Gaussian likelihood or restricted likelihood corresponding to a general linear mixed model. Included are arbitrary covariance structures for both the random effects and errors. Formulas are also given for the first and second derivatives of the likelihoods, thus enabling a Newton–Raphson implementation. The algorithms make heavy use of the Cholesky decomposition, the sweep operator, and the W-transformation. Also described are the modifications needed for variance profiling, Fisher scoring, and MIVQUE(0), as well as the computational order of the procedures.
Related Papers
- → Multilevel Models for Survival Analysis with Random Effects(2001)86 cited
- → Selecting the covariance structure for the seemingly unrelated regression models(2022)2 cited
- A mixed model for repeated split-plot data(2010)
- → Extending generalized linear models with random effects and components of dispersion = [Gegeneraliseerde lineaire modellen met extra stochastische termen en bijbehorende variantiecomponenten](1997)
- → STARTING VALUES FOR PROC MIXED WITH REPEATED MEASURES DATA(1999)