An empirical comparison of GARCH option pricing models

Recent empirical studies have shown that GARCH models can be successfully used to describe option prices. Pricing such contracts requires knowledge of the risk neutral cumulative return distribution. Since the analytical forms of these distributions are generally unknown, computationally intensive numerical schemes are required for pricing to proceed. Heston and Nandi (2000) consider a particular GARCH structure that permits analytical solutions for pricing European options and they provide empirical support for their model. The analytical tractability comes at a potential cost of realism in the underlying GARCH dynamics. In particular, their model falls in the affine family, whereas most GARCH models that have been examined fall in the non-affine family. This article takes a closer look at this model with the objective of establishing whether there is a cost to restricting focus to models in the affine family. We confirm Heston and Nandi's findings, namely that their model can explain a significant portion of the volatility smile. However, we show that a simple non affine NGARCH option model is superior in removing biases from pricing residuals for all moneyness and maturity categories especially for out-the-money contracts. The implications of this finding are examined. JEL Classification G13     

A General Affine Earnings Valuation Model

We introduce a methodology, with two applications, that incorporates stochastic interest rates, heteroskedasticity and risk aversion into the residual income model. In the first application, goodwill is an affine (constant plus linear term) function where the constant and linear coefficients are time-varying. Homoskedastic risk gives rise to a constant risk premium, while heteroskedastic risk gives rise to linear state-dependent risk premiums. In the second application, we present a class of models where a non-linear function for the price-to-book ratio can be derived. We show how interest rates, risk, profitability and growth affect the price-to-book ratio.     

American option pricing under the double Heston model based on asymptotic expansion

This paper focuses on pricing American put options under the double Heston model proposed by Christoffersen et al. By introducing an explicit exercise rule, we obtain the asymptotic expansion of the solution to the partial differential equation for pricing American put options. We calculate American option price by the sum of the European option price and the early exercise premium. The early exercise premium is calculated by the difference between the American and European option prices based on asymptotic expansions. The European option price is obtained by the efficient COS method. Based on the obtained American option price, the double Heston model is calibrated by minimizing the distance between model and market prices, which yields an optimization problem that is solved by a differential evolution algorithm combined with the Matlab function fmincon.m. Numerical results show that the pricing approach is fast and accurate. Empirical results show that the double Heston model has better performance in pricing short-maturity American put options and capturing the volatility term structure of American put options than the Heston model.     

The Jacobi stochastic volatility model

We introduce a novel stochastic volatility model where the squared volatility of the asset return follows a Jacobi process. It contains the Heston model as a limit case. We show that the joint density of any finite sequence of log-returns admits a Gram–Charlier A expansion with closed-form coefficients. We derive closed-form series representations for option prices whose discounted payoffs are functions of the asset price trajectory at finitely many time points. This includes European call, put and digital options, forward start options, and can be applied to discretely monitored Asian options. In a numerical study, we show that option prices can be accurately and efficiently approximated by truncating their series representations.     

A multivariate stochastic volatility model with applications in the foreign exchange market

The main objective of this paper is to study the behavior of a daily calibration of a multivariate stochastic volatility model, namely the principal component stochastic volatility (PCSV) model, to market data of plain vanilla options on foreign exchange rates. To this end, a general setting describing a foreign exchange market is introduced. Two adequate models—PCSV and a simpler multivariate Heston model—are adjusted to suit the foreign exchange setting. For both models, characteristic functions are found which allow for an almost instantaneous calculation of option prices using Fourier techniques. After presenting the general calibration procedure, both the multivariate Heston and the PCSV models are calibrated to a time series of option data on three exchange rates—USD-SEK, EUR-SEK, and EUR-USD—spanning more than 11 years. Finally, the benefits of the PCSV model which we find to be superior to the multivariate extension of the Heston model in replicating the dynamics of these options are highlighted.     

Maximum likelihood estimation of stochastic volatility models

We develop and implement a method for maximum likelihood estimation in closed-form of stochastic volatility models. Using Monte Carlo simulations, we compare a full likelihood procedure, where an option price is inverted into the unobservable volatility state, to an approximate likelihood procedure where the volatility state is replaced by proxies based on the implied volatility of a short-dated at-the-money option. The approximation results in a small loss of accuracy relative to the standard errors due to sampling noise. We apply this method to market prices of index options for several stochastic volatility models, and compare the characteristics of the estimated models. The evidence for a general CEV model, which nests both the affine Heston model and a GARCH model, suggests that the elasticity of variance of volatility lies between that assumed by the two nested models.     

