You are not currently logged in.
Access JSTOR through your library or other institution:
Newton-Raphson and EM Algorithms for Linear Mixed-Effects Models for Repeated-Measures Data
Mary J. Lindstrom and Douglas M. Bates
Journal of the American Statistical Association
Vol. 83, No. 404 (Dec., 1988), pp. 1014-1022
Stable URL: http://www.jstor.org/stable/2290128
Page Count: 9
You can always find the topics here!Topics: Matrices, Algorithms, Covariance, Statistical models, Maximum likelihood estimation, Statistical estimation, Statistical variance, Data models, Estimation methods, Estimators
Were these topics helpful?See something inaccurate? Let us know!
Select the topics that are inaccurate.
Preview not available
We develop an efficient and effective implementation of the Newton-Raphson (NR) algorithm for estimating the parameters in mixed-effects models for repeated-measures data. We formulate the derivatives for both maximum likelihood and restricted maximum likelihood estimation and propose improvements to the algorithm discussed by Jennrich and Schluchter (1986) to speed convergence and ensure a positive-definite covariance matrix for the random effects at each iteration. We use matrix decompositions to develop efficient and computationally stable implementations of both the NR algorithm and an EM algorithm (Laird and Ware 1982) for this model. We compare the two methods (EM vs. NR) in terms of computational order and performance on two sample data sets and conclude that in most situations a well-implemented NR algorithm is preferable to the EM algorithm or EM algorithm with Aitken's acceleration. The term repeated measures refers to experimental designs where there are several individuals and several measurements taken on each individual. In the mixed-effects model each individual's vector of responses is modeled as a parametric function, where some of the parameters or "effects" are random variables with a multivariate normal distribution. This model has been successful because it can handle unbalanced data (different designs for different individuals), missing data (observations on all individuals are taken at the same design points, but some individuals have missing data), and jointly dependent random effects. The price for this flexibility is that the parameter estimates may be difficult to compute. We propose some new methods for implementing the EM and NR algorithms and draw conclusions about their performance. We also discuss extensions of the mixed-effects model to incorporate nonindependent conditional error structure and nested-type designs.
Journal of the American Statistical Association © 1988 American Statistical Association