Hostname: page-component-848d4c4894-pftt2 Total loading time: 0 Render date: 2024-06-07T07:56:10.626Z Has data issue: false hasContentIssue false

Exact simulation for multivariate Itô diffusions

Published online by Cambridge University Press:  03 December 2020

Jose Blanchet*
Affiliation:
Stanford University
Fan Zhang*
Affiliation:
Stanford University
*
*Postal address: Huang Engineering Center, 475 Via Ortega, Stanford, CA 94305, United States.
*Postal address: Huang Engineering Center, 475 Via Ortega, Stanford, CA 94305, United States.

Abstract

We provide the first generic exact simulation algorithm for multivariate diffusions. Current exact sampling algorithms for diffusions require the existence of a transformation which can be used to reduce the sampling problem to the case of a constant diffusion matrix and a drift which is the gradient of some function. Such a transformation, called the Lamperti transformation, can be applied in general only in one dimension. So, completely different ideas are required for the exact sampling of generic multivariate diffusions. The development of these ideas is the main contribution of this paper. Our strategy combines techniques borrowed from the theory of rough paths, on the one hand, and multilevel Monte Carlo on the other.

Type
Original Article
Copyright
© Applied Probability Trust 2020

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Ait-Sahalia, Y. (2008). Closed-form likelihood expansions for multivariate diffusions. Ann. Statist. 36, 906937.CrossRefGoogle Scholar
Aronson, D. (1967). Bounds for the fundamental solution of a parabolic equation. Bull. Amer. Math. Soc. 73, 890896.CrossRefGoogle Scholar
Bally, V. (2006). Lower bounds for the density of locally elliptic Itô processes. Ann. Prob. 34, 24062440.CrossRefGoogle Scholar
Beskos, A., Papaspiliopoulos, O. and Roberts, G. O. (2006). Retrospective exact simulation of diffusion sample paths with applications. Bernoulli 12, 10771098.CrossRefGoogle Scholar
Beskos, A. and Roberts, G. O. (2005). Exact simulation of diffusions. Ann. Appl. Prob. 15, 24222444.CrossRefGoogle Scholar
Blanchet, J., Chen, X. and Dong, J. (2017). $\epsilon$-strong simulation for multidimensional stochastic differential equations via rough path analysis. Ann. Appl. Prob. 27, 275336.CrossRefGoogle Scholar
Chen, N. and Huang, Z. (2013). Localization and exact simulation of Brownian motion-driven stochastic differential equations. Math. Operat. Res. 38, 591616.CrossRefGoogle Scholar
Fearnhead, P., Papaspiliopoulos, O., Roberts, G. O. and Stuart, A. (2010). Random-weight particle filtering of continuous time processes. J. R. Statist. Soc. B [Statist. Methodology] 72, 497512.CrossRefGoogle Scholar
Fleming, W. H. and Rishel, R. W. (2012). Deterministic and Stochastic Optimal Control. Springer Science & Business Media, New York.Google Scholar
Friedman, A. (2013). Partial Differential Equations of Parabolic Type. Courier Corporation, Chelmsford, MA.Google Scholar
Giles, M. B. (2008). Multilevel Monte Carlo path simulation. Operat. Res. 56, 607617.CrossRefGoogle Scholar
Giles, M. B. and Szpruch, L. (2014). Antithetic multilevel Monte Carlo estimation for multi-dimensional SDEs without Lévy area simulation. Ann. Appl. Prob. 24, 15851620.CrossRefGoogle Scholar
Henry-Labordère, P., Tan, X., Touzi, N. (2017). Unbiased simulation of stochastic differential equations. Ann. Appl. Prob. 27, 33053341.CrossRefGoogle Scholar
Huber, M. (2016). Nearly optimal Bernoulli factories for linear functions. Combinatorics Prob. Comput. 25, 577591.CrossRefGoogle Scholar
Karatzas, I. and Shreve, S. (2012). Brownian Motion and Stochastic Calculus. Springer Science & Business Media, New York.Google Scholar
Kusuoka, S. and Stroock, D. (1987). Applications of the Malliavin calculus, part 3. J. Fac. Sci. Univ. Tokyo Sect. IA Math. 34, 397442.Google Scholar
Łatuszyńki, K., Kosmidis, I., Papaspiliopoulos, O. and Roberts, G. O. (2011). Simulating events of unknown probabilities via reverse time martingales. Random Structures Algorithms 38, 441452.CrossRefGoogle Scholar
McLeish, D. (2011). A general method for debiasing a Monte Carlo estimator. Monte Carlo Meth. Appl. 17, 301315.CrossRefGoogle Scholar
Nacu, Ş. and Peres, Y. (2005). Fast simulation of new coins from old. Ann. Appl. Prob. 15, 93115.CrossRefGoogle Scholar
Revuz, D. and Yor, M. (2013). Continuous Martingales and Brownian Motion. Springer Science & Business Media, New York.Google Scholar
Rhee, C.-H. and Glynn, P. W. (2012). A new approach to unbiased estimation for SDE’s. In WSC ’12: Proceedings of the Winter Simulation Conference, Berlin, pp. 17.CrossRefGoogle Scholar
Rhee, C.-H. and Glynn, P. W. (2015). Unbiased estimation with square root convergence for SDE models. Operat. Res. 63, 10261043.CrossRefGoogle Scholar
Sheu, S.-J. (1991). Some estimates of the transition density of a nondegenerate diffusion Markov process. Ann. Prob. 19, 538561.CrossRefGoogle Scholar