English

Conditional sequential Monte Carlo in high dimensions

Computation 2021-08-24 v1

Abstract

The iterated conditional sequential Monte Carlo (i-CSMC) algorithm from Andrieu, Doucet and Holenstein (2010) is an MCMC approach for efficiently sampling from the joint posterior distribution of the TT latent states in challenging time-series models, e.g. in non-linear or non-Gaussian state-space models. It is also the main ingredient in particle Gibbs samplers which infer unknown model parameters alongside the latent states. In this work, we first prove that the i-CSMC algorithm suffers from a curse of dimension in the dimension of the states, DD: it breaks down unless the number of samples ("particles"), NN, proposed by the algorithm grows exponentially with DD. Then, we present a novel "local" version of the algorithm which proposes particles using Gaussian random-walk moves that are suitably scaled with DD. We prove that this iterated random-walk conditional sequential Monte Carlo (i-RW-CSMC) algorithm avoids the curse of dimension: for arbitrary NN, its acceptance rates and expected squared jumping distance converge to non-trivial limits as DD \to \infty. If T=N=1T = N = 1, our proposed algorithm reduces to a Metropolis--Hastings or Barker's algorithm with Gaussian random-walk moves and we recover the well known scaling limits for such algorithms.

Keywords

Cite

@article{arxiv.2108.10277,
  title  = {Conditional sequential Monte Carlo in high dimensions},
  author = {Axel Finke and Alexandre H. Thiery},
  journal= {arXiv preprint arXiv:2108.10277},
  year   = {2021}
}

Comments

47 pages, 5 figures

R2 v1 2026-06-24T05:21:13.262Z