Scalable Bayesian Structure Learning for Gaussian Graphical Models Using Marginal Pseudo-likelihood

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Mohammadi, Reza, Schoonhoven, Marit, Vogels, Lucas, Birbil, S. Ilker
Natura: Preprint
Pubblicazione: 2023
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866916921494994944
author Mohammadi, Reza
Schoonhoven, Marit
Vogels, Lucas
Birbil, S. Ilker
author_facet Mohammadi, Reza
Schoonhoven, Marit
Vogels, Lucas
Birbil, S. Ilker
contents Bayesian methods for learning Gaussian graphical models offer a principled framework for quantifying model uncertainty and incorporating prior knowledge. However, their scalability is constrained by the computational cost of jointly exploring graph structures and precision matrices. To address this challenge, we perform inference directly on the graph by integrating out the precision matrix. We adopt a marginal pseudo-likelihood approach, eliminating the need to compute intractable normalizing constants and perform computationally intensive precision matrix sampling. Building on this framework, we develop continuous-time (birth-death) and discrete-time (reversible jump) Markov chain Monte Carlo (MCMC) algorithms that efficiently explore the posterior over graph space. We establish theoretical guarantees for posterior contraction, convergence, and graph selection consistency. The algorithms scale to large graph spaces, enabling parallel exploration for graphs with over 1,000 nodes, while providing uncertainty quantification and supporting flexible prior specification over the graph space. Extensive simulations show substantial computational gains over state-of-the-art Bayesian approaches without sacrificing graph recovery accuracy. Applications to human and mouse gene expression datasets demonstrate the ability of our approach to recover biologically meaningful structures and quantify uncertainty in complex networks. An implementation is available in the R package BDgraph.
format Preprint
id arxiv_https___arxiv_org_abs_2307_00127
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Scalable Bayesian Structure Learning for Gaussian Graphical Models Using Marginal Pseudo-likelihood
Mohammadi, Reza
Schoonhoven, Marit
Vogels, Lucas
Birbil, S. Ilker
Methodology
Computation
Machine Learning
Bayesian methods for learning Gaussian graphical models offer a principled framework for quantifying model uncertainty and incorporating prior knowledge. However, their scalability is constrained by the computational cost of jointly exploring graph structures and precision matrices. To address this challenge, we perform inference directly on the graph by integrating out the precision matrix. We adopt a marginal pseudo-likelihood approach, eliminating the need to compute intractable normalizing constants and perform computationally intensive precision matrix sampling. Building on this framework, we develop continuous-time (birth-death) and discrete-time (reversible jump) Markov chain Monte Carlo (MCMC) algorithms that efficiently explore the posterior over graph space. We establish theoretical guarantees for posterior contraction, convergence, and graph selection consistency. The algorithms scale to large graph spaces, enabling parallel exploration for graphs with over 1,000 nodes, while providing uncertainty quantification and supporting flexible prior specification over the graph space. Extensive simulations show substantial computational gains over state-of-the-art Bayesian approaches without sacrificing graph recovery accuracy. Applications to human and mouse gene expression datasets demonstrate the ability of our approach to recover biologically meaningful structures and quantify uncertainty in complex networks. An implementation is available in the R package BDgraph.
title Scalable Bayesian Structure Learning for Gaussian Graphical Models Using Marginal Pseudo-likelihood
topic Methodology
Computation
Machine Learning
url https://arxiv.org/abs/2307.00127