Stochastic volatility options pricing with wavelets and artificial neural networks

Abstract

The paper describes an alternative options pricing method which uses a binomial tree linked to an innovative stochastic volatility model. The volatility model is based on wavelets and artificial neural networks. Wavelets provide a convenient signal/noise decomposition of the volatility in the nonlinear feature space. Neural networks are used to infer future volatility from the wavelets feature space in an iterative manner. The bootstrap method provides the 95% confidence intervals for the options prices. Market options prices as quoted on the Chicago Board Options Exchange are used for performance comparison between the Black‐Scholes model and a new options pricing scheme. The proposed dynamic volatility model produces as good as and often better options prices than the conventional Black‐Scholes formulae.     

Maximum likelihood estimation of non-affine volatility processes

In this paper we develop a new estimation method for extracting non-affine latent stochastic volatility and risk premia from measures of model-free realized and risk-neutral integrated volatility. We estimate non-affine models with nonlinear drift and constant elasticity of variance and we compare them to the popular square-root stochastic volatility model. Our empirical findings are: (1) the square-root model is misspecified; (2) the inclusion of constant elasticity of variance and nonlinear drift captures stylized facts of volatility dynamics and (3) the square-root stochastic volatility model is explosive under the risk-neutral probability measure.     

On VIX futures in the rough Bergomi model

The rough Bergomi model introduced by Bayer et al. [Quant. Finance, 2015, 1–18] has been outperforming conventional Markovian stochastic volatility models by reproducing implied volatility smiles in a very realistic manner, in particular for short maturities. We investigate here the dynamics of the VIX and the forward variance curve generated by this model, and develop efficient pricing algorithms for VIX futures and options. We further analyse the validity of the rough Bergomi model to jointly describe the VIX and the SPX, and present a joint calibration algorithm based on the hybrid scheme by Bennedsen et al. [Finance Stoch., forthcoming].     

Affine fractional stochastic volatility models

By fractional integration of a square root volatility process, we propose in this paper a long memory extension of the Heston (Rev Financ Stud 6:327–343, 1993) option pricing model. Long memory in the volatility process allows us to explain some option pricing puzzles as steep volatility smiles in long term options and co-movements between implied and realized volatility. Moreover, we take advantage of the analytical tractability of affine diffusion models to clearly disentangle long term components and short term variations in the term structure of volatility smiles. In addition, we provide a recursive algorithm of discretization of fractional integrals in order to be able to implement a method of moments based estimation procedure from the high frequency observation of realized volatilities.     

Orthogonal expansions for VIX options under affine jump diffusions

In this work we derive new closed-form pricing formulas for VIX options in the jump-diffusion SVJJ model proposed by Duffie et al. [Econometrica, 2000, 68, 1343–1376]. Our approach is based on the classic methodology of approximating a density function with an orthogonal expansion of polynomials weighted by a kernel. Orthogonal expansions based on the Gaussian distribution, such as Edgeworth or Gram–Charlier expansions, have been successfully employed by a number of authors in the context of equity options. However, these expansions are not quite suitable for volatility or variance densities as they inherently assign positive mass to the negative real line. Here we approximate option prices via expansions that instead are based on kernels defined on the positive real line. Specifically, we consider a flexible family of distributions, which generalizes the gamma kernel associated with the classic Laguerre expansions. The method can be employed whenever the moments of the underlying variance distribution are known. It provides fast and accurate price computations, and therefore it represents a valid and possibly more robust alternative to pricing techniques based on Fourier transform inversions.     

Keep on smiling? The pricing of Quanto options when all covariances are stochastic

The paper introduces a model for the joint dynamics of asset prices which can capture both a stochastic correlation between stock returns as well as between stock returns and volatilities (stochastic leverage). By relying on two factors for stochastic volatility, the model allows for stochastic leverage and is thus able to explain time-varying slopes of the smiles. The use of Wishart processes for the covariance matrix of returns enables the model to also capture stochastic correlations between the assets. Our model offers an integrated pricing approach for both Quanto and plain-vanilla options on the stock as well as the foreign exchange rate. We derive semi-closed form solutions for option prices and analyze the impact of state variables. Quanto options offer a significant exposure to the stochastic covariance between stock prices and exchange rates. In contrast to standard models, the smile of stock options, the smile of currency options, and the price differences between Quanto options and plain-vanilla options can change independently of each other.     

