Extended Libor Models and Their Calibration
|
|
- Ezra Barnett
- 5 years ago
- Views:
Transcription
1 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 / 43
2 Overview 1 Introduction Forward Libor Models 2 Modelling Modelling under Terminal Measure Modelling under Forward Measures 3 Pricing and Calibration Pricing of Caplets Specification Analysis Calibration Procedure 4 Calibration in work Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
3 Forward Libor Models Introduction Forward Libor Models Tenor structure: 0 = T 0 < T 1 <... < T M < T M+1 with accrual periods δ i := T i+1 T i Zero coupon bonds: B k (t), t [0, T k ] with B k (T k ) = 1 Forward Libor rates: L 1 (t),..., L M (t) L k (t) = 1 ( ) Bk (t) δ k B k+1 (t) 1, t [0, T k ], k = 1,..., M Remark L 1,..., L M are based on simple compounding that is an investor receives 1$ at T k and pays 1 + δ k L k (t) at T k+1 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
4 Forward Libor Models Introduction Forward Libor Models Tenor structure: 0 = T 0 < T 1 <... < T M < T M+1 with accrual periods δ i := T i+1 T i Zero coupon bonds: B k (t), t [0, T k ] with B k (T k ) = 1 Forward Libor rates: L 1 (t),..., L M (t) L k (t) = 1 ( ) Bk (t) δ k B k+1 (t) 1, t [0, T k ], k = 1,..., M Remark L 1,..., L M are based on simple compounding that is an investor receives 1$ at T k and pays 1 + δ k L k (t) at T k+1 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
5 Modelling Modelling under Terminal Measure The BGM/Jamshidian LIBOR market model dl i = M j=i+1 δ j L i L j γ i γ j 1 + δ j L j dt + L i γ i dw (M+1), where (W (M+1) (t) 0 t T M ) is a D-dimensional Wiener process under P M+1 and γ i = (γ i,1,..., γ i,d ) R D, i = 1,..., M are deterministic volatility vector functions, called factor loadings, to be determined via calibration. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
6 Modelling Modelling under Terminal Measure Scalar volatility and correlation functions A scalar volatility functions t σ i (t) := γ i (t) = d γi,k 2 (t), 0 t T i, 1 i M, k=1 and a (local) correlation functions t ϱ ij (t) := γ i(t) γ j (t) γ i (t) γ j (t), 0 t min(t i, T j ), 1 i, j M, are the only economically relevant objects! Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
7 Modelling Modelling under Terminal Measure Scalar volatility and correlation functions For any deterministic orthogonal matrix valued map Q : t Q(t) R D D the process W t = t 0 Q(u)dW u is again a standard Brownian motion in R D. Hence, taking dw t := Q d W t, yields the same Libor process with γ := γq. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
8 Modelling The Hull-White Parametrization Modelling under Terminal Measure For a fixed correlation matrix ϱ (0) := [ϱ (0) ij ] 0 i,j<m of rank D and a fixed non-negative vector Λ := [Λ i ] 0 i<m take σ i (t) := Λ i m(t), 0 t T i, ϱ ij (t) := ϱ (0) i m(t),j m(t), 0 t min(t i, T j ), 1 i, j M with m(t) := min{m : T m t}. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
9 Modelling Modelling under Terminal Measure The Hull-White Parametrization Find unit vectors e (0) i R D satisfying ϱ (0) ij = e (0) i e (0) j and set e i (t) := e (0) i m(t). The dynamics of the Hull-White LIBOR model is determined by the volatility structure γ i (t) := σ i (t)e i (t) = Λ i m(t) e (0) i m(t), 1 i M. Observation Good feature: A kind of time-shift homogeneity which is economically sensible. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
10 Modelling Modelling under Terminal Measure The Hull-White Parametrization Find unit vectors e (0) i R D satisfying ϱ (0) ij = e (0) i e (0) j and set e i (t) := e (0) i m(t). The dynamics of the Hull-White LIBOR model is determined by the volatility structure γ i (t) := σ i (t)e i (t) = Λ i m(t) e (0) i m(t), 1 i M. Observation Good feature: A kind of time-shift homogeneity which is economically sensible. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
11 Modelling Modelling under Terminal Measure More flexible volatility structure Define γ i (t) = c i g(t i t)e (0) i m(t), 0 t min(t i, T j ), 1 i, j M, ϱ (0) kl = e (0) k e (0) l, 0 k, l < M with suitably parameterized function g and correlation matrix ϱ (0). Observation For c i c, model is basically time-shift homogeneous. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
12 Modelling Modelling under Terminal Measure More flexible volatility structure Define γ i (t) = c i g(t i t)e (0) i m(t), 0 t min(t i, T j ), 1 i, j M, ϱ (0) kl = e (0) k e (0) l, 0 k, l < M with suitably parameterized function g and correlation matrix ϱ (0). Observation For c i c, model is basically time-shift homogeneous. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
13 Modelling Modelling under Terminal Measure Simple parametrization of scalar volatilities g(x) := g a,b,g (x) := g + (1 g + ax)e bx, a, b, g > 0 A typical shape of g : x g 0.5,0.4,0.6 (x) g (x) x Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
14 Modelling Modelling under Terminal Measure Semi-parametric full rank correlation structures Consider initial correlation structures of the following form ϱ (0) ij :== min(b i, b j ), 1 i, j < M, max(b i, b j ) where the sequences i b i and i b i /b i+1 strictly increasing. are positive and Observation Implicit consequence: D = M, w.l.g. b 1 = 1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
15 Modelling Modelling under Terminal Measure Semi-parametric full rank correlation structures The assumption that the sequence i b i b i+1 is increasing forces that for fixed p is increasing. i ϱ i,i+p = Cor( L i, L i+p ) Example Correlation between a seven and a nine year forward is higher than the correlation between a three and a five year forward. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
16 Modelling Modelling under Terminal Measure Semi-parametric full rank correlation structures The assumption that the sequence i b i b i+1 is increasing forces that for fixed p is increasing. i ϱ i,i+p = Cor( L i, L i+p ) Example Correlation between a seven and a nine year forward is higher than the correlation between a three and a five year forward. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
17 Modelling Modelling under Terminal Measure Semi-parametric full rank correlation structures Question Is min(b i, b j )/ max(b i, b j ) a correlation structure at all? Take i b i increasing with b 0 = 0, b 1 = 1 and set a i := Consider standard normal i.i.d r.v. Z k and set Then, for i j Y i := i a k Z k. k=1 Cov(Y i, Y j ) = and the correlations are ϱ Yi,Y j = b i /b j. i ak 2 = b2 i k=1 b 2 i b 2 i 1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
18 Modelling Modelling under Terminal Measure Semi-parametric full rank correlation structures Question Is min(b i, b j )/ max(b i, b j ) a correlation structure at all? Take i b i increasing with b 0 = 0, b 1 = 1 and set a i := Consider standard normal i.i.d r.v. Z k and set Then, for i j Y i := i a k Z k. k=1 Cov(Y i, Y j ) = and the correlations are ϱ Yi,Y j = b i /b j. i ak 2 = b2 i k=1 b 2 i b 2 i 1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
19 Modelling Modelling under Terminal Measure Semi-parametric full rank correlation structures For every correlation structure b i /b j, i < j, with increasing b i and b i /b i+1 it holds [ M ] b i = exp min(l 1, i 1) l l=2 for a nonnegative numbers i, i 0, 2 i M. In particular, i := (ln b i 1 + ln b i+1 2 ln b i ), 2 i M 1, M := ln b m ln b m 1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
20 Modelling Modelling under Terminal Measure Semi-parametric full rank correlation structures The following representation thus holds ϱ ij = min(b [ i, b j ) max(b i, b j ) = exp l 0, 2 l M. M l=i+1 min(l i, j i ) l ], From this representation we may derive conveniently various low parametric structures consistent with the structure ϱ ij = min(b i, b j ) max(b i, b j ). Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
21 Modelling Modelling under Terminal Measure Example: a two parametric structure The choice 2 =... = M 1 =: α > 0 and M =: β > 0 gives the correlation structure i+j+1 i j (β+α(m ϱ ij = e 2 )), i, j = 1,..., m. To gain parameter stability we set, η := α(m 1)(m 2)/2, ϱ := ϱ 1m = exp[ α(m 1)(m 2)/2 (m 1)β] and get a two parametric structure: ϱ ij = e i j M 1( ln ϱ +η M i j+1 M 2 ), 0 < η < ln ϱ with i, j = 1,..., M. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
22 Modelling Modelling under Terminal Measure General Libor Model under P M+1 For i = 1,..., M dl i (t) L i (t) = A(M+1) i (dt)+ Γ i dw (M+1) (t) }{{} Continuous Part + E ( ψ i (t, u) µ ν (M+1)) (dt, du) } {{ } Jump Part A i are predictable drift processes W (M+1) is a D-dimensional Brownian motion under P M+1 Γ i are predictable D-dimensional volatility processes ω µ(dt, du, ω) is a random point measure on R + E with P M+1 -compensator ν (M+1) (dt, du) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
23 Modelling Modelling under Terminal Measure General Libor Model under P M+1 For i = 1,..., M dl i (t) L i (t) = A(M+1) i (dt)+ Γ i dw (M+1) (t) }{{} Continuous Part + E ( ψ i (t, u) µ ν (M+1)) (dt, du) } {{ } Jump Part A i are predictable drift processes W (M+1) is a D-dimensional Brownian motion under P M+1 Γ i are predictable D-dimensional volatility processes ω µ(dt, du, ω) is a random point measure on R + E with P M+1 -compensator ν (M+1) (dt, du) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
24 Modelling Modelling under Terminal Measure General Libor Model under P M+1 For i = 1,..., M dl i (t) L i (t) = A(M+1) i (dt)+ Γ i dw (M+1) (t) }{{} Continuous Part + E ( ψ i (t, u) µ ν (M+1)) (dt, du) } {{ } Jump Part A i are predictable drift processes W (M+1) is a D-dimensional Brownian motion under P M+1 Γ i are predictable D-dimensional volatility processes ω µ(dt, du, ω) is a random point measure on R + E with P M+1 -compensator ν (M+1) (dt, du) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
25 Modelling Modelling under Terminal Measure General Libor Model under P M+1 For i = 1,..., M dl i (t) L i (t) = A(M+1) i (dt)+ Γ i dw (M+1) (t) }{{} Continuous Part + E ( ψ i (t, u) µ ν (M+1)) (dt, du) } {{ } Jump Part A i are predictable drift processes W (M+1) is a D-dimensional Brownian motion under P M+1 Γ i are predictable D-dimensional volatility processes ω µ(dt, du, ω) is a random point measure on R + E with P M+1 -compensator ν (M+1) (dt, du) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
26 Modelling Modelling under Terminal Measure General Libor Model under P M+1 For i = 1,..., M dl i (t) L i (t) = A(M+1) i (dt)+ Γ i dw (M+1) (t) }{{} Continuous Part + E ( ψ i (t, u) µ ν (M+1)) (dt, du) } {{ } Jump Part A i are predictable drift processes W (M+1) is a D-dimensional Brownian motion under P M+1 Γ i are predictable D-dimensional volatility processes ω µ(dt, du, ω) is a random point measure on R + E with P M+1 -compensator ν (M+1) (dt, du) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
27 Modelling Modelling under Terminal Measure Libor Model under P M+1 We consider random point measures µ of finite activity satisfying E [ ( ψ i (t, u) µ ν (M+1)) Nt ] (dt, du) = d ψ i (s l, u l ), where (s l, u l ) R + E are jumps of µ and l=1 N t is a Poisson process with intensity λ E = R... R }{{} m u l R m is distributed with p 1 (dx 1 )... p m (dx m ) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
28 Modelling Modelling under Terminal Measure Libor Model under P M+1 We consider random point measures µ of finite activity satisfying E [ ( ψ i (t, u) µ ν (M+1)) Nt ] (dt, du) = d ψ i (s l, u l ), where (s l, u l ) R + E are jumps of µ and l=1 N t is a Poisson process with intensity λ E = R... R }{{} m u l R m is distributed with p 1 (dx 1 )... p m (dx m ) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
29 Modelling Modelling under Terminal Measure Libor Model under P M+1 We consider random point measures µ of finite activity satisfying E [ ( ψ i (t, u) µ ν (M+1)) Nt ] (dt, du) = d ψ i (s l, u l ), where (s l, u l ) R + E are jumps of µ and l=1 N t is a Poisson process with intensity λ E = R... R }{{} m u l R m is distributed with p 1 (dx 1 )... p m (dx m ) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
30 Modelling Modelling under Terminal Measure Libor Model under P M+1 We consider random point measures µ of finite activity satisfying E [ ( ψ i (t, u) µ ν (M+1)) Nt ] (dt, du) = d ψ i (s l, u l ), where (s l, u l ) R + E are jumps of µ and l=1 N t is a Poisson process with intensity λ E = R... R }{{} m u l R m is distributed with p 1 (dx 1 )... p m (dx m ) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
31 Modelling Modelling under Terminal Measure Drift term under P M+1 The requirement that L i is a martingale under P (M+1) implies A (M+1) i M δ j L j (dt) = Γ i Γ j dt+ 1 + δ j L j j=i+1 + λ(t)dt ψ i (u, t)p(du) R m M j=i+1 ( 1 + δ ) jl j ψ i (t, u). 1 + δ j L j Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
32 Modelling Modelling under Forward Measures Dynamic of L i under P i+1 Since L i is a martingale under P i+1 dl i = Γ i dw (i+1) + ψ i (t, u)(µ ν (i+1) )(dt, du), L i where dw (i+1) = M j=i+1 E δ j L j 1 + δ j L j Γ i dt + dw (M+1) is a standard Brownian motion under P i+1 and M ( ν (i+1) (dt, du) = ν (M+1) (dt, du) 1 + δ ) jl j ψ j (t, u). 1 + δ j L j j=i+1 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
33 Modelling Modelling under Forward Measures Dynamic of L i under P i+1 The logarithmic version reads as [ Nt ] d ln(l i ) = A (i+1) (dt) + Γ i dw (i+1) + d φ i (s l, u l ) i=1 with φ i = ln(1 + ψ i ) and A (i+1) (dt) = 1 2 Γ i 2 dt ψ i (t, u)ν (i+1) (dt, du) R m Observation For i < M the new compensator ν (i+1) is not deterministic and ln(l i ) is generally not affine under P i+1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
34 Modelling Modelling under Forward Measures Dynamic of L i under P i+1 The logarithmic version reads as [ Nt ] d ln(l i ) = A (i+1) (dt) + Γ i dw (i+1) + d φ i (s l, u l ) i=1 with φ i = ln(1 + ψ i ) and A (i+1) (dt) = 1 2 Γ i 2 dt ψ i (t, u)ν (i+1) (dt, du) R m Observation For i < M the new compensator ν (i+1) is not deterministic and ln(l i ) is generally not affine under P i+1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
35 Caplet Volas Caplets Pricing and Calibration Pricing of Caplets The price of j-th caplet at time zero is given by C j (K ) = δ j B j+1 (0)E Pj+1 [(L j (T j ) K ) + ] Tenors Strikes Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
36 Pricing and Calibration Pricing of Caplets Pricing Caplets under P i+1 In terms of log-forward moneyness v = ln(k /L j (0)) C j (v) := δ j B j+1 (0)L j (0)E Pj+1 [(e X j (t) e v ) + ], with X j (t) = log(l j (t)) log(l j (0)). Define then O j (v) = E Pj+1 [(e X j (t) e v ) + ] (1 e v ) +, F{O j }(z) := R O j (v)e ivz dz = 1 Φ P j+1 (z i; T j ). z(z i) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
37 Pricing and Calibration Pricing of Caplets Pricing Caplets under P i+1 In terms of log-forward moneyness v = ln(k /L j (0)) C j (v) := δ j B j+1 (0)L j (0)E Pj+1 [(e X j (t) e v ) + ], with X j (t) = log(l j (t)) log(l j (0)). Define then O j (v) = E Pj+1 [(e X j (t) e v ) + ] (1 e v ) +, F{O j }(z) := R O j (v)e ivz dz = 1 Φ P j+1 (z i; T j ). z(z i) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
38 Pricing and Calibration Specification Analysis Characteristic Function of L M under P M+1 Since d ln(l M ) = 1 2 Γ i 2 dt + Γ i dw (M+1) (t) + d [ Nt ] φ i (s l, u l ) l=1 and N t, W (M+1) and u l are mutually independent Φ PM+1 (z; T ) = Φ C P M+1 (z; T )Φ J P M+1 (z; T ), where Φ C P M+1 (z; T ) ( Φ J P M+1 (z; T )) is the c.f. of continuous (jump) part. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
39 Pricing and Calibration Specification Analysis Specification Analysis: Continuous Part For some predictable vector volatility process (v 1 (t),..., v d (t)) define Γ i = 1 r 2 SV γ i1 1 r 2 SV γ i2.. 1 r 2 SV γ id r SV β i1 v1 (t).. r SV β id vd (t), W (M+1) = W (M+1) 1 W (M+1) 2.. W (M+1) d W (M+1) 1.. W (M+1) d with mutually independent d-dimensional Brownian motions W (M+1) and W (M+1). Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
40 Pricing and Calibration Specification Analysis Specification Analysis: Continuous Part Let γ i (t) = c i g(t i t)e i, e i R d, where c i > 0 are loading factors g i ( ) is a scalar volatility function e i are unit vectors coming from the decomposition of the correlation matrix ζ ζ ij = e i e j, 1 i, j M, be a deterministic volatility structure of the input Libor market model calibrated to ATM caps and ATM swaptions. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
41 Pricing and Calibration Specification Analysis Specification Analysis: Continuous Part Define a new time independent volatility structure via β i β j = min(i,j) 1 min(i, j) k=1 1 T k Tk 0 γ i (t)γ j (t) dt. Remark The covariance of the process ξ i (t) := t 0 Γ i (t)dw (M+1) satisfies cov(ξ i (t), ξ j (t)) t 0 γ i (t)γ j (t)ds and is approximately the same as in the input LMM. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
42 Pricing and Calibration Specification Analysis Specification Analysis: Continuous Part Define a new time independent volatility structure via β i β j = min(i,j) 1 min(i, j) k=1 1 T k Tk 0 γ i (t)γ j (t) dt. Remark The covariance of the process ξ i (t) := t 0 Γ i (t)dw (M+1) satisfies cov(ξ i (t), ξ j (t)) t 0 γ i (t)γ j (t)ds and is approximately the same as in the input LMM. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
43 Pricing and Calibration Specification Analysis Specification Analysis: Continuous Part Two possible specifications for the volatility process v Stochastic Volatility Heston Model dv k = κ k (1 v k )dt +σ k ϱ k vk dw (M+1) k +σ k (1 ϱ 2 k ) v k dv (M+1) k, Stochastic Volatility BN Model dv k = κ k v k dt + σ k ϱ k dw (M+1) k + σ k (1 ϱ 2 (M+1) k )dv k. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
44 Pricing and Calibration Specification Analysis Specification Analysis: Continuous Part Two possible specifications for the volatility process v Stochastic Volatility Heston Model dv k = κ k (1 v k )dt +σ k ϱ k vk dw (M+1) k +σ k (1 ϱ 2 k ) v k dv (M+1) k, Stochastic Volatility BN Model dv k = κ k v k dt + σ k ϱ k dw (M+1) k + σ k (1 ϱ 2 (M+1) k )dv k. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
45 Pricing and Calibration Specification Analysis Specification Analysis: Continuous Part Two possible specifications for the volatility process v Stochastic Volatility Heston Model dv k = κ k (1 v k )dt +σ k ϱ k vk dw (M+1) k +σ k (1 ϱ 2 k ) v k dv (M+1) k, Stochastic Volatility BN Model dv k = κ k v k dt + σ k ϱ k dw (M+1) k + σ k (1 ϱ 2 (M+1) k )dv k. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
46 Pricing and Calibration Specification Analysis Specification Analysis It holds Φ C P M+1 (z; T ) = Φ C D,P M+1 (z; T ) Φ C SV,P M+1 (z; T ), where Φ C D,P M+1 (z; T ) = exp ( 12 ( ) ) T θ2m (T ) z 2 + iz, θm 2 (T ) = γ M 2 dt 0 and Φ C SV,P M+1 (z; T ) = exp (A M (z; T ) + B M (z; T )) Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
47 Pricing and Calibration Specification Analysis Specification Analysis In particular and A M (z; T ) = κ { [ M 1 gm e d ]} MT σm 2 (a M + d M )T 2 ln 1 g M B M (z; T ) = (a M + d M )(1 e d MT ) σ 2 M (1 g Me d MT ) a M = κ M iϱ M ω M z d M = am 2 + ω2 M (z2 + iz) g M = a M + d M a M d M, ω M = r SV β MM σ M Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
48 Pricing and Calibration Specification Analysis Specification Analysis As can be easily seen ) A M (z; T ) lim = α M ω M (iϱ M + 1 ϱ 2 z z M T and with B M (z; T ) 1 ϱ 2 M lim = + iϱ M z z α M := κ M σ 2 M σ M Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
49 Pricing and Calibration Specification Analysis Specification Analysis Let us take φ i (u, t) = u β i, then the characteristic function of the jump part is given by ( ) Φ J P M+1 (z; T ) = exp λt (e izv 1)µ M (v) dv, R where µ M is the density of u β M (t). Observation Due to the Riemann-Lebesgue theorem Φ J P M+1 (z; T ) exp( λt ), z. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
50 Pricing and Calibration Specification Analysis Specification Analysis Let us take φ i (u, t) = u β i, then the characteristic function of the jump part is given by ( ) Φ J P M+1 (z; T ) = exp λt (e izv 1)µ M (v) dv, R where µ M is the density of u β M (t). Observation Due to the Riemann-Lebesgue theorem Φ J P M+1 (z; T ) exp( λt ), z. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
51 Pricing and Calibration Specification Analysis Specification Analysis: Asymptotic Properties Computing sequentially L 2 := lim z log(φ PM+1 (z; T ))/z 2, we get L 1 := lim z [ log(φpm+1 (z; T ))/z (z + i)l 2 ], [ ] L 0 := lim log(φ PM+1 (z; T )) (z 2 + iz)l 2 zl 1, z L 0 = λ, L 2 = 1 2 θ2 M (T ) and Re L 1 = 1 ϱ 2 M σ M α M ω M 1 ϱ 2 M T, Im L 1 = ϱ M σ M α M ω M ϱ M T Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
52 Pricing and Calibration Specification Analysis Specification Analysis: Asymptotic Properties Computing sequentially L 2 := lim z log(φ PM+1 (z; T ))/z 2, we get L 1 := lim z [ log(φpm+1 (z; T ))/z (z + i)l 2 ], [ ] L 0 := lim log(φ PM+1 (z; T )) (z 2 + iz)l 2 zl 1, z L 0 = λ, L 2 = 1 2 θ2 M (T ) and Re L 1 = 1 ϱ 2 M σ M α M ω M 1 ϱ 2 M T, Im L 1 = ϱ M σ M α M ω M ϱ M T Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
53 Pricing and Calibration Specification Analysis Specification Analysis: Asymptotic Properties Computing sequentially L 2 := lim z log(φ PM+1 (z; T ))/z 2, we get L 1 := lim z [ log(φpm+1 (z; T ))/z (z + i)l 2 ], [ ] L 0 := lim log(φ PM+1 (z; T )) (z 2 + iz)l 2 zl 1, z L 0 = λ, L 2 = 1 2 θ2 M (T ) and Re L 1 = 1 ϱ 2 M σ M α M ω M 1 ϱ 2 M T, Im L 1 = ϱ M σ M α M ω M ϱ M T Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
54 Pricing and Calibration Calibration Procedure Parameters Estimation: Linearization Observation From the knowledge of L 1 (T ) for two different T one can reconstruct all parameters of the SV process Theorem Ψ PM+1 (z; T ) := log(φ PM+1 (z; T )) = L 2 (z 2 + iz) + L 1 z + L 0 + R 0 + R 1 (z), where R 0 = R 0 (α M, κ M, ϱ M, ω M ) is a constant not depending on λ and R 1 (z) 0, z. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
55 Pricing and Calibration Calibration Procedure Parameters Estimation: Linearization Observation From the knowledge of L 1 (T ) for two different T one can reconstruct all parameters of the SV process Theorem Ψ PM+1 (z; T ) := log(φ PM+1 (z; T )) = L 2 (z 2 + iz) + L 1 z + L 0 + R 0 + R 1 (z), where R 0 = R 0 (α M, κ M, ϱ M, ω M ) is a constant not depending on λ and R 1 (z) 0, z. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
56 Pricing and Calibration Calibration Procedure Parameters Estimation: Projection Estimators We find estimates for L 2, L 1 and L 0 in the form of weighted averages L 2,U := L 1,U := L 0,U := Re( Ψ PM+1 (u))w2 U (u) du, Im( Ψ PM+1 (u))w1 U (u) du i L 2,U, Re( Ψ PM+1 (u))w0 U (u) du R 0 with ( ) Ψ PM+1 (u) := ln 1 u(u + i)f{õm}(u + i). Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
57 Pricing and Calibration Parameters Estimation Calibration Procedure The weights are given by w U 2 = U 3 w 2 (u/u), w U 1 = U 2 w 1 (u/u), w U 0 = U 1 w 0 (u/u), where 1 1 w 2 (u)du = 0, 1 1 uw 2 (u)du = 0, 1 1 u 2 w 2 (u)du = 1, w 0 (u)du = 1, w 1 (u) du = 0, uw 0 (u)du = 0, uw 1 (u) du = 1, 1 1 u 2 w 0 (u)du = 0. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
58 Pricing and Calibration Calibration Procedure Parameters Estimation: Jump distribution Define F{ µ M }(z) = Ψ PM+1 (z; T ) L 2 (z 2 + iz) L 1 z L 0 R 0 or equivalently µ M := F 1 [( ΨPM+1 ( ; T ) L 2 ( 2 + i ) L 1 L 0 R ) ] 0 1 [ U,U] ( ) Remark Due to lack of data and numerical errors µ M may not be a density and needs to be corrected. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
59 Pricing and Calibration Calibration Procedure Parameters Estimation: Jump distribution Define F{ µ M }(z) = Ψ PM+1 (z; T ) L 2 (z 2 + iz) L 1 z L 0 R 0 or equivalently µ M := F 1 [( ΨPM+1 ( ; T ) L 2 ( 2 + i ) L 1 L 0 R ) ] 0 1 [ U,U] ( ) Remark Due to lack of data and numerical errors µ M may not be a density and needs to be corrected. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
60 Pricing and Calibration Calibration Procedure Parameters Estimation: Further optimization Upon finding ( L0,U, L 1,U, L ) 2,U T := we may consider T as a final set of parameters or consider nonlinear least-squares ( σ M, ϱ M, κ M, λ ) J (T ) = N w i CM T (K i) C M (K i ) 2 i=1 and minimize J (T ) over the parametric set S R 4 taking as initial value T. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
61 Pricing and Calibration Calibration Procedure Approximative dynamics of L i under P i+1 It holds approximately dl i Γ i dw (i+1) + e u β j (µ ν (i+1) )(dt, du), L i E where dw (i+1) is a standard Brownian motion under P i+1 and ν (i+1) (dt, du) = ν (M+1) (dt, du) M j=i+1 ( 1 + δ jl j (0)e u β j 1 + δ j L j (0) ). Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
62 Pricing and Calibration Calibration Procedure Approximative dynamics of v k under P i+1 By freezing the Libors at their initial values we obtain an approximative v k dynamics ( ) ( dv k κ (i+1) k θ (i+1) ) (i+1) k v k dt + σ k vk ϱ k d W k + 1 ϱ 2 (i+1) kdw k with reversion speed parameter κ (i+1) k := and mean reversion level κ k r SV σ k ϱ k M j=i+1 θ (i+1) k := κ k κ (i+1) k δ j L j (0) 1 + δ j L j (0) β jk,. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
63 Pricing and Calibration Calibration Procedure Pricing Caplets under P M+1 The price of j-th caplet at time zero can be alternatively written as [ ] Bj+1 (T j ) C j (K ) = δ j B M+1 (0)E PM+1 B M+1 (T j ) (L j(t j ) K ) + Note that B j+1 (T j ) B M+1 (T j ) = M k=j+1 (1 + δ k L k (T j )) = M k=j+1 (1 + δ k )E ξ exp M k=j+1 ξ k ln(l k (T j )), where {ξ k } M k=j+1 are independent random variables and each ξ k takes two values 0 and 1 with probabilities 1/(1 + δ k ) and δ k /(1 + δ k ). Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
64 Pricing and Calibration Calibration Procedure Pricing Caplets under P M+1 Thus, F{O j }(z) = 1 E ξφ M+1 (z i, ξ j+1,..., ξ M ), z(z i) where Φ M+1 (z j, z j+1,..., z M ) is the joint characteristic function of (ln(l j (T j )),..., ln(l M (T j ))) under P M+1. Remark Instead of terminal measure P M+1 we could consider P l+1 with 1 < l < M + 1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
65 Pricing and Calibration Calibration Procedure Pricing Caplets under P M+1 Thus, F{O j }(z) = 1 E ξφ M+1 (z i, ξ j+1,..., ξ M ), z(z i) where Φ M+1 (z j, z j+1,..., z M ) is the joint characteristic function of (ln(l j (T j )),..., ln(l M (T j ))) under P M+1. Remark Instead of terminal measure P M+1 we could consider P l+1 with 1 < l < M + 1. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
66 Calibration in work Calibration results for Caplet volas for different caplet periods Caplet Volatilities [ 17.5, 18 ] [ 15, 15.5 ] Caplet Volatilities [ 12.5, 13 ] [ 10, 10.5 ] Strikes Strikes Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
67 Bibliography Belomestny, D. and Spokoiny, V. Spatial aggregation of local likelihood estimates with applications to classification, Annals of Statistics, 2007, 35(5), Belomestny, D. and Reiß, M. Spectral calibration of exponential Lévy models, Finance and Stochastics, 2006, 10(4), Belomestny, D. and Schoenmakers, J. A jump-diffusion Libor model and its robust calibration, SFB649 Discussion Paper, 2006, 037. Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar / 43
Extended 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 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 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 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 informationQuadratic hedging in affine stochastic volatility models
Quadratic hedging in affine stochastic volatility models Jan Kallsen TU München Pittsburgh, February 20, 2006 (based on joint work with F. Hubalek, L. Krawczyk, A. Pauwels) 1 Hedging problem S t = S 0
More informationINTEREST RATES AND FX MODELS
INTEREST RATES AND FX MODELS 6. LIBOR Market Model Andrew Lesniewski Courant Institute of Mathematical Sciences New York University New York March 6, 2013 2 Interest Rates & FX Models Contents 1 Introduction
More informationParametric Inference and Dynamic State Recovery from Option Panels. Torben G. Andersen
Parametric Inference and Dynamic State Recovery from Option Panels Torben G. Andersen Joint work with Nicola Fusari and Viktor Todorov The Third International Conference High-Frequency Data Analysis in
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 informationA Hybrid Commodity and Interest Rate Market Model
A Hybrid Commodity and Interest Rate Market Model University of Technology, Sydney June 1 Literature A Hybrid Market Model Recall: The basic LIBOR Market Model The cross currency LIBOR Market Model LIBOR
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 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 informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More informationOptimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models
Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University University of Missouri-Kansas City Department of Mathematics
More informationIMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS
IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS DENIS BELOMESTNY AND MARKUS REISS 1. Introduction The aim of this report is to describe more precisely how the spectral calibration method
More informationParameters Estimation in Stochastic Process Model
Parameters Estimation in Stochastic Process Model A Quasi-Likelihood Approach Ziliang Li University of Maryland, College Park GEE RIT, Spring 28 Outline 1 Model Review The Big Model in Mind: Signal + Noise
More informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
More informationCredit Risk in Lévy Libor Modeling: Rating Based Approach
Credit Risk in Lévy Libor Modeling: Rating Based Approach Zorana Grbac Department of Math. Stochastics, University of Freiburg Joint work with Ernst Eberlein Croatian Quants Day University of Zagreb, 9th
More informationHedging of swaptions in a Lévy driven Heath-Jarrow-Morton framework
Hedging of swaptions in a Lévy driven Heath-Jarrow-Morton framework Kathrin Glau, Nele Vandaele, Michèle Vanmaele Bachelier Finance Society World Congress 2010 June 22-26, 2010 Nele Vandaele Hedging of
More informationHedging under Arbitrage
Hedging under Arbitrage Johannes Ruf Columbia University, Department of Statistics Modeling and Managing Financial Risks January 12, 2011 Motivation Given: a frictionless market of stocks with continuous
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 informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
More informationThe Valuation of Bermudan Guaranteed Return Contracts
The Valuation of Bermudan Guaranteed Return Contracts Steven Simon 1 November 2003 1 K.U.Leuven and Ente Luigi Einaudi Abstract A guaranteed or minimum return can be found in different financial products,
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 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 informationStochastic modelling of electricity markets Pricing Forwards and Swaps
Stochastic modelling of electricity markets Pricing Forwards and Swaps Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Clip for this slide Pricing
More informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More informationSupplementary Appendix to The Risk Premia Embedded in Index Options
Supplementary Appendix to The Risk Premia Embedded in Index Options Torben G. Andersen Nicola Fusari Viktor Todorov December 214 Contents A The Non-Linear Factor Structure of Option Surfaces 2 B Additional
More informationModern Methods of Option Pricing
Modern Methods of Option Pricing Denis Belomestny Weierstraß Institute Berlin Motzen, 14 June 2007 Denis Belomestny (WIAS) Modern Methods of Option Pricing Motzen, 14 June 2007 1 / 30 Overview 1 Introduction
More informationParametric Inference and Dynamic State Recovery from Option Panels. Nicola Fusari
Parametric Inference and Dynamic State Recovery from Option Panels Nicola Fusari Joint work with Torben G. Andersen and Viktor Todorov July 2012 Motivation Under realistic assumptions derivatives are nonredundant
More informationPolicy iterated lower bounds and linear MC upper bounds for Bermudan style derivatives
Finance Winterschool 2007, Lunteren NL Policy iterated lower bounds and linear MC upper bounds for Bermudan style derivatives Pricing complex structured products Mohrenstr 39 10117 Berlin schoenma@wias-berlin.de
More informationInterest Rate Curves Calibration with Monte-Carlo Simulatio
Interest Rate Curves Calibration with Monte-Carlo Simulation 24 june 2008 Participants A. Baena (UCM) Y. Borhani (Univ. of Oxford) E. Leoncini (Univ. of Florence) R. Minguez (UCM) J.M. Nkhaso (UCM) A.
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 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 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 informationA Closed-form Solution for Outperfomance Options with Stochastic Correlation and Stochastic Volatility
A Closed-form Solution for Outperfomance Options with Stochastic Correlation and Stochastic Volatility Jacinto Marabel Romo Email: jacinto.marabel@grupobbva.com November 2011 Abstract This article introduces
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 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 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 informationDiscrete time interest rate models
slides for the course Interest rate theory, University of Ljubljana, 2012-13/I, part II József Gáll University of Debrecen, Faculty of Economics Nov. 2012 Jan. 2013, Ljubljana Introduction to discrete
More informationPolynomial processes in stochastic portofolio theory
Polynomial processes in stochastic portofolio theory Christa Cuchiero University of Vienna 9 th Bachelier World Congress July 15, 2016 Christa Cuchiero (University of Vienna) Polynomial processes in SPT
More informationStochastic volatility modeling in energy markets
Stochastic volatility modeling in energy markets Fred Espen Benth Centre of Mathematics for Applications (CMA) University of Oslo, Norway Joint work with Linda Vos, CMA Energy Finance Seminar, Essen 18
More informationGRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS
GRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS Patrick GAGLIARDINI and Christian GOURIÉROUX INTRODUCTION Risk measures such as Value-at-Risk (VaR) Expected
More informationInsiders Hedging in a Stochastic Volatility Model with Informed Traders of Multiple Levels
Insiders Hedging in a Stochastic Volatility Model with Informed Traders of Multiple Levels Kiseop Lee Department of Statistics, Purdue University Mathematical Finance Seminar University of Southern California
More informationSLE and CFT. Mitsuhiro QFT2005
SLE and CFT Mitsuhiro Kato @ QFT2005 1. Introduction Critical phenomena Conformal Field Theory (CFT) Algebraic approach, field theory, BPZ(1984) Stochastic Loewner Evolution (SLE) Geometrical approach,
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 informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationArbitrage of the first kind and filtration enlargements in semimartingale financial models. Beatrice Acciaio
Arbitrage of the first kind and filtration enlargements in semimartingale financial models Beatrice Acciaio the London School of Economics and Political Science (based on a joint work with C. Fontana and
More informationStochastic Volatility and Jump Modeling in Finance
Stochastic Volatility and Jump Modeling in Finance HPCFinance 1st kick-off meeting Elisa Nicolato Aarhus University Department of Economics and Business January 21, 2013 Elisa Nicolato (Aarhus University
More informationSensitivity Analysis on Long-term Cash flows
Sensitivity Analysis on Long-term Cash flows Hyungbin Park Worcester Polytechnic Institute 19 March 2016 Eastern Conference on Mathematical Finance Worcester Polytechnic Institute, Worceseter, MA 1 / 49
More informationPolynomial Models in Finance
Polynomial Models in Finance Martin Larsson Department of Mathematics, ETH Zürich based on joint work with Damir Filipović, Anders Trolle, Tony Ware Risk Day Zurich, 11 September 2015 Flexibility Tractability
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 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 informationThe Term Structure of Interest Rates under Regime Shifts and Jumps
The Term Structure of Interest Rates under Regime Shifts and Jumps Shu Wu and Yong Zeng September 2005 Abstract This paper develops a tractable dynamic term structure models under jump-diffusion and regime
More informationHedging under arbitrage
Hedging under arbitrage Johannes Ruf Columbia University, Department of Statistics AnStAp10 August 12, 2010 Motivation Usually, there are several trading strategies at one s disposal to obtain a given
More informationLongevity risk: past, present and future
Longevity risk: past, present and future Xiaoming Liu Department of Statistical & Actuarial Sciences Western University Longevity risk: past, present and future Xiaoming Liu Department of Statistical &
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 informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
More informationA Correlation-sensitive Calibration of a Stochastic Volatility LIBOR Market Model
A Correlation-sensitive Calibration of a Stochastic Volatility LIBOR Market Model Man Kuan Wong Lady Margaret Hall University of Oxford A thesis submitted in partial fulfillment of the MSc in Mathematical
More informationGirsanov s Theorem. Bernardo D Auria web: July 5, 2017 ICMAT / UC3M
Girsanov s Theorem Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 5, 2017 ICMAT / UC3M Girsanov s Theorem Decomposition of P-Martingales as Q-semi-martingales Theorem
More informationGenetics and/of basket options
Genetics and/of basket options Wolfgang Karl Härdle Elena Silyakova Ladislaus von Bortkiewicz Chair of Statistics Humboldt-Universität zu Berlin http://lvb.wiwi.hu-berlin.de Motivation 1-1 Basket derivatives
More informationMultiple stochastic volatility extension of the Libor market model and its implementation
Multiple stochastic volatility extension of the Libor maret model and its implementation Denis Belomestny 1, Stanley Mathew 2, and John Schoenmaers 1 April 9, 29 Keywords: Libor modeling, stochastic volatility,
More informationA More Detailed and Complete Appendix for Macroeconomic Volatilities and Long-run Risks of Asset Prices
A More Detailed and Complete Appendix for Macroeconomic Volatilities and Long-run Risks of Asset Prices This is an on-line appendix with more details and analysis for the readers of the paper. B.1 Derivation
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 informationExponential utility maximization under partial information
Exponential utility maximization under partial information Marina Santacroce Politecnico di Torino Joint work with M. Mania AMaMeF 5-1 May, 28 Pitesti, May 1th, 28 Outline Expected utility maximization
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 informationNormal Inverse Gaussian (NIG) Process
With Applications in Mathematical Finance The Mathematical and Computational Finance Laboratory - Lunch at the Lab March 26, 2009 1 Limitations of Gaussian Driven Processes Background and Definition IG
More informationTesting for non-correlation between price and volatility jumps and ramifications
Testing for non-correlation between price and volatility jumps and ramifications Claudia Klüppelberg Technische Universität München cklu@ma.tum.de www-m4.ma.tum.de Joint work with Jean Jacod, Gernot Müller,
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 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 informationWeierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions
Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions Hilmar Mai Mohrenstrasse 39 1117 Berlin Germany Tel. +49 3 2372 www.wias-berlin.de Haindorf
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 informationEnlargement of filtration
Enlargement of filtration Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 6, 2017 ICMAT / UC3M Enlargement of Filtration Enlargement of Filtration ([1] 5.9) If G is a
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 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 informationBIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS
BIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS PRICING EMMS014S7 Tuesday, May 31 2011, 10:00am-13.15pm
More informationCalibration Risk for Exotic Options
SFB 649 Discussion Paper 2006-001 Calibration Risk for Exotic Options Kai Detlefsen* Wolfgang K. Härdle** * CASE - Center for Applied Statistics and Economics, Humboldt-Universität zu Berlin, Germany SFB
More informationOn the Ross recovery under the single-factor spot rate model
.... On the Ross recovery under the single-factor spot rate model M. Kijima Tokyo Metropolitan University 11/08/2016 Kijima (TMU) Ross Recovery SMU @ August 11, 2016 1 / 35 Plan of My Talk..1 Introduction:
More informationL 2 -theoretical study of the relation between the LIBOR market model and the HJM model Takashi Yasuoka
Journal of Math-for-Industry, Vol. 5 (213A-2), pp. 11 16 L 2 -theoretical study of the relation between the LIBOR market model and the HJM model Takashi Yasuoka Received on November 2, 212 / Revised on
More informationUniversität Regensburg Mathematik
Universität Regensburg Mathematik Modeling financial markets with extreme risk Tobias Kusche Preprint Nr. 04/2008 Modeling financial markets with extreme risk Dr. Tobias Kusche 11. January 2008 1 Introduction
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 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 information16. Inflation-Indexed Swaps
6. Inflation-Indexed Swaps Given a set of dates T,...,T M, an Inflation-Indexed Swap (IIS) is a swap where, on each payment date, Party A pays Party B the inflation rate over a predefined period, while
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 informationOptimal trading strategies under arbitrage
Optimal trading strategies under arbitrage Johannes Ruf Columbia University, Department of Statistics The Third Western Conference in Mathematical Finance November 14, 2009 How should an investor trade
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 informationRMSC 4005 Stochastic Calculus for Finance and Risk. 1 Exercises. (c) Let X = {X n } n=0 be a {F n }-supermartingale. Show that.
1. EXERCISES RMSC 45 Stochastic Calculus for Finance and Risk Exercises 1 Exercises 1. (a) Let X = {X n } n= be a {F n }-martingale. Show that E(X n ) = E(X ) n N (b) Let X = {X n } n= be a {F n }-submartingale.
More informationIlliquidity, Credit risk and Merton s model
Illiquidity, Credit risk and Merton s model (joint work with J. Dong and L. Korobenko) A. Deniz Sezer University of Calgary April 28, 2016 Merton s model of corporate debt A corporate bond is a contingent
More informationLecture 1: Lévy processes
Lecture 1: Lévy processes A. E. Kyprianou Department of Mathematical Sciences, University of Bath 1/ 22 Lévy processes 2/ 22 Lévy processes A process X = {X t : t 0} defined on a probability space (Ω,
More informationModeling the dependence between a Poisson process and a continuous semimartingale
1 / 28 Modeling the dependence between a Poisson process and a continuous semimartingale Application to electricity spot prices and wind production modeling Thomas Deschatre 1,2 1 CEREMADE, University
More informationHedging of Contingent Claims under Incomplete Information
Projektbereich B Discussion Paper No. B 166 Hedging of Contingent Claims under Incomplete Information by Hans Föllmer ) Martin Schweizer ) October 199 ) Financial support by Deutsche Forschungsgemeinschaft,
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 informationPDE Approach to Credit Derivatives
PDE Approach to Credit Derivatives Marek Rutkowski School of Mathematics and Statistics University of New South Wales Joint work with T. Bielecki, M. Jeanblanc and K. Yousiph Seminar 26 September, 2007
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 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 informationChapter 14. The Multi-Underlying Black-Scholes Model and Correlation
Chapter 4 The Multi-Underlying Black-Scholes Model and Correlation So far we have discussed single asset options, the payoff function depended only on one underlying. Now we want to allow multiple underlyings.
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 informationNonlinear Filtering in Models for Interest-Rate and Credit Risk
Nonlinear Filtering in Models for Interest-Rate and Credit Risk Rüdiger Frey 1 and Wolfgang Runggaldier 2 June 23, 29 3 Abstract We consider filtering problems that arise in Markovian factor models for
More informationA Multifrequency Theory of the Interest Rate Term Structure
A Multifrequency Theory of the Interest Rate Term Structure Laurent Calvet, Adlai Fisher, and Liuren Wu HEC, UBC, & Baruch College Chicago University February 26, 2010 Liuren Wu (Baruch) Cascade Dynamics
More information