Fast Ninomiya-Victoir Calibration of the Double-Mean-Reverting Model

We consider the three-factor double mean reverting (DMR) option pricing model of Gatheral [ Consistent Modelling of SPX and VIX Options , 2008], a model which can be successfully calibrated to both VIX options and SPX options simultaneously. One drawback of this model is that calibration may be slow because no closed form solution for European options exists. In this paper, we apply modified versions of the second-order Monte Carlo scheme of Ninomiya and Victoir [ Appl. Math. Finance , 2008, 15 , 107--121], and compare these to the Euler--Maruyama scheme with full truncation of Lord et al. [ Quant. Finance , 2010, 10 (2), 177--194], demonstrating on the one hand that fast calibration of the DMR model is practical, and on the other that suitably modified Ninomiya--Victoir schemes are applicable to the simulation of much more complicated time-homogeneous models than may have been thought previously.

[1]  L. Bergomi.  Smile Dynamics IV , 2009 .

[2]  Terry Lyons,et al.  Cubature on Wiener space , 2004, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[3]  Jorge J. Moré,et al.  User guide for MINPACK-1. [In FORTRAN] , 1980 .

[4]  P. Hagan,et al.  MANAGING SMILE RISK , 2002 .

[5]  Ioannis Karatzas,et al.  Brownian Motion and Stochastic Calculus , 1987 .

[6]  Aurélien Alfonsi,et al.  High order discretization schemes for the CIR process: Application to affine term structure and Heston models , 2010, Math. Comput..

[7]  Hans Buehler,et al.  Consistent Variance Curve Models , 2006, Finance Stochastics.

[8]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[9]  S. Ninomiya,et al.  Weak Approximation of Stochastic Differential Equations and Application to Derivative Pricing , 2006, math/0605361.

[10]  G. Strang.  Accurate partial difference methods I: Linear cauchy problems , 1963 .

[11]  Frances Y. Kuo,et al.  Constructing Sobol Sequences with Better Two-Dimensional Projections , 2008, SIAM J. Sci. Comput..

[12]  Jim Gatheral.  Consistent Modeling of SPX and VIX options , 2008 .

[13]  Jim Gatheral,et al.  Arbitrage-free SVI volatility surfaces , 2012, 1204.0646.

[14]  Jim Gatheral.  The Volatility Surface: A Practitioner's Guide , 2006 .

[15]  Gilbert Strang,et al.  Accurate partial difference methods , 1964 .

[16]  Marcos Eugênio da Silva,et al.  Quasi-Monte Carlo in finance: extending for problems of high effective dimension , 2005 .

[17]  Jorge J. Moré,et al.  User Guide for Minpack-1 , 1980 .

[18]  L. Bergomi.  Smile Dynamics II , 2005 .

[19]  Paul Glasserman,et al.  Monte Carlo Methods in Financial Engineering , 2003 .

[20]  Peter K. Friz,et al.  Semi-closed form cubature and applications to financial diffusion models , 2010 .

[21]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[22]  F. Diebold,et al.  The distribution of realized stock return volatility , 2001 .

[23]  Syoiti Ninomiya,et al.  A new higher-order weak approximation scheme for stochastic differential equations and the Runge–Kutta method , 2009, Finance Stochastics.

[24]  D. Talay,et al.  Expansion of the global error for numerical schemes solving stochastic differential equations , 1990 .

[25]  D. Dijk,et al.  A comparison of biased simulation schemes for stochastic volatility models , 2008 .