Normal approximation and concentration of spectral projectors of sample covariance

preprint2015arXivOpen access

Abstract

Let X,X1,…,XnX,X_1,\dots, X_n be i.i.d. Gaussian random variables in a separable Hilbert space H{\mathbb H} with zero mean and covariance operator Σ=E(X⊗X),Σ={\mathbb E}(X\otimes X), and let Σ^:=n−1∑j=1n(Xj⊗Xj)\hat Σ:=n^{-1}\sum_{j=1}^n (X_j\otimes X_j) be the sample (empirical) covariance operator based on (X1,…,Xn).(X_1,\dots, X_n). Denote by PrP_r the spectral projector of ΣΣ corresponding to its rr-th eigenvalue μrμ_r and by P^r\hat P_r the empirical counterpart of Pr.P_r. The main goal of the paper is to obtain tight bounds on sup⁡x∈R∣P{∥P^r−Pr∥22−E∥P^r−Pr∥22Var1/2(∥P^r−Pr∥22)≤x}−Φ(x)∣, \sup_{x\in {\mathbb R}} \left|{\mathbb P}\left\{\frac{\|\hat P_r-P_r\|_2^2-{\mathbb E}\|\hat P_r-P_r\|_2^2}{{\rm Var}^{1/2}(\|\hat P_r-P_r\|_2^2)}\leq x\right\}-Φ(x)\right|, where ∥⋅∥2\|\cdot\|_2 denotes the Hilbert--Schmidt norm and ΦΦ is the standard normal distribution function. Such accuracy of normal approximation of the distribution of squared Hilbert--Schmidt error is characterized in terms of so called effective rank of ΣΣ defined as r(Σ)=tr(Σ)∥Σ∥∞,{\bf r}(Σ)=\frac{{\rm tr}(Σ)}{\|Σ\|_{\infty}}, where tr(Σ){\rm tr}(Σ) is the trace of ΣΣ and ∥Σ∥∞\|Σ\|_{\infty} is its operator norm, as well as another parameter characterizing the size of Var(∥P^r−Pr∥22).{\rm Var}(\|\hat P_r-P_r\|_2^2). Other results include non-asymptotic bounds and asymptotic representations for the mean squared Hilbert--Schmidt norm error E∥P^r−Pr∥22{\mathbb E}\|\hat P_r-P_r\|_2^2 and the variance Var(∥P^r−Pr∥22),{\rm Var}(\|\hat P_r-P_r\|_2^2), and concentration inequalities for ∥P^r−Pr∥22\|\hat P_r-P_r\|_2^2 around its expectation.

Explore this paper’s authors, topics and related work.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Reviews 0

Write a reviewWrite

No reviews yet.

Discussion 0

Add a commentWrite

No comments yet.