Abstract:Across many disciplines spanning from neuroscience and genomics to machine learning, atmospheric science, and finance, the problems of denoising large data matrices to recover hidden signals obscured by noise, and of estimating the structure of these signals, is of fundamental importance. A key to solving these problems lies in understanding how the singular value structure of a signal is deformed by noise. This question has been thoroughly studied in the well-known spiked matrix model, in which data matrices originate from low-rank signal matrices perturbed by additive noise matrices, in an asymptotic limit where matrix size tends to infinity but the signal rank remains finite. We first show, strikingly, that the singular value structure of large finite matrices (of size ∼1000) with even moderate-rank signals, as low as 10, is not accurately predicted by the finite-rank theory, thereby limiting the application of this theory to real data. To address these deficiencies, we analytically compute how the singular values and vectors of an arbitrary high-rank signal matrix are deformed by additive noise. We focus on an asymptotic limit corresponding to an extensive spike model, in which both the signal rank and the size of the data matrix tend to infinity at a constant ratio. We map out the phase diagram of the singular value structure of the extensive spike model as a joint function of signal strength and rank. We further exploit these analytics to derive optimal rotationally invariant denoisers to recover the hidden high-rank signal from the data, as well as optimal invariant estimators of the signal covariance structure. Our extensive-rank results yield several conceptual differences compared to the finite-rank case: (1) as signal strength increases, the singular value spectrum does not directly transition from a unimodal bulk phase to a disconnected phase, but instead there is a bimodal connected regime separating them; (2) the signal singular vectors can be partially estimated even in the unimodal bulk regime, and thus the transitions in the data singular value spectrum do not coincide with a detectability threshold for the signal singular vectors, unlike in the finite-rank theory; (3) signal singular values interact nontrivially to generate data singular values in the extensive-rank model, whereas they are noninteracting in the finite-rank theory; and (4) as a result, the more sophisticated data denoisers and signal covariance estimators we derive, which take into account these nontrivial extensive-rank interactions, significantly outperform their simpler, noninteracting, finite-rank counterparts, even on data matrices of only moderate rank. Overall, our results provide fundamental theory governing how high-dimensional signals are deformed by additive noise, together with practical formulas for optimal denoising and covariance estimation.

Fundamental limits of symmetric low-rank matrix estimation

Fundamental limits of Non-Linear Low-Rank Matrix Estimation

A multiscale cavity method for sublinear-rank symmetric matrix factorization

Low Rank Matrix Recovery with Adversarial Sparse Noise

Rank-one matrix estimation: analysis of algorithmic and information theoretic limits by the spatial coupling method

Matrix Inference in Growing Rank Regimes

Estimation of high-dimensional low-rank matrices

A New Statistic for Testing Covariance Equality in High-Dimensional Gaussian Low-Rank Models

Large information plus noise random matrix models and consistent subspace estimation in large sensor networks

Statistical and Computational Limits for Sparse Matrix Detection

Coherence-free Entrywise Estimation of Eigenvectors in Low-rank Signal-plus-noise Matrix Models

Information limits and Thouless-Anderson-Palmer equations for spiked matrix models with structured noise

Singular Vectors of Sums of Rectangular Random Matrices and Optimal Estimators of High-Rank Signals: The Extensive Spike Model

A Random Matrix Approach to Low-Multilinear-Rank Tensor Approximation

Universal Rank Inference via Residual Subsampling with Application to Large Networks

Completing Low-Rank Matrices with Corrupted Samples from Few Coefficients in General Basis

Singular vectors of sums of rectangular random matrices and optimal estimation of high-rank signals: The extensive spike model

Fundamental limits in structured principal component analysis and how to reach them

Optimal low-rank approximations of posteriors for linear Gaussian inverse problems on Hilbert spaces

Statistical limits of dictionary learning: random matrix theory and the spectral replica method

Simple bounds for recovering low-complexity models