Popular Methods for Pricing American Options

Popular Methods for Pricing Ameri bum OptionsChapter 1 IntroductionAmerican alternatives atomic number 18 financial derivatives, an instrument whose nurture is derived from an be asset, usually a stock. bare and Scholes (1973) exposit an pick as a security giving the right to buy or sell an asset, subject to sealed conditions, within a specified period of m.The main question of this speech is how American preferences can be valued. The option value is still whapn with proof when the option is exercised, either at maturity or not. When the owner decides to exercise the option or it is the option maturity time, it is possible to determine the determine of the option as the strike will be exchanged by the asset in the case that the conditions be favourable for the owner of the option. When the one buys the option, she does not roll in the hay what will be the future toll of the rudimentary asset, and assuming it follows a random transit it is hard to put a price on much(prenominal)(prenominal) contract without knowing what will be the price change. This non linear feature of the option makes calculating the price to pay for such contracts a challenging process and has been the focus of a large number of financial studies and publications.This dissertation deals with the nigh popular methods for determine American options and their implementation in MatLab, including a graphic exploiter interface.The methods studied include the Black and Scholes (1973) European option price as the starting point, followed by the Barone Adesi and Whaley (1987) analytical approximation. Then the binomial and trinomial fretwork methods presented in Cox, Ross and Rubinstein (1979) are considered in like manner as the Finite difference approximations positions AAA. The most sophisticated method is the Least Squares Monte Carlo pretension presented in Longstaff and Schwartz (2001).The analysis of the diametric option price methods in this dissertation follo w most of the assumptions made by Black and Scholes (1973), the short term interest rate and the dividend are faux to be cognize and constant, the cardinal stock follows a log normal distributed geometric Brownian motion, the markets are frictionless and finally it exists the possibility of forming a riskless portfolio, consisting of the option and underlying stock.The dissertation is organised as follows a brief literature survey is provided in the near Chapter. The analytical approximation method and the numerical methods employ are described on Chapter 3 and their implementation in Matlab environment is go forthn in chapter 4. Numerical results are given in Chapter 5. The conclusion and future developments are presented in Chapter 6.Chapter 2 provides a survey of some of the most relevant publications in American Option Pricing, with focus on analytical approximations, grillwork and finite difference methods, more precisely, binomial and trinomial trees, explicit, implicit and Crank Nicolson Scheme, and in like manner on Monte Carlo Simulation.Chapter 3 provides a description of the methods used, their advantages, disadvantages and limitations. Here the required equivalences will be derived and the theme for the determine of American options will be provided.Chapter 4 focus on the algorithms used and their implementation on the MatLab environment, besides as the procedures for the development of the GUI for easier exploiter interface.On Chapter 5 results and their comparison are shown for the different methods used, with the required figures to support the numerical answers.In the final chapter the dissertation is concluded and a summary of the findings is provided, also as with further work on this subject.Chapter 2 Literature SurveyBlack and Scholes (1973) and Merton (1973) developed the eldest analytical closed form solution for the pricing of European type options and certain types of American options, such as American call options on non dividend paying stocks. The option pricing fashion model developed by Black and Scholes and extended by Merton gives rise to partial derivative instrument equations governing the value of an option Schwartz (1976).Black and Scholes (1973) develop their model on the basis of the no arbitrage theory, If options are correctly priced in the market, it should not be possible to make sure profits by creating portfolios of long and short positions in options and their underlying stocks Black and Scholes (1973).The Black and Scholes (1973) model valued European options on non dividend paying stocks, and with a number of quite restrictive assumptions, constant and cognize interest rates, the markets are frictionless with no transaction costs and penalties for short selling. The Black and Scholes (1973) model also assumes that the underlying stocks follow a random walk. Due to all this assumptions the pricing model Black and Scholes (1973) proposed was of easy use, and there is only the n eed to input the required values on the proposed pricing equation. The model they have proposed does not sign on into consideration early on exercise of the option so it is inaccurate for pricing American Options.One of the most popular analytical approximation models that starts from the Black and Scholes (1973) model and adjusts it to consider the scenario of early exercise strategies is the work by Baron Adesi and Whaley (1987) which was based on the paper by MacMillan (1986).Baron Adesi and Whaley (1987) consider that the Black and Scholes (1973) partial differential equation must apply to the early exercise premium as this is just the difference between the American and the European option prices, which are also priced by the same partial differential equation. After some transformation they end with an easily solvable through an interactive process second outrank differential equation.When closed form solutions, like the Black and Scholes (1973) valuation model cannot be d erived, numerical methods must be developed. These are computational methods where the values for the underlying assets are modelled up to maturity and the price of the options is derived from them. In the case of American options this is a complex process, as the modelled price changes may have to be set to include dividend payments and the derivation of the option price must also include the possibility of early exercise.Cox, Ross and Rubinstein (1979) developed a simple discrete time lattice model to deal with the complexity of option valuation, as they considered the methods of Black and Scholes (1973) quite advanced and have tended to obscure the underlying economics Cos, Ross and Rubinstein (1979). The use of lattice models such as the one by Cox, Ross and Rubinstein (1979) is the simplicity of its application.The most significant drawback of the Cox, Ross and Rubinstein (1979) model, is to increase its accuracy the number of time intervals must increase, in order to approach a continuous time model, which will significantly increase the computational time, needed for affect the entire tree in order to derive the option value.Others such as Hull and White (1988), (1993) and Trigeorgis (1991) have extended the model of Cox, Ross and Rubinstein (1979).Hull and White (1988) present a theater of operations of the use of lattice models for underlying assets with known dividends instead of known divided yields. They also consider the use of a control variate to price a option numerically, by a the lattice model, exploitation the price of a similar option calculated analytically. While Trigeorgis (1991) proposes a log transformed variation of binomial option pricing designed to overcome problems of consistency, stability and efficiency encountered in the Cox, Ross and Rubinstein (1979) foc apply on the pricing of foreign options. Hull and White (1993) also present an application of binomial and trinomial procedures for exotic path dependent options, where t hey developed a model faster than Monte Carlo simulation and faster than other numerical methods. commonly the analytical procedures are applicable to simple payoffs of the American Options, but in the cases where this is not possible numerical solutions must be developed. Geske and Shastri (1985) give a detailed comparison of the lattice methods to the different numerical methods, finite difference methods and other simulation methods.The model proposed by Brennan and Schwartz (1978) for valuing options was the first approach that used the finite difference method. This approach was used callable to the fact that most of the times an analytical solution for the option pricing problem does not exist. The finite difference method uses the heat equation derived from the Black and Sholes PDE to obtain an approximation of the option price. Courtadon (1998) goes further to reduce the approximation error of the Brennan and Schwartz (1978) model but only applies his findings only to simpl e option pay offs.Geske and Shastri (1985) give a good description of the finite difference method The finite difference technique analyze the partial differential equation () by using discrete themes of the changes in the options value for small changes in time or the underlying stock price to form equations as approximations to the continuous partial derivatives. Usually the approximations is done using frontwards, backward or central difference theorem, which respectively result in the explicit, implicit and Crank Nicolson schemes, the procedure used in this study will be shown further in the paper.In this case as with most of the methods for pricing options, the most significant drawback is the duality between accuracy and processing time. In order to increase accuracy the time and stock change steps must be smaller, increasing their number and the number of computations to make, this issue also affects the stability and convergence of the methods.Another approach used for solv ing the option pricing problem, especially for path dependent American options is the use of simulation. This means that the option price is derived from a simulated underlying asset price, usually using a Monte Carlo simulation method. Boyle (1977) and Schwartz (1977) pioneered the use of Monte Carlo simulation which is nowadays used to price complex options contracts. The Monte Carlo simulation method is very powerful in terms of its flexibility to generate the returns of the underlying asset of the options, by ever-changing the random variables used to generate the process a new returns distribution may be easily obtained, Boyle (1977).Boyle (1977) introduces the Monte Carlo technique for pricing European option where there is a dividend payment, but Schwartz (1977) was the true pioneer, pricing American options, with the underlying asset paying discrete dividends, and also deriving an optimal strategy for early exercise of the option, which is the crucial point for pricing Ameri can type options. Schwartz (1997) focused on a particular type of contract, warrants, so in fairness his first model is not exactly on an American type option.Tilley (1993) was one of the first to fully focus on the pricing of American option using a Monte Carlo simulation method as he mentioned that simulation methods were reserved for exotic options or other complex debt products. His findings are only utilize to American options on non dividend paying stocks, but he develops an important part of the model which is the optimal early exercise option.Carriere (1996) presents a development of the Monte Carlo simulation method presented by Tilley (1993). The paper by Carriere (1996) presents a model where the optima early exercise strategy is based on conditional expectations of Markov processes by carrying a nonparametric regression on the simulated underlying asset return paths.Brodie and Glasserman (1997) extended the previous studies by considering an upper and lower converging move of the option price. These estimated bounds are calculated using a high and a low bias, which Combining the two estimators yields a confidence interval for the true price. Brodie and Glasserman (1997)One of the most important papers, and probably one of the most used ones, is the paper by Longstaff Schwartz (2001). Their Least Squares Monte Carlo (LSM) valuation model is very simple and straight forward which combined with the accuracy of the method made it famous. Their big(p)est advance can be described as The key to this approach is the use of least squares to estimate the conditional expected payoff to the option holder from continuation Longstaff Schwartz (2001). They applied their model to a series of exotic path dependent American options with great success.Chapter 3 Pricing American Options Methods3.1 Asset Prices ModelsThe Black and Scholes (1973) and Merton(1973) pricing methods which are the basis for most of this paper assume that the stock returns follow a geo metric Brownian motions, with the stock prices log normally distributed.The stock returns can be represent by the undermentioned stochastic differential equation,(3.1.1)Where St is the asset price at time t, is the assets expected return, is the assets instantaneous volatility and Wt is a Wiener process.3.2 Analytical Approximation by Barone Adesi and Whaley (1987)Barone Adesi and Whaley (1987) developed a method to approximate analytically and easily the price of American options. They considered that the American and European option pricing equation is represented by the partial differential equation (3.2.1) developed by Black and Scholes (1987) and Merton (1987),(3.2.1)Barone Adesi and Whaley (1987) assumed that if this is true, then the early exercise premium of the American option, which is the price difference between the American and the European call option prices (3.2.2), can be represented by the same partial differential equation (3.2.3).(3.2.2)(3.2.3)The above equation after some transformation, shown on Barone Adesi and Whaley (1987) paper, and applying an approximation of a term tending to zero, yields the following quadratic equation,(3.2.4)Where (3.2.5), (3.2.6) and (3.2.7). Equation (3.2.4) is a second order ordinary differential equation with two linearly independent solutions of the form . They can be prove by substituting (3.2.8) into equation (3.2.4) Barone Adesi and Whaley (1987),(3.2.9)With a general solution of the form, (3.2.10)When the American option boundary conditions are applied to the above solution and considering , then must be equal to 0 as when the asset price tends to zero so does the option price, resulting in the following American call option pricing equation, Barone Adesi and Whaley (1987),(3.2.11)From (3.2.9) we have the value for so the only value missing is . This can be calculated interactively considering another boundary condition of American call options. We know that in early exercise the payoff will never be higher than S X, so from a critical underlying asset value the option payoff curve must be tangent to the S X curve, which means that below the critical asset value the pricing equation is represented by (3.2.11), Barone Adesi and Whaley (1987).The algorithm presented by Barone Adesi and Whaley (1987) for the above pricing problem is presented further in the paper in the section dedicated to the implementation of the American option pricing models.3.3 Lattice MethodsCox, Ross and Rubinstein (1979) proposed a model where the underlying asset would go up or down from one time step to the next by a certain proportional amount and with a certain probability until maturity. Due to the up and down characteristic of the asset price model these type of models are characterised by a binomial tree or, in the cases of the existence of a third possible movement, they are characterised by a trinomial tree, therefrom named as Binomial or Trinomial modelsThe price of the option would be recursi vely derived from maturity, due to the boundary condition as has been referenced before that the price of the option is only known with certainty at maturity.This means that the price of the option is calculated at maturity and recursively at to each one node up to the initial value, by discounting backwards at the risk free rate and respective probabilities. Due to the characteristic of American options, the model has to check if it is optimal to exercise the option at each node or if it has the advantage to continue to the next one, for example on the case of dividend payments.In the case that it is optimal to exercise the option at a certain node, its price will be equal to the intrinsic value at that same node. Every node will be check over for the optimality of exercising the option or not, until we have reached the initial point where we want to price the option.3.3.1 Binomial Tree ModelThe model starts being built for a American option of a non dividend paying stock and aft er that the scenario of dividend payments and optimal early exercise strategy is considered.As referenced before the stock goes up and down by a certain amount form one period to the next, if u is the up movement and d the down movement, then they can be calculated as, (3.3.1.1) and (3.3.1.2) as in Cox, Ross and Rubinstein (1979). In no arbitrage conditions it is possible to calculate the probability of the up and down movements, with the up being specify as, (3.3.1.3) where from the definition of probability and the down movement as (3.3.1.4).The tree formed using these specifications from Cox, Ross and Rubinstein (1979), can have the following graphical representationThe option is price is calculated from the asset price binomial tree. The maturity boundary condition for an American option, is that the payoff is equal to , we already have S at each maturity node from the asset price model, so we can calculate backwards the price of the option as the expectation of the future payo ff of the option.At each node we calculate the expectation of the future payoffs, where the price of the option will be a compound of expectations. These can be represented by the multi period case for a call as in Cox, Ross and Rubinstein (1979),The option prices are calculated as the expectation of the options future payoffs using their respective leaden risk neutral probabilities of an up movement and a down movement and then discounted at the risk free rate r. The Binomial value is fix for each node, starting at the final time step, and working backwards to the