Abstract:Many Markov chain Monte Carlo techniques currently available rely on discrete-time reversible Markov processes whose transition kernels are variations of the Metropolis–Hastings algorithm. We explore and generalize an alternative scheme recently introduced in the physics literature (Peters and de With 2012) where the target distribution is explored using a continuous-time nonreversible piecewise-deterministic Markov process. In the Metropolis–Hastings algorithm, a trial move to a region of lower target density, equivalently of higher “energy,” than the current state can be rejected with positive probability. In this alternative approach, a particle moves along straight lines around the space and, when facing a high energy barrier, it is not rejected but its path is modified by bouncing against this barrier. By reformulating this algorithm using inhomogeneous Poisson processes, we exploit standard sampling techniques to simulate exactly this Markov process in a wide range of scenarios of interest. Additionally, when the target distribution is given by a product of factors dependent only on subsets of the state variables, such as the posterior distribution associated with a probabilistic graphical model, this method can be modified to take advantage of this structure by allowing computationally cheaper “local” bounces, which only involve the state variables associated with a factor, while the other state variables keep on evolving. In this context, by leveraging techniques from chemical kinetics, we propose several computationally efficient implementations. Experimentally, this new class of Markov chain Monte Carlo schemes compares favorably to state-of-the-art methods on various Bayesian inference tasks, including for high-dimensional models and large datasets. Supplementary materials for this article are available online.

On sequential Monte Carlo, partial rejection control and approximate Bayesian computation

An invitation to sequential Monte Carlo samplers

Multicanonical sequential Monte Carlo sampler for uncertainty quantification

Stratified sampling and bootstrapping for approximate Bayesian computation

Sequential Monte Carlo for cut-Bayesian posterior computation

Reducing rejection exponentially improves Markov chain Monte Carlo sampling

The Bouncy Particle Sampler: A Nonreversible Rejection-Free Markov Chain Monte Carlo Method

On sequential Monte Carlo sampling methods for Bayesian filtering

Sequential monte carlo samplers

Improving approximate Bayesian computation via quasi-Monte Carlo

Sequential Monte Carlo with Highly Informative Observations

Convergence of sequential Markov Chain Monte Carlo methods: I. Nonlinear flow of probability measures

Using deterministic approximations to accelerate SMC for posterior sampling

Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems

Monte Carlo Approximation of Bayes Factors via Mixing with Surrogate Distributions

Approximating multivariate posterior distribution functions from Monte Carlo samples for sequential Bayesian inference

A Sequential Monte Carlo Algorithm to Incorporate Model Uncertainty in Bayesian Sequential Design

Simple rejection Monte Carlo algorithm and its application to multivariate statistical inference

Sequential Monte Carlo samplers: error bounds and insensitivity to initial conditions

Stochastic Gradient Piecewise Deterministic Monte Carlo Samplers

No Free Lunch for Approximate MCMC