Confined exponential approximations for the valuation of American options

We provide an alternative analytic approximation for the value of an American option using a confined exponential distribution with tight upper bounds. This is an extension of the Geske and Johnson compound option approach and the Ho et al. exponential extrapolation method. Use of a perpetual American put value, and then a European put with high input volatility is suggested in order to provide a tighter upper bound for an American put price than simply the exercise price. Numerical results show that the new method not only overcomes the deficiencies in existing two-point extrapolation methods for long-term options but also further improves pricing accuracy for short-term options, which may substitute adequately for numerical solutions. As an extension, an analytic approximation is presented for a two-factor American call option.     

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.     

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.     

Valuation of American partial barrier options

This paper concerns barrier options of American type where the underlying asset price is monitored for barrier hits during a part of the option’s lifetime. Analytic valuation formulas of the American partial barrier options are provided as the finite sum of bivariate normal distribution functions. This approximation method is based on barrier options along with constant early exercise policies. In addition, numerical results are given to show the accuracy of the approximating price. Our explicit formulas provide a very tight lower bound for the option values, and moreover, this method is superior in speed and its simplicity.     

An alternative approach to the valuation of American options and applications

In this paper we examine the structure of American option valuation problems and derive the analytic valuation formulas under general underlying security price processes by an alternative but intuitive method. For alternative diffusion processes, we derive closed-form analytic valuation formulas and analyze the implications of asset price dynamics on the early exercise premiums of American options. In this regard, we introduce useful and interesting diffusion processes into American option-pricing literature, thus providing a wide range of choices of pricing models for various American-type derivative assets. This work offers a useful analytic framework for empirical testing and practical applications such as the valuation of corporate securities and examining the impact of options trading on market micro-structure.     

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.     

A closed-form solution to American options under general diffusion processes

This paper investigates American option pricing under general diffusion processes. Specifically, the underlying asset price is assumed to follow a diffusion process in which both the dividend yield and volatility are functions of time and the underlying asset price. Using the generalized homotopy analysis method, the determination of the early exercise boundary is separated from the valuation procedure of American options. Then, an exact and explicit solution for American options on a dividend-paying stock is derived as a Maclaurin series. In addition, the corresponding optimal early exercise boundary and the Greeks are obtained in closed-form solutions. A nonlinear sequence transformation, the Padé technique, is used to effectively accelerate the convergence of the partial sums of the infinite series. As the homotopy constructed in this paper is based on a generalized deformation with a shape parameter and kernel function, the error of the homotopic approximation could be reduced further for a fixed order. Numerical examples demonstrate the validity, effectiveness, and flexibility of the proposed approach.     

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

Efficient estimation of lower and upper bounds for pricing higher-dimensional American arithmetic average options by approximating their payoff functions

In this paper, we develop an efficient payoff function approximation approach to estimating lower and upper bounds for pricing American arithmetic average options with a large number of underlying assets. The crucial step in the approach is to find a geometric mean which is more tractable than and highly correlated with a given arithmetic mean. Then the optimal exercise strategy for the resultant American geometric average option is used to obtain a low-biased estimator for the corresponding American arithmetic average option. This method is particularly efficient for asset prices modeled by jump-diffusion processes with deterministic volatilities because the geometric mean is always a one-dimensional Markov process regardless of the number of underlying assets and thus is free from the curse of dimensionality. Another appealing feature of our method is that it provides an extremely efficient way to obtain tight upper bounds with no nested simulation involved as opposed to some existing duality approaches. Various numerical examples with up to 50 underlying stocks suggest that our algorithm is able to produce computationally efficient results.     

Pricing average options under time-changed Lévy processes

This paper presents an approximate formula for pricing average options when the underlying asset price is driven by time-changed Lévy processes. Time-changed Lévy processes are attractive to use for a driving factor of underlying prices because the processes provide a flexible framework for generating jumps, capturing stochastic volatility as the random time change, and introducing the leverage effect. There have been very few studies dealing with pricing problems of exotic derivatives on time-changed Lévy processes in contrast to standard European derivatives. Our pricing formula is based on the Gram–Charlier expansion and the key of the formula is to find analytic treatments for computing the moments of the normalized average asset price. In numerical examples, we demonstrate that our formula give accurate values of average call options when adopting Heston’s stochastic volatility model, VG-CIR, and NIG-CIR models.     

Analysis of the Nonlinear Option Pricing Model Under Variable Transaction Costs

In this paper we analyze a nonlinear Black–Scholes model for option pricing under variable transaction costs. The diffusion coefficient of the nonlinear parabolic equation for the price V is assumed to be a function of the underlying asset price and the Gamma of the option. We show that the generalizations of the classical Black–Scholes model can be analyzed by means of transformation of the fully nonlinear parabolic equation into a quasilinear parabolic equation for the second derivative of the option price. We show existence of a classical smooth solution and prove useful bounds on the option prices. Furthermore, we construct an effective numerical scheme for approximation of the solution. The solutions are obtained by means of the efficient numerical discretization scheme of the Gamma equation. Several computational examples are presented.     

