Emergent Mind

Estimation of distributions via multilevel Monte Carlo with stratified sampling

(1906.00126)
Published Jun 1, 2019 in math.NA and cs.NA

Abstract

We design and implement a novel algorithm for computing a multilevel Monte Carlo (MLMC) estimator of the cumulative distribution function of a quantity of interest in problems with random input parameters or initial conditions. Our approach combines a standard MLMC method with stratified sampling by replacing standard Monte Carlo at each level with stratified Monte Carlo with proportional allocation. We show that the resulting stratified MLMC algorithm is more efficient than its standard MLMC counterpart, due to the reduction in variance at each level provided by the stratification of the random parameter's domain. A smoothing approximation for the indicator function based on kernel density estimation yields a more efficient algorithm compared to the typically used polynomial smoothing. The difference in computational cost between the smoothing methods depends on the required error tolerance.

We're not able to analyze this paper right now due to high demand.

Please check back later (sorry!).

Generate a summary of this paper on our Pro plan:

We ran into a problem analyzing this paper.

Newsletter

Get summaries of trending comp sci papers delivered straight to your inbox:

Unsubscribe anytime.