Optimization Online


Generalized power method for sparse principal component analysis

Michel Journee(m.journee***at***ulg.ac.be)
Yurii Nesterov(nesterov***at***core.ucl.ac.be)
Peter Richtarik(richtarik***at***core.ucl.ac.be)
Rodolphe Sepulchre(r.sepulchre***at***ulg.ac.be)

Abstract: In this paper we develop a new approach to sparse principal component analysis (sparse PCA). We propose two single-unit and two block optimization formulations of the sparse PCA problem, aimed at extracting a single sparse dominant principal component of a data matrix, or more components at once, respectively. While the initial formulations involve nonconvex functions, and are therefore computationally intractable, we rewrite them into the form of an optimization program involving maximization of a convex function on a compact set. The dimension of the search space is decreased enormously if the data matrix has many more columns (variables) than rows. We then propose and analyze a simple gradient method suited for the task. It appears that our algorithm has best convergence properties in the case when either the objective function or the feasible set are strongly convex, which is the case with our single-unit formulations and can be enforced in the block case. Finally, we demonstrate numerically on a set of random and gene expression test problems that our approach outperforms existing algorithms both in quality of the obtained solution and in computational speed.

Keywords: sparse PCA, power method, gradient ascent, strongly convex sets, block algorithms

Category 1: Applications -- Science and Engineering

Category 2: Convex and Nonsmooth Optimization (Nonsmooth Optimization )

Category 3: Nonlinear Optimization


Download: [PDF]

Entry Submitted: 11/28/2008
Entry Accepted: 11/28/2008
Entry Last Modified: 11/28/2008

Modify/Update this entry

  Visitors Authors More about us Links
  Subscribe, Unsubscribe
Digest Archive
Search, Browse the Repository


Coordinator's Board
Classification Scheme
Give us feedback
Optimization Journals, Sites, Societies
Mathematical Programming Society