Stability of Conditional Sequential Monte Carlo

06/18/2018
by   Bernd Kuhlenschmidt, et al.
0

The particle Gibbs (PG) sampler is a Markov Chain Monte Carlo (MCMC) algorithm, which uses an interacting particle system to perform the Gibbs steps. Each Gibbs step consists of simulating a particle system conditioned on one particle path. It relies on a conditional Sequential Monte Carlo (cSMC) method to create the particle system. We propose a novel interpretation of the cSMC algorithm as a perturbed Sequential Monte Carlo (SMC) method and apply telescopic decompositions developed for the analysis of SMC algorithms delmoral2004 to derive a bound for the distance between the expected sampled path from cSMC and the target distribution of the MCMC algorithm. This can be used to get a uniform ergodicity result. In particular, we can show that the mixing rate of cSMC can be kept constant by increasing the number of particles linearly with the number of observations. Based on our decomposition, we also prove a central limit theorem for the cSMC Algorithm, which cannot be done using the approaches in Andrieu2013 and Lindsten2014.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/14/2020

Particle Gibbs Sampling for Bayesian Phylogenetic inference

The combinatorial sequential Monte Carlo (CSMC) has been demonstrated to...
research
02/04/2022

De-Sequentialized Monte Carlo: a parallel-in-time particle smoother

Particle smoothers are SMC (Sequential Monte Carlo) algorithms designed ...
research
11/04/2020

Waste-free Sequential Monte Carlo

A standard way to move particles in a SMC sampler is to apply several st...
research
08/30/2023

Sequential Bayesian Predictive Synthesis

Dynamic Bayesian predictive synthesis is a formal approach to coherently...
research
04/21/2019

Particle filter efficiency under limited communication

Sequential Monte Carlo (SMC) methods are typically not straightforward t...
research
09/21/2022

Particle-based, rapid incremental smoother meets particle Gibbs

The particle-based, rapid incremental smoother (PARIS) is a sequential M...
research
05/31/2018

Bayesian inference in decomposable graphical models using sequential Monte Carlo methods

InthisstudywepresentasequentialsamplingmethodologyforBayesian inference ...

Please sign up or login with your details

Forgot password? Click here to reset