Closed-form approximations for operational VaR
Lorenzo Hernández[2], Jorge Tejero[2], Alberto Suárez[1,2,3] , Santiago Carrillo-Menéndez[2,3,4]
[1] Computer Science Department, Universidad Autónoma de Madrid [2] Quantitative Risk Research (QRR) [3] RiskLab Madrid[4] Mathematics Department, Universidad Autónoma de Madrid
2
Extreme events in finance
Sums of heavy-tailed random variables
Rare events Large impact Difficult to predict their specific characteristics
(but can be explained after the fact). Amenable to modeling
Modeling of extreme events is possibleand provides insight
3
Basel II
http://www.bis.org/publ/bcbs128.htmSums of heavy-tailed random variables
4
The three pillars approach
First pillar: Minimum capital requirements(quantification of risk)
• Specifies the guiding principles for the estimation of regulatory/economic capital.
• Operational risk is included as a new type of risk.
Second pillar: Supervisory review process.
Third pillar: Market discipline (+ public disclosure)
Sums of heavy-tailed random variables
5
Operational Risk: Definition
[source: Basel II]644. Operational risk is defined as the risk of
loss resulting from inadequate or failedinternal processes, people and systemsor from external events. This definition includes legal risk, but excludes strategic and reputational risk.
Sums of heavy-tailed random variables
6
Operational Risk: Measurement
Sums of heavy-tailed random variables
7
Loss distribution approach
Model the distribution of the aggregate losses for a given business line & risk type
Calculate Capital at Risk (99.9% percentile) of the aggregate loss distribution per business line & risk typeand add them.
. risk type and line businessfor year in losses ofnumber theis
;
],[1
],[],[],[
jitN
LLoss
jit
N
n
jint
jit
jit
8
1
7
1
],[
i j
jiCaRCaR
Sums of heavy-tailed random variables
8
Actuarial models: Frequency + Severity
Hypothesis Severities of individual losses are independent. Severities and frequencies are independent.
Model separately Frequency {Nt}
E.g. Poisson, negative binomial, Cox process,… Severity {Lnt}
E. g. Lognormal, Pareto, g-and-h, … Approximate the distribution of aggregate losses by
combining these distributions. [Panjer, FFT, MC,single-loss approximations and beyond]
Sums of heavy-tailed random variables
9
Aggregation of frequency and severity dists.
Sums of heavy-tailed random variables
10
Expected & unexpected loss
Sums of heavy-tailed random variables
11
Risk analysis
Calculate aggregate yearly loss distribution from the frequency and severity distributions.
Compute risk measures• Expected loss• Capital at Risk (CaR)
e.g. 99.9% percentile of the aggregate loss distribution.• Conditional CaR (Expected shortfall)
Expected loss, given that the loss is above CaR
11Sums of heavy-tailed random variables
12
Computational issues in risk analysis
Algorithms to compute risk measures Deterministic algorithms
• Discretized approximation (FFT, Panjer).• Piecewise approximation for the distribution
• Body approximation: Not very important• Tail approximation. E.g. Single loss
Monte Carlo algorithmsEmpirical compound distribution obtained by simulation.Computationally costly.
• Use variance reduction techniques• Hardware solutions: Grid computing, GPU’s, …
12Sums of heavy-tailed random variables
13
Distribution of aggregate losses
Consider L1 ,…, LN iid samples from the density f(x)For the individual losses
The aggregate loss Z = L1 +…+ LN is a random variable with density g(x)
13Sums of heavy-tailed random variables
on]distributi [Tail )(1Prob)(
')'(Prob)(0
xFxLxF
dxxfxLxFx
on]distributi [Tail )(1Prob)(
)(')'(Prob)( *
0
zGzZzG
zFdzzgzZzG Nz
14
The lognormal distribution
exp() scale tails
)1,0(~ );exp(~0
log2
1exp2
1),;( 222
NZZXx
xx
xLNpdf
Sums of heavy-tailed random variables
15
Sums of iid lognormal rvs
Z = L1 +L2 +...+LN {Li} ~ iid LN(,) For sufficiently large , and specially at the tails,
the sum is dominated by the maximum.
Sums of heavy-tailed random variables
Ncentral limit
0.5 21.0 81.5 912.0 29812.5 2683383.0 656599703.5 4.37E+104.0 7.90E+134.5 3.88E+175.0 5.18E+21
16
Heavy-tailed distributions
A heavy-tailed distribution is such that the tails of the distribution decay more slowly than the exponential.
16Sums of heavy-tailed random variables
0,)(lim xx e
xF
17
Subexponential distributions
Subexponential distributions are a special class of heavy-tailed distributions.
Definition: F(x) is a subexponential distribution if
where L1 ,…, LN an iid sample from F(x)
17Sums of heavy-tailed random variables
2 ,lim
)(1)(1lim
)()(lim
1
1
**
NNxLP
xLLPxFxF
xFxF
N
x
N
x
N
x
18
One loss causes ruin [subexponential distributions]
One loss causes ruin: If the distribution of individual losses is subexponential, large values of the aggregate loss from N events are dominated by the maximum loss in one of the N events.
18Sums of heavy-tailed random variables
xLPNxFNxFN
xFxFxFxLP
xLLPxLLP
x
N
i
iNN
ii
NN
1
1
01
11
)()(1
)()(1)(11
),,max(1),,max(
2,1
,,maxlim
1
1
NxLLP
xLLP
N
N
x
19
Single loss approximation [Böcker + Klüppelberg, 2005]
From the definition of subexponential distribution
From the definition of VaRα
19Sums of heavy-tailed random variables
xxFNxGxLPNxLLP
xxLPNxLLP
N
N
as )(1)(111
as
11
11
)()(1)( 11 VaRGVaRLPVaRLP
NFVaRVaRFN
11)(11 1
20
Single loss approximation + mean correction
Single loss approximation for the aggregate loss distribution[Böcker + Sprittulla, 2006]
Single loss approximation + mean correction for VaRα
20Sums of heavy-tailed random variables
1][][
11)( 11LNE
NEFGVaR
LENENEpFpG
xNExFNExG
LL
L
;1][ ][
11~)(
as 1][1][~1
11
21
Beyond the single-loss approximation
For the Pareto distribution Based on [Omey & Willekens, 1986, 1987],
[Sahay et al., 2007] A. Sahay, Z. Wan, B. Keller, Operational risk capital: asymptoticsin the case of heavy-tailed severity. The journal of operational risk 2(2):61–72 (2007)[Suárez, 2010] http://www.fields.utoronto.ca/audio/09-10/forum-oprisk/suarez/[Degen, 2010] Degen, M. The calculation of minimum regulatory capital using single-loss approximations. Journal of Operational Risk 5(4):3–17 (2010)
For rapidly varying subexponential distributionsBased on [Barbe, McCormick & Zhang, 2007][Suárez & Fernández, 2011] http://www.risklab.es/es/jornadas/2011/index.html
Based on asymptotics with a shifted argument[Albrecher et al., 2010] H. Albrecher, C. Hipp, D. Kortschak, Higher-order expansions
for compound distributions and ruin probabilities with subexponential claims. Scandinavian Actuarial Journal 2010(2):105–135 (2010)
21Sums of heavy-tailed random variables
22
Correction by the mean from the perturbative expansion [Hernández et al. 2011]
Using the first order approximation and assuming that Q0 >> X
22Sums of heavy-tailed random variables
L
N
NQQLLENQQ
NEFFQ
)1( |)1(
][11
)0(0
)0()1(
1/11)0(
, |)1( ),( with
,2,1,0,!
121
!1
01/11
0
10
)(
22
101
00
QLLENQFQ
KQk
QQVaR
QQQQk
QQQVaR
N
K
k
kk
K
k
kk
23
Estimation of the mean can be difficult
23Sums of heavy-tailed random variables
24
Error in the sample estimate of the mean
24Sums of heavy-tailed random variables
25
Pre-CLT behavior in convergence to mean
25Sums of heavy-tailed random variables
= 0; = 3 mean = 90.0171
N = 104 N = 105N = 103N = 102N = 10
26
Perturbative expansion I
Consider a random variable Z = X + Y
Let Q0 = FX-1() be the percentile of X at probability level
Let Q be the percentile of Z at probability level :
From these two equations, one obtains
26Sums of heavy-tailed random variables
),(, yxfdxdy YX
yQ
),()( ,
00 yxfdxdyxfdx YX
X
),(0 ,
0
yxfdxdy YX
yQ
Q
27
Perturbative expansion II
Expressing Q as a perturbative series in
27Sums of heavy-tailed random variables
),(0 ,
0
yxfdxdy YX
yQ
Q
1
0 !1,
n
nnQ
nQQQQ
0
0
0
),(10
),(0
,1
,
QxYXxyQ
YX
yQQ
Q
yxfedy
yxfdxdy
x
28
Perturbative expansion III
Taylor expansion
In terms of Bell polynomials
28Sums of heavy-tailed random variables
n
nnx
QxYX
nx
n
n xyxfyQ
ndy
;),(
!10
0
,1
1
0
),(,,,!
10 ,211
1QxYXxnxxnx
n
n
yxfQQyQBn
dy
29
Perturbative expansion IV
Therefore, for the nth term in the series
29Sums of heavy-tailed random variables
,3,2;~;~
),(,,,~11
)(1
),,(11
),,(1
|),(0
;)(1
),(,,,0
11
,21
1
10
1
1
11
010,1
111
,211
0
0
mQQyQQ
yxfQQyQBQdymn
QfQ
xxBxmn
xxxBn
QXYEQyQfyQdy
xxBn
yxfQQyQBdy
mm
QxYXxmnxxmnm
n
mXn
mnmnm
n
mnnn
YX
QxYXxnxxnx
30
Perturbative expansion V
30Sums of heavy-tailed random variables
31
Perturbative expansion VI
31Sums of heavy-tailed random variables
32
Expansion around the maximum I
Consider a random variable ZN = L1+ L2 + ... + LN = L[1]+ L[2] + ... + L[N]
where the order statistics are L[1] L[2] ... L[N]
Identify XN with L[N] = max {L1, L2, ... , LN}YN with L[1]+ L[2] + ... + L[N-1]
Define ZN() = XN + YN [ZN = ZN(=1)]
32Sums of heavy-tailed random variables
33
Expansion around the maximum II
33Sums of heavy-tailed random variables
)()()( )(ProbProb)( 1][][][ xfxFNxfxFxLxLxF N
NNN
NN
34
Expansion around the maximum III
34Sums of heavy-tailed random variables
35
Stochastic frequency
35Sums of heavy-tailed random variables
36
Benchmarks for comparison (finite mean)
36Sums of heavy-tailed random variables
37
Benchmarks for comparison (divergent mean)
37Sums of heavy-tailed random variables
38
Lévy distribution
38Sums of heavy-tailed random variables
39
Lévy distribution (deterministic N)
39Sums of heavy-tailed random variables
40
Lognormal distribution (stochastic N)
40Sums of heavy-tailed random variables
41
Pareto distribution: dependence on
41Sums of heavy-tailed random variables
42
Pareto distribution (finite mean)
42Sums of heavy-tailed random variables
43
Pareto distribution (divergent mean)
43Sums of heavy-tailed random variables
44
Pareto distribution (asymptotic analysis)
44Sums of heavy-tailed random variables
Dominant for a > 1 up to order k a
/N: Expansion parameter
=1- 0+
45
Pareto distribution (asymptotic analysis)
45Sums of heavy-tailed random variables
Dominant for a < 1and for a >1 (k > a)
No recognizableexpansion parameter
=1- 0+
46
Asymptotic formulas to operational VaR
The single-loss approximation is insufficient, specially with Lower percentiles. Less heavy tails. Higher frequencies.
The single loss formula corrected by the mean Can be derived from the second order asymptotics. Accurate in a wide range of cases Estimation of the mean can be inaccurate if tails are
heavy The following orders in the perturbative approximation
Improve estimate and are easy to compute.46Sums of heavy-tailed random variables
47
Open questions
The perturbative approximation is a formal expansion with good numerical accuracy (at least for low orders)
Convergence properties?General asymptotic behavior?Resummation of the series?Transition to CLT regime?
Extend the analysis to Non-identically distributed random variables Dependent variables
47Sums of heavy-tailed random variables
48
Thanks for your attention!
48Sums of heavy-tailed random variables