Options Arbitrage in Imperfect Markets

Option valuation models are based on an arbitrage strategy—hedging the option against the underlying asset and rebalancing continuously until expiration—that is only possible in a frictionless market. This paper simulates the impact of market imperfections and other problems with the “standard” arbitrage trade, including uncertain volatility, transactions costs, indivisibilities, and rebalancing only at discrete intervals. We find that, in an actual market such as that for stock index options, the standard arbitrage is exposed to such large risk and transactions costs that it can only establish very wide bounds on equilibrium options prices. This has important implications for price determination in options markets, as well as for testing of valuation models.     

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.     

General equilibrium pricing of options on the market portfolio with discontinuous returns

When the price process for a long-lived asset is of a mixedjump-diffusion type, pricing of options on that asset by arbitrageis not possible if trading is allowed only in the underlyingasset and a risk-less bond. Using a general equilibrium framework,we derive and analyse option prices when the underlying assetis the market portfolio with discontinuous returns. The premiumfor the risk of jumps and the diffusions risk forms a significantpart of the prices of the options. In this economy, an attemptedreplication of call and put options by the Black-Scholes typeof trading strategies may require substantial infusion of fundswhen jumps occur. We study the cost and risk implications ofsuch dynamic hedging plans.     

Financial market models with Lévy processes and time-varying volatility

Asset management and pricing models require the proper modeling of the return distribution of financial assets. While the return distribution used in the traditional theories of asset pricing and portfolio selection is the normal distribution, numerous studies that have investigated the empirical behavior of asset returns in financial markets throughout the world reject the hypothesis that asset return distributions are normally distribution. Alternative models for describing return distributions have been proposed since the 1960s, with the strongest empirical and theoretical support being provided for the family of stable distributions (with the normal distribution being a special case of this distribution). Since the turn of the century, specific forms of the stable distribution have been proposed and tested that better fit the observed behavior of historical return distributions. More specifically, subclasses of the tempered stable distribution have been proposed. In this paper, we propose one such subclass of the tempered stable distribution which we refer to as the “KR distribution”. We empirically test this distribution as well as two other recently proposed subclasses of the tempered stable distribution: the Carr–Geman–Madan–Yor (CGMY) distribution and the modified tempered stable (MTS) distribution. The advantage of the KR distribution over the other two distributions is that it has more flexible tail parameters. For these three subclasses of the tempered stable distribution, which are infinitely divisible and have exponential moments for some neighborhood of zero, we generate the exponential Lévy market models induced from them. We then construct a new GARCH model with the infinitely divisible distributed innovation and three subclasses of that GARCH model that incorporates three observed properties of asset returns: volatility clustering, fat tails, and skewness. We formulate the algorithm to find the risk-neutral return processes for those GARCH models using the “change of measure” for the tempered stable distributions. To compare the performance of those exponential Lévy models and the GARCH models, we report the results of the parameters estimated for the S&P 500 index and investigate the out-of-sample forecasting performance for those GARCH models for the S&P 500 option prices.     

Generalized DEA model of fundamental analysis and its application to portfolio optimization

Fundamental analysis is used in asset selection for equity portfolio management. In this paper, a generalized data envelopment analysis (DEA) model is developed to analyze a firm’s financial statements over time in order to determine a relative financial strength indicator (RFSI) that is predictive of firm’s stock price returns. RFSI is based on maximizing the correlation between the DEA-based score of financial strength and the stock market performance. This maximization involves a difficult binary nonlinear program that requires iterative re-configuration of parameters of financial statements as inputs and outputs. We utilize a two-step heuristic algorithm that combines random sampling and local search optimization. The proposed approach is tested with 230 firms from various US technology-industries to determine optimized RFSI indicators for stock selection. Then, those selected stocks are used within portfolio optimization models to demonstrate the usefulness of the scheme for portfolio risk management.     

Recovering the real-world density and liquidity premia from option data

In this paper, we develop a methodology for simultaneous recovery of the real-world probability density and liquidity premia from observed S&P 500 index option prices. Assuming the existence of a numéraire portfolio for the US equity market, fair prices of derivatives under the benchmark approach can be obtained directly under the real-world measure. Under this modelling framework, there exists a direct link between observed call option prices on the index and the real-world density for the underlying index. We use a novel method for the estimation of option-implied volatility surfaces of high quality, which enables the subsequent analysis. We show that the real-world density that we recover is consistent with the observed realized dynamics of the underlying index. This admits the identification of liquidity premia embedded in option price data. We identify and estimate two separate liquidity premia embedded in S&P 500 index options that are consistent with previous findings in the literature.     

