The Coordinate Sampler: A Non-Reversible Gibbs-like MCMC Sampler

09/10/2018
by   Changye Wu, et al.
0

In this article, we derive a novel non-reversible, continuous-time Markov chain Monte Carlo (MCMC) sampler, called Coordinate Sampler, based on a piecewise deterministic Markov process (PDMP), which can be seen as a variant of the Zigzag sampler. In addition to proving a theoretical validation for this new sampling algorithm, we show that the Markov chain it induces exhibits geometrical ergodicity convergence, for distributions whose tails decay at least as fast as an exponential distribution and at most as fast as a Gaussian distribution. Several numerical examples highlight that our coordinate sampler is more efficient than the Zigzag sampler, in terms of effective sample size.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/12/2018

Bouncy Hybrid Sampler as a Unifying Device

This work introduces a class of rejection-free Markov chain Monte Carlo ...
research
01/16/2020

A piecewise deterministic Monte Carlo method for diffusion bridges

We introduce the use of the Zig-Zag sampler to the problem of sampling o...
research
08/26/2019

MALA-within-Gibbs samplers for high-dimensional distributions with sparse conditional structure

Markov chain Monte Carlo (MCMC) samplers are numerical methods for drawi...
research
10/22/2020

Reversible Jump PDMP Samplers for Variable Selection

A new class of Markov chain Monte Carlo (MCMC) algorithms, based on simu...
research
02/20/2023

Gibbs Sampler for Matrix Generalized Inverse Gaussian Distributions

Sampling from matrix generalized inverse Gaussian (MGIG) distributions i...
research
09/28/2021

Reversible Gromov-Monge Sampler for Simulation-Based Inference

This paper introduces a new simulation-based inference procedure to mode...
research
05/30/2017

Zonotope hit-and-run for efficient sampling from projection DPPs

Determinantal point processes (DPPs) are distributions over sets of item...

Please sign up or login with your details

Forgot password? Click here to reset