Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

Multilevel branching splitting algorithm for estimating rare event probabilities

Abstract : We analyse the splitting algorithm performance in the estimation of rare event probabilities and this in a discrete multidimensional framework. For this we assume that each threshold is partitioned into disjoint subsets and the probability for a particle to reach the next threshold will depend on the starting subset. A straightforward estimator of the rare event probability is given by the proportion of simulated particles for which the rare event occurs. The variance of this estimator we get is the sum of two parts: one part resuming the variability due to each threshold and a second part resuming the variability due to the thresholds number. This decomposition is analogous to that of the continuous case. The optimal algorithm is then derived by cancelling the first term leading to optimal thresholds. Then we compare this variance with that of the algorithm in which one of the threshold has been deleted. Finally, we investigate the sensitivity of the variance of the estimator with respect to a shape deformation of an optimal threshold. As an example, we consider a two-dimensional Ornstein-Uhlenbeck process with conformal maps for shape deformation.
Document type :
Preprints, Working Papers, ...
Complete list of metadata

https://hal.archives-ouvertes.fr/hal-01269766
Contributor : Agnes Lagnoux Connect in order to contact the contributor
Submitted on : Monday, February 8, 2016 - 10:03:36 PM
Last modification on : Friday, December 24, 2021 - 3:18:04 PM
Long-term archiving on: : Saturday, November 12, 2016 - 11:13:07 AM

Files

simus_multidim_4.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-01269766, version 1
  • ARXIV : 1602.02893

Collections

Citation

Agnès Lagnoux, Pascal Lezaud. Multilevel branching splitting algorithm for estimating rare event probabilities. 2016. ⟨hal-01269766v1⟩

Share

Metrics

Les métriques sont temporairement indisponibles