Statistical Romberg extrapolation: A new variance reduction method and applications to option pricing

Abstract

We study the approximation of $\mathbb{E}f(X_T)$ by a Monte Carlo algorithm, where $X$ is the solution of a stochastic differential equation and $f$ is a given function. We introduce a new variance reduction method, which can be viewed as a statistical analogue of Romberg extrapolation method. Namely, we use two Euler schemes with steps $\delta$ and ${\delta }^{\beta },0<\beta <1$. This leads to an algorithm which, for a given level of the statistical error, has a complexity significantly lower than the complexity of the standard Monte Carlo method. We analyze the asymptotic error of this algorithm in the context of general (possibly degenerate) diffusions. In order to find the optimal $\beta$ (which turns out to be $\beta =1/2$), we establish a central limit type theorem, based on a result of Jacod and Protter for the asymptotic distribution of the error in the Euler scheme. We test our method on various examples. In particular, we adapt it to Asian options. In this setting, we have a CLT and, as a by-product, an explicit expansion of the discretization error.

Published at http://dx.doi.org/10.1214/105051605000000511 in the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

References in corpus (1)

• Statistical Romberg extrapolation: A new variance reduction method and applications to option pricing

Cited by in corpus (17)

• Statistical Romberg extrapolation: A new variance reduction method and applications to option pricing
• A Continuation Multilevel Monte Carlo algorithm
• On a perturbation theory and on strong convergence rates for stochastic ordinary and partial differential equations with non-globally monotone coefficients
• Divergence of the multilevel Monte Carlo Euler method for nonlinear stochastic differential equations
• Multilevel Picard iterations for solving smooth semilinear parabolic heat equations
• Optimization of mesh hierarchies in Multilevel Monte Carlo samplers
• A multilevel Monte Carlo method for computing failure probabilities
• Multilevel Richardson-Romberg extrapolation
• Central limit theorem for the multilevel Monte Carlo Euler method
• Iterative Particle Approximation for McKean-Vlasov SDEs with application to Multilevel Monte Carlo estimation
• Asymptotic error for the Milstein scheme for SDEs driven by continuous semimartingales
• Consistency and Stability of a Milstein-Galerkin Finite Element Scheme for Semilinear SPDE
• Antithetic multilevel particle system sampling method for McKean-Vlasov SDEs
• Antithetic multilevel Monte Carlo method for approximations of SDEs with non-globally Lipschitz continuous coefficients
• On the Acceleration of the Multi-Level Monte Carlo Method
• Central Limit Theorem for the $\sigma$-antithetic multilevel Monte Carlo method
• Importance sampling for a robust and efficient multilevel Monte Carlo estimator for stochastic reaction networks