Advanced Search
MyIDEAS: Login to save this article or follow this journal

A note on the lack of symmetry in the graphical lasso

Contents:

Author Info

  • Rolfs, Benjamin T.
  • Rajaratnam, Bala
Registered author(s):

    Abstract

    The graphical lasso (glasso) is a widely-used fast algorithm for estimating sparse inverse covariance matrices. The glasso solves an ℓ1 penalized maximum likelihood problem and is available as an R library on CRAN. The output from the glasso, a regularized covariance matrix estimate Σˆglasso and a sparse inverse covariance matrix estimate Ωˆglasso, not only identify a graphical model but can also serve as intermediate inputs into multivariate procedures such as PCA, LDA, MANOVA, and others. The glasso indeed produces a covariance matrix estimate Σˆglasso which solves the ℓ1 penalized optimization problem in a dual sense; however, the method for producing Ωˆglasso after this optimization is inexact and may produce asymmetric estimates. This problem is exacerbated when the amount of ℓ1 regularization that is applied is small, which in turn is more likely to occur if the true underlying inverse covariance matrix is not sparse. The lack of symmetry can potentially have consequences. First, it implies that Σˆglasso−1≠Ωˆglasso and, second, asymmetry can possibly lead to negative or complex eigenvalues, rendering many multivariate procedures which may depend on Ωˆglasso unusable. We demonstrate this problem, explain its causes, and propose possible remedies.

    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.sciencedirect.com/science/article/pii/S0167947312002873
    Download Restriction: Full text for ScienceDirect subscribers only.

    As the access to this document is restricted, you may want to look for a different version under "Related research" (further below) or search for a different version of it.

    Bibliographic Info

    Article provided by Elsevier in its journal Computational Statistics & Data Analysis.

    Volume (Year): 57 (2013)
    Issue (Month): 1 ()
    Pages: 429-434

    as in new window
    Handle: RePEc:eee:csdana:v:57:y:2013:i:1:p:429-434

    Contact details of provider:
    Web page: http://www.elsevier.com/locate/csda

    Related research

    Keywords: Concentration model selection; Glasso; Graphical Gaussian models; Graphical lasso; ℓ1 regularization;

    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:eee:csdana:v:57:y:2013:i:1:p:429-434. 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: (Zhang, Lei).

    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.