Fast and numerically stable particle-based online additive smoothing: the AdaSmooth algorithm
Abstract
We present a novel sequential Monte Carlo approach to online smoothing of additive functionals in a very general class of path-space models. Hitherto, the solutions proposed in the literature suffer from either long-term numerical instability due to particle-path degeneracy or, in the case that degeneracy is remedied by particle approximation of the so-called backward kernel, high computational demands. In order to balance optimally computational speed against numerical stability, we propose to furnish a (fast) naive particle smoother, propagating recursively a sample of particles and associated smoothing statistics, with an adaptive backward-sampling-based updating rule which allows the number of (costly) backward samples to be kept at a minimum. This yields a new, function-specific additive smoothing algorithm, AdaSmooth, which is computationally fast, numerically stable and easy to implement. The algorithm is provided with rigorous theoretical results guaranteeing its consistency, asymptotic normality and long-term stability as well as numerical results demonstrating empirically the clear superiority of AdaSmooth to existing algorithms.
Cite
@article{arxiv.2108.00432,
title = {Fast and numerically stable particle-based online additive smoothing: the AdaSmooth algorithm},
author = {Alessandro Mastrototaro and Jimmy Olsson and Johan Alenlöv},
journal= {arXiv preprint arXiv:2108.00432},
year = {2022}
}
Comments
67 pages, 5 figures, 2 tables. Added initial note: "This is an original manuscript of an article published by Taylor & Francis in the Journal of the American Statistical Association (JASA) on 10 October 2022, available online: https://www.tandfonline.com/doi/full/10.1080/01621459.2022.2118602"