Abstract:We derive a formula for optimal hard thresholding of the singular value decomposition in the presence of correlated additive noise; although it nominally involves unobservables, we show how to apply it even where the noise covariance structure is not a-priori known or is not independently estimable. The proposed method, which we call ScreeNOT, is a mathematically solid alternative to Cattell's ever-popular but vague Scree Plot heuristic from 1966. ScreeNOT has a surprising oracle property: it typically achieves exactly, in large finite samples, the lowest possible MSE for matrix recovery, on each given problem instance - i.e. the specific threshold it selects gives exactly the smallest achievable MSE loss among all possible threshold choices for that noisy dataset and that unknown underlying true low rank model. The method is computationally efficient and robust against perturbations of the underlying covariance structure. Our results depend on the assumption that the singular values of the noise have a limiting empirical distribution of compact support; this model, which is standard in random matrix theory, is satisfied by many models exhibiting either cross-row correlation structure or cross-column correlation structure, and also by many situations where there is inter-element correlation structure. Simulations demonstrate the effectiveness of the method even at moderate matrix sizes. The paper is supplemented by ready-to-use software packages implementing the proposed algorithm: package ScreeNOT in Python (via PyPI) and R (via CRAN).

Supplementary Material of Generalized Singular Value Thresholding

Generalized Singular Value Thresholding

Optimized Number of the Primary Singular Values for Image Reconstruction in Reflection Matrix Based Optical Coherence Tomography

A Hybrid Algorithm for Computing a Partial Singular Value Decomposition Satisfying a Given Threshold

Adaptive Hybrid Threshold Shrinkage Using Singular Value Decomposition

Accelerated singular value thresholding for matrix completion.

ScreeNOT: Exact MSE-Optimal Singular Value Thresholding in Correlated Noise

A Contour Integral-Based Algorithm for Computing Generalized Singular Values

Efficient Thresholded Correlation using Truncated Singular Value Decomposition

Generalized Singular Value Thresholding Operator to Affine Matrix Rank Minimization Problem

Comments on Efficient Singular Value Thresholding Computation

Improved upper bounds for the structured singular value

Target Recognition Based on Singular Value Decomposition in a Single-Pixel Non-Imaging System

A General-Thresholding Solution for ℓp (0 < p <1) Regularized CT Reconstruction

Some Software Packages for Partial SVD Computation

Further lower bounds for the smallest singular value and applications

On the Distribution of the Squared Generalized Singular Values and Its Applications.

Fast Singular Value Shrinkage with Chebyshev Polynomial Approximation Based on Signal Sparsity

New discovery of a conversion of generalized singular values

Supplementary Material to “ Convolutional Sparse Coding for Image Super-resolution ”

Singularity, Wielandt's Lemma and Singular Values