Multilevel sequential Monte Carlo: Mean square error bounds under verifiable conditions

ABSTRACT In this article, we consider the multilevel sequential Monte Carlo (MLSMC) method of Beskos et al. (Stoch. Proc. Appl. [to appear]). This is a technique designed to approximate expectations w.r.t. probability laws associated to a discretization. For instance, in the context of inverse problems, where one discretizes the solution of a partial differential equation. The MLSMC approach is especially useful when independent, coupled sampling is not possible. Beskos et al. show that for MLSMC the computational effort to achieve a given error, can be less than independent sampling. In this article we significantly weaken the assumptions of Beskos et al., extending the proofs to non-compact state-spaces. The assumptions are based upon multiplicative drift conditions as in Kontoyiannis and Meyn (Electron. J. Probab. 10 [2005]: 61–123). The assumptions are verified for an example.

[1]  N. Whiteley.  Stability properties of some particle filters , 2011, 1109.6779.

[2]  A. Beskos,et al.  Multilevel sequential Monte Carlo samplers , 2015, 1503.07259.

[3]  Ajay Jasra,et al.  On the Behaviour of the Backward Interpretation of Feynman-Kac Formulae Under Verifiable Conditions , 2013, J. Appl. Probab..

[4]  Stefan Heinrich,et al.  Monte Carlo Complexity of Global Solution of Integral Equations , 1998, J. Complex..

[5]  Stefan Heinrich,et al.  Monte Carlo Complexity of Parametric Integration , 1999, J. Complex..

[6]  P. Moral.  Feynman-Kac Formulae: Genealogical and Interacting Particle Systems with Applications , 2004 .

[7]  P. Moral,et al.  Sequential Monte Carlo samplers , 2002, cond-mat/0212648.

[8]  Andrew M. Stuart,et al.  Complexity analysis of accelerated MCMC methods for Bayesian inversion , 2012, 1207.2411.

[9]  Ajay Jasra,et al.  Forward and Inverse Uncertainty Quantification using Multilevel Monte Carlo Algorithms for an Elliptic Nonlocal Equation , 2016, 1603.06381.

[10]  S. Meyn,et al.  Large Deviations Asymptotics and the Spectral Theory of Multiplicatively Regular Markov Processes , 2005, math/0509310.

[11]  Nick Whiteley,et al.  Linear Variance Bounds for Particle Approximations of Time-Homogeneous Feynman-Kac Formulae , 2012 .

[12]  Michael B. Giles,et al.  Multilevel Monte Carlo Path Simulation , 2008, Oper. Res..

[13]  P. Moral,et al.  On adaptive resampling strategies for sequential Monte Carlo methods , 2012, 1203.0464.

[14]  N. Whiteley.  Sequential Monte Carlo Samplers: Error Bounds and Insensitivity to Initial Conditions , 2011, 1103.3970.