Pricing variance and volatility swaps in a stochastic volatility model with regime switching: discrete observations case

This study presents a set of closed-form exact solutions for pricing discretely sampled variance swaps and volatility swaps, based on the Heston stochastic volatility model with regime switching. In comparison with all the previous studies in the literature, this research, which obtains closed-form exact solutions for variance and volatility swaps with discrete sampling times, serves several purposes. (1) It verifies the degree of validity of Elliott et al.'s [Appl. Math. Finance, 2007, 14(1), 41–62] continuous-sampling-time approximation for variance and volatility swaps of relatively short sampling periods. (2) It examines the effect of ignoring regime switching on pricing variance and volatility swaps. (3) It contributes to bridging the gap between Zhu and Lian's [Math. Finance, 2011, 21(2), 233–256] approach and Elliott et al.'s framework. (4) Finally, it presents a semi-Monte-Carlo simulation for the pricing of other important realized variance based derivatives.     

A Comparison of Option Prices Under Different Pricing Measures in a Stochastic Volatility Model with Correlation

This paper investigates option prices in an incomplete stochastic volatility model with correlation. In a general setting, we prove an ordering result which says that prices for European options with convex payoffs are decreasing in the market price of volatility risk.As an example, and as our main motivation, we investigate option pricing under the class of q-optimal pricing measures. The q-optimal pricing measure is related to the marginal utility indifference price of an agent with constant relative risk aversion. Using the ordering result, we prove comparison theorems between option prices under the minimal martingale, minimal entropy and variance-optimal pricing measures. If the Sharpe ratio is deterministic, the comparison collapses to the well known result that option prices computed under these three pricing measures are the same.As a concrete example, we specialize to a variant of the Hull-White or Heston model for which the Sharpe ratio is increasing in volatility. For this example we are able to deduce option prices are decreasing in the parameter q. Numerical solution of the pricing pde corroborates the theory and shows the magnitude of the differences in option price due to varying q.JEL Classification: D52, G13     

Generic pricing of FX,inflation and stock options under stochastic interest rates and stochastic volatility

We consider the pricing of FX, inflation and stock options under stochastic interest rates and stochastic volatility, for which we use a generic multi-currency framework. We allow for a general correlation structure between the drivers of the volatility, the inflation index, the domestic (nominal) and the foreign (real) rates. Having the flexibility to correlate the underlying FX/inflation/stock index with both stochastic volatility and stochastic interest rates yields a realistic model that is of practical importance for the pricing and hedging of options with a long-term exposure. We derive explicit valuation formulas for various securities, such as vanilla call/put options, forward starting options, inflation-indexed swaps and inflation caps/floors. These vanilla derivatives can be valued in closed form under Schöbel and Zhu [Eur. Finance Rev., 1999, 4, 23–46] stochastic volatility, whereas we devise an (Monte Carlo) approximation in the form of a very effective control variate for the general Heston [Rev. Financial Stud., 1993, 6, 327–343] model. Finally, we investigate the quality of this approximation numerically and consider a calibration example to FX and inflation market data.     

The representation of American options prices under stochastic volatility and jump-diffusion dynamics

This paper considers the problem of pricing American options when the dynamics of the underlying are driven by both stochastic volatility following a square-root process as used by Heston [Rev. Financial Stud., 1993, 6, 327–343], and by a Poisson jump process as introduced by Merton [J. Financial Econ., 1976, 3, 125–144]. Probability arguments are invoked to find a representation of the solution in terms of expectations over the joint distribution of the underlying process. A combination of Fourier transform in the log stock price and Laplace transform in the volatility is then applied to find the transition probability density function of the underlying process. It turns out that the price is given by an integral dependent upon the early exercise surface, for which a corresponding integral equation is obtained. The solution generalizes in an intuitive way the structure of the solution to the corresponding European option pricing problem obtained by Scott [Math. Finance, 1997, 7(4), 413–426], but here in the case of a call option and constant interest rates.     

