In this paper, we consider high-dimensional Gaussian graphical models where the true underlying graph is decomposable. A hierarchical G-Wishart prior is proposed to conduct a Bayesian inference for the precision matrix and its graph structure. Although the posterior asymptotics using the G-Wishart prior has received increasing attention in recent years, most of the results assume moderate high-dimensional settings, where the number of variables p is smaller than the sample size n. However, this assumption might not hold in many real applications such as genomics, speech recognition and climatology. Motivated by this gap, we investigate asymptotic properties of posteriors under the high-dimensional setting where p can be much larger than n. The pairwise Bayes factor consistency, posterior ratio consistency and graph selection consistency are obtained in this high-dimensional setting. Furthermore, the posterior convergence rate for precision matrices under the matrix -norm is derived, which is faster than posterior convergence rates obtained in existing literature. A simulation study confirms that the proposed Bayesian procedure outperforms competitors.
KL was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2020R1A4A1018207) and INHA UNIVERSITY Research Grant. XC was supported by Simons Foundation’s collaboration grant (No.635213) and Taft Summer Research Fellowship at University of Cincinnati.
We would like to thank the Editor, the Associate Editor and two reviewers for their insightful comments which have led to improvements of an earlier version of this paper.
"Bayesian inference for high-dimensional decomposable graphs." Electron. J. Statist. 15 (1) 1549 - 1582, 2021. https://doi.org/10.1214/21-EJS1822