Merton s Jump-Diffusion Model
|
|
- Sylvia Patrick
- 6 years ago
- Views:
Transcription
1 Merton s Jump-Diffusion Model Empirically, stock returns tend to have fat tails, inconsistent with the Black-Scholes model s assumptions. Stochastic volatility and jump processes have been proposed to address this problem. Merton s jump-diffusion model is our focus. a a Merton (1976). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 698
2 Merton s Jump-Diffusion Model (continued) This model superimposes a jump component on a diffusion component. The diffusion component is the familiar geometric Brownian motion. The jump component is composed of lognormal jumps driven by a Poisson process. It models the sudden changes in the stock price because of the arrival of important new information. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 699
3 Merton s Jump-Diffusion Model (continued) Let S t be the stock price at time t. The risk-neutral jump-diffusion process for the stock price follows ds t S t =(r λ k) dt + σdw t + kdq t. (81) Above, σ denotes the volatility of the diffusion component. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 700
4 Merton s Jump-Diffusion Model (continued) The jump event is governed by a compound Poisson process q t with intensity λ, where k denotes the magnitude of the random jump. The distribution of k obeys ln(1 + k) N ( γ,δ 2) with mean k E (k) =e γ+δ2 /2 1. The model with λ = 0 reduces to the Black-Scholes model. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 701
5 Merton s Jump-Diffusion Model (continued) The solution to Eq. (81) on p. 700 is where S t = S 0 e (r λ k σ 2 /2) t+σw t U(n(t)), (82) U(n(t)) = n(t) i=0 (1 + k i ). k i is the magnitude of the ith jump with ln(1 + k i ) N(γ,δ 2 ). k 0 =0. n(t) is a Poisson process with intensity λ. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 702
6 Merton s Jump-Diffusion Model (concluded) Recall that n(t) denotes the number of jumps that occur up to time t. As k> 1, stock prices will stay positive. The geometric Brownian motion, the lognormal jumps, and the Poisson process are assumed to be independent. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 703
7 Tree for Merton s Jump-Diffusion Model a Define the S-logarithmic return of the stock price S as ln(s /S). Define the logarithmic distance between stock prices S and S as ln(s ) ln(s) = ln(s /S). a Dai (R , D ), Wang (F ), Lyuu, and Liu (2010). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 704
8 Tree for Merton s Jump-Diffusion Model (continued) Take the logarithm of Eq. (82) on p. 702: ( ) St M t ln = X t + Y t, (83) S 0 where X t ( r λ k σ 2 /2 ) t + σw t, (84) Y t n(t) i=0 ln (1 + k i ). (85) It decomposes the S 0 -logarithmic return of S t into the diffusion component X t and the jump component Y t. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 705
9 Tree for Merton s Jump-Diffusion Model (continued) Motivated by decomposition (83) on p. 705, the tree construction divides each period into a diffusion phase followed by a jump phase. In the diffusion phase, X t BOPM. is approximated by the Hence X t canmakeanupmoveto X t + σ Δt with probability p u or a down move to X t σ Δt with probability p d. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 706
10 Tree for Merton s Jump-Diffusion Model (continued) According to BOPM, p u = eμδt d u d, p d = 1 p u, except that μ = r λ k here. The diffusion component gives rise to diffusion nodes. They are spaced at 2σ Δt apart such as the white nodes A, B, C, D, E, F, and G on p c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 707
11 ( l 1) Δt l Δt ( l + 1) Δt p u p d q 1 q 0 q h = γ + δ 2σ Δt White nodes are diffusion nodes. Gray nodes are jump nodes. In the diffusion phase, the solid black lines denote the binomial structure of BOPM, whereas the dashed lines denote the trinomial structure. Here m is set to one here for simplicity. Only the doublecircled nodes will remain after the construction. Note that a and b are diffusion nodes because no jump occurs in the jump phase. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 708
12 Tree for Merton s Jump-Diffusion Model (concluded) In the jump phase, Y t+δt is approximated by moves from each diffusion node to 2m jump nodes that match the first 2m moments of the lognormal jump. The m jump nodes above the diffusion node are spaced at h apart. Thesameholdsforthe m jump nodes below the diffusion node. The gray nodes at time lδt on p. 708 are jump nodes. After some work, the size of the tree is O(n 2.5 ). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 709
13 Multivariate Contingent Claims They depend on two or more underlying assets. The basket call on m assets has the terminal payoff ( m ) max α i S i (τ) X, 0, i=1 where α i is the percentage of asset i. Basket options are essentially options on a portfolio of stocks or index options. Option on the best of two risky assets and cash has a terminal payoff of max(s 1 (τ),s 2 (τ),x). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 710
14 Multivariate Contingent Claims (concluded) a Name Payoff Exchange option max(s 1 (τ) S 2 (τ), 0) Better-off option max(s 1 (τ),...,s k (τ), 0) Worst-off option min(s 1 (τ),...,s k (τ), 0) Binary maximum option I{ max(s 1 (τ),...,s k (τ)) >X} Maximum option max(max(s 1 (τ),...,s k (τ)) X, 0) Minimum option max(min(s 1 (τ),...,s k (τ)) X, 0) Spread option max(s 1 (τ) S 2 (τ) X, 0) Basket average option max((s 1 (τ),...,s k (τ))/k X, 0) Multi-strike option max(s 1 (τ) X 1,...,S k (τ) X k, 0) Pyramid rainbow option max( S 1 (τ) X S k (τ) X k X, 0) Madonna option max( (S 1 (τ) X 1 ) 2 + +(S k (τ) X k ) 2 X, 0) a Lyuu and Teng (R ) (2011). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 711
15 Correlated Trinomial Model a Two risky assets S 1 and S 2 follow ds i S i = rdt+ σ i dw i in a risk-neutral economy, i =1, 2. Let M i e rδt, V i M 2 i (e σ2 i Δt 1). S i M i is the mean of S i at time Δt. S 2 i V i the variance of S i at time Δt. a Boyle, Evnine, and Gibbs (1989). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 712
16 Correlated Trinomial Model (continued) The value of S 1 S 2 at time Δt has a joint lognormal distribution with mean S 1 S 2 M 1 M 2 e ρσ 1σ 2 Δt,where ρ is the correlation between dw 1 and dw 2. Next match the 1st and 2nd moments of the approximating discrete distribution to those of the continuous counterpart. At time Δt from now, there are five distinct outcomes. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 713
17 Correlated Trinomial Model (continued) The five-point probability distribution of the asset prices is (as usual, we impose u i d i =1) Probability Asset 1 Asset 2 p 1 S 1 u 1 S 2 u 2 p 2 S 1 u 1 S 2 d 2 p 3 S 1 d 1 S 2 d 2 p 4 S 1 d 1 S 2 u 2 p 5 S 1 S 2 c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 714
18 Correlated Trinomial Model (continued) The probabilities must sum to one, and the means must be matched: 1 = p 1 + p 2 + p 3 + p 4 + p 5, S 1 M 1 = (p 1 + p 2 ) S 1 u 1 + p 5 S 1 +(p 3 + p 4 ) S 1 d 1, S 2 M 2 = (p 1 + p 4 ) S 2 u 2 + p 5 S 2 +(p 2 + p 3 ) S 2 d 2. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 715
19 Correlated Trinomial Model (concluded) Let R M 1 M 2 e ρσ 1σ 2 Δt. Match the variances and covariance: S 2 1 V 1 = (p 1 + p 2 )((S 1 u 1 ) 2 (S 1 M 1 ) 2 )+p 5 (S 2 1 (S 1M 1 ) 2 ) +(p 3 + p 4 )((S 1 d 1 ) 2 (S 1 M 1 ) 2 ), S 2 2 V 2 = (p 1 + p 4 )((S 2 u 2 ) 2 (S 2 M 2 ) 2 )+p 5 (S 2 2 (S 2M 2 ) 2 ) +(p 2 + p 3 )((S 2 d 2 ) 2 (S 2 M 2 ) 2 ), S 1 S 2 R = (p 1 u 1 u 2 + p 2 u 1 d 2 + p 3 d 1 d 2 + p 4 d 1 u 2 + p 5 ) S 1 S 2. The solutions are complex (see text). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 716
20 Correlated Trinomial Model Simplified a Let μ i r σ2 i /2andu i e λσ i Δt for i =1, 2. The following simpler scheme is good enough: p 1 = p 2 = p 3 = p 4 = p 5 = 1 1 λ 2. a Madan, Milne, and Shefrin (1989). [ ( 1 1 Δt μ 4 λ μ ) 2 + ρ ] λ σ 1 σ 2 λ 2, [ ( 1 1 Δt μ 4 λ μ ) 2 ρ ] λ σ 1 σ 2 λ 2, [ ( 1 1 Δt 4 λ 2 + μ 1 μ ) 2 + ρ ] λ σ 1 σ 2 λ 2, [ ( 1 1 Δt 4 λ 2 + μ 1 + μ ) 2 ρ ] λ σ 1 σ 2 λ 2, c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 717
21 Correlated Trinomial Model Simplified (continued) All of the probabilities lie between 0 and 1 if and only if 1+λ Δt μ 1 σ 1 + μ 2 σ 2 ρ 1 λ Δt μ 1 μ 2 σ 1 σ, (86) 2 1 λ (87) We call a multivariate tree (correlation-) optimal if it guarantees valid probabilities as long as 1+O( Δt) <ρ<1 O( Δt), such as the above one. a a Kao (R ) (2011) and Kao (R ), Lyuu, and Wen (D ) (2014). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 718
22 Correlated Trinomial Model Simplified (concluded) But this model cannot price 2-asset 2-barrier options accurately. a Few multivariate trees are both optimal and able to handle multiple barriers. b An alternative is to use orthogonalization. c a See Chang (B , R ), Hsu (R , D ), and Lyuu (2006) and Kao (R ), Lyuu and Wen (D ) (2014) for solutions. b See Kao (R ), Lyuu, and Wen (D ) (2014) for one. c Hull and White (1990) and Dai (R , D ), Lyuu, and Wang (F ) (2012). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 719
23 Extrapolation It is a method to speed up numerical convergence. Say f(n) converges to an unknown limit f at rate of 1/n: f(n) =f + c ( ) 1 n + o. (88) n Assume c is an unknown constant independent of n. Convergence is basically monotonic and smooth. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 720
24 Extrapolation (concluded) From two approximations f(n 1 ) and f(n 2 ) and ignoring the smaller terms, f(n 1 ) = f + c, n 1 f(n 2 ) = f + c. n 2 A better approximation to the desired f is f = n 1f(n 1 ) n 2 f(n 2 ) n 1 n 2. (89) This estimate should converge faster than 1/n. The Richardson extrapolation uses n 2 =2n 1. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 721
25 Improving BOPM with Extrapolation Consider standard European options. Denote the option value under BOPM using n time periods by f(n). It is known that BOPM convergences at the rate of 1/n, consistent with Eq. (88) on p But the plots on p. 282 (redrawn on next page) demonstrate that convergence to the true option value oscillates with n. Extrapolation is inapplicable at this stage. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 722
26 Call value n Call value n c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 723
27 Improving BOPM with Extrapolation (concluded) Take the at-the-money option in the left plot on p The sequence with odd n turns out to be monotonic and smooth (see the left plot on p. 725). a Apply extrapolation (89) on p. 721 with n 2 = n 1 +2, where n 1 is odd. Result is shown in the right plot on p The convergence rate is amazing. See Exercise of the text (p. 111) for ideas in the general case. a This can be proved; see Chang and Palmer (2007). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 724
28 Call value n Call value n c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 725
29 Numerical Methods c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 726
30 All science is dominated by the idea of approximation. Bertrand Russell c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 727
31 Finite-Difference Methods Place a grid of points on the space over which the desired function takes value. Then approximate the function value at each of these points (p. 729). Solve the equation numerically by introducing difference equations in place of derivatives. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 728
32 c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 729
33 Example: Poisson s Equation It is 2 θ/ x θ/ y 2 = ρ(x, y). Replace second derivatives with finite differences through central difference. Introduce evenly spaced grid points with distance of Δx along the x axis and Δy along the y axis. The finite difference form is ρ(x i,y j )= θ(x i+1,y j ) 2θ(x i,y j )+θ(x i 1,y j ) (Δx) 2 + θ(x i,y j+1 ) 2θ(x i,y j )+θ(x i,y j 1 ) (Δy) 2. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 730
34 Example: Poisson s Equation (concluded) In the above, Δx x i x i 1 and Δy y j y j 1 for i, j =1, 2,... When the grid points are evenly spaced in both axes so that Δx =Δy = h, the difference equation becomes h 2 ρ(x i,y j )=θ(x i+1,y j )+θ(x i 1,y j ) +θ(x i,y j+1 )+θ(x i,y j 1 ) 4θ(x i,y j ). Given boundary values, we can solve for the x i sandthe y j s within the square [ ±L, ±L ]. From now on, θ i,j will denote the finite-difference approximation to the exact θ(x i,y j ). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 731
35 Explicit Methods Consider the diffusion equation D( 2 θ/ x 2 ) ( θ/ t)=0,d>0. Use evenly spaced grid points (x i,t j ) with distances Δx and Δt, where Δx x i+1 x i and Δt t j+1 t j. Employ central difference for the second derivative and forward difference for the time derivative to obtain θ(x, t) t = θ(x, t j+1) θ(x, t j ) +, (90) t=tj Δt 2 θ(x, t) x 2 = θ(x i+1,t) 2θ(x i,t)+θ(x i 1,t) x=xi (Δx) 2 +. (91) c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 732
36 Explicit Methods (continued) Next, assemble Eqs. (90) and (91) into a single equation at (x i,t j ). But we need to decide how to evaluate x in the first equation and t in the second. Since central difference around x i is used in Eq. (91), we might as well use x i for x in Eq. (90). Two choices are possible for t in Eq. (91). The first choice uses t = t j finite-difference equation, to yield the following θ i,j+1 θ i,j Δt = D θ i+1,j 2θ i,j + θ i 1,j (Δx) 2. (92) c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 733
37 Explicit Methods (continued) The stencil of grid points involves four values, θ i,j+1, θ i,j, θ i+1,j,andθ i 1,j. Rearrange Eq. (92) on p. 733 as θ i,j+1 = DΔt (Δx) 2 θ i+1,j + ( 1 2DΔt (Δx) 2 ) θ i,j + DΔt (Δx) 2 θ i 1,j. We can calculate θ i,j+1 from θ i,j,θ i+1,j,θ i 1,j,atthe previous time t j (see exhibit (a) on next page). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 734
38 Stencils c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 735
39 Explicit Methods (concluded) Starting from the initial conditions at t 0,thatis, θ i,0 = θ(x i,t 0 ), i =1, 2,...,wecalculate θ i,1, i =1, 2,.... And then θ i,2, i =1, 2,.... And so on. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 736
40 Stability The explicit method is numerically unstable unless Δt (Δx) 2 /(2D). A numerical method is unstable if the solution is highly sensitive to changes in initial conditions. The stability condition may lead to high running times and memory requirements. For instance, halving Δx would imply quadrupling (Δt) 1, resulting in a running time 8 times as much. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 737
41 Recall that Explicit Method and Trinomial Tree θ i,j+1 = DΔt (Δx) 2 θ i+1,j + ( 1 2DΔt ) θ (Δx) 2 i,j + DΔt (Δx) θ i 1,j. 2 When the stability condition is satisfied, the three coefficients for θ i+1,j, θ i,j,and θ i 1,j all lie between zero and one and sum to one. They can be interpreted as probabilities. So the finite-difference equation becomes identical to backward induction on trinomial trees! The freedom in choosing Δx corresponds to similar freedom in the construction of trinomial trees. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 738
42 Implicit Methods Suppose we use t = t j+1 in Eq. (91) on p. 732 instead. The finite-difference equation becomes θ i,j+1 θ i,j Δt = D θ i+1,j+1 2θ i,j+1 + θ i 1,j+1 (Δx) 2. (93) The stencil involves θ i,j, θ i,j+1, θ i+1,j+1,andθ i 1,j+1. This method is implicit: The value of any one of the three quantities at t j+1 cannot be calculated unless the other two are known. See exhibit (b) on p c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 739
43 Implicit Methods (continued) Equation (93) can be rearranged as θ i 1,j+1 (2 + γ) θ i,j+1 + θ i+1,j+1 = γθ i,j, where γ (Δx) 2 /(DΔt). This equation is unconditionally stable. Suppose the boundary conditions are given at x = x 0 and x = x N+1. After θ i,j has been calculated for i =1, 2,...,N,the values of θ i,j+1 at time t j+1 can be computed as the solution to the following tridiagonal linear system, c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 740
44 Implicit Methods (continued) a a a a a θ 1,j+1 θ 2,j+1 θ 3,j θ N,j+1 = γθ 1,j θ 0,j+1 γθ 2,j γθ 3,j γθ N 1,j γθ N,j θ N+1,j+1, where a 2 γ. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 741
45 Implicit Methods (concluded) Tridiagonal systems can be solved in O(N) time and O(N) space. Never invert a matrix to solve a tridiagonal system. The matrix above is nonsingular when γ 0. A square matrix is nonsingular if its inverse exists. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 742
46 Crank-Nicolson Method Take the average of explicit method (92) on p. 733 and implicit method (93) on p. 739: θ i,j+1 θ i,j = 1 2 Δt ( D θ i+1,j 2θ i,j + θ i 1,j (Δx) 2 + D θ i+1,j+1 2θ i,j+1 + θ i 1,j+1 (Δx) 2 ). After rearrangement, γθ i,j+1 θ i+1,j+1 2θ i,j+1 + θ i 1,j+1 2 = γθ i,j + θ i+1,j 2θ i,j + θ i 1,j 2. This is an unconditionally stable implicit method with excellent rates of convergence. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 743
47 Stencil x i+1 x i x i+1 t j t j+1 c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 744
48 Numerically Solving the Black-Scholes PDE (65) on p. 583 See text. Brennan and Schwartz (1978) analyze the stability of the implicit method. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 745
49 Monte Carlo Simulation a Monte Carlo simulation is a sampling scheme. In many important applications within finance and without, Monte Carlo is one of the few feasible tools. When the time evolution of a stochastic process is not easy to describe analytically, Monte Carlo may very well be the only strategy that succeeds consistently. a A top 10 algorithm according to Dongarra and Sullivan (2000). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 746
50 The Big Idea Assume X 1,X 2,...,X n have a joint distribution. θ E[ g(x 1,X 2,...,X n ) ] for some function g is desired. We generate ( x (i) 1,x(i) 2,...,x(i) n ), 1 i N independently with the same joint distribution as (X 1,X 2,...,X n ). Set Y i g ( x (i) 1,x(i) 2,...,x(i) n ). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 747
51 The Big Idea (concluded) Y 1,Y 2,...,Y N are independent and identically distributed random variables. Each Y i has the same distribution as Y g(x 1,X 2,...,X n ). Since the average of these N random variables, Y, satisfies E[ Y ]=θ, it can be used to estimate θ. The strong law of large numbers says that this procedure converges almost surely. The number of replications (or independent trials), N, is called the sample size. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 748
52 Accuracy The Monte Carlo estimate and true value may differ owing to two reasons: 1. Sampling variation. 2. The discreteness of the sample paths. a The first can be controlled by the number of replications. The second can be controlled by the number of observations along the sample path. a This may not be an issue if the financial derivative only requires discrete sampling along the time dimension, such as the discrete barrier option. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 749
53 Accuracy and Number of Replications The statistical error of the sample mean Y of the random variable Y grows as 1/ N. Because Var[ Y ]=Var[Y ]/N. In fact, this convergence rate is asymptotically optimal. a So the variance of the estimator Y can be reduced by a factor of 1/N by doing N times as much work. This is amazing because the same order of convergence holds independently of the dimension n. a The Berry-Esseen theorem. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 750
54 Accuracy and Number of Replications (concluded) In contrast, classic numerical integration schemes have an error bound of O(N c/n ) for some constant c>0. n is the dimension. The required number of evaluations thus grows exponentially in n to achieve a given level of accuracy. The curse of dimensionality. The Monte Carlo method is more efficient than alternative procedures for multivariate derivatives when n is large. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 751
55 Monte Carlo Option Pricing For the pricing of European options on a dividend-paying stock, we may proceed as follows. Assume ds = μdt+ σdw. S Stock prices S 1,S 2,S 3,... at times Δt, 2Δt, 3Δt,... can be generated via S i+1 = S i e (μ σ2 /2) Δt+σ Δtξ, ξ N(0, 1). (94) c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 752
56 Monte Carlo Option Pricing (continued) If we discretize ds/s = μdt+ σdw directly, we will obtain S i+1 = S i + S i μ Δt + S i σ Δtξ. But this is locally normally distributed, not lognormally, hence biased. a In practice, this is not expected to be a major problem as long as Δt is sufficiently small. a Contributed by Mr. Tai, Hui-Chin (R ) on April 22, c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 753
57 Monte Carlo Option Pricing (continued) Non-dividend-paying stock prices in a risk-neutral economy can be generated by setting μ = r and Δt = T. 1: C := 0; {Accumulated terminal option value.} 2: for i =1, 2, 3,...,N do 3: P := S e (r σ2 /2) T +σ T ξ, ξ N(0, 1); 4: C := C +max(p X, 0); 5: end for 6: return Ce rt /N ; c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 754
58 Monte Carlo Option Pricing (concluded) Pricing Asian options is also easy. 1: C := 0; 2: for i =1, 2, 3,...,N do 3: P := S; M := S; 4: for j =1, 2, 3,...,n do 5: P := P e (r σ2 /2)(T/n)+σ T/n ξ ; 6: M := M + P ; 7: end for 8: C := C +max(m/(n +1) X, 0); 9: end for 10: return Ce rt /N ; c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 755
59 How about American Options? Standard Monte Carlo simulation is inappropriate for American options because of early exercise (why?). It is difficult to determine the early-exercise point based on one single path. But Monte Carlo simulation can be modified to price American options with small biases (pp. 807ff). a a Longstaff and Schwartz (2001). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 756
60 Delta and Common Random Numbers In estimating delta, it is natural to start with the finite-difference estimate rτ E[ P (S + ɛ)] E[ P (S ɛ)] e. 2ɛ P (x) is the terminal payoff of the derivative security when the underlying asset s initial price equals x. Use simulation to estimate E[ P (S + ɛ)] first. Use another simulation to estimate E[ P (S ɛ)]. Finally, apply the formula to approximate the delta. This is also called the bump-and-revalue method. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 757
61 Delta and Common Random Numbers (concluded) This method is not recommended because of its high variance. A much better approach is to use common random numbers to lower the variance: [ ] P (S + ɛ) P (S ɛ) e rτ E. 2ɛ Here, the same random numbers are used for P (S + ɛ) and P (S ɛ). This holds for gamma and cross gammas (for multivariate derivatives). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 758
62 Problems with the Bump-and-Revalue Method Consider the binary option with payoff 1, if S(T ) >X, 0, otherwise. Then P (S + ɛ) P (S ɛ) = 1, if P (S + ɛ) >X and P (S ɛ)]<x, 0, otherwise. So the finite-difference estimate per run for the (undiscounted) delta is 0 or O(1/ɛ). This means high variance. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 759
63 Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S)+P(S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gamma 2 P (S 1,S 2,...)/( S 1 S 2 ) is: [ P e rτ (S1 + ɛ 1,S 2 + ɛ 2 ) P (S 1 ɛ 1,S 2 + ɛ 2 ) E 4ɛ 1 ɛ 2 P (S 1 + ɛ 1,S 2 ɛ 2 )+P (S 1 ɛ 1,S 2 ɛ 2 ) ]. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 760
64 Gamma (continued) Choosing an ɛ of the right magnitude can be challenging. If ɛ is too large, inaccurate Greeks result. If ɛ is too small, unstable Greeks result. This phenomenon is sometimes called the curse of differentiation. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 761
65 Gamma (continued) In general, suppose i θ i e rτ E[ P (S)]=e rτ E [ i P (S) θ i ] holds for all i>0, where θ is a parameter of interest. Then formulas for the Greeks become integrals. As a result, we avoid ɛ, finite differences, and resimulation. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 762
66 Gamma (continued) This is indeed possible for a broad class of payoff functions. a Roughly speaking, any payoff function that is equal to a sum of products of differentiable functions and indicator functions with the right kind of support. For example, the payoff of a call is max(s(t ) X, 0) = (S(T ) X)I { S(T ) X 0 }. The results are too technical to cover here (see next page). a Teng (R ) (2004) and Lyuu and Teng (R ) (2011). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 763
67 Gamma (continued) Suppose h(θ, x) Hwith pdf f(x) forx and g j (θ, x) G for j B, a finite set of natural numbers. Then = h(θ, x) θ R j B R h θ (θ, x) j B + l B h(θ, x)j l (θ, x) 1 {gj (θ,x)>0} (x) f(x) dx 1 {gj (θ,x)>0} (x) f(x) dx 1 {gj (θ, x)>0} (x) f(x), j B\l x=χ l (θ) where ( ) gl (θ, x) gl (θ, x)/ θ J l (θ, x) =sign x k g l (θ, x)/ x for l B. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 764
68 Gamma (concluded) Similar results have been derived for Levy processes. a Formulas are also recently obtained for credit derivatives. b In queueing networks, this is called infinitesimal perturbation analysis (IPA). c a Lyuu, Teng (R ), and Wang (2013). b Lyuu, Teng (R ), and Tzeng (2014). c Cao (1985); Ho and Cao (1985). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 765
69 Biases in Pricing Continuously Monitored Options with Monte Carlo We are asked to price a continuously monitored up-and-out call with barrier H. The Monte Carlo method samples the stock price at n discrete time points t 1,t 2,...,t n. A sample path is produced. S(t 0 ),S(t 1 ),...,S(t n ) Here, t 0 = 0 is the current time, and t n = T is the expiration time of the option. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 766
70 Biases in Pricing Continuously Monitored Options with Monte Carlo (continued) If all of the sampled prices are below the barrier, this sample path pays max(s(t n ) X, 0). Repeating these steps and averaging the payoffs yield a Monte Carlo estimate. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 767
71 1: C := 0; 2: for i =1, 2, 3,...,N do 3: P := S; hit := 0; 4: for j =1, 2, 3,...,n do 5: P := P e (r σ2 /2) (T/n)+σ (T/n) ξ ; 6: if P H then 7: hit := 1; 8: break; 9: end if 10: end for 11: if hit = 0 then 12: C := C +max(p X, 0); 13: end if 14: end for 15: return Ce rt /N ; c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 768
72 Biases in Pricing Continuously Monitored Options with Monte Carlo (continued) This estimate is biased. a Suppose none of the sampled prices on a sample path equals or exceeds the barrier H. It remains possible for the continuous sample path that passes through them to hit the barrier between sampled time points (see plot on next page). a Shevchenko (2003). c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 769
73 H c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 770
74 Biases in Pricing Continuously Monitored Options with Monte Carlo (concluded) The bias can certainly be lowered by increasing the number of observations along the sample path. However, even daily sampling may not suffice. The computational cost also rises as a result. c 2015 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 771
Trinomial Tree. Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a
Trinomial Tree Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a The three stock prices at time t are S, Su, and Sd, where ud = 1. Impose the matching of mean and
More informationTrinomial Tree. Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a
Trinomial Tree Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a The three stock prices at time t are S, Su, and Sd, where ud = 1. Impose the matching of mean and
More informationGamma. The finite-difference formula for gamma is
Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S) + P (S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gammas
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationBrownian Bridge Approach to Pricing Barrier Options
Brownian Bridge Approach to Pricing Barrier Options We desire an unbiased estimate which can be calculated efficiently. The above-mentioned payoff should be multiplied by the probability p that a continuous
More informationEquilibrium Term Structure Models. c 2008 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 854
Equilibrium Term Structure Models c 2008 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 854 8. What s your problem? Any moron can understand bond pricing models. Top Ten Lies Finance Professors Tell
More informationMAFS Computational Methods for Pricing Structured Products
MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 )
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More informationMonte Carlo Methods in Structuring and Derivatives Pricing
Monte Carlo Methods in Structuring and Derivatives Pricing Prof. Manuela Pedio (guest) 20263 Advanced Tools for Risk Management and Pricing Spring 2017 Outline and objectives The basic Monte Carlo algorithm
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets (Hull chapter: 12, 13, 14) Liuren Wu ( c ) The Black-Scholes Model colorhmoptions Markets 1 / 17 The Black-Scholes-Merton (BSM) model Black and Scholes
More informationFIN FINANCIAL INSTRUMENTS SPRING 2008
FIN-40008 FINANCIAL INSTRUMENTS SPRING 2008 The Greeks Introduction We have studied how to price an option using the Black-Scholes formula. Now we wish to consider how the option price changes, either
More informationFixed-Income Options
Fixed-Income Options Consider a two-year 99 European call on the three-year, 5% Treasury. Assume the Treasury pays annual interest. From p. 852 the three-year Treasury s price minus the $5 interest could
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationOption Pricing Models. c 2013 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 205
Option Pricing Models c 2013 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 205 If the world of sense does not fit mathematics, so much the worse for the world of sense. Bertrand Russell (1872 1970)
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets Liuren Wu ( c ) The Black-Merton-Scholes Model colorhmoptions Markets 1 / 18 The Black-Merton-Scholes-Merton (BMS) model Black and Scholes (1973) and Merton
More informationEstimating the Greeks
IEOR E4703: Monte-Carlo Simulation Columbia University Estimating the Greeks c 207 by Martin Haugh In these lecture notes we discuss the use of Monte-Carlo simulation for the estimation of sensitivities
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationPricing Convertible Bonds under the First-Passage Credit Risk Model
Pricing Convertible Bonds under the First-Passage Credit Risk Model Prof. Tian-Shyr Dai Department of Information Management and Finance National Chiao Tung University Joint work with Prof. Chuan-Ju Wang
More informationChapter 5 Finite Difference Methods. Math6911 W07, HM Zhu
Chapter 5 Finite Difference Methods Math69 W07, HM Zhu References. Chapters 5 and 9, Brandimarte. Section 7.8, Hull 3. Chapter 7, Numerical analysis, Burden and Faires Outline Finite difference (FD) approximation
More informationProblems; the Smile. Options written on the same underlying asset usually do not produce the same implied volatility.
Problems; the Smile Options written on the same underlying asset usually do not produce the same implied volatility. A typical pattern is a smile in relation to the strike price. The implied volatility
More information2 f. f t S 2. Delta measures the sensitivityof the portfolio value to changes in the price of the underlying
Sensitivity analysis Simulating the Greeks Meet the Greeks he value of a derivative on a single underlying asset depends upon the current asset price S and its volatility Σ, the risk-free interest rate
More informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationToward the Black-Scholes Formula
Toward the Black-Scholes Formula The binomial model seems to suffer from two unrealistic assumptions. The stock price takes on only two values in a period. Trading occurs at discrete points in time. As
More informationPricing with a Smile. Bruno Dupire. Bloomberg
CP-Bruno Dupire.qxd 10/08/04 6:38 PM Page 1 11 Pricing with a Smile Bruno Dupire Bloomberg The Black Scholes model (see Black and Scholes, 1973) gives options prices as a function of volatility. If an
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationFINITE DIFFERENCE METHODS
FINITE DIFFERENCE METHODS School of Mathematics 2013 OUTLINE Review 1 REVIEW Last time Today s Lecture OUTLINE Review 1 REVIEW Last time Today s Lecture 2 DISCRETISING THE PROBLEM Finite-difference approximations
More informationAn Adjusted Trinomial Lattice for Pricing Arithmetic Average Based Asian Option
American Journal of Applied Mathematics 2018; 6(2): 28-33 http://www.sciencepublishinggroup.com/j/ajam doi: 10.11648/j.ajam.20180602.11 ISSN: 2330-0043 (Print); ISSN: 2330-006X (Online) An Adjusted Trinomial
More informationANALYSIS OF THE BINOMIAL METHOD
ANALYSIS OF THE BINOMIAL METHOD School of Mathematics 2013 OUTLINE 1 CONVERGENCE AND ERRORS OUTLINE 1 CONVERGENCE AND ERRORS 2 EXOTIC OPTIONS American Options Computational Effort OUTLINE 1 CONVERGENCE
More informationThe Binomial Model. The analytical framework can be nicely illustrated with the binomial model.
The Binomial Model The analytical framework can be nicely illustrated with the binomial model. Suppose the bond price P can move with probability q to P u and probability 1 q to P d, where u > d: P 1 q
More informationComputational Finance Binomial Trees Analysis
Computational Finance Binomial Trees Analysis School of Mathematics 2018 Review - Binomial Trees Developed a multistep binomial lattice which will approximate the value of a European option Extended the
More informationNumerical Evaluation of Multivariate Contingent Claims
Numerical Evaluation of Multivariate Contingent Claims Phelim P. Boyle University of California, Berkeley and University of Waterloo Jeremy Evnine Wells Fargo Investment Advisers Stephen Gibbs University
More informationLecture Quantitative Finance Spring Term 2015
and Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 06: March 26, 2015 1 / 47 Remember and Previous chapters: introduction to the theory of options put-call parity fundamentals
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 13 Lecture 13 November 15, 217 Derivation of the Black-Scholes-Merton
More informationHandbook of Financial Risk Management
Handbook of Financial Risk Management Simulations and Case Studies N.H. Chan H.Y. Wong The Chinese University of Hong Kong WILEY Contents Preface xi 1 An Introduction to Excel VBA 1 1.1 How to Start Excel
More informationMSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013
MSc Financial Engineering 2012-13 CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL To be handed in by monday January 28, 2013 Department EMS, Birkbeck Introduction The assignment consists of Reading
More informationAN IMPROVED BINOMIAL METHOD FOR PRICING ASIAN OPTIONS
Commun. Korean Math. Soc. 28 (2013), No. 2, pp. 397 406 http://dx.doi.org/10.4134/ckms.2013.28.2.397 AN IMPROVED BINOMIAL METHOD FOR PRICING ASIAN OPTIONS Kyoung-Sook Moon and Hongjoong Kim Abstract. We
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 13. The Black-Scholes PDE Steve Yang Stevens Institute of Technology 04/25/2013 Outline 1 The Black-Scholes PDE 2 PDEs in Asset Pricing 3 Exotic
More informationThe Black-Scholes Model
IEOR E4706: Foundations of Financial Engineering c 2016 by Martin Haugh The Black-Scholes Model In these notes we will use Itô s Lemma and a replicating argument to derive the famous Black-Scholes formula
More informationUtility Indifference Pricing and Dynamic Programming Algorithm
Chapter 8 Utility Indifference ricing and Dynamic rogramming Algorithm In the Black-Scholes framework, we can perfectly replicate an option s payoff. However, it may not be true beyond the Black-Scholes
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationThe Black-Scholes PDE from Scratch
The Black-Scholes PDE from Scratch chris bemis November 27, 2006 0-0 Goal: Derive the Black-Scholes PDE To do this, we will need to: Come up with some dynamics for the stock returns Discuss Brownian motion
More informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
More informationRisk Neutral Valuation
copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential
More informationComputational Finance Improving Monte Carlo
Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal
More informationMATH 476/567 ACTUARIAL RISK THEORY FALL 2016 PROFESSOR WANG
MATH 476/567 ACTUARIAL RISK THEORY FALL 206 PROFESSOR WANG Homework 5 (max. points = 00) Due at the beginning of class on Tuesday, November 8, 206 You are encouraged to work on these problems in groups
More informationComputational Efficiency and Accuracy in the Valuation of Basket Options. Pengguo Wang 1
Computational Efficiency and Accuracy in the Valuation of Basket Options Pengguo Wang 1 Abstract The complexity involved in the pricing of American style basket options requires careful consideration of
More informationDRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics
Chapter 12 American Put Option Recall that the American option has strike K and maturity T and gives the holder the right to exercise at any time in [0, T ]. The American option is not straightforward
More informationAD in Monte Carlo for finance
AD in Monte Carlo for finance Mike Giles giles@comlab.ox.ac.uk Oxford University Computing Laboratory AD & Monte Carlo p. 1/30 Overview overview of computational finance stochastic o.d.e. s Monte Carlo
More informationBinomial model: numerical algorithm
Binomial model: numerical algorithm S / 0 C \ 0 S0 u / C \ 1,1 S0 d / S u 0 /, S u 3 0 / 3,3 C \ S0 u d /,1 S u 5 0 4 0 / C 5 5,5 max X S0 u,0 S u C \ 4 4,4 C \ 3 S u d / 0 3, C \ S u d 0 S u d 0 / C 4
More informationMATH6911: Numerical Methods in Finance. Final exam Time: 2:00pm - 5:00pm, April 11, Student Name (print): Student Signature: Student ID:
MATH6911 Page 1 of 16 Winter 2007 MATH6911: Numerical Methods in Finance Final exam Time: 2:00pm - 5:00pm, April 11, 2007 Student Name (print): Student Signature: Student ID: Question Full Mark Mark 1
More informationComputational Finance Finite Difference Methods
Explicit finite difference method Computational Finance Finite Difference Methods School of Mathematics 2018 Today s Lecture We now introduce the final numerical scheme which is related to the PDE solution.
More informationTEACHING NOTE 98-04: EXCHANGE OPTION PRICING
TEACHING NOTE 98-04: EXCHANGE OPTION PRICING Version date: June 3, 017 C:\CLASSES\TEACHING NOTES\TN98-04.WPD The exchange option, first developed by Margrabe (1978), has proven to be an extremely powerful
More information1 The continuous time limit
Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1
More informationThe Pennsylvania State University. The Graduate School. Department of Industrial Engineering AMERICAN-ASIAN OPTION PRICING BASED ON MONTE CARLO
The Pennsylvania State University The Graduate School Department of Industrial Engineering AMERICAN-ASIAN OPTION PRICING BASED ON MONTE CARLO SIMULATION METHOD A Thesis in Industrial Engineering and Operations
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationAdvanced Numerical Methods
Advanced Numerical Methods Solution to Homework One Course instructor: Prof. Y.K. Kwok. When the asset pays continuous dividend yield at the rate q the expected rate of return of the asset is r q under
More informationInterest-Sensitive Financial Instruments
Interest-Sensitive Financial Instruments Valuing fixed cash flows Two basic rules: - Value additivity: Find the portfolio of zero-coupon bonds which replicates the cash flows of the security, the price
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationFINANCIAL OPTION ANALYSIS HANDOUTS
FINANCIAL OPTION ANALYSIS HANDOUTS 1 2 FAIR PRICING There is a market for an object called S. The prevailing price today is S 0 = 100. At this price the object S can be bought or sold by anyone for any
More informationChapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets
Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,
More informationDefinition Pricing Risk management Second generation barrier options. Barrier Options. Arfima Financial Solutions
Arfima Financial Solutions Contents Definition 1 Definition 2 3 4 Contenido Definition 1 Definition 2 3 4 Definition Definition: A barrier option is an option on the underlying asset that is activated
More informationLecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.
Lecture 7 Overture to continuous models Before rigorously deriving the acclaimed Black-Scholes pricing formula for the value of a European option, we developed a substantial body of material, in continuous
More informationFast Convergence of Regress-later Series Estimators
Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser
More information( ) since this is the benefit of buying the asset at the strike price rather
Review of some financial models for MAT 483 Parity and Other Option Relationships The basic parity relationship for European options with the same strike price and the same time to expiration is: C( KT
More informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
More informationHEDGING RAINBOW OPTIONS IN DISCRETE TIME
Journal of the Chinese Statistical Association Vol. 50, (2012) 1 20 HEDGING RAINBOW OPTIONS IN DISCRETE TIME Shih-Feng Huang and Jia-Fang Yu Department of Applied Mathematics, National University of Kaohsiung
More informationMathematics of Finance Final Preparation December 19. To be thoroughly prepared for the final exam, you should
Mathematics of Finance Final Preparation December 19 To be thoroughly prepared for the final exam, you should 1. know how to do the homework problems. 2. be able to provide (correct and complete!) definitions
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationComputational Finance
Path Dependent Options Computational Finance School of Mathematics 2018 The Random Walk One of the main assumption of the Black-Scholes framework is that the underlying stock price follows a random walk
More informationNumerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationStochastic Processes and Brownian Motion
A stochastic process Stochastic Processes X = { X(t) } Stochastic Processes and Brownian Motion is a time series of random variables. X(t) (or X t ) is a random variable for each time t and is usually
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose
More informationPricing Dynamic Guaranteed Funds Under a Double Exponential. Jump Diffusion Process. Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay
Pricing Dynamic Guaranteed Funds Under a Double Exponential Jump Diffusion Process Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay ABSTRACT This paper complements the extant literature to evaluate the
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationBinomial Model for Forward and Futures Options
Binomial Model for Forward and Futures Options Futures price behaves like a stock paying a continuous dividend yield of r. The futures price at time 0 is (p. 437) F = Se rt. From Lemma 10 (p. 275), the
More informationComputational Finance. Computational Finance p. 1
Computational Finance Computational Finance p. 1 Outline Binomial model: option pricing and optimal investment Monte Carlo techniques for pricing of options pricing of non-standard options improving accuracy
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationNUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE
Trends in Mathematics - New Series Information Center for Mathematical Sciences Volume 13, Number 1, 011, pages 1 5 NUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE YONGHOON
More informationApproximation Methods in Derivatives Pricing
Approximation Methods in Derivatives Pricing Minqiang Li Bloomberg LP September 24, 2013 1 / 27 Outline of the talk A brief overview of approximation methods Timer option price approximation Perpetual
More informationDepartment of Mathematics. Mathematics of Financial Derivatives
Department of Mathematics MA408 Mathematics of Financial Derivatives Thursday 15th January, 2009 2pm 4pm Duration: 2 hours Attempt THREE questions MA408 Page 1 of 5 1. (a) Suppose 0 < E 1 < E 3 and E 2
More informationPricing Barrier Options under Local Volatility
Abstract Pricing Barrier Options under Local Volatility Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 16 November 2002 We study pricing under the local volatility. Our research is mainly
More informationA new PDE approach for pricing arithmetic average Asian options
A new PDE approach for pricing arithmetic average Asian options Jan Večeř Department of Mathematical Sciences, Carnegie Mellon University, Pittsburgh, PA 15213. Email: vecer@andrew.cmu.edu. May 15, 21
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationMATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, Student Name (print):
MATH4143 Page 1 of 17 Winter 2007 MATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, 2007 Student Name (print): Student Signature: Student ID: Question
More informationHedging Credit Derivatives in Intensity Based Models
Hedging Credit Derivatives in Intensity Based Models PETER CARR Head of Quantitative Financial Research, Bloomberg LP, New York Director of the Masters Program in Math Finance, Courant Institute, NYU Stanford
More informationMarket risk measurement in practice
Lecture notes on risk management, public policy, and the financial system Allan M. Malz Columbia University 2018 Allan M. Malz Last updated: October 23, 2018 2/32 Outline Nonlinearity in market risk Market
More informationDerivative Securities Fall 2012 Final Exam Guidance Extended version includes full semester
Derivative Securities Fall 2012 Final Exam Guidance Extended version includes full semester Our exam is Wednesday, December 19, at the normal class place and time. You may bring two sheets of notes (8.5
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationMASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.
MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS May/June 2006 Time allowed: 2 HOURS. Examiner: Dr N.P. Byott This is a CLOSED
More informationCalibration Lecture 4: LSV and Model Uncertainty
Calibration Lecture 4: LSV and Model Uncertainty March 2017 Recap: Heston model Recall the Heston stochastic volatility model ds t = rs t dt + Y t S t dw 1 t, dy t = κ(θ Y t ) dt + ξ Y t dw 2 t, where
More information