Statistical estimation of Lévy-type stochastic volatility models

Volatility clustering and leverage are two of the most prominent stylized features of the dynamics of asset prices. In order to incorporate these features as well as the typical fat-tails of the log return distributions, several types of exponential Lévy models driven by random clocks have been proposed in the literature. These models constitute a viable alternative to the classical stochastic volatility approach based on SDEs driven by Wiener processes. This paper has two main objectives. First, using threshold type estimators based on high-frequency discrete observations of the process, we consider the recovery problem of the underlying random clock of the process. We show consistency of our estimator in the mean-square sense, extending former results in the literature for more general Lévy processes and for irregular sampling schemes. Secondly, we illustrate empirically the estimation of the random clock, the Blumenthal-Geetor index of jump activity, and the spectral Lévy measure of the process using real intraday high-frequency data.     

Calibrating the exponential Ornstein–Uhlenbeck multiscale stochastic volatility model

This paper demonstrates a tractable and efficient way of calibrating a multiscale exponential Ornstein–Uhlenbeck stochastic volatility model including a correlation between the asset return and its volatility. As opposed to many contributions where this correlation is assumed to be null, this framework allows one to describe the leverage effect widely observed in equity markets. The resulting model is non-exponential and driven by a degenerate noise, thus requiring a high level of care in designing the estimation algorithm. The way this difficulty is overcome provides guidelines concerning the development of an estimation algorithm in a non-standard framework. The authors propose using a block-type expectation maximization algorithm along with particle smoothing. This method results in an accurate calibration process able to identify up to three timescale factors. Furthermore, we introduce an intuitive heuristic which can be used to choose the number of factors.     

Multiple time scales and the exponential Ornstein–Uhlenbeck stochastic volatility model

We study the exponential Ornstein–Uhlenbeck stochastic volatility model and observe that the model shows a multiscale behaviour in the volatility autocorrelation. It also exhibits a leverage correlation and a probability profile for the stationary volatility which are consistent with market observations. All these features make the model quite appealing since it appears to be more complete than other stochastic volatility models also based on a two-dimensional diffusion. We finally present an approximate solution for the return probability density designed to capture the kurtosis and skewness effects.     

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/     

Comparison results for stochastic volatility models via coupling

The aim of this paper is to investigate the properties of stochastic volatility models, and to discuss to what extent, and with regard to which models, properties of the classical exponential Brownian motion model carry over to a stochastic volatility setting. The properties of the classical model of interest include the fact that the discounted stock price is positive for all t but converges to zero almost surely, the fact that it is a martingale but not a uniformly integrable martingale, and the fact that European option prices (with convex payoff functions) are convex in the initial stock price and increasing in volatility. We explain why these properties are significant economically, and give examples of stochastic volatility models where these properties continue to hold, and other examples where they fail. The main tool is a construction of a time-homogeneous autonomous volatility model via a time-change.     

Short-term asymptotics for the implied volatility skew under a stochastic volatility model with Lévy jumps

The implied volatility skew has received relatively little attention in the literature on short-term asymptotics for financial models with jumps, despite its importance in model selection and calibration. We rectify this by providing high order asymptotic expansions for the at-the-money implied volatility skew, under a rich class of stochastic volatility models with independent stable-like jumps of infinite variation. The case of a pure-jump stable-like Lévy model is also considered under the minimal possible conditions for the resulting expansion to be well defined. Unlike recent results for “near-the-money” option prices and implied volatility, the results herein aid in understanding how the implied volatility smile near expiry is affected by important features of the continuous component, such as the leverage and vol-of-vol parameters. As intermediary results, we obtain high order expansions for at-the-money digital call option prices, which furthermore allow us to infer analogous results for the delta of at-the-money options. Simulation results indicate that our asymptotic expansions give good fits for options with maturities up to one month, underpinning their relevance in practical applications, and an analysis of the implied volatility skew in recent S&P 500 options data shows it to be consistent with the infinite variation jump component of our models.     

Analysis of variance based instruments for Ornstein–Uhlenbeck type models: swap and price index

In this paper a couple of variance dependent instruments in the financial market are studied. Firstly, a number of aspects of the variance swap in connection to the Barndorff-Nielsen and Shephard model are studied. A partial integro-differential equation that describes the dynamics of the arbitrage-free price of the variance swap is formulated. Under appropriate assumptions for the first four cumulants of the driving subordinator, a Ve?e?-type theorem is proved. The bounds of the arbitrage-free variance swap price are also found. Finally, a price-weighted index modulated by market variance is introduced. The large-basket limit dynamics of the price index and the “error term” are derived. Empirical data driven numerical examples are provided in support of the proposed price index.     

Option pricing under stochastic volatility and tempered stable Lévy jumps

The purpose of this paper is to introduce a stochastic volatility model for option pricing that exhibits Lévy jump behavior. For this model, we derive the general formula for a European call option. A well known particular case of this class of models is the Bates model, for which the jumps are modeled by a compound Poisson process with normally distributed jumps. Alternatively, we turn our attention to infinite activity jumps produced by a tempered stable process. Then we empirically compare the estimated log-return probability density and the option prices produced from this model to both the Bates model and the Black–Scholes model. We find that the tempered stable jumps describe more precisely market prices than compound Poisson jumps assumed in the Bates model.     

Well-posedness and invariant measures for HJM models with deterministic volatility and Lévy noise

We give sufficient conditions for the existence, uniqueness and ergodicity of invariant measures for Musiela's stochastic partial differential equation with deterministic volatility and a Hilbert space valued driving Lévy noise. Conditions for the absence of arbitrage and for the existence of mild solutions are also discussed.     

