IDEAS home Printed from https://ideas.repec.org/a/eee/csdana/v56y2012i8p2430-2441.html
   My bibliography  Save this article

The cost of using decomposable Gaussian graphical models for computational convenience

Author

Listed:
  • Fitch, A. Marie
  • Jones, Beatrix

Abstract

Graphical models are a powerful tool for describing patterns of conditional independence, and can also be used to regularize the covariance matrix. Vertices in the graph represent variables, and in the Gaussian setting, edges between vertices are equivalent to non-zero elements in the inverse covariance matrix. Models that can be represented as a decomposable (triangulated) graph are more computationally tractable; in fact, in the high-dimensional Bayesian setting it is common to restrict model selection procedures to decomposable models. We consider estimation of the covariance and inverse covariance matrix where the true model forms a cycle, but estimation is performed supposing that the pattern of zeros is a decomposable graphical model, where the elements restricted to zero are a subset of those in the true matrix. The variance of the maximum likelihood estimator based on the decomposable model is demonstrably larger than for the true non-decomposable model, and which decomposable model is selected affects the variance of particular elements of the matrix. When estimating the inverse covariance matrix the cost in terms of accuracy for using the decomposable model is fairly small, even when the difference in sparsity is large and the sample size is fairly small (e.g., the true model is a cycle of size 50, and the sample size is 51). However, when estimating the covariance matrix, the estimators for most elements had a dramatic increase in variance (200-fold in some cases) when a decomposable model was substituted. These increases become more pronounced as the difference in sparsity between models increases.

Suggested Citation

  • Fitch, A. Marie & Jones, Beatrix, 2012. "The cost of using decomposable Gaussian graphical models for computational convenience," Computational Statistics & Data Analysis, Elsevier, vol. 56(8), pages 2430-2441.
  • Handle: RePEc:eee:csdana:v:56:y:2012:i:8:p:2430-2441
    DOI: 10.1016/j.csda.2012.01.020
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0167947312000400
    Download Restriction: Full text for ScienceDirect subscribers only.

    As the access to this document is restricted, you may want to search for a different version of it.

    References listed on IDEAS

    as
    1. Aliye Atay-Kayis & Helène Massam, 2005. "A Monte Carlo method for computing the marginal likelihood in nondecomposable Gaussian graphical models," Biometrika, Biometrika Trust, vol. 92(2), pages 317-335, June.
    2. Alberto Roverato, 2002. "Hyper Inverse Wishart Distribution for Non-decomposable Graphs and its Application to Bayesian Inference for Gaussian Graphical Models," Scandinavian Journal of Statistics, Danish Society for Theoretical Statistics;Finnish Statistical Society;Norwegian Statistical Association;Swedish Statistical Association, vol. 29(3), pages 391-411.
    3. Frederick Wong, 2003. "Efficient estimation of covariance selection models," Biometrika, Biometrika Trust, vol. 90(4), pages 809-830, December.
    4. Mathias Drton & Michael Eichler, 2006. "Maximum Likelihood Estimation in Gaussian Chain Graph Models under the Alternative Markov Property," Scandinavian Journal of Statistics, Danish Society for Theoretical Statistics;Finnish Statistical Society;Norwegian Statistical Association;Swedish Statistical Association, vol. 33(2), pages 247-257.
    Full references (including those not matched with items on IDEAS)

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:csdana:v:56:y:2012:i:8:p:2430-2441. 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: (Dana Niculescu). General contact details of provider: http://www.elsevier.com/locate/csda .

    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 CitEc recognized a reference but did not link an item in RePEc 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 RePEc Author Service 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.

    IDEAS is a RePEc service hosted by the Research Division of the Federal Reserve Bank of St. Louis . RePEc uses bibliographic data supplied by the respective publishers.