Implied integrated variance and hedging

If the volatility is stochastic, stock price returns and European option prices depend on the time average of the variance, i.e. the integrated variance, not on the path of the volatility. Applying a Bayesian statistical approach, we compute a forward-looking estimate of this variance, an option-implied integrated variance. Simultaneously, we obtain estimates of the correlation coefficient between stock price and volatility shocks, and of the parameters of the volatility process. Due to the convexity of the Black–Scholes formula with respect to the volatility, pricing and hedging with Black–Scholes-type formulas and the implied volatility often lead to inaccuracies if the volatility is stochastic. Theoretically, this problem can be avoided by using Hull–White-type option pricing and hedging formulas and the integrated variance. We use the implied integrated variance and Hull–White-type formulas to hedge European options and certain volatility derivatives.     

Are classical option pricing models consistent with observed option second-order moments? Evidence from high-frequency data

As a means of validating an option pricing model, we compare the ex-post intra-day realized variance of options with the realized variance of the associated underlying asset that would be implied using assumptions as in the Black and Scholes (BS) model, the Heston, and the Bates model. Based on data for the S&P 500 index, we find that the BS model is strongly directionally biased due to the presence of stochastic volatility. The Heston model reduces the mismatch in realized variance between the two markets, but deviations are still significant. With the exception of short-dated options, we achieve best approximations after controlling for the presence of jumps in the underlying dynamics. Finally, we provide evidence that, although heavily biased, the realized variance based on the BS model contains relevant predictive information that can be exploited when option high-frequency data is not available.     

An Asymptotic Expansion Approach to Currency Options with a Market Model of Interest Rates under Stochastic Volatility Processes of Spot Exchange Rates

This paper proposes an asymptotic expansion scheme of currency options with a libor market model of interest rates and stochastic volatility models of spot exchange rates. In particular, we derive closed-form approximation formulas for the density functions of the underlying assets and for pricing currency options based on a third order asymptotic expansion scheme; we do not model a foreign exchange rate’s variance such as in Heston [(1993) The Review of Financial studies, 6, 327–343], but its volatility that follows a general time-inhomogeneous Markovian process. Further, the correlations among all the factors such as domestic and foreign interest rates, a spot foreign exchange rate and its volatility, are allowed. Finally, numerical examples are provided and the pricing formula are applied to the calibration of volatility surfaces in the JPY/USD option market.     

Pricing double barrier options under a volatility regime-switching model with psychological barriers

The prices of lots of assets have been proved in literature to exhibit special behaviors around psychological barriers, which is an important fact needed to be considered when pricing derivatives. In this paper, we discuss the valuation problem of double barrier options under a volatility regime-switching model where there exist psychological barriers in the prices of underlying assets. The volatility can shift between two regimes, that is to say, when the asset price rises up or falls down through the psychological barrier, the volatility takes two different values. Using the Laplace transform approach, we obtain the price of the double barrier knock-out call option as well as its delta. We also provide the eigenfunction expansion pricing formula and examine the effect of the psychological barrier on the option price and delta, finding that the gamma of the option is discontinuous at such barriers.     

A numerical method to price exotic path-dependent options on an underlying described by the Heston stochastic volatility model

We consider the problem of pricing European exotic path-dependent derivatives on an underlying described by the Heston stochastic volatility model. Lipton has found a closed form integral representation of the joint transition probability density function of underlying price and variance in the Heston model. We give a convenient numerical approximation of this formula and we use the obtained approximated transition probability density function to price discrete path-dependent options as discounted expectations. The expected value of the payoff is calculated evaluating an integral with the Monte Carlo method using a variance reduction technique based on a suitable approximation of the transition probability density function of the Heston model. As a test case, we evaluate the price of a discrete arithmetic average Asian option, when the average over n = 12 prices is considered, that is when the integral to evaluate is a 2n = 24 dimensional integral. We show that the method proposed is computationally efficient and gives accurate results.     

