Importance sampling and statistical Romberg method for Lévy processes

详细信息    查看全文

• 作者：Mohamed Ben Alaya ^{mba@math.univ-paris13.fr" class="auth_mail" title="E-mail the corresponding author} ; Kaouther Hajji ^{hajji@math.univ-paris13.fr" class="auth_mail" title="E-mail the corresponding author} ; Ahmed Kebaier ; ^{kebaier@math.univ-paris13.fr" class="auth_mail" title="E-mail the corresponding author}
• 关键词：60E07 ; 60G51 ; 60F05 ; 62L20 ; 65C05 ; 60H35
• 刊名：Stochastic Processes and their Applications
• 出版年：2016
• 出版时间：July 2016
• 年：2016
• 卷：126
• 期：7
• 页码：1901-1931
• 全文大小：1286 K

文摘

An important family of stochastic processes arising in many areas of applied probability is the class of Lévy processes. Generally, such processes are not simulatable especially for those with infinite activity. In practice, it is common to approximate them by truncating the jumps at some cut-off size k to view the MathML source">ε$\epsilon$ (k to view the MathML source">ε↘0$\epsilon \searrow 0$). This procedure leads us to consider a simulatable compound Poisson process. This paper first introduces, for this setting, the statistical Romberg method to improve the complexity of the classical Monte Carlo method. Roughly speaking, we use many sample paths with a coarse cut-off k to view the MathML source">ε^β${\epsilon }^{\beta }$, k to view the MathML source">β∈(0,1)$\beta \in (0,1)$, and few additional sample paths with a fine cut-off k to view the MathML source">ε$\epsilon$. Central limit theorems of Lindeberg–Feller type for both Monte Carlo and statistical Romberg method for the inferred errors depending on the parameter k to view the MathML source">ε$\epsilon$ are proved with explicit formulas for the limit variances. This leads to an accurate description of the optimal choice of parameters. Afterwards, the authors propose a stochastic approximation method in order to find the optimal measure change by Esscher transform for Lévy processes with Monte Carlo and statistical Romberg importance sampling variance reduction. Furthermore, we develop new adaptive Monte Carlo and statistical Romberg algorithms and prove the associated central limit theorems. Finally, numerical simulations are processed to illustrate the efficiency of the adaptive statistical Romberg method that reduces at the same time the variance and the computational effort associated to the effective computation of option prices when the underlying asset process follows an exponential pure jump CGMY model.