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 δ and δβ,0<β<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 β (which turns out to be β=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.
References
Aldous, D. J. and Eagleson, G. K. (1978). On mixing and stability of limit theorems. Ann. Probab. 6 325–331. MR517416 0376.60026 10.1214/aop/1176995577 euclid.aop/1176995577 Aldous, D. J. and Eagleson, G. K. (1978). On mixing and stability of limit theorems. Ann. Probab. 6 325–331. MR517416 0376.60026 10.1214/aop/1176995577 euclid.aop/1176995577
Bally, V. and Talay, D. (1996). The law of the Euler scheme for stochastic differential equations. I. Convergence rate of the distribution function. Probab. Theory Related Fields 104 43–60. MR1367666 0838.60051 10.1007/BF01303802 Bally, V. and Talay, D. (1996). The law of the Euler scheme for stochastic differential equations. I. Convergence rate of the distribution function. Probab. Theory Related Fields 104 43–60. MR1367666 0838.60051 10.1007/BF01303802
Boyle, P., Broadie, M. and Glasserman, P. (1997). Monte Carlo methods for security pricing. J. Econom. Dynam. Control 21 1267–1321. MR1470283 0901.90007 10.1016/S0165-1889(97)00028-6 Boyle, P., Broadie, M. and Glasserman, P. (1997). Monte Carlo methods for security pricing. J. Econom. Dynam. Control 21 1267–1321. MR1470283 0901.90007 10.1016/S0165-1889(97)00028-6
Broadie, M. and Detemple, J. (1997). Recent advances in numerical methods for pricing derivative securities. In Numerical Methods in Finance 43–66. Cambridge Univ. Press. MR1470508 0898.90029 Broadie, M. and Detemple, J. (1997). Recent advances in numerical methods for pricing derivative securities. In Numerical Methods in Finance 43–66. Cambridge Univ. Press. MR1470508 0898.90029
Duffie, D. and Glynn, P. (1995). Efficient Monte Carlo simulation of security prices. Ann. Appl. Probab. 5 897–905. MR1384358 0877.65099 10.1214/aoap/1177004598 euclid.aoap/1177004598 Duffie, D. and Glynn, P. (1995). Efficient Monte Carlo simulation of security prices. Ann. Appl. Probab. 5 897–905. MR1384358 0877.65099 10.1214/aoap/1177004598 euclid.aoap/1177004598
Jacod, J. (1997). On continuous conditional Gaussian martingales and stable convergence in law. Séminaire de Probabilités XXXI. Lecture Notes in Math. 1655 232–246. Springer, Berlin. MR1478732 Jacod, J. (1997). On continuous conditional Gaussian martingales and stable convergence in law. Séminaire de Probabilités XXXI. Lecture Notes in Math. 1655 232–246. Springer, Berlin. MR1478732
Jacod, J., Kurtz, T. G., Méléard, S. and Protter, P. (2005). The approximate Euler method for Lévy driven stochastic differential equations. Ann. Inst. H. Poincaré. Special Issue Devoted to the Memory of P. A. Meyer 41 523–558. MR2139032 10.1016/j.anihpb.2004.01.007 Jacod, J., Kurtz, T. G., Méléard, S. and Protter, P. (2005). The approximate Euler method for Lévy driven stochastic differential equations. Ann. Inst. H. Poincaré. Special Issue Devoted to the Memory of P. A. Meyer 41 523–558. MR2139032 10.1016/j.anihpb.2004.01.007
Jacod, J. and Protter, P. (1998). Asymptotic error distributions for the Euler method for stochastic differential equations. Ann. Probab. 26 267–307. MR1617049 0937.60060 10.1214/aop/1022855419 euclid.aop/1022855419 Jacod, J. and Protter, P. (1998). Asymptotic error distributions for the Euler method for stochastic differential equations. Ann. Probab. 26 267–307. MR1617049 0937.60060 10.1214/aop/1022855419 euclid.aop/1022855419
Kebaier, A. (2004). Romberg extrapolation, variance reduction and applications to option pricing. Preprint. 1099.65011 10.1214/105051605000000511 euclid.aoap/1133965776 Kebaier, A. (2004). Romberg extrapolation, variance reduction and applications to option pricing. Preprint. 1099.65011 10.1214/105051605000000511 euclid.aoap/1133965776
Kurtz, T. G. and Protter, P. (1991). Wong–Zakai corrections, random evolutions, and simulation schemes for SDEs. In Stochastic Analysis 331–346. Academic Press, Boston, MA. MR1119837 0762.60047 Kurtz, T. G. and Protter, P. (1991). Wong–Zakai corrections, random evolutions, and simulation schemes for SDEs. In Stochastic Analysis 331–346. Academic Press, Boston, MA. MR1119837 0762.60047
Lamberton, D. and Lapeyre, B. (1997). Introduction au calcul stochastique appliqué à la finance, 2nd ed. Ellipses, Paris. MR1607509 0949.60005 Lamberton, D. and Lapeyre, B. (1997). Introduction au calcul stochastique appliqué à la finance, 2nd ed. Ellipses, Paris. MR1607509 0949.60005
Protter, P. (1990). Stochastic Integration and Differential Equations. Springer, Berlin. MR1037262 0694.60047 Protter, P. (1990). Stochastic Integration and Differential Equations. Springer, Berlin. MR1037262 0694.60047
Talay, D. and Tubaro, L. (1990). Expansion of the global error for numerical schemes solving stochastic differential equations. Stochastic Anal. Appl. 8 483–509. MR1091544 0718.60058 10.1080/07362999008809220 Talay, D. and Tubaro, L. (1990). Expansion of the global error for numerical schemes solving stochastic differential equations. Stochastic Anal. Appl. 8 483–509. MR1091544 0718.60058 10.1080/07362999008809220
