Constrained S-estimators for linear mixed effects models with covariance components.

Linear mixed effects (LME) models are increasingly used for analyses of biological and biomedical data. When the multivariate normal assumption is not adequate for an LME model, then a robust estimation approach is preferable to the maximum likelihood one. M-estimators were considered before for robust estimation of the LME models, and recently a constrained S-estimator was proposed. This S-estimator cannot be applied directly to LME models with correlated error terms and vector random effects with correlated dimensions. Therefore, a modification is proposed, which extends application of the constrained S-estimator to the LME models for multivariate responses with correlated dimensions and to longitudinal data. Also, a new computational algorithm is developed for computing constrained S-estimators. Performance of the S-estimators based on the original Tukey's biweight and translated biweight is evaluated in a small simulation study with repeated multivariate responses with correlated dimensions. The proposed methodology is applied to jointly analyze repeated measures on three cholesterol components, HDL, LDL, and triglycerides.

[1]  David M. Rocke.  Robustness properties of S-estimators of multivariate location and shape in high dimension , 1996 .

[2]  R. Maronna.  Robust $M$-Estimators of Multivariate Location and Scatter , 1976 .

[3]  J. Tukey,et al.  The Fitting of Power Series, Meaning Polynomials, Illustrated on Band-Spectroscopic Data , 1974 .

[4]  Pierre Duchesne,et al.  Robust estimation of the SUR model , 2000 .

[5]  V. Yohai,et al.  Robust regression with both continuous and categorical predictors , 2000 .

[6]  D. Ruppert.  Computing S Estimators for Regression and Multivariate Location/Dispersion , 1992 .

[7]  Peter J. Rousseeuw,et al.  An Algorithm for Positive-Breakdown Regression Based on Concentration Steps , 2000 .

[8]  E. Vonesh,et al.  Efficient inference for random-coefficient growth curve models with unbalanced data. , 1987, Biometrics.

[9]  R. Huggins,et al.  A Robust Approach to the Analysis of Repeated Measures , 1993 .

[10]  Peter J. Rousseeuw,et al.  ROBUST REGRESSION BY MEANS OF S-ESTIMATORS , 1984 .

[11]  Stephane Heritier,et al.  Robust Alternatives to the F‐Test in Mixed Linear Models Based on MM‐Estimates , 2007, Biometrics.

[12]  Maria-Pia Victoria-Feser,et al.  High-Breakdown Inference for Mixed Linear Models , 2006 .

[13]  Stefan Van Aelst,et al.  MULTIVARIATE REGRESSION S-ESTIMATORS FOR ROBUST ESTIMATION AND INFERENCE , 2005 .

[14]  M. Sperling,et al.  Effects of antiepileptic drugs on lipids, homocysteine, and C‐reactive protein , 2009, Annals of neurology.

[15]  H. P. Lopuhaä.  On the relation between S-estimators and M-estimators of multivariate location and covariance , 1989 .

[16]  E. Ronchetti,et al.  Robust Bounded-Influence Tests in General Parametric Models , 1994 .

[17]  V. Yohai,et al.  A Fast Algorithm for S-Regression Estimates , 2006 .

[18]  W. Fung,et al.  High Breakdown Estimation for Multiple Populations with Applications to Discriminant Analysis , 2000 .

[19]  A. Welsh,et al.  Robust Restricted Maximum Likelihood in Mixed Linear Models , 1995 .

[20]  Alice Richardson,et al.  Bounded Influence Estimation in the Mixed Linear Model , 1997 .