Fixed and Random Effects Selection in Mixed Effects Models

We consider selecting both fixed and random effects in a general class of mixed effects models using maximum penalized likelihood (MPL) estimation along with the smoothly clipped absolute deviation (SCAD) and adaptive least absolute shrinkage and selection operator (ALASSO) penalty functions. The MPL estimates are shown to possess consistency and sparsity properties and asymptotic normality. A model selection criterion, called the IC_Qstatistic, is proposed for selecting the penalty parameters (Ibrahim, Zhu, and Tang, 2008,Journal of the American Statistical Association103, 1648-1658). The variable selection procedure based on IC_Qis shown to consistently select important fixed and random effects. The methodology is very general and can be applied to numerous situations involving random effects, including generalized linear mixed models. Simulation studies and a real data set from a Yale infant growth study are used to illustrate the proposed methodology.