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.

The Mixture of Markov Jump Processes: Monte Carlo Method and the EM Estimation

Statistical inference for a mixture of Markov jump processes

Mixing it up: A general framework for Markovian statistics

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

An EM algorithm for continuous-time bivariate Markov chains

Entropic Matching for Expectation Propagation of Markov Jump Processes

Pathwise and distributional approximations of semi-Markov processes

A Multi-Index Markov Chain Monte Carlo Method

Strongly Convergent Homogeneous Approximations to Inhomogeneous Markov Jump Processes and Applications

Bayesian Static Parameter Estimation for Partially Observed Diffusions Via Multilevel Monte Carlo.

Parameter Estimation and Model Testing for Markov Processes Via Conditional Characteristic Functions

Markov Chain Monte Carlo for Continuous-Time Switching Dynamical Systems

Parameter Estimation for Partially Observed McKean-Vlasov Diffusions

Estimating the Mixing Coefficients of Geometrically Ergodic Markov Processes

Dynamical mixture modeling with fast, automatic determination of Markov chains

Mixed memory markov models: Decomposing complex stochastic processes as mixtures of simpler ones

Interacting Markov chain Monte Carlo methods for solving nonlinear measure-valued equations

Likelihood Estimation of Jump-Diffusions. Extensions from Diffusions to Jump-Diffusions, Implementation with Automatic Differentiation, and Applications

Empirical and Instance-Dependent Estimation of Markov Chain and Mixing Time

Neural Markov Jump Processes