ON A GENERALIZED SPLITTING METHOD FOR SAMPLING FROM A CONDITIONAL DISTRIBUTION

We study the behavior of a splitting method for sampling from a given distribution conditional on the occurrence of a rare event. The method returns a random-sized sample of points such that unconditionally on the sample size, each point is distributed exactly according to the original distribution conditional on the rare event. For a cost function which is nonzero only when the rare event occurs, the method provides an unbiased estimator of the expected cost, but if we select at random one of the returned points, its distribution differs in general from the exact conditional distribution given the rare event. However, we prove that if we repeat the algorithm n times and select one of the returned points at random, the distribution of the selected point converges to the exact one in total variation when n increases.

[1]  Pierre L'Ecuyer,et al.  Reliability Estimation for Networks with Minimal Flow Demand and Random Link Capacities , 2018, 1805.03326.

[2]  Peter W. Glynn,et al.  Stochastic Simulation: Algorithms and Analysis , 2007 .

[3]  Gerardo Rubino,et al.  Static Network Reliability Estimation via Generalized Splitting , 2013, INFORMS J. Comput..

[4]  Bruno Tuffin,et al.  Splitting for rare-event simulation , 2006, WSC.

[5]  Paul Glasserman,et al.  Multilevel Splitting for Estimating Rare Event Probabilities , 1999, Oper. Res..

[6]  Bruno Tuffin,et al.  Static Network Reliability Estimation under the Marshall-Olkin Copula , 2015, ACM Trans. Model. Comput. Simul..

[7]  Yasuhiro Yamai,et al.  Value-at-risk versus expected shortfall: A practical perspective , 2005 .

[8]  Dirk P. Kroese,et al.  Efficient Monte Carlo simulation via the generalized splitting method , 2012, Stat. Comput..

[9]  Robert E. Shannon,et al.  Design and analysis of simulation experiments , 1978, WSC '78.

[10]  Bruno Tuffin,et al.  Rare events, splitting, and quasi-Monte Carlo , 2007, TOMC.

[11]  W. Hoeffding Probability Inequalities for sums of Bounded Random Variables , 1963 .

[12]  Pierre L'Ecuyer,et al.  Splitting Techniques , 2009, Rare Event Simulation using Monte Carlo Methods.

[13]  James G. Scott,et al.  The Bayesian bridge , 2011, 1109.2279.

[14]  Ludovic Goudenège,et al.  Unbiasedness of some generalized Adaptive Multilevel Splitting algorithms , 2016 .

[15]  Bruno Tuffin,et al.  Markov chain importance sampling with applications to rare event probability estimation , 2011, Stat. Comput..

[16]  Christian Léger,et al.  Bootstrap confidence intervals for ratios of expectations , 1999, TOMC.

[17]  Jan-Kees C. W. van Ommeren,et al.  On the importance function in splitting simulation , 2002, Eur. Trans. Telecommun..