A Benchmark Approach To Quantitative Finance - Lecture Notes
A Benchmark Approach To Quantitative Finance - Lecture Notes
Eckhard Platen
School of Finance and Economics and Department of Mathematical Sciences
University of Technology, Sydney
Platen, E.: A benchmark approach to finance. Mathematical Finance 16(1), 131-151 (2006).
Bühlmann, H. & Platen, E.: A discrete time benchmark approach for insurance and finance
ASTIN Bulletin 33(2), 153-172 (2003).
Contents
1 Asset Pricing 1
3 Portfolio Optimization 91
References 337
1 Asset Pricing
Cochrane (2001), . . .
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.1: Logarithms of savings bond, fair zero coupon bond and savings
account.
c Copyright Eckhard Platen 10 BA - Toronto C1 2
Financial Market
Stj
j ∈ {0, 1, . . . , d}
• savings account
St0
t≥0
predictable
• portfolio
d
X
Stδ = δtj Stj
j=0
• self-financing
d
X
dStδ = δtj dStj
j=0
S0δ = x > 0
=⇒ S δ ∈ Vx+
Bühlmann & Pl. (2003), Goll & Kallsen (2003), Karatzas & Kardaras (2007)
7,000
6,000
5,000
DWI
4,000
MCI
3,000
2,000
1,000
0
1/1/73 28/9/76 25/6/80 22/3/84 18/12/87 15/9/91 12/6/95 9/3/99 4/12/02 31/8/06
Figure 1.2: Constructed MCI, DWI, EWI and WSI. In the long term the WSI
and EWI outperform all other indices, Le & Platen (2006).
Supermartingale Property
0≤t≤s<∞
Ŝ δ supermartingale
g δ ≤ g δ∗ .
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.3: Logarithms of fair zero coupon bond and savings account.
δ
St+h
Aδt,h =
Stδ
• expected growth
δ
gt,h = Et ln Aδt,h
t, h ≥ 0
• expected log-utility
Kelly (1956)
Hakansson (1971a)
Merton (1973a)
Roll (1973)
Markowitz (1976)
Stδ
P > 0 > 0.
0.8
0.6
0.4
0.2
-0.2
1930 1940 1950 1960 1970 1980 1990 2000
time
0.0016
0.0014
0.0012
0.001
0.0008
0.0006
0.0004
0.0002
-0.0002
-0.0004
1920 1922 1924 1926 1928
time
(no trend)
Ŝtδ = Et Ŝsδ
0 ≤ t ≤ s < ∞.
0.1
0.08
0.06
0.04
0.02
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.6: Benchmarked savings bond and benchmarked fair zero coupon
bond.
c Copyright Eckhard Platen 10 BA - Toronto C1 27
Law of the Minimal Price
Pl. (2008)
• least expensive
• minimal hedge
HT !
HT
E0 <∞
STδ∗
• StδH fair if
!
StδH HT
ŜtδH = = Et
Stδ∗ STδ∗
Corollary 1.14
=⇒
actuarial pricing
risk neutral pricing
pricing with stochastic discount factor
pricing with numeraire change
pricing with deflator
pricing with state pricing density
pricing with pricing kernel
pricing with numeraire portfolio
=⇒ rigorous derivation of
S00
S0δH = E0 ΛT HT
ST0
Ŝt0
Λt = supermartingale
Ŝ00
1 = Λ0 ≥ E0 (ΛT )
=⇒
S00
E0 ΛT S 0 HT
S0δH ≤ T
E0 (ΛT )
90
80
70
60
50
40
30
20
10
0
1930 1940 1950 1960 1970 1980 1990 2000
time
1.4
1.2
0.8
0.6
0.4
0.2
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.8: Radon-Nikodym derivative and total mass of putative risk neutral
measure.
c Copyright Eckhard Platen 10 BA - Toronto C1 35
• special case when savings account is fair:
dQ
=⇒ ΛT = dP
forms martingale; E0 (ΛT ) = 1;
Bayes’ formula =⇒
risk neutral pricing formula
S00
S0δH = E0Q HT
ST0
-1
1930 1940 1950 1960 1970 1980 1990 2000
time
• savings bond
St0
P ∗ (t, T ) =
ST0
P (t, T )
P̂ (t, T ) = martingale
Stδ∗
Stδ∗
S̄tδ∗ =
St0
αt = α exp{η t}
=⇒
-1
-2
1930 1940 1950 1960 1970 1980 1990 2000
time
√
• quadratic variation of Yt
p 1
d Yt = . . . + dWt
2
it h√ i
X p q 2 t
Ytℓ − Ytℓ−1 ≈ Y =
ℓ=1
t 4
25
20
15
10
0
1930 1940 1950 1960 1970 1980 1990 2000
time
√
Figure 1.11: Quadratic Variation of Yt .
( )!
∗ 2η S̄tδ∗
P (t, T ) = P (t, T ) 1 − exp −
α (exp{η T } − exp{η t})
• initial prices
P ∗ (0, T ) = 0.0335
P (0, T ) = 0.00077
P (0,T )
P ∗ (0,T )
< 0.023 =⇒ 2.3%
no strong arbitrage
0.12
0.1
0.08
0.06
0.04
0.02
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.12: Benchmarked savings bond and benchmarked fair zero coupon
bond.
c Copyright Eckhard Platen 10 BA - Toronto C1 46
1
savings bond
fair zero coupon bond
0.9 savings account
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.13: Savings bond, fair zero coupon bond and savings account.
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.14: Logarithms of fair zero coupon bond and savings account.
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.15: Logarithms of fair zero coupon bond and savings account.
∂ P̄ (t, T )
δt∗ =
∂ S̄tδ∗
( )
∗ −2 η S̄tδ∗
= P (0, T ) exp
α (exp{η T } − exp{η t})
2η
×
α (exp{η T } − exp{η t})
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 1.16: Logarithm of zero coupon bond and self-financing hedge port-
folio.
c Copyright Eckhard Platen 10 BA - Toronto C1 52
7e-005
6e-005
5e-005
4e-005
3e-005
2e-005
1e-005
-1e-005
-2e-005
1930 1940 1950 1960 1970 1980 1990 2000
time
0.8
0.6
0.4
0.2
W 1, W 2, . . . , W d, d ∈ {1, 2, . . .}
j ∈ {0, 1, . . . , d}
k=1
Z t
St0 = exp rs ds
0
solution of relation
bt θ t = at − rt 1
where
at = (a1t , a2t , . . . , adt )⊤
1 = (1, 1, . . . , 1)⊤
=⇒
θ t = b−1
t (at − rt 1)
( d
)
X
dStj = Stj rt dt + bj,k k k
t (θt dt + dWt )
k=1
• portfolio
d
X
Stδ = δtj Stj
j=0
d
X
dStδ = δtj dStj
j=0
j Stj
πδ,t = δtj
Stδ
j ∈ {0, 1, . . . , d}
d
X j
πδ,t =1
j=0
d
!
X
dStδ Stδ bkδ,t θtk dWtk
= rt dt + dt +
k=1
d
X j
bkδ,t = πδ,t bj,k
t
j=1
d
X
d ln(Stδ ) = gtδ dt + bkδ,t dWtk
k=1
d
1
X
gtδ = rt + bkδ,t θtk − bkδ,t
k=1
2
d X
d d
!
X j 1X
gtδ = rt + πδ,t bj,k
t θtk − ℓ
πδ,t bℓ,k
t
k=1 j=1
2 ℓ=1
d d
!
X X
0= bj,k
t θtk − πδℓ∗ ,t bℓ,k
t
k=1 ℓ=1
=⇒
θ⊤ ⊤
t = π δ∗ ,t bt
π δ∗ ,t = (πδ1∗ ,t , . . . , πδd∗ ,t )⊤
= (b−1 ⊤
t ) θt
=⇒ NP
" d
# d
!
X X
dStδ∗ = Stδ∗ rt + (θtk )2 dt + θtk dWtk
k=1 k=1
finite NP exists
Stδ
Ŝtδ =
Stδ∗
d X
X d
j
dStδ Stδ bj,k θtk dWtk
= rt dt + πδ,t t dt +
k=1 j=0
d
!
X
dStδ∗ Stδ∗ θtk θtk dWtk
= rt dt + dt +
k=1
d X
X d
j
dŜtδ = −Ŝtδ πδ,t sj,k k
t dWt
k=1 j=0
" d
#
1 σ
X
dStj = Stj r+σ 2
1+ √ dt + √ dWtk + σ dWtj
d d k=1
j ∈ {1, 2, . . . , d}
√ −1
πδ0∗ ,t = 1+ d
• NP SDE
d
!
2 σ X
dStδ∗ Stδ∗ dWtk
= r+σ dt + √
d k=1
σ = 0.15, r = 0.05
20
15
10
0
0 5 10 15 20 25 30 35
time
0
0 5 10 15 20 25 30 35
time
j K2
πδ,t ≤ 1
d 2 +K1
a.s. for all j ∈ {0, 1, . . . , d} and t ∈ [0, T ] for K1 , K2 ∈ (0, ∞) and
K3 ∈ {1, 2, . . .}, independent of d ∈ {K3 , K3 + 1, . . .}.
Stj
Ŝtj =
Stδ∗
d
X
dŜtj = −Ŝtj sj,k k
t dWt
k=1
sj,k
t = θt
k
− bj,k
t
d
X
ŝkt = |sj,k
t |
j=0
2
k
E ŝt ≤ K5 .
2
d
X d
X j
Rδd (t) = πδ,t sj,k
t
k=1 j=0
benchmarked NP constant =⇒
Rδd∗ (t) = 0
Rδd (t)
lim P > ε = 0.
d→∞
• model independent
• similarity to CLT
j ∈ {0, 1, . . . , d}
0
0 5 10 15 20 25 30 35
time
0
0 5 10 15 20 25 30 35
time
100,000
80,000
60,000
40,000
20,000
0
1/1/73 28/9/76 25/6/80 22/3/84 18/12/87 15/9/91 12/6/95 9/3/99 4/12/02 31/8/06
Figure 2.5: World industry sector stock market indices as primary security
accounts.
7,000
6,000
5,000
DWI
4,000
MCI
3,000
2,000
1,000
0
1/1/73 28/9/76 25/6/80 22/3/84 18/12/87 15/9/91 12/6/95 9/3/99 4/12/02 31/8/06
Figure 2.6: Constructed MCI, DWI, EWI and WSI. In the long term the WSI
and EWI outperform all other indices, Le & Platen (2006).
400
350
300
250
200
150
100
50
0
0 5 9 14 18 23 27 32
50
40
30
20
10
0
0 5 9 14 18 23 27 32
EWI
90 GOP
80
70
60
50
40
30
20
10
0
0 5 9 14 18 23 27 32
Market index
90
GOP
80
70
60
50
40
30
20
10
0
0 5 9 14 18 23 27 32
Markowitz (1959)
Merton (1973a)
• examples
power utility
1
U (x) = xγ
γ
for γ 6= 0 and γ < 1
log-utility
U (x) = ln(x)
x
5 10 15 20
-5
-10
-15
Figure 3.1: Examples for power utility (upper graph) and log-utility (lower
graph).
c Copyright Eckhard Platen 10 BA - Toronto C3 94
• assume Ŝ 0 is scalar diffusion
S̄0δ = S0 > 0
Ŝtδ̃ = û(t, Ŝt0 ) = Et U ′−1 λ ŜT0 ŜT0 ,
1
=⇒ invest of wealth in NP and remainder in savings account
Jtδ̃
′ 1
U (x) =
x
′−1 1
U (y) =
y
′′ 1
U (x) = −
x2
utility concave
!
′−1
1 1
û(t, Ŝt0 ) = Et U λ ŜT0 ŜT0 = Et ŜT
0
=
λ ŜT0 λ
Lagrange multiplier
• expected log-utility
δ̃ δ∗
v = E0 ln S̄T
T
1
Z
δ∗ 2
= ln(λ) + ln(S0 ) + E0 θ(s, S̄s ) ds
2 0
1
U (x) = xγ
γ
for γ < 1, γ 6= 0
U ′ (x) = xγ−1
1
U ′−1 (y) = y γ−1
U ′′ (x) = (γ − 1) xγ−2
concave
1
! γ−1
γ
λ 1
1−γ
1
û(t, Ŝt0 ) = Et = λ γ−1 Et S̄Tδ∗
S̄Tδ∗ S̄T δ∗
1 γ
û(t, Ŝt0 ) S̄tδ∗
1−γ
= λ γ−1
2
γ θ
× Et exp (T − t) + θ (WT − Wt )
1−γ 2
2
θ γ
1 γ
δ∗ 1−γ
= λ γ−1 S̄t exp (T − t)
2 (1 − γ)2
Jtδ̃ = 1 − γ
expected utility
1 θ2 γ
γ
v δ̃ = E0 S̄Tδ̃ = exp − T (S0 )γ
γ 2 1−γ
H
discounted payoff H̄ = 0
ST
(not perfectly hedgable)
!!
δ̃ S̄Tδ̃
vε,V = E0 U (S0 − ε V ) + ε H̄
S0
δ̃ δ̃
vε,V − v0,V
lim =0
ε→0 ε
Davis (1997)
!!
δ̃
′
V S̄Tδ̃
vε,V ≈ E0 U S̄Tδ̃ + εU S̄Tδ̃ H̄ −
S0
!
S̄ δ̃
δ̃
= v0,V + ε E0 U ′ S̄Tδ̃ H̄ − V ε E0 U ′ S̄T δ̃ T
S0
′−1
S̄Tδ̃ = ŜTδ̃ /ŜT0 =U λ ŜT0 =⇒
!!
δ̃
S̄T
δ̃
vε,V δ̃
− v0,V = ε E0 λ ŜT0 H̄ − V E0 λ ŜT0
S0
d Z T
H X
= ÛH (t) + xkH (s) dWsk + MH (t)
STδ∗ k=1 t
MH -orthogonal martingale
k
E MH , W t = 0
fractions
⊤
π δH (t) = bδH (t) ⊤
b−1
t
d
1 X j xkH (t)
bkδH (t) = δH (t) Ŝtj bj,k
t = + θtk .
ÛH (t) j=0
ÛH (t)
1
δ∗
|θt |2 dt + |θt | dŴt
d ln St = rt dt +
2
• volatility
r
d
|θt | = [ln (S δ∗ )]t
dt
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
1982 1992 2002
time
-10 -5 0 5 10
cτ,K (0, S, σ)
maturity τ ∈ [0, T ]
volatility σ
Vc,τ,K (0, S0 )
• implied volatility
call
σBS (0, S0 , τ, K)
call
Vc,τ,K (0, S0 ) = cτ,K 0, S0 , σBS (0, S0 , τ, K)
time 200
0.6
300
0.26
0.24
0.22
0.5
0.75
1
1.25
K S
1.5
• CEV Model
Lewis (2000)
Beckers (1980)
Schroder (1989)
• savings account
dSt0 = r St0 dt
{θt ≥ 0, t ∈ [0, T ]}
• NP volatility
θt = (Stδ∗ )a−1 ψ
exponent a ∈ (−∞, ∞)
stochastic for a 6= 1
• If a = 1 =⇒ BS model
δ∗ 2a−1
2 δ∗ a
dStδ∗ Stδ∗
= rt + St ψ dt + St ψ dŴt
δ∗ 2 (1−a)
Xt = St
2
dXt = 2(1 − a) rt Xt + ψ (1 − a) (3 − 2a) dt
p
+2 ψ (1 − a) Xt dŴt
and dimension
3 − 2a
δ=
1−a
dPθ
= Λθ (T )
dP
S0δ∗
q
X0
Λθ (t) = St0 = St0
Stδ∗ Xt
1
q=
2 (1 − a)
2
dXt = 2 (1 − a) rt Xt + ψ (1 − a)(1 − 2 a) dt
p
+ 2 ψ (1 − a) Xt dŴθ (t)
dimension
1 − 2a
δθ =
1−a
risk neutral
1 − 2a
δθ =
1−a
real world
3 − 2a
δ=
1−a
-2
-4
-1 0 1 2 3
dPθ (ω)
Z Z
Pθ (A) = dPθ (ω) = dP (w)
A A dP (ω)
Z
= Λθ (T ) dP (ω)
A
=⇒
Z
Pθ (Ω) = Λθ (T ) dP (ω) = E0 (Λθ (T ))
Ω
=⇒ Λθ is strict supermartingale =⇒
• contingent claim
Hτ
Hτ = Hτ (Sτδ∗ ), E <∞
Sτδ∗
• fair price
benchmarked price
!
Hτ (Sτδ∗ )
ûHτ (t, Stδ∗ ) = Et
Sτδ∗
Hτ
= ûHτ (τ, Sτδ∗ )
Sτδ∗
τ
∂ ûHτ (s, Ssδ∗ )
Z
δ∗ a
ûHτ (t, Stδ∗ )
= + Ss ψ dŴs
t ∂S δ∗
−2
10
−3
10
−4
10
−10 −5 0 5 10
−2.84
−2.86
−2.88
−2.9
2
−2.92
1.5
−2.94
−5
−4 1
−3
−2.15
−1
0 0.5
Estimated λ 1
2
3
4 0 α
5
λ
• clustering
θt2
∼ N h, θt2 h
2
• discounted NP
Stδ∗
S̄tδ∗ =
St0
=⇒
κ 1
θ02 > 0, γ2
≥ 2
Scott (1987)
Stein & Stein (1991)
Ornstein-Uhlenbeck process
θ0 ∈ ℜ
Wiggins (1987)
Chesney & Scott (1989)
Melino & Turnbull (1990)
d ln(θt ) = k ln(θ̄) − ln(θt ) dt + γ dW̄t
θ0 > 0
1 1 − 12 ν−1 (
1 2
)
( 2 ν) 2 ν y νε
2
p̄ (y) =
θ2 exp −
ε2 Γ( 12 ν) ε2 y
4(ξ−1)
dθt2 = k θt (θ̄ 2 − θt2 ) dt + γ θt2 ξ dW̄t
degree of freedom
2 (2 ξ − 1)
ν= ε2
1− θ̄ 2
ξ=1
two-factor model
0.2 0.4
0.18 0.2
80 0
rho
90
100 -0.2
K 110
-0.4
120
GLMR
Heston
Black-Scholes
-0.10 -0.15
Figure 4.12: Option price differences between the ARCH diffusion (GLMR),
Heston and Black-Scholes models, with time to maturity of three months.
GLMR
Heston
Black-Scholes
0.14
• Diversification Theorem =⇒
• discounted NP
where
d
1 X
dWt = θtk dWtk
|θt | k=1
90
80
70
60
50
40
30
20
10
0
1930 1940 1950 1960 1970 1980 1990 2000
time
-1
-2
1930 1940 1950 1960 1970 1980 1990 2000
time
• discounted NP drift
αt = S̄tδ∗ |θt |2
• volatility
s
αt
|θt | =
S̄tδ∗
fundamental value
• transformed time
t
1
Z
ϕt = αs ds
4 0
Xϕt = S̄tδ∗
√
dW (ϕt ) = αt dWt
p
dXϕ = 4 dϕ + 2 Xϕ dW (ϕ)
3 αt 1√
q
d S̄tδ∗ = q dt + αt dWt
8 S̄tδ∗ 2
• quadratic variation
Z t
hp i 1
S̄ δ∗ = αs ds
t 4 0
• transformed time hp i
ϕt = S̄ δ∗
t
observable
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 5.3: Empirical quadratic variation of the square root of the discounted
NP.
c Copyright Eckhard Platen 10 BA - Toronto C5 162
Stylized Minimal Market Model
αt = α exp {η t}
t
α
Z
ϕt = exp {η z} dz
4 0
α
= (exp {η t} − 1)
4η
0
1930 1940 1950 1960 1970 1980 1990 2000
time
S̄tδ∗
Yt =
αt
p
dYt = (1 − η Yt ) dt + Yt dWt
• discounted NP
S̄tδ∗ = Yt αt
• NP
Stδ∗ = St0 S̄tδ∗ = St0 Yt αt
• speed of adjustment η
ln(2)
• half life time of major displacement η
≈ 13 years
80
70
60
50
40
30
20
10
0
1930 1940 1950 1960 1970 1980 1990 2000
time
1
|θt | = √
Yt
• squared volatility
1
2 2 2
32
d|θt | = d = |θt | η dt − |θt | dWt
Yt
0.3
0.25
0.2
0.15
0.1
1930 1940 1950 1960 1970 1980 1990 2000
time
12
1 y x+y
p(s, x; t, y) = exp −
2 (ϕt − ϕs ) x 2 (ϕt − ϕs )
√
xy
× I1
ϕt − ϕ s
α
ϕt = (exp{η t} − 1)
4η
I1 (·) modified Bessel function of the first kind
y S̄tδ∗
=
ϕt − ϕ s ϕt − ϕ s
x S̄sδ∗
=
ϕt − ϕ s ϕt − ϕ s
100
0.1
y 110
120
P (0, T ) = 0.00077
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 5.8: Savings bond, fair zero coupon bond and savings account.
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 5.9: Logarithms of savings bond and fair zero coupon bond.
• forward rate
∂
f (t, T ) = − ln(P (t, T ))
∂T
= rT + n(t, T )
( )!
∂ S̄tδ∗
n(t, T ) = − ln 1 − exp −
∂T 2(ϕT − ϕt )
1 S̄tδ∗ αT
= 2
S̄tδ∗
n o
exp −1 (ϕ T − ϕ t ) 8
2(ϕT −ϕt )
lim n(t, T ) = η
T →∞
lower bound
Stδ ≥ −P (0, T ) exp{r t}
=⇒
there is a free lunch with vanishing risk
• fair zero coupon bond has a lower price than the savings bond
( Z )
T 0
S
P (t, T ) < PT∗ (t) = exp − rs ds = 0 t
t ST
• Radon-Nikodym derivative
Ŝt0 S̄0δ∗
Λt = =
Ŝ00 S̄tδ∗
1.4
1.2
0.8
0.6
0.4
0.2
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 5.11: Radon-Nikodym derivative and total mass of putative risk neu-
tral measure.
c Copyright Eckhard Platen 10 BA - Toronto C5 185
• hypothetical risk neutral measure
total mass:
( )
S̄0δ∗
Pθ,T (Ω) = E0 (ΛT ) = 1 − exp − < Λ0 = 1
2 ϕT
!
(STδ∗ − K) +
cT ,K (t, Stδ∗ ) = Stδ∗ Et
STδ∗
!+
K Stδ∗
= Et Stδ∗ −
STδ∗
2
Stδ∗
= 1 − χ (d1 ; 4, ℓ2 )
for t ∈ [0, T )
European calls
European puts
barrier options
local martingale
=⇒
c∗T ,K (t, Stδ∗ ) = cT ,K (t, Stδ∗ )
for t ∈ [0, T )
a.s. St0
= K >0
ST0
dramatic differences can arise
Figure 5.13: Logarithms of savings bond times K, risk neutral put and fair
put.
c Copyright Eckhard Platen 10 BA - Toronto C5 199
MMM with Random Market Activity
• volatility
s
αt
|θt | =
S̄tδ∗
• discounted NP drift
αt to be modeled
0.26
0.24
0.22
0.5
0.75
1
1.25
K S
1.5
time 200
0.6
300
Figure 5.16: Implied volatilities for S&P500 three month options.
time 200
0.6
300
Figure 5.17: Implied volatilities for S&P500 one year options.
p
2
dmt = k(mt ) β dt + β mt ̺ dWt + 1 − ̺2 dW̃t
̺ - correlation parameter
• drift
mt
k(mt ) = (p − g mt )
2
• stationary density
g p−1
pm (y) = y p−2 exp{−g y}
Γ(p − 1)
p−1 1
gamma density with mean g
and variance g
for p > 1 and g > 0
1 40
0 30
0.5 g
20
1
10
y
1.5
PT (t, S̄tδ∗ , mt ) = Stδ∗ P̂T (t, S̄tδ∗ , mt ) = St0 S̄tδ∗ P̂T (t, S̄tδ∗ , mt )
Figure 5.20: Implied volatilities for put options as a function of strike K and
speed of adjustment g.
c Copyright Eckhard Platen 10 BA - Toronto C5 212
0.4
0.375 10
0.35
8
0.325
0.3 6
80 T
4
100
K 2
120
Figure 5.21: Implied volatilities for long dated put options as a function of
strike K and maturity T .
Pl. (2001)
Heath & Pl. (2005a)
d + 1 currencies
• ith normalized NP
q d+1
X
dYti = 1 − η i Yt i i
q i,k dWtk
dt + Yt
k=1
!
d+1 i,k j,k i,k
X q q q
dXti,j = Xti,j (r i − r j ) dt + p −q p dt + dWtk
k=1 Yti Ytj Yti
!
d+1
X q i,k q j,k q i,k
dSij (t) = Sij (t) r i dt + p −q p dt + dWtk
i i
k=1 Y t Yt j Y t
i,k
q
θik (t) = p
Yti
bj,k
i (t) = θi
k
(t) − θ k
j (t)
equity prices
commodity prices
roll-ups:
original investment accrued at a pre-defined interest rate
ratchets:
death benefit based upon the highest anniversary account value
embedded options:
hedging against market downturn occurrence of death
IFRS Phase II
CFO-Forum (2008)
CRO-Forum (2008)
Solvency II
Verheugen & Hines (2008)
benchmark approach in
Wt - Brownian motion
• savings account
dBt = rt Bt dt
µt − γ − rt
θt =
σt
• risky asset
dSt
= rt dt + σt θt dt + σt dWt
St
Vt = δt0 Bt + δt1 St
• fractions
0 Bt 1 St
πt0 = δt , πt1 = δt
Vt Vt
πt0 + πt1 = 1
Long (1990)
V ∗ is best performing in several ways
Kelly (1956)
V ∗ = max E(log VT )
µt − γ − rt θt
πt1∗ = =
σt2 σt
dVt∗
= rt dt + θt (θt dt + dWt )
Vt∗
Merton (1992)
1 VT∗ 1 VT
lim sup log ≥ lim sup log
T →∞ T V0∗ T →∞ T V0
Pl. (2005b)
90
80
70
60
50
40
30
20
10
0
1930 1940 1950 1960 1970 1980 1990 2000
time
-1
-2
1930 1940 1950 1960 1970 1980 1990 2000
time
Ût ≥ Et Ûs
t≤s
Pl. (2002)
no strong arbitrage
Ût = Et Ûs
t≤s
Pl. (2008)
VT = VT′
supermartingale property
=⇒
Vt ≤ Vt′
t ∈ [0, T ],
0.12
0.1
0.08
0.06
0.04
0.02
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 7.3: Benchmarked savings bond and benchmarked fair zero coupon
bond.
c Copyright Eckhard Platen 10 BA - Toronto C7 233
• real world pricing formula
for
HT
Et <∞
VT∗
HT
UH (t) = Vt∗ Et
VT∗
t ∈ [0, T ]
=⇒
UH (t) = P (t, T ) Et (HT )
complete market
dQ Bt V0∗
Λt = =
dP At B0 Vt∗
supermartingale =⇒
1 = Λ0 ≥ E0 (ΛT )
B0
UH (0) = EQ HT
BT
NP as Long (1990)
P (t, T ) 1
P̂ (t, T ) = = Et
Vt∗ VT∗
martingale (no trend)
for rt - deterministic
T !
1 V̄t∗
Z
P (t, T ) = Vt∗ Et ∗
= exp − rs ds Et ,
VT V̄T∗
t
Vt∗
V̄t∗ = Bt
- discounted NP
0.12
0.1
0.08
0.06
0.04
0.02
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 7.4: Benchmarked savings bond and benchmarked fair zero coupon
bond.
c Copyright Eckhard Platen 10 BA - Toronto C7 241
=⇒ downward trend reflects equity premium
T
Bt
Z
P (t, T ) < exp − rs ds =
BT
t
1.4
1.2
0.8
0.6
0.4
0.2
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 7.5: Radon-Nikodym derivative and total mass of putative risk neutral
measure.
c Copyright Eckhard Platen 10 BA - Toronto C7 243
1
savings bond
fair zero coupon bond
0.9 savings account
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 7.6: Savings bond, fair zero coupon bond and savings account.
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 7.7: Logarithms of savings bond and fair zero coupon bond.
-1
-2
1930 1940 1950 1960 1970 1980 1990 2000
time
=⇒ volatility
αt
r
θt =
V̄t∗
reflects leverage effect
MMM
discounted NP drift
assume
αt = α0 exp{ηt}
α0 > 0
=⇒ MMM
-1
-2
1930 1940 1950 1960 1970 1980 1990 2000
time
V̄t∗
Yt =
αt
p
dYt = (1 − ηYt ) dt + Yt dWt
1 αt
r
θt = √ =
Yt V̄t∗
rt - deterministic
T ( )!
∗
Z V̄t
P (t, T ) = exp − rs ds 1 − exp −
2(ϕ(T ) − ϕ(t))
t
Explicit formula !
-1
-2
-3
-4
-5
-6
-7
-8
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 7.10: Logarithms of savings bond and fair zero coupon bond.
0.8
0.6
0.4
0.2
0
1930 1940 1950 1960 1970 1980 1990 2000
time
6e-005
5e-005
4e-005
3e-005
2e-005
1e-005
-1e-005
-2e-005
1930 1940 1950 1960 1970 1980 1990 2000
time
! ( )!
ΛT V̄t∗ V̄t∗
Et = Et = 1 − exp −
Λt V̄T∗ 2(ϕ(T ) − ϕ(t))
Bt
=⇒ P (t, T ) < BT
1.4
1.2
0.8
0.6
0.4
0.2
0
1930 1940 1950 1960 1970 1980 1990 2000
time
Figure 7.13: Radon-Nikodym derivative and total mass of putative risk neu-
tral measure.
c Copyright Eckhard Platen 10 BA - Toronto C7 258
• European put option under the MMM
!
(K − Vτ∗ )+
p(t, Vt∗ , T, K, r) = Vt∗ Et
Vτ∗
−r(T −t) 2
+ Ke χ (d1 ; 0, l2 ) − exp {−l2 /2}
and
4ηVt∗
l2 =
Bt αt (exp{η(T − t)} − 1)
n ≥ 0 degrees of freedom
∞
exp − 2l l k Γ x2 ; n+2k
X 2
χ2 (x; n, l) = 2 1 −
k=0
k! Γ n+2k2
and
p̃(t, Vt∗ , T, K, r) → K e−r(T −t) > 0
Figure 7.14: Logarithms of savings bond times K, risk neutral put and fair
put.
c Copyright Eckhard Platen 10 BA - Toronto C7 263
Guaranteed Minimum Death Benefit (GMDB)
max(egτ V0 , Vτ )
time of death τ
Vt = e−ξt Vt∗
=⇒
=⇒ payoff
h i
HT = GM DBT = e−ξT (e(g+ξ)T V0∗ − VT∗ )+ + VT∗
GM DBT
GM DB0 = V0∗ E
VT∗
= e−ξT p(0, V0∗ , T, e(g+ξ)T V0∗ , r) + V0∗
GMDB0
0
1 1 1
GMDB
g=0
0.95 0.95 0.95
GMDB0
0
GMDB
1 1 1
Figure 7.15: Present value of the GMDB under the real world pricing formula (left),
the risk neutral pricing formula (middle) and the Black Scholes formula (right) for
η = 0.05, α0 = 0.05, r = 0.05, ξ = 0.01 and Y0 = 20.
c Copyright Eckhard Platen 10 BA - Toronto C7 268
• lifetime τ is stochastic
GM DBτ
GM DB0 = V0∗ E E Fτ
Vτ∗
Z T
GM DB0 = p(0, V0∗ , t, e(g+ξ)t V0∗ , r) + V0∗ e−ξt fτ (t) dt
0
male
−1
female
10
−2
10
log(qx)
−3
10
−4
10
−5
10
0 20 40 60 80 100 120
Age x
surrender charge βt
(8 − ⌈t⌉)%, t ≤ 7,
βt =
0, t>7
no credit risk
no accumulation phase
male
female
1.1
GMDB0
1.05
1
20 30 40 50 60 70 80
Age x
Figure 7.17: Value of the GMDB under the MMM for male and female pol-
icyholders aged x, assuming an irrational lapsation of l = 1%.
m ∈ {1, 2, . . . , d}
hkt > 0
Z t
hks ds < ∞
0
• jump martingale
k − 12
dqtk = dpkt − hkt dt ht
k ∈ {1, 2, . . . , d−m}
• trading uncertainty
d
!
X
dStj = St−
j
ajt dt + bj,k
t dWt
k
k=1
j ∈ {1, 2, . . . , d}
at = (a1t , . . . , adt )⊤
1 = (1, . . . , 1)⊤
=⇒
d
!
X
dStj = St−
j
rt dt + bj,k
t (θt
k
dt + dW k
t )
k=1
j ∈ {0, 1, . . . , d}
• portfolio
d
X
Stδ = δtj Stj
j=0
• self-financing if
d
X
dStδ = δtj dStj
j=0
• fraction
j Stj
πδ,t = δtj
Stδ
1 d ⊤
π δ,t = (πδ,t , . . . , πδ,t )
n o
⊤
dStδ = δ
St− rt dt + π δ,t− bt (θ t dt + dW t )
strictly positive ⇐⇒
d
X q
j
πδ,t bj,k
t >− hk−m
t
j=1
k ∈ {m + 1, . . . , d}
k=1 j=1
d d
bj,k
q
t
X X j
+ ln1 + πδ,t− q hk−m
t dWtk
k=m+1 j=1 hk−m
t
• growth rate
2
m d d
X X j 1 X
gtδ = rt + πδ,t bj,k
t θtk − j
πδ,t bj,k
t
2
k=1 j =1 j =1
d d d j,k
b
X X q X j
j j,k k k−m hk−m
+ πδ,t bt θt − ht +ln1+ πδ,t q t t
k=m+1 j =1 j =1 hk−m
t
k ∈ {m + 1, . . . , d}
⊤
π δ∗ ,t = (πδ1∗ ,t , . . . , πδd∗ ,t )⊤ = c⊤
t b−1
t
θtk for k ∈ {1, 2, . . . , m}
ckt =
θtk
k−m − 1
for k ∈ {m + 1, . . . , d}
1−θtk (ht ) 2
θtk
if q
k−m
≪1
ht
=⇒ ckt ≈ θtk
dStδ∗ = δ∗
St− rt dt + c⊤
t (θ t dt + dW t )
m
X
δ∗
= St− rt dt + θtk (θtk dt + dWtk )
k=1
d
X θtk
+ (θtk dt + dWtk )
−1
k=m+1 1 − θtk (hk−m
t ) 2
=⇒
S δ∗ is a NP.
=⇒
m d
1 X X θtk θtk
gtδ∗ = rt + (θtk )2 − hk−m
t
ln1 +q + q
2 k=1 k=m+1 hk−m − θtk hk−m
t t
d
1 X |θ t |2
gtδ∗ ≈ rt + (θtk )2 = rt +
2 k=1
2
d
X
d ln(Stδ∗ ) ≈ gtδ∗ dt + θtk dWtk
k=1
d
!
X
dStδ∗ Stδ∗ θtk θtk dWtk
≈ rt + dt +
k=1
• benchmarked portfolio
Stδ
Ŝtδ =
Stδ∗
m
X d
X
dŜtδ = δtj Ŝtj bj,k
t − Ŝtδ θtk dWtk
k=1 j=1
d d
X X θtk
+ δtj j
Ŝt− bj,k
t
1 −q − Ŝ δ θ kdW k
t− t t
k=m+1 j=1 hk−m
t
driftless
d
X
dŜt0 = −Ŝt−
0
θtk dWtk
k=1
driftless
Corollary 8.6
A JDM does not allow nonnegative portfolios that permit strong arbitrage.
H
Corollary 8.7 Aτ -measurable payoff H with E( |A0 ) < ∞. If
Sτδ∗
there exists a fair nonnegative portfolio S δ that replicates the payoff
Sτδ = H,
=⇒
=⇒
Ŝt0
derivative Ŝ00
forms (A, P )-martingale
d
X
dP̂ (t, T ) = −P̂ (t−, T ) σ k (t, T ) dWtk
k=1
d
! !
t t k
σ (s, T )
X Z q Z
+ σ k (s, T ) hk−m
s ds + ln 1 − p dpks
k=m+1 0 0 hk−m
s
=⇒
• forward rate equation
m Z t
∂
X
k k
dWsk
f (t, T ) = f (0, T ) + σ (s, T ) σ (s, T ) ds +
k=1 0 ∂T
d t
1 ∂
X Z
k k
dWsk
+ σ k (s,T )
σ (s, T ) σ (s, T ) ds +
k=m+1 0 1− √ ∂T
k−m
hs
δ
Sequence (S(d) )d∈N is a sequence of DPs if
j K2
πδ,t ≤ 1
d 2 +K1
d
X
j j j,k
dŜ(d) (t) = −Ŝ(d) (t−) σ(d) (t) dWtk
k=1
d
X
k j,k
σ̂(d) (t) = |σ(d) (t)|
j=0
2
k
E σ̂(d) (t) ≤ K5
• tracking rate
2
d
X d
X
δ j j,k
R(d) (t) = πδ,t σ(d) (t)
k=1 j=0
δ δ
S(d) is NP ⇐⇒ Ŝ(d) ≡ constant
⇐⇒
δ
R(d) (t) = 0
δ P
lim R(d) (t) = 0
d→∞
δ
(S(d) )d∈N has vanishing expected tracking rate, if
=⇒
Platen (2005a)
JD δ
For a regular sequence of JDMs (S(d) )d∈N , each sequence (S(d) )d∈N
(K2 )2 K5
eδ(d) (t) ≤ .
d2K1
• savings account
0
S(d) (t) = exp{r t}
• discounted NP drift
αδt ∗ = α0 exp{η t}
• NP
0
δ∗
S(d) (t)
S(d) (t) = 0
Ŝ(d) (t)
j j δ∗
S(d) (t) = Ŝ(d) (t) S(d) (t)
400
350
300
250
200
150
100
50
0
0 5 9 14 18 23 27 32
50
40
30
20
10
0
0 5 9 14 18 23 27 32
EWI
90 GOP
80
70
60
50
40
30
20
10
0
0 5 9 14 18 23 27 32
Market index
90
GOP
80
70
60
50
40
30
20
10
0
0 5 9 14 18 23 27 32
Specifying a Continuous NP
θtk = 0
NP !
m
X
dStδ∗ Stδ∗ θtk θtk dWtk
= rt dt + dt +
k=1
S0δ∗ = 1
d
X
dŜtj = −Ŝt−
j
σtj,k dWtk
k=1
( m m
)
t X t
1
Z Z
j,k 2
X
Ŝtj S0j σsj,k dWsk
= exp − σs ds −
2 0 k=1 k=1 0
k−m
Z t d q d pt στj,k
k
l −
X Y Y
× exp σsj,k hk−m
s ds 1 − q
0 k=m+1
k=m+1 l=1 hk−m
τ k−
l
k=1
m Z t
X σsj,k
Ŵtj = dWsk
k=1 0 |σsj |
bj,k
t = θtk − σtj,k
for k ∈ {1, 2, . . . , m}
bj,k
t = −σtj,k
k ∈ {m + 1, . . . , d}
√
σtj,k =σ j,k
≤ hk−m
t t
1
Z Z
Ŝtj,c = S0j exp − |σsj |2 ds − |σsj | dŴsj
2 0 0
!pk−m
d
X √ d
Y σ j,k
t
rt = r, σtj,k = σ j,k
1
Ŝtj,c = S0j exp − |σ j |2 t − |σ j | Ŵtj
2
• discounted NP drift
q
dYtj = 1 − η j Ytj dt + Ytj dŴtj
• continuous part
1
Ŝtj,c =
αj (t)Ytj
j 1
Xϕ j (t) =α j
(t)Ytj =
Ŝtj,c
! ( ) !
T
1 1
Z
P (t, T ) = Stδ∗ E At = E exp − rs ds ŜT0 At
STδ∗ Ŝt0 t
• MM case
1
P (t, T ) = exp{−r(T −t)} E ŜT0 At = exp{−r(T −t)}
Ŝt0
1
λjt =
Ŝtj (ϕj (t) − ϕj (T ))
ŜT0 = ŜT0,c
( ) !
T
1
Z
P (t, T ) = E exp − rs ds At E ŜT0 At
t Ŝt0
( ) !
T
1
Z
= E exp − rs ds At 1 − exp − λ0t
t 2
!
j
F (t, T ) − STj
Stδ∗ E At =0
STδ∗
forward price
j
S
j
Stδ∗ E ŜTj At t
P (t,T 1
) Ŝ j
E ŜT At if Stj > 0
j t
F (t, T ) = =
Stδ∗ E 1
δ∗ At
0 if Stj = 0
ST
• MMM case
1 1 1
E ŜTj At = E ŜTj,c At E STj,d At
Ŝtj Ŝtj,c Stj,d
1
= 1 − exp − λjt
2
forward price
n o !!−1
1 j
1 − exp − 2 λt
( Z )
T
F j (t, T ) = Stj 1 0 E exp − rs ds At
1 − exp − 2 λt t
Stj
= E 1{Ŝ j ≥K(S 0 )−1 Ŝ 0 } ŜTj At
Ŝtj T T T
Stj
= E 1{Ŝ j,c ≥g(pk−m −pk−m )Ŝ 0 }
Ŝtj,c T T t T
!pk−m k−m
−pt
n √ o σ j,k T
∞ √
X k−m (hk (T − t))n n
j,k
o
= exp −h (T − t) exp σ hk−m (T − t)
n=0
n!
!n
σ j,k
Stj
× 1− √ E 1{Ŝ j,c ≥g(n)Ŝ 0 } ŜTj,c At
hk−m Ŝtj,c T T
for all n ∈ N
∞
j,k
X k−m (hk−m (T − t))n
A (t, T, K) = exp −h (T − t)
n=0
n!
!n
n √ o σ j,k
× exp σj,k hk−m (T − t) 1 − √ Stj N (d1 (n))
hk−m
σ j,k
Stj
√ ln 1− √
j,k hk−m 1 0,j 2
ln K
+r + σ hk−m + n T −t
+ 2
σ̂ (T − t)
d1 (n) = √
σ̂ 0,j T −t
̺i,j - correlation
∞
X (hk−m (T − t))n
Aj,k (t, T, K) = exp −hk−m (T − t)
n=0
n!
!n
n √ o σ j,k
j,k
× exp σ hk−m (T − t) 1− √
hk−m
ϕ0 (T ) − ϕ0 (t)
1 j
j ′′ j 0
× St G0,4 ; λt , λt − exp − λt
g(n) 2
MM case
∞
X (hk−m (T − t))n
B j,k (t, T, K) = exp{−hk−m (T − t)}
n=0
n!
B j,k (t, T, K)
∞
X (hk−m (T − t))n
= exp{−hk−m (T − t)}
n=0
n!
× K exp{−r(T − t)} 1 − G′′
0,4 (ϕ j
(T ) − ϕ j
(t))g(n); λ0
t , λ j
t
+
j
ST −K
j,k
cT ,K (t) = Stδ∗ E At
STδ∗
!
STj −K
= Stδ∗ E 1{S j ≥K} At
T
STδ∗
∞
X (hk−m (T − t))n
cj,k
T ,K (t) = exp{−hk−m (T − t)}
n=0
n!
!n
n √ o σ j,k
× exp σ j,k
hk−m (T − t) 1 − √ Stj N (d1 (n))
hk−m
!
− K exp{−r(T − t)}N (d2 (n))
∞
X (hk−m (T − t))n
cj,k
T ,K (t) = exp{−h k−m
(T − t)}
n=0
n!
" !n
n √ o σ j,k
× exp σ j,k hk−m (T − t) 1 − √
hk−m
0 0 !
ϕ (T ) − ϕ (t) 1
× Stj G′′
0,4 ; λjt , λ0t − exp − λjt
g(n) 2
#
− K exp{−r(T − t)} 1 − G′′ j j 0 j
0,4 (ϕ (T ) − ϕ (t))g(n); λt , λt
maturity T
zero recovery
!
1{τ k−m >T }
k−m
P̃ (t, T ) = Stδ∗ E 1
At
STδ∗
!
1
= Stδ∗ E At E 1{τ k−m >T } At
STδ∗ 1
= P (t, T )P (pk−m
T = 0 At )
P pk−m
T = 0 At = E 1{pk−m =0} 1{pk−m −pk−m =0} At
t T t
= 1{pk−m =0} P pk−m
T − pk−m
t = 0 At
t
( ) !
Z T
= 1{pk−m =0} E exp − hk−m
s ds At
t
t