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 条
[11]   Feynman-Kac Particle Integration with Geometric Interacting Jumps [J].
Del Moral, Pierre ;
Jacob, Pierre E. ;
Lee, Anthony ;
Murray, Lawrence ;
Peters, Gareth W. .
STOCHASTIC ANALYSIS AND APPLICATIONS, 2013, 31 (05) :830-871
[12]  
DelMoral P, 2013, MONOGR STAT APPL PRO, V126, P1
[13]   Comparison of resampling schemes for particle filtering [J].
Douc, R ;
Cappé, O ;
Moulines, E .
ISPA 2005: Proceedings of the 4th International Symposium on Image and Signal Processing and Analysis, 2005, :64-69
[14]  
Ethier S. N., 2005, WILEY SERIES PROBABI, DOI [10.1002/9780470316658, DOI 10.1002/9780470316658]
[15]   BATCH MEANS AND SPECTRAL VARIANCE ESTIMATORS IN MARKOV CHAIN MONTE CARLO [J].
Flegal, James M. ;
Jones, Galin L. .
ANNALS OF STATISTICS, 2010, 38 (02) :1034-1070
[16]   NEGATIVE ASSOCIATION, ORDERING AND CONVERGENCE OF RESAMPLING METHODS [J].
Gerber, Mathieu ;
Chopin, Nicolas ;
Whiteley, Nick .
ANNALS OF STATISTICS, 2019, 47 (04) :2236-2260
[17]   THE ASYMPTOTIC EFFICIENCY OF SIMULATION ESTIMATORS [J].
GLYNN, PW ;
WHITT, W .
OPERATIONS RESEARCH, 1992, 40 (03) :505-520
[18]   NOVEL-APPROACH TO NONLINEAR NON-GAUSSIAN BAYESIAN STATE ESTIMATION [J].
GORDON, NJ ;
SALMOND, DJ ;
SMITH, AFM .
IEE PROCEEDINGS-F RADAR AND SIGNAL PROCESSING, 1993, 140 (02) :107-113
[19]   An adaptive Metropolis algorithm [J].
Haario, H ;
Saksman, E ;
Tamminen, J .
BERNOULLI, 2001, 7 (02) :223-242
[20]   Monte Carlo filter using the genetic algorithm operators [J].
Higuchi, T .
JOURNAL OF STATISTICAL COMPUTATION AND SIMULATION, 1997, 59 (01) :1-23