Mixed Models: Modelling Covariance Structure in the Analysis of Repeated Measures Data

The term 'repeated measures' refers to data with multiple observations on the same sampling unit. In most cases, the multiple observations are taken over time, but they could be over space. It is usually plausible to assume that observations on the same unit are correlated. Hence, statistical analysis of repeated measures data must address the issue of covariation between measures on the same unit. Until recently, analysis techniques available in computer software only offered the user limited and inadequate choices. One choice was to ignore covariance structure and make invalid assumptions. Another was to avoid the covariance structure issue by analysing transformed data or making adjustments to otherwise inadequate analyses. Ignoring covariance structure may result in erroneous inference, and avoiding it may result in inefficient inference. Recently available mixed model methodology permits the covariance structure to be incorporated into the statistical model. The MIXED procedure of the SAS((R)) System provides a rich selection of covariance structures through the RANDOM and REPEATED statements. Modelling the covariance structure is a major hurdle in the use of PROC MIXED. However, once the covariance structure is modelled, inference about fixed effects proceeds essentially as when using PROC GLM. An example from the pharmaceutical industry is used to illustrate how to choose a covariance structure. The example also illustrates the effects of choice of covariance structure on tests and estimates of fixed effects. In many situations, estimates of linear combinations are invariant with respect to covariance structure, yet standard errors of the estimates may still depend on the covariance structure.

[1]  D. E. Johnson,et al.  Analysis of Messy Data Volume I: Designed Experiments , 1985 .

[2]  P. Diggle Analysis of Longitudinal Data , 1995 .

[3]  P. Cornelius,et al.  Approximate F-tests of multiple degree of freedom hypotheses in generalized least squares analyses of unbalanced split-plot experiments , 1996 .

[4]  G. Robinson That BLUP is a Good Thing: The Estimation of Random Effects , 1991 .

[5]  P. Diggle An approach to the analysis of repeated measurements. , 1988, Biometrics.

[6]  Simo Puntanen,et al.  The Equality of the Ordinary Least Squares Estimator and the Best Linear Unbiased Estimator , 1989 .

[7]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[8]  H. Huynh,et al.  Estimation of the Box Correction for Degrees of Freedom from Sample Data in Randomized Block and Split-Plot Designs , 1976 .

[9]  R. Wolfinger Covariance structure selection in general mixed models , 1993 .

[10]  Franklin A. Graybill,et al.  Theory and Application of the Linear Model , 1976 .

[11]  Kathryn S. Dawson,et al.  Two Graphical Techniques Useful in Detecting Correlation Structure in Repeated Measures Data , 1997 .

[12]  L. Wasserman,et al.  A Reference Bayesian Test for Nested Hypotheses and its Relationship to the Schwarz Criterion , 1995 .

[13]  D. Harville,et al.  Some Bayesian and Non-Bayesian Procedures for the Analysis of Comparative Experiments and for Small-Area Estimation: Computational Aspects, Frequentist Properties, and Relationships , 1991 .

[14]  Nan M. Laird,et al.  Using the General Linear Mixed Model to Analyse Unbalanced Repeated Measures and Longitudinal Data , 1997 .

[15]  Robert Haining,et al.  Statistics for spatial data: by Noel Cressie, 1991, John Wiley & Sons, New York, 900 p., ISBN 0-471-84336-9, US $89.95 , 1993 .

[16]  Geert Molenberghs,et al.  Linear Mixed Models in Practice , 1997 .

[17]  H. Akaike A new look at the statistical model identification , 1974 .