A piecewise deterministic Monte Carlo method for diffusion bridges

Abstract

We introduce the use of the Zig-Zag sampler to the problem of sampling conditional diffusion processes (diffusion bridges). The Zig-Zag sampler is a rejection-free sampling scheme based on a non-reversible continuous piecewise deterministic Markov process. Similar to the Lévy-Ciesielski construction of a Brownian motion, we expand the diffusion path in a truncated Faber-Schauder basis. The coefficients within the basis are sampled using a Zig-Zag sampler. A key innovation is the use of the fully local Algorithm for the Zig-Zag sampler that allows to exploit the sparsity structure implied by the dependency graph of the coefficients and by the subsampling technique to reduce the complexity of the algorithm. We illustrate the performance of the proposed methods in a number of examples.

arXiv comment31 pages, 11 figures. Implementation at https://github.com/SebaGraz/ZZDiffusionBridge/

References in corpus (12)

• A Conceptual Introduction to Hamiltonian Monte Carlo
• The Zig-Zag Process and Super-Efficient Sampling for Bayesian Analysis of Big Data
• Rejection-free Monte-Carlo sampling for general potentials
• Piecewise Deterministic Markov Processes for Continuous-Time Monte Carlo
• Simple simulation of diffusion bridges with application to likelihood inference for diffusions
• Bayesian estimation of discretely observed multi-dimensional diffusion processes using guided proposals
• All-atom computations with irreversible Markov chains
• Continuous-discrete smoothing of diffusions
• Clock Monte Carlo methods
• Simulation of elliptic and hypo-elliptic conditional diffusions
• Adaptive nonparametric drift estimation for diffusion processes using Faber-Schauder expansions
• Simulating bridges using confluent diffusions