Sequential Monte Carlo for rare event estimation

  • Authors:
  • F. Cérou;P. Moral;T. Furon;A. Guyader

  • Affiliations:
  • INRIA Rennes - Bretagne Atlantique, Rennes Cedex, France 35042;INRIA Bordeaux Sud-Ouest & Institut de Mathéématiques de Bordeaux, Université Bordeaux 1, Talence Cedex, France 33405;INRIA Rennes - Bretagne Atlantique, Rennes Cedex, France 35042;Equipe de Statistique, Université de Haute Bretagne, Rennes Cedex, France 35043

  • Venue:
  • Statistics and Computing
  • Year:
  • 2012

Quantified Score

Hi-index 0.00

Visualization

Abstract

This paper discusses a novel strategy for simulating rare events and an associated Monte Carlo estimation of tail probabilities. Our method uses a system of interacting particles and exploits a Feynman-Kac representation of that system to analyze their fluctuations. Our precise analysis of the variance of a standard multilevel splitting algorithm reveals an opportunity for improvement. This leads to a novel method that relies on adaptive levels and produces, in the limit of an idealized version of the algorithm, estimates with optimal variance. The motivation for this theoretical work comes from problems occurring in watermarking and fingerprinting of digital contents, which represents a new field of applications of rare event simulation techniques. Some numerical results show performance close to the idealized version of our technique for these practical applications.