Pricing exotic options using MSL-MC

Today, better numerical approximations are required for multi-dimensional SDEs to improve on the poor performance of the standard Monte Carlo pricing method. With this aim in mind, this paper presents a method (MSL-MC) to price exotic options using multi-dimensional SDEs (e.g. stochastic volatility models). Usually, it is the weak convergence property of numerical discretizations that is most important, because, in financial applications, one is mostly concerned with the accurate estimation of expected payoffs. However, in the recently developed Multilevel Monte Carlo path simulation method (ML-MC), the strong convergence property plays a crucial role. We present a modification to the ML-MC algorithm that can be used to achieve better savings. To illustrate these, various examples of exotic options are given using a wide variety of payoffs, stochastic volatility models and the new Multischeme Multilevel Monte Carlo method (MSL-MC). For standard payoffs, both European and Digital options are presented. Examples are also given for complex payoffs, such as combinations of European options (Butterfly Spread, Strip and Strap options). Finally, for path-dependent payoffs, both Asian and Variance Swap options are demonstrated. This research shows how the use of stochastic volatility models and the θ scheme can improve the convergence of the MSL-MC so that the computational cost to achieve an accuracy of O(ε) is reduced from O(ε^?3) to O(ε^?2) for a payoff under global and non-global Lipschitz conditions.     

An efficient algorithm for pricing barrier options in arbitrage-free binomial models with calibrated drift terms

The interrelation between the drift coefficient of price processes on arbitrage-free financial markets and the corresponding transition probabilities induced by a martingale measure is analysed in a discrete setup. As a result, we obtain a flexible setting that encompasses most arbitrage-free binomial models. It is argued that knowledge of the link between drift and transition probabilities may be useful for pricing derivatives such as barrier options. The idea is illustrated in a simple example and later extended to a general numerical procedure. The results indicate that the option values in our fitted drift model converge much faster to closed-form solutions of continuous models for a wider range of contract specifications than those of conventional binomial models.     

Pseudospectral methods for pricing options

Models with two or more risk sources have been widely applied in option pricing in order to capture volatility smiles and skews. However, the computational cost of implementing these models can be large—especially for American-style options. This paper illustrates how numerical techniques called ‘pseudospectral’ methods can be used to solve the partial differential and partial integro-differential equations that apply to these multifactor models. The method offers significant advantages over finite-difference and Monte Carlo simulation schemes in terms of accuracy and computational cost.     

Identifying,valuing and hedging of embedded options in non-maturity deposits

Non-maturity deposits like savings accounts or demand deposits contain significant option risks caused by the bank’s discretionary pricing and the customers’ withdrawal right. Option risks follow from inherent non-linear factor exposures. I propose an ordinal response model for deposit rate jumps to identify non-linear factor exposures and a discrete-time term structure model to value the resulting option risks and to derive hedge measures “outside the model”. My delta profile resembles a constant maturity swap, but vega and gamma are more pronounced, which demonstrates that the widespread practice of static hedging with zero bonds is inadequate.     

Some calculations for Israeli options

Recently Kifer (2000) introduced the concept of an Israeli (or Game) option. That is a general American-type option with the added possibility that the writer may terminate the contract early inducing a payment exceeding the holders claim had they exercised at that moment. Kifer shows that pricing and hedging of these options reduces to evaluating a saddle point problem associated with Dynkin games. In this short text we give two examples of perpetual Israeli options where the solutions are explicit.Received: December 2002, Mathematics Subject Classification: 90A09, 60J40, 90D15JEL Classification: G13, C73I would like to express thanks to Chris Rogers for a valuable conversation.     

Generative Bayesian neural network model for risk-neutral pricing of American index options

Financial models with stochastic volatility or jumps play a critical role as alternative option pricing models for the classical Black–Scholes model, which have the ability to fit different market volatility structures. Recently, machine learning models have elicited considerable attention from researchers because of their improved prediction accuracy in pricing financial derivatives. We propose a generative Bayesian learning model that incorporates a prior reflecting a risk-neutral pricing structure to provide fair prices for the deep ITM and the deep OTM options that are rarely traded. We conduct a comprehensive empirical study to compare classical financial option models with machine learning models in terms of model estimation and prediction using S&P 100 American put options from 2003 to 2012. Results indicate that machine learning models demonstrate better prediction performance than the classical financial option models. Especially, we observe that the generative Bayesian neural network model demonstrates the best overall prediction performance.     

A one-factor volatility smile model with closed-form solutions for European options

The common practice of using different volatilities for options of different strikes in the Black-Scholes (1973) model imposes inconsistent assumptions on underlying securities. The phenomenon is referred to as the volatility smile. This paper addresses this problem by replacing the Brownian motion or, alternatively, the Geometric Brownian motion in the Black-Scholes model with a two-piece quadratic or linear function of the Brownian motion. By selecting appropriate parameters of this function we obtain a wide range of shapes of implied volatility curves with respect to option strikes. The model has closed-form solutions for European options, which enables fast calibration of the model to market option prices. The model can also be efficiently implemented in discrete time for pricing complex options.
G1     

Primal–dual quasi-Monte Carlo simulation with dimension reduction for pricing American options

