Fast Gibbs sampling for high-dimensional Bayesian inversion
From MaRDI portal
Abstract: Solving ill-posed inverse problems by Bayesian inference has recently attracted considerable attention. Compared to deterministic approaches, the probabilistic representation of the solution by the posterior distribution can be exploited to explore and quantify its uncertainties. In applications where the inverse solution is subject to further analysis procedures, this can be a significant advantage. Alongside theoretical progress, various new computational techniques allow to sample very high dimensional posterior distributions: In [Lucka2012], a Markov chain Monte Carlo (MCMC) posterior sampler was developed for linear inverse problems with -type priors. In this article, we extend this single component Gibbs-type sampler to a wide range of priors used in Bayesian inversion, such as general priors with additional hard constraints. Besides a fast computation of the conditional, single component densities in an explicit, parameterized form, a fast, robust and exact sampling from these one-dimensional densities is key to obtain an efficient algorithm. We demonstrate that a generalization of slice sampling can utilize their specific structure for this task and illustrate the performance of the resulting slice-within-Gibbs samplers by different computed examples. These new samplers allow us to perform sample-based Bayesian inference in high-dimensional scenarios with certain priors for the first time, including the inversion of computed tomography (CT) data with the popular isotropic total variation (TV) prior.
Recommendations
- Fast Markov chain Monte Carlo sampling for sparse Bayesian inference in high-dimensional inverse problems using L1-type priors
- Low-rank independence samplers in hierarchical Bayesian inverse problems
- Fast Bayesian inversion for high dimensional inverse problems
- Analysis of the Gibbs Sampler for Hierarchical Inverse Problems
- Bayesian inverse problems with l₁ priors: a randomize-then-optimize approach
Cited in
(23)- Gibbs sampling for parameters estimation and change-points detection in inverse problems. Application to electromagnetic imaging.
- A method for estimating dominant acoustic backscatter mechanism of water-seabed interface via relative entropy estimation
- Proximal nested sampling for high-dimensional Bayesian model selection
- Fast Bayesian inversion for high dimensional inverse problems
- Fast automatic Bayesian cubature using lattice sampling
- Analysis of the Gibbs Sampler for Hierarchical Inverse Problems
- Hybrid samplers for ill-posed inverse problems
- Fast sampling in a linear-Gaussian inverse problem
- Recent advances in inferential solutions to inverse problems
- A hierarchical Bayesian perspective on majorization-minimization for non-convex sparse regression: application to M/EEG source imaging
- Efficient Bayesian computation by proximal Markov chain Monte Carlo: when Langevin meets Moreau
- Low-rank independence samplers in hierarchical Bayesian inverse problems
- Fast Markov chain Monte Carlo sampling for sparse Bayesian inference in high-dimensional inverse problems using L1-type priors
- A proximal Markov chain Monte Carlo method for Bayesian inference in imaging inverse problems: when Langevin meets Moreau
- A hybrid Gibbs sampler for edge-preserving tomographic reconstruction with uncertain view angles
- MCMC Algorithms for Computational UQ of Nonnegativity Constrained Linear Inverse Problems
- Two Metropolis--Hastings Algorithms for Posterior Measures with Non-Gaussian Priors in Infinite Dimensions
- Bayesian inverse problems with l₁ priors: a randomize-then-optimize approach
- A CVAE-within-Gibbs sampler for Bayesian linear inverse problems with hyperparameters
- The Split Gibbs Sampler Revisited: Improvements to Its Algorithmic Structure and Augmented Target Distribution
- Certified coordinate selection for high-dimensional Bayesian inversion with Laplace prior
- Uncertainty quantification for linear inverse problems with Besov prior: a randomize-then-optimize method
- Fused \(L_{1/2}\) prior for large scale linear inverse problem with Gibbs bouncy particle sampler
This page was built for publication: Fast Gibbs sampling for high-dimensional Bayesian inversion
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2835454)
