American stochastic volatility call option pricing: A lattice based approach

This study presents a new method of pricing options on assets with stochastic volatility that is lattice based, and can easily accommodate early exercise for American options. Unlike traditional lattice methods, recombination is not a problem in the new model, and it is easily adapted to alternative volatility processes. Approximations are developed for European C.E.V. calls and American stochastic volatility calls. The application of the pricing model to exchange traded calls is also illustrated using a sample of market prices. Modifying the model to price American puts is straightforward, and the approach can easily be extended to other non-recombining lattices.     

American options under stochastic volatility: control variates,maturity randomization & multiscale asymptotics

American options are actively traded worldwide on exchanges, thus making their accurate and efficient pricing an important problem. As most financial markets exhibit randomly varying volatility, in this paper we introduce an approximation of an American option price under stochastic volatility models. We achieve this by using the maturity randomization method known as Canadization. The volatility process is characterized by fast and slow-scale fluctuating factors. In particular, we study the case of an American put with a single underlying asset and use perturbative expansion techniques to approximate its price as well as the optimal exercise boundary up to the first order. We then use the approximate optimal exercise boundary formula to price an American put via Monte Carlo. We also develop efficient control variates for our simulation method using martingales resulting from the approximate price formula. A numerical study is conducted to demonstrate that the proposed method performs better than the least squares regression method popular in the financial industry, in typical settings where values of the scaling parameters are small. Further, it is empirically observed that in the regimes where the scaling parameter value is equal to unity, fast and slow-scale approximations are equally accurate.     

A new approach for option pricing under stochastic volatility

We develop a new approach for pricing European-style contingent claims written on the time T spot price of an underlying asset whose volatility is stochastic. Like most of the stochastic volatility literature, we assume continuous dynamics for the price of the underlying asset. In contrast to most of the stochastic volatility literature, we do not directly model the dynamics of the instantaneous volatility. Instead, taking advantage of the recent rise of the variance swap market, we directly assume continuous dynamics for the time T variance swap rate. The initial value of this variance swap rate can either be directly observed, or inferred from option prices. We make no assumption concerning the real world drift of this process. We assume that the ratio of the volatility of the variance swap rate to the instantaneous volatility of the underlying asset just depends on the variance swap rate and on the variance swap maturity. Since this ratio is assumed to be independent of calendar time, we term this key assumption the stationary volatility ratio hypothesis (SVRH). The instantaneous volatility of the futures follows an unspecified stochastic process, so both the underlying futures price and the variance swap rate have unspecified stochastic volatility. Despite this, we show that the payoff to a path-independent contingent claim can be perfectly replicated by dynamic trading in futures contracts and variance swaps of the same maturity. As a result, the contingent claim is uniquely valued relative to its underlying’s futures price and the assumed observable variance swap rate. In contrast to standard models of stochastic volatility, our approach does not require specifying the market price of volatility risk or observing the initial level of instantaneous volatility. As a consequence of our SVRH, the partial differential equation (PDE) governing the arbitrage-free value of the contingent claim just depends on two state variables rather than the usual three. We then focus on the consistency of our SVRH with the standard assumption that the risk-neutral process for the instantaneous variance is a diffusion whose coefficients are independent of the variance swap maturity. We show that the combination of this maturity independent diffusion hypothesis (MIDH) and our SVRH implies a very special form of the risk-neutral diffusion process for the instantaneous variance. Fortunately, this process is tractable, well-behaved, and enjoys empirical support. Finally, we show that our model can also be used to robustly price and hedge volatility derivatives.     

Stock price distributions with stochastic volatility: an analytic approach

We study the stock price distributions that arise when pricesfollow a diffusion process with a stochastically varying volatilityparameters. We use analytic techniques to derive an explicitclosed-form solution for the case when volatility is drivenby an arithmetic Ornstein-Uhlenbeck (or AR1) process. We thenapply our results to two related problems in the finance literature:(i) options pricing in a world of stochastic volatility, and(ii) the relationship between stochastic volatility and thenature of 'fat tails' in stock price distributions.     

Using Richardson extrapolation techniques to price American options with alternative stochastic processes