How Integrated are Credit and Equity Markets? Evidence from Index Options

We study the extent to which credit index (CDX) options are priced consistent with S&P 500 (SPX) equity index options. We derive analytical expressions for CDX and SPX options within a structural credit-risk model with stochastic volatility and jumps using new results for pricing compound options via multivariate affine transform analysis. The model captures many aspects of the joint dynamics of CDX and SPX options. However, it cannot reconcile the relative levels of option prices, suggesting that credit and equity markets are not fully integrated. A strategy of selling CDX volatility yields significantly higher excess returns than selling SPX volatility.     

A transform approach to compute prices and Greeks of barrier options driven by a class of Lévy processes

In this paper we propose a transform method to compute the prices and Greeks of barrier options driven by a class of Lévy processes. We derive analytical expressions for the Laplace transforms in time of the prices and sensitivities of single barrier options in an exponential Lévy model with hyper-exponential jumps. Inversion of these single Laplace transforms yields rapid, accurate results. These results are employed to construct an approximation of the prices and sensitivities of barrier options in exponential generalized hyper-exponential Lévy models. The latter class includes many of the Lévy models employed in quantitative finance such as the variance gamma (VG), KoBoL, generalized hyperbolic, and the normal inverse Gaussian (NIG) models. Convergence of the approximating prices and sensitivities is proved. To provide a numerical illustration, this transform approach is compared with Monte Carlo simulation in cases where the driving process is a VG and a NIG Lévy process. Parameters are calibrated to Stoxx50E call options.     

Lévy processes driven by stochastic volatility

In this paper we extend option pricing under Lévy dynamics, by assuming that the volatility of the Lévy process is stochastic. We, therefore, develop the analog of the standard stochastic volatility models, when the underlying process is not a standard (unit variance) Brownian motion, but rather a standardized Lévy process. We present a methodology that allows one to compute option prices, under virtually any set of diffusive dynamics for the parameters of the volatility process. First, we use ‘local consistency’ arguments to approximate the volatility process with a finite, but sufficiently dense Markov chain; we then use this regime switching approximation to efficiently compute option prices using Fourier inversion. A detailed example, based on a generalization of the popular stochastic volatility model of Heston (Rev Financial Stud 6 (1993) 327), is used to illustrate the implementation of the algorithms. Computer code is available at www.theponytail.net/     

Estimation and pricing under long-memory stochastic volatility

We treat the problem of option pricing under a stochastic volatility model that exhibits long-range dependence. We model the price process as a Geometric Brownian Motion with volatility evolving as a fractional Ornstein–Uhlenbeck process. We assume that the model has long-memory, thus the memory parameter H in the volatility is greater than 0.5. Although the price process evolves in continuous time, the reality is that observations can only be collected in discrete time. Using historical stock price information we adapt an interacting particle stochastic filtering algorithm to estimate the stochastic volatility empirical distribution. In order to deal with the pricing problem we construct a multinomial recombining tree using sampled values of the volatility from the stochastic volatility empirical measure. Moreover, we describe how to estimate the parameters of our model, including the long-memory parameter of the fractional Brownian motion that drives the volatility process using an implied method. Finally, we compute option prices on the S&P 500 index and we compare our estimated prices with the market option prices.     

Implied volatilities,stochastic interest rates,and currency futures options valuation: an empirical investigation

Different models of pricing currency call and put options on futures are empirically tested. Option prices are determined using different models and compared to actual market prices. Option prices are determined using historical as well as implied volatility. The different models tested include both constant and stochastic interest rate models. To determine if the model prices are different from the market prices, regression analysis and paired t-tests are performed. To see which model misprices the least, root mean square errors are determined. It is found that better results are obtained when implied volatility is used. Stochastic interest rate models perform better than constant interest rate models.