The performance of model based option trading strategies

This paper analyzes returns to trading strategies in options markets that exploit information given by a theoretical asset pricing model. We examine trading strategies in which a positive portfolio weight is assigned to assets which market prices exceed the price of a theoretical asset pricing model. We investigate portfolio rules which mimic standard mean-variance analysis is used to construct optimal model based portfolio weights. In essence, these portfolio rules allow estimation risk, as well as price risk to be approximately hedged. An empirical exercise shows that the portfolio rules give out-of-sample Sharpe ratios exceeding unity for S&P 500 options. Portfolio returns have no discernible correlation with systematic risk factors, which is troubling for traditional risk based asset pricing explanations.     

Generalized parameter functions for option pricing

We extend the benchmark nonlinear deterministic volatility regression functions of Dumas et al. (1998) to provide a semi-parametric method where an enhancement of the implied parameter values is used in the parametric option pricing models. Besides volatility, skewness and kurtosis of the asset return distribution can also be enhanced. Empirical results, using closing prices of the S&P 500 index call options (in one day ahead out-of-sample pricing tests), strongly support our method that compares favorably with a model that admits stochastic volatility and random jumps. Moreover, it is found to be superior in various robustness tests. Our semi-parametric approach is an effective remedy to the curse of dimensionality presented in nonparametric estimation and its main advantage is that it delivers theoretically consistent option prices and hedging parameters. The economic significance of the approach is tested in terms of hedging, where the evaluation and estimation loss functions are aligned.     

Modelling the implied volatility surface: Does market efficiency matter?: An application to MIB30 index options

We analyze the volatility surface vs. moneyness and time-to-expiration implied by MIBO options written on the MIB30, the most important Italian stock index. We specify and fit a number of models of the implied volatility surface and find that it has a rich and interesting structure that strongly departs from a constant volatility, Black-Scholes benchmark. This result is robust to alternative econometric approaches, including generalized least squares approaches that take into account both the panel structure of the data and the likely presence of heteroskedasticity and serial correlation in the random disturbances. Finally we show that the degree of pricing efficiency of this options market can strongly condition the results of the econometric analysis and therefore our understanding of the pricing mechanism underlying observed MIBO option prices. Applications to value-at-risk and portfolio choice calculations illustrate the importance of using arbitrage-free data only.     

Two-dimensional risk-neutral valuation relationships for the pricing of options

The Black–Scholes model is based on a one-parameter pricing kernel with constant elasticity. Theoretical and empirical results suggest declining elasticity and, hence, a pricing kernel with at least two parameters. We price European-style options on assets whose probability distributions have two unknown parameters. We assume a pricing kernel which also has two unknown parameters. When certain conditions are met, a two-dimensional risk-neutral valuation relationship exists for the pricing of these options: i.e. the relationship between the price of the option and the prices of the underlying asset and one other option on the asset is the same as it would be under risk neutrality. In this class of models, the price of the underlying asset and that of one other option take the place of the unknown parameters.      

Price discovery in the U.S. stock and stock options markets: A portfolio approach

Option prices vary with not only the underlying asset price, but also volatilities and higher moments. In this paper, we use a portfolio of options to seclude the value change of the portfolio from the impact of volatility and higher moments. We apply this portfolio approach to the price discovery analysis in the U.S. stock and stock options markets. We find that the price discovery on the directional movement of the stock price mainly occurs in the stock market, more so now than before as an increasing proportion of options market makers adopt automated quoting algorithms. Nevertheless, the options market becomes more informative during periods of significant options trading activities. The informativeness of the options quotes increases further when the options trading activity generates net sell or buy pressure on the underlying stock price, even more so when the pressure is consistent with deviations between the stock and the options market quotes. JEL Classification C52, G10, G13, G14     

Multivariate option pricing with time varying volatility and correlations

In this paper we consider option pricing using multivariate models for asset returns. Specifically, we demonstrate the existence of an equivalent martingale measure, we characterize the risk neutral dynamics, and we provide a feasible way for pricing options in this framework. Our application confirms the importance of allowing for dynamic correlation, and it shows that accommodating correlation risk and modeling non-Gaussian features with multivariate mixtures of normals substantially changes the estimated option prices.     