In this paper the authors investigate the performance of the original and repeated Richardson extrapolation methods for American option pricing by implementing both the original and modified Geske?CJohnson approximation formulae. A comprehensive numerical comparison includes alternative stochastic processes of the underlying asset price. The numerical results show that whether the original or modified formula is implemented, the Richardson extrapolation techniques work very well. The repeated Richardson extrapolation strongly outperforms the original, especially when the underlying asset price follows a stochastic volatility process. Moreover, this study verifies the feasibility of the estimated error bounds of the American option prices under alternative stochastic processes by applying the repeated Richardson extrapolation method and estimating the interval of true American option values, as well as determining the number of options needed for an approximation to achieve a desired accuracy level.     

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.     

Pricing via recursive quantization in stochastic volatility models

We provide the first recursive quantization-based approach for pricing options in the presence of stochastic volatility. This method can be applied to any model for which an Euler scheme is available for the underlying price process and it allows one to price vanillas, as well as exotics, thanks to the knowledge of the transition probabilities for the discretized stock process. We apply the methodology to some celebrated stochastic volatility models, including the Stein and Stein [Rev. Financ. Stud. 1991, (4), 727–752] model and the SABR model introduced in Hagan et al. [Wilmott Mag., 2002, 84–108]. A numerical exercise shows that the pricing of vanillas turns out to be accurate; in addition, when applied to some exotics like equity-volatility options, the quantization-based method overperforms by far the Monte Carlo simulation.     

Pricing Options under Generalized GARCH and Stochastic Volatility Processes

In this paper, we develop an efficient lattice algorithm to price European and American options under discrete time GARCH processes. We show that this algorithm is easily extended to price options under generalized GARCH processes, with many of the existing stochastic volatility bivariate diffusion models appearing as limiting cases. We establish one unifying algorithm that can price options under almost all existing GARCH specifications as well as under a large family of bivariate diffusions in which volatility follows its own, perhaps correlated, process.     

Smile and default: the role of stochastic volatility and interest rates in counterparty credit risk

In this research, we investigate the impact of stochastic volatility and interest rates on counterparty credit risk (CCR) for FX derivatives. To achieve this we analyse two real-life cases in which the market conditions are different, namely during the 2008 credit crisis where risks are high and a period after the crisis in 2014, where volatility levels are low. The Heston model is extended by adding two Hull–White components which are calibrated to fit the EURUSD volatility surfaces. We then present future exposure profiles and credit value adjustments (CVAs) for plain vanilla cross-currency swaps (CCYS), barrier and American options and compare the different results when Heston-Hull–White-Hull–White or Black–Scholes dynamics are assumed. It is observed that the stochastic volatility has a significant impact on all the derivatives. For CCYS, some of the impact can be reduced by allowing for time-dependent variance. We further confirmed that Barrier options exposure and CVA is highly sensitive to volatility dynamics and that American options’ risk dynamics are significantly affected by the uncertainty in the interest rates.     

Pricing Asian options with stochastic volatility

Abstract

In this paper, we generalize the recently developed dimension reduction technique of Vecer for pricing arithmetic average Asian options. The assumption of constant volatility in Vecer's method will be relaxed to the case that volatility is randomly fluctuating and is driven by a mean-reverting (or ergodic) process. We then use the fast mean-reverting stochastic volatility asymptotic analysis introduced by Fouque, Papanicolaou and Sircar to derive an approximation to the option price which takes into account the skew of the implied volatility surface. This approximation is obtained by solving a pair of one-dimensional partial differential equations.     

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 fast Fourier transform technique for pricing American options under stochastic volatility

This paper develops a non-finite-difference-based method of American option pricing under stochastic volatility by extending the Geske-Johnson compound option scheme. The characteristic function of the underlying state vector is inverted to obtain the vector’s density using a kernel-smoothed fast Fourier transform technique. The method produces option values that are closely in line with the values obtained by finite-difference schemes. It also performs well in an empirical application with traded S&P 100 index options. The method is especially well suited to price a set of options with different strikes on the same underlying asset, which is a task often encountered by practitioners.     

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 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.     

Pricing swaps and options on quadratic variation under stochastic time change models—discrete observations case

