"Pricing Exotic Options using Strong Convergence Properties
|
|
- Claude Bryant
- 5 years ago
- Views:
Transcription
1 Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk Prof. Mike Giles Mike.Giles@comlab.ox.ac.uk
2 Outline Introduction (purpose of the talk) Strong and weak order of convergence [3] Orthogonal transformation [],[4] Pricing exotic options using ML-MC [1] Conclusions and comments [1] Michael Giles: Multi-level Monte Carlo path simulation (ML-MC). Technical Report No. NA06/03, Oxford University Computing Laboratory, Parks Road, Oxford, U.K., (006). [] A. B. Cruzeiro, P. Malliavin and A. Thalmaier: Geometrization of Monte-Carlo numerical analysis of an elliptic operator: strong approximation. C. R. Acad. Sci. Paris, Ser. I, 338, , (004). [3] K. Schmitz-Abe and W. T. Shaw: Measure Order of Convergence without an Exact Solution, Euler vs Milstein Scheme. International Journal of Pure and Applied Mathematics, Vol. 4, No 3, , (005). [4] K. Schmitz-Abe and M. Giles: Pricing Exotic Options using Strong Convergence Properties, ECMI 006, conference proceedings, to be published by Springer Verlag (September 007). "Pricing Exotic Options using Strong Schemes, Oxford, UK 01-0
3 Purpose of this talk Speed-up the computation of financial instruments TIME + accuracy(time) = = $ 0-0
4 Standard Euler & Milstein approximation Stochastic variance model: 1, (, ) (, ), t ds = µ S, ν dt + σ S, ν dw t dν = ϖ Sν dt+ ξ Sν dw Euler approximation: S µ σ 0 = t + W + W ν ϖ ρξ 1, t ρξ, t Milstein approximation: S ( ) t+ t = EULER + ϕ W 1, t, W, t + A1, A L ( 1, ν ) t Where the Lie brackets is defined by: σ ρξ ν A, A = 1 ( A A A A 1 1 ) = ξ ρσ S and the Lévy area is defined by: t+ t t+ t k t+ t k E dw 1, t, dw, t = ρdt E dw, dw ρ L = dw dw dw dw 1, t, t = 0 = 1 ρ ( 1,) 1,,, 1, t u k u k t t t t In practice: What is the main difference between both schemes? 03-0
5 Standard Euler & Milstein approximation Stochastic variance model: 1, (, ) (, ), t ds = µ S, ν dt + σ S, ν dw t dν = ϖ Sν dt+ ξ Sν dw E dw 1, t, dw, t = ρdt E dw, dw 1, t, t = 0 Simulation of the asset price S for a Monte Carlo path ( t=1/0) Is Euler scheme a good approximation? 04-0
6 Standard Euler & Milstein approximation Strong convergence Weak convergence 1 (, ) = E ST ST (, t) = C t γ E g ST E g ST t C t β 05-0
7 Standard Euler & Milstein approximation Strong convergence Weak convergence 1 (, ) = E ST ST (, t) = C t γ E g ST E g ST t C t β Euler: γ=0.51 Milstein: γ=1.0 Euler: β =1.05 Milstein: β =1.06 Can we take advantage of the strong convergence properties of Milstein scheme to price options? 05-0
8 Convergence without an Exact Solution Theorem 1: If a time discrete approximation Ŝ(T) converges strongly with order γ>0 at time T as t 0 to the exact solution S(T), i.e. exist a positive constant "C", that does not depend on γ such that: E S( T) S ( T, t) = C t γ Then using the same Wiener process, there exists a positive constant "C1", that is not dependent of γ, such that: t γ E S(T, t) S T, =C1 ( t) Surface error using the exact solution. Surface error using Theorem 1. [3] K. Schmitz-Abe and W. T. Shaw: Measure Order of Convergence without an Exact Solution, Euler vs Milstein Scheme. International Journal of Pure and Applied Mathematics, Vol. 4, No 3, , (005). 07-0
9 Problem with Milstein scheme Stochastic variance model: Milstein approximation: S ( ) t+ t = EULER + ϕ W 1, t, W, t + A1, A L ( 1, ν ) t Where the Lie brackets is defined by: σ ρξ ν A, A = 1 ξ ρσ and the Lévy area is defined by: S t+ t t+ t k t+ t k L = dw dw dw dw ( 1,) 1,,, 1, t 1, (, ) (, ), t ds = µ S, ν dt + σ S, ν dw t dν = ϖ Sν dt+ ξ Sν dw u k u k t t t t E dw 1, t, dw, t = ρdt E dw, dw ρ 1, t, t = 0 = 1 ρ Can we avoid the simulation of the Lévy area? 08-0
10 Orthogonal Milstein approximation If we replace the Wiener process dw by an orthogonal transform, the probability distribution does not change and we obtain the set of all orthogonal transforms: The Lie brackets for the new orthogonal process: dw 1, t cos θ sin θ dw1, t = dw, t sin ( θ) cos ( θ) dw, t ρξ σ S σ θ S ρσξ θ A, A = ν 1 ξ θ θ ρσ ξ ρσξ S ν S If we do not want to deal with the Lévy area, we need the Lie brackets to be identically zero, i.e.: A, A = 0 1 So then the rotation θ has to satisfy: dw θ 1 ξ σ ρ ξ θ 1 σ ξ ρ σ = and S + ρ = + ν ξ S ν σ ρ ξ S σ ν [4] K. Schmitz-Abe and M. Giles: Pricing Exotic Options using Strong Convergence Properties, ECMI 006, conference proceedings, to be published by Springer Verlag (September 007). 09-0
11 Example 1: D -- θ scheme Stochastic volatility model: ds = S µ dt + Sσ dw 1, t dσ = κ ( ϖ σ) dt + βσ λ dw, t where: µ, κ, ϖ, β, λ are constants E dw, 1, t dw, t = ρ dt ρ = 1 ρ If we do not want to deal with the Lévy area, θ has to satisfy: If we assume that θ is only a function of S and σ, we can solve it if λ= and the solution is: Remark: θ λ θ βσ θ ρ = and = S ρs σ ρσ ρ () t = log ( t ) log( S() t ) σ () ρ If λ = A, A = 0 A, A L = ( 1,) β ρ t+ t t [4] K. Schmitz-Abe and M. Giles: Pricing Exotic Options using Strong Convergence Properties, ECMI 006, conference proceedings, to be published by Springer Verlag (September 007). 10-0
12 Example : 3D -- θ scheme Stochastic variance model: Orthogonal conditions: α 1.5 = and = S S ν θ βν θ ρ ρ ρν Using the chain rule, our SVM (1) becomes a 3-dimensional Itô process: ds S µ S 0 ν dν κ ( ϖ ν ) dt dw βν α dw = + + θ θ dθ θ θ S S ν βν α µ + κ ϖ ν S ν S ν The Lie brackets using independent Wiener processes: A, A 1 = ds = Sµ dt + S νdw 1, t dν = κ ϖ ν dt + βν α dw, t 0 0 β S ν (1) where: µ, κ, ϖ, β, α are constants ( α 1.5) E dw, 1, t dw, t = ρ dt 0 1, t, t Note that: If α = 1.5 A, A = 0 1 () t ρ β = log ( v() t ) log( S() t ) How accuracy have to be θ to obtain 1 strong order convergence in S and v? θ ρ σ = ρ = 1 ρ ρ ( ν ) max,0 11-0
13 Strong convergence test ds = S ( µ dt + σdw 1, t ) Case 1, Case 1: Quadratic Volatility Model dσ = κ ϖ σ dt+ βσ dw, Case : 3/ Model (A. Lewis) 3/ dν = κ ϖ ν dt + βν dw 1, t Case 3: GARCH Diffusion Model dν = κ ( ϖ ν ) dt + β νdw, t Case 4: Heston Model 1 dν = κ ϖ ν dt + β νdw 3, t t Case 3 Case 4 S(0) = 1, σ(0) = 0.1, ν(0) = 0.1 Scheme Case 1, Case 3 Case 4 Euler scheme Milstein scheme (L=0) Milstein scheme D θ scheme 0.98 N/A N/A 3D θ scheme (L=0) D θ scheme Convergence orders for all cases T = 1, ρ = 0.3, µ = 0.05, κ = 1.8, ϖ = 0.6 S(0) = 1, ν(0) = 0.1, ν(0) = ν(0) β = 1, β = 0.5, β =
14 Theorem1: D Orthogonal Milstein scheme (Exact solution) If we have a -D Itô stochastic differential equation (SDE) with a -D Wiener process: X1 a1 b1,1 a1, W1 d dt d (1) X = a + b,1 a, W where we suppose ai,bi,j are sufficiently smooth functions of X & Y in [t₀...t] and satisfy the following condition: Ψ Φ = () X1 X where Hi are the coefficients of the Levy area (Lie bracket) of (1) and: H1( b,1+ b,) H( b1,1b,1+ b1,b,) H b1,1 b1, H1 b1,1b,1 b1,b, Ψ= Φ= b b b b ( 1,1, 1,,1 ) then; if we apply an orthogonal transformation to (1) described by: dw 1, t cos( θt) sin ( θt) dw1, t = where: dw, t sin ( θt) cos( θt) dw ±, t θ = Ψ dx +ΦdX t ( X, X ) 1 ( + ) ( + ) ( bb 1,1, bb 1,,1 ) 1 the new orthogonal process has 1 strong order convergence using Milstein scheme neglecting the simulation of the Lévy Area. Conversely, for Hi 0, the strong order Milstein scheme of (1) has 0.5 strong order convergence. Proof: Ph.D. Thesis; Schmitz, Klaus 13-0
15 Theorem: θ scheme If we have a -D Itô stochastic differential equation (SDE) with a -D Wiener process: X1 a1 b1,1 a1, W1 d dt d X = b a +,1 a, W (1) then, if θ is described using a third SDE: where: X a b b d a dt d 1 1 1,1 1, W1 X b,1 b, = + W θ Φ a1 + Ψa bθ,1 b θ, b b = Φ b + Ψb Φ b + Ψb θ,1 θ,1 1,1,1 1,, (3) then, the 3-dimensional SDE (3) can have better strong convergence than (1) using Milstein scheme neglecting the simulation of the Lévy Area. The accuracy of θ and hence in Xi depends directly in the value of the Lie bracket RL of the system: R L 0 = 0 Ψ Φ ( bb 11, bb 1,1 ) X1 X Proof: Ph.D. Thesis; Schmitz, Klaus 14-0
16 Pricing Exotic Options using ML-MC MC The expectation of a payoff P with maturity T is calculated by: L F E PL = E P F 0 + E PL P L 1 L = 1 [ ] where L is the level in the algorithm that simulates the scheme or time approximation with different time steps t: TT - t = 0 for: M L M L Z For each level L, we have a repeated cycle where we calculate the option price using NL paths, i.e. simulate extra samples at each level as needed for new NL. L N = ε V t Vl L L L l = 1 t l For a given ε = exact solution approximation, the algorithm will stop when it converges: { M 1 P P 1 } < ( M ) max, 1 ε L 1 L Computation cost Standard Method M SE c N 1 + c dt 1 = O ( ε -3 ) -3 ( ε ) - ε ( ε) -3 ( ε ) - ( ε ) Euler : O O log Milstein : O O ML-MC [1] Michael Giles: Multi-level Monte Carlo path simulation. Technical Report No. NA06/03, Oxford University Computing Laboratory, Parks Road, Oxford, U.K., (006). 15-0
17 Pricing Exotic Options using ML-MC MC Case 1: Quadratic Volatility Model ds = S ( µ dt + σ dw 1, t ) dσ = κ( ϖ σ) dt + β σ dw 1, t Monte Carlo paths Cost of simulation European Option: PCall = max S( T ) K,0 PPut = max K S( T ),0 For ε = , the ML-MC method is more efficient than the standard Euler method for: 50 times for Euler scheme. 150 times for Milstein scheme. 300 times for θ scheme. Value of θ for one random path Top left: convergence in option value with grid level. Bottom left: reduction in ML-MC variance. Top right: # of Monte Carlo paths required on each level. Bottom right: overall computational cost. θ ρ () t = log ( t ) log( S() t ) σ () ρ T = 1, ρ = 0.3, µ = 0.05, κ = 1.8, ϖ = 0.6 S(0) = 1, σ(0) = 0.1, β = 1, K = 1, put option 1 t+ t L ( 1,) = Lévy area = 0 t 16-0 β ρ
18 Pricing Exotic Options using ML-MC MC Case : GARCH Diffusion Model ds = S ( µ dt + σdw 1, t ) ( ) dν = κ ϖ ν dt + β νdw, t σ = ( ν ) max,0 Monte Carlo paths Digital Options P = H ST K ( ) H x = 1 if x>0 Heaviside function = H x = 0 if x 0 Cost of simulation For ε =0.0001, the ML-MC method is more efficient than the standard Euler method for: 3 times for Euler scheme. 60 times for Milstein scheme. 90 times for θ scheme. Top left: convergence in option value with grid level. Bottom left: reduction in ML-MC variance. Top right: # of Monte Carlo paths required on each level. Bottom right: overall computational cost. T = 1, ρ = 0.3, µ = 0.05, κ = 1.8, ϖ = 0.6 (0) = 1, (0) = 0.1, = 0.5 S ν β K = 1 t+ t Lévy area 0 1, t L = = 17-0
19 Pricing Exotic Options using ML-MC MC Case : GARCH Diffusion Model ds = S ( µ dt + σdw 1, t ) ( ) dν = κ ϖ ν dt + β νdw, t σ = ( ν ) max,0 Monte Carlo paths Asian Options P = max ST K,0 Arithmetic average: T N dt dt S = S ( t) dt S n + S n 0 ( ) n = 1 1 Cost of simulation For ε = , the ML-MC method is more efficient than the standard Euler method for: 50 times for Euler scheme. 80 times for Milstein scheme. 110 times for θ scheme. Top left: convergence in option value with grid level. Bottom left: reduction in ML-MC variance. Top right: # of Monte Carlo paths required on each level. Bottom right: overall computational cost. T = 1, ρ = 0.3, µ = 0.05, κ = 1.8, ϖ = 0.6 (0) = 1, (0) = 0.1, = 0.5 S ν β K = 1 t+ t Lévy area 0 1, t L = = 18-0
20 Pricing Exotic Options using ML-MC MC Case : GARCH Diffusion Model ds = S ( µ dt + σdw 1, t ) ( ) dν = κ ϖ ν dt + β νdw, t σ = ( ν ) max,0 Monte Carlo paths Variance Swaps P ( ν ) = N T K Arithmetic average: T N dt dt ν = ν ( t) dt ν n + ν n 0 Var n = 1 1 Cost of simulation For ε = , the ML-MC method is more efficient than the standard Euler method for: 150 times for Euler scheme. 380 times for Milstein scheme. 360 times for θ scheme. Top left: convergence in option value with grid level. Bottom left: reduction in ML-MC variance. Top right: # of Monte Carlo paths required on each level. Bottom right: overall computational cost. T ρ µ κ ϖ = 1, = 0.3, = 0.05, = 1.8, = 0.6 S(0) = 1, ν(0) = 0.1, β = 0.5 N = 10, K = 0.6 Var t+ t Lévy area 0 1, t L = = 19-0
21 Conclusions 1. θ scheme: When a specific orthogonal transformation is applied to a -Dimensional SDE and the conditions for integrability are satisfied, we can use the formula for θ to obtain the value of the rotation angle and obtain first order strong convergence without the simulation of expensive Lévy areas. Otherwise, we can use the 3-Dimensional transformation and check the magnitude of the Lie brackets to decide if it is likely to give computational savings in the solution of our system. The numerical results demonstrate considerable computational savings when the orthogonal transformation is applied to the quadratic volatility model or the 3/ model or the GARCH model. Unfortunately, similar savings are not achieved with the Heston model.. Multilevel Monte Carlo method: The ML-MC works without any problems with all schemes and does not depend on the value of the parameters of the system. In combination with θ scheme, it can reduce substantially the computational cost in pricing exotic options, reducing the cost to achieve an r.m.s. error of size ε from O(ε ³) to O(ε ²). Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk "Pricing Exotic Options using Strong Schemes, Oxford, UK 0-0
22 Thank you for your attention? Klaus E. Schmitz Abe "Pricing Exotic Options using Strong Schemes, Oxford, UK
23 Strong Convergence DEFINITION: We shall state that a time discrete approximation Ŝ(T) converges strongly with order γ>0 at time T as t 0 to the exact solution S(T) if there exists a positive constant "C", that does not depend on γ, such that: E S( T) S ( T, t) =C t γ (1) Kloeden & Platen If (1) is true, it can be implied that the order of strong convergence is not only in the last point T, but also uniformly within the whole time interval [0,T]. REMARK: When we talk about strong convergence, we are referring to how fast our time discrete approximation converges to the exact solution.
24 Weak Convergence DEFINITION: We shall state that a time discrete approximation S(T) converges weakly with order β>0 at time T as t 0 to the exact solution S(T) if there exists a positive constant "C", that does not depend on β, such that: ( ) ( (, )) = ( ) E g S T E g S T t C t β Monte Carlo Expectation: ( 1) p d (, ) g C β + M 1 E g S(T) - g S k (T, ) = M ( t ) O( t) k=1 REMARK: When we talk about weak convergence, we want to know how fast the expectation of our time discrete approximation converges to the exact expectation of our system.
25 Convergence without an Exact Solution Theorem 1: If a time discrete approximation Ŝ(T) converges strongly with order γ>0 at time T as t 0 to the exact solution S(T), i.e. exist a positive constant "C", that does not depend on γ such that: E S( T) S ( T, t) = C t γ Then using the same Wiener process, there exists a positive constant "C1", that is not dependent of γ, such that: t γ E S(T, t) S T, =C1 ( t) Surface error using the exact solution. Surface error using Theorem 1. [3] K. Schmitz-Abe and W. T. Shaw: Measure Order of Convergence without an Exact Solution, Euler vs Milstein Scheme. International Journal of Pure and Applied Mathematics, Vol. 4, No 3, , (005).
26 Conclusions (Paper) 1. In finance, stochastic variance and volatility models are very important for the valuation of exotic options. We have shown that the use of the orthogonal θ scheme can achieve the first order strong convergence properties of the Milstein numerical discretisation without the expensive simulation of Lévy areas. In combination with the recently introduced Multilevel Monte Carlo method it can reduce substantially the computational cost in pricing exotic options, reducing the cost to achieve an r.m.s. error of size ε from O(ε ³) to O(e ²).. The ML-MC works without any problems with all schemes and does not depend on the value of the parameters of the system. However, when a specific orthogonal transformation (θ scheme) is applied to a -Dimensional SDE it is only possible under certain conditions to avoid calculation of the Lévy area. When the conditions for integrability are satisfied, we can use the formula for θ to obtain the value of the rotation angle and obtain first order strong convergence. Otherwise, we have to use the 3-Dimensional transformation (θ scheme) and check the magnitude of the Lie brackets to decide if it is likely to give computational savings in the solution of our system. 3. The numerical results demonstrate considerable computational savings when the orthogonal transformation is applied to either the quadratic volatility model or the stochastic variance model. Unfortunately, similar savings are not achieved with the Heston model Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk "Pricing Exotic Options using Strong Schemes London, UK (007)
Monte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
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 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 informationMultilevel Monte Carlo Simulation
Multilevel Monte Carlo p. 1/48 Multilevel Monte Carlo Simulation Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance Workshop on Computational
More informationMultilevel quasi-monte Carlo path simulation
Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein
More informationMultilevel Monte Carlo for Basket Options
MLMC for basket options p. 1/26 Multilevel Monte Carlo for Basket Options Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance WSC09,
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 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 informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
More information"Vibrato" Monte Carlo evaluation of Greeks
"Vibrato" Monte Carlo evaluation of Greeks (Smoking Adjoints: part 3) Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance MCQMC 2008,
More informationLogarithmic derivatives of densities for jump processes
Logarithmic derivatives of densities for jump processes Atsushi AKEUCHI Osaka City University (JAPAN) June 3, 29 City University of Hong Kong Workshop on Stochastic Analysis and Finance (June 29 - July
More informationAnalytical formulas for local volatility model with stochastic. Mohammed Miri
Analytical formulas for local volatility model with stochastic rates Mohammed Miri Joint work with Eric Benhamou (Pricing Partners) and Emmanuel Gobet (Ecole Polytechnique Modeling and Managing Financial
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 informationModule 4: Monte Carlo path simulation
Module 4: Monte Carlo path simulation Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Module 4: Monte Carlo p. 1 SDE Path Simulation In Module 2, looked at the case
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 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 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 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 informationSample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models
Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Scott Robertson Carnegie Mellon University scottrob@andrew.cmu.edu http://www.math.cmu.edu/users/scottrob June
More informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
More informationPricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case
Pricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case Guang-Hua Lian Collaboration with Robert Elliott University of Adelaide Feb. 2, 2011 Robert Elliott,
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 informationMonte Carlo Methods. Prof. Mike Giles. Oxford University Mathematical Institute. Lecture 1 p. 1.
Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Lecture 1 p. 1 Geometric Brownian Motion In the case of Geometric Brownian Motion ds t = rs t dt+σs
More informationComputing Greeks with Multilevel Monte Carlo Methods using Importance Sampling
Computing Greeks with Multilevel Monte Carlo Methods using Importance Sampling Supervisor - Dr Lukas Szpruch Candidate Number - 605148 Dissertation for MSc Mathematical & Computational Finance Trinity
More informationTwo-dimensional COS method
Two-dimensional COS method Marjon Ruijter Winterschool Lunteren 22 January 2013 1/29 Introduction PhD student since October 2010 Prof.dr.ir. C.W. Oosterlee). CWI national research center for mathematics
More informationFinancial Mathematics and Supercomputing
GPU acceleration in early-exercise option valuation Álvaro Leitao and Cornelis W. Oosterlee Financial Mathematics and Supercomputing A Coruña - September 26, 2018 Á. Leitao & Kees Oosterlee SGBM on GPU
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 informationDr. Maddah ENMG 625 Financial Eng g II 10/16/06
Dr. Maddah ENMG 65 Financial Eng g II 10/16/06 Chapter 11 Models of Asset Dynamics () Random Walk A random process, z, is an additive process defined over times t 0, t 1,, t k, t k+1,, such that z( t )
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 informationMultilevel Monte Carlo for VaR
Multilevel Monte Carlo for VaR Mike Giles, Wenhui Gou, Abdul-Lateef Haji-Ali Mathematical Institute, University of Oxford (BNP Paribas, Hong Kong) (also discussions with Ralf Korn, Klaus Ritter) Advances
More informationAn overview of some financial models using BSDE with enlarged filtrations
An overview of some financial models using BSDE with enlarged filtrations Anne EYRAUD-LOISEL Workshop : Enlargement of Filtrations and Applications to Finance and Insurance May 31st - June 4th, 2010, Jena
More informationStochastic Runge Kutta Methods with the Constant Elasticity of Variance (CEV) Diffusion Model for Pricing Option
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 18, 849-856 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4381 Stochastic Runge Kutta Methods with the Constant Elasticity of Variance
More informationModule 2: Monte Carlo Methods
Module 2: Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute MC Lecture 2 p. 1 Greeks In Monte Carlo applications we don t just want to know the expected
More informationExact Sampling of Jump-Diffusion Processes
1 Exact Sampling of Jump-Diffusion Processes and Dmitry Smelov Management Science & Engineering Stanford University Exact Sampling of Jump-Diffusion Processes 2 Jump-Diffusion Processes Ubiquitous in finance
More informationAn Efficient Numerical Scheme for Simulation of Mean-reverting Square-root Diffusions
Journal of Numerical Mathematics and Stochastics,1 (1) : 45-55, 2009 http://www.jnmas.org/jnmas1-5.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 An Efficient Numerical Scheme for Simulation of
More informationMultilevel path simulation for jump-diffusion SDEs
Multilevel path simulation for jump-diffusion SDEs Yuan Xia, Michael B. Giles Abstract We investigate the extension of the multilevel Monte Carlo path simulation method to jump-diffusion SDEs. We consider
More informationNumerical Solution of Stochastic Differential Equations with Jumps in Finance
Numerical Solution of Stochastic Differential Equations with Jumps in Finance Eckhard Platen School of Finance and Economics and School of Mathematical Sciences University of Technology, Sydney Kloeden,
More informationMonte Carlo Simulation of Stochastic Processes
Monte Carlo Simulation of Stochastic Processes Last update: January 10th, 2004. In this section is presented the steps to perform the simulation of the main stochastic processes used in real options applications,
More informationA Moment Matching Approach To The Valuation Of A Volume Weighted Average Price Option
A Moment Matching Approach To The Valuation Of A Volume Weighted Average Price Option Antony Stace Department of Mathematics and MASCOS University of Queensland 15th October 2004 AUSTRALIAN RESEARCH COUNCIL
More information3.1 Itô s Lemma for Continuous Stochastic Variables
Lecture 3 Log Normal Distribution 3.1 Itô s Lemma for Continuous Stochastic Variables Mathematical Finance is about pricing (or valuing) financial contracts, and in particular those contracts which depend
More informationApplication of Stochastic Calculus to Price a Quanto Spread
Application of Stochastic Calculus to Price a Quanto Spread Christopher Ting http://www.mysmu.edu/faculty/christophert/ Algorithmic Quantitative Finance July 15, 2017 Christopher Ting July 15, 2017 1/33
More informationConvergence Analysis of Monte Carlo Calibration of Financial Market Models
Analysis of Monte Carlo Calibration of Financial Market Models Christoph Käbe Universität Trier Workshop on PDE Constrained Optimization of Certain and Uncertain Processes June 03, 2009 Monte Carlo Calibration
More informationAnalysing multi-level Monte Carlo for options with non-globally Lipschitz payoff
Finance Stoch 2009 13: 403 413 DOI 10.1007/s00780-009-0092-1 Analysing multi-level Monte Carlo for options with non-globally Lipschitz payoff Michael B. Giles Desmond J. Higham Xuerong Mao Received: 1
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Spring 2010 Computer Exercise 2 Simulation This lab deals with
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationOptimal Stopping for American Type Options
Optimal Stopping for Department of Mathematics Stockholm University Sweden E-mail: silvestrov@math.su.se ISI 2011, Dublin, 21-26 August 2011 Outline of communication Multivariate Modulated Markov price
More informationThe Evaluation of American Compound Option Prices under Stochastic Volatility. Carl Chiarella and Boda Kang
The Evaluation of American Compound Option Prices under Stochastic Volatility Carl Chiarella and Boda Kang School of Finance and Economics University of Technology, Sydney CNR-IMATI Finance Day Wednesday,
More informationValuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model
Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility
More informationAdvanced computational methods X SDE Lecture 5
Advanced computational methods X071521-SDE Lecture 5 1 Other weak schemes Here, I list out some typical weak schemes. If you are interested in them, you can read in details. 1.1 Weak Taylor approximations
More informationSimulating more interesting stochastic processes
Chapter 7 Simulating more interesting stochastic processes 7. Generating correlated random variables The lectures contained a lot of motivation and pictures. We'll boil everything down to pure algebra
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
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 informationTime-changed Brownian motion and option pricing
Time-changed Brownian motion and option pricing Peter Hieber Chair of Mathematical Finance, TU Munich 6th AMaMeF Warsaw, June 13th 2013 Partially joint with Marcos Escobar (RU Toronto), Matthias Scherer
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 informationHeston Stochastic Local Volatility Model
Heston Stochastic Local Volatility Model Klaus Spanderen 1 R/Finance 2016 University of Illinois, Chicago May 20-21, 2016 1 Joint work with Johannes Göttker-Schnetmann Klaus Spanderen Heston Stochastic
More informationParallel Multilevel Monte Carlo Simulation
Parallel Simulation Mathematisches Institut Goethe-Universität Frankfurt am Main Advances in Financial Mathematics Paris January 7-10, 2014 Simulation Outline 1 Monte Carlo 2 3 4 Algorithm Numerical Results
More informationA Continuity Correction under Jump-Diffusion Models with Applications in Finance
A Continuity Correction under Jump-Diffusion Models with Applications in Finance Cheng-Der Fuh 1, Sheng-Feng Luo 2 and Ju-Fang Yen 3 1 Institute of Statistical Science, Academia Sinica, and Graduate Institute
More informationUniversity of Oxford. Robust hedging of digital double touch barrier options. Ni Hao
University of Oxford Robust hedging of digital double touch barrier options Ni Hao Lady Margaret Hall MSc in Mathematical and Computational Finance Supervisor: Dr Jan Ob lój Oxford, June of 2009 Contents
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 informationLocal Volatility Dynamic Models
René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding
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 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 informationRough volatility models
Mohrenstrasse 39 10117 Berlin Germany Tel. +49 30 20372 0 www.wias-berlin.de October 18, 2018 Weierstrass Institute for Applied Analysis and Stochastics Rough volatility models Christian Bayer EMEA Quant
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 informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More informationInvestigation of Dependency between Short Rate and Transition Rate on Pension Buy-outs. Arık, A. 1 Yolcu-Okur, Y. 2 Uğur Ö. 2
Investigation of Dependency between Short Rate and Transition Rate on Pension Buy-outs Arık, A. 1 Yolcu-Okur, Y. 2 Uğur Ö. 2 1 Hacettepe University Department of Actuarial Sciences 06800, TURKEY 2 Middle
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationPricing and hedging with rough-heston models
Pricing and hedging with rough-heston models Omar El Euch, Mathieu Rosenbaum Ecole Polytechnique 1 January 216 El Euch, Rosenbaum Pricing and hedging with rough-heston models 1 Table of contents Introduction
More informationMultiname and Multiscale Default Modeling
Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis
More informationSTOCHASTIC INTEGRALS
Stat 391/FinMath 346 Lecture 8 STOCHASTIC INTEGRALS X t = CONTINUOUS PROCESS θ t = PORTFOLIO: #X t HELD AT t { St : STOCK PRICE M t : MG W t : BROWNIAN MOTION DISCRETE TIME: = t < t 1
More informationAn Asymptotic Expansion Formula for Up-and-Out Barrier Option Price under Stochastic Volatility Model
CIRJE-F-873 An Asymptotic Expansion Formula for Up-and-Out Option Price under Stochastic Volatility Model Takashi Kato Osaka University Akihiko Takahashi University of Tokyo Toshihiro Yamada Graduate School
More informationCounterparty Credit Risk Simulation
Counterparty Credit Risk Simulation Alex Yang FinPricing http://www.finpricing.com Summary Counterparty Credit Risk Definition Counterparty Credit Risk Measures Monte Carlo Simulation Interest Rate Curve
More informationLecture 3. Sergei Fedotov Introduction to Financial Mathematics. Sergei Fedotov (University of Manchester) / 6
Lecture 3 Sergei Fedotov 091 - Introduction to Financial Mathematics Sergei Fedotov (University of Manchester) 091 010 1 / 6 Lecture 3 1 Distribution for lns(t) Solution to Stochastic Differential Equation
More informationOn VIX Futures in the rough Bergomi model
On VIX Futures in the rough Bergomi model Oberwolfach Research Institute for Mathematics, February 28, 2017 joint work with Antoine Jacquier and Claude Martini Contents VIX future dynamics under rbergomi
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 informationValue at Risk Ch.12. PAK Study Manual
Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and
More informationStochastic Volatility (Working Draft I)
Stochastic Volatility (Working Draft I) Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu 1 Introduction When using the Black-Scholes-Merton model to price derivative
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 informationAnumericalalgorithm for general HJB equations : a jump-constrained BSDE approach
Anumericalalgorithm for general HJB equations : a jump-constrained BSDE approach Nicolas Langrené Univ. Paris Diderot - Sorbonne Paris Cité, LPMA, FiME Joint work with Idris Kharroubi (Paris Dauphine),
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 informationApplication of Moment Expansion Method to Option Square Root Model
Application of Moment Expansion Method to Option Square Root Model Yun Zhou Advisor: Professor Steve Heston University of Maryland May 5, 2009 1 / 19 Motivation Black-Scholes Model successfully explain
More informationOption Pricing Models for European Options
Chapter 2 Option Pricing Models for European Options 2.1 Continuous-time Model: Black-Scholes Model 2.1.1 Black-Scholes Assumptions We list the assumptions that we make for most of this notes. 1. The underlying
More informationNumerical Methods for Stochastic Differential Equations with Applications to Finance
Numerical Methods for Stochastic Differential Equations with Applications to Finance Matilde Lopes Rosa Instituto Superior Técnico University of Lisbon, Portugal May 2016 Abstract The pricing of financial
More informationA GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies
A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL Stephen Chin and Daniel Dufresne Centre for Actuarial Studies University of Melbourne Paper: http://mercury.ecom.unimelb.edu.au/site/actwww/wps2009/no181.pdf
More informationMath 623 (IOE 623), Winter 2008: Final exam
Math 623 (IOE 623), Winter 2008: Final exam Name: Student ID: This is a closed book exam. You may bring up to ten one sided A4 pages of notes to the exam. You may also use a calculator but not its memory
More informationMachine Learning for Quantitative Finance
Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY
ESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY Kai Detlefsen Wolfgang K. Härdle Rouslan A. Moro, Deutsches Institut für Wirtschaftsforschung (DIW) Center for Applied Statistics
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 informationVariance Reduction Through Multilevel Monte Carlo Path Calculations
Variance Reduction Through Mutieve Monte Caro Path Cacuations Mike Gies gies@comab.ox.ac.uk Oxford University Computing Laboratory Mutieve Monte Caro p. 1/30 Mutigrid A powerfu technique for soving PDE
More informationReplication and Absence of Arbitrage in Non-Semimartingale Models
Replication and Absence of Arbitrage in Non-Semimartingale Models Matematiikan päivät, Tampere, 4-5. January 2006 Tommi Sottinen University of Helsinki 4.1.2006 Outline 1. The classical pricing model:
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 informationMultiscale Stochastic Volatility Models
Multiscale Stochastic Volatility Models Jean-Pierre Fouque University of California Santa Barbara 6th World Congress of the Bachelier Finance Society Toronto, June 25, 2010 Multiscale Stochastic Volatility
More informationValuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005
Valuation of Volatility Derivatives Jim Gatheral Global Derivatives & Risk Management 005 Paris May 4, 005 he opinions expressed in this presentation are those of the author alone, and do not necessarily
More informationOptimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University
Optimal Hedging of Variance Derivatives John Crosby Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation at Baruch College, in New York, 16th November 2010
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 informationNumerical Solution of Stochastic Differential Equations with Jumps in Finance
Numerical Solution of Stochastic Differential Equations with Jumps in Finance Eckhard Platen School of Finance and Economics and School of Mathematical Sciences University of Technology, Sydney Kloeden,
More informationOptimal robust bounds for variance options and asymptotically extreme models
Optimal robust bounds for variance options and asymptotically extreme models Alexander Cox 1 Jiajie Wang 2 1 University of Bath 2 Università di Roma La Sapienza Advances in Financial Mathematics, 9th January,
More information