A unified approach for the pricing of options relating to averages

In this paper, we consider generalized Asian options and propose a unified approximation method for the pricing of such options when the underlying process is a diffusion. Through numerical examples, we show that our approximation method is accurate enough to be used in practice for the pricing of any type of Asian options that has been treated separately in the literature. Comparisons are made with the existing methods in the literature to support the usefulness of our method.     

Analytical approximations for the critical stock prices of American options: a performance comparison

Many efficient and accurate analytical methods for pricing American options now exist. However, while they can produce accurate option prices, they often do not give accurate critical stock prices. In this paper, we propose two new analytical approximations for American options based on the quadratic approximation. We compare our methods with existing analytical methods including the quadratic approximations in Barone-Adesi and Whaley (J Finance 42:301–320, 1987) and Barone-Adesi and Elliott (Stoch Anal Appl 9(2):115–131, 1991), the lower bound approximation in Broadie and Detemple (Rev Financial Stud 9:1211–1250, 1996), the tangent approximation in Bunch and Johnson (J Finance 55(5):2333–2356, 2000), the Laplace inversion method in Zhu (Int J Theor Appl Finance 9(7):1141–1177, 2006b), and the interpolation method in Li (Working paper, 2008). Both of our methods give much more accurate critical stock prices than all the existing methods above.     

Pricing bounds and approximations for discrete arithmetic Asian options under time-changed Lévy processes

We derive efficient and accurate analytical pricing bounds and approximations for discrete arithmetic Asian options under time-changed Lévy processes. By extending the conditioning variable approach, we derive the lower bound on the Asian option price and construct an upper bound based on the sharp lower bound. We also consider the general partially exact and bounded (PEB) approximations, which include the sharp lower bound and partially conditional moment matching approximation as special cases. The PEB approximations are known to lie between a sharp lower bound and an upper bound. Our numerical tests show that the PEB approximations to discrete arithmetic Asian option prices can produce highly accurate approximations when compared to other approximation methods. Our proposed approximation methods can be readily applied to pricing Asian options under most common types of underlying asset price processes, like the Heston stochastic volatility model nested in the class of time-changed Lévy processes with the leverage effect.     

A quasi-analytical interpolation method for pricing American options under general multi-dimensional diffusion processes

We present a quasi-analytical method for pricing multi-dimensional American options based on interpolating two arbitrage bounds, along the lines of Johnson in J Financ Quant Anal 18(1):141–148 (1983). Our method allows for the close examination of the interpolation parameter on a rigorous theoretical footing instead of empirical regression. The method can be adapted to general diffusion processes as long as quick and accurate pricing methods exist for the corresponding European and perpetual American options. The American option price is shown to be approximately equal to an interpolation of two European option prices with the interpolation weight proportional to a perpetual American option. In the Black-Scholes model, our method achieves the same efficiency as the quadratic approximation of Barone-Adesi and Whaley in J Financ 42:301–320 (1987), with our method being generally more accurate for out-of-the-money and long-maturity options. When applied to Heston’s stochastic volatility model, our method is shown to be extremely efficient and fairly accurate.     

Pricing options on discrete realized variance with partially exact and bounded approximations

We derive efficient and accurate analytic approximation formulas for pricing options on discrete realized variance (DRV) under affine stochastic volatility models with jumps using the partially exact and bounded (PEB) approximations. The PEB method is an enhanced extension of the conditioning variable approach commonly used in deriving analytic approximation formulas for pricing discrete Asian style options. By adopting either the conditional normal or gamma distribution approximation based on some asymptotic behaviour of the DRV of the underlying asset price process, we manage to obtain PEB approximation formulas that achieve a high level of numerical accuracy in option values even for short-maturity options on DRV.     

A general closed form option pricing formula

A new method to retrieve the risk-neutral probability measure from observed option prices is developed and a closed form pricing formula for European options is obtained by employing a modified Gram–Charlier series expansion, known as the Gauss–Hermite expansion. This expansion converges for fat-tailed distributions commonly encountered in the study of financial returns. The expansion coefficients can be calibrated from observed option prices and can also be computed, for example, in models with the probability density function or the characteristic function known in closed form. We investigate the properties of the new option pricing model by calibrating it to both real-world and simulated option prices and find that the resulting implied volatility curves provide an accurate approximation for a wide range of strike prices. Based on an extensive empirical study, we conclude that the new approximation method outperforms other methods both in-sample and out-of-sample.

     

The Valuation of Options on Futures Contracts

