Lecture 8 - SABR Model
Lecture 8 - SABR Model
Lecture 8 - SABR Model
Daniel Brødsgaard
14 November 2022
Overview of today
▶ In a casino you realise that the roulette is biased, so that red comes
up 70% of the time.
▶ There is no green, ie. the casino earns nothing in expectation.
▶ The odds are 50/50.
▶ This means that if you bet 50 on black or red and you colour comes
up, you gain 100.
▶ You decide to start a derivatives market where you, the market
maker, quote prices on tickets that pay 100 if red comes up, zero if
not.
▶ Q: Quote a two-way price for the tickets.
▶ Task illustrates again the difference between the P and Q measures:
speculators vs hedgers.
Motivation
"Z #
Z T T
1
Γ(u)S 2 (u) σr2 (u) − σBS
2
Et [ PnL(u)du] = Et du
t 2 t
Z + inf
C (K , T ) = (x − K )+ ϕT (x)dx
− inf
Z + inf 2
∂ 2 C (K , T ) ∂ (x − K )+
2
= ϕT (x)dx
∂K − inf ∂K 2
Z + inf
= δ(K )ϕT (x)dx = ϕT (K )
− inf
From the last lecture we know how we can use local vols to generate
implied BS volatilities:
hR i
T
Et t
Γ(u)S 2 (u)σ 2 (S, u)du
2
σBS = hR
T
i
Et t Γ(u)S 2 (u)du
RT R
t
ϕ(S, u)Γ(u)S 2 (u)σ 2 (S, u)dSdu
= RT R
t
ϕ(S, u)Γ(u)S 2 (u)dSdu
Figure 4: Left: Weight large due to density. Right: Weight large due to
gamma. Source: Antoine Savine.
Some useful remarks on local vol
Bruno Dupire (1992) showed that if the local volatility function is
2 2∂
C (K ,T )
∂K 2
σ (K , T ) = ∂C (K ,T )
K 2
∂T
then the local volatility model will be able to fit all European option
prices of all strike as well as maturities. Furthermore, two useful
approximations follow in the case where local vol is either independent of
time or describing short maturities:
σBS (K ) ≈ σ(S(t))
K =S(t)
∂σBS (K ) 1 ∂σ(S(t))
≈
∂K K =S(t) 2 ∂S
Ie. implied ATM vol is equal to local vol, and the implied ATM vol skew
(in the strike dimension) is roughly half of the local vol skew (in the spot
dimension).
Some drawbacks of local vol
▶ Given the formulas above we can infer some dynamics. Assume for
an example that 3m expiry options are currently pricing a smile as
shown below (blue line). Note that the local vol function must be a
more convex function in this case, since implied vol is approximately
a weighted average of local vols between spot and strike.
▶ In this case a local vol model predicts that implied ATM vol is
assumed to increase no matter what direction the underlying moves!
▶ Furthermore, if the underling declines then the ATM skew increases
and vice versa. This is not in line with how the market trades.
▶ ATM vol is typically higher with lower swap rates. Note, however,
that normal volatilities are closer to constant since σnormal ≈ σBS S
▶ The market inhibits a skew/smile at all times, although it can be
steeper/flatter.
Stylized facts on smile/skew: Backbone and skew effect
Figure 7: Implied vol for given strike.
dC 1 d 2C 2 1 d 2C 2 d 2C
dt = − (dS) − (dσ) − (dSdσ)
dt 2 dS 2 2 |dσ 2
{z } |dσdS
{z }
volga vanna
where we have assumed that the quadratic variation terms: (dS)2 and
(dσ)2 , are functions of dt only.
▶ Thus, we have two new greeks, vanna and volga, to take into
account, which the BS model was not able to.
▶ Signs depends on the vol-of-vol (dσ)2 and on correlation between
spot and vol (dSdσ) besides the original gamma which depends on
the vol of spot (dS)2 .
What does volga and vanna look like?
Figure 8: Top-L: Volga. Top-R: Vanna. Low-L: volga
vega . Low-R:
vanna
vega .
Functions of strike (S=100). See appendix for intuition.
More on vanna and volga
We will see now why we above represented volga and vanna in ratios
relative to vega as well.
▶ We have seen that with stochastic vol there are two more greeks to
take into account.
▶ Thus, in a stochastic volatility model we can, relative to the
constant vol in Black-Sholes, write the excess value (E ) for a strike
K as (where SV (K ) the stochastic vol value of the option)
dC 1 d 2C 1 d 2C d 2C
dt = − 2
(dS)2 − 2
(dσ)2 − (dSdσ)
dt 2 dS |2 dσ
{z } |dσdS
{z }
volga vanna
Thus, we have a local vol function of the form σloc (f , t) = σt ftα−1 From
our local vol approximations we know that
∂σBS
K =f (0)
σBS ≈ σ0 f0α−1 and ≈ (α − 1)σ0 f0α−2
K =f (0) ∂f0
∂σBS
K =f (0) σBS (ATM)
=⇒ = (α − 1) <0 if α < 1
∂f0 f0
Thus, if ATM rates move higher the backbone effect implies a lower
implied vol for the new ATM vol when α < 1.
How α controls backbone continued
The effect on risk from a choice of α can be seen by noting that for an
option with strike K , changing the ATM forward rate
∂σBS (K , f0 )
σBS (K , f0 ) ≈ σBS (f0 , f0 ) + (K − f0 )
| {z } ∂K
ATM imp. vol
| {z }
Skew given ATM vol
∂σBS (K , f0 ) ∂σBS (f0 , f0 ) ∂σBS (K , f0 )
=⇒ ≈ −
∂f0 ∂f0 ∂K
| {z } | {z }
backbone skew at strike K
last line follows from differentiating with respect to the ATM rate f0 .
▶ So the change in implied vol for a given strike can be decomposed
into a backbone effect and a skew effect.
▶ If for an example α is close to zero and skew at strike K is close to
flat, then we should expect implied vol to drop at strike K, when
ATMF increases.
▶ If on the other hand the skew at strike K is is highly negative, then
we cannot a priori say if implied vol decreases or increases.
α effect on change in implied vol
Figure 16: Change in smile for different α. Blue line is initial smile. Red
and grey lines after +100bp bump.
SABR effect on delta
▶ When talking greeks in the SABR model one can either differentiate
directly wrt. to the SABR parameters or alternatively adjust the
greeks in the BS model. We will do the latter as is market practice.
▶ Implied vol in the BS model is now a function of the spot, meaning
that for delta we get
Figure 17: Left: SABR vs BS delta. Right: Smile after underlying being
bumped more than 100bp.
A closed form solution still exists for european options in the BS model in
this framework since implied vol can be found to equal
s
Z T
1
σBS = σ 2 (u)du
T t
Appendix: Local vol w. time dependance continued
Pj
2 i=1 σi2 · (Ti − Ti−1 )
σBS (Tj ) =
Tj
2 2
σBS (Tj ) · Tj − σBS (Tj−1 ) · Tj−1 = σj · (Tj − Tj−1 )
2 2
σBS (Tj ) · Tj − σBS (Tj−1 ) · Tj−1
=⇒ σj2 =
(Tj − Tj−1 )
Appendix: Other drawbacks of local vol
▶ Unless very simply specified (as the CEV case) the local vol models
a la Dupire requires an entire continuous surface of options prices in
both strike and expiry dimension.
▶ This makes the model somewhat difficult to calibrate, since first we
need to interpolate in both strike and expiry dimension (not as
simple as you might think - we need to keep the model arbitrage
free).
▶ Somewhat of a black box since non-parametric (Dupire calibration).
Lacking intuition as to what factors are relevant for implied vol.
▶ We need frequent re-calibration since, as we have seen, the local
volatility model implies implied vol moves which are different,
sometimes even opposite, of empirics. Thus, the local volatility
model quickly ”moves away” from markets.
▶ Instability of the model is bad for traders and hedge ratios as they
will first of all be wrong, and secondly change often.
Appendix: A simple local vol specification: CEV
A simple example of a local vol model is the CEV (contant elasticity of
volatility) model, which is specified as
Thus, the CEV specification has a local vol function of the form
σ(S; λ, α) = λS(t)α−1 . Following our approximations above
=⇒ λ = σBS (ATM)S(t)1−α
In order to hit ATM vol. Further we have that this model will only be
able to generate a downward sloping skew. Flat if α = 1 (BS) and steep
if α = 0 (Bachelier/Normal). Distribution is chi-suared for α ∈ (0, 1). No
close-form solution exist (numerical implementation needed but can be
approximated).
Appendix: The CEV model in charts
Figure 19: Left: Model vs. market. Right: Implied skew for different α.
Appendix: A note on the difficulty of calibration and jumps
▶ Not all models are able to calibrate to all smiles/skews. Opposite to
local volatility models stochastic vol models cannot be guaranteed
to hit all implied volatilities for every expiry. This is especially the
case for short expiries, where the market can trade with a highly
pronounced smile to which the parameters of volatility models is
unable to calibrate.
▶ One reason for a highly pronounced skew/smile for short expiry
options is the existence of jumps.
▶ Jumps are per definition a dis-continuity and renders everything
we’ve done so far wrong. In the case of jumps we won’t be able to
re-adjust our delta hedges during the jump. In other words: We will
have unhedged risks on our books as the underlying moves
(unexplained PnL).
▶ Overall the existence of jumps should shift upwards all implied vols
at each expiry since jumps will increase average realized vol.
▶ However, it should lift wings more than ATM since jumps towards
the strike level of an option is valuable due to gamma being largest
close to the strike. As OTM or ITM options are more likely to
observe jumps towards them this can explain highly pronounced
Appendix: Vanna and volga intuition
Volga
▶ Volga is intuitively zero ATM. Why? Implied vol increases when the
market expects a higher average volatility to prevail over some time
period. That is a positive for all options. However, it also changes
the market implied probability distribution at the same time.
▶ When implied vol increases, for at the money options the probability
that the spot is gonna be fairly close to the current spot (ATM) for
at least some time is close to unchanged.
▶ However, for OTM/ITM options the market implied probability that
the spot will at some future point in time be close to the strike level
of these options increases. Thereby it increases the buyer of the
option’s ability to monetize moves in the underlying at the same
time as the market expects a higher average vol. Thus, owners of
OTM/ITM options are exactly exposed towards the high implied vol
when it increases and vice versa!
▶ For VERY OTM and ITM the market implied probability changes by
a very small amount.
Appendix: Vanna and volga intuition continued
Vanna
▶ Vanna is the change in delta given a change in implied vol. Vanna is
positive for high strike options (OTM payers and ITM receivers).
Why? When implied vol increases it increases the probability of
exercise, and thereby results in a more positive delta, for OTM
payers. However, it lowers the probability of exercise for ITM
receivers, and thereby results in a less negative delta.
▶ ATM vanna is close to zero since the probability of exercise doesn’t
change for ATM options, given the higher implied vol gives rise to
symmetric changes in the market implied density.
▶ Vanna is negative for low strike options (ITM payers and OTM
receivers). Why? When implied vol increases the probability of
exercise increases for OTM receivers (more negative delta), while the
probability of exercise decreases for ITM payers (less positive delta).
Skew and backbone decomposed in a drawing