Pricing American options written on two underlying assets

This paper extends the integral transform approach of McKean [Ind. Manage. Rev., 1965, 6, 32–39] and Chiarella and Ziogas [J. Econ. Dyn. Control, 2005, 29, 229–263] to the pricing of American options written on more than one underlying asset under the Black and Scholes [J. Polit. Econ., 1973, 81, 637–659] framework. A bivariate transition density function of the two underlying stochastic processes is derived by solving the associated backward Kolmogorov partial differential equation. Fourier transform techniques are used to transform the partial differential equation to a corresponding ordinary differential equation whose solution can be readily found by using the integrating factor method. An integral expression of the American option written on any two assets is then obtained by applying Duhamel’s principle. A numerical algorithm for calculating American spread call option prices is given as an example, with the corresponding early exercise boundaries approximated by linear functions. Numerical results are presented and comparisons made with other alternative approaches.     

When are Options Overpriced? The Black--Scholes Model and Alternative Characterisations of the Pricing Kernel

An important determinant of option prices is the elasticityof the pricing kernel used to price all claims in the economy.In this paper, we first show that for a given forward priceof the underlying asset, option prices are higher when the elasticityof the pricing kernel is declining than when it is constant.We then investigate the implications of the elasticity of thepricing kernel for the stochastic process followed by the underlyingasset. Given that the underlying information process followsa geometric Brownian motion, we demonstrate that constant elasticityof the pricing kernel is equivalent to a Brownian motion forthe forward price of the underlying asset, so that the Black–Scholesformula correctly prices options on the asset. In contrast,declining elasticity implies that the forward price processis no longer a Brownian motion: it has higher volatility andexhibits autocorrelation. In this case, the Black–Scholesformula underprices all options.     

Algorithmic market making for options

In this article, we tackle the problem of a market maker in charge of a book of options on a single liquid underlying asset. By using an approximation of the portfolio in terms of its vega, we show that the seemingly high-dimensional stochastic optimal control problem of an option market maker is in fact tractable. More precisely, when volatility is modeled using a classical stochastic volatility model—e.g. the Heston model—the problem faced by an option market maker is characterized by a low-dimensional functional equation that can be solved numerically using a Euler scheme along with interpolation techniques, even for large portfolios. In order to illustrate our findings, numerical examples are provided.     

Pricing by American option by approximating its early exercise boundary as a multipiece exponential function

This article proposes to price an American option by approximatingits early exercise boundary as a multipiece exponential function.Closed form formulas are obtained in terms of the bases andexponents of the multipiece exponential function. It is demonstratedthat a three-point extrapolation scheme has the accuracy ofan 800-time-step binomial tree, but is about 130 times faster.An intuitive argument is given to indicate why this seeminglycrude approximation works so well. Our method is very simpleand easy to implement. Comparisons with other leading competingmethods are also included.     

Saddlepoint approximations to option price in a regime-switching model

In this paper we consider the saddlepoint approximation for the valuation of a European-style call option in a Markovian, regime-switching, Black–Scholes–Merton economy, where the price process of an underlying risky asset is assumed to follow a Markov-modulated geometric Brownian motion. The standard option pricing procedure under this model becomes problematic as the occupation time of chains for a given state cannot be evaluated easily. In the case of two-state Markov chains, we present an explicit analytic formula of the cumulant generating function (CGF). When the process has more than two states, an approximate formula of the CGF is provided. We adopt a splitting method to reduce the complexity of computing the matrix exponential function. Then we use these CGFs to develop the use of the saddlepoint approximations. The numerical results show that the saddlepoint approximation is an efficient and reliable approach for option pricing under a multi-state regime-switching model.     

Skew Brownian Motion and Pricing European Options

Abstract

The volatility smile and systematic mispricing of the Black–Scholes option pricing model are the typical motivation for examining stochastic processes other than geometric Brownian motion to describe the underlying stock price. In this paper a new stochastic process is presented, which is a special case of the skew-Brownian motion of Itô and McKean. The process in question is the sum of a standard Brownian motion and an independent reflecting Brownian motion that is similar in construction to the stochastic representation of a skew-normal random variable. This stochastic process is taken in its exponential form to price European options. The derived option price nests the Black–Scholes equation as a special case and is flexible enough to accommodate stochastic volatility as well as stochastic skewness.     

Lapse rate modeling: a rational expectation approach

The surrender option embedded in many life insurance products is a clause that allows policyholders to terminate the contract early. Pricing techniques based on the American Contingent Claim (ACC) theory are often used, though the actual policyholders' behavior is far from optimal. Inspired by many prepayment models for mortgage backed securities, this paper builds a Rational Expectation (RE) model describing the policyholders' behavior in lapsing the contract. A market model with stochastic interest rates is considered, and the pricing is carried out through numerical approximation of the corresponding two-space-dimensional parabolic partial differential equation. Extensive numerical experiments show the differences in terms of pricing and interest rate elasticity between the ACC and RE approaches as well as the sensitivity of the contract price with respect to changes in the policyholders' behavior.