Polynomial processes and their applications to mathematical finance

We introduce a class of Markov processes, called m-polynomial, for which the calculation of (mixed) moments up to order m only requires the computation of matrix exponentials. This class contains affine processes, processes with quadratic diffusion coefficients, as well as Lévy-driven SDEs with affine vector fields. Thus, many popular models such as exponential Lévy models or affine models are covered by this setting. The applications range from statistical GMM estimation procedures to new techniques for option pricing and hedging. For instance, the efficient and easy computation of moments can be used for variance reduction techniques in Monte Carlo methods.     

Unspanned stochastic volatility and fixed income derivatives pricing

We propose a parsimonious ‘unspanned stochastic volatility’ model of the term structure and study its implications for fixed-income option prices. The drift and quadratic variation of the short rate are affine in three state variables (the short rate, its long-term mean and variance) which follow a joint Markov (vector) process. Yet, bond prices are exponential affine functions of only two state variables, independent of the current interest rate volatility level. Because this result holds for an arbitrary volatility process, such a process can be calibrated to match fixed income derivative prices. Furthermore, this model can be ‘extended’ (by relaxing the time-homogeneity) to fit any arbitrary term structure. In its ‘HJM’ form, this model nests the analogous stochastic equity volatility model of Heston (1993) [Heston, S.L., 1993. A closed form solution for options with stochastic volatility. Review of Financial Studies 6, 327–343]. In particular, if the volatility process is specified to be affine, closed-form solutions for interest rate options obtain. We propose an efficient algorithm to compute these prices. An application using data on caps and floors shows that the model can capture very well the implied Black spot volatility surface, while simultaneously fitting the observed term structure.     

Risk-Neutral Parameter Shifts and Derivatives Pricing in Discrete Time

We obtain a large class of discrete‐time risk‐neutral valuation relationships, or “preference‐free” derivatives pricing models, by imposing a simple restriction on the state‐price density process. The risk‐neutral stock‐return and forward‐rate dynamics are obtained by changing only a location parameter, which can be determined independent of the preference and true location parameters. The Gaussian models of Rubinstein (1976) , Brennan (1979) , and Câmera (2003) , and the gamma model of Heston (1993) are all special cases. The model provides simple relationships between expected returns and state‐price density parameters analogous to the diffusion case.     

Option pricing in a Garch model with tempered stable innovations

The key problem for option pricing in Garch models is that the risk-neutral distribution of the underlying at maturity is unknown. Heston and Nandi solved this problem by computing the characteristic function of the underlying by a recursive procedure. Following the same idea, Christoffersen, Heston and Jacobs proposed a Garch-like model with inverse Gaussian innovations and recently Bellini and Mercuri obtained a similar procedure in a model with Gamma innovations. We present a model with tempered stable innovations that encompasses both the CHJ and the BM models as special cases. The proposed model is calibrated on S&P500 closing option prices and its performance is compared with the CHJ, the BM and the Heston–Nandi models.     

Variance swap dynamics

We compare several parametric and non-parametric approaches for modelling variance swap curves by conducting an in-sample and an out-of-sample analysis using market prices. The forecasted Heston model gives the best overall performance. Moreover, the static Heston model highlights some problems of stochastic volatility models in option pricing of forward starting products.     

A comparison principle between rough and non-rough Heston models—with applications to the volatility surface

We present a number of related comparison results, which allow one to compare moment explosion times, moment generating functions and critical moments between rough and non-rough Heston models of stochastic volatility. All results are based on a comparison principle for certain non-linear Volterra integral equations. Our upper bound for the moment explosion time is different from the bound introduced by Gerhold, Gerstenecker and Pinter [Moment explosions in the rough Heston model. Decisions in Economics and Finance, 2019, 42, 575–608] and tighter for typical parameter values. The results can be directly transferred to a comparison principle for the asymptotic slope of implied variance between rough and non-rough Heston models. This principle shows that the ratio of implied variance slopes in the rough versus non-rough Heston model increases at least with power-law behavior for small maturities.     

International money and stock market contingent claims

