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.     

Option pricing under a stressed-beta model

Empirical studies have concluded that stochastic volatility is an important component of option prices. We introduce a regime-switching mechanism into a continuous-time Capital Asset Pricing Model which naturally induces stochastic volatility in the asset price. Under this Stressed-Beta model, the mechanism is relatively simple: the slope coefficient—which measures asset returns relative to market returns—switches between two values, depending on the market being above or below a given level. After specifying the model, we use it to price European options on the asset. Interestingly, these option prices are given explicitly as integrals with respect to known densities. We find that the model is able to produce a volatility skew, which is a prominent feature in option markets. This opens the possibility of forward-looking calibration of the slope coefficients, using option data, as illustrated in the paper.     

Option market making under inventory risk

We propose a mean-variance framework to analyze the optimal quoting policy of an option market maker. The market maker’s profits come from the bid-ask spreads received over the course of a trading day, while the risk comes from uncertainty in the value of his portfolio, or inventory. Within this framework, we study the impact of liquidity and market incompleteness on the optimal bid and ask prices of the option. First, we consider a market maker in a complete market, where continuous trading in a perfectly liquid underlying stock is allowed. In this setting, the market maker may remove all risk by Delta hedging, and the optimal quotes will depend on the option’s liquidity, but not on the inventory. Second, we model a market maker who may not trade continuously in the underlying stock, but rather sets bid and ask quotes in the option and this illiquid stock. We find that the optimal stock and option quotes depend on the relative liquidity of both instruments as well as on the net Delta of the inventory. Third, we consider an incomplete market with residual risks due to stochastic volatility and large overnight moves in the stock price. In this setting, the optimal quotes depend on the liquidity of the option and on the net Vega and Gamma of the inventory.      

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.     

Asymmetric information about volatility: How does it affect implied volatility,option prices and market liquidity?

This paper develops a model of asymmetric information in which an investor has information regarding the future volatility of the price process of an asset and trades an option on the asset. The model relates the level and curvature of the smile in implied volatilities as well as mispricing by the Black-Scholes model to net options order flows (to the market maker). It is found that an increase in net options order flows (to the market maker) increases the level of implied volatilities and results in greater mispricing by the Black-Scholes model, besides impacting the curvature of the smile. The liquidity of the option market is found to be decreasing in the amount of uncertainty about future volatility that is consistent with existing evidence. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.     

Realizing smiles: Options pricing with realized volatility

We develop a discrete-time stochastic volatility option pricing model exploiting the information contained in the Realized Volatility (RV), which is used as a proxy of the unobservable log-return volatility. We model the RV dynamics by a simple and effective long-memory process, whose parameters can be easily estimated using historical data. Assuming an exponentially affine stochastic discount factor, we obtain a fully analytic change of measure. An empirical analysis of Standard and Poor's 500 index options illustrates that our model outperforms competing time-varying and stochastic volatility option pricing models.     

Portfolio optimization for student t and skewed t returns

We consider the problem of valuing a European option written on an asset whose dynamics are described by an exponential Lévy-type model. In our framework, both the volatility and jump-intensity are allowed to vary stochastically in time through common driving factors—one fast-varying and one slow-varying. Using Fourier analysis we derive an explicit formula for the approximate price of any European-style derivative whose payoff has a generalized Fourier transform; in particular, this includes European calls and puts. From a theoretical perspective, our results extend the class of multiscale stochastic volatility models of Fouque et al. [Multiscale Stochastic Volatility for Equity, Interest Rate, and Credit Derivatives, 2011] to models of the exponential Lévy type. From a financial perspective, the inclusion of jumps and stochastic volatility allow us to capture the term-structure of implied volatility, as demonstrated in a calibration to S&;P500 options data.     

Understanding Delta-Hedged Option Returns in Stochastic Volatility Environments

In this paper, we provide a novel representation of delta-hedged option returns in a stochastic volatility environment. The representation of delta-hedged option returns provided in this paper consists of two terms: volatility risk premium and parameter estimation risk. In an empirical analysis, we examine delta-hedged option returns based on the result of a historical simulation with the USD-JPY currency option market data from October 2003 to June 2010. We find that the delta-hedged option returns for OTM put options are strongly affected by parameter estimation risk as well as the volatility risk premium, especially in the post-Lehman shock period.     

A chaos expansion approach under hybrid volatility models

In this paper, we propose an approximation method based on the Wiener–Ito chaos expansion for the pricing of European contingent claims. Our method is applicable to widely used option pricing models such as local volatility models, stochastic volatility models and their combinations. This method is useful in practice since the resulting approximation formula is not computationally expensive, hence it is suitable for calibration purposes. We will show through some numerical examples that our approximation remains quite good even for the long maturity and/or the high volatility cases, which is a desired feature. As an example, we propose a hybrid volatility model and apply our approximation formula to the JPY/USD currency option market obtaining very accurate results.     

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.     

Testing the local volatility assumption: a statistical approach

