Maximum likelihood computations with repeated measures: application of the EM algorithm

Abstract The purpose of this article is to consider the use of the EM algorithm (Dempster, Laird, and Rubin 1977) for both maximum likelihood (ML) and restricted maximum likelihood (REML) estimation in a general repeated measures setting using a multivariate normal data model with linear mean and covariance structure (Anderson 1973). Several models and methods of analysis have been proposed in recent years for repeated measures data; Ware (1985) presented an overview. Because the EM algorithm is a general-purpose, iterative method for computing ML estimates with incomplete data, it has often been used in this particular setting (Dempster et al. 1977; Andrade and Helms 1984; Jennrich and Schluchter 1985). There are two apparently different approaches to using the EM algorithm in this setting. In one application, each experimental unit is observed under a standard protocol specifying measurements at each of n occasions (or under n conditions), and incompleteness implies that the number of measurements actua...

[1]  R. Jennrich,et al.  Unbalanced repeated-measures models with structured covariance matrices. , 1986, Biometrics.

[2]  N. Laird,et al.  An algorithmic approach for the fitting of a general mixed ANOVA model appropriate in longitudinal settings , 1986 .

[3]  Gregory C. Reinsel,et al.  Mean Squared Error Properties of Empirical Bayes Estimators in a Multivariate Random Effects General Linear Model , 1985 .

[4]  J. Ware.  Linear Models for the Analysis of Longitudinal Studies , 1985 .

[5]  Gregory C. Reinsel,et al.  Estimation and Prediction in a Multivariate Random Effects Generalized Linear Model , 1984 .

[6]  Ronald W. Helms,et al.  Maximum likelihood estimates in the multivariate normal with patterned mean and covariance via the em algorithm , 1984 .

[7]  J. Berger,et al.  Empirical Bayes Estimation of Rates in Longitudinal Studies , 1983 .

[8]  J. Ware,et al.  Random-effects models for longitudinal data. , 1982, Biometrics.

[9]  G. Reinsel,et al.  Multivariate Repeated-Measurement or Growth Curve Models with Multivariate Random-Effects Covariance Structure , 1982 .

[10]  T. Louis.  Finding the Observed Information Matrix When Using the EM Algorithm , 1982 .

[11]  D. Rubin,et al.  Estimation in Covariance Components Models , 1981 .

[12]  J. Miller,et al.  Explicit Maximum Likelihood Estimates from Balanced Data in the Mixed Model of the Analysis of Variance , 1980 .

[13]  Ted H. Szatrowski,et al.  Necessary and Sufficient Conditions for Explicit Solutions in the Multivariate Normal Estimation Problem for Patterned Means and Covariances , 1980 .

[14]  V. Klema.  LINPACK user's guide , 1980 .

[15]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[16]  D. Harville.  Maximum Likelihood Approaches to Variance Component Estimation and to Related Problems , 1977 .

[17]  E. Beale,et al.  Missing Values in Multivariate Analysis , 1975 .

[18]  T. W. Anderson.  Asymptotically Efficient Estimation of Covariance Matrices with Linear Structure , 1973 .

[19]  M. Woodbury.  A missing information principle: theory and applications , 1972 .

[20]  T. W. Anderson,et al.  Statistical analysis of time series , 1972 .

[21]  Curtis F. Gerald,et al.  APPLIED NUMERICAL ANALYSIS , 1972, The Mathematical Gazette.

[22]  D. M. Allen,et al.  Analysis of growth and dose response curves. , 1969, Biometrics.

[23]  C. G. Khatri,et al.  A note on a manova model applied to problems in growth curve , 1966 .

[24]  R. Potthoff,et al.  A generalized multivariate analysis of variance model useful especially for growth curve problems , 1964 .