ON RESAMPLING SCHEMES FOR PARTICLE FILTERS WITH WEAKLY INFORMATIVE OBSERVATIONS

被引:4
作者
Chopin, Nicolas [1 ]
Singh, Sumeetpal S. [2 ]
Soto, Tomas [3 ]
Vihola, Matti [4 ]
机构
[1] Inst Polytech Paris, ENSAE, Palaiseau, France
[2] Univ Cambridge, Dept Engn, Cambridge, England
[3] LUT Univ, Sch Engn Sci, Lappeenranta, Finland
[4] Univ Jyvaskyla, Dept Math & Stat, Jyvaskyla, Finland
基金
芬兰科学院;
关键词
Feynman-Kac model; hidden Markov model; particle filter; path integral; resampling;
D O I
10.1214/22-AOS2222
中图分类号
O21 [概率论与数理统计]; C8 [统计学];
学科分类号
020208 ; 070103 ; 0714 ;
摘要
We consider particle filters with weakly informative observations (or 'po-tentials') relative to the latent state dynamics. The particular focus of this work is on particle filters to approximate time-discretisations of continuous -time Feynman-Kac path integral models-a scenario that naturally arises when addressing filtering and smoothing problems in continuous time-but our findings are indicative about weakly informative settings beyond this con-text too. We study the performance of different resampling schemes, such as systematic resampling, SSP (Srinivasan sampling process) and stratified resampling, as the time-discretisation becomes finer and also identify their continuous-time limit, which is expressed as a suitably defined 'infinitesimal generator.' By contrasting these generators, we find that (certain modifica-tions of) systematic and SSP resampling 'dominate' stratified and indepen-dent 'killing' resampling in terms of their limiting overall resampling rate. The reduced intensity of resampling manifests itself in lower variance in our numerical experiment. This efficiency result, through an ordering of the re -sampling rate, is new to the literature. The second major contribution of this work concerns the analysis of the limiting behaviour of the entire population of particles of the particle filter as the time discretisation becomes finer. We provide the first proof, under general conditions, that the particle approxima-tion of the discretised continuous-time Feynman-Kac path integral models converges to a (uniformly weighted) continuous-time particle system.
引用
收藏
页码:3197 / 3222
页数:26
相关论文
共 28 条
[1]   Particle Markov chain Monte Carlo methods [J].
Andrieu, Christophe ;
Doucet, Arnaud ;
Holenstein, Roman .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-STATISTICAL METHODOLOGY, 2010, 72 :269-342
[2]   A duality formula and a particle Gibbs sampler for continuous time Feynman-Kac measures on path spaces [J].
Arnaudon, Marc ;
Del Moral, Pierre .
ELECTRONIC JOURNAL OF PROBABILITY, 2020, 25 :1-54
[3]   Julia: A Fresh Approach to Numerical Computing [J].
Bezanson, Jeff ;
Edelman, Alan ;
Karpinski, Stefan ;
Shah, Viral B. .
SIAM REVIEW, 2017, 59 (01) :65-98
[4]   A nonasymptotic theorem for unnormalized Feynman-Kac particle models [J].
Cerou, F. ;
Del Moral, P. ;
Guyader, A. .
ANNALES DE L INSTITUT HENRI POINCARE-PROBABILITES ET STATISTIQUES, 2011, 47 (03) :629-649
[5]  
Chopin N., 2020, An introduction to sequential Monte Carlo, DOI [10.1007/978-3-030-47845-2, DOI 10.1007/978-3-030-47845-2]
[6]  
CHOPIN N., 2022, SUPPLEMENT RESAMPLIN, DOI [10.1214/22-AOS2222SUPPA, DOI 10.1214/22-AOS2222SUPPA]
[7]  
CRISAN D., 1999, MARKOV PROCESS RELAT, V5, P293
[8]  
Del Moral P, 2000, STOCH PROC APPL, V86, P193
[9]  
Del Moral P, 2000, LECT NOTES MATH, V1729, P2
[10]  
DEL MORAL P., 2004, PROB APPL S, DOI 10 . 1007 / 978 - 1 - 4684 - 9393 - 1