A multifactor, stochastic volatility HJM model in a low dimensional markov representation: theory overview and implementation details
|
|
- Dwayne Peters
- 6 years ago
- Views:
Transcription
1 A multifactor, stochastic volatility HJM model in a low dimensional markov representation: theory overview and implementation details Dr. Michael Dirkmann
2 A stochastic volatility extended Cheyette Model: theory overview and implementation details Dr. Michael Dirkmann
3 YASM Yet Another Shortrate Model Dr. Michael Dirkmann
4 Contents 1. Relationship between Shortrate Model and HJM 2. 1 factor Cheyette Model 3. Swaptions in the 1 factor Cheyette Model 4. 1 factor Stochastic Volatility Cheyette Model 5. Swaptions in the 1 factor Stochastic Volatility Cheyette Model 6. Calibration issues 7. Simulation issues 8. PDE issues 9. Multifactor SV Cheyette Modell
5 Introduction We have the LMMs why bother with shortrate models anyway? looking for a slim alternative to LMM try to capture as many features as possible with only few state variables few state variables will be faster in MC simulation few state variables will provide chance for PDE solution So give it a try...
6 setting and integrating yields df(t, T ) = σ f (t, T ) ( T f(t, T ) = f(0, t) + g(t ) g(t) HJM (I) t σ f (t, s)ds ) σ f (t, T ) = h(t)g(t ) ( x(t) + y(t) 1 g(t) dt + σ f (t, T )dw (t) T t g(s)ds ) with dx(t) = dy(t) = ( g ) (t) x(t) + y(t) dt + g(t)h(t)dw (t) g(t) ( ) g 2 (t)h 2 (t) + 2 g (t) g(t) y(t) dt 1/33
7 HJM (II) dx(t) = dy(t) = ( g ) (t) x(t) + y(t) dt + g(t)h(t)dw (t) g(t) ( ) g 2 (t)h 2 (t) + 2 g (t) g(t) y(t) dt with g (t) = k(t) g(t) g(t)h(t) = η ( t, x(t), y(t) ) yields dx(t) = ( y(t) k(t)x(t) ) dt + η (t, x, (t), y(t)) dw (t) dy(t) = ( η 2 (t, x(t), y(t)) 2k(t)y(t) ) dt 2/33
8 See e.g. [1], [2] Cheyette 1f (I) dx(t) = ( y(t) k(t)x(t) ) dt + η(t, x(t), y(t))dw (t) dy(t) = ( η 2 (t, x(t), y(t)) 2k(t)y(t) ) dt r(t) = x(t) + f(0, t) Zero Coupon Bond: P ( t, T ; x(t), y(t) ) = D(T ) ( D(t) exp G(t, T )x(t) + 1 ) 2 G2 (t, T )y(t) G(t, T ) = T t e R u t For any choice of η ( t, x(t), y(t) )! So, what is the best choice of η ( t, x(t), y(t) )? k(s)ds du = 1 k k(t (1 e t)) 3/33
9 Cheyette 1f (II) Displaced diffusion (DD) η(t, x) = ( m r(x, t) + (1 m)l ) σ(t) η(t, x, y) = ( m S nn (x, y, t) + (1 m)l ) σ(t) Constant Elasticity of Variance (CEV) η(t, x) = r(x, t) α σ(t) ; η(t, x) = S nn (x, y, t) α σ(t) DD with stochastic volatility (SV) η(t, x) = ( m r(t) + (1 m)l ) V (t)σ(t) dv = β(1 V (t))dt + V (t)ɛdw V ; V (0) = 1 4/33
10 Cheyette 1f (III) Displaced Diffusion: η(t, x) = ( m r(t) + (1 m)l ) σ(t) with m = 0 η(t, x) = ( mr(t) + (1 m)l ) σ(t) = Lσ(t) dx(t) = ( y(t) kx(t) ) dt + Lσ(t)dW (t) dy(t) = ( L 2 σ 2 (t) 2ky(t) ) dt yields extended Vasicek model t y(t) = e 2kt L 2 σ 2 (s)e 2ks ds 0 dx(t) = ( y(t) kx(t) ) dt + Lσ(t)dW (t) 5/33
11 Cheyette 1f (IV) Displaced Diffusion: η(t, x) = ( m r(t) + (1 m)l ) σ(t) with m = 1 η(x, t) = ( mr(t) + (1 m)l ) σ(t) = r(x, t)σ(t) dx(t) = ( y(t) kx(t) ) dt + η(x, t)dw (t) dy(t) = ( η 2 (x, t) 2ky(t) ) dt = ( x 2 (t)σ 2 (t) 2ky(t) ) dt we get a kind of lognormal model dx(t) = ( y(t) kx(t) ) dt + r(x, t)σ(t)dw (t) dy(t) = ( η 2 (x, t) 2ky(t) ) dt = ( σ 2 (t)x 2 (t) 2ky(t) ) dt 6/33
12 What do we have? Intermediate Summary (I) shortrate model with only few state variables variable volatility specification (including stochastic volatility) analytic function for zero coupon bond: P ( t, T, x(t), y(t) ) What do we need? fast calculation of instruments to calibrate to (e.g. swaptions) simulation of process scheme for PDE solving 7/33
13 Cheyette 1f Swaptions (I) Displaced Diffusion: η(t, x) = ( mr(t) + (1 m)l ) σ(t) Swaption in swap-measure (annuity measure) [3]: S nn (t) = P (t,t n) P (t,t N ) P N1 α i P (t,t i ) ds nn = S nn x dx = S ( ) nn mr(t) + (1 m)l σ(t)dw x [ ] SnN (mr(t) + (1 m)l) ( ) ms nn + (1 m)l σ(t)dw x (ms nn (t) + (1 m)l) ds nn ( ms nn + (1 m)l ) λ(t)dw x=y=0 Displaced Diffusion (lognormal in S): SWP = P V 01 E [ (S K) +] = P V [ 01 m E ( S K) +] = P V 01 BS( S, K) m with S = ms + (1 m)l ; K = mk + (1 m)l 8/33
14 Comments on approximation [3]: Cheyette 1f Swaptions (II) S nn (t) = P (t, t n) P (t, t N ) P N1 α i P (t, t i ) ; P (t, T ) = D(T ) D(t) exp G(t, T )x + 1 «2 G2 (t, T )y Errors in approximation freezing S x [ ] S x x=y=0 since E [ ] [ [ S ] S x E x removes variance ] expectation is biased x=y=0 The above effects are small and partially offsetting. (Unfortunately true only for m 0.) no drift correction for x(t) in annuity meassure residual skew and biased expectation 9/33
15 Cheyette 1f Swaptions (III) impl.vol log (m = 1) norm (m = 0) rel.strike ds nn = ( ms nn + (1 m)l ) λ(t)dw 10/33
16 Cheyette 1f Swaptions (IV) impl.vol approximation simulation normal (m = 0) 10y into 10y rel.strike 11/33
17 Cheyette 1f Swaptions (V) impl.vol lognormal (m = 1) 10y into 10y approximation simulation rel.strike 12/33
18 Cheyette 1f stochastic volatility (SV) SDE (see e.g. [2]): dx(t) = ( y(t) kx(t) ) dt + η(x, t) V (t)dw x (t) dy(t) = ( η 2 (x, t) 2ky(t) ) dt η(x, t) = ( m(t)x(t) + (1 m(t))l ) σ(t) dv (t) = β ( 1 V (t) ) dt + ɛ(t) V (t)dw V (t) dw x dw V = 0 Zero Coupon Bond: like in Cheyette without SV (Unspanned Stochastic Volatility) Zero Coupon Bond: P (t, T ; x, y) = D(T ) D(t) exp G(t, T ) = 1 k k(t (1 e t)) ( G(t, T )x(t) G2 (t, T )y(t) ) 13/33
19 Cheyette 1f SV Swaption (I) Start with time constant displacement (m): η(x, t) = ( mr(t) + (1 m)l ) σ(t) dv (t) = β ( 1 V (t) ) dt + ɛ(t) V (t)dw V (t) change of measure into swap-measure (like in non-sv case) due to 0-correlation between dw x and dw V no change of drift in dv Swaption in swap-measure (annuity measure) plus approximations: ( ) ds nn = ms nn + (1 m)l λ(t) V (t)dw dv (t) = β ( 1 V (t) ) dt + ɛ V (t)dw V (t) well known Heston SDE (+ Displaced Diffusion) 14/33
20 Cheyette 1f SV Swaption (II) ds nn = ( ) ms nn + (1 m)l λ(t) V (t)dw dv (t) = β ( 1 V (t) ) dt + ɛ(t) V (t)dw V (t) SWP = P V 01 E [ (S K) +] = P V [ 01 m E ( S K) +] = P V 01 Heston( S, K) m with S = ms + (1 m)l ; K = mk + (1 m)l The Heston part can be solved by fundamental transform (see below) volatility expansion 15/33
21 Heston( S, K) = S(0) + K 2π R da dt db dt Cheyette 1f SV Swaption (III) k i =0.5 = S(0) + K π iki + ik i S(0) K = βb ; a(t, ω) = 0 0 S(0) iω ln e K ω 2 iω ea(0,ω)+b(0,ω) dω cos ( k ln S(0) ) K k 2 e a(0,k+0.5i)+b(0,k+0.5i) dk = βb 1 2 ɛ2 (t)b λ2 (t)m 2 k 2 ; b(t, ω) = 0 Riccati ordinary differential equation has analytical solution for ɛ and λ constant, piecewise analytical solution for ɛ and λ piecewise constant and numerical solution for ɛ(t) and λ(t) (which is slower) Unfortunately λ(t) is not piecewise constant; (whereas σ(t) is) 16/33
22 Cheyette 1f SV Swaption (IV) can solve model with λ(t) and ɛ(t) and m cannot use lognormal property of Displaced Diffusion anymore when m becomes m(t) ds nn = ( m(t)s nn + (1 m(t))l ) λ(t) V (t)dw S dv (t) = β ( 1 V (t) ) dt + ɛ(t) V (t)dw V find time averaged variables to approximate dynamics: m for solvebility ; λ and ɛ for speed up d S nn = ( m S nn + (1 m)l ) λ V (t)dw S d V (t) = β ( 1 V (t) ) dt + ɛ V (t)dw V 17/33
23 Cheyette 1f SV Swaption (V) Average volatility of variance ɛ (V.Piterbarg[5]): Equate variance of realized volatility [( T E 0 λ 2 (t) V ) 2 ] (t)dt [( T = E 0 ) 2 ] λ 2 (t)v (t)dt yielding with w(t) = T t T s ɛ 2 = T 0 T 0 ɛ2 (t)w(t)dt w(t)dt λ 2 (u)λ 2 (s)e β(u s) e 2β(s t) du ds 18/33
24 Cheyette 1f SV Swaption (VI) impl.vol[%] y into 10y averaging: ε slow sol. to ODE: ε(t) rel.strike L = 4% k = 2.4% m = 0.2 σ = 18.5% ɛ(0) = 80%, ɛ(10) = 20% 19/33
25 Cheyette 1f SV Swaption (VII) Effective skew parameter m (V.Piterbarg[5]): q ds nn = η (S(t), t) λ(t) V (t)dw S with η (S(t), t) = m(t)s nn + (1 m(t))l q dv (t) = β (1 V (t)) dt + ɛ(t) V (t)dw V V.P. shows that η 2 (S(t)) = T 0 w(t)η(s(t), t)λ2 (t)dt T 0 w(t)λ2 (t)dt with w(t) = E [V (t) ( S(t) S(0) ) ] 2 minimizes the difference of the second and third moment between the time averaged and timedependent SDE. For the given η (S(t), t) this yields m = T 0 T 0 m(t)w(t)λ2 (t)dt w(t)λ2 (t)dt with w(t) = t 0 λ 2 (s)ds+ ɛ 2 e βt t 0 λ 2 (s)( e βs e βs) ds 2β 20/33
26 Cheyette 1f SV Swaption (VIII) impl.vol[%] swaption (20y into 10y) approximation (avg. m) simulation (annuity measure) rel.strike L = 4% k = 2.4% m(0) = 0.8, m(15) = 0.3 σ = 18.5% ɛ = 60% 21/33
27 Cheyette 1f SV Swaption (IX) Effective volatility λ (V.Piterbarg[5]): step 1: Black-Scholes ATM-Swaption formula as function of variance is approximated with by setting g(x) = S nn m ( 2N( 1 ) 2 m x) 1 g(x) = a + be cx. g(ξ) = g(ξ) ; g (ξ) = g (ξ) and g (ξ) = g (ξ) step 2: Equate expectations for time dependent and time averaged λ a + be E [ g ( T 0 λ 2 (t)v (t) )] [ = E g ( λ2 T [e c R T 0 λ2 (t)v (t)dt ] = a + be V (t) )] 0 [e c λ 2 R T 0 V (t)dt] 22/33
28 Cheyette 1f SV Swaption (X) still step 2: a + be [e c R ] T 0 λ2 (t)v (t)dt = a + be [e c λ 2 R T 0 V (t)dt] E [e c R ] T 0 λ2 (t)v (t)dt = E [e c λ 2 R T 0 V (t)dt] e A(0,T )+B(0,T ) = eā(0,t )+ B(0,T ) a() and a() follow Riccati ODE. For ā() and ā() anayltic solutions exist. time dependent part: once solve for a() and b() numerically time independent part: solve for λ numerically (cheap since analytical solutions to ā() and b() exist 23/33
29 Cheyette 1f SV Swaption (XI) Some empirical observations for swaptions in the Euro market in June: the skew parameter m tends to be well below 1 and not strongly varying the smile parameter ɛ is constant around 50% for larger t, a bit larger (60%) at small t (t < 2years) the volatiliy σ(t) is about 18% and not strongly fluctuating 24/33
30 Cheyette 1f SV Swaption (XIIa) impl.vol market (10y into 10y) approximation simulation simulation tuned rel.strike L = 4% k = 2.4% m 0.2[0.1] σ 18% ɛ = 60%[50%] 25/33
31 impl.vol.diff.[%] Cheyette 1f SV Swaption (XIIb) diff to market (10y into 10y) diff to approximation diff to simulation diff to tuned simulation rel.strike L = 4% k = 2.4% m 0.2[0.1] σ 18% ɛ = 60%[50%] 26/33
32 Intermediate Summary (II) by transformation into swap measure freezing of ds nn dx η(r) [ dsnn dx η(r) η(s nn ) ] η(s nn ) x=y=0 using fact that displaced diffusion leads to lognormal distribution in substitute variable time averaging of σ(t), λ(t) and ɛ(t) we get more (m small) or less (m large) accurate numerical solution can (largely) remove residual skew of approximation by correting for deviation between approximation and simulation once (denoted by tuned simulation) hence we are prepared for calibration of the model 27/33
33 Calibration to Swaptions minimize sum of weighted squared relative deviations of model price from market price global: for all model parameters at once : slow and unstable stepwise: assume knowledge of time independent parameters k and β; fit time dependent parameters m(t), σ(t) and ɛ(t) for subsequent swaption expiries iterative: optimize time independent parameters combined with a stepwise calibration for each valuation of the cost function depending of the number of factors, time intervals and market quotes this takes from fractions of a second to several minutes for very accurate calibration one intermediate simulation step is helpful (removing residual skew in swaption price approximation for m > 0) 28/33
34 Simulation of Variance Process (CIR) dv (t) = β ( 1 V (t) ) dt + ɛ V (t)dw V (t) V (0) = 1 V (t + t) = V (t) + β ( 1 V (t) ) t + ɛ V (t) tz process itself cannot yield V (t) < 0 for large β and / or large ɛ : Euler discretization can yield negative V (t) if negative, set to 0 : does not work at all if negative, mirror to postive value : does not work well sample from non-central χ 2 [6]: works well, no antithetics (fast implementation in GSL (not inverting distribution function)) moment matched log-normal[7] : works well, antithetics possible V (t + t) = ( 1 + (V (t) 1)e β ) e 0.5Γ2 (t)+γ(t)z Γ 2 (t) = ln ( 1 + ɛ 2 V (t)(1 e 2β t ) ) 2β(1 + (V (t) 1)e β t ) 2 29/33
35 PDE solution (I) PDE for 1factor without stochastic volatility V t = (D x + D y r)v ( ) 2 2 mr(t) (1 m)l x 2 D x = (y kx) x ( ( ) ) 2 D y = mr(t) (1 m)l 2ky y two dimensional PDE; y-dimension is convection only! ADI Craig Sneyd (splitting scheme) IMEX (see below) 30/33
36 PDE solution (II) IMEX-scheme [8]: V (t + ) V (t) t = µ y y (V (y+, t) V (t)) h µy alt. y 1 6 V (y, t) + V (y, t) 1 2 V (t) 1 «i 3 V (y+, t) 1 V (x +, t) + V (x, t) 2V (t) V (x +, t) V (x, t) 2 σ2 x x 2 + µ x 2 x rv (t) 1 V (x +, t + ) + V (x, t + ) 2V (t + ) V (x +, t + ) V (x, t + ) 2 σ2 x x 2 + µ x 2 x! rv (t + )! purely explicit in y use upwind scheme in y, not central differences : oscillations use third order upwind in y to increase accuracy to O( y 3 ): no penalty in solving system of linear equations since explicit can get accurate solution of swaption with 15 points in y (number of necesarry gridpoints in y depends on m; for m 0 we can get away with even less than 15) 31/33
37 Multifactor One multidimensional extension of the model is straight forward dx i (t) = ( y ij (t) k i x i (t) ) f dt + η ij (t)dw j (t) ; dw i dw j = 0 j=1 dy ij (t) = ( f η ik (t)η jk (t) (k i + k j )y ij (t) ) dt k=1 η ij (t, x i, y ij ) = σ ij (t) ( m(t)r(t) + (1 m(t))l ) σ ij (t) = lower triangular matrix r(t) = i x i (t) + f(0, t) f stochastic and f 3+f 2 total number of state variables analytic formula for P(t, x, ȳ) similar approximations for swaption 32/33
38 Summary / Outlook What we have: class of shortrate models with flexible volatility function flexible smile and / or skew stochastic volatility with corresponding dynamics low dimensional Markow and analytical solution for ZCB: relatively cheap in MC-simulation and PDE 1 factor version is capable of reproducing the smile for one underlying per expiry nicely (e.g coterminal swaptions); even swaptions not used in calibration process often still ok; further improvements by multi factor version expected Future plans: improve swaption approximation in annuity measure for m 0 extensions and checks of multi factor model include jump diffusion 33/33
39 References [1] Markov Representation of the HJM Model; O. Cheyette; Working paper: Barra [2] Volatile Volatilities; L. Andersen, J.Andreasen; Risk December 2002 [3] Pricing Swaptions and Coupon Bond Options in Affine TS Models D. Schrager, A. Pelsser [4] Option Valuation under Stochastic Volatility; A. Lewis [5] Time to smile and Stochastic Volatility Model with Timedependent Skew; V. Piterbarg; Risk May 2005 [6] Exact Simulation of SV and other Affine Jump Diffusion Processes M. Broadie, Ö. Kaya [7] Extended Libor Market Models with Stochastic Volatility L. Andersen, R. Brotherton-Ratcliffe
40 [8] Numerical Solution of Time-Dep. Advection-Diffusion-Reaction Eqs. W. Hundsdorfer, J. Verwer 35/33
Risk 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 informationLecture 5: Review of interest rate models
Lecture 5: Review of interest rate models Xiaoguang Wang STAT 598W January 30th, 2014 (STAT 598W) Lecture 5 1 / 46 Outline 1 Bonds and Interest Rates 2 Short Rate Models 3 Forward Rate Models 4 LIBOR and
More informationLIBOR models, multi-curve extensions, and the pricing of callable structured derivatives
Weierstrass Institute for Applied Analysis and Stochastics LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives John Schoenmakers 9th Summer School in Mathematical Finance
More informationInterest Rate Volatility
Interest Rate Volatility III. Working with SABR Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline Arbitrage free SABR 1 Arbitrage free
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 informationInterest rate models and Solvency II
www.nr.no Outline Desired properties of interest rate models in a Solvency II setting. A review of three well-known interest rate models A real example from a Norwegian insurance company 2 Interest rate
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 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 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 informationMarket interest-rate models
Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations
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 informationEuropean call option with inflation-linked strike
Mathematical Statistics Stockholm University European call option with inflation-linked strike Ola Hammarlid Research Report 2010:2 ISSN 1650-0377 Postal address: Mathematical Statistics Dept. of Mathematics
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 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 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 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 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 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 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 informationMarkovian Projection, Heston Model and Pricing of European Basket Optio
Markovian Projection, Heston Model and Pricing of European Basket Options with Smile July 7, 2009 European Options on geometric baskets European Options on arithmetic baskets Pricing via Moment Matching
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 informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
More informationMulti-dimensional Term Structure Models
Multi-dimensional Term Structure Models We will focus on the affine class. But first some motivation. A generic one-dimensional model for zero-coupon yields, y(t; τ), looks like this dy(t; τ) =... dt +
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 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 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 informationOne-Factor Models { 1 Key features of one-factor (equilibrium) models: { All bond prices are a function of a single state variable, the short rate. {
Fixed Income Analysis Term-Structure Models in Continuous Time Multi-factor equilibrium models (general theory) The Brennan and Schwartz model Exponential-ane models Jesper Lund April 14, 1998 1 Outline
More informationarxiv: v1 [q-fin.pr] 23 Feb 2014
Time-dependent Heston model. G. S. Vasilev, Department of Physics, Sofia University, James Bourchier 5 blvd, 64 Sofia, Bulgaria CloudRisk Ltd (Dated: February 5, 04) This work presents an exact solution
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 informationUnified Credit-Equity Modeling
Unified Credit-Equity Modeling Rafael Mendoza-Arriaga Based on joint research with: Vadim Linetsky and Peter Carr The University of Texas at Austin McCombs School of Business (IROM) Recent Advancements
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 informationAN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL
AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it 1 Stylized facts Traders use the Black-Scholes formula to price plain-vanilla options. An
More informationHeston Model Version 1.0.9
Heston Model Version 1.0.9 1 Introduction This plug-in implements the Heston model. Once installed the plug-in offers the possibility of using two new processes, the Heston process and the Heston time
More informationCredit Risk using Time Changed Brownian Motions
Credit Risk using Time Changed Brownian Motions Tom Hurd Mathematics and Statistics McMaster University Joint work with Alexey Kuznetsov (New Brunswick) and Zhuowei Zhou (Mac) 2nd Princeton Credit Conference
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 informationMultiscale Stochastic Volatility Models Heston 1.5
Multiscale Stochastic Volatility Models Heston 1.5 Jean-Pierre Fouque Department of Statistics & Applied Probability University of California Santa Barbara Modeling and Managing Financial Risks Paris,
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 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 information1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:
1 Mathematics in a Pill The purpose of this chapter is to give a brief outline of the probability theory underlying the mathematics inside the book, and to introduce necessary notation and conventions
More informationOptimal asset allocation under forward performance criteria Oberwolfach, February 2007
Optimal asset allocation under forward performance criteria Oberwolfach, February 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 References Indifference valuation in binomial models (with
More informationOPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF FINITE
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 WeA11.6 OPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF
More informationRough volatility models: When population processes become a new tool for trading and risk management
Rough volatility models: When population processes become a new tool for trading and risk management Omar El Euch and Mathieu Rosenbaum École Polytechnique 4 October 2017 Omar El Euch and Mathieu Rosenbaum
More informationGeneralized Affine Transform Formulae and Exact Simulation of the WMSV Model
On of Affine Processes on S + d Generalized Affine and Exact Simulation of the WMSV Model Department of Mathematical Science, KAIST, Republic of Korea 2012 SIAM Financial Math and Engineering joint work
More information(1) Consider a European call option and a European put option on a nondividend-paying stock. You are given:
(1) Consider a European call option and a European put option on a nondividend-paying stock. You are given: (i) The current price of the stock is $60. (ii) The call option currently sells for $0.15 more
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 informationLecture 11: Ito Calculus. Tuesday, October 23, 12
Lecture 11: Ito Calculus Continuous time models We start with the model from Chapter 3 log S j log S j 1 = µ t + p tz j Sum it over j: log S N log S 0 = NX µ t + NX p tzj j=1 j=1 Can we take the limit
More informationFINANCIAL PRICING MODELS
Page 1-22 like equions FINANCIAL PRICING MODELS 20 de Setembro de 2013 PhD Page 1- Student 22 Contents Page 2-22 1 2 3 4 5 PhD Page 2- Student 22 Page 3-22 In 1973, Fischer Black and Myron Scholes presented
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 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 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 informationCrashcourse Interest Rate Models
Crashcourse Interest Rate Models Stefan Gerhold August 30, 2006 Interest Rate Models Model the evolution of the yield curve Can be used for forecasting the future yield curve or for pricing interest rate
More informationPortability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans
Portability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans An Chen University of Ulm joint with Filip Uzelac (University of Bonn) Seminar at SWUFE,
More informationManaging Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives
Managing Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives Simon Man Chung Fung, Katja Ignatieva and Michael Sherris School of Risk & Actuarial Studies University of
More informationMarket models for the smile Local volatility, local-stochastic volatility
Market models for the smile Local volatility, local-stochastic volatility Lorenzo Bergomi lorenzo.bergomi@sgcib.com Global Markets Quantitative Research European Summer School in Financial Mathematics
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 informationLinear-Rational Term-Structure Models
Linear-Rational Term-Structure Models Anders Trolle (joint with Damir Filipović and Martin Larsson) Ecole Polytechnique Fédérale de Lausanne Swiss Finance Institute AMaMeF and Swissquote Conference, September
More informationMarkovian projection to a Displaced Volatility Heston model
Markovian projection to a Displaced Volatility Heston model A. Antonov, M. Arneguy and N. Audet March 12, 2008 Abstract Markovian Projection is an optimal approximation of a complex underlying process
More information18. Diffusion processes for stocks and interest rates. MA6622, Ernesto Mordecki, CityU, HK, References for this Lecture:
18. Diffusion processes for stocks and interest rates MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: P. Willmot, Paul Willmot on Quantitative Finance. Volume 1, Wiley, (2000) A.
More informationHedging under Model Uncertainty
Hedging under Model Uncertainty Efficient Computation of the Hedging Error using the POD 6th World Congress of the Bachelier Finance Society June, 24th 2010 M. Monoyios, T. Schröter, Oxford University
More informationLOGNORMAL MIXTURE SMILE CONSISTENT OPTION PRICING
LOGNORMAL MIXTURE SMILE CONSISTENT OPTION PRICING FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it Daiwa International Workshop on Financial Engineering, Tokyo, 26-27 August 2004 1 Stylized
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 informationThings You Have To Have Heard About (In Double-Quick Time) The LIBOR market model: Björk 27. Swaption pricing too.
Things You Have To Have Heard About (In Double-Quick Time) LIBORs, floating rate bonds, swaps.: Björk 22.3 Caps: Björk 26.8. Fun with caps. The LIBOR market model: Björk 27. Swaption pricing too. 1 Simple
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Haindorf, 7 Februar 2008 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar
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 informationA Brief Introduction to Stochastic Volatility Modeling
A Brief Introduction to Stochastic Volatility Modeling Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction When using the Black-Scholes-Merton model to
More informationLibor Market Model Version 1.0
Libor Market Model Version.0 Introduction This plug-in implements the Libor Market Model (also know as BGM Model, from the authors Brace Gatarek Musiela). For a general reference on this model see [, [2
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 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 informationStochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models
Stochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models Eni Musta Università degli studi di Pisa San Miniato - 16 September 2016 Overview 1 Self-financing portfolio 2 Complete
More informationOptimal investments under dynamic performance critria. Lecture IV
Optimal investments under dynamic performance critria Lecture IV 1 Utility-based measurement of performance 2 Deterministic environment Utility traits u(x, t) : x wealth and t time Monotonicity u x (x,
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 informationLocal Volatility Pricing Models for Long-Dated FX Derivatives
Local Volatility Pricing Models for Long-Dated FX Derivatives G. Deelstra, G. Rayee Université Libre de Bruxelles grayee@ulb.ac.be Gregory Rayee (ULB) Bachelier Congress 2010 Toronto, 22-26 june 2010 1
More informationθ(t ) = T f(0, T ) + σ2 T
1 Derivatives Pricing and Financial Modelling Andrew Cairns: room M3.08 E-mail: A.Cairns@ma.hw.ac.uk Tutorial 10 1. (Ho-Lee) Let X(T ) = T 0 W t dt. (a) What is the distribution of X(T )? (b) Find E[exp(
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 informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
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 informationExploring Volatility Derivatives: New Advances in Modelling. Bruno Dupire Bloomberg L.P. NY
Exploring Volatility Derivatives: New Advances in Modelling Bruno Dupire Bloomberg L.P. NY bdupire@bloomberg.net Global Derivatives 2005, Paris May 25, 2005 1. Volatility Products Historical Volatility
More informationA new approach to LIBOR modeling
A new approach to LIBOR modeling Antonis Papapantoleon FAM TU Vienna Based on joint work with Martin Keller-Ressel and Josef Teichmann Istanbul Workshop on Mathematical Finance Istanbul, Turkey, 18 May
More informationLecture Notes for 5765/6895, Part II
Lecture Notes for 5765/6895, Part II The choices of material we cover after Chapter 5 have been more flexible, reflecting recent developments and our own interest. In the following, we discuss and summarize
More informationZABR -- Expansions for the Masses
ZABR -- Expansions for the Masses Preliminary Version December 011 Jesper Andreasen and Brian Huge Danse Marets, Copenhagen want.daddy@danseban.com brno@danseban.com 1 Electronic copy available at: http://ssrn.com/abstract=198076
More informationOption Pricing for a Stochastic-Volatility Jump-Diffusion Model
Option Pricing for a Stochastic-Volatility Jump-Diffusion Model Guoqing Yan and Floyd B. Hanson Department of Mathematics, Statistics, and Computer Science University of Illinois at Chicago Conference
More informationForward Dynamic Utility
Forward Dynamic Utility El Karoui Nicole & M RAD Mohamed UnivParis VI / École Polytechnique,CMAP elkaroui@cmapx.polytechnique.fr with the financial support of the "Fondation du Risque" and the Fédération
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 informationMarket Price of Longevity Risk for A Multi-Cohort Mortality Model with Application to Longevity Bond Option Pricing
1/51 Market Price of Longevity Risk for A Multi-Cohort Mortality Model with Application to Longevity Bond Option Pricing Yajing Xu, Michael Sherris and Jonathan Ziveyi School of Risk & Actuarial Studies,
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 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 informationFourier Space Time-stepping Method for Option Pricing with Lévy Processes
FST method Extensions Indifference pricing Fourier Space Time-stepping Method for Option Pricing with Lévy Processes Vladimir Surkov University of Toronto Computational Methods in Finance Conference University
More informationReading: You should read Hull chapter 12 and perhaps the very first part of chapter 13.
FIN-40008 FINANCIAL INSTRUMENTS SPRING 2008 Asset Price Dynamics Introduction These notes give assumptions of asset price returns that are derived from the efficient markets hypothesis. Although a hypothesis,
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 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 informationdt+ ρσ 2 1 ρ2 σ 2 κ i and that A is a rather lengthy expression that we may or may not need. (Brigo & Mercurio Lemma Thm , p. 135.
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where ( κ1 0 dx(t) = 0 κ 2 r(t) = δ 0 +X 1 (t)+x 2 (t) )( X1 (t) X 2 (t) ) ( σ1 0 dt+ ρσ 2 1 ρ2 σ 2 )( dw Q 1 (t) dw Q 2 (t) ) In this
More information4. Black-Scholes Models and PDEs. Math6911 S08, HM Zhu
4. Black-Scholes Models and PDEs Math6911 S08, HM Zhu References 1. Chapter 13, J. Hull. Section.6, P. Brandimarte Outline Derivation of Black-Scholes equation Black-Scholes models for options Implied
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 informationdt + ρσ 2 1 ρ2 σ 2 B i (τ) = 1 e κ iτ κ i
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where dx(t) = ( κ1 0 0 κ 2 ) ( X1 (t) X 2 (t) In this case we find (BLACKBOARD) that r(t) = δ 0 + X 1 (t) + X 2 (t) ) ( σ1 0 dt + ρσ 2
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU IFS, Chengdu, China, July 30, 2018 Peter Carr (NYU) Volatility Smiles and Yield Frowns 7/30/2018 1 / 35 Interest Rates and Volatility Practitioners and
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 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 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 informationKØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours
This question paper consists of 3 printed pages FinKont KØBENHAVNS UNIVERSITET (Blok 2, 211/212) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours This exam paper
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More information