Robust and Scalable Sure Screening of Fixed effects in Ultrahigh-dimensional Linear Mixed Models

2026-06-26

AI Summary

This study addresses the challenges posed by exponentially growing fixed-effect covariates, dependent random effects, and data contamination in ultra-high-dimensional linear mixed models. The authors propose a robust and efficient variable screening method that decouples random effects via a proxy transformation and constructs marginal utilities based on minimum density power divergence, enabling robust marginal screening in the transformed model. This work represents the first integration of robust statistics with ultra-high-dimensional mixed-model screening, offering sure screening consistency, a high breakdown point, and bounded influence functions. It accommodates non-Gaussian errors and non-polynomial dimensionality growth, and allows for iterative refinement through prior information. Theoretical analysis establishes that all relevant variables are retained with exponentially high probability, and both simulations and real-data analysis of the ADNI2 dataset demonstrate superior performance over existing methods under data contamination.

Abstract

In modern applications of linear mixed models, the number of candidate fixed-effects covariates can grow exponentially with the sample size, while dependence induced by random effects and possible data contamination pose substantial challenges for existing variable screening methods. We propose a robust and computationally efficient sure screening procedure for identifying relevant fixed-effects covariates in ultrahigh-dimensional linear mixed models with known random effects. The proposed method leverages a proxy-based transformation to decouple dependence induced by random effects, enabling screening via marginal analysis in a transformed regression model. Robustness is achieved by constructing marginal utilities based on minimum density power divergence, yielding stability under data contamination and model misspecification without sacrificing scalability. The resulting procedure, termed DPD-SISP, is shown to retain all relevant covariates (sure screening property) with exponentially high probability under general conditions, allowing for non-Gaussian errors and nonpolynomial growth of dimensionality. In addition, DPD-SISP exhibits strong robustness properties supported by influence function and breakdown point analyses. The framework is further extended to incorporate prior information through conditional screening, mitigate correlation-induced masking via iterative refinement, and enable robust post-screening estimation of fixed effects. Extensive simulation studies demonstrate competitive performance of DPD-SISP under ideal settings and substantial gains in stability under data contamination. Its practical utility is illustrated through an application to high-dimensional longitudinal data from the ADNI2 study. The proposed framework thus provides a unified, robust, and scalable approach for variable screening in ultrahigh-dimensional linear mixed models.

Problem

• ultrahigh-dimensional
• linear mixed models
• fixed-effects screening
• data contamination
• dependence

Innovation

• sure screening
• linear mixed models
• robustness
• minimum density power divergence
• ultrahigh-dimensional