Convergence Analysis of Monte Carlo Calibration of Financial Market Models
|
|
- Denis Ball
- 6 years ago
- Views:
Transcription
1 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
2 Monte Carlo Calibration Fundamentals Focus will be on calibration of European call options Definition 1 (European Call Option) A European call option is the right to buy a predetermined underlying (e.g. stock) at a certain time T (maturity) for a certain price K (strike).
3 Monte Carlo Calibration Fundamentals Focus will be on calibration of European call options Definition 1 (European Call Option) A European call option is the right to buy a predetermined underlying (e.g. stock) at a certain time T (maturity) for a certain price K (strike). Definition 2 (Price of a Call Option) The price of a call option C in t = 0 can be calculated through C = e rt E (max(s T K, 0)) where r is the risk free rate and S T the value of the underlying at future time T.
4 Monte Carlo Calibration Stochastic Differential Equation L-dimensional system of stochastic differential equations (SDE): dy t (x) = a(x, Y t (x))dt + b(x, Y t (x))dw t where x R P Y t = [S t, Yt 2,..., Yt L ] R L W t = (Wt 1,..., Wt L ) R L a : R P R L R L b : R P R L R L R L vector of parameters Solution of SDE Vector of Brownian motions a l (x, Y t (x))dt L ν=1 bl,ν (x, Y t (x))dwt ν, l = 1,..., L
5 Monte Carlo Calibration Least Squares Problem Continuous Optimization Problem (True Problem) min f(x) := I ( C i (x) C i 2 x X obs) i=1 where C i (x) = e rt i E (max(s Ti (x) K i, 0)) s.t. dy t (x) = a(x, Y t (x))dt + b(x, Y t (x))dw t, Y 0 > 0 X R P convex and compact
6 Monte Carlo Calibration Least Squares Problem Continuous Optimization Problem (True Problem) min f(x) := I ( C i (x) C i 2 x X obs) i=1 where C i (x) = e rt i E (max(s Ti (x) K i, 0)) s.t. dy t (x) = a(x, Y t (x))dt + b(x, Y t (x))dw t, Y 0 > 0 X R P convex and compact Discretized Optimization Problem (SAA Problem) min f M, t,ɛ := I ( ) 2 CM, t,ɛ i (x) x X Ci obs i=1 M where C i M, t,ɛ (x) := e rt i 1 M s.t. m=1 ( ) π ɛ (s m N i,ɛ (x) K i) y m n+1,ɛ (x) = ym n,ɛ(x) + a ɛ (x, y m n,ɛ(x)) t n + b ɛ (x, y m n,ɛ(x)) W m n
7 Monte Carlo Calibration Smoothing Non-differentiabilities Consider Heston s Model: ds t,ɛ = (r δ)s t,ɛ dt + v + t,ɛ S t,ɛdw 1 t dv t,ɛ = κ(θ v + t,ɛ )dt + σ v + t,ɛ (ρdw 1 t + 1 ρ 2 dw 2 t )
8 Overview Table of Contents 1 Monte Carlo Calibration 2 Convergence Overview Pathwise Uniqueness Uniform Convergence First Order Optimality Condition 3 Conclusions
9 Overview Convergence True Problem min f(x) := I ( C i (x) C i 2 x X obs) i=1 SAA Problem min f M k, t k,ɛ k := I x X i=1 Increase number of simulations: ( ) 2 CM i k, t k,ɛ k (x) Cobs i M k Decrease discretization step size: t k 0 Decrease smoothing parameter: ɛ k 0 x k X solutions
10 Overview Convergence True Problem min f(x) := I ( C i (x) C i 2 x X obs) i=1 SAA Problem min f M k, t k,ɛ k := I x X i=1 Increase number of simulations: ( ) 2 CM i k, t k,ɛ k (x) Cobs i M k Decrease discretization step size: t k 0 Decrease smoothing parameter: ɛ k 0 X compact x kl x with x in X. x k X solutions
11 Overview Convergence True Problem min f(x) := I ( C i (x) C i 2 x X obs) i=1 SAA Problem min f M k, t k,ɛ k := I x X i=1 Increase number of simulations: ( ) 2 CM i k, t k,ɛ k (x) Cobs i M k Decrease discretization step size: t k 0 Decrease smoothing parameter: ɛ k 0 X compact x kl x with x in X. Question x solution of the true problem? x k X solutions
12 Overview Local Minima min f(x) := x2 x [ 1;1] min f M(x) := x 2 M 1 sin(mx 2 ) x [ 1;1]
13 Overview Local Minima min f(x) := x2 x [ 1;1] min f M(x) := x 2 M 1 sin(mx 2 ) x [ 1;1] Local minima might lead to problems:
14 Overview Literature Review True Problem: min h(x) := E(H(x, ω)) min x X SAA Problem x X h M(x) := 1 M M m=1 H(x, ω m) Shapiro (2000): Convergence if min h(x) produces global minimum Rubinstein & Shapiro (1993): Convergence to a critical first order point under assumption that H(x, ω) is dominated integrable and continuous Bastin et al. (2006): Additionally second order convergence even for stochastic constraints Dependence on three error sources: Monte Carlo, discretization and smoothing!
15 Overview Goal: First Order Optimality Steps to be taken: 1 Pathwise Uniqueness of SDE 2 Uniform Convergence: lim sup f Mk, t k,ɛ k (x) f(x) = 0 k x X lim sup f Mk, t k,ɛ k (x) f(x) = 0 k x X 3 First Order Optimality Condition: f(x ) T (x x ) 0 x X
16 Pathwise Uniqueness Table of Contents 1 Monte Carlo Calibration 2 Convergence Overview Pathwise Uniqueness Uniform Convergence First Order Optimality Condition 3 Conclusions
17 Pathwise Uniqueness Pathwise Uniqueness under Lipschitz Continuity Theorem 3 (Kloeden & Platen) Under the assumptions that There exists a constant K Lip > 0 such that t [0, T ] and y R L a(t, y) a(t, z) + b(t, y) b(t, z) K Lip y z There exists a constant K Grow > 0 such that t [0, T ] and y R L a(t, y) + b(t, y) K Grow (1 + y ) the stochastic differential equation dy t = a(t, Y t )dt + b(t, Y t )dw t, Y 0 (0, ). has a pathwise unique strong solution Y t on [0, T ].
18 Pathwise Uniqueness Problem: Lipschitz Continuity Consider Heston s model ds t,ɛ = (r δ)s t,ɛ dt + π ɛ (v t,ɛ )S t,ɛ dwt 1 dv t,ɛ = κ(θ π ɛ (v t,ɛ ))dt + σ π ɛ (v t,ɛ )(ρdwt ρ 2 dwt 2 ) Lipschitz continuity for ɛ > 0
19 Pathwise Uniqueness Yamada Condition Theorem 4 Let with dy t,ɛ = a ɛ (t, Y t,ɛ )dt + b ɛ (t, Y t,ɛ )dw t. a ɛ (t, Y t,ɛ ) = (a 1 ɛ(t, Y 1 t,ɛ),..., a L ɛ (t, Y L t,ɛ)) T b ɛ (t, Y t,ɛ ) = diag(b 1 ɛ(t, Y 1 t,ɛ),..., b L ɛ (t, Y L t,ɛ)) If there exists a positive increasing function β : [0, ) [0, ) with and b i (t, x) b i (t, y) β( x y ) x, y R, i = 1,..., L δ with an arbitrarily small δ > β 2 (z)dz =.
20 Pathwise Uniqueness Yamada Condition (2)... and a positive increasing concave function α : [0, ) [0, ) such that with a i (t, x) a i (t, y) α( x y ) x, y R, i = 1,..., L δ 0 α 1 (z)dz =. with an arbitrarily small δ > 0, the SDE has a pathwise unique solution. Proof: Yamada (1971)
21 Pathwise Uniqueness Yamada Condition (3) Reconsider Heston s model ds t,ɛ = (r δ)s t,ɛ dt + π ɛ (v t,ɛ )S t,ɛ dwt 1 dv t,ɛ = κ(θ π ɛ (v t,ɛ ))dt + σ π ɛ (v t,ɛ )(ρdwt ρ 2 dwt 2 ) The drift is Lipschitz continuous: a i (t, x) a i (t, y) K Lip x y x, y R, i = 1, 2 and the diffusion is Hölder continuous: b i (t, x) b i (t, y) x y x, y R, i = 1, 2 with δ 0 1 dz = ; K Lip z δ 0 1 dz =. z
22 Pathwise Uniqueness Problem: Independent components required Heston s model: ds t,ɛ = (r δ)s t,ɛ dt + π ɛ (v t,ɛ )S t,ɛ dwt 1 dv t,ɛ = κ(θ π ɛ (v t,ɛ ))dt + σ π ɛ (v t,ɛ )dwt 2
23 Pathwise Uniqueness Problem: Independent components required Heston s model: ds t,ɛ = (r δ)s t,ɛ dt + π ɛ (v t,ɛ )S t,ɛ dwt 1 dv t,ɛ = κ(θ π ɛ (v t,ɛ ))dt + σ π ɛ (v t,ɛ )dwt 2 Solution: Process v t,ɛ has pathwise unique solution following Yamada s Theorem Insert this unique solution in process S t,ɛ Process S t,ɛ has pathwise unique solution following Yamada s Theorem Pathwise unique solution via Yamada s Theorem
24 Uniform Convergence Table of Contents 1 Monte Carlo Calibration 2 Convergence Overview Pathwise Uniqueness Uniform Convergence First Order Optimality Condition 3 Conclusions
25 Uniform Convergence Convergence of the Problem Reconsider: f M, t,ɛ (x) f(x) f M, t,ɛ (x) f t,ɛ (x) (1) + f t,ɛ (x) f ɛ (x) (2) + f ɛ (x) f(x) (3)
26 Uniform Convergence Convergence of the Problem Reconsider: f M, t,ɛ (x) f(x) f M, t,ɛ (x) f t,ɛ (x) (1) + f t,ɛ (x) f ɛ (x) (2) + f ɛ (x) f(x) (3) Assumption: There exists a constant K Grow > 0 such that t [0, T ] and y R L a ɛ (t, y) + b ɛ (t, y) K Grow (1 + y ).
27 Uniform Convergence Convergence of Smoothed and Discretized SDE Theorem 5 Consider the SDE dy t,ɛ = a ɛ (t, Y t,ɛ )dt + b ɛ (t, Y t,ɛ )dw t. and the continuously interpolated process t y t,ɛ = Y t a ɛ (x, y τ(s),ɛ )ds + 0 b ɛ (x, y τ(s),ɛ )dw s where τ(s) = n, s [τ n, τ n+1 ) and n = 0,..., N 1. Assuming that the growth condition holds and the SDE has a pathwise unique solution it holds ( lim sup E y t,ɛ Y T,ɛ 2) = 0. t 0 x X Proof: Kaneko & Nakao (1988)
28 Uniform Convergence Convergence of Smoothed SDE Theorem 6 Assume that the growth condition and the pathwise uniqueness holds for a solution of dy t = a(t, Y t )dt + b(t, Y t )dw t and let Y t,ɛ be a solution of dy t,ɛ = a ɛ (t, Y t,ɛ )dt + b ɛ (t, Y t,ɛ )dw t. If a ɛ and b ɛ converge uniformly to a and b for ɛ 0, i.e. lim sup ɛ 0 t [0,T ] x X sup ( a ɛ (t, x) a(t, x) + b ɛ (t, x) b(t, x) ) = 0. where is a matrix norm, it holds lim sup E ɛ 0 x X Proof: Kaneko & Nakao (1988) ( Y t,ɛ Y t 2) = 0.
29 Uniform Convergence Dominated Integrability & Continuity Lemma 7 Assume that the families {π(s T (x, ω) K), x X} are dominated by a Q-integrable function P (ω). Then there exist t > 0 and ɛ > 0 such that {π ɛ (s N,ɛ (x, ω) K), x X} is dominated by a Q-integrable function for all t [0, t] and ɛ [0, ɛ]. Lemma 8 If the functions π(s T (, ω) K) are continuous on X for Q almost every ω, the functions π ɛ (s N,ɛ (x, ω) K) are continuous on X for 0 < t < and 0 < ɛ <.
30 Uniform Convergence Uniform Convergence Theorem 9 Assume that the families {π(s T (x, ω) K), x X} are dominated by a Q-integrable function P (ω) and furthermore the functions π(s T (, ω) K) are continuous on X for Q almost every ω. If additionally X is compact, then f(x) is continuous on X. Furthermore f M, t,ɛ converges uniformly to f on X, i.e. for given sequences (M k ) k IN, ( t k ) k R + and (ɛ k ) k R + satisfying M k, t k 0, ɛ k 0 it holds lim sup f Mk, t k,ɛ k (x) f(x) = 0. k x X Note that the same can be shown for the gradients!
31 First Order Optimality Condition Table of Contents 1 Monte Carlo Calibration 2 Convergence Overview Pathwise Uniqueness Uniform Convergence First Order Optimality Condition 3 Conclusions
32 First Order Optimality Condition First Order Optimality Condition Theorem 10 Assume that the families {π(s T (x, ω) K), x X} and { x p π(s T (, ω) K), x X}, i = 1,..., I are dominated by a Q-integrable function P (ω) and furthermore the functions π(s T (, ω) K) and x p π(s T (, ω) K), i = 1,..., I are continuous on X for Q almost every ω and additionally that X is compact. Further let (M k ) k N +, ( t k ) k R +, (ɛ k ) k R + and (γ k ) k R + with M k, t k 0, ɛ k 0 and γ k 0 be given sequences and assume that (x k ) k IN X is a sequence of points satisfying f(x k ) T (x x k ) γ k x X. Then every limit point x X of (x k ) k almost surely satisfies the first order optimality condition f(x ) T (x x ) 0 x X.
33 First Order Optimality Condition Convergence: Graphical Illustration
34 First Order Optimality Condition Convergence: Graphical Illustration
35 First Order Optimality Condition Convergence: Graphical Illustration
36 Conclusions Table of Contents 1 Monte Carlo Calibration 2 Convergence Overview Pathwise Uniqueness Uniform Convergence First Order Optimality Condition 3 Conclusions
37 Conclusions Conclusions Set up calibration problem Discretized via Monte Carlo, Euler-Maruyama and smoothing Pathwise Uniqueness for resulting SDE under Yamada Condition Uniform convergence of objectives under unrestrictive assumptions First order optimality condition satisfied for limit point x
38 Conclusions Bibliography Bastin,F., Cirillo,C. and Toint,P.L.: Convergence Theory for Nonconvex Stochastic Programming with an Application to Mixed Logit, Mathematical Programming Series B, Vol. 108, 2006, Rubinstein,R.Y. and Shapiro,A.: Discrete Event Systems, John Wiley, 1993 Shapiro, A.: Stochastic Programming by Monte Carlo Simulation Methods, Stochastic Programming E-Print Series 2000, Kaneko,H. and Nakao,S.: A Note on Approximation for Stochastic Differential Equations, Seminaire de Probabilites, XXII, Lecture Notes in Mathematics, Vol. 1321, 1988, Yamada, T. and Watanabe, S.: On the Uniqueness of Solutions of Stochastic Differential Equations, Journal of Mathematics of Kyoto University, Vol 11, 1971 Kaebe, C., Maruhn, J. and Sachs, E.W.: Adjoint Based Monte Carlo Calibration of Financial Market Models, Journal of Finance and Stochastics (to appear)
AMH4 - 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 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 informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationMLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th 2017 1 https://papers.ssrn.com/sol3/papers.cfm?abstract
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 informationCalibration Lecture 1: Background and Parametric Models
Calibration Lecture 1: Background and Parametric Models March 2016 Motivation What is calibration? Derivative pricing models depend on parameters: Black-Scholes σ, interest rate r, Heston reversion speed
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 informationMonte 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 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 informationLecture 8: The Black-Scholes theory
Lecture 8: The Black-Scholes theory Dr. Roman V Belavkin MSO4112 Contents 1 Geometric Brownian motion 1 2 The Black-Scholes pricing 2 3 The Black-Scholes equation 3 References 5 1 Geometric Brownian motion
More informationConvergence of Discretized Stochastic (Interest Rate) Processes with Stochastic Drift Term.
Convergence of Discretized Stochastic (Interest Rate) Processes with Stochastic Drift Term. G. Deelstra F. Delbaen Free University of Brussels, Department of Mathematics, Pleinlaan 2, B-15 Brussels, Belgium
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 informationROBUST STATIC HEDGING OF BARRIER OPTIONS IN STOCHASTIC VOLATILITY MODELS
ROBUST STATIC HEDGING OF BARRIER OPTIONS IN STOCHASTIC VOLATILITY MODELS J. H. MARUHN AND E. W. SACHS Abstract. Static hedge portfolios for barrier options are extremely sensitive with respect to changes
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 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 informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
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
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 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 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 informationIntroduction to Affine Processes. Applications to Mathematical Finance
and Its Applications to Mathematical Finance Department of Mathematical Science, KAIST Workshop for Young Mathematicians in Korea, 2010 Outline Motivation 1 Motivation 2 Preliminary : Stochastic Calculus
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 informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationThe Correlation Smile Recovery
Fortis Bank Equity & Credit Derivatives Quantitative Research The Correlation Smile Recovery E. Vandenbrande, A. Vandendorpe, Y. Nesterov, P. Van Dooren draft version : March 2, 2009 1 Introduction Pricing
More informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More informationEquivalence between Semimartingales and Itô Processes
International Journal of Mathematical Analysis Vol. 9, 215, no. 16, 787-791 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.215.411358 Equivalence between Semimartingales and Itô Processes
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
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 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 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 informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
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 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 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 informationVariance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment
Variance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment Jean-Pierre Fouque Tracey Andrew Tullie December 11, 21 Abstract We propose a variance reduction method for Monte Carlo
More informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
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 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 informationCS 774 Project: Fall 2009 Version: November 27, 2009
CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS
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 informationAmerican Foreign Exchange Options and some Continuity Estimates of the Optimal Exercise Boundary with respect to Volatility
American Foreign Exchange Options and some Continuity Estimates of the Optimal Exercise Boundary with respect to Volatility Nasir Rehman Allam Iqbal Open University Islamabad, Pakistan. Outline Mathematical
More informationParameter sensitivity of CIR process
Parameter sensitivity of CIR process Sidi Mohamed Ould Aly To cite this version: Sidi Mohamed Ould Aly. Parameter sensitivity of CIR process. Electronic Communications in Probability, Institute of Mathematical
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 informationStochastic Partial Differential Equations and Portfolio Choice. Crete, May Thaleia Zariphopoulou
Stochastic Partial Differential Equations and Portfolio Choice Crete, May 2011 Thaleia Zariphopoulou Oxford-Man Institute and Mathematical Institute University of Oxford and Mathematics and IROM, The University
More informationA new approach for scenario generation in risk management
A new approach for scenario generation in risk management Josef Teichmann TU Wien Vienna, March 2009 Scenario generators Scenarios of risk factors are needed for the daily risk analysis (1D and 10D ahead)
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 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 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 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 informationRegression estimation in continuous time with a view towards pricing Bermudan options
with a view towards pricing Bermudan options Tagung des SFB 649 Ökonomisches Risiko in Motzen 04.-06.06.2009 Financial engineering in times of financial crisis Derivate... süßes Gift für die Spekulanten
More informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
More informationPricing theory of financial derivatives
Pricing theory of financial derivatives One-period securities model S denotes the price process {S(t) : t = 0, 1}, where S(t) = (S 1 (t) S 2 (t) S M (t)). Here, M is the number of securities. At t = 1,
More informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
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 informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
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 informationSPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin
SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market
More informationLast Time. Martingale inequalities Martingale convergence theorem Uniformly integrable martingales. Today s lecture: Sections 4.4.1, 5.
MATH136/STAT219 Lecture 21, November 12, 2008 p. 1/11 Last Time Martingale inequalities Martingale convergence theorem Uniformly integrable martingales Today s lecture: Sections 4.4.1, 5.3 MATH136/STAT219
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 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 informationNo-arbitrage theorem for multi-factor uncertain stock model with floating interest rate
Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer
More informationAsset Pricing Models with Underlying Time-varying Lévy Processes
Asset Pricing Models with Underlying Time-varying Lévy Processes Stochastics & Computational Finance 2015 Xuecan CUI Jang SCHILTZ University of Luxembourg July 9, 2015 Xuecan CUI, Jang SCHILTZ University
More informationRobust Pricing and Hedging of Options on Variance
Robust Pricing and Hedging of Options on Variance Alexander Cox Jiajie Wang University of Bath Bachelier 21, Toronto Financial Setting Option priced on an underlying asset S t Dynamics of S t unspecified,
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 informationExam Quantitative Finance (35V5A1)
Exam Quantitative Finance (35V5A1) Part I: Discrete-time finance Exercise 1 (20 points) a. Provide the definition of the pricing kernel k q. Relate this pricing kernel to the set of discount factors D
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 informationRichardson Extrapolation Techniques for the Pricing of American-style Options
Richardson Extrapolation Techniques for the Pricing of American-style Options June 1, 2005 Abstract Richardson Extrapolation Techniques for the Pricing of American-style Options In this paper we re-examine
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 informationOutline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0.
Outline Coordinate Minimization Daniel P. Robinson Department of Applied Mathematics and Statistics Johns Hopkins University November 27, 208 Introduction 2 Algorithms Cyclic order with exact minimization
More informationWhat can we do with numerical optimization?
Optimization motivation and background Eddie Wadbro Introduction to PDE Constrained Optimization, 2016 February 15 16, 2016 Eddie Wadbro, Introduction to PDE Constrained Optimization, February 15 16, 2016
More informationYoungrok Lee and Jaesung Lee
orean J. Math. 3 015, No. 1, pp. 81 91 http://dx.doi.org/10.11568/kjm.015.3.1.81 LOCAL VOLATILITY FOR QUANTO OPTION PRICES WITH STOCHASTIC INTEREST RATES Youngrok Lee and Jaesung Lee Abstract. This paper
More informationA discretionary stopping problem with applications to the optimal timing of investment decisions.
A discretionary stopping problem with applications to the optimal timing of investment decisions. Timothy Johnson Department of Mathematics King s College London The Strand London WC2R 2LS, UK Tuesday,
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 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 of Interest Rates
WDS'12 Proceedings of Contributed Papers, Part I, 25 30, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Calibration of Interest Rates J. Černý Charles University, Faculty of Mathematics and Physics, Prague,
More informationMartingales & Strict Local Martingales PDE & Probability Methods INRIA, Sophia-Antipolis
Martingales & Strict Local Martingales PDE & Probability Methods INRIA, Sophia-Antipolis Philip Protter, Columbia University Based on work with Aditi Dandapani, 2016 Columbia PhD, now at ETH, Zurich March
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 information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More informationOptimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error
Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error José E. Figueroa-López Department of Mathematics Washington University in St. Louis Spring Central Sectional Meeting
More informationA Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 129 136 A Numerical Approach to the Estimation of Search Effort in a Search for a Moving
More informationThe Capital Asset Pricing Model as a corollary of the Black Scholes model
he Capital Asset Pricing Model as a corollary of the Black Scholes model Vladimir Vovk he Game-heoretic Probability and Finance Project Working Paper #39 September 6, 011 Project web site: http://www.probabilityandfinance.com
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 1, Part 1. Overview of Lecture 1, Part 1: Background Mater.
Numerical Simulation of Stochastic Differential Equations: Lecture, Part Des Higham Department of Mathematics University of Strathclyde Course Aim: Give an accessible intro. to SDEs and their numerical
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 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 informationOn Complexity of Multistage Stochastic Programs
On Complexity of Multistage Stochastic Programs Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA e-mail: ashapiro@isye.gatech.edu
More informationMartingale representation theorem
Martingale representation theorem Ω = C[, T ], F T = smallest σ-field with respect to which B s are all measurable, s T, P the Wiener measure, B t = Brownian motion M t square integrable martingale with
More informationAnalysis of pricing American options on the maximum (minimum) of two risk assets
Interfaces Free Boundaries 4, (00) 7 46 Analysis of pricing American options on the maximum (minimum) of two risk assets LISHANG JIANG Institute of Mathematics, Tongji University, People s Republic of
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 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 informationLocally risk-minimizing vs. -hedging in stochastic vola
Locally risk-minimizing vs. -hedging in stochastic volatility models University of St. Andrews School of Economics and Finance August 29, 2007 joint work with R. Poulsen ( Kopenhagen )and K.R.Schenk-Hoppe
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 informationRemarks on stochastic automatic adjoint differentiation and financial models calibration
arxiv:1901.04200v1 [q-fin.cp] 14 Jan 2019 Remarks on stochastic automatic adjoint differentiation and financial models calibration Dmitri Goloubentcev, Evgeny Lakshtanov Abstract In this work, we discuss
More informationEuropean option pricing under parameter uncertainty
European option pricing under parameter uncertainty Martin Jönsson (joint work with Samuel Cohen) University of Oxford Workshop on BSDEs, SPDEs and their Applications July 4, 2017 Introduction 2/29 Introduction
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 informationContinous time models and realized variance: Simulations
Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU CBOE Conference on Derivatives and Volatility, Chicago, Nov. 10, 2017 Peter Carr (NYU) Volatility Smiles and Yield Frowns 11/10/2017 1 / 33 Interest Rates
More informationConvergence of trust-region methods based on probabilistic models
Convergence of trust-region methods based on probabilistic models A. S. Bandeira K. Scheinberg L. N. Vicente October 24, 2013 Abstract In this paper we consider the use of probabilistic or random models
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 informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
More informationOn Using Shadow Prices in Portfolio optimization with Transaction Costs
On Using Shadow Prices in Portfolio optimization with Transaction Costs Johannes Muhle-Karbe Universität Wien Joint work with Jan Kallsen Universidad de Murcia 12.03.2010 Outline The Merton problem The
More information