We use a forward characteristic function approach to price variance and volatility swaps and options on swaps. The swaps are defined via contingent claims whose payoffs depend on the terminal level of a discretely monitored version of the quadratic variation of some observable reference process. As such a process we consider a class of Levy models with stochastic time change. Our analysis reveals a natural small parameter of the problem which allows a general asymptotic method to be developed in order to obtain a closed-form expression for the fair price of the above products. As examples, we consider the CIR clock change, general affine models of activity rates and the 3/2 power clock change, and give an analytical expression of the swap price. Comparison of the results obtained with a familiar log-contract approach is provided.     

The Valuation of Volatility Options

This paper examines the valuation of European- and American-style volatilityoptions based on a general equilibrium stochastic volatility framework.Properties of the optimal exercise region and of the option price areprovided when volatility follows a general diffusion process. Explicitvaluation formulas are derived in four particular cases. Emphasis is placedon the MRLP (mean-reverting in the log) volatility model which has receivedconsiderable empirical support. In this context we examine the propertiesand hedging behavior of volatility options. Unlike American options,European call options on volatility are found to display concavity at highlevels of volatility.     

American options: the EPV pricing model

Summary We explicitly solve the pricing problem for perpetual American puts and calls, and provide an efficient semi-explicit pricing procedure for options with finite time horizon. Contrary to the standard approach, which uses the price process as a primitive, we model the price process as the expected present value of a stream, which is a monotone function of a Lévy process. Certain processes exhibiting mean-reverting, stochastic volatility and/or switching features can be modeled this way. This specification allows us to consider assets that pay no dividends at all when the level of the underlying stochastic factor is too low, assets that pay dividends at a fixed rate when the underlying stochastic process remains in some range, or capped dividends.The authors are grateful to the anonymous referees for valuable comments and suggestions.     

American option valuation under stochastic interest rates

By applying Ho, Stapleton and Subrahmanyam's (1997, hereafter HSS) generalised Geske–Johnson (1984, hereafter GJ) method, this paper provides analytic solutions for the valuation and hedging of American options in a stochastic interest rate economy. The proposed method simplifies HSS's three-dimensional solution to a one-dimensional solution. The simulations verify that the proposed method is more efficient and accurate than the HSS (1997) method. We illustrate how the price, the delta, and the rho of an American option vary between the stochastic and non-stochastic interest rate models. The magnitude of this effect depends on the moneyness of the option, interest rates, volatilities of the underlying asset price and the bond price, as well as the correlation between them. This revised version was published online in June 2006 with corrections to the Cover Date.     

Stochastic Dominance Bounds on American Option Prices in Markets with Frictions

We derive equilibrium restrictions on the range of the transactionprices of American options on the stock market index and indexfutures. Trading over the lifetime of the options is accountedfor, in contrast to earlier single-period results. The boundson the reservation purchase price of American puts and the reservationwrite price of American calls are tight. We allow the marketto be incomplete and imperfect due to the presence of proportionaltransaction costs in trading the underlying security and dueto bid-ask spreads in option prices. The bounds may be derivedfor any given probability distribution of the return of theunderlying security and admit price jumps and stochastic volatility.We assume that at least some of the traders maximize a time-separable utility function. The bounds are derived by applyingthe weak notion of stochastic dominance and are independentof a trader's particular utility function and initial portfolioposition.     

The evaluation of European compound option prices under stochastic volatility using Fourier transform techniques

Compound options are not only sensitive to future movements of the underlying asset price, but also to future changes in volatility levels. Because the Black–Scholes analytical valuation formula for compound options is not able to incorporate the sensitivity to volatility, the aim of this paper is to develop a numerical pricing procedure for this type of option in stochastic volatility models, specifically focusing on the model of Heston. For this, the compound option value is represented as the difference of its exercise probabilities, which depend on three random variables through a complex functional form. Then the joint distribution of these random variables is uniquely determined by their characteristic function and therefore the probabilities can each be expressed as a multiple inverse Fourier transform. Solving the inverse Fourier transform with respect to volatility, we can reduce the pricing problem from three to two dimensions. This reduced dimensionality simplifies the application of the fast Fourier transform (FFT) method developed by Dempster and Hong when transferred to our stochastic volatility framework. After combining their approach with a new extension of the fractional FFT technique for option pricing to the two-dimensional case, it is possible to obtain good approximations to the exercise probabilities. The resulting upper and lower bounds are then compared with other numerical methods such as Monte Carlo simulations and show promising results.