Fast Ninomiya–Victoir calibration of the double-mean-reverting model

Abstract

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.     

Real-world jump-diffusion term structure models

This paper considers interest rate term structure models in a market attracting both continuous and discrete types of uncertainty. The event-driven noise is modelled by a Poisson random measure. Using as numeraire the growth optimal portfolio, interest rate derivatives are priced under the real-world probability measure. In particular, the real-world dynamics of the forward rates are derived and, for specific volatility structures, finite-dimensional Markovian representations are obtained. Furthermore, allowing for a stochastic short rate in a non-Markovian setting, a class of tractable affine term structures is derived where an equivalent risk-neutral probability measure may not exist.     

A new closed-form solution as an extension of the Black–Scholes formula allowing smile curve plotting

In this paper, as a generalization of the Black–Scholes (BS) model, we elaborate a new closed-form solution for a uni-dimensional European option pricing model called the J-model. This closed-form solution is based on a new stochastic process, called the J-process, which is an extension of the Wiener process satisfying the martingale property. The J-process is based on a new statistical law called the J-law, which is an extension of the normal law. The J-law relies on four parameters in its general form. It has interesting asymmetry and tail properties, allowing it to fit the reality of financial markets with good accuracy, which is not the case for the normal law. Despite the use of one state variable, we find results similar to those of Heston dealing with the bi-dimensional stochastic volatility problem for pricing European calls. Inverting the BS formula, we plot the smile curve related to this closed-form solution. The J-model can also serve to determine the implied volatility by inverting the J-formula and can be used to price other kinds of options such as American options.     

Pricing equity options everywhere

Finite difference methods are a popular technique for pricing American options. Since their introduction to finance by Brennan and Schwartz their use has spread from vanilla calls and puts on one stock to path-dependent and exotic options on multiple assets. Despite the breadth of the problems they have been applied to, and the increased sophistication of some of the newer techniques, most approaches to pricing equity options have not adequately addressed the issues of unbounded computational domains and divergent diffusion coefficients. In this article it is shown that these two problems are related and can be overcome using multiple grids. This new technique allows options to be priced for all values of the underlying, and is illustrated using standard put options and the call on the maximum of two stocks. For the latter contract, I also derive a characterization of the asymptotic continuation region in terms of a one-dimensional option pricing problem, and give analytic formulae for the perpetual case.     

Turbocharging Monte Carlo pricing for the rough Bergomi model

The rough Bergomi model, introduced by Bayer et al. [Quant. Finance, 2016, 16(6), 887–904], is one of the recent rough volatility models that are consistent with the stylised fact of implied volatility surfaces being essentially time-invariant, and are able to capture the term structure of skew observed in equity markets. In the absence of analytical European option pricing methods for the model, we focus on reducing the runtime-adjusted variance of Monte Carlo implied volatilities, thereby contributing to the model’s calibration by simulation. We employ a novel composition of variance reduction methods, immediately applicable to any conditionally log-normal stochastic volatility model. Assuming one targets implied volatility estimates with a given degree of confidence, thus calibration RMSE, the results we demonstrate equate to significant runtime reductions—roughly 20 times on average, across different correlation regimes.     

A multivariate Lévy process model with linear correlation

In this paper, we develop a multivariate risk-neutral Lévy process model and discuss its applicability in the context of the volatility smile of multiple assets. Our formulation is based upon a linear combination of independent univariate Lévy processes and can easily be calibrated to a set of one-dimensional marginal distributions and a given linear correlation matrix. We derive conditions for our formulation and the associated calibration procedure to be well-defined and provide some examples associated with particular Lévy processes permitting a closed-form characteristic function. Numerical results of the option premiums on three currencies are presented to illustrate the effectiveness of our formulation with different linear correlation structures.     

Regression-based Monte Carlo methods for stochastic control models: variable annuities with lifelong guarantees

We present regression-based Monte Carlo simulation algorithm for solving the stochastic control models associated with pricing and hedging of the guaranteed lifelong withdrawal benefit (GLWB) in variable annuities, where the dynamics of the underlying fund value is assumed to evolve according to the stochastic volatility model. The GLWB offers a lifelong withdrawal benefit, even when the policy account value becomes zero, while the policyholder remains alive. Upon death, the remaining account value will be paid to the beneficiary as a death benefit. The bang-bang control strategy analysed under the assumption of maximization of the policyholder’s expected cash flow reduces the strategy space of optimal withdrawal policies to three choices: zero withdrawal, withdrawal at the contractual amount or complete surrender. The impact on the GLWB value under various withdrawal behaviours of the policyholder is examined. We also analyse the pricing properties of GLWB subject to different model parameter values and structural features.     

Pricing via recursive quantization in stochastic volatility models

We provide the first recursive quantization-based approach for pricing options in the presence of stochastic volatility. This method can be applied to any model for which an Euler scheme is available for the underlying price process and it allows one to price vanillas, as well as exotics, thanks to the knowledge of the transition probabilities for the discretized stock process. We apply the methodology to some celebrated stochastic volatility models, including the Stein and Stein [Rev. Financ. Stud. 1991, (4), 727–752] model and the SABR model introduced in Hagan et al. [Wilmott Mag., 2002, 84–108]. A numerical exercise shows that the pricing of vanillas turns out to be accurate; in addition, when applied to some exotics like equity-volatility options, the quantization-based method overperforms by far the Monte Carlo simulation.     

A PDE method for estimation of implied volatility

In this paper it is proved that the Black–Scholes implied volatility satisfies a second order non-linear partial differential equation. The obtained PDE is then used to construct an algorithm for fast and accurate polynomial approximation for Black–Scholes implied volatility that improves on the existing numerical schemes from literature, both in speed and parallelizability. We also show that the method is applicable to other problems, such as approximation of implied Bachelier volatility.     

American option pricing with discrete and continuous time models: An empirical comparison

This paper considers discrete time GARCH and continuous time SV models and uses these for American option pricing. We first of all show that with a particular choice of framework the parameters of the SV models can be estimated using simple maximum likelihood techniques. We then perform a Monte Carlo study to examine their differences in terms of option pricing, and we study the convergence of the discrete time option prices to their implied continuous time values. Finally, a large scale empirical analysis using individual stock options and options on an index is performed comparing the estimated prices from discrete time models to the corresponding continuous time model prices. The results show that, while the overall differences in performance are small, for the in the money put options on individual stocks the continuous time SV models do generally perform better than the discrete time GARCH specifications.