Rational restrictions are derived for the values of American options on futures contracts. For these options, the optimal policy, in general, involves premature exercise. A model is developed for valuing options on futures contracts in a constant interest rate setting. Despite the fact that premature exercise may be optimal, the value of this American feature appears to be small and a European formula due to Black serves as a useful approximation. Finally, a model is developed to value these options in a world with stochastic interest rates. It is shown that the pricing errors caused by ignoring the location of the interest rate (relative to its long-run mean) range from ?5% to 7%, when the current rate is ±200 basis points from its long-run value. The role of interest rate expectations is, therefore, crucial to the valuation. Optimal exercise policies are found from numerical methods for both models.     

Determining the implied volatility for American options using the QAM

This paper describes an efficient numerical procedure which may be used to determine implied volatilities for American options using the quadratic approximation method. Simulation results are presented. The procedure usually converges in five or six iterations with extreme accuracy under a wide variety of option market conditions. A comparison of American implied volatilities with European model implied volatilities indicates that significant differences may arise. This suggests that reliance on European model volatilities estimates may lead to significant pricing errors.     

Efficient,exact algorithms for asian options with multiresolution lattices

Asian options are a kind of path-dependent derivative. How to price such derivatives efficiently and accurately has been a long-standing research and practical problem. This paper proposes a novel multiresolution (MR) trinomial lattice for pricing European- and American-style arithmetic Asian options. Extensive experimental work suggests that this new approach is both efficient and more accurate than existing methods. It also computes the numerical delta accurately. The MR algorithm is exact as no errors are introduced during backward induction. In fact, it may be the first exact discrete-time algorithm to break the exponential-time barrier. The MR algorithm is guaranteed to converge to the continuous-time value. This revised version was published online in June 2006 with corrections to the Cover Date.     

Calibration to American options: numerical investigation of the de-Americanization method

American options are the reference instruments for the model calibration of a large and important class of single stocks. For this task, a fast and accurate pricing algorithm is indispensable. The literature mainly discusses pricing methods for American options that are based on Monte Carlo, tree and partial differential equation methods. We present an alternative approach that has become popular under the name de-Americanization in the financial industry. The method is easy to implement and enjoys fast run-times (compared to a direct calibration to American options). Since it is based on ad hoc simplifications, however, theoretical results guaranteeing reliability are not available. To quantify the resulting methodological risk, we empirically test the performance of the de-Americanization method for calibration. We classify the scenarios in which de-Americanization performs very well. However, we also identify the cases where de-Americanization oversimplifies and can result in large errors.     

Option pricing with quadratic volatility: a revisit

This paper considers the pricing of European options on assets that follow a stochastic differential equation with a quadratic volatility term. We correct several errors in the existing literature, extend the pricing formulas to arbitrary root configurations, and list alternative representations of option pricing formulas to improve computational performance. Our exposition is based entirely on probabilistic arguments, adding a fresh perspective and new intuition to the existing PDE-dominated literature on the subject. Our main tools are martingale methods and shifts of probability measures; the fact that the underlying process is typically a strict local martingale is carefully considered throughout the paper.     

Pricing Discrete Barrier Options Under Stochastic Volatility

This paper proposes a new approximation method for pricing barrier options with discrete monitoring under stochastic volatility environment. In particular, the integration-by-parts formula and the duality formula in Malliavin calculus are effectively applied in pricing barrier options with discrete monitoring. To the best of our knowledge, this paper is the first one that shows an analytical approximation for pricing discrete barrier options with stochastic volatility models. Furthermore, it provides numerical examples for pricing double barrier call options with discrete monitoring under Heston and λ-SABR models.     

Analytical pricing of American options

By using the homotopy analysis method, we derive a new explicit approximate formula for the optimal exercise boundary of American options on an underlying asset with dividend yields. Compared with highly accurate numerical values, the new formula is shown to be valid for up to 2?years of time to maturity, which is ten times longer than existing explicit approximate formulas. The option price errors computed with our formula are within a few cents for American options that have moneyness (the ratio between stock and strike prices) from 0.8 to 1.2, strike prices of 100 dollars and 2?years to maturity.     

Asymptotic and exact pricing of options on variance

We consider the pricing of derivatives written on the discretely sampled realized variance of an underlying security. In the literature, the realized variance is usually approximated by its continuous-time limit, the quadratic variation of the underlying log-price. Here, we characterize the small-time limits of options on both objects. We find that the difference between them strongly depends on whether or not the stock price process has jumps. Subsequently, we propose two new methods to evaluate the prices of options on the discretely sampled realized variance. One of the methods is approximative; it is based on correcting prices of options on quadratic variation by our asymptotic results. The other method is exact; it uses a novel randomization approach and applies Fourier?CLaplace techniques. We compare the methods and illustrate our results by some numerical examples.     

