Gradient-Based Markov Chain Monte Carlo for Bayesian Inference With Non-differentiable Priors

From MaRDI portal
Revision as of 06:18, 10 July 2024 by Import240710060729 (talk | contribs) (Created automatically from import240710060729)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Publication:6110732

DOI10.1080/01621459.2021.1909600zbMATH Open1524.62473arXiv2103.09017OpenAlexW3152141212MaRDI QIDQ6110732

Author name not available (Why is that?)

Publication date: 6 July 2023

Published in: (Search for Journal in Brave)

Abstract: The use of non-differentiable priors in Bayesian statistics has become increasingly popular, in particular in Bayesian imaging analysis. Current state of the art methods are approximate in the sense that they replace the posterior with a smooth approximation via Moreau-Yosida envelopes, and apply gradient-based discretized diffusions to sample from the resulting distribution. We characterize the error of the Moreau-Yosida approximation and propose a novel implementation using underdamped Langevin dynamics. In misson-critical cases, however, replacing the posterior with an approximation may not be a viable option. Instead, we show that Piecewise-Deterministic Markov Processes (PDMP) can be utilized for exact posterior inference from distributions satisfying almost everywhere differentiability. Furthermore, in contrast with diffusion-based methods, the suggested PDMP-based samplers place no assumptions on the prior shape, nor require access to a computationally cheap proximal operator, and consequently have a much broader scope of application. Through detailed numerical examples, including a non-differentiable circular distribution and a non-convex genomics model, we elucidate the relative strengths of these sampling methods on problems of moderate to high dimensions, underlining the benefits of PDMP-based methods when accurate sampling is decisive.


Full work available at URL: https://arxiv.org/abs/2103.09017



No records found.


No records found.








This page was built for publication: Gradient-Based Markov Chain Monte Carlo for Bayesian Inference With Non-differentiable Priors

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6110732)