Scholion.ai

Testing in high-dimensional spiked models

Iain M. Johnstone, Alexei Onatski

The Annals of Statistics · 2020 · Volume 48 · Issue 3

DOI: 10.1214/18-AOS1697

Abstract

We consider the five classes of multivariate statistical problems identified by James (Ann. Math. Stat. 35 (1964) 475–501), which together cover much of classical multivariate analysis, plus a simpler limiting case, symmetric matrix denoising. Each of James’ problems involves the eigenvalues of $E^{-1}H$ where $H$ and $E$ are proportional to high-dimensional Wishart matrices. Under the null hypothesis, both Wisharts are central with identity covariance. Under the alternative, the noncentrality or the covariance parameter of $H$ has a single eigenvalue, a spike, that stands alone. When the spike is smaller than a case-specific phase transition threshold, none of the sample eigenvalues separate from the bulk, making the testing problem challenging. Using a unified strategy for the six cases, we show that the log likelihood ratio processes parameterized by the value of the subcritical spike converge to Gaussian processes with logarithmic correlation. We then derive asymptotic power envelopes for tests for the presence of a spike.

Keywords

62E20, 62H15, canonical correlations, hypergeometric function, likelihood ratio test, Matrix denoising, multiple response regression, principal components analysis

Enable JavaScript for citation tools and additional paper interactions.