Advanced Search
MyIDEAS: Login

Generalized power method for sparse principal component analysis

Contents:

Author Info

  • Journée, Michel

    (---)

  • Nesterov, Yurii

    (Université catholique de Louvain (UCL). Center for Operations Research and Econometrics (CORE))

  • Richtarik, Peter

    (Université catholique de Louvain (UCL). Center for Operations Research and Econometrics (CORE))

  • Sepulchre, Rodolphe
Registered author(s):

    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.

    Download Info

    If you experience problems downloading a file, check if you have the proper application to view it first. In case of further problems read the IDEAS help page. Note that these files are not on the IDEAS site. Please be patient as the files may be large.
    File URL: http://www.uclouvain.be/cps/ucl/doc/core/documents/coredp2008_70.pdf
    Download Restriction: no

    Bibliographic Info

    Paper provided by Université catholique de Louvain, Center for Operations Research and Econometrics (CORE) in its series CORE Discussion Papers with number 2008070.

    as in new window
    Length:
    Date of creation: 01 Nov 2008
    Date of revision:
    Handle: RePEc:cor:louvco:2008070

    Contact details of provider:
    Postal: Voie du Roman Pays 34, 1348 Louvain-la-Neuve (Belgium)
    Phone: 32(10)474321
    Fax: +32 10474304
    Email:
    Web page: http://www.uclouvain.be/core
    More information through EDIRC

    Related research

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

    This paper has been announced in the following NEP Reports:

    References

    No references listed on IDEAS
    You can help add them by filling out this form.

    Citations

    Lists

    This item is not listed on Wikipedia, on a reading list or among the top items on IDEAS.

    Statistics

    Access and download statistics

    Corrections

    When requesting a correction, please mention this item's handle: RePEc:cor:louvco:2008070. See general information about how to correct material in RePEc.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Alain GILLIS).

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If references are entirely missing, you can add them using this form.

    If the full references list an item that is present in RePEc, but the system did not link to it, you can help with this form.

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your profile, as there may be some citations waiting for confirmation.

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.