670
92
v1v2v3v4 (latest)

On the sample covariance matrix estimator of reduced effective rank population matrices, with applications to fPCA

Abstract

This work provides a unified analysis of the properties of the sample covariance matrix Σn\Sigma_n over the class of p×pp\times p population covariance matrices Σ\Sigma of reduced effective rank re(Σ)r_e(\Sigma). This class includes scaled factor models and covariance matrices with decaying spectrum. We consider re(Σ)r_e(\Sigma) as a measure of matrix complexity, and obtain sharp minimax rates on the operator and Frobenius norm of ΣnΣ\Sigma_n-\Sigma, as a function of re(Σ)r_e(\Sigma) and Σ2\|\Sigma\|_2, the operator norm of Σ\Sigma. With guidelines offered by the optimal rates, we define classes of matrices of reduced effective rank over which Σn\Sigma_n is an accurate estimator. Within the framework of these classes, we perform a detailed finite sample theoretical analysis of the merits and limitations of the empirical scree plot procedure routinely used in PCA. We show that identifying jumps in the empirical spectrum that consistently estimate jumps in the spectrum of Σ\Sigma is not necessarily informative for other goals, for instance for the selection of those sample eigenvalues and eigenvectors that are consistent estimates of their population counterparts. The scree plot method can still be used for selecting consistent eigenvalues, for appropriate threshold levels. We provide a threshold construction and also give a rule for checking the consistency of the corresponding sample eigenvectors. We specialize these results and analysis to population covariance matrices with polynomially decaying spectra, and extend it to covariance operators with polynomially decaying spectra. An application to fPCA illustrates how our results can be used in functional data analysis.

View on arXiv
Comments on this paper

We use cookies and other tracking technologies to improve your browsing experience on our website, to show you personalized content and targeted ads, to analyze our website traffic, and to understand where our visitors are coming from. See our policy.