Pricing Bermudan options by nonparametric regression: optimal rates of convergence for lower estimates

The problem of pricing Bermudan options using simulations and nonparametric regression is considered. We derive optimal nonasymptotic bounds for the low biased estimate based on a suboptimal stopping rule constructed from some estimates of the optimal continuation values. These estimates may be of different nature, local or global, with the only requirement being that the deviations of these estimates from the true continuation values can be uniformly bounded in probability. As an illustration, we discuss a class of local polynomial estimates which, under some regularity conditions, yield continuation values estimates possessing the required property.     

Pricing methods for α-quantile and perpetual early exercise options based on Spitzer identities

We present new numerical schemes for pricing perpetual Bermudan and American options as well as α-quantile options. This includes a new direct calculation of the optimal exercise boundary for early-exercise options. Our approach is based on the Spitzer identities for general Lévy processes and on the Wiener–Hopf method. Our direct calculation of the price of α-quantile options combines for the first time the Dassios–Port–Wendel identity and the Spitzer identities for the extrema of processes. Our results show that the new pricing methods provide excellent error convergence with respect to computational time when implemented with a range of Lévy processes.     

A perturbative approach to Bermudan options pricing with applications

In this paper we address the problem of the valuation of Bermudan option derivatives in the framework of multi-factor interest rate models. We propose a solution in which the exercise decision entails a properly defined series expansion. The method allows for the fast computation of both a lower and an upper bound for the option price, and a tight control of its accuracy, for a generic Markovian interest rate model. In particular, we show detailed computations in the case of the Bond Market Model. As examples we consider the case of a zero coupon Bermudan option and a coupon bearing Bermudan option; in order to demonstrate the wide applicability of the proposed methodology we also consider the case of a last generation payoff, a Bermudan option on a CMS spread bond.     

Pricing equity options everywhere

Finite difference methods are a popular technique for pricing American options. Since their introduction to finance by Brennan and Schwartz their use has spread from vanilla calls and puts on one stock to path-dependent and exotic options on multiple assets. Despite the breadth of the problems they have been applied to, and the increased sophistication of some of the newer techniques, most approaches to pricing equity options have not adequately addressed the issues of unbounded computational domains and divergent diffusion coefficients. In this article it is shown that these two problems are related and can be overcome using multiple grids. This new technique allows options to be priced for all values of the underlying, and is illustrated using standard put options and the call on the maximum of two stocks. For the latter contract, I also derive a characterization of the asymptotic continuation region in terms of a one-dimensional option pricing problem, and give analytic formulae for the perpetual case.     

An Explicit Finite Difference Approach to the Pricing Problems of Perpetual Bermudan Options

The pricing problem of options with an early exercise feature, such as American options, is one of the important topics in mathematical finance. Pricing formulas for options with the early exercise feature, however, are not easy to obtain and the numerical methods are thus frequently required to derive the price of these options. The value function of perpetual Bermudan options is characterized with the partial differential equation and this is solved by the finite difference method in this article.     

Pricing and hedging power options

In this paper we study the pricing and hedging of options whose payoff is a polynomial function of the underlying price at expiration; so-called ‘power options’. Working in the well-known Black and Scholes (1973) framework we derive closed-form formulas for the prices of general power calls and puts. Parabola options are studied as a special case. Power options can be hedged by statically combining ordinary options in such a way that their payoffs form a piecewise linear function which approximates the power option's payoff. Traditional delta hedging may subsequently be used to reduce any residual risk.     

Pricing options on scenario trees

We examine valuation procedures that can be applied to incorporate options in scenario-based portfolio optimization models. Stochastic programming models use discrete scenarios to represent the stochastic evolution of asset prices. At issue is the adoption of suitable procedures to price options on the basis of the postulated discrete distributions of asset prices so as to ensure internally consistent portfolio optimization models. We adapt and implement two methods to price European options in accordance with discrete distributions represented by scenario trees and assess their performance with numerical tests. We consider features of option prices that are observed in practice. We find that asymmetries and/or leptokurtic features in the distribution of the underlying materially affect option prices; we quantify the impact of higher moments (skewness and excess kurtosis) on option prices. We demonstrate through empirical tests using market prices of the S&P500 stock index and options on the index that the proposed procedures consistently approximate the observed prices of options under different market regimes, especially for deep out-of-the-money options.     

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.     

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.     

Pricing options on mean reverting underliers

For mean reverting base probabilities, option pricing models are developed, using an explicit measure change induced by the selection of a terminal time and a terminal random variable. The models employed are the square root process and an OU equation driven by centred variance gamma shocks. VIX options are calibrated using the square root process. The OU equation driven by centred variance gamma shocks is applied in pricing options on the ratio of the stock price for J. P. Morgan Chase (JPM) to the Exchange Traded Fund for the financial sector with ticker XLF. For the purposes of calibrating the ratio option pricing model to market data, we indirectly infer the prices for stock options on JPM from the prices for options on the ratio, by hedging the conditional value of JPM options given XLF, using options on XLF. The implied volatilities for the options on the ratio are then indirectly observed to be fairly flat. This suggests that for JPM, the use XLF as a benchmark is a possibly good choice. It is shown to perform better than the use of the S&P 500 index. Furthermore, though the use of an unrelated stock price like Johnson and Johnson as a benchmark for JPM provides as a good fit as does the use of XLF, this comes at the cost of requiring a considerable smile for the implied volatilities on the ratio options and hence a more complex model for the implied distribution on the ratio.     

Pricing Asian options in a semimartingale model

In this paper we studyy arithmetic Asian options when the underlying stock is driven by special semimartingale processes. We show that the inherently path dependent problem of pricing Asian options can be transformed into a problem without path dependence in the payoff function. We also show that the price is driven by a process with independent increments, Levy processes being a special case. This approach applies for both discretely or continuously options.     

Pricing Black-Scholes options with correlated credit risk

This paper presents an improved method of pricing vulnerable Black-Scholes options under assumptions which are appropriate in many business situations. An analytic pricing formula is derived which allows not only for correlation between the option's underlying asset and the credit risk of the counterparty, but also for the option writer to have other liabilities. Further, the proportion of nominal claims paid out in default is endogenous to the model and is based on the terminal value of the assets of the counterparty and the amount of other equally ranking claims. Numerical examples compare the results of this model with those of other pricing formulas based on alternative assumptions, and illustrate how the model can be calibrated using market data.     

Pricing exotic options in a path integral approach

In the framework of the Black–Scholes–Merton model of financial derivatives, a path integral approach to option pricing is presented. A general formula to price European path-dependent options on multidimensional assets is obtained and implemented by means of various flexible and efficient algorithms. As an example, we detail the cases of Asian, barrier knock out, reverse cliquet and basket call options, evaluating prices and Greeks. The numerical results are compared with those obtained with other procedures used in quantitative finance and found to be in good agreement. In particular, when pricing at-the-money and out-of-the-money options, the path integral approach exhibits competitive performances.     

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.