# Pricing of American Options using Simulation Models

**Info:** 12822 words (51 pages) Dissertation

**Published:** 16th Dec 2019

**Tagged:**
Finance

__Abstract__

**Notations -----------------------------------------------------------------------------------------**4

- Introduction to Financial Derivatives ------------------------------------------------------5
- European and American Options-------------------------------------------------5-6
- Options Position: Long and short------------------------------------------------ 6-7
- Factors affecting Option prices---------------------------------------------------7-9
- Put-Call Parity for American Options-------------------------------------------9-10
- Bounds on Options ---------------------------------------------------------------10-12
- Arbitrage on American Options --------------------------------------------------12

- Pricing Approaches to American Options ------------------------------------------------13
- Numerical Pricing Model
- Binomial Model --------------------------------------------------13-15
- Trinomial Tree Model--------------------------------------------15-16
- Risk-Neutral Valuation------------------------------------------16-17
- Monte Carlo-------------------------------------------------------17-18

- Numerical Pricing Model

- Black Scholes Equations--------------------------------------------------------18-19
- Partial Differential Equations (PDE)------------------------19-20
- American Puts--------------------------------------------------20-21

- Monte Carlo Simulation---------------------------------------------------------------------26-28
- Least Square Monte Carlo Method (LSM)------------------------------------28
- Quasi Monte Carlo Methods---------------------------------------------------29-30

- Comparisons of Pricing Models-----------------------------------------------------------30-31

- Conclusion -------------------------------------------------------------------------------------31-32

- Acknowledgements ----------------------------------------------------------------------------33

- References---------------------------------------------------------------------------------------34-36

__Notations:__

- t – time

- T - the expiry (expiration date)

- E - the exercise price of the option

- K-Strike price

- S - market price of the underlying asset; S = S(t)

- V- the value of the option; V = V (S, t)

- C - the value of a call option; C = C (S, t)

- P - the value of a put option; P = P (S, t)

- r- the interest rate

- σ- the volatility of the underlying asset (measure of the standard deviation of the returns of the asset)

- μ- the drift of the underlying asset (measure of the average rate of growth of the asset)

- D- Dividend rate

*“The modelling of financial derivative products is a fast-growing area of applied mathematics with ‘real-world’ applications to problems originating from the industrial revolution.”*Such instruments were developed to control the spread of risk caused by adverse changes in the market. They also provide the opportunity to make significant profit or loss for those prepared to accept risk. As we will be mainly discussing pricing methods for American options, it is rightly to include this quote from the New Palgrave Dictionary of Economics where Ross (1987) says:

*This does not mean, however, that there are no important gaps in the (option pricing) theory. Perhaps of most importance, beyond numerical results…, very little is known about most American options which expire in finite time…Despite such gaps, when judged by its ability to explain the empirical data, option pricing theory is the most successful theory not only in finance, but in all of economics.*The structure of this dissertation will be split into different chapters. In Chapter 1, we would give a brief introduction on Options and find out about the factors that affect them. We would also look into the various positions and types of options mainly focusing on American options. Chapter 2 would then introduce us to the various pricing approaches and their validity to the modern trading structure while again emphasising on American options. In Chapter 3, we would further discuss the Monte Carlo Simulation, its properties as well as drawbacks. Chapter 4 will provide us with some insight concerning numerical solutions for American puts. Finally, in Chapter 5, we would attempt to compare and differentiate between the different approaches. We would then conclude this thesis by summarising what we have achieved throughout our research and give a final verdict on pricing American options.

- Options Contract: European and American

- What should be the correct price of this option? i.e., what amount of premium would the writer be ready to pay?

- Will the writer be willing to undertake the risk associated with his obligation?

- Options Position: Long and Short

*There are two sides to every option contract. On one side is the holder who has taken a long position (i.e., has bought the option) while on the other side is the writer who has taken a short position (i.e., has sold the option). There are four types of option positions:*

- A long position in a call option.

- A short position in a call option.

- A long position in a put option.

- A short position in a put option.

- Factors affecting option prices

Factors affecting Option Prices |
Notation |

Price of underlying asset/stock | S |

Strike price of the option | K |

Time to Maturity | T-t |

Volatility | σ |

Riskless interest rate | R |

Present value of dividends from t-T | D |

Options |
European |
American |
||

Variable |
Call |
Put |
Call |
Put |

Stock Price | + | - | + | - |

Strike Price | - | + | - | + |

Time to maturity | ? | ? | + | + |

Volatility | + | + | + | + |

Riskless Interest rate | + | - | + | - |

Dividends | - | + | - | + |

- Put-Call parity on American options

*.*The exercise policy of the option holder is necessarily based on the information accumulated to date and not on the future prices of the stock. American put-call parity is a pair of inequalities for the same strike price options of the same maturity as shown in equations below: P+S-D≤C+X (1.1) C+Xe-r(T-t)≤P+S (1.2), where P stands for Put options and C for call options. If equation 1.1 does not hold, then P+S-D>C+X . In other words, P+S>D+C+X . Therefore, in order to generate an arbitrage profit, one would need to short the put and the stock and use more resources to buy the call, invest in the strike, and invest in the present value of the dividends D. If equation 1.2 does not hold, then C+Xe-r(T-t)>P+S . In such a situation, we would sell the call, borrow the present value (PV) of the strike, buy the put and the stock, and reduce the surplus. Per se, we cannot derive an equality for the American put-call parity as half of the arbitrage argument would break down. For example, if the call option ‘C’ is very high as compared to the put option ‘P’ in (1.1), then we cannot generate an arbitrage by shorting the call, borrowing the strike, borrowing the present value of the dividend, and using the proceeds to buy the put and the stock. The reason is that if we are assigned an exercise on the short call, although we have the stock to deliver, the receipt of the strike combined with liquidation of the put will not necessarily cover your borrowing. A similar argument would apply to equation 1.2

