Sufficient Ensemble Size for Random Matrix Theory based Handling of Singular Covariance Matrices

Research output: Contribution to journalArticle

Authors

Colleges, School and Institutes

Abstract

Singular covariance matrices are frequently encountered in both machine learning and optimization problems, most commonly due to high dimensionality of data and insufficient sample sizes. Among many methods of regularization, here we focus on a relatively recent random matrix theoretic approach, the idea of which is to create well-conditioned approximations of a singular covariance matrix and its inverse by taking the expectation of its random projections. We are interested in the error of a Monte Carlo implementation of this approach, which allows subsequent parallel processing in low dimensions in practice. We find that O(d) random projections, where d is the size of the original matrix,
are sufficient for the Monte Carlo error to become negligible, in the sense of expected spectral norm difference, for both covariance and inverse covariance approximation, in the latter case under mild assumptions.

Details

Original languageEnglish
Number of pages21
JournalAnalysis and Applications
Publication statusAccepted/In press - 15 Apr 2020

Keywords

  • Singular covariance, precision matrix, curse of dimensionality, random projections;, Monte Carlo error