A QUANTIZATION TREE METHOD FOR PRICING AND HEDGING MULTIDIMENSIONAL AMERICAN OPTIONS

We present here the quantization method which is well-adapted for the pricing and hedging of American options on a basket of assets. Its purpose is to compute a large number of conditional expectations by projection of the diffusion on optimal grids designed to minimize the (square mean) projection error ( Graf and Luschgy 2000 ). An algorithm to compute such grids is described. We provide results concerning the orders of the approximation with respect to the regularity of the payoff function and the global size of the grids. Numerical tests are performed in dimensions 2, 4, 5, 6, 10 with American style exchange options. They show that theoretical orders are probably pessimistic.     

PRICING AND HEDGING DOUBLE-BARRIER OPTIONS: A PROBABILISTIC APPROACH

Barrier options have become increasingly popular over the last few years. Less expensive than standard options, they may provide the appropriate hedge in a number of risk management strategies. In the case of a single-barrier option, the valuation problem is not very difficult (see Merton 1973 and Goldman, Sosin, and Gatto 1979). the situation where the option gets knocked out when the underlying instrument hits either of two well-defined boundaries is less straightforward. Kunitomo and Ikeda (1992) provide a pricing formula expressed as the sum of an infinite series whose convergence is studied through numerical procedures and suggested to be rapid. We follow a methodology which proved quite successful in the case of Asian options (see Geman and Yor 1992,1993) and which has its roots in some fundamental properties of Brownian motion. This methodology permits the derivation of a simple expression of the Laplace transform of the double-barrir price with respect to its maturity date. the inversion of the Laplace transform using techniques developed by Geman and Eydeland (1995), is then fairly easy to perform.     

A NEW METHOD OF PRICING LOOKBACK OPTIONS

A new method for pricing lookback options (a.k.a. hindsight options) is presented, which simplifies the derivation of analytical formulas for this class of exotics in the Black-Scholes framework. Underlying the method is the observation that a lookback option can be considered as an integrated form of a related barrier option. The integrations with respect to the barrier price are evaluated at the expiry date to derive the payoff of an equivalent portfolio of European-type binary options. The arbitrage-free price of the lookback option can then be evaluated by static replication as the present value of this portfolio. We illustrate the method by deriving expressions for generic, standard floating-, fixed-, and reverse-strike lookbacks, and then show how the method can be used to price the more complex partial-price and partial-time lookback options. The method is in principle applicable to frameworks with alternative asset-price dynamics to the Black-Scholes world.     

ANAYTICAL SOLUTIONS FOR THE PRICING OF AMERICAN BOND AND YIELD OPTIONS1

In this paper we use the Cox, Ingersoll, and Ross (1985b) single-factor, term structure model and extend it to the pricing of American default-free bond puts. We provide a quasi-analytical formula for these option prices based on recently established mathematical results for Bessel bridges, coupled with the optimal stopping time method. We extend our results to another interest rate contingent claim and provide a quasi-analytical solution for American yield option prices which illustrates the flexibility of our framework.     

A METHOD FOR PRICING AMERICAN OPTIONS USING SEMI‐INFINITE LINEAR PROGRAMMING

We introduce a new approach for the numerical pricing of American options. The main idea is to choose a finite number of suitable excessive functions (randomly) and to find the smallest majorant of the gain function in the span of these functions. The resulting problem is a linear semi‐infinite programming problem, that can be solved using standard algorithms. This leads to good upper bounds for the original problem. For our algorithms no discretization of space and time and no simulation is necessary. Furthermore it is applicable even for high‐dimensional problems. The algorithm provides an approximation of the value not only for one starting point, but for the complete value function on the continuation set, so that the optimal exercise region and, for example, the Greeks can be calculated. We apply the algorithm to (one‐ and) multidimensional diffusions and show it to be fast and accurate.     

A STATE-SPACE PARTITIONING METHOD FOR PRICING HIGH-DIMENSIONAL AMERICAN-STYLE OPTIONS

The pricing of American-style options by simulation-based methods is an important but difficult task primarily due to the feature of early exercise, particularly for high-dimensional derivatives. In this paper, a bundling method based on quasi-Monte Carlo sequences is proposed to price high-dimensional American-style options. The proposed method substantially extends Tilley's bundling algorithm to higher-dimensional situations. By using low-discrepancy points, this approach partitions the state space and forms bundles. A dynamic programming algorithm is then applied to the bundles to estimate the continuation value of an American-style option. A convergence proof of the algorithm is provided. A variety of examples with up to 15 dimensions are investigated numerically and the algorithm is able to produce computationally efficient results with good accuracy.     