Volatility dynamics for the S&P 500: Further evidence from non-affine,multi-factor jump diffusions

We apply Markov chain Monte Carlo methods to time series data on S&P 500 index returns, and to its option prices via a term structure of VIX indices, to estimate 18 different affine and non-affine stochastic volatility models with one or two variance factors, and where jumps are allowed in both the price and the instantaneous volatility. The in-sample fit to the VIX term structure shows that the second (stochastic long-term volatility) factor is required to fit the VIX term structure. Out-of-sample tests on the fit to individual option prices, as well as in-sample tests, show that the inclusion of jumps is less important than allowing for non-affine dynamics. The estimation and testing periods together cover more than 21 years of daily data.     

Pure jump Lévy processes for asset price modelling

The goal of the paper is to show that some types of Lévy processes such as the hyperbolic motion and the CGMY are particularly suitable for asset price modelling and option pricing. We wish to review some fundamental mathematic properties of Lévy distributions, such as the one of infinite divisibility, and how they translate observed features of asset price returns. We explain how these processes are related to Brownian motion, the central process in finance, through stochastic time changes which can in turn be interpreted as a measure of the economic activity. Lastly, we focus on two particular classes of pure jump Lévy processes, the generalized hyperbolic model and the CGMY models, and report on the goodness of fit obtained both on stock prices and option prices.     

The information content of implied volatilities and model-free volatility expectations: Evidence from options written on individual stocks

We measure the volatility information content of stock options for individual firms using option prices for 149 US firms and the S&P 100 index. We use ARCH and regression models to compare volatility forecasts defined by historical stock returns, at-the-money implied volatilities and model-free volatility expectations for every firm. For 1-day-ahead estimation, a historical ARCH model outperforms both of the volatility estimates extracted from option prices for 36% of the firms, but the option forecasts are nearly always more informative for those firms that have the more actively traded options. When the prediction horizon extends until the expiry date of the options, the option forecasts are more informative than the historical volatility for 85% of the firms. However, at-the-money implied volatilities generally outperform the model-free volatility expectations.     

Recovering Probability Distributions from Option Prices

This article derives underlying asset risk-neutral probability distributions of European options on the S&P 500 index. Nonparametric methods are used to choose probabilities that minimize an objective function subject to requiring that the probabilities are consistent with observed option and underlying asset prices. Alternative optimization specifications produce approximately the same implied distributions. A new and fast optimization technique for estimating probability distributions based on maximizing the smoothness of the resulting distribution is proposed. Since the crash, the risk-neutral probability of a three (four) standard deviation decline in the index (about ?36 percent (?46 percent) over a year) is about 10 (100) times more likely than under the assumption of lognormality.     

Is there price discovery in equity options?

We use tick-by-tick quote data for 39 liquid US stocks and options on them, and we focus on events when the two markets disagree about the stock price in the sense that the option-implied stock price obtained from the put-call parity relation is inconsistent with the actual stock price. Option market quotes adjust to eliminate the disagreement, while the stock market quotes behave normally, as if there were no disagreement. The disagreement events are typically precipitated by stock price movements and display signed option volume in the direction that tends to eliminate the disagreements. These results show that option price quotes do not contain economically significant information about future stock prices beyond what is already reflected in current stock prices, i.e., no economically significant price discovery occurs in the option market. We also find no option market price discovery using a much larger sample of disagreement events based on a weaker definition of a disagreement, which verifies that the findings for the primary sample are not due to unusual or unrepresentative market behavior during the put-call parity violations.     

Predictable dynamics in implied volatility surfaces from OTC currency options

Recent empirical studies report predictable dynamics in the volatility surfaces that are implied by observed index option prices, such as those prescribed by general equilibrium models. Using an extensive data set from the over-the-counter options market, we document similar predictability in the factors that capture the daily variation of surfaces implied by options on 25 different foreign exchange rates. We proceed to demonstrate that simple vector autoregressive specifications for the factors can help produce accurate out-of-sample forecasts of the systematic component of the surface at short horizons. Profitable delta-hedged positions can be set up based on these forecasts; however, profits disappear when typical transaction costs are taken into account and when trading rules on wide segments of the surface are sought.