ExSIS: Extended Sure Independence Screening for Ultrahigh-dimensional Linear Models

TL;DR

This paper extends correlation-based variable screening methods to arbitrary linear models, providing theoretical conditions for successful screening and demonstrating the potential to reduce ultrahigh-dimensional models to near-sample size dimensions.

Contribution

It introduces a generalized screening condition applicable to any linear model and post-screening inference, broadening the scope beyond prior assumptions.

Findings

The screening condition is sufficient for successful variable selection.

Under certain conditions, models can be reduced to nearly the sample size.

Specialization to sub-Gaussian models validates the generalization.

Abstract

Statistical inference can be computationally prohibitive in ultrahigh-dimensional linear models. Correlation-based variable screening, in which one leverages marginal correlations for removal of irrelevant variables from the model prior to statistical inference, can be used to overcome this challenge. Prior works on correlation-based variable screening either impose statistical priors on the linear model or assume specific post-screening inference methods. This paper first extends the analysis of correlation-based variable screening to arbitrary linear models and post-screening inference techniques. In particular, (i) it shows that a condition---termed the screening condition---is sufficient for successful correlation-based screening of linear models, and (ii) it provides insights into the dependence of marginal correlation-based screening on different problem parameters. Numerical experiments confirm that these insights are not mere artifacts of analysis; rather, they are reflective of the challenges associated with marginal correlation-based variable screening. Second, the paper explicitly derives the screening condition for arbitrary (random or deterministic) linear models and, in the process, it establishes that---under appropriate conditions---it is possible to reduce the dimension of an ultrahigh-dimensional, arbitrary linear model to almost the sample size even when the number of active variables scales almost linearly with the sample size. Third, it specializes the screening condition to sub-Gaussian linear models and contrasts the final results to those existing in the literature. This specialization formally validates the claim that the main result of this paper generalizes existing ones on correlation-based screening.