A RANDOMIZED ALGORITHM FOR PRINCIPAL COMPONENT ANALYSIS

Principal component analysis (PCA) requires the computation of a low-rank approximation to a matrix containing the data being analyzed. In many applications of PCA, the best possible accuracy of any rank-deficient approximation is at most a few digits (measured in the spectral norm, relative to the...

Full description

Saved in:
Bibliographic Details
Published inSIAM journal on matrix analysis and applications Vol. 31; no. 3; pp. 1100 - 1124
Main Authors ROKHLIN, Vladimir, SZLAM, Arthur, TYGERT, Mark
Format Journal Article
LanguageEnglish
Published Philadelphia, PA Society for Industrial and Applied Mathematics 01.01.2009
Subjects
Online AccessGet full text

Cover

Loading…
More Information
Summary:Principal component analysis (PCA) requires the computation of a low-rank approximation to a matrix containing the data being analyzed. In many applications of PCA, the best possible accuracy of any rank-deficient approximation is at most a few digits (measured in the spectral norm, relative to the spectral norm of the matrix being approximated). In such circumstances, efficient algorithms have not come with guarantees of good accuracy, unless one or both dimensions of the matrix being approximated are small. We describe an efficient algorithm for the low-rank approximation of matrices that produces accuracy that is very close to the best possible accuracy, for matrices of arbitrary sizes. We illustrate our theoretical results via several numerical examples.
Bibliography:ObjectType-Article-2
SourceType-Scholarly Journals-1
ObjectType-Feature-1
content type line 23
ISSN:0895-4798
1095-7162
DOI:10.1137/080736417