Using SAS PROC MIXED to Fit Multilevel Models, Hierarchical Models, and Individual Growth Models

SAS PROC MIXED is a flexible program suitable for fitting multilevel models, hierarchical linear models, and individual growth models. Its position as an integrated program within the SAS statistical package makes it an ideal choice for empirical researchers and applied statisticians seeking to do data reduction, management, and analysis within a single statistical package. Because the program was developed from the perspective of a “mixed” statistical model with both random and fixed effects, its syntax and programming logic may appear unfamiliar to users in education and the social and behavioral sciences who tend to express these models as multilevel or hierarchical models. The purpose of this paper is to help users familiar with fitting multilevel models using other statistical packages (e.g., HLM, MLwiN, MIXREG) add SAS PROC MIXED to their array of analytic options. The paper is written as a step-by-step tutorial that shows how to fit the two most common multilevel models: (a) school effects models, designed for data on individuals nested within naturally occurring hierarchies (e.g., students within classes); and (b) individual growth models, designed for exploring longitudinal data (on individuals) over time. The conclusion discusses how these ideas can be extended straighforwardly to the case of three level models. An appendix presents general strategies for working with multilevel data in SAS and for creating data sets at several levels.

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

[2]  I G Kreft,et al.  The Effect of Different Forms of Centering in Hierarchical Linear Models. , 1995, Multivariate behavioral research.

[3]  Russell D. Wolfinger,et al.  Analysis Issues in the Evaluation of Community Trials - Progress Toward Solutions in Sas/Stat Mixed , 1994 .

[4]  Jan de Leeuw,et al.  Review of Five Multilevel Analysis Programs: BMDP-5V, GENMOD, HLM, ML3, VARCL , 1994 .

[5]  John B. Willett,et al.  Questions and Answers in the Measurement of Change , 1988 .

[6]  John M. Ferron,et al.  Moving between Hierarchical Modeling Notations , 1997 .

[7]  H Goldstein,et al.  Multilevel time series models with applications to repeated measures data. , 1994, Statistics in medicine.

[8]  Joop J. Hox,et al.  Applied Multilevel Analysis. , 1995 .

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

[10]  R. Wolfinger Heterogeneous Variance-Covariance Structures for Repeated Measures , 1996 .

[11]  Ita G. G. Kreft,et al.  Multilevel Analysis Methods , 1994 .

[12]  R. Littell SAS System for Mixed Models , 1996 .

[13]  Anthony S. Bryk,et al.  Hierarchical Linear Models: Applications and Data Analysis Methods , 1992 .

[14]  Roel Bosker,et al.  Modeled Variance in Two-Level Models , 1994 .

[15]  John M. Ferron Teacher’s Corner: Moving Between Hierarchical Modeling Notations , 1997 .

[16]  Jan de Leeuw,et al.  Review of Five Multilevel Analysis Programs: BMDP-5V, GENMOD, HLM, ML3, VARCL@@@BMDP-5V@@@GENMOD@@@HLM, Version 2.1@@@ML3, Version 2.2@@@VARCL , 1994 .

[17]  R. Mclean,et al.  A Unified Approach to Mixed Linear Models , 1991 .

[18]  D. Hedeker,et al.  MIXREG: a computer program for mixed-effects regression analysis with autocorrelated errors. , 1996, Computer methods and programs in biomedicine.

[19]  Ita G. G. Kreft,et al.  Special Issue: Hierarchical Linear Models: Problems and Prospects , 1995 .

[20]  D. VanLeeuwen,et al.  A note on the covariance structure in a linear model , 1997 .

[21]  H. Goldstein Multilevel Statistical Models , 2006 .