Bootstrapping the Operator Norm in High Dimensions: Error Estimation for Covariance Matrices and Sketching

09/13/2019
by   Miles E. Lopes, et al.
0

Although the operator (spectral) norm is one of the most widely used metrics for covariance estimation, comparatively little is known about the fluctuations of error in this norm. To be specific, let Σ̂ denote the sample covariance matrix of n observations in R^p that arise from a population matrix Σ, and let T_n=√(n)Σ̂-Σ_op. In the setting where the eigenvalues of Σ have a decay profile of the form λ_j(Σ) j^-2β, we analyze how well the bootstrap can approximate the distribution of T_n. Our main result shows that up to factors of log(n), the bootstrap can approximate the distribution of T_n at the dimension-free rate of n^-β-1/2/6β+4, with respect to the Kolmogorov metric. Perhaps surprisingly, a result of this type appears to be new even in settings where p< n. More generally, we discuss the consequences of this result beyond covariance matrices and show how the bootstrap can be used to estimate the errors of sketching algorithms in randomized numerical linear algebra (RandNLA). An illustration of these ideas is also provided with a climate data example.

READ FULL TEXT

page 19

page 20

research
08/05/2022

Improved Rates of Bootstrap Approximation for the Operator Norm: A Coordinate-Free Approach

Let Σ̂=1/n∑_i=1^n X_i⊗ X_i denote the sample covariance operator of cent...
research
04/15/2021

Rates of Bootstrap Approximation for Eigenvalues in High-Dimensional PCA

In the context of principal components analysis (PCA), the bootstrap is ...
research
02/23/2020

Hypothesis testing for eigenspaces of covariance matrix

Eigenspaces of covariance matrices play an important role in statistical...
research
08/26/2022

Quantitative limit theorems and bootstrap approximations for empirical spectral projectors

Given finite i.i.d. samples in a Hilbert space with zero mean and trace-...
research
07/27/2022

Exact spectral norm error of sample covariance

Let X_1,…,X_n be i.i.d. centered Gaussian vectors in ℝ^p with covariance...
research
06/25/2021

The conjugate gradient algorithm on a general class of spiked covariance matrices

We consider the conjugate gradient algorithm applied to a general class ...
research
09/08/2022

A Bootstrap Method for Spectral Statistics in High-Dimensional Elliptical Models

Although there is an extensive literature on the eigenvalues of high-dim...

Please sign up or login with your details

Forgot password? Click here to reset