Consumption-investment problem with transaction costs for Lévy-driven price processes

We consider an optimal control problem for a linear stochastic integro-differential equation with conic constraints on the phase variable and with the control of singular–regular type. Our setting includes consumption-investment problems for models of financial markets in the presence of proportional transaction costs, where the prices of the assets are given by a geometric Lévy process, and the investor is allowed to take short positions. We prove that the Bellman function of the problem is a viscosity solution of an HJB equation. A uniqueness theorem for the solution of the latter is established. Special attention is paid to the dynamic programming principle.     

Option pricing and hedging under a stochastic volatility Lévy process model

In this paper, we discuss a stochastic volatility model with a Lévy driving process and then apply the model to option pricing and hedging. The stochastic volatility in our model is defined by the continuous Markov chain. The risk-neutral measure is obtained by applying the Esscher transform. The option price using this model is computed by the Fourier transform method. We obtain the closed-form solution for the hedge ratio by applying locally risk-minimizing hedging.     

Multilevel Monte Carlo for exponential Lévy models

We apply the multilevel Monte Carlo method for option pricing problems using exponential Lévy models with a uniform timestep discretisation. For lookback and barrier options, we derive estimates of the convergence rate of the error introduced by the discrete monitoring of the running supremum of a broad class of Lévy processes. We then use these to obtain upper bounds on the multilevel Monte Carlo variance convergence rate for the variance gamma, NIG and \(\alpha\)-stable processes. We also provide an analysis of a trapezoidal approximation for Asian options. Our method is illustrated by numerical experiments.     

Leveraged Lévy processes as models for stock prices

Adopting a constant elasticity of variance formulation in the context of a general Lévy process as the driving uncertainty we show that the presence of the leverage effect? ?One explanation of the documented negative relation between market volatilities and the level of asset prices (the ‘smile’ or ‘skew’), we term the ‘leverage effect’, argues that this negative relation reflects greater risk taking by the management, induced by a fall in the asset price, with a view of maximizing the option value of equity shareholders. in this form has the implication that asset price processes satisfy a scaling hypothesis. We develop forward partial integro-differential equations under a general Markovian setup, and show in two examples (both continuous and pure-jump Lévy) how to use them for option pricing when stock prices follow our leveraged Lévy processes. Using calibrated models we then show an example of simulation-based pricing and report on the adequacy of using leveraged Lévy models to value equity structured products.     

The maximum surplus before ruin for dependent risk models through Farlie–Gumbel–Morgenstern copula

We extend the classical compound Poisson risk model to consider the distribution of the maximum surplus before ruin where the claim sizes depend on inter-claim times via the Farlie–Gumbel–Morgenstern copula. We derive an integro-differential equation with certain boundary conditions for this distribution, of which the Laplace transform is provided. We obtain the renewal equation and explicit expressions for this distribution are derived when the claim amounts are exponentially distributed. Finally, we present numerical examples.     

Have the GIPSI settled down? Breaks and multivariate stochastic volatility models for,and not against,the European financial integration

We investigate the integration of the European peripheral financial markets with Germany, France, and the UK using a combination of tests for structural breaks and return correlations derived from several multivariate stochastic volatility models. Our findings suggest that financial integration intensified in anticipation of the Euro, further strengthened by the EMU inception, and amplified in response to the 2007/2008 financial crisis. Hence, no evidence is found of decoupling of the equity markets in more troubled European countries from the core. Interestingly, the UK, despite staying outside the EMU, is not worse integrated with the GIPSI than Germany or France.     

Esscher transforms and the minimal entropy martingale measure for exponential Lévy models

In this paper we offer a systematic survey and comparison of the Esscher martingale transform for linear processes, the Esscher martingale transform for exponential processes, and the minimal entropy martingale measure for exponential Lévy models, and present some new results in order to give a complete characterization of those classes of measures. We illustrate the results with several concrete examples in detail.     

Comparing alternative Lévy base correlation models for pricing and hedging CDO tranches

In this paper we investigate alternative Lévy base correlation models that arise from the Gamma, Inverse Gaussian and CMY distribution classes. We compare these models with the basic (exponential) Lévy base correlation model and the classical Gaussian base correlation model. For all investigated models, the Lévy base correlation curve is significantly flatter than the corresponding Gaussian curve, which indicates better correspondence of the Lévy models with reality. Furthermore, we present the results of pricing bespoke tranchlets and comparing deltas of both standard and custom-made tranches under all the considered models. We focus on deltas with respect to the CDS index and individual CDSs, and the hedge ratio for hedging the equity tranche with the junior mezzanine.     

A Feynman–Kac-type formula for Lévy processes with discontinuous killing rates

The challenge to fruitfully merge state-of-the-art techniques from mathematical finance and numerical analysis has inspired researchers to develop fast deterministic option pricing methods. As a result, highly efficient algorithms to compute option prices in Lévy models by solving partial integro-differential equations have been developed. In order to provide a solid mathematical foundation for these methods, we derive a Feynman–Kac representation of variational solutions to partial integro-differential equations that characterize conditional expectations of functionals of killed time-inhomogeneous Lévy processes. We allow a wide range of underlying stochastic processes, comprising processes with Brownian part as well as a broad class of pure jump processes such as generalized hyperbolic, multivariate normal inverse Gaussian, tempered stable, and \(\alpha\)-semistable Lévy processes. By virtue of our mild regularity assumptions as to the killing rate and the initial condition of the partial integro-differential equation, our results provide a rigorous basis for numerous applications in financial mathematics and in probability theory. We implement a Galerkin scheme to solve the corresponding pricing equation numerically and illustrate the effect of a killing rate.