We develop a unified approach with closed-form solutions for pricing bonds, stocks, currencies and their derivatives. The specification assumes a fundamental risk factor represented by a stochastic positive definite matrix following a Wishart autoregressive (WAR) process. By assuming a volatility-in-mean specification for the domestic stock returns and the relative changes of exchange rates, and a domestic stochastic discount factor exponential affine with respect to the fundamental risk, it is possible to derive closed form solutions for the term structures of interest rates and for the risk-neutral probabilities while keeping the flexibility of the model. In particular:

i) The domestic and foreign term structures are jointly affine and correspond to Wishart quadratic term structures, which can ensure the positivity of interest rates;

ii) In this framework where the stock price follows a model with stochastic volatility, we obtain explicit or quasi-explicit formulas for futures and forward contracts, swaps and options. This extends results by

Heston (1993)

and

Ball and Roma (1994)

.

A variance reduction technique based on integral representations

Abstract

Standard Monte Carlo methods can often be significantly improved with the addition of appropriate variance reduction techniques. In this paper a new and powerful variance reduction technique is presented. The method is based directly on the Itô calculus and is used to find unbiased variance-reduced estimators for the expectation of functionals of Itô diffusion processes. The approach considered has wide applicability: for instance, it can be used as a means of approximating solutions of parabolic partial differential equations or applied to valuation problems that arise in mathematical finance. We illustrate how the method can be applied by considering the pricing of European-style derivative securities for a class of stochastic volatility models, including the Heston model.     

Solvable local and stochastic volatility models: supersymmetric methods in option pricing

In this paper we provide an extensive classification of one- and two-dimensional diffusion processes which admit an exact solution to the Kolmogorov (and hence Black–Scholes) equation (in terms of hypergeometric functions). By identifying the one-dimensional solvable processes with the class of integrable superpotentials introduced recently in supersymmetric quantum mechanics, we obtain new analytical solutions. In particular, by applying supersymmetric transformations on a known solvable diffusion process (such as the Natanzon process for which the solution is given by a hypergeometric function), we obtain a hierarchy of new solutions. These solutions are given by a sum of hypergeometric functions, generalizing the results obtained in a paper by Albanese et al. (Albanese, C., Campolieti, G., Carr, P. and Lipton, A., Black–Scholes goes hypergeometric. Risk Mag., 2001, 14, 99–103). For two-dimensional processes, more precisely stochastic volatility models, the classification is achieved for a specific class called gauge-free models including the Heston model, the 3?/?2-model and the geometric Brownian model. We then present a new exact stochastic volatility model belonging to this class.     

Probability distribution of returns in the Heston model with stochastic volatility*

Abstract

We study the Heston model, where the stock price dynamics is governed by a geometrical (multiplicative) Brownian motion with stochastic variance. We solve the corresponding Fokker‐Planck equation exactly and, after integrating out the variance, find an analytic formula for the time‐dependent probability distribution of stock price changes (returns). The formula is in excellent agreement with the Dow‐Jones index for time lags from 1 to 250 trading days. For large returns, the distribution is exponential in log‐returns with a time‐dependent exponent, whereas for small returns it is Gaussian. For time lags longer than the relaxation time of variance, the probability distribution can be expressed in a scaling form using a Bessel function. The Dow‐Jones data for 1982–2001 follow the scaling function for seven orders of magnitude.     

Regime switching affine processes with applications to finance

We introduce the notion of a regime switching affine process. Informally this is a Markov process that behaves conditionally on each regime as an affine process with specific parameters. To facilitate our analysis, specific restrictions are imposed on these parameters. The regime switches are driven by a Markov chain. We prove that the joint process of the Markov chain and the conditionally affine part is a process with an affine structure on an enlarged state space, conditionally on the starting state of the Markov chain. Like for affine processes, the characteristic function can be expressed in a set of ordinary differential equations that can sometimes be solved analytically. This result unifies several semi-analytical solutions found in the literature for pricing derivatives of specific regime switching processes on smaller state spaces. It also provides a unifying theory that allows us to introduce regime switching to the pricing of many derivatives within the broad class of affine processes. Examples include European options and term structure derivatives with stochastic volatility and default. Essentially, whenever there is a pricing solution based on an affine process, we can extend this to a regime switching affine process without sacrificing the analytical tractability of the affine process.