arrow
Return

Robust Linear Mixed Models Using Hierarchical Gamma-Divergence

delete2025-11-01
delete0
PRE
AI
S
Shonosuke Sugasawa *
F
Francis K. C. Hui
A
A. H. Welsh
DOI:10.1080/10618600.2025.2565235delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
Linear mixed models (LMMs) are a popular class of methods for analyzing longitudinal and clustered data. However, such models can be sensitive to outliers, and this can lead to biased inference on model parameters and inaccurate prediction of random effects if the data are contaminated. We propose a new approach to robust estimation and inference for LMMs using a hierarchical gamma-divergence, which offers an automated, data-driven approach to downweight the effects of outliers occurring in both the error and the random effects, using normalized powered density weights. For estimation and inference, we develop a computationally scalable minorization-maximization algorithm for the resulting objective function, along with a clustered bootstrap method for uncertainty quantification and a Hyvarinen score criterion for selecting a tuning parameter controlling the degree of robustness. Under suitable regularity conditions, we show the resulting robust estimates can be asymptotically controlled even under a heavy level of (covariate-dependent) contamination. Simulation studies demonstrate hierarchical gamma-divergence consistently outperforms several currently available methods for robustifying LMMs. We also illustrate the proposed method using data from a multi-center AIDS cohort study. Supplementary materials for this article are available online.
Keywords:
Clustered data
Divergence
Maximum likelihood estimation
MM-algorithm
Outliers
Random effects

Journal

J
Journal of Computational and Graphical Statistics
IF:
1.8
Papers:
138
Citations:
6.4K

Organization

A
Australian National University
Scholars:
2.1W
Papers: 2.3W
Citations: 3.9W
K
keio university
Scholars:
3.5K
Papers: 1.4K
Citations: 0