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

DOI10.3150/22-BEJ1463arXiv1909.06120MaRDI QIDQ6325231FDOQ6325231

N. Benjamin Erichson, Michael W. Mahoney, Miles E. Lopes

Publication date: 13 September 2019

Abstract: 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 $hatSigma$ denote the sample covariance matrix of $n$ observations in $mathbb{R}^p$ that arise from a population matrix $Sigma$, and let $T_n=sqrtn|hatSigma-Sigma{|}_{extop}$. In the setting where the eigenvalues of $Sigma$ have a decay profile of the form $$, 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 $$, 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.

Mathematics Subject Classification ID

Nonparametric inference (62Gxx) Multivariate analysis (62Hxx) Numerical linear algebra (65Fxx)

This page was built for publication: Bootstrapping the Operator Norm in High Dimensions: Error Estimation for Covariance Matrices and Sketching

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6325231)