In practice, the choice of using a local volatility model or a stochastic volatility model is made according to their respective ability to fit implied volatility surfaces. In this paper, we adopt a different point of view. Indeed, using a purely statistical methodology, we design new procedures aiming at testing the assumption of a local volatility model for the price dynamics, against the alternative of a stochastic volatility model. These test procedures are based only on historical data and do not require any calibration procedures via option prices. We also provide a convincing simulation study and an empirical analysis on future contracts on interest rates.     

Displaced lognormal volatility skews: analysis and applications to stochastic volatility simulations

We analyze the implied volatility skews generated by displaced lognormal diffusions. In particular, we prove the global monotonicity of implied volatility, and an at-the-money bound on the steepness of downward volatility skews, under displaced lognormal dynamics, which therefore cannot reproduce some features observed in equity markets. A variant, the displaced anti-lognormal, overcomes the steepness constraint, but its state space is bounded above and unbounded below. In light of these limitations on what features the displaced (anti-)lognormal (DL) can model, we exploit the DL, not as a model, but as a control variate, to reduce variance in Monte Carlo simulations of the CEV and SABR local/stochastic volatility models. For either use—as model, or as control variate—the DL’s parameters require estimation. We find an explicit formula for the DL’s short-expiry limiting volatility skew, which allows direct calibration of its parameters to volatility skews implied by market data or by other models.     

Hedging volatility risk

Volatility risk plays an important role in the management of portfolios of derivative assets as well as portfolios of basic assets. This risk is currently managed by volatility “swaps” or futures. However, this risk could be managed more efficiently using options on volatility that were proposed in the past but were never introduced mainly due to the lack of a cost efficient tradable underlying asset.The objective of this paper is to introduce a new volatility instrument, an option on a straddle, which can be used to hedge volatility risk. The design and valuation of such an instrument are the basic ingredients of a successful financial product. In order to value these options, we combine the approaches of compound options and stochastic volatility. Our numerical results show that the straddle option is a powerful instrument to hedge volatility risk. An additional benefit of such an innovation is that it will provide a direct estimate of the market price for volatility risk.     

Modeling the bid/ask spread: measuring the inventory-holding premium

The need to understand and measure the determinants of market maker bid/ask spreads is crucial in evaluating the merits of competing market structures and the fairness of market maker rents. This study develops a simple, parsimonious model for the market maker's spread that accounts for the effects of price discreteness induced by minimum tick size, order-processing costs, inventory-holding costs, adverse selection, and competition. The inventory-holding and adverse selection cost components of spread are modeled as an option with a stochastic time to expiration. This inventory-holding premium embedded in the spread represents compensation for the price risk borne by the market maker while the security is held in inventory. The premium is partitioned in such a way that the inventory-holding and adverse selection cost components, as well as the probability of an informed trade, are identified. The model is tested empirically using Nasdaq stocks in three distinct minimum tick size regimes and is shown to perform well both in an absolute sense and relative to competing specifications.     

Maturity cycles in implied volatility

The skew effect in market implied volatility can be reproduced by option pricing theory based on stochastic volatility models for the price of the underlying asset. Here we study the performance of the calibration of the S&P 500 implied volatility surface using the asymptotic pricing theory under fast mean-reverting stochastic volatility described in [8]. The time-variation of the fitted skew-slope parameter shows a periodic behaviour that depends on the option maturity dates in the future, which are known in advance. By extending the mathematical analysis to incorporate model parameters which are time-varying, we show this behaviour can be explained in a manner consistent with a large model class for the underlying price dynamics with time-periodic volatility coefficients.Received: December 2003, Mathematics Subject Classification (2000): 91B70, 60F05, 60H30JEL Classification: C13, G13Jean-Pierre Fouque: Work partially supported by NSF grant DMS-0071744.Ronnie Sircar: Work supported by NSF grant DMS-0090067. We are grateful to Peter Thurston for research assistance.We thank a referee for his/her comments which improved the paper.     

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/     

Optimal portfolios and Heston's stochastic volatility model: an explicit solution for power utility

Given an investor maximizing utility from terminal wealth with respect to a power utility function, we present a verification result for portfolio problems with stochastic volatility. Applying this result, we solve the portfolio problem for Heston's stochastic volatility model. We find that only under a specific condition on the model parameters does the problem possess a unique solution leading to a partial equilibrium. Finally, it is demonstrated that the results critically hinge upon the specification of the market price of risk. We conclude that, in applications, one has to be very careful when exogenously specifying the form of the market price of risk.     

Application of Heston’s Model to the Chinese Stock Market

This article applies Heston’s (1993) stochastic volatility model to the Chinese stock market indices and subsequently assesses its pricing performance. A two-step estimation procedure is adopted to calibrate Heston’s model. First, we find that the option price is affected by both the moneyness and the maturity. Second, Heston’s model is more likely to overprice options, whereas the BS model tends to underestimate options. Finally, Heston’s model, by employing volatility as a random process, significantly improves the pricing accuracy compared to the BS model. Therefore, Heston’s model is tractable to analyze the Chinese stock market indices, and there is volatility risk that must not be overlooked in the Chinese stock market.     

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.