Operational risk Dependencies and the Determination of Risk Capital Stefan Mittnik Chair of Financial Econometrics, LMU Munich & CEQURA finmetrics@stat.uni-muenchen.de Sandra Paterlini EBS Universität für Wirtschaft und Recht & CEQURA sandra.paterlini@ebs.edu Tina Yener Linde AG & CEQURA tina.yener@linde.com March 22, 2013 S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 1 / 25
Operational Risk: Heterogeneity of Events The heterogeneity of Operational Risk leads to the regulatory requirement of a separate modeling within 56 event type/business-line combinations. Business Lines Corporate Finance... Retail Brokerage Event Types Internal Fraud L 1,1... L 1,8........ Execution, Delivery & L 7,1... L 7,8 Process Management S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 2 / 25
Operational Risk: Total Risk Capital The quantity of interest is ( 56 ) VaR.999 (L) = VaR.999 L i ; (1) clearly, it is influenced by dependencies among cells i and j. However, Basel II prescribes to calculate Total Risk Capital as i=1 56 TRC = VaR.999 (L i ) ; (2) i=1 Only under certain qualifying conditions, banks may explicitly model dependencies. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 3 / 25
VaR and Subadditivity It can be shown (Frachot et al., 2004) that for the case of comonotonic risks, VaR co α (L i + L j ) = VaR α (L i ) + VaR α (L j ). (3) Comonotonicity translates into perfect positive correlation in the elliptical (and thus, also the popular Gaussian) world. For elliptical distributions, the sum of the single VaRs provides an upper bound and thus a worst case scenario for VaR α (L), VaR α (L i + L j ) VaR α (L i ) + VaR α (L j ). (4) S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 4 / 25
VaR and Subadditivity However, for non elliptical distributions, it may happen that VaR α (L 1 + L 2 ) > VaR α (L 1 ) + VaR α (L 2 ), (5) the reason being the lack of subadditivity of the VaR measure (Artzner et al., 1999). Does this mean that banks may not be rewarded for a more realistic dependency modeling by a decrease in risk capital, but instead be punished by an increase? Yes! (see, e.g., Embrechts et al., 2002) But, is this practically relevant? S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 5 / 25
Aim of our Analyses Based on n = 60 observations of monthly aggregate losses from the Italian DIPO 1 database, we aim at evaluating VaR (L i + L j ) (VaR (L i ) + VaR (L j )) }{{} =TRC (6) for different cells i and j of the event type/business line matrix. This task is non trivial, because it means analyzing the 99.9% quantile of a distribution estimated from a small sample with extreme data under consideration of dependencies. To model dependencies, we focus on Correlation, Copulas and Nonparametric Tail Dependence measures. 1 http://www.dipo-operationalrisk.it S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 6 / 25
Linear (Pearson) Correlation The well known fact that linear correlation is prone to extremes is quickly revealed by the data. For example, for event type combination (2. External Fraud; 5. Damage to Physical Assets), as the two most extreme observations drop out of the sample, correlation becomes negative. 01/03 12/07 (entire sample): 01/03 04/06 (2/3 of sample): l5 l5 l 2 ρ 2,5 = 0.0258 ρ 2,5 = 0.1144 l 2 S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 7 / 25
Linear (Pearson) Correlation Similarly, for event type combination (3. Employment Practices & Workplace Safety; 4. Clients, Products & Business Practice): 01/03 12/07 (entire sample): 01/03 04/06 (2/3 of sample): l4 l4 l 3 ρ 3,4 = 0.5284 ρ 3,4 = 0.5882 l 3 If we further remove the most extreme observation between 01/03 and 04/06, the correlation decreases to ρ 3,4 = 0.3232. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 8 / 25
Correlation: Results The well known sensitivity of linear correlation with respect to extremes leads to substantial variations, depending on the sample size. Its inability to capture possible nonlinear dependency structures provides another important reason for discarding linear correlation as a reliable measure of dependency. Rank correlations were also considered but not found to lead to considerably more stable results. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 9 / 25
Copulas Instead of boiling down dependency into one single number, copulas contain the dependence structure of a joint distribution. The central theorem of copula theory can be traced back to Sklar (1959) and summed up by C(u1, u2) 1 0.5 1 0 0.5 u 2 0 0 u 0.5 1 1 F i,j (l i, l j ) = C(F i (l i ), F }{{} j (l j )), (7) }{{} u i u j where C denotes the copula of L i and L j and U i, U j Unif(0, 1). S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 10 / 25
Copulas and Tail Dependence Tail dependence accounts for possibly nonlinear dependence among extremes and thus overcomes one drawback of correlation. ℓi Gaussian copula: no tail dependence S. Mittnik, S. Paterlini and T. Yener ℓj ℓj ℓj Different copulas imply different tail dependence structures. ℓi Gumbel copula: upper tail dependence ℓi Clayton copula: lower tail dependence CFS, 22.03.2013 11 / 25
Copulas Estimation We fit alternative parametric copulas (i.e.: Gaussian, Student t, Gumbel and Clayton). ET 2 ET 3 ET 4 ET 5 ET 6 ET 7 ET 1 0.136 0.275 0.243 0.147 0.182 0.251 ET 2 0.176 0.235 0.020 0.000 0.115 ET 3 0.574 0.318 0.272 0.297 ET 4 0.359 0.353 0.154 ET 5 0.037 0.000 ET 6 0.218 Upper tail dependence coefficient implied by the Gumbel copula parameter estimates. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 12 / 25
Copulas: Results Copula fitting suggests that some event type combinations are characterized by tail dependence, while others are not; i.e. (3;4) exhibits tail dependence, while (2;5) do not. However, we do neither find an overall best fitting copula, nor can we exclude any copula family considered. Again, the availability of a small data set affects the stability of estimation results. We also consider Nonparametric Tail Dependence measures and empirical results support the presence of quantile dependence for (3;4) and its absence for (2;5). S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 13 / 25
Range of Risk Capital Estimates Now, we want to assess the effects of realistic dependency structures on risk capital estimates. To this end, we estimate 250 99.9% VaR figures per model and event type combination, using different numbers of replications. For each event type combination, we use the copula parameter values obtained from Maximum Likelihood estimation. For the margins, a lognormal distribution was fitted and is used here to derive risk capital estimates. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 14 / 25
Range of Risk Capital Estimates We consider VaR.999 (L i + L j ) (VaR.999 (L i ) + VaR.999 (L j )) (VaR.999 (L i ) + VaR.999 (L j )) (8) under two different assumptions: 1 the Gaussian copula for all event type combinations, 2 the worst case copula, i.e., that copula yielding the highest tail/quantile dependence coefficient for the respective event type combination. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 15 / 25
Range of Risk Capital Estimates: Boxplots 30 30 30 15 15 15 Rel. Diff. (%) 0 Rel. Diff. (%) 0 Rel. Diff. (%) 0 15 15 15 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination B rc = 10,000 B rc = 50,000 B rc = 100,000 Range of simulated risk capital changes with Gaussian copula. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 16 / 25
Range of Risk Capital Estimates: Boxplots 30 30 30 15 15 15 Rel. Diff. (%) 0 Rel. Diff. (%) 0 Rel. Diff. (%) 0 15 15 15 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination B rc = 10,000 B rc = 50,000 B rc = 100,000 Range of simulated risk capital changes with worst case copula. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 17 / 25
Bounds on Risk Capital Estimates Obviously, increasing the number of replications for VaR calculations narrows the range of possible risk capital estimates. It is, thus, not clear which part of a change is due to the subadditivity problem, and which one is due to computational issues. A natural question is then: What could be the worst capital estimate? Statistically, this means to study whether there are theoretical bounds on VaR. This topic has been treated, for example, by Makarov (1981) and Frank et al. (1987), and recently by Embrechts and Puccetti (2006). S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 18 / 25
Bounds on Risk Capital Estimates The Fréchet Höffding bounds (Fréchet, 1951; Höffding, 1940) apply to any n dimensional copula, i.e., max(u 1 +... + u n n + 1, 0) }{{} C l (u) C(u) min(u). (9) }{{} C u(u) 1 1 1 C(u1, u2) 0.5 C(u1, u2) 0.5 C(u1, u2) 0.5 1 0 1 0 1 0 0.5 1 0.5 1 0.5 1 u 2 0 0.5 0 u 0.5 1 u 2 0 0 u 0.5 1 u 2 0 0 u 1 lower bound C l (u 1, u 2) Gaussian copula upper bound C u(u 1, u 2) S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 19 / 25
Bounds on Risk Capital Estimates: Assumptions The tightness of the bounds on VaR depends on the dependence assumption. We evaluate upper and lower bounds for three scenarios. 1 C 0 = C 1 = C l : We do not use any restriction on the dependence structure and thus use the lower Fréchet bound, C l. 2 C 0 = C 1 = u i u j : We assume that C u i u j, that is, we have positive quadrant dependence (PQD). 3 C 0 = C θ i,j CS, Ĉ1 = C θ i,j C : We take the Clayton Survival copula as lower bound, using the parameter values estimated for the DIPO data. For the survival copula of C 1, we accordingly assume the Clayton copula with respective parameter values. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 20 / 25
Bounds on Risk Capital Estimates: Boxplots 45 45 45 C0 = Cl C0 = Cl C0 = Cl C0 = C1 C0 = C1 C0 = C1 30 C0 = C γs CS, C1 = C γ C 30 C0 = C γs CS, C1 = C γ C 30 C0 = C γs CS, C1 = C γ C Rel. Diff. (%) 15 0 Rel. Diff. (%) 15 0 Rel. Diff. (%) 15 0 15 15 15 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination B rc = 10,000 B rc = 50,000 B rc = 100,000 Relative variations in simulated risk capital and theoretical bounds based on a Gaussian copula. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 21 / 25
Bounds on Risk Capital Estimates: Boxplots 45 45 45 C0 = Cl C0 = Cl C0 = Cl C0 = C1 C0 = C1 C0 = C1 30 C0 = C γs CS, C1 = C γ C 30 C0 = C γs CS, C1 = C γ C 30 C0 = C γs CS, C1 = C γ C Rel. Diff. (%) 15 0 Rel. Diff. (%) 15 0 Rel. Diff. (%) 15 0 15 15 15 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination 30 2/5 3/4 Event Type Combination B rc = 10,000 B rc = 50,000 B rc = 100,000 Relative variations in simulated risk capital and theoretical bounds based on a worst case copula. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 22 / 25
Bounds on Risk Capital Estimates: Results Risk capital estimates may increase when departing from the comonotonicity assumption. However, this effect depends on the presence of extremal (tail/quantile) dependence; such an increase may as well be caused by an insufficient number of replications in the simulation of losses. Theoretical bounds may help to assess which part of the change in risk capital is due to computational effects. The more restrictive the dependence assumptions used in deriving these bounds the more helpful they are. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 23 / 25
Conclusion The question whether risk capital estimates may increase/decrease compared to the comonotonicity case crucially depends on the presence of tail dependence and the ellipticity of the multivariate distribution. Simple methods, as correlations, may not lead to a complete and/or reliable picture of dependencies in operational risk losses. More sophisticated methods, such as copulas, could provide relevant information about dependencies and their effect on risk capital estimates. Risk capital estimates may increase when departing from the comonotonicity assumption. Theoretical bounds on VaR may help to assess which part of the change in risk capital stems from effects due to the computational setup. Serious efforts towards improving database for operational risk losses should be undertaken. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 24 / 25
Further Research What is the effect of modelling multivariate dependence (beyond the bivariate case) on Total Risk Capital? Robust and stable estimation of risk capital by maximum entropy methods S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 25 / 25
Contacts Prof. Stefan Mittnik, PhD Chair of Financial Econometrics, LMU Munich Akademiestr. 1/I, 80799 Munich, Germany. Tel. +49 (0) 8921803224, Fax. +49 (0) 8921805044 Email: finmetrics@stat.uni-muenchen.de Prof. Sandra Paterlini, PhD Chair of Financial Econometrics and Asset Management EBS Universität für Wirtschaft und Recht, Gustav-Stresemann-Ring 3 65189 Wiesbaden, Germany. Tel.: +49 (0) 611 7102 1217 ; fax: +49 (0) 611 7102 1217 Email: sandra.paterlini@ebs.edu Dr. Tina Yener Linde AG, Klosterhofstrasse 1, 80331 Munich, Germany Tel.: +49 (0) 89 35757 1614, Fax: +49 (0) 89 35757 1605, Email: tina.yener@linde.com S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 26 / 25
Acknowledgements & Copyright We are thankful to Andrea Resti, Claudia Pasquini, Claudia Capobianco, Marco Belluomini, and Vincenzo Bugge for helpful comments, and the Database Italiano delle Perdite Operative (DIPO) and its Statistical Committee for their support. The views expressed in this paper are those of the authors and do not necessarily reflect the viewpoints of DIPO or the DIPO Statistical Committee. The material cannot be copied, modified, distributed or displayed without the authors explicit written permission S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 27 / 25
References I Artzner, P., Delbaen, F., Eber, J., and Heath, D. (1999). Coherent Measures of Risk. Mathematical Finance, 9:203 228. Embrechts, P., McNeil, A., and Straumann, D. (2002). Correlation and dependence in risk management: properties and pitfalls. In Dempster, M., editor, Risk Management: Value at Risk and Beyond. Cambridge University Press. Embrechts, P. and Puccetti, G. (2006). Bounds for functions of multivariate risks. Journal of Multivariate Analysis, 97(2):526 547. Frachot, A., Roncalli, T., and Salomon, E. (2004). The correlation problem in operational risk. Working paper, Crédit Lyonnais. Frank, M. J., Nelsen, R. B., and Schweizer, B. (1987). Best possible bounds for the distribution of a sum a problem of kolmogorov. Probability Theory and Related Fields, 74(2):199 211. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 28 / 25
References II Fréchet, M. (1951). Sur les tableaux de corrélation dont les marges sont donnés. Annales de l Université de Lyon, 3(14):53 77. Höffding, W. (1940). Masstabinvariante korrelationstheorie. Schriften des Mathematischen Instituts und des Instituts fur Angewandte Mathematik der Universität Berlin, 5:179 233. Makarov, G. (1981). Estimates for the distribution function of a sum of two random variables when the marginal distributions are fixed. Theory of Probability and its Applications, 26:803 806. Sklar, A. (1959). Fonctions de répartition a n dimensions et leurs marges. Publications de l Institut de Statistique de L Université de Paris, 8:229 231. S. Mittnik, S. Paterlini and T. Yener CFS, 22.03.2013 29 / 25