The pricing of American options is one of the most challenging problems in financial engineering due to the involved optimal stopping time problem, which can be solved by using dynamic programming (DP). But applying DP is not always practical, especially when the state space is high dimensional. However, the curse of dimensionality can be overcome by Monte Carlo (MC) simulation. We can get lower and upper bounds by MC to ensure that the true price falls into a valid confidence interval. During the recent decades, progress has been made in using MC simulation to obtain both the lower bound by least-squares Monte Carlo method (LSM) and the upper bound by duality approach. However, there are few works on pricing American options using quasi-Monte Carlo (QMC) methods, especially to compute the upper bound. For comparing the sample variances and standard errors in the numerical experiments, randomized QMC (RQMC) methods are usually used. In this paper, we propose to use RQMC to replace MC simulation to compute both the lower bound (by the LSM) and the upper bound (by the duality approach). Moreover, we propose to use dimension reduction techniques, such as the Brownian bridge, principal component analysis, linear transformation and the gradients based principle component analysis. We perform numerical experiments on American–Asian options and American max-call options under the Black–Scholes model and the variance gamma model, in which the options have the path-dependent feature or are written on multiple underlying assets. We find that RQMC in combination with dimension reduction techniques can significantly increase the efficiency in computing both the lower and upper bounds, resulting in better estimates and tighter confidence intervals of the true price than pure MC simulation.     

Efficient and accurate quadratic approximation methods for pricing Asian strike options

This study is on valuing Asian strike options and presents efficient and accurate quadratic approximation methods that work extremely well, both with regard to the volatility of a wide range of underlying assets, and longer average time windows. We demonstrate that most of the well-known quadratic approximation methods used in the literature for pricing Asian strike options are special cases of our model, with the numerical results demonstrating that our method significantly outperforms the other quadratic approximation methods examined here. Using our method for the calculation of hundreds of Asian strike options, the pricing errors (in terms of the root mean square errors) are reasonably small. Compared with the Monte Carlo benchmark method, our method is shown to be rapid and accurate. We further extend our method to the valuing of quanto forward-starting Asian strike options, with the pricing accuracy of these options being largely the same as the pricing of plain vanilla Asian strike options.     

Pricing of vanilla and first-generation exotic options in the local stochastic volatility framework: survey and new results

Stochastic volatility (SV) and local stochastic volatility (LSV) processes can be used to model the evolution of various financial variables such as FX rates, stock prices and so on. Considerable efforts have been devoted to pricing derivatives written on underliers governed by such processes. Many issues remain, though, including the efficacy of the standard alternating direction implicit (ADI) numerical methods for solving SV and LSV pricing problems. In general, the amount of required computations for these methods is very substantial. In this paper, we address some of these issues and propose a viable alternative to the standard ADI methods based on Galerkin-Ritz ideas. We also discuss various approaches to solving the corresponding pricing problems in a semi-analytical fashion. We use the fact that in the zero correlation case some of the pricing problems can be solved analytically, and develop a closed-form series expansion in powers of correlation. We perform a thorough benchmarking of various numerical solutions by using analytical and semi-analytical solutions derived in the paper.     

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.      

Numerical solutions to an integro-differential parabolic problem arising in the pricing of financial options in a Levy market

We study the numerical solutions for an integro-differential parabolic problem modeling a process with jumps and stochastic volatility in financial mathematics. We present two general algorithms to calculate numerical solutions. The algorithms are implemented in PDE2D, a general-purpose, partial differential equation solver.     

A semi-analytic method for valuing high-dimensional options on the maximum and minimum of multiple assets

Valuing high-dimensional options has many important applications in finance but when the true distributions are unknown or complex, numerical approximations must be used. Approximation methods based on Monte-Carlo simulation show a steep trade-off between estimation accuracy and computational efficiency. This article presents an alternative semi-analytic approximation method for pricing options on the maximum or minimum of multiple assets with unknown distributions. Computational efficiency is shown to improve significantly without sacrificing estimation accuracy. The method is illustrated with applications to options on underlying assets with mean-reverting prices, time-dependent correlations, and stochastic volatility The authors would like to thank the two anonymous referees, the associate editor, and Dr. Jess H. Chua at the University of Calgary for valuable comments and insights on this research. This research was partly supported by NUS grant R-146-000-059-112     

On the Robustness of Least-Squares Monte Carlo (LSM) for Pricing American Derivatives

This paper analyses the robustness of Least-Squares Monte Carlo, a technique proposed by Longstaff and Schwartz (2001) for pricing American options. This method is based on least-squares regressions in which the explanatory variables are certain polynomial functions. We analyze the impact of different basis functions on option prices. Numerical results for American put options show that this approach is quite robust to the choice of basis functions. For more complex derivatives, this choice can slightly affect option prices. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.     

Optimal exercise of executive stock options and implications for firm cost

This paper conducts a comprehensive study of the optimal exercise policy for an executive stock option and its implications for option cost, average life, and alternative valuation concepts. The paper is the first to provide analytical results for an executive with general concave utility. Wealthier or less risk-averse executives exercise later and create greater option cost. However, option cost can decline with volatility. We show when there exists a single exercise boundary, yet demonstrate the possibility of a split continuation region. We also show that, for constant relative risk averse utility, the option value does not converge to the Black and Scholes value as the correlation between the stock and the market portfolio converges to one. We compare our model's option cost with the modified Black and Scholes approximation typically used in practice and show that the approximation error can be large or small, positive or negative, depending on firm characteristics.     

Truncation and acceleration of the Tian tree for the pricing of American put options

We present a new method for truncating binomial trees based on using a tolerance to control truncation errors and apply it to the Tian tree together with acceleration techniques of smoothing and Richardson extrapolation. For both the current (based on standard deviations) and the new (based on tolerance) truncation methods, we test different truncation criteria, levels and replacement values to obtain the best combination for each required level of accuracy. We also provide numerical results demonstrating that the new method can be 50% faster than previously presented methods when pricing American put options in the Black–Scholes model.