Variance Reduction through Multilevel Monte Carlo Path Calculations

Author(s):  
Michael Giles
2019 ◽  
Vol 24 (1) ◽  
pp. 311-331 ◽  
Author(s):  
Prashant Kumar ◽  
Carmen Rodrigo ◽  
Francisco J. Gaspar ◽  
Cornelis W. Oosterlee

AbstractWe present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that is modeled using the Richards equation. We propose a stochastic extension for the empirical models that are typically employed to close the Richards equations. This is achieved by treating the soil parameters in these models as spatially correlated random fields with appropriately defined marginal distributions. As some of these parameters can only take values in a specific range, non-Gaussian models are utilized. The randomness in these parameters may result in path-wise highly nonlinear systems, so that a robust solver with respect to the random input is required. For this purpose, a solution method based on a combination of the modified Picard iteration and a cell-centered multigrid method for heterogeneous diffusion coefficients is utilized. Moreover, we propose a non-standard MLMC estimator to solve the resulting high-dimensional stochastic Richards equation. The improved efficiency of this multilevel estimator is achieved by parametric continuation that allows us to incorporate simpler nonlinear problems on coarser levels for variance reduction while the target strongly nonlinear problem is solved only on the finest level. Several numerical experiments are presented showing computational savings obtained by the new estimator compared with the original MC estimator.


2017 ◽  
Vol 6 (4) ◽  
pp. 121
Author(s):  
Hamid Reza Erfanian ◽  
Seyed Jaliledin Ghaznavi Bidgoli ◽  
Parvin Shakibaei

Monte Carlo simulation is one of the most common and popular method of options pricing. The advantages of this method are being easy to use, suitable for all kinds of standard and exotic options and also are suitable for higher dimensional problems. But on the other hand Monte Carlo variance convergence rate is which due to that it will have relatively slow convergence rate to answer the problems, as to achieve  accuracy when it has been d-dimensions, complexity is . For this purpose, several methods are provided in quasi Monte Carlo simulation to increase variance convergence rate as variance reduction techniques, so far. One of the latest presented methods is multilevel Monte Carlo that is introduced by Giles in 2008. This method not only reduces the complexity of computing amount  in use of Euler discretization scheme and the amount  in use of Milstein discretization scheme, but also has the ability to combine with other variance reduction techniques. In this paper, using Multilevel Monte Carlo method by taking Milstein discretization scheme, pricing spread option and compared complexity of computing with standard Monte Carlo method. The results of Multilevel Monte Carlo method in pricing spread options are better than standard Monte Carlo simulation.


2018 ◽  
Vol 482 (6) ◽  
pp. 627-630
Author(s):  
D. Belomestny ◽  
◽  
L. Iosipoi ◽  
N. Zhivotovskiy ◽  
◽  
...  

2020 ◽  
Vol 26 (2) ◽  
pp. 131-161
Author(s):  
Florian Bourgey ◽  
Stefano De Marco ◽  
Emmanuel Gobet ◽  
Alexandre Zhou

AbstractThe multilevel Monte Carlo (MLMC) method developed by M. B. Giles [Multilevel Monte Carlo path simulation, Oper. Res. 56 2008, 3, 607–617] has a natural application to the evaluation of nested expectations {\mathbb{E}[g(\mathbb{E}[f(X,Y)|X])]}, where {f,g} are functions and {(X,Y)} a couple of independent random variables. Apart from the pricing of American-type derivatives, such computations arise in a large variety of risk valuations (VaR or CVaR of a portfolio, CVA), and in the assessment of margin costs for centrally cleared portfolios. In this work, we focus on the computation of initial margin. We analyze the properties of corresponding MLMC estimators, for which we provide results of asymptotic optimality; at the technical level, we have to deal with limited regularity of the outer function g (which might fail to be everywhere differentiable). Parallel to this, we investigate upper and lower bounds for nested expectations as above, in the spirit of primal-dual algorithms for stochastic control problems.


2021 ◽  
Vol 433 ◽  
pp. 110164
Author(s):  
S. Ben Bader ◽  
P. Benedusi ◽  
A. Quaglino ◽  
P. Zulian ◽  
R. Krause

2020 ◽  
Vol 26 (3) ◽  
pp. 193-203
Author(s):  
Shady Ahmed Nagy ◽  
Mohamed A. El-Beltagy ◽  
Mohamed Wafa

AbstractMonte Carlo (MC) simulation depends on pseudo-random numbers. The generation of these numbers is examined in connection with the Brownian motion. We present the low discrepancy sequence known as Halton sequence that generates different stochastic samples in an equally distributed form. This will increase the convergence and accuracy using the generated different samples in the Multilevel Monte Carlo method (MLMC). We compare algorithms by using a pseudo-random generator and a random generator depending on a Halton sequence. The computational cost for different stochastic differential equations increases in a standard MC technique. It will be highly reduced using a Halton sequence, especially in multiplicative stochastic differential equations.


Sign in / Sign up

Export Citation Format

Share Document