THE PRICING AND HEDGING OF LIMITED EXERCISE CAPS AND SPREADS

Capped options are barrier option spreads that automatically create simultaneous long and short positions. Exchange-traded capped options were introduced in 1991, though with limited volume. Such options, however, have traded on the over-the-counter markets for several years. Most of these options have the unusual feature that they automatically exercise when the underlying asset closes beyond a critical strike, making them a hybrid of European and American options. In this paper I present their boundary conditions and examine the prices, deltas, gammas, and thetas of caps as well as spreads constructed with European and American options. I also examine the effect of permitting exercise based only on the closing price as opposed to exercise at any time the critical strike is reached. I show that assuming that exercise can occur at any time can lead to serious pricing errors. The results have implications for the pricing of barrier options in general, which nearly always exercise early based only on the closing price.     

Convergence of numerical methods for valuing path-dependent options using interpolation

One method for valuing path-dependent options is the augmented state space approach described in Hull and White (1993) and Barraquand and Pudet (1996), among others. In certain cases, interpolation is required because the number of possible values of the additional state variable grows exponentially. We provide a detailed analysis of the convergence of these algorithms. We show that it is possible for the algorithm to be non-convergent, or to converge to an incorrect answer, if the interpolation scheme is selected in appropriately. We concentrate on Asian options, due to their popularity and because of some errors in the previous literature. This revised version was published online in June 2006 with corrections to the Cover Date.     

Assessing the performance of symmetric and asymmetric implied volatility functions

This study examines several alternative symmetric and asymmetric model specifications of regression-based deterministic volatility models to identify the one that best characterizes the implied volatility functions of S&P 500 Index options in the period 1996–2009. We find that estimating the models with nonlinear least squares, instead of ordinary least squares, always results in lower pricing errors in both in- and out-of-sample comparisons. In-sample, asymmetric models of the moneyness ratio estimated separately on calls and puts provide the overall best performance. However, separating calls from puts violates the put-call-parity and leads to severe model mis-specification problems. Out-of-sample, symmetric models that use the logarithmic transformation of the strike price are the overall best ones. The lowest out-of-sample pricing errors are observed when implied volatility models are estimated consistently to the put-call-parity using the joint data set of out-of-the-money options. The out-of-sample pricing performance of the overall best model is shown to be resilient to extreme market conditions and compares quite favorably with continuous-time option pricing models that admit stochastic volatility and random jump risk factors.     

High-order accurate implicit finite difference method for evaluating American options

A numerical method is presented for valuing vanilla American options on a single asset that is up to fourth-order accurate in the log of the asset price, and second-order accurate in time. The method overcomes the standard difficulty encountered in developing high-order accurate finite difference schemes for valuing American options; that is, the lack of smoothness in the option price at the critical boundary. This is done by making special corrections to the right-hand side of the differnce equations near the boundary, so they retain their level of accuracy. These corrections are easily evaluated using estimates of the boundary location and jump in the gamma that occurs there, such as those developed by Carr and Eaguet. Results of numerical experiments are presented comparing the method with more standard finite difference methods.     

Interest Rate Futures Options and Interest Rate Options

This paper derives pricing models of interest rate options and interest rate futures options. The models utilize the arbitrage-free interest rate movements model of Ho and Lee. In their model, they take the initial term structure as given, and for the subsequent periods, they only require that the bond prices move relative to each other in an arbitrage-free manner. Viewing the interest rate options as contingent claims to the underlying bonds, we derive the closed-form solutions to the options. Since these models are sufficiently simple, they can be used to investigate empirically the pricing of bond options. We also empirically examine the pricing of Eurodollar futures options. The results show that the model has significant explanatory power and, on average, has smaller estimation errors than Black's model. The results suggest that the model can be used to price options relative to each other, even though they may have different expiration dates and strike prices.     

Efficient Pricing of Ratchet Equity-Indexed Annuities in a Variance-Gamma Economy

Abstract

In this paper we propose a new method for approximating the price of arithmetic Asian options in a Variance-Gamma (VG) economy, which is then applied to the problem of pricing equityindexed annuity contracts. The proposed procedure is an extension to the case of a VG-based model of the moment-matching method developed by Turnbull and Wakeman and Levy for the pricing of this class of path-dependent options in the traditional Black-Scholes setting. The accuracy of the approximation is analyzed against RQMC estimates for the case of ratchet equityindexed annuities with index averaging.