Proximal nested sampling for high-dimensional Bayesian model selection
From MaRDI portal
Abstract: Bayesian model selection provides a powerful framework for objectively comparing models directly from observed data, without reference to ground truth data. However, Bayesian model selection requires the computation of the marginal likelihood (model evidence), which is computationally challenging, prohibiting its use in many high-dimensional Bayesian inverse problems. With Bayesian imaging applications in mind, in this work we present the proximal nested sampling methodology to objectively compare alternative Bayesian imaging models for applications that use images to inform decisions under uncertainty. The methodology is based on nested sampling, a Monte Carlo approach specialised for model comparison, and exploits proximal Markov chain Monte Carlo techniques to scale efficiently to large problems and to tackle models that are log-concave and not necessarily smooth (e.g., involving l_1 or total-variation priors). The proposed approach can be applied computationally to problems of dimension O(10^6) and beyond, making it suitable for high-dimensional inverse imaging problems. It is validated on large Gaussian models, for which the likelihood is available analytically, and subsequently illustrated on a range of imaging problems where it is used to analyse different choices of dictionary and measurement model.
Recommendations
- Efficient Bayesian computation by proximal Markov chain Monte Carlo: when Langevin meets Moreau
- A proximal Markov chain Monte Carlo method for Bayesian inference in imaging inverse problems: when Langevin meets Moreau
- Scalable Bayesian uncertainty quantification in imaging inverse problems via convex optimization
- A Bayesian model selection method with applications
- Editorial: Bayesian model selection: some thoughts on future directions
Cites work
- scientific article; zbMATH DE number 509150 (Why is no real title available?)
- scientific article; zbMATH DE number 2117879 (Why is no real title available?)
- scientific article; zbMATH DE number 5043156 (Why is no real title available?)
- Accelerating Proximal Markov Chain Monte Carlo by Using an Explicit Stabilized Method
- Accurate Approximations for Posterior Moments and Marginal Densities
- Bayesian computation: a summary of the current state, and samples backwards and forwards
- Bayesian inference and uncertainty quantification for medical image reconstruction with Poisson data
- Convex analysis and monotone operator theory in Hilbert spaces
- Diffusive nested sampling
- Efficient Bayesian computation by proximal Markov chain Monte Carlo: when Langevin meets Moreau
- Fast Gibbs sampling for high-dimensional Bayesian inversion
- Layered adaptive importance sampling
- Marginal Likelihood From the Metropolis–Hastings Output
- Marginal Likelihood from the Gibbs Output
- Maximum likelihood estimation of regularization parameters in high-dimensional inverse problems: an empirical Bayesian approach. I: Methodology and experiments
- Nested sampling for general Bayesian computation
- Numerical Bayesian Methods Applied to Signal Processing
- On learned operator correction in inverse problems
- Properties of nested sampling
- Proximal Markov chain Monte Carlo algorithms
- Proximal splitting methods in signal processing
- Reversible jump Markov chain Monte Carlo computation and Bayesian model determination
- Statistical and computational inverse problems.
- The Bayesian Choice
Cited in
(9)- A proximal Markov chain Monte Carlo method for Bayesian inference in imaging inverse problems: when Langevin meets Moreau
- Accelerated Bayesian imaging by relaxed proximal-point Langevin sampling
- Marginal likelihood estimation in semiblind image deconvolution: a stochastic approximation approach
- An optimal Bayesian strategy for comparing Wiener-hunt deconvolution models in the absence of ground truth
- Subgradient Langevin methods for sampling from nonsmooth potentials
- The Split Gibbs Sampler Revisited: Improvements to Its Algorithmic Structure and Augmented Target Distribution
- Efficient Bayesian computation by proximal Markov chain Monte Carlo: when Langevin meets Moreau
- Smoothing unadjusted Langevin algorithms for nonsmooth composite potential functions
- The method of forced probabilities: a computation trick for Bayesian model evidence
This page was built for publication: Proximal nested sampling for high-dimensional Bayesian model selection
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2084321)