A UNIFIED FRAMEWORK FOR PRICING CREDIT AND EQUITY DERIVATIVES

We propose a model which can be jointly calibrated to the corporate bond term structure and equity option volatility surface of the same company. Our purpose is to obtain explicit bond and equity option pricing formulas that can be calibrated to find a risk neutral model that matches a set of observed market prices. This risk neutral model can then be used to price more exotic, illiquid, or over‐the‐counter derivatives. We observe that our model matches the equity option implied volatility surface well since we properly account for the default risk in the implied volatility surface. We demonstrate the importance of accounting for the default risk and stochastic interest rate in equity option pricing by comparing our results to Fouque et al., which only accounts for stochastic volatility.     

HEDGING STRATEGIES AND MINIMAL VARIANCE PORTFOLIOS FOR EUROPEAN AND EXOTIC OPTIONS IN A LÉVY MARKET

This paper presents hedging strategies for European and exotic options in a Lévy market. By applying Taylor’s theorem, dynamic hedging portfolios are constructed under different market assumptions, such as the existence of power jump assets or moment swaps. In the case of European options or baskets of European options, static hedging is implemented. It is shown that perfect hedging can be achieved. Delta and gamma hedging strategies are extended to higher moment hedging by investing in other traded derivatives depending on the same underlying asset. This development is of practical importance as such other derivatives might be readily available. Moment swaps or power jump assets are not typically liquidly traded. It is shown how minimal variance portfolios can be used to hedge the higher order terms in a Taylor expansion of the pricing function, investing only in a risk‐free bank account, the underlying asset, and potentially variance swaps. The numerical algorithms and performance of the hedging strategies are presented, showing the practical utility of the derived results.     

A TEST OF A GENERAL EQUILIBRIUM STOCK OPTION PRICING MODEL

An empirical version of the Cox, Ingersoll, and Ross (1985a) call option pricing model is derived, assuming execution price uncertainty in the options market. the pricing restrictions come in the form of moment conditions in the option pricing error. These can be estimated and tested using a version of the method of simulated moments (MSM). Simulation estimates, obtained by discretely approximating the risk-neutral processes of the underlying stock price and the interest rate, are substituted for analytically unknown call prices. the asymptotics and other aspects of the MSM estimator are discussed. the model is tested on transaction prices at 15-minute intervals. It substantially outperforms the Black-Scholes model. the empirical success of the Cox-Ingersoll-Ross model implies that the continuous-time interest rate implicit in synchronous transaction quotes of 90-day Treasury-bill futures contracts is an-albeit noisy-proxy for the instantaneous volatility on common stock. the process of the instantaneous volatility is found to be close to nonstationary. It is well approximated by a heteroskedastic unit-root process. With this approximation, the Cox-Ingersoll-Ross model only slightly overprices long-maturity options.     

ON THE AMERICAN OPTION PROBLEM

We show how the change-of-variable formula with local time on curves derived recently in Peskir (2002) can be used to prove that the optimal stopping boundary for the American put option can be characterized as the unique solution of a nonlinear integral equation arising from the early exercise premium representation. This settles the question raised in Myneni (1992) and dating back to McKean (1965) .     

NONCONVEXITY OF THE OPTIMAL EXERCISE BOUNDARY FOR AN AMERICAN PUT OPTION ON A DIVIDEND‐PAYING ASSET

We prove that when the dividend rate of the underlying asset following a geometric Brownian motion is slightly larger than the risk‐free interest rate, the optimal exercise boundary of the American put option is not convex.     

PRICING DISCRETELY MONITORED BARRIER OPTIONS AND DEFAULTABLE BONDS IN LÉVY PROCESS MODELS: A FAST HILBERT TRANSFORM APPROACH

This paper presents a novel method to price discretely monitored single- and double-barrier options in Lévy process-based models. The method involves a sequential evaluation of Hilbert transforms of the product of the Fourier transform of the value function at the previous barrier monitoring date and the characteristic function of the (Esscher transformed) Lévy process. A discrete approximation with exponentially decaying errors is developed based on the Whittaker cardinal series (Sinc expansion) in Hardy spaces of functions analytic in a strip. An efficient computational algorithm is developed based on the fast Hilbert transform that, in turn, relies on the FFT-based Toeplitz matrix–vector multiplication. Our method also provides a natural framework for credit risk applications, where the firm value follows an exponential Lévy process and default occurs at the first time the firm value is below the default barrier on one of a discrete set of monitoring dates.