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.

Twisting the Alive Particle Filter

The Alive Particle Filter and Its Use in Particle Markov Chain Monte Carlo

New Real-Time Resampling Algorithm For Particle Filters

Monte Carlo twisting for particle filters

Persistent Sampling: Unleashing the Potential of Sequential Monte Carlo

A General-Purpose Fixed-Lag No-U-Turn Sampler for Nonlinear Non-Gaussian State Space Models

The Application of Zig-Zag Sampler in Sequential Markov Chain Monte Carlo

Reversed particle filtering for hidden markov models

A Vine-Copula Based Adaptive Mcmc Sampler for Efficient Inference of Dynamical Systems

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

Particle Gibbs for Likelihood-Free Inference of State Space Models with Application to Stochastic Volatility

Sequential Monte Carlo with Highly Informative Observations

Numerically stable online estimation of variance in particle filters

Guidance for twisted particle filter: a continuous-time perspective

Ecological non-linear state space model selection via adaptive particle Markov chain Monte Carlo (AdPMCMC)

Adaptive sequential Monte Carlo by means of mixture of experts

Auxiliary MCMC and particle Gibbs samplers for parallelisable inference in latent dynamical systems

On Disturbance State-Space Models and the Particle Marginal Metropolis-Hastings Sampler

Augmentation Schemes for Particle MCMC

Bayesian Conditional Monte Carlo Algorithms for Sequential Single and Multi-Object filtering

A Novel SMC-PHD Filter Based on Particle Compensation