Bayesian Covariance Matrix Estimation using a Mixture of Decomposable Graphical Models

Helen Armstrong (), Christopher Carter (), Kevin Wong and Robert Kohn ()
Additional contact information
Helen Armstrong: School of Mathematics, University of New South Wales
Christopher Carter: School of Economics, University of New South Wales
Kevin Wong: Graduate University for Advanced Studies, Tokyo, Japan

No 2007-13, Discussion Papers from School of Economics, The University of New South Wales

Abstract: Estimating a covariance matrix efficiently and discovering its structure are important statistical problems with applications in many fields. This article takes a Bayesian approach to estimate the covariance matrix of Gaussian data. We use ideas from Gaussian graphical models and model selection to construct a prior for the covariance matrix that is a mixture over all decomposable graphs, where a graph means the configuration of nonzero offdiagonal elements in the inverse of the covariance matrix. Our prior for the covariance matrix is such that the probability of each graph size is specified by the user and graphs of equal size are assigned equal probability. Most previous approaches assume that all graphs are equally probable. We give empirical results that show the prior that assigns equal probability over graph sizes outperforms the prior that assigns equal probability over all graphs, both in identifying the correct decomposable graph and in more efficiently estimating the covariance matrix. The advantage is greatest when the number of observations is small relative to the dimension of the covariance matrix. The article also shows empirically that there is minimal change in statistical efficiency in using the mixture over decomposable graphs prior for estimating a general covariance compared to the Bayesian estimator by Wong et al. (2003), even when the graph of the covariance matrix is nondecomposable. However, our approach has some important advantages over that of Wong et al. (2003). Our method requires the number of decomposable graphs for each graph size. We show how to estimate these numbers using simulation and that the simulation results agree with analytic results when such results are known. We also show how to estimate the posterior distribution of the covariance matrix using Markov chain Monte Carlo with the elements of the covariance matrix integrated out and give empirical results that show the sampler is computationally efficient and converges rapidly. Finally, we note that both the prior and the simulation method to evaluate the prior apply generally to any decomposable graphical model.

Keywords: Covariance selection; Graphical models; Reduced conditional sampling; Variable selection (search for similar items in EconPapers)
Pages: 27 pages
Date: 2007-04
New Economics Papers: this item is included in nep-ecm and nep-ore
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
http://wwwdocs.fce.unsw.edu.au/economics/Research/WorkingPapers/2007_13.pdf (application/pdf)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:swe:wpaper:2007-13

Access Statistics for this paper

More papers in Discussion Papers from School of Economics, The University of New South Wales Contact information at EDIRC.
Bibliographic data for series maintained by Hongyi Li ().