*.*Let’s assume that the equation holds with strict inequality and we short the put and stock, using the proceeds to buy a call plus invest the present value of the strike. If we are assigned an exercise on the short put, then the receipt of the stock can be used to cover your short stock, but the present value of the strike together with the liquidation value of the call will not necessarily suffice to pay the strike. Furthermore, we would not be able to pay the dividends on the short stock in the interim.

- Bounds on Option Prices

*“If this relationship was not true, then an arbitrageur could make a risk-less profit by buying the underlying asset and selling the call option.”*There is a lower bound for the price of a European call option which is derived by considering two alternative investment strategies:

- Buy a share today at a price St, or

- Buy a European call option on the share at a price of CEtoday and, at the same time, deposit enough money at the risk-free interest rate to yield the exercise cost Kat time T. The sum to deposit now is the present value of K, PV(K)which is known from the payoff equation: PV(K)=Ke-r(T-t).

**(1.3) The worst that can happen to an option is that it expires worthless, hence its value cannot be negative. This gives the following equation,**

**CE≥max(St-PV(K), 0) (1.4) Similarly, put options have the following price bounds PE≥0 and PA≥0 (1.5) PE≤K and PA≤K (1.6) where PE is the value of a European put option and PA is the value of an American put option, respectively. The bounds imply that a put cannot have a negative value (1.5), and the price of a put cannot exceed the exercise price (1.6). Similar to the call equation (1.3), a European put option has a lower bound given as PE≥PV K-St**

**(1.7) Also, similar to (1.4), the following holds: PE≥ max(PV(K)-St,0)**

**(1.8) In the case of early exercise, CA ≥ CE PA≥PE since American options provide more exercise opportunities than the corresponding European options. It follows from (1.4) and (1.8) that CE≥St-PV(K)**

**PE≥PV (K)–St Given r>0 , it follows that CE>St-K . It is thus never best to exercise an American call option on a non-dividend-paying stock before the expiry date. (If it was best to exercise early, CE would be equal St-K ). However, it may be optimal to exercise an American put option on a non-dividend-paying stock early. At any given time during a contract, a put option should always be exercised early if it is sufficiently deeply in the money. For an American put option, the stronger condition PE≥ K-St must always hold since immediate exercise is always possible. In general, the early exercise of a put option becomes more attractive as St decreases, as r increases and as the volatility decreases. If a put-call parity doesn’t hold up, then this gives rise to arbitrage opportunities.**

- Arbitrage on American Options

*.*Due to this convention, we may assume that an American call or put option should be exercised by its holder either prior to or at maturity date. Chapter 2: Pricing Approaches to American Options There have been many approaches put forward in order to price financial derivatives. While it is easy to price European Options, it is not the case for American options due to its ability to exercise any time before and on the maturity date. In this chapter, we will be reviewing the various pricing approaches put forward. While some of these approaches are purely numerical in nature, others are more or less analytical. 2.1 Numerical Pricing Methods 2.1.1. Binomial option pricing Simple yet powerful, such is the binomial option pricing method which is used to solve a lot of complex option pricing problems. It was first introduced in 1979 by John C. Cox, Stephen Ross and Mark E. Rubinstein. As such the binomial option pricing method is often denoted as the CRR model. The plus point about this model is that it covers an underlying asset over a period of time rather than just a point. Because of this, this model mostly useful to American options. Also, this method also takes into account the effects of dividends paid during the life of an option. It uses an iterative procedure, allowing for the specification of nodes, or points in time, during the time period between the valuation date and the option's expiry date. The model reduces possibilities of price changes as well as removes the possibility for arbitrage.

*Fig2.1: Two-step binomial model with period 1*

*The binomial tree consists of N time steps corresponding to times k=0, 1, . . . , N , and models an asset price Sk . The price is normalized to S0=1 . Every outgoing branch in the tree has the price either moving up to Sk+1=uSk or down to Sk+1=dSk where u=eσt*

**and**d=e-σt=1u , where u>1 If we define the ‘up’ probability p (the same at every node) then Sk is a discrete-time stochastic process. We let Fk be the σ -field generated by {S0,S1,…,Sk} . At time k the possible price values are specified by Sk=Sk0,…, Sk[k] with Sjj=ukd2j=d2j-k (2.1) (as such Sk0=uk, Skk=dk , as such the prices are listed in a decreasing order) For instance, there is a riskless savings account with £1 invested at time k and is worth £ R at time k+1 . The condition for no arbitrage is as such: d<r<u (2.2) The ‘up’ probability p is risk neutral if the discounted price Sk/Rk is a martingale. At time 0, this requires that: S0=1=1Rpu+1-pd=ES1R And it is easily seen that this is the same condition as that required at any node in the tree. Therefore, p=R-du-d We note that 0< p<1 if and only if (2.2) holds. Our interest lies more in identifying the node the price process is at, rather than the actual value Sk of the price. If we define Jk=12k+logSklogd (2.3) Then from (2.1) we see that if Sk=d2j-k then Jk=j , so Jk is the process followed by the index of the price. Since there is at each k , a 1-1 correspondence between Sk and Jk , the σ - fields these processes generate are the same: Fk=σS0,…,Sk= σ{J0,…, Jk} . In the classic American put option with strike K , the holder has the right to exercise at any time k=1,…,N and receives at that time the exercise value [K-Sk]+ . If Sk=Sk [j] then then exercise value will be [K-ukd2j]+ . The binomial option pricing makes the following assumptions:

- That there is a perfectly efficient market due to which it is able to provide a mathematical solution for an option at each point in the time frame specified.

- That the spot price can only go in two direction hence the name binomial at each node and that there can be only two possible prices (up and down).

- That the underlying asset does not pay any dividends at any point.

- That the rate of interest ( r) is constant throughout the life of the option.

- That there are no taxes and no transaction cost in motion on the market.

- That the risk factor is not really problem for the investors.

**(2.4)**In the above equation, μ is the measure of the average rate of growth of the underlying asset price which is also known as the drift and σ is the volatility of the asset, which measures the standard deviation of the returns while X=X(t) is the Brownian motion. Both the drift and the volatility are considered to be constant. The option payout function is u(S,t)

*.*A path dependent option is one for which u(S,t) depends on the entire path {S(t):0<t<t} ; whereas a simple (non-path dependent) option has uS,t=u(St,t) . For a European option the payout may only be collected at the final time so that it is f(T)=u(S(T),T) . For an American option, since we can exercise it at any time before T , the payout is fτ=u(Sτ,τ) in which τ is an optimally chosen stopping time. In some papers, Black and Scholes (1973) and Merton (1973) both described two methods of valuation of derivatives. The first one being risk-neutral valuation, which is applicable to path dependent options and other derivatives for which the PDE is either unavailable or intractable, is the formula below: FS,t=maxE'[e-rτ-tu(Sτ,τ)St=S]

**(2.5)**In the above equation, E' is the risk-neutral expectation, for which the growth rate μ in (2.4) is replaced by r . The maximum is taken all over all stopping times τ with t<τ<T . The risk-neutral method offers a near true tool to find the value of an option in the real work. The fact that many real projects deal with a combination of embedded options, this method make it easier to price and incorporate multiple options into analysis. The fact that the risk-neutral valuation model requires thousands of calculations in order to get the best approximate value of the option makes it a great setback. Also, this method requires that the current value and variance of the growth rate of the underlying asset. Another drawback is that risk is definitely not an asset, so this method cannot really be applied. There is much technological advances nowadays, so it is a belief that in the near future, risk-neutral techniques would be widely used as a tool to price real options. 2.1.4 Monte Carlo Monte-Carlo simulation is one of the most powerful and most easily implemented methods for the calculation of option values. However, the procedure is essentially only useful for European-style options, where no decisions are made until expiration. Monte Carlo has increasingly become more attractive as the problem of pricing has become even bigger.

*As such, there are several Monte Carlo simulation models are known and used in practice.*The Monte Carlo technique for valuing derivatives is mainly based on the probability distribution of complete histories of the underlying derivative. The difficulty that occurs when pricing American options through Monte Carlo valuation is that Monte Carlo requires options that depend on multiple underlying derivatives or those that

*involve path dependent features. Since the determination of the optimal exercise time depends on an average over future events, Monte Carlo simulation for an American option has a “Monte Carlo on Monte Carlo” feature that makes it computationally very complex.*For a Monte-Carlo simulation of an American option, one performs the same type of Monte-Carlo simulation as for European options but then takes the Sup for all times τ ∈ [t, T] to obtain a realization of Pt . Then by doing this N times and taking the average, one obtains Pt for one value of S . In practice this is quite expensive and one may prefer binomial trees (see above) or ﬁnite difference or ﬁnite element methods A disadvantage of this method is that suitable accuracy may require a very large number of simulation trials. In general, the expected error decreases with the number of trials n by factor 1/n ; so, one more digit of accuracy requires 100 times as many trials. Often tens of thousands of trials are required to obtain two-place accuracy. 2.1.4.1 The Feynman-Kac Theorem Named after Richard Feynman and Mark Kac, this theorem establishes a possible relationship between solutions of parabolic partial differential equations (PDEs) to an

*expectation which establishes a mathematical link between PDE formulation of the diffusion complications in finance and the Monte Carlo simulations.*Given the stochastic processes: dXi=bidt+∑j=1naijdWj For i =1…n , With formal solution XiT=Xit+∫tTbidt+∫tT∑j=1naijdWj

