Island filters for partially observed spatiotemporal systems
share
Statistical inference for high-dimensional partially observed, nonlinear, stochastic processes is a methodological challenge with applications including spatiotemporal analysis of epidemiological and ecological systems. Standard particle filter algorithms, which provide an effective approach for general low-dimensional partially observed Markov processes, suffer from a curse of dimensionality that limits their applicability beyond low-dimensional systems. We show that many independent Monte Carlo calculations, none of which, by itself, attempts to solve the filtering problem, can be combined to give a global filtering solution that theoretically beats the curse of dimensionality under weak coupling conditions. The independent Monte Carlo calculations are called islands, and the corresponding algorithm is called a basic island filter (BIF). Carrying out an operation called adapted simulation on each island results in a variant called an adapted simulation island filter (ASIF). Adapted simulation can be implemented using a Monte Carlo technique called intermediate resampling to give the ASIF-IR algorithm which has improved theoretical and empirical scaling. Our focus is on evaluation of the likelihood function, but our algorithms and theory are also pertinent to latent state estimation. We demonstrate our methodology on coupled population dynamics arising in the epidemiology of infectious disease. In this context, our methods retain the generality of particle filter algorithms while scaling to systems an order of magnitude larger.
READ FULL TEXT