Multivariate MM-estimators with auxiliary Scale for Linear Models with Structured Covariance Matrices
Hendrik Paul Lopuhaa
TL;DR
This work develops a unified MM-estimation framework with auxiliary scale for balanced linear models with structured covariance, extending to linear mixed effects and several multivariate models. By introducing MM-functionals based on two losses $\rho_0$ and $\rho_1$ and an auxiliary scale, it establishes existence, continuity, consistency, and asymptotic normality under mild distributional assumptions, far beyond elliptically contoured densities. The authors derive global breakdown points and local influence functions, including explicit forms in the general case and under elliptical models, linking to existing multivariate MM/S-functionals and CM-functionals. A comprehensive application demonstrates practical tuning via Tukey biweight losses, quantifying the trade-off between efficiency and robustness through scalar indices, and providing guidance for robust covariance estimation in complex multivariate settings.
Abstract
We provide a unified approach to MM-estimation with auxiliary scale for balanced linear models with structured covariance matrices. This approach leads to estimators that are highly robust against outliers and highly efficient for normal data. These properties not only hold for estimators of the regression parameter, but also for estimators of scale invariant transformations of the variance parameters. Of main interest are MM-estimators for linear mixed effects models, but our approach also includes MM-estimators in several other standard multivariate models. We provide sufficient conditions for the existence of MM-functionals and MM-estimators, establish asymptotic properties such as consistency and asymptotic normality, and derive their robustness properties in terms of breakdown point and influence function. All the results are obtained for general identifiable covariance structures and are established under mild conditions on the distribution of the observations, which goes far beyond models with elliptically contoured densities.
