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.

Novel sampling method for the von Mises–Fisher distribution

A fast algorithm for sampling from the posterior of a von Mises distribution

Penalized Maximum Likelihood Estimator for Mixture of von Mises-Fisher Distributions

An importance sampling method for Feldman-Cousins confidence intervals

IID Sampling from Intractable Distributions

A Novel Quota Sampling Algorithm for Generating Representative Random Samples given Small Sample Size

Sampling Approximately Low-Rank Ising Models: MCMC meets Variational Methods

Rejection sampling for generating random numbers from weighted distributions

The Coupled Rejection Sampler

Importance sampling method for Feldman-Cousins confidence intervals

New algorithms for sampling and diffusion models

Ensemble Rejection Sampling

Efficient Importance Sampling for Rare Event Simulation with Applications

Schr{\"o}dinger-F{\"o}llmer Sampler: Sampling without Ergodicity

Sampling-Based Approaches to Calculating Marginal Densities

Markov Chain Sampling Methods for Dirichlet Process Mixture Models

Cross-entropy-based directional importance sampling with von Mises-Fisher mixture model for reliability analysis

Schr{ö}dinger-F{ö}llmer Sampler: Sampling without Ergodicity

In search of rogue waves: a novel proposal distribution for parallelized rejection sampling of the truncated KdV Gibbs measure

The distribution of order statistics under sampling without replacement

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