Skip to content
DraconiticResearch
Search Sources API
197,693 studies31 sourcesUpdated 4 min ago
PreprintCC BY 4.0arXiv

Metropolized backward sampling for conditional sequential Monte Carlo

Abstract

Particle Gibbs samplers with backward sampling are widely used for sampling the smoothing distribution of state-space models, but, with few particles, the new trajectory often coincides with the current one at many time steps. A known but rarely used remedy replaces the random selection of the final particle by a Metropolis--Hastings move which avoids the current particle. We show that, with backward sampling, this remedy only helps at the last time steps, because its effect on the update rates decays geometrically with the distance to the final time. We propose instead Metropolized backward sampling, which makes such a move at every step of the backward pass, forcibly switching trajectories whenever possible. Additionally, we introduce a coupled version which prevents proposing the current trajectory, maximizing the chances of an actual switch versus the current trajectory and performing perfect sampling when the state-space model has no observations, a feature not shared by backward sampling. We prove that, at every time step, both methods dominate backward sampling as they update the state more often and have smaller lag-one autocorrelations than it does. Our methods are drop-in replacements for backward sampling, and can be used to improve any particle Gibbs sampler.

Related

  1. Parameter estimation in Conditional Sequential Monte Carlo algorithms through Particle Learning
    Alfonso Diz-Lois Palomares, Geir Storvik
    arXivPreprintCC BY-NC-ND 4.0Full textStatisticsarXiv
  2. Sequential Monte Carlo Sampling of Feasible Regions Represented by R-Functions
    Sinan Yıldırım, Sergei Kucherenko
    arXivPreprintarXiv licenceFull textNumerical analysis & scientific computingStatisticsarXiv
  3. Conditional Reflected Backward Stochastic Differential Equations
    Hanwu Li
    arXivPreprintarXiv licenceFull textProbabilityarXiv
  4. Hamiltonian Monte Carlo vs. event-chain Monte Carlo: An appraisal of sampling strategies beyond the diffusive regime
    Werner Krauth
    Journal of Computational PhysicsPeer-reviewedLicence not stated10 citedStatisticsOpenAlex
  5. Differentiable Systematic Resampling for Variational Sequential Monte Carlo
    Fredrik Cumlin, Saikat Chatterjee
    arXivPreprintCC BY 4.0Full textStatisticsarXiv
BibTeX
@misc{corenflos2026metropolized,
  title = {{Metropolized backward sampling for conditional sequential Monte Carlo}},
  author = {Adrien Corenflos},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2610.38753},
  eprint = {2610.38753},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2610.38753},
  note = {CC BY 4.0},
}
RIS
TY  - UNPB
TI  - Metropolized backward sampling for conditional sequential Monte Carlo
AU  - Adrien Corenflos
PY  - 2026
DA  - 2026/10/28
PB  - arXiv
DO  - 10.48550/arxiv.2610.38753
UR  - https://arxiv.org/abs/2610.38753
AB  - Particle Gibbs samplers with backward sampling are widely used for sampling the smoothing distribution of state-space models, but, with few particles, the new trajectory often coincides with the current one at many time steps. A known but rarely used remedy replaces the random selection of the final particle by a Metropolis--Hastings move which avoids the current particle. We show that, with backward sampling, this remedy only helps at the last time steps, because its effect on the update rates decays geometrically with the distance to the final time. We propose instead Metropolized backward sampling, which makes such a move at every step of the backward pass, forcibly switching trajectories whenever possible. Additionally, we introduce a coupled version which prevents proposing the current trajectory, maximizing the chances of an actual switch versus the current trajectory and performing perfect sampling when the state-space model has no observations, a feature not shared by backward sampling. We prove that, at every time step, both methods dominate backward sampling as they update the state more often and have smaller lag-one autocorrelations than it does. Our methods are drop-in replacements for backward sampling, and can be used to improve any particle Gibbs sampler.
ER  -