On sequential Monte Carlo, partial rejection control and approximate Bayesian computation
Computation
2009-11-11 v2 Statistics Theory
Statistics Theory
Abstract
We present a sequential Monte Carlo sampler variant of the partial rejection control algorithm, and show that this variant can be considered as a sequential Monte Carlo sampler with a modified mutation kernel. We prove that the new sampler can reduce the variance of the incremental importance weights when compared with standard sequential Monte Carlo samplers. We provide a study of theoretical properties of the new algorithm, and make connections with some existing algorithms. Finally, the sampler is adapted for application under the challenging "likelihood free," approximate Bayesian computation modelling framework, where we demonstrate superior performance over existing likelihood-free samplers.
Cite
@article{arxiv.0808.3466,
title = {On sequential Monte Carlo, partial rejection control and approximate Bayesian computation},
author = {G. W. Peters and Y. Fan and S. A. Sisson},
journal= {arXiv preprint arXiv:0808.3466},
year = {2009}
}