Linear Variance Bounds for Particle Approximations of Time-Homogeneous Feynman-Kac Formulae

This article establishes sufficient conditions for a linear-in-time bound on the non-asymptotic variance for particle approximations of time-homogeneous Feynman–Kac formulae. These formulae appear in a wide variety of applications including option pricing in finance and risk sensitive control in engineering. In direct Monte Carlo approximation of these formulae, the non-asymptotic variance typically increases at an exponential rate in the time parameter. It is shown that a linear bound holds when a non-negative kernel, defined by the logarithmic potential function and Markov kernel which specify the Feynman–Kac model, satisfies a type of multiplicative drift condition and other regularity assumptions. Examples illustrate that these conditions are general and flexible enough to accommodate two rather extreme cases, which can occur in the context of a non-compact state space: (1) when the potential function is bounded above, not bounded below and the Markov kernel is not ergodic; and (2) when the potential function is not bounded above, but the Markov kernel itself satisfies a multiplicative drift condition.

[1]  F. Gland,et al.  STABILITY AND UNIFORM APPROXIMATION OF NONLINEAR FILTERS USING THE HILBERT METRIC AND APPLICATION TO PARTICLE FILTERS1 , 2004 .

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

[3]  G. Roberts,et al.  Polynomial convergence rates of Markov chains. , 2002 .

[4]  Pierre Del Moral,et al.  Stability of Feynman-Kac formulae with path-dependent potentials , 2009, 0910.4870.

[5]  A. Doucet,et al.  Particle Motions in Absorbing Medium with Hard and Soft Obstacles , 2004 .

[6]  P. Ney.  GENERAL IRREDUCIBLE MARKOV CHAINS AND NON‐NEGATIVE OPERATORS (Cambridge Tracts in Mathematics, 83) , 1986 .

[7]  E. Nummelin.  General irreducible Markov chains and non-negative operators: Preface , 1984 .

[8]  P. Moral,et al.  On the stability of interacting processes with applications to filtering and genetic algorithms , 2001 .

[9]  D. Crisan,et al.  Fundamentals of Stochastic Filtering , 2008 .

[10]  Nando de Freitas,et al.  Sequential Monte Carlo Methods in Practice , 2001, Statistics for Engineering and Information Science.

[11]  Pierre Del Moral,et al.  Particle approximations of Lyapunov exponents connected to Schrödinger operators and Feynman–Kac semigroups , 2003 .

[12]  P. Moral,et al.  Tree based functional expansions for Feynman–Kac particle models , 2009, 0906.4249.

[13]  Neil J. Gordon,et al.  Editors: Sequential Monte Carlo Methods in Practice , 2001 .

[14]  Lukasz Stettner,et al.  Risk-Sensitive Control of Discrete-Time Markov Processes with Infinite Horizon , 1999, SIAM J. Control. Optim..

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

[16]  S. Meyn.  Large deviation asymptotics and control variates for simulating large functions , 2006, math/0603328.

[17]  Ajay Jasra,et al.  Sequential Monte Carlo Methods for Option Pricing , 2010, 1005.4797.

[18]  A. Doucet,et al.  Sequential Monte Carlo methods for diffusion processes , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[19]  R. Handel.  Uniform time average consistency of Monte Carlo particle filters , 2008, 0812.0350.

[20]  Richard L. Tweedie,et al.  Markov Chains and Stochastic Stability , 1993, Communications and Control Engineering Series.

[21]  Pierre Del Moral,et al.  Feynman-Kac formulae , 2004 .

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

[23]  P. Moral,et al.  A nonasymptotic theorem for unnormalized Feynman-Kac particle models , 2011 .

[24]  A. Beskos,et al.  On the stability of sequential Monte Carlo methods in high dimensions , 2011, 1103.3965.

[25]  P. Whittle.  Risk-Sensitive Optimal Control , 1990 .

[26]  S. Meyn,et al.  Spectral theory and limit theorems for geometrically ergodic Markov processes , 2002, math/0209200.

[27]  Tze Leung Lai,et al.  A sequential Monte Carlo approach to computing tail probabilities in stochastic models , 2011 .

[28]  S. Ross,et al.  A theory of the term structure of interest rates'', Econometrica 53, 385-407 , 1985 .