**,**Any function V(t,X) with boundary conditions VT,X=f(X) That satisfies the partial differential equation ∂V∂t+g+12∑i,j=1ncij∂2V∂Xi∂Xj+∑i=1nbi∂V∂XI=kV with cij∶= ∑k=1naikajk Can be represented as the expectation Vt,X=E[fXTe-∫tTkdu+∫tTge-∫tskduds

*Hereby, all the coefficients*aij, bi, k and g

*can be functions both of time*t

*and the state vector*X(t)

*.*

- Black-Scholes equation [Analytical]

*to calculate the theoretical value of European-style options using current stock prices, expected dividends, the option's strike price, expected interest rates, time to expiration and expected volatility.*The formula was formerly developed by three economists namely; Fischer Black, Myron Scholes and Robert Merton and was first introduced in their 1973 paper, "The Pricing of Options and Corporate Liabilities" published in the Journal of Political Economy. The Black-Scholes model makes the following assumptions:

- Only European options can be exercised at expiration
- During the life span of the option, no dividends is paid out.
- The need for efficient markets (i.e., market movements cannot be predicted)
- There are no transaction costs in buying the option
- The risk-free rate and volatility of the underlying are known and constant
- The returns on the underlying assets are normally distributed

- The underlying asset price follows the lognormal random walk: dSt=μSt+σStdWt.

- The price of the option
uis a continuous, differentiable function of the asset price
Stand time
t. i.e.,
u ≡ u(St, t)
**.**

- The risk-less rate of interest is a known function of time during the life of the option, investors can borrow and lend at this rate.

- The underlying asset pays no dividends during the life of the option.

**(2.6). The change in the value of the portfolio in time [0, t] is therefore ∏t-∏0= ∫0tduSτ,τ-∫0t∆τdSτ Note that Δt has not changed during the time-step. Ito’s lemma implies that the portfolio follows the random walk: d∏t=∂u∂tdt+∂u∂SdSt+12σ2St2∂2u∂S2dt-∆tdSt 2.2.2 American Put Option The American put option gives the owner the right to sell the underlying asset at a fixed strike price K∈R+ at any time before the expiry date of the option T>0 . Let’s assume that the spot price of the American put is St and strike price K and time expiration T . If the spot price is St* at time t* is 0, then St=0 holds for t≥t* since the price process follows a geometric Brownian motion. At this point, it is not recommended to exercise the option at a later date. If St*=0 , then the value of the put at t* is K , i.e. the initial value. The no-arbitrage condition demands that the price of the option is never less than the payoff value, g(S)=max(K-S, 0) . This condition leads to a partial differential inequality equation satisfied by the value of the American. To derive this, a riskless portfolio is constructed. The possibility of early exercise demands that the appreciation in the value of the portfolio is bounded above by the risk-free rate of interest. The pricing function u satisfies the inequality: -∂u∂t≥LBSu In Ω={S: S≥0} × [0, T] where LBS is the Black-Scholes differential operator of equation: -∂u∂t=12σ2S2∂2u∂S2+rS∂u∂S-ru=LBS u . The valuation of American options is a free boundary problem. Typically, at each time t there is a value of the spot which marks the boundary between two regions; to one side one should hold the option and to the other side one should exercise it. That is, the domain Ω is demarcated into two regions C and ℇ where C={(S, t)∈Ω:u(S, t)>g(S)} , is the continuation or holding region. Its complement is given by ℇ={(S, t)∈Ω:u(S, t)=g(S)} Which is the exercise region. If the option is exercised, then its value is determined by the payoff. If it is held then its value must be greater than the immediate payoff. Since the option is either exercised or held, the pricing function satisfies the linear complementarity problem given by, -∂uS,t∂t≥LBSuS,t∈R+×[o,T] uS,t≥gS∀S,t∈R+×[o,T] -∂u∂t-LBSu.u-g=0∀(S,t) ∈R+×[o,T] and final condition u(S, T)=g(S) . The complementarity formulation can be generalised to multi-factor financial derivatives. While it initially came into existence in 1973, the Black-Scholes formula was initially devised to price European options on stocks that do not pay any dividends. The main advantage of the Black-Scholes is that it allows one to calculate large amount of option prices in a very short time. Its major disadvantage however comes when pricing the American options. Black-Scholes cannot accurately price American option as the formula can only calculate the option at the expirations time. It ignores all possibility of an early exercise of the American option. However, a loophole can be found in the case of American calls where the American call is worth the same as a European call making it unsuitable for exercising it early. 2.3. Other Pricing Models 2.3.1 Optimal-stopping formulation In mathematical theory, optional stopping is one where the problem arises when choosing the time to exercise an option in order to maximise profits or minimize loss. Due to the fact that the American option can be exercised on or before its maturity date sort of makes it an optimal-stopping problem. Following a classic set up under the Black-Scholes. The stock price would normally follow a geometric Brownian motion: St=S0expr-δ-σ22t+σBt {\displaystyle r}{\displaystyle S_{t}=S_{0}\exp \left\{\left(r-\delta -{\frac {\sigma ^{2}}{2}}\right)t+\sigma B_{t}\right\}} Let us find the fair value of an American option (Vam) contract which is equivalent to solving the following optimal-stopping-time problem**

**Vam=supt∈T E[e-∫0τruduPτSτ] (2.1) In the above equation, the expectation is taken under the risk-neutral martingale measure. The instantaneous short rate process r(t) is used to discount the payoff Pτ at time τ , and the supremum in (2.1) is taken over all admissible stopping times τ ∈ T between now and the final expiry of the option. The rate r is constant. Let us consider the case of an American put option with strike price K on a single underlying S with price process St . Assuming that the process follows a geometric Brownian motion with interest rate r under the risk-neutral martingale measure, the equation 2.1 can be further converted into the following; supτ∈T E[e-τrmaxK-Sτ,0] (2.2) In case of a put option on a single underlying the hold and exercise regions are half planes divided by the optimal exercise boundary B*(t) and the supremum in (2.2) is attained for an optimal stopping time τ* of the form τ*=inft≥0:St≤B*t (2.3) Figure 2.1 illustrates the situation for a sample path which hits the optimal exercise boundary. Fig 2.1: hold and exercise region under a given exercise rule The difference in value between an American derivative and its European equivalent is referred to as the early exercise premium. In contrast to European options there are generally no closed-form solutions for American option even for plain vanilla instruments. The pricing of an American option is based on the assumption that a rational option holder will choose to exercise his right optimally in some sense. Thus, finding an optimal exercise strategy is at the core of pricing an American option. Mathematically, this translates into the optimal stopping formulation. 2.3.2 Free-boundary formulation Another plausible method on how we can possibly price the American options is the free boundary problem. Let us consider the case of American put option. The decision to whether exercise an option or to hold it for a little further depends naturally on the current of the underlying asset basically at every point in time. Mathematically, the (t, S) plane can be divided into a hold region (where it is optimal to hold the option) and an exercise region (where it is optimal to exercise it). The barrier between those two regions is called the optimal exercise boundary B*(t) (the blue line). Fig2.3: (t, S) plane showing the hold and exercise region. Let the fair value of the American put option be Pam . In this situation, the hold region lies above the exercise boundary B*(t) in the (t, S) plane. Within the hold region the value function Pam satisfies the Black-Scholes equation, i.e. LBSPam∶=∂Pam∂t+12σ2S2∂2Pam∂S2+rS∂Pam∂S-rPam=0 (2.4) where r and σ are the usual drift and volatility parameters from the geometric Brownian motion SDE. At expiry T and along the optimal exercise boundary B*(t) the option is worth its payoff, i.e. Pam(T, S)=max (K-S, 0) and Pam(t, B*(t))=max (K-B*(t), 0). The transition between the two regions is smooth, in the sense of the so-called smooth pasting condition ∂Pam∂S t,B*t=-1 One of the drawbacks of the optimal exercise boundary is that the boundary is not known but must in fact be computed as part of the solution, making the problem is much harder to solve than a simple boundary value problem. We note that the smooth pasting condition only holds for options whose payoff is sufficiently smooth. It does not hold in the case of an American digital put for example. Another point worth mentioning is that for more complicated American derivatives the hold region might not be simply connected at all times. Overall, free boundary problems are hard to solve numerically. 2.3.3 Linear complementarity formulation Another way to characterize the price of an American option is as the solution of a linear complementarity formulation. In this formulation, we are looking for a function V (S, t) such that LBSV≤0, VS,t≥P(S,t) VS,t-PS,tLBSV=0 And VS,T=P(S,T)**

**Here, P(S, t) denotes the payoff function of the option and the Black-Scholes operator LBS is defined as in (2.4). The linear complementarity formulation is the clever combination of several features of an American option. Some of these can easily be argued heuristically. We had already argued before that based on the no-arbitrage assumption the option must be worth no less than its immediate payoff, i.e. V(S, t)≥P(S, t) . Moreover, it is clear that the option value coincides with its payoff value whenever it is optimal to exercise. It can be shown that within the hold region of the option its price function V(S,t) satisfies the Black Scholes equation. Combining these observations one ends up with the linear complementarity formulation above. 2.3.3 Stopping rules and parametric approximations The connection in (2.3) allows us to interchangeably formulate the pricing problem either in terms of stopping times or in terms of exercise regions. Note that any choice of a stopping time τ̂∈T will produce a (usually suboptimal) value in the expectation from (2.1), i.e. Vam(τ̂)=Ee-τ̂rPτ̂Sτ̂≤Vam (2.5) If the optimal exercise region ε* is known, then the optimal stopping time τ* is the first time that St enters the exercise region τ*=infτ∈T{t≥0:St∈ε* } (2.6) and the fair value of the option can be computed as a simple expected value according to Vam(τ*)=Ee-τ*rPτ*Sτ*=Vam (2.7) The mapping implicit in (2.6) can be used to assign a stopping time to any given “exercise region”. As such there is a possibility in terms of optimizing over exercise regions (or boundaries) and find the value of an American option. The interesting consequence is that given an exercise strategy (i.e. an exercise region), a forward Monte Carlo simulation of the underlying path process can be used, compute the payoff for each path and take an average to find the value of the option under the chosen exercise strategy. The emphasis here is that most numerical pricing algorithms for American options including the famous Longstaff- Schwartz algorithm are based on backward iterations and that forward algorithms are usually unsuited for the problem. Suppose now that we can parametrize all feasible exercise regions εθ by some parameter θ ∈ Θ and denote by τ(θ) the stopping time associated with εθ according to the mapping presented in (2.6). The option value is then approximately equal to Ṽam=Supθ∈ΘE[e-τθrPτθSτθ] (2.8) however, the value computed in (2.8) is in general biased low because we might not hit all feasible stopping times in our approximation. In other words, the set τ(Θ)⊂T will in general be a proper subset of T and hence Ṽam=Supθ∈ΘEe-τθrPτθSτθ≤Supτ∈TE[e-τrPτSτ=Vam Nonetheless, equation (2.8) is the fundamental base of our numerical approach. For an American put option a lot is known about the optimal exercise boundary. One of its main characteristics is its asymptotic behaviour close to the expiry date of the option, where it clearly hits the strike. In fact, the optimal exercise boundary of an American put option asymptotically approaches its strike level K displaying a square root i.e. B*t~K-cT-t for T-t≪1 (2.9) with an appropriate constant c . Even more importantly, it influences our choice of a parametric family of feasible exercise boundaries, which we use in the optimization problem (2.8). Chapter 3: The Monte Carlo Simulation (MC) The Monte Carlo method is well suited once it is agreed that security pricing can be represented by expectations. Monte Carlo seeks to solve the pricing equation satisfied by the price differential equation by representing the solution as an expectation, then using numerical simulation to evaluate the expectation. The approach consists of the following steps:**

*Sample paths of the underlying state variables (underlying asset prices and interest rates) are**simulated over the relevant time frame.*

*The discounted cash flows of a security on each sample path, as determined by the structure of the security, is evaluated. Let’s call it*uj*.*

- The discounted cash flows are averaged over the sample paths. That is,

*method for pricing complex financial derivatives, such as contracts whose payoff depends on several correlated assets or on the entire sample path of an asset price. We can simulate as many sample paths as desired according to the underlying stochastic differential equation that describes the stock process. For each sample path, the option value is determined and the average from all paths is the estimated option price. The option price μ is written as an integral that represents the mathematical expectation of the discounted payoff under a so-called risk-neutral probability measure. This expectation is usually with respect to a non-uniform density over the real space, but with a change of variables, it can be rewritten as an integral over the s-dimensional unit hypercube:*

*0,1)t=u=u0,…, ut-1:0<uj<1 for all j: μ=μf=∫01…∫01fu0,…, ut-1du0…dut-1=∫[0,1)t=EfU, (3.1) For some function f:[0,1)t→R , where u represents a point in [0,1)t , and U~U[0,1)t is a random point with the uniform distribution over the unit hypercube. Let’s assume that the integral is pre-written as in (3.1) for a fixed non-negative integer t . As such, we would also want to estimate the value of μ . The integer t represents the number of calls to the underlying random number generator that would be used in our simulation. When the number of calls is random and limitless, we can denote t as an infinite while assuming that the probability is 1. Only a finite number of coordinates of U need to be plainly generated. In the case of European option, the MC method works fine. As a matter of fact, we even have an analytical solution in the form of the Black-Scholes formula. More importantly, the value is determined only by*

*the final stock price if one assumes a given starting point, time, constant interest rate and volatility.*It is much easier to see that the Monte Carlo simulation would work in a forward direction. However, for the American option, it might not be the case. Because of early exercise condition,

*in contrast to a partial differential equation, we would also need to know the option value at the in-between times between the simulation start time and the option expiry time. In Monte-Carlo this information is harder to obtain, therefore, even though it is simple and capable of handling multi-factor problems, we would*need to solve the problem by going backwards. As such the Monte Carlo simulation becomes difficult to implement. Like all the methods described above, the Monte Carlo simulation also comes with its advantages and disadvantages. To begin with, the probability distributions within the model can be easily implemented without the need to approximate them. Furthermore, the correlations of the model can be modelled with ease. The method itself is quite easy to implement. Moreover, the simulation can be easily mounted on various software packages such as Matlab and Excel. The main drawback of this simulation however is that it is more of an approximate way to price options than a sure way although it is possible to reduce the percentage of deviation by simply increasing the number of simulations. That said, this set new limitations to the model such as how much time will take for an analyst to generate the simulations. 3.1 Least Square Monte Carlo method (LSM) In 2001, Longstaff and Schwartz introduced the Monte Carlo simulation and least squares algorithm of Carriere in order to price American options. Option holders would compare the payoff between immediate exercise and exercise it at a later date (continuation date) depending on which one is higher. Else, the options would remain as it is. The expected payoff for a later exercise is conditional on the information available at that time point. The main emphasis for this approach is that the conditional expectation can be estimated from the cross-sectional information in the simulation by using least squares. This makes this approach readily applicable in path-dependent and multifactor situations where traditional finite difference techniques cannot be used. In order to find out the conditional expectation function, we revert back the realized payoffs from continuation on a set of basic functions in the underlying asset prices. The fitted values are chosen as the expected continuation values. We would compare these continuation values with the immediate exercise values and make the optimal exercise decisions. We would then obtain a complete specification of the optimal exercise strategy along each path. We can use this algorithm and discount the optimal payoffs to time zero. That would give the option price. The method starts with N random paths (Skn, tn) for 1≤k≤N and tn=ndt . Valuation is performed by rolling-back on these paths. Suppose that Fn+1k=FSkn+1, tn+1 is known. For points Stn, tn set X= Skn the current equity value and Y=e-rdtF(Skn+1,tn+1) the value of deferred exercise. Then perform regression of Y as a function of the polynomials X,X2,…, Xm for some small value of m which is called basic function; i.e. approximate Yk by a least squares fit of these polynomials in X . Hence we use this regressed value in deciding whether to exercise early. 3.2 Quasi-Monte Carlo Methods-QMC (low-discrepancy)

*Instead of generating sample paths randomly, it is possible to systematically select points in a probability spaces so as to optimally "fill up" the space. The selection of points is a low-discrepancy sequence. Taking averages of derivative payoffs at points in a low-discrepancy sequence is often more efficient than taking averages of payoffs at random points. Because of the multidimensional numerical integration problem, the quasi-Monte Carlo methods was proposed in lieu of the Monte Carlo methods. In order to resolve this problem, we would construct the following equation:*μ=∫[0,1)tf(u)du

**where f is a square-integral function. To approximate μ , both MC and QMC proceed by choosing a point set Pn=u0, ⋯, un-1∈[0,1)t , and then the average value of f over Pn is computed, i.e., we get the estimator Qn=1n∑i=0n-1f(ui) In Monte Carlo method, the points u0,…, un-1 are independent random vectors and uniformly distributed over [0,1)t and n is the number of replications (a constant). This estimator is unbiased and has variance σ2/n , and define σ2=VarfUi, where VarfUi=∫[0,1)t f2udu-μ2**

**(3.2)**If σ2<∞ , then Qn obeys a central-limit theorem, and we can rely on it to compute a confidence interval on μ , whose width converges roughly as Oσn-1/2. The use of MC for pricing financial options was first proposed by Boyle. Since then the approach have evolved quite greatly.

*In practice, one uses a pseudorandom number generator to choose these points. The idea of quasi-Monte Carlo methods is to use a more regularly distributed point set, so that a better sampling of the function can be achieved. An important difference with Monte Carlo is that the set*Pn

*is typically deterministic when a quasi-Monte Carlo method is applied. In other words, Quasi-Monte Carlo replaces the independent random points*Ui

*in (3.2) by a set of n deterministic points,*Pn=U0,⋯, Un-1,

*which cover the unit hypercube*[0,1]t

*more evenly (uniformly) than a typical set of random points. The point set*Pn

*is called a design by some statisticians. Niederreiter presents these methods in detail in his book, and describes different ways of measuring the quality of the point sets*Pn

*on which quasi-Monte Carlo methods rely. More specifically, the goal is to measure how far is the empirical distribution induced by*Pn

*from the uniform distribution over*[0,1]t

*. Such measures can be useful for providing upper bounds on the deterministic integration error*∣Qn–μ∣

*.*Typically, a point set Pn is called a low-discrepancy point set if D*(Pn)=O(n-1 logtn)

*.*For a function of bounded variation in the sense of Hardy and Krause, the integration error ∣Qn –μ∣is in O(n-1logtn) when Pn is a low-discrepancy point set The upper bound states that the advantage that QMC methods have over MC is that QMC has a probable error in O(n-1/2) which will ultimately be lost as the dimension t increases. Consequently, it suggests that quasi-Monte Carlo methods will require a sample size n too large, for practical purposes, to improve upon Monte Carlo when t is large. In this context, the programming codes for QMC simulation on several assets are too much complicated, then numerical results showing an improvement of QMC over MC in high dimensions and using a relatively small sample size n seem hard to explain

*.*Chapter 5: Comparisons of Pricing Models Throughout this thesis, we have reviewed the various pricing methods that can be used on American options. The pricing of European options is rather deemed as being easy and accurate due to the fact that they expire only at maturity and are hence done through the Black-Scholes equation. This is not the case for the American option due to its ability to be exercised any time before and on the expiration date. As such, numerous mathematicians and analyst have come forward to present multiple models to price the latter, the most famous being the Black-Scholes, Boundary formulation and the Monte Carlo simulation. First to be analysed was the binomial model. The binomial model is rather very easy and simple to use. With its various conditions like the price would move in only two directions; up and down, it can be assumed that it is very easy. However, the problem arises when we have to price multiple assets over periods more than 1, then it becomes much complicated and time-consuming to do so. The same would go for the trinomial model, where instead of two, the price would move in either three directions; up, down or no change. This model however is more adapted to real life situations. The binomial model calculates the option prices at each next step. As such if there are any adjustments to be done such as an early exercise of the American option, it can still be adjusted and worked through at that specific time. The binomial is among one of those who calculates by working backwards, till we are left with only one option price. It should be noted that in the case there is no alteration to the model and that the same inputs are used as in a Black-Scholes model, with the help of sufficient steps, both the Black-Scholes and the binomial model would look alike. The main difference between these two is that the binomial model offers a greater flexibility than the Black-Scholes as the inputs can be changed at each step in order to process for differences arising while exercising a particular option. Longstaff and Schwartz introduced their well-known pricing algorithm for American options in [Fa01] based on a dynamic programming formulation. The Longstaff- Schwartz algorithm uses a backward iteration where at every timestep the continuation value of the option is approximated with least squares fit against a set of basic functions. Most tree based methods to approximate the value of an American option are also based on backward iterations. Such methods are simple yet powerful pricing tools. With a Monte Carlo approach pricing and managing the risks of American derivatives is far more involved than pricing and managing the risks of European options. The added difficulty stems from the fact that it is a priori unclear when the option holder will choose to exercise the option. Among all the pricing methods that work backwards, the Monte Carlo is the only one to be working forward in order to price options. As such, it uses computerised modelling to predict outcomes. At first the model generates a random number based on a probability distribution. The number then uses additional inputs of volatility and time to expiration to generate a stock price. The latter is then used to calculate the value of the option. In turn the model will keep generating results each time using a different set of random values from the probability functions. Depending on what model of Monte Carlo we are using, the number of uncertainties and probability distribution used, the simulation can generate a thousand of recalculations before it is completed. Monte Carlo simulation generally would take the average of the all the results as the option price. Conclusion The purpose of this thesis was to attempt to price American options using various pricing approaches and to compare them to the Monte Carlo Simulation. Also, we concluded that among the many Monte Carlo approaches, using the Least Squares Monte Carlo method is more suitable in solving problems of bigger size. While the pricing of European options have been quite easy due to the fact that it can only be exercised at maturity time, it has been rather difficult to do the same with American options. In fact, over the years many mathematicians have devised many strategies for pricing the American option. Due to the fact that underlying dynamics can be very complicated, the analytical models for option pricing involve a lot of restrictive assumptions. Coming to real-world applications, we use numerical methods that can approximate options, notably for American options. These include the valuation of options, the estimation of their sensitivities, risk analysis, and stress testing of portfolios. However, in recent years numerical experiments have become a lot complicated in financial theory and practice has increased enormously. putting more demands on computational speed and efficiency, giving way to further develop and enhance the methods. We have seen as well how the risk-neutral valuation model can be a valuable tool in the future to accurately price options due to the rapid advancement in computational technology.

*The pricing of American options remains a very interesting research area, especially while using Monte Carlo techniques. The Monte Carlo method is best used in handling complex and big options. In the near future, we would be able to value American options with much ease as the simulation techniques of the Monte Carlo method are becoming broader and promising, particularly while dealing with multiple factors in the market. When options are traded, they are often done so in bulk rather than single. As such, to price each and every one of them can be quite tiring and time-consuming. As far as calculations are concerned, I personally feel the use of Binomial model would be a perfect tool for pricing as it is simple and easy to use. The only problem though would come with multi assets over multiple periods of time. As such one can hope to use the Monte Carlo to calculate heavy portfolios and assets Its ability to calculate hefty input at high speed makes it one of the most sort out model for option pricing. So much that there have been a lot of derivatives from the original simulation. In case of heavy inputs, I would say that the Monte Carlo model is the perfect model to price options However.in order to correctly execute a Monte Carlo simulation, one must have more than basic knowledge about unlike the binomial model which require very basic knowledge. Although, the Monte Carlo is becoming more a computerised system and is advancing increasingly in terms of technology. The world has become a very unstable financial and economical place. Over years, there have been some deep impacting events that have had damaging effects on the global financial health. Events such as the great recession, multiple terrorists’ attacks over the world, and the recent Brexit situation have cause major fluctuations. As such the market would be fluctuating a lot. One can only rely on the effectiveness of risk analysts and market watchers as to not further get in to financial debts. These analysts and watcher in turn depend on the reliability of such models to correctly price options and a make a reasonable gain from it and limit the loss to the minimum. As such, in my opinion, due to the fact that there are a lot of software being developed and these models being implemented on them, in the coming decade the pricing of American option is particular would be made quite accurate and quick. One such model would be the risk neutral valuation and the Monte Carlo simulation.*

__References__

**[And09]**A. Pascucci February 2009:

*“A short course on American options”*, web:http://www.dm.unibo.it/~pascucci/web/Ricerca/PDF/american.pdf

**[Boy77]**P. P Boyle. Options: A Monte Carlo approach. Journal of Financial Economics, 4:323–338, 1977.

**[Boy97]**P. P Boyle, M Broadie, and P Glasserman: Monte Carlo Methods for Security Pricing, pages 1267–1321. 1997

**[Carr98]**P. Carr, “

*Randomization and the American Put*,” 1998,

*Review of Financial Studies*, 11, 3, 597–626.

**[Cle02]**E, Clement. D, Lamberton. P, Potter: 2002,

*“An analysis of a least squares regression method for American option pricing”.*

**[Dam09]**D. lamberton September 2009: “Optimal stopping and American options”,

**[Dew93]**Dewynne, J. N., Howison, S. D., Rupf, I. and Wilmott, P.: 1993,

*“Some mathematical results in the pricing of American options”*, European J. Appl. Math. 4(4), 381–398.

**[Dup98]**B Dupire, editor.

*“Monte Carlo: Methodologies and Applications for Pricing and Risk Management”*. Risk Publications, London, 1998.

**[Fra01]**Francis A. Longstaff and Eduardo S. Schwartz. Valuing American options by simulation: A simple least-squares approach. Review of Financial Studies, 14(1):113–147, 2001.

**[Gla04]**Paul Glasserman:

*“Monte Carlo Methods in Financial Engineering”*. Springer, 2004.

**[Hir03]**A. Hirsa:

*“Why Be Backward: Forward Equations for American Options*

*”*, January 2003, 103–107

**[Hull03]**J Hull:

*“Options, Futures and Other Derivatives*”. Prentice Hall International, fifth edition, 2003.

**[Jar92]**R.Jarrow, R. Myneni: “

*Alternative Characterizations of American Put Options*,” (with R. Jarrow and R. Myneni), Mathematical Finance, April 1992, 87–105.

**[Joo07]**Joonhee Huh.:

*“Computation of Multivariate Barrier Crossing Probability, and Its Applications in Finance”*. PhD thesis, University of Waterloo, 2007.

**[Kar88]**Karatzas, I.: 1988,

*“On the pricing American options”*, Appl. Math. Optim. 17, 37–60.

**[Kie52]**J. Kiefer and J. Wolfowitz.:

*“Stochastic estimation of the maximum of a regression function”*, Annals of Mathematical Statistics, 23(3), 1952.

**[Kim90]**Kim, I. J.: 1990,

*“The Analytic Valuation of American Options”*, The Review of Financial Studies 3(4), 547–572.

**[Kir83]**S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi.:

*“Optimization by Simulated Annealing. Science”*, 220(4598):671–680, 1983.

**[Kus98]**Kuske, R. A. and Keller, J. B.: 1998,

*“Optimal exercise boundary for an American put option”*, Applied Mathematical Finance 5, 107–116

**[Mus98]**M, Musiela. and M, Rutkowski.:1998,

*“Martingale Methods in Financial Modelling*”, 19, 48-50

**[Qui09]**Q.Jia June 2009: “

*Pricing American Options using the Monte Carlo Methods”*web: http://www.diva-portal.org/smash/get/diva2:301061/FULLTEXT01.pdf

**[Rac98]**Rachel A. Kuske and Joseph B. Keller:

*“Optimal exercise boundary for an American put option”*. Applied Mathematical Finance, 5(2):107–116, 1998.

**[Rob51]**H. Robbins and S. Monro:

*“A Stochastic Approximation Method. The Annals of Mathematical Statistics”*, 22(3), 1951.

**[Rog02]**L. C. G. Rogers:

*“Monte Carlo valuation of American options”*. Mathematical Finance, 12:271–286, July 2002.

**[Rog97]**L.C.G. Rogers and D. Talay:

*“*

*American Options: A Comparison of Numerical Methods*

*”*, in Numerical Methods in Finance, ed. Cambridge University Press, 1997, 67–87.

**[Sin10]**Euan Sinclair:

*“Option Trading, Pricing and Volatility Strategies and Techniques”*, Wiley, 2010.

**[Tho00]**Thomas Little, Vijay Pant, and Chunli Hou:

*“A new integral representation of the early exercise boundary for American put options”*, Journal of Computational Finance, pages 73–96, 2000.

**[Wil93]**P Wilmott, J Dewynne, and S Howison:

*“Option Pricing: Mathematical Models and Computation”*, Oxford Financial Press, 1993

**[Wu97]**Wu, L. and Kwok, Y. K.: 1997:

*“A front-fixing finite difference method for the valuation of American options*”, Journal of Financial Engineering 6, 83–97. Web: http://fedc.wiwi.hu-berlin.de/xplore/tutorials/sfehtmlnode40.html

## Cite This Work

To export a reference to this article please select a referencing stye below:

## Related Services

View allRelated Content

All Tags**Content relating to: "Finance"**

Finance is a field of study involving matters of the management, and creation, of money and investments including the dynamics of assets and liabilities, under conditions of uncertainty and risk.

**Related Articles**

### DMCA / Removal Request

If you are the original writer of this dissertation and no longer wish to have your work published on the UKDiss.com website then please: