Correction des exercices du livre
La Gestion des Risques Financiers
Thierry Roncalli
4 Décembre 2011
Ce document présente les corrections des questions de cours et des exercices des pages 535 à 552 du livre : Thierry Roncalli, La Gestion des Risques Financiers, Economica, deuxième édition, 2009. Notons que ce livre est cité sous la référence [TR-GDR] par la suite. Ce document présente aussi de nouveaux exercices (à partir de la page 46) et la correction de ceux-ci (de la page 49 à la page 61). Pour chaque exercice, nous indiquons à quelle page du livre le lecteur peut trouver l’énoncé correspondant. Par exemple, l’énoncé du premier exercice “valeur en risque d’un portefeuille long/short” se trouve à la page 538 du livre “La Gestion des Risques Financiers”.
1 Correction des exercices 2
1.1 Valeur en risque d’un portefeuille long/short (TR-GDR, page 538) . . . . . . . . . . . . . 2
1.2 Valeur en risque d’un portefeuille de gestion (TR-GDR, page 538) . . . . . . . . . . . . . 3
1.3 Construction d’un scénario de stress à partir de la théorie des valeurs extrêmes (TR-GDR,
page 542) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4 Stress-testing et statistiques d’ordre (TR-GDR, page 545) . . . . . . . . . . . . . . . . . . 7
1.5 Les mesures de risque (TR-GDR, page 549) . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.6 Valeur en risque crédit d’un portefeuille non granulaire (TR-GDR, page 542) . . . . . . . 13
1.7 Contribution en risque dans le modèle Bâle II (TR-GDR, page 541) . . . . . . . . . . . . 15
1.8 Le risque de contrepartie sur opérations de marché (TR-GDR, page 551) . . . . . . . . . . 17
1.9 Le risque opérationnel (TR-GDR, page 550) . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.10 Valorisation d’un CDS (TR-GDR, page 539) . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.11 Calibration du paramètre LGD (TR-GDR, page 543) . . . . . . . . . . . . . . . . . . . . . 28
1.12 Modélisation du temps de défaut (TR-GDR, page 544) . . . . . . . . . . . . . . . . . . . . 30
1.13 Les spreads de crédit (TR-GDR, page 546) . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.14 Valeur en risque d’une position optionnelle (TR-GDR, page 551) . . . . . . . . . . . . . . 34
1.15 Corrélation et copules (TR-GDR, page 547) . . . . . . . . . . . . . . . . . . . . . . . . . . 37
1.16 La distribution exponentielle (TR-GDR, page 548) . . . . . . . . . . . . . . . . . . . . . . 39
1.17 Les fonctions copules (TR-GDR, page 540) . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2 De nouveaux exercices 46
2.1 Calculs statistiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.2 Contribution en risque . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.3 Dérivés de crédit et corrélation de défaut . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.4 Modélisation de la perte en cas de défaut . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.5 Produits exotiques et gestion des risques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.6 Calcul des bornes de corrélation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.7 Le modèle exponentiel généralisé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3 Correction des exercices supplémentaires 49
3.1 Calculs statistiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.2 Contribution en risque . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.3 Dérivés de crédit et corrélation de défaut . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.4 Modélisation de la perte en cas de défaut . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.5 Produits exotiques et gestion des risques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.6 Calcul des bornes de corrélation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.7 Le modèle exponentiel généralisé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4 Éléments de correction des questions de cours 61
4.1 La réglementation Bâle II (TR-GDR, page 535) . . . . . . . . . . . . . . . . . . . . . . . . 61
4.2 Le risque de marché (TR-GDR, page 535) . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
4.3 Le risque de crédit (TR-GDR, page 536) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
4.4 Le risque opérationnel (TR-GDR, page 536) . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.6 Les dérivés de crédit (TR-GDR, page 538) . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
4.7 La gestion du risque de crédit (TR-GDR, page 538) . . . . . . . . . . . . . . . . . . . . . 64
La référence principale de cet exercice est (TR-GDR, pages 61-63). On note PA (t) (resp. PB (t)) la valeur de l’action A (resp. B) à la date t. La valeur du portefeuille est :
P (t) = xA · PA (t) + xB · PB (t)
avec xA et xB le nombre d’actions A et B. On en déduit que le PnL entre les dates t et t + 1 est:
PnL(t;t + 1) = P (t + 1) ? P (t)
avec RA (t;t + 1) et RB (t;t + 1) les rendements des actions A et B entre les dates t et t + 1.
1. On a xA = +1, xB = ?1 et PA (t) = PB (t) = 100. On en déduit que :
PnL(t;t + 1) = 100 × (RA ? RB)
On a RA ? RB ? N (0,?A?B) avec :
?A?B = ?0,20+ (?0,20)2 + 2 × 0,5 × 0,20 × (?0,20)
= 20%
La volatilité annuelle de cette position long/short est donc de 20%. Pour calculer la VaR pour un horizon 1 jour, on utilise la méthode du scaling (TR-GDR, page 74). On obtient2 :
La probabilité de perdre 2,89 euros par jour est de 1%.
2. On a PnL(t;t + 1) = 100×(RA ? RB). On utilise les données historiques pour calculer les scénarios de rendements joints (RA,RB). On en déduit ensuite la distribution empirique du PnL. Enfin, on calcule le quantile empirique correspondant. Avec 250 scénarios, le quantile empirique 1% est à mi-chemin entre la deuxième et troisième pires pertes :
3. L’expression du PnL devient (TR-GDR, pages 91-95) :
PnL(
avec CA (t) la valeur de l’option d’achat. On a :
CA (t + 1) ? CA (t) ? ? × (PA (t + 1) ? PA (t))
avec ? le delta de l’option. On en déduit que :
PnL(t;t + 1) = 50 × RA ? 100 × RB
À titre d’illustration, on obtient pour la VaR analytique :
car :
La VaR 1 jour 99% est donc passée de 2,89 euros à 2,50 euros.
1.2 Valeur en risque d’un portefeuille de gestion (TR-GDR, page 538)
1. On a (TR-GDR, page 63) :
PnL(t;t + 1) = P (t + 1) ? P (t)
avec P (t) la valeur du portefeuille à la date t et R(t;t + 1) le rendement du portefeuille pour la période [t;t + 1]. À la date t, P (t) est une variable certaine – P (t) est égal à 20000 euros – tandis que R(t;t + 1) est aléatoire. Par construction, R(t?;t + 1) suit une distribution N (µp,?p)
avec µp = x?µ le rendement espéré du portefeuille et ?p = x??x la volatilité du portefeuille. Soit VaR? la valeur en risque du portefeuille. Celle-ci vérifie :
Pr{PnL(t;t + 1) ? ?VaR} = 1 ? ?
ou de façon équivalente :
Pr{L(t;t + 1) ? VaR} = 1 ? ?
avec L(t;t + 1) = ?PnL(t;t + 1) la fonction de perte. On en déduit que :
PnL(
Pr{
On obtient finalement que :
Pour déterminer la VaR gaussienne du portefeuille, il suffit de calculer l’expression précédente avec les valeurs données de P (t), x, µ et ?. Dans cet exercice, on a µp = ?34% et ?p = 40%. La VaR 99% pour une période de détention d’un an est donc :
VaR1Y = 20000 × (2,33 × 0,40 + 0,34) = 25440
On applique la méthode du scaling pour calculer la VaR mensuelle (TR-GDR, page 74) :
2. Le quantile empirique 1% correspond à la statistique d’ordre 2,6 : 260. On procède par interpolation et on a (TR-GDR, page 66) :
VaR1D = ?(?1124 + 0,6 × (?1012 ? (?1124))) = 1056.8
Par la méthode du scaling, on obtient (TR-GDR, page 74) :
VaR1M = ?20VaR1D = 4726
La VaR historique 99% pour une période de détention d’un mois est égale à 4726 euros.
3. On écrit la relation de façon vectorielle (TR-GDR, pages 119-121) :
rt = BFt + ?t
On en déduit que :
var(rt) = B?B? + D
avec B la matrice n×m de sensibilité aux facteurs et D une matrice diagonale d’éléments.
(a) On a :
?p = ?x?B?B?x + x?Dx
Il vient donc que :
VaR = P (t) × ??1 (?) × ?x?B?B?x + x?Dx
(b) On en déduit que :
VaR = P (t) × ??1 (?) × ?(??F)2x?11?x + x?Dx
Si ?i2 = ?j2 = ?2, alors le portefeuille de VaR minimale satisfait le programme d’optimisation suivant :
x? = argmin
s.c. 1?x = 1
Le Lagrangien est :
avec Q = (??F)211? + ?2In. Les conditions de premier ordre sont :
On obtient x = ?Q?11. On en déduit aussi que :
d’où :
La solution optimale est donc :
La matrice( Q?1 est particulière. C’est une matrice telle que (Q?1)i,i = (Q?1)j,j et (Q?1)i,j =
Q?1)j,i car les termes de variance sont les mêmes ainsi que les termes de covariance. Le portefeuille de valeur en risque minimale est donc le portefeuille équipondéré (ou portefeuille 1/n) :
alors G est une distribution de Gumbel, de Fréchet ou de Weibull (TR-GDR, page 139).
2. On a (TR-GDR, page 140) :
Pr
(a) On a (TR-GDR, page 141) :
Fn (anx + bn) = (1 ? e??(??1x+??1 lnn))n
On en déduit que :
G
(b) On a :
F
On en déduit que :
G
(c) On a :
n
F
On en déduit que :
G
3. La distribution GEV permet de caractériser les trois types de distribution par une distribution unique. On a (TR-GDR, page 146) :
On en déduit que :
4. On a :
On en déduit que :
5. Le temps de retour est la durée moyenne entre deux dépassements successifs de la VaR (TR-GDR, pages 79-81). On a :
1
t?(VaR?) = × p 1 ? ?
avec ? le seuil de confiance et p la période de détention. Si ? est égal à 99% (resp. 95%) et si p est 1 jour (resp. 1 an), t?(VaR?) est égal à 100 jours (resp. 20 ans).
7. On a (TR-GDR, page 146) :
G?1 (?) = µ ? ???1 [1 ? (?ln?)??]
Pour ? = 1, on obtient :
G?1 (?) = µ ? ? (1 ? (?ln?)?1)
Par définition, on a t? = (1 ? ?)?1n. Le temps de retour t? est donc associé au seuil de confiance ? = 1 ? n/t?. On en déduit que:
On remplace ensuite µ et ? par les estimateurs du maximum de vraisemblance µˆ et ?ˆ.
8. On a : | ||
r(P1) | = | ?(0,01 + (12 ? 1)0,03) = ?34% |
r(P2) | = | ?(0,10 + (12 ? 1)0,02) = ?32% |
Le scénario de stress est plus sévère pour le portefeuille P1. Si on considère cette mesure de risque, le portefeuille P1 est donc plus risqué que le portefeuille P2.
1. Le stress-testing permet de calculer la perte maximale d’un portefeuille pour un scénario extrême donné (TR-GDR, pages 121-129). Il permet donc de compléter la valeur en risque. En effet, celleci est généralement basée sur les 260 derniers scénarios historiques journaliers. Ces scénarios ne correspondent pas forcément à des scénarios extrêmes ou adverses. Dans la réglementation, le stresstesting est utilisé pour identifier les scénarios les plus risqués pour la banque. Ils ne donnent pas lieu à une exigence de fonds propres. Néanmoins, ils doivent être réalisés de façon périodique (par exemple tous les trimestres) et les résultats doivent être communiqués aux autorités réglementaires.
Pour le risque de marché, un exemple de scénario de crise est une baisse des actions de 50% en un mois. On utilise aussi le stress-testing pour mesurer les risques de marché liés à des paramètres inobservables (pour le pricing d’options exotiques par exemple). Pour le risque de crédit, un exemple de scénario de crise est un défaut de plusieurs pays de la zone Euro.
2. On a (TR-GDR, pages 131-133) :
Gn (x) = Pr{max(X1, ,Xn) ? x}
La log-vraisemblance associée à un échantillon (x1, ,xm) de la statistique d’ordre Xn:n est donc (TR-GDR, page 135) :
Pour chaque statistique d’ordre Xn:n, nous pouvons donc calculer la log-vraisemblance ?(n) et déterminer les estimations µˆ(n) et ?ˆ (n). Il suffit ensuite de tester les hypothèses µˆ(1) = µˆ(2) = = µˆ(n) et ?ˆ (1) = ?ˆ (2) = = ?ˆ (n) (TR-GDR, pages 136-137).
4. Le temps de retour est la durée moyenne entre deux dépassements successifs de la VaR. On a (TR-GDR, page 81) :
avec ? le seuil de confiance et p la période de détention. Si ? = 99%, t?(VaR?) est égal à 100 jours (resp. 4 ans et 100 ans) si la période de détention est 1 jour (resp. 10 jours et 1 an). On a :
Le temps de retour associé aux quantiles G?11 (99%), G?51 (99%) et G jours, 500 jours et 2200 jours.
5. On cherche le seuil de confiance ? tel que (TR-GDR, page 150) :
On a donc :
ou encore :
? = 1 ? 20 × 0,001 = 98%
est donc égal à 100
1.5 Les mesures de risque (TR-GDR, page 549)
1. (a) On a (TR-GDR, page 29) :
VaR(?) = inf {x : Pr{L ? x} ? ?} ES(?) = E[L|L ? VaR(?)]
(b) On suppose que F est continue. On a donc VaR(?) = F?1 (?). On en déduit que :
ES(
On considère le changement de variable t = F(x). Comme dt = f (x) dx et F(?) = 1, on obtient :
ES(
(c) On a :
On a donc :
On en déduit que :
et :
On remarque aussi que :
x permet de localiser la distribution tandis que ? contrôle la queue de distribution. On a : ?
On en déduit que :
F?1 (?) = x? (1 ? ?)???1
On obtient donc :
ES(
Comme ? > 1, on a alors et :
ES( VaR( )
(d) On a :
ES(
On considère le changement de variable t = ??1 (x ? µ). On obtient (TR-GDR, page 498) :
ES(
Comme ?? (x) = ?x?(x), on a :
On considère l’intégration par partie avec u(t) = ?t?1 et v? (t) = ?(t) et on obtient :
On considère de nouveau l’intégration par partie avec u(t) = t?3 et v? (t) = ?(t) et on obtient :
d’où on en déduit que :
On reprend l’expression de ES(?) et on pose x = ??1 (?) :
ES(
On en déduit que ES(?) ? VaR(?) lorsque ? ? 1 car :
(e) Pour la distribution gaussienne, l’expected shortfall et la valeur en risque coincident pour des niveaux élevés de probabilité, ce qui n’est pas le cas de la distribution de Pareto, qui présente une queue épaisse. L’utilisation de la distribution de Pareto peut donc produire des mesures de risque beaucoup plus élevées que celles basées sur la distribution gaussienne.
2. On a (TR-GDR, page 29) :
R(L1 + L2) = E[L1 + L2] = E[L1] + E[L2] = R(L1) + R(L2)
R(?L) = E[?L] = ?E[L] = ?R(L)
R(L + m) = E[L ? m] = E[L] ? m = R(L) ? m
On remarque que :
On en déduit que si F1 (x) ? F2 (x), alors F?11 (t) ? F?21 (t) et E[L1] ? E[L2]. R est donc une mesure de risque cohérente.
3. On procède de la même façon pour montrer que R(L) = E[L] + ? (L) est une mesure de risque cohérente.
4. On a :
?i | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
Pr{L = ?i} Pr L ?i | 0,2 0,2 | 0,1 0,3 | 0,1 0,4 | 0,1 0,5 | 0,1 0,6 | 0,1 0,7 | 0,1 0,8 | 0,1 0,9 | 0,1 1,0 |
(a) On a VaR(50%) = 3, VaR(75%) = 6 et VaR(90%) = 7. On en déduit que :
3 × 10% + + 8 × 10%
ES(50%) == 5,5
60%
6 × 10% + + 8 × 10%
ES(75%) == 7,0
30%
7 × 10% + 8 × 10%
20%
(b) On doit construire une distribution bivariée telle que (TR-GDR, page 30) :
F
Pour cela, on peut s’aider des différents résultats sur les inégalités de Makarov (TR-GDR, pages 344-350). Par exemple, il suffit de considérer la somme ordinale de la copule C+ pour u1 ? ? et u2 ? ? et d’une autre copule C? pour u1> ? et u2> ? pour produire dans la plupart des cas une distribution bivariée qui ne vérifie pas la propriété de sous-additivité pour la valeur en risque. Prenons par exemple ? = 70% et C? = C?, on obtient la distribution bivariée suivante :
On a donc :
?i | 0 | 2 | 4 | 6 | 8 | 10 | 14 |
Pr{L1 + L2 = ?i} Pr L1 + L2 ?i | 0,2 0,2 | 0,1 0,3 | 0,1 0,4 | 0,1 0,5 | 0,1 0,6 | 0,1 0,7 | 0,3 1,0 |
Comme F?11 (80%) = F?21 (80%) = 6 et F, on en déduit que :
F
5. (a) On a (TR-GDR, page 498) :
(b) On a :
L ? N (x?µ,?x??x)
On en déduit que :
VaR(?) = x?µ + ??1 (?)?x??x
et :
On peut aussi utiliser le résultat précédent pour retrouver cette expression (TR-GDR, page 499).
On rappelle le résultat suivant. Soient A et B deux événements indépendants. On a :
Pr{A et B} | = | Pr{A} × Pr{B} |
Pr{A ou B} | = | Pr{A} + Pr{B} ? Pr{A} × Pr{B} |
La référence principale de cet exercice est (TR-GDR, pages 30 et 197-200). La perte aléatoire L a pour expression :
L = 1000 × LGD×D
avec LGD la perte en cas de défaut et D la variable aléatoire représentant le défaut. D est une variable de Bernoulli de paramètre PD.
1. On a L = 1000 × D. L peut donc prendre deux valeurs 0 et 1000. On a :
=
Pr{1000 × D = 0}
=
Pr{D = 0}
=
1 ? PD
On obtient finalement :
x | 0 | 1000 |
Pr{L = x} Pr L x | 99,5% 99,5% | 0,5% 100% |
{ ? }
On en déduit que la VaR 1 an à 99% est nulle.
2. Dans ce cas, L peut prendre trois valeurs 0, 500 et 1000. On a :
Pr{L = 0} | = | Pr{1000 × LGD×D = 0} |
= | Pr{D = 0} | |
et : | = | 99,5% |
Pr{L = 500} | = | Pr{1000 × LGD×D = 500} |
= | Pr{LGD = 50% et D = 1} | |
= | Pr{LGD = 50%} × Pr{D = 1} | |
= | 25% × 0,5% |
On obtient finalement :
x | 0 | 500 | 1000 |
Pr{L = x} Pr L x | 99,5% 99,5% | 0,125% 99,625% | 0,375% 100% |
On en déduit que la VaR 1 an à 99% est nulle.
3. Dans le premier cas, on obtient :
x | 0 | 1000 |
Pr{L = x} Pr L x | 90% 90% | 10% 100% |
{ ? }
Dans le deuxième cas, les résultats deviennent :
x | 0 | 500 | 1000 |
Pr{L = x} Pr L x | 90% 90% | 2,5% 92,5% | 7,5% 100% |
La VaR 1 an à 99% est égale à 1000 euros dans les 2 cas.
4. L’expression de la perte aléatoire L devient :
L peut donc prendre trois valeurs 0 (si A et B ne font pas défaut), 1000 (si A ou B fait défaut) et
2000 (si A et B font défaut). On a : | ||
Pr{L = 0} | = | Pr{DA = 0 et DB = 0} |
= | Pr{DA = 0} × Pr{DB = 0} | |
et : | = | 0.9 × 0.9 |
Pr{L = 1000} = Pr{DA = 0 et DB = 1} + Pr{DA = 1 et DB = 0}
= 2 × 0,9 × 0,1
On obtient :
x | 0 | 1000 | 2000 |
Pr{L = x} Pr L x | 81% 81% | 18% 99% | 1% 100% |
La VaR 1 an à 99% est égale à 1000 euros.
5. Dans ce cas, L prend cinq valeurs 0, 500, 1000, 1500 et 2000. On a :
Pr{L = 0} | = | Pr{DA = 0 et DB = 0} = 0,9 × 0,9 |
Pr{L = 500} | = | Pr{DA = 1, LGDA = 50% et DB = 0} + Pr{DA = 0, DB = 1 et LGDB = 50%} |
= | 2 × 0,9 × 0,25 × 0,1 | |
Pr{L = 1000} | = | Pr{DA = 1, LGDA = 100% et DB = 0} + Pr{DA = 0, DB = 1 et LGDB = 100%} + Pr{DA = 1, LGDA = 50%, DB = 1 et LGDB = 50%} |
= | 2 × 0,9 × 0,75 × 0,1 + 0,1 × 0,25 × 0,1 × 0,25 | |
Pr{L = 1500} | = | Pr{DA = 1, LGDA = 100%, DB = 1 et LGDB = 50%} + Pr{DA = 1, LGDA = 50%, DB = 1 et LGDB = 100%} |
= | 2 × 0,1 × 0,75 × 0,1 × 0,25 | |
Pr{L = 2000} | = | Pr{DA = 1, LGDA = 50%, DB = 1 et LGDB = 50%} |
= | 0,1 × 0,75 × 0,1 × 0,75 |
On obtient :
x | 0 | 500 | 1000 | 1500 | 2000 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Pr{L = x} Pr L x | 81% 81% | 4,5% 85,5% | 13,5625% 99,0625% | 0,375% 99,4375% | 0,5625% 100% 1.7 Contribution en risque dans le modèle Bâle II (TR-GDR, page 541)1. On suppose que (TR-GDR, page 179) : (a) les pertes en cas de défaut LGDi sont indépendantes des temps de défaut ?i ; (b) le défaut de la i-ième créance dépend d’un ensemble de facteurs communs X1, ,Xm ; (c) le portefeuille est infiniment granulaire; il n’y a donc pas de concentration d’exposition sur une contrepartie, ce qui se peut se traduire par la propriété mathématique suivante : EAD 2. Soit R une mesure de risque. La contribution en risque de la créance i est le produit de l’exposition de la créance i et du risque marginal (TR-GDR, page 497) : RCi = EAD ? EADi Dans le cas d’une mesure de risque convexe, on a : La mesure de risque est alors égale à la somme des différentes contributions en risque. 3. On a : EL = E[L] UL = F?1 (?) ? EL avec F la fonction de répartition de L. Si les temps de défaut sont indépendants, on obtient : I E[L | X1, ,Xm] = ?EADi ×E[LGDi] × PDi = EL i=1 L n’est donc plus aléatoire, on a Pr{L = EL} = 1. Il vient que F?1 (?) = EL et UL = 0. 4. On note (TR-GDR, page 180) : I g (x) = ?EADi ×E[LGDi] × PDi (x) i=1 d’où : E[L | X = x] = g (x) On a donc : F(?) = Pr{L ? ?} = Pr{g (X) ? ?} Comme EADi ? 0 et E[LGDi] ? 0, g (x) est une fonction monotone croissante (resp. décroissante) si PDi (x) est une fonction croissante (resp. décroissante) de x. Dans le cas où g (x) est une fonction
On en déduit donc que : F[LGDi] × PDi (H?1 (?)) Dans le cas où g (x) est une fonction décroissante, on a :
F[LGDi] × PDi (H?1 (1 ? ?)) Si au moins une des expositions EADi est négative, la fonction g (x) n’est plus forcément monotone. On ne vérifie plus g (X) ? ? ? X ? g?1 (?) (cas croissant) ou g (X) ? ? ? X ? g?1 (?) (cas décroissant) et les expressions de F?1 (?) ne sont plus valides. On ne peut donc pas utiliser le modèle de capital réglementaire si le portefeuille de crédit contient une (ou plusieurs) exposition nette négative sur une contrepartie. Pour la gestion de portefeuille de crédit, cela implique que l’achat d’une protection de type CDS sur une contrepartie ne doit se faire que pour réduire l’exposition sur cette contrepartie et non pour prendre une position short sur celle-ci. 5. On a (TR-GDR, pages 181-182) : PDi = Pr{?i ? Mi} = Pr{Zi ? Bi} = ?(Bi) d’où Bi = ??1 (PDi). On a aussi : Pr ? En utilisant la paramétrisation de ?(x1,x2;?) donnée à la page 296 de TR-GDR, on a : 6. La probabilité de défaut conditionnelle étant une fonction décroissante de x, on a : F Comme fF?1 (?)) = (??F?1 (?))?1, on en déduit que (TR-GDR, page 505) : 1 f (?) = avec ? = F?1 (?). 7. ? est la corrélation constante des actifs : cor(Zi,Zj) = E[ZiZj] On peut estimer ce paramètre par maximum de vraisemblance en considérant un historique de taux de défaut (TR-GDR, pages 474-477). 1.8 Le risque de contrepartie sur opérations de marché (TR-GDR, page 551)1. (a) Notons MtMA (C) et MTMB (C) le mark-to-market des banques A et B pour le contrat C. On doit vérifier théoriquement que : MtMA+B (C) = MTMA (C) + MTMB (C) = 0 Dans le cas de contrats listés, la relation précédente est généralement vérifiée. Dans le cas de contrats OTC, il n’y a pas de prix de marché. Dans ce cas, la banque utilise un modèle pour valoriser le contrat OTC. On parle d’une valorisation en mark-to-model. La relation MTMA+B (C) = 0 (ou MTMA+B (C) ? 0) n’est plus forcément vérifiée pour plusieurs raisons : les deux contreparties n’utilisent pas forcément le même modèle de valorisation, les paramètres utilisés ne sont pas exactement les mêmes, les systèmes de provisionnement pour le risque de modèle (appelés aussi réfactions) ne sont pas comparables, etc. Dans cet exemple, on obtient : MTMA+B (C1) = 10 ? 11 = ?1 MTMA+B (C2) = ?5 + 6 = 1 MTMA+B (C3) = 6 ? 3 = 3 MTMA+B (C4) = 17 ? 12 = 5 MTMA+B (C6) = ?5 + 5 = 0 MTMA+B (C7) = 1 + 1 = 2 Il y a donc un contrat qui vérifie exactement la relation (c’est le contrat C6). La relation est presque vérifiée pour deux autres contrats (ce sont les contrats C1 et C2). Pour le contrat C7, chaque contrepartie pense être gagnante et affiche un mark-to-market positif! Enfin, il y a trois contrats qui présentent de grosses différences (ce sont les contrats C3, C4 et C5). (b) On a (TR-GDR, pages 216-217) : EAD =max(MTM(Ci),0) On obtient donc : EADA = 10 + 6 + 17 + 1 = 34 EADB = 6 + 9 + 5 + 1 = 21 (c) On a (TR-GDR, page 217) : EAD = maxMTM(Ci),0) On obtient donc : EADA = max(10 ? 5 + 6 + 17 ? 5 ? 5 + 1,0) = max(19,0) = 19 EADB = max(?11 + 6 ? 3 ? 12 + 9 + 5 + 1,0) = max(?5,0) = 0 (d) On a (TR-GDR, page 217) : EADA = max(10 ? 5 + 6,0) + 17 + 1 = 29 EADB = max(?11 + 6 ? 3,0) + 9 + 5 + 1 = 15 (a) On a (TR-GDR, pages 220-221) : – L’exposition au défaut d’un contrat OTC est égale à : e(t) = max(MTM(t),0) avec MTM(t) le mark-to-market du contrat à la date future t. e(t) est donc une variable aléatoire. – La distribution de l’exposition au défaut e(t) pour la date future t est notée F[0,t]. L’exposition future potentielle (potential future exposure) est le quantile de la distribution F[0,t] pour un seuil de confiance ? donné : PFE? (0;t) = F?[01,t] (?) – L’exposition maximale (peak exposure) est le maximum des expositions futures potentielles : PEPFE? (0;t) Cette quantité est aussi connue sous le nom d’exposition future potentielle maximale (maximum potential future exposure ou MPFE). EE(0;t) = E[e(t)] = ? xdF[0,t] (x) – L’exposition positive attendue (expected positive exposure) est la moyenne pondérée dans le temps des expositions attendues : EPE(0; EE(0,t) dt – L’exposition attendue effective (effective expected exposure) est l’exposition maximale attendue jusqu’à la date t : EEE(0;EE(0;?) = max(EEE(0;t?),EE(0;t)) – L’exposition positive attendue effective (effective expected positive exposure) est la moyenne pondérée dans le temps des expositions attendues : EEPE(0;EEE(0;t) dt – Si e(t) = ??t? avec ? ? U[0,1], on obtient : F avec x ? [0,??t]. On a donc : 1 PFE? (0;t) = F?[0,t] (?) = ???t PE? (0) = ???T EE(0; EPE(0; EEE(0; EEPE(0; – Si e(t) = exp(??t?) avec ? ? N (0,1), on obtient : F avec x ? [0,?]. On a donc : PFE PE? (0) = exp(??T??1 (?)) EE(0; EPE(0; EEE(0; EEPE(0; – Si, on obtient : F avec. On a donc : PFE PE? (0) = 1{t < t?} · PFE? (0;t) + 1{t ? t?} · PFE? (0;t?) EE(0; EPE(0; EEPE(0;EEE(0;t) dt avec : Cette exercice est plus difficile que le précédent, car la fonction e(t) n’est pas monotone croissante. Elle est croissante pour et ensuite décroissante. Ceci explique que l’on doit déterminer pour PE? (0) et EEE(0;t) l’optimum t?. À titre de comparaison, nous reportons les différentes fonctions EE(0;t), EPE(0;h), EEE(0;t) et EEPE(0;h) pour ces deux exercices sur les graphiques 1 et 2. Figure 1 – Exposition au défaut si e(t) = exp(??tN (0,1)) Figure 2 – Exposition au défaut si (b) On a (TR-GDR, pages 224-225) : e(t) = max(x1 + ?1W1 (t),0) On en déduit que : avec f (x) la fonction de densité de e1 (t). Comme e(t) ? N (x1,?1?t), on obtient : On considère le changement de variable y = ?1?1t?1/2 (x ? x1). On obtient : car ?(?x) = ?(x) et ?(?x) = 1 ? ?(x). (c) On a : e(t) = max(x1 + x2 + ?1W1 (t) + ?2W2 (t),0) On en déduit que : En utilisant le résultat de la question (b), il vient que : (d) On a représenté la fonction E[e(t)] sur le graphique 3 dans le cas x1 = x2 = 0 et ?1 = ?2. On note que la moyenne de l’exposition au défaut est une fonction croissante du temps et de la volatilité et que l’impact d’un contrat de compensation globale peut être très important, surtout dans le cas d’une corrélation faible ou négative. Figure 3 – Espérance de l’exposition au défaut E[e(t)] de la banque A 1.9 Le risque opérationnel (TR-GDR, page 550)1. (a) On rappelle que la fonction de densité de la distribution LN (µ,?) est (TR-GDR, pages 239) : On en déduit que la fonction de log-vraisemblance est : Pour calculer les expressions des estimateurs µˆ et ?ˆ, on peut procéder de deux façons. #1 On sait que si Y ? LN (µ,?), alors X = lnY ? N (µ,?) et que les estimateurs µˆ et ?ˆ du maximum de vraisemblance de X ? N (µ,?) sont : #2 On maximize la log-vraisemblance et on en déduit les estimateurs µˆ et ?ˆ : (µ,ˆ ?ˆ) = argmax?(µ,?) Cela revient à résoudre les conditions du 1er ordre : On a : et donc : On a : et donc : Pour m ? 1, on a : On considère les changements de variable y = ??1 (lnx ? µ) et z = y ? ?. On obtient : On en déduit que : On peut donc estimer les paramètres µ et ? par la méthode généralisée des moments en considérant les moments empiriques suivantes (TR-GDR, pages 239) : (b) On a : On en déduit que la fonction de log-vraisemblance est : La condition du 1er ordre est : On en déduit que : L’estimateur du maximum de vraisemblance est donc : Pour m ? 1, on a : On en déduit que : et : On remarque aussi que : On peut donc estimer les paramètres µ et ? par la méthode généralisée des moments en considérant les moments empiriques suivantes : La mise en place de la méthode GMM peut se faire en considérant soit le premier moment hi,1 (?), soit le second moment hi,2 (?) soit les deux moments joints (hi,1 (?),hi,2 (?)). (c) L’expression de la densité conditionnelle dans le cas (ii) est (TR-GDR, page 242) : La distribution conditionnelle est donc une distribution de Pareto de même paramètre ? avec x = H. L’estimateur ?ˆ du maximum de vraisemblance devient donc : ? Pour la distribution (i), on obtient (TR-GDR, page 243) : Les deux premiers moments de Li | Li ? H sont donc : hi,1 (µ,?) = Li ? m(µ,?) hi,2 (µ,?) = (Li ? m(µ,?))2 ? (v (µ,?) ? m2 (µ,?)) (d) On peut utiliser un test dérivé de la représentation QQ-plot ou un test basé sur les statistiques d’ordre (et plus particulièrement la statistique du maximum). 2. On considère un échantillon de T nombres de pertes {N1, ,NT}. On suppose que le nombre de pertes suit une distribution de Poisson P (?). (a) On a : n Pr ! On en déduit que la log-vraisemblance a pour expression : ! La condition du premier ordre est : L’estimateur du maximum de vraisemblance est donc : 1Y (b) On sait que : N1Y = NQ1 + NQ2 + NQ3 + NQ4 D’après la propriété de la distribution de Poisson, si N1Y ? P (?) alors N1Q ? P (?/4). On en déduit que : 1Q 3. (a) On sait que la durée d entre deux pertes consécutives supérieures à ? est une distribution exponentielle de paramètre ?(1 ? F (?)) (TR-GDR, page 251). On en déduit que : (b) On peut demander aux experts d’évaluer les temps de retour dj de différents scénarios ?j et on peut calibrer les paramètres ? et ? par la méthode des moments en remarquant que (TR-GDR, pages 252-253) : 1.10 Valorisation d’un CDS (TR-GDR, page 539)1. Si on note N le notionnel et s le spread, la jambe fixe (ou de prime) annuelle payée par l’acheteur B de la protection est (TR-GDR, page 410) : avec R le taux de recouvrement. On en déduit le diagramme des flux suivant (TR-GDR, page 411) : Figure'4 – Diagramme des flux du point de vue de l’acheteur de la protection$ S’il n’y a pas de défaut avant les 3 ans, l’acheteur B paye donc au vendeur A 10000 euros pendant 6 semestres. Le PnL du vendeur A est donc égal à 60000 euros. Si la contrepartie X fait défaut dans 2 ans et 2 mois, le PnL de l’acheteur B est égal à : PnLB = 600000 ? 4 × 10000 = 560000 2. Si le spread de la contrepartie X augmente sept mois plus tard et vaut 1000 pb, alors l’acheteur B de la protection gagne de l’argent s’il retourne sa position dans le marché puisque la probabilité de défaut de X a augmenté. La nouvelle valeur de la jambe de prime semestrielle est : Comme il reste encore 5 versements semestriels de la jambe de prime, le PnL de l’acheteur B est égal à (si on néglige l’effet d’actualisation) (TR-GDR, page 410) : PnLB = 5 × (50000 ? 10000) ? 10000 = 190000 On a la relation suivante (TR-GDR, page 411) : PD La probabilité de défaut initiale est donc : PD pb La nouvelle probabilité de défaut devient : PD ? 1.11 Calibration du paramètre LGD (TR-GDR, page 543)1. On calcule la perte en cas de défaut moyenne empirique de la classe de risque C : 100 × 0% + 100 × 25% + 600 × 50% + 100 × 75% + 100 × 100% mˆ LGD = 50% Soit L la perte aléatoire du portefeuille. On a (TR-GDR, page 457) : 100 L = ?EADi ×LGDi ×1{?i ? 1} i=1 avec EADi l’exposition au défaut de la créance i qui vaut 10000 euros, LGDi la perte en cas de défaut aléatoire de la contrepartie i et ?i le temps de défaut aléatoire de la contrepartie i. Comme les variables aléatoires LGDi et ?i sont indépendantes, on en déduit que : EL = LGDi ×1{?i ? 1}] 100 = ?E[EADi ×LGDi ×1{?i ? 1}] = ?EADi ×E[LGDi] × E[1{?i ? 1}] = ?EADi ×E[LGDi] × PDi i=1 avec PDi la probabilité de défaut annuelle de la contrepartie i. On obtient donc : = 5000 euros 2. L’espérance de perte par créance est donc égale à 50 euros. Si on suppose que l’exigence des fonds propres porte aussi sur la perte espérée, alors il convient d’impacter le bilan en considérant une charge supplémentaire de capital par créance égale à 50 euros. Si ce n’est pas le cas, alors il convient de provisionner la perte espérée au niveau du compte de résultat. Dans ce cas, cela revient à considérer que la marge commerciale est de 600 euros par créance. 3. On a : mˆ LGD = 50% et : ?ˆLGD = 25% On rappelle que : µLGD ?LGD2 On en déduit que (TR-GDR, pages 452-453) : (1 ? µLGD) b = a µLGD et : (a + b)2 (a + b + 1)?LGD2 (1 2 ? LGD = a µLGDµLGD On obtient finalement que : LGD LGD b = ?LGD2 (1 µLGD) En utilisant les expressions numériques de µLGD et ?LGD, on a : 4. On peut considérer le portefeuille précédent comme infiniment granulaire puisque c’est un portefeuille homogène avec 100 créances. On sait dans ce cas que l’unexpected loss dépend seulement de l’espérance de la perte en cas de défaut. Comme l’espérance de la distribution uniforme est 50% (comme la distribution Bêta calibrée), alors il n’y a pas de différence dans ce cas entre utiliser une distribution uniforme ou la distribution Bêta calibrée. Ce résultat est valable pour le modèle Bâle II à un seul facteur, mais aussi dans le cas où il y a plusieurs facteurs (TR-GDR, page 458). 1.12 Modélisation du temps de défaut (TR-GDR, page 544)1. Si t? [t?m?1,t?m[, on a (TR-GDR, page 427) : S Par itération, on obtient : S On en déduit que : f (t) = ??tS(t) ? = ??t (S(t?m?1)e?m?m(t?m?tm1 ?1)) = ?mS(tm? ?1)e?? (t?t? ) = ?mS(t) = ?me??mi=1?1 ?i(t?i?ti??1)e??m(t?t?m?1) 2. On a donc : S e? 15e? 2( ?5) si t >? 5 avec ?1 = 0,01 et ?1 = 0,02. On vérifie que : S car F(?) ? U[0,1]. Soit U ? U[0,1]. On en déduit que : ? = S?1 (u) Si u ? S(5), on a e??1t = u ou encore t = ???11 lnu. Si u >S(5), on a e??15e??2(t?5) = u ou encore t = 5 ? ??21 (lnu + ?15). On simule donc ? de la façon suivante : ???1 1 lnu si u ? e??115 ? =5 ? ??1 (lnu + ? 5) si u > e?? 5 2 1 Comme e??15 = e?0,01×5 = 0,951, les valeurs simulées de ? sont donc : 163,44 pour u = 0,04 ? =75,98 pour u = 0,23 ? 3,05 pour u = 0,97 3. On a : S On en déduit que (TR-GDR, page 431) : avec S. On considère la classe de risque A. On en déduit que : On obtient finalement les résultats suivantes :
? 4. Soit ? (t) la matrice de transition entre les dates 0 et t. On appelle générateur de ? (t) la matrice ? = (?i,j) définie par (TR-GDR, page 441) : ? (t) = et? où eM est l’exponentielle matricielle de M. On en déduit que (TR-GDR, page 443) : On a donc ? = (ˆ ln?)/2. Le générateur est Markovien si ?4j=1?i,j = 0 et ?i,j ? 0 si i ?= j. Ces conditions sont vérifiées, on en déduit que ?ˆ est une générateur Markovien. Si (t) = 1 ? e?i ? (t)e4 = 1 ? e?iet?ˆe4 On en déduit que (TR-GDR, page 467) : e?i ?ˆet?ˆe4 = 1 ? e?iet?ˆe4 On utilise donc la formule précédente pour calculer ?(t) directement à partir du générateur Markovien. À long terme, le système de transition est stationnaire. La probabilité de défaut d’une contrepartie appartenant à la classe de risque A est donc à long terme la même que celle d’une contrepartie appartenant à la classe de risque B ou C. On a donc : ?A (?) = ?B (?) = ?C (?) = 147,6 bp A court terme, ce n’est pas le cas et les probabilités de défaut sont ordonnées selon la classe de risque : ?A (0) < ?B (0) < ?C (0) Il est donc tout à fait normal que l’on obtienne une courbe croissante pour le rating A et une courbe décroissante pour le rating C (TR-GDR, page 448). 1.13 Les spreads de crédit (TR-GDR, page 546) 1. On a (TR-GDR, page 427) : F Sf (t) = ?e??t On note X = S(?). On a X ? [0,1] et : Pr{X ? x} = Pr{S(?) ? x} = Pr ? ? S?1 (x)} = S({S?1 (x)) 2. On a (TR-GDR, pages 409-411) : avec s le spread du CDS et N le notionnel. Le paiement trimestriel de la prime de protection explique le facteur 1/4 dans la formule de JF. On en déduit le diagramme des flux suivant : 'Figure 5 – Diagramme des flux du point de vue de l’acheteur de la protection $ 3. La marge à la monnaie ou le spread du CDS est la valeur de s telle que le prix du CDS soit nul (TR-GDR, page 410) : P (t) = E[JF ? JV ] On a la relation triangulaire suivante (TR-GDR, page 410) : s ? ? × (1 ? R) 4. Soit PD la probabilité de défaut annuelle. On a : PD = 1 car ? est généralement petit (? ? 10%). On en déduit donc que : PD On a : 2% PD = = 267 pb 1 ? 25% Une stratégie de valeur relative consiste à construire son propre modèle théorique pour valoriser le spread de CDS et à acheter ou vendre la protection en fonction de la position du spread estimé par rapport au spread de marché (TR-GDR, page 413). 1.14 Valeur en risque d’une position optionnelle (TR-GDR, page 551)1. L’expression du PnL est (TR-GDR, page 91) : PnL = C (St+1) ? C (St) avec C (St) la valeur de l’option à la date t. On a : ? = ?SC = 50% = 4% et St+1 = (1 + Rt+1)St avec Rt+1 le rendement journalier: Rt+1 ? N (0,2%) (a) On a (TR-GDR, page 93) : PnL On en déduit que : PnL ? N (0,??St)
(b) On a (TR-GDR, page 94) : PnL On en déduit que : E[PnL] = et : E[PnL On a Rt+1 = ?X avec X ? N (0,1). Il vient que : car E[X] = 0, E[X2] = 1, E[X3] = 0 et E[X4] = 3. L’expression de l’écart-type du PnL est donc : ? (PnL) = L’approximation gaussienne du PnL est alors la suivante : PnL On en déduit que la VaR gaussienne est : En utilisant les valeurs numériques des différents paramètres, on obtient : VaR(99%) = 2,29 euros (c) Soit L = ?PnL la perte. La VaR Cornish-Fisher est égale à (TR-GDR, page 82) : VaR(?) = z? (?1,?2) × ? (L) avec et z? = ??1 (?). ?1 et ?2 sont respectivement le coefficient de skewness et l’excès de kurtosis associés à la distribution de la perte L. On a : PnL = ? avec X ? N (0,1). On rappelle que :
On en déduit que : et : L’expression des moments centrés est donc : E[(PnL et : L’expression de ?1 et ?2 est donc (TR-GDR, page 82) : ?1 (L) = ??1 (PnL) E[(PnL?E[PnL])3] E[(PnL?E[PnL])4] En utilisant les valeurs numériques des différents paramètres, on obtient ? (L) = 1,0016, ?1 (L) = ?0,2394, ?2 (L) = 0,0764, z? (?1,?2) = 2,1466 et VaR(?) = 2,15 euros. (d) La VaR est réduite avec la méthode Cornish-Fisher à cause de l’effet de skewness . Comme celle-ci est négative, cela revient à réduire la queue de distribution de probabilité de la perte et à augmenter la masse de probabilité des pertes plus petites que la perte moyenne. Il est donc normal que la valeur en risque soit plus petite que celle obtenue avec l’approximation gaussienne (car on rappelle que l’excès de kurtosis est très faible). 2. On a : PnL = ? En utilisant un historique de rendement, on peut calculer la valeur du PnL avec la formule précédente en remplaçant Rt+1 par les scénarios historiques. Il suffit ensuite de classer les valeurs simulées du PnL et de prendre la statistique d’ordre correspondant au quantile 1 ? ? (TR-GDR, pages 94-95). Pour la méthode de Monte Carlo, on remplace les scénarios historiques par des scénarios simulés de Rt+1 ? N (0,2%). 1.15 Corrélation et copules (TR-GDR, page 547)1. D’après le thèorème de Sklar, on a (TR-GDR, page 261) : ?2 (x1,x2) = C?X1,X2? (?(x1),?(x2)) Soit F la distribution bivariée du vecteur aléatoire (Y1,Y2) avec Y1 = ?(X1) et Y2 = ?(X2). On a
F1 (y1) = Pr{Y1 ? y1} = y1 On en déduit que Y1 ? U[0,1] et Y2 ? U[0,1] et que la décomposition canonique de F est : F(y1,y2) = C?X1,X2? (y1,y2) On a donc montré que : C?Y1,Y2? = C?X1,X2? Simuler la copule gaussienne de paramètre ? revient à simuler le vecteur aléatoire (Y1,Y2). Pour cela, on simule (x1,x2) une réalisation du vecteur gaussien (X1,X2) et on applique la transformation (?(x1),?(x2)). 2. On a : On a (TR-GDR, page 484) : ??S1,S2? = 1 ? ?(a,b) ? R × R+ : S1 = a + bS2 On en déduit que : Comme (X1,X2) est un vecteur gaussien, on en déduit que a = 0 et : Comme E[X1] = 0, b = 0. La condition implique que ?1 = ?2. On obtient donc le résultat suivant : ??S1,S2? = 1 ? ?1 = ?2 Corréler parfaitement les mouvements browniens W1 (t) et W2 (t) dans le modèle Black-Scholes n’est pas suffisant pour corréler parfaitement les prix des actifs S1 (t) et S2 (t). Il faut aussi imposer que ?1 = ?2. Si ?1 ?= ?2 et W1 (t) = W2 (t), on obtient alors que (TR-GDR, page 484) : ??S1,S2? < 1 3. On a : C?X1,X2? = ?C? + (1 ? ?)C+ (a) On a (TR-GDR, page 269) : F(x1,x2) = ? max(?(x1) + ?(x2) ? 1) + (1 ? ?)min(?(x1),?(x2)) On en déduit que : E[X1X2] = ?? x1x2 dF(x1,x2) = ? × (?1) + (1 ? ?) × (+1) On obtient donc : ??X1,X2? = E[X1X2] La corrélation linéaire de X1 et X2 est nulle lorsque : (b) Soient. On a (TR-GDR, page 284) : (c) On a : On calcule donc la corrélation empirique ?ˆ et on en déduit que l’estimateur des moments est (TR-GDR, page 328) : ?ˆ= 1 ? ?ˆ 2 1.16 La distribution exponentielle (TR-GDR, page 548)1. D’après la formule de Bayes, on a : Pr ? > t ? > s Pr ? > s Pr{? > t | ? > s} = { { ? } } Pr{? > t} = Pr{? > s} S(s) On dit que la distribution exponentielle possède la propriété d’absence de mémoire (TR-GDR, page 427). Pr{? > t | ? > s} = Pr{? > t ? s} Cela veut dire que la distribution de ? reste inchangée quelque soit la date d’origine pour mesurer ? et n’est pas influencée par le passé. Cette propriété est très intéressante dans la modélisation du risque de crédit car on utilise généralement la date présente pour mesurer ?. Si la propriété d’absence de mémoire n’est pas vérifiée, cela voudrait dire que modéliser le temps de défaut d’une entreprise en octobre 2011 dépendrait par exemple de sa date de création. La propriété Markovienne de la distribution exponentielle permet ainsi d’oublier tout ce qui s’est passé avant octobre 2011. 2. (a) On a : Pr et : Pr On en déduit que (TR-GDR, page 427) : (b) ?1 et ?i sont co-monotones si et seulement si (TR-GDR, page 274) : U1 = Ui S1 (?1) = Si (?i) Soit ?+ = max(?1, ,?n) et ?? = min(?1, ,?n). On en déduit que : Pr et : Pr On obtient finalement que :
Pr Pr 3. (a) Si la fonction de dépendance de (?1,?2) est C+, alors ?1 et ?2 sont co-monotones. En utilisant les résultats de la question 2b, on obtient : (b) Si la fonction de dépendance de (?1,?2) est C?, on a (TR-GDR, page 274) : U1 = 1 ? U2 S1 (?1) = 1 ? S2 (?2) ? ? ? ? Il existe bien une fonction f telle que ?1 = f (?2) avec : (c) On sait que ? ? [?min,?max] avec ?min et ?max les coefficients de corrélation linéaire de (?1,?2) pour C = C? et C = C+ (TR-GDR, page 283). On sait aussi que ? = 1 (resp. ? = 1) si et seulement si il existe une fonction f linéaire croissante (resp. décroissante) telle que ?1 = f (?2) (TR-GDR, page 284). D’après les résultats de la question 3a, on a ?1 = f (?2) pour C = C+ avec : La fonction f (t) est linéaire et croissante, on en déduit que ?max = 1. D’après les résultats de la question 3b, on a ?1 = f (?2) pour C = C? avec : La fonction f (t) est bien décroissante, mais elle n’est pas linéaire. On en déduit donc que ?min = 1? . On obtient finalement que : ?1 < ? ? 1 (d) Pour C = C?, on sait qu’il existe une fonction décroissante f telle que (TR-GDR, page 274) : X2 = f (X1) On sait aussi que ? atteint sa borne basse si C = C? (TR-GDR, page 284) et que ? = ?1 s’il existe une fonction linéaire décroissante f telle que : X2 = a + bX1 On doit donc vérifier que X2 = a + bX1 avec b < 0. Dans le cas où X1 = +?, on obtient : Comme X2 est borné par 0, on en déduit donc que : a = +? Si le support des variables aléatoires X1 et X2 est [0,+?], il n’existe donc pas de fonction linéaire décroissante f telle que X2 = f (X1). On en déduit que la borne ? = ?1 ne peut pas être atteinte. 1.17 Les fonctions copules (TR-GDR, page 540)1. (a) On a : C(u,0) = C(0,u) = 0 C(u,1) = C(1,u) = u On en déduit que C est une fonction copule (TR-GDR, page 260). (b) On a (TR-GDR, page 289) :
On a : on obtient le résultat suivant (TR-GDR, page 281) : (c) On simule tout d’abord (U1,U2) en utilisant la méthode des distributions conditionnelles (TRGDR, pages 315 et 318) et on applique la méthode de l’inversion (TR-GDR, page 314) : (d) On a (TR-GDR, page 261) : On en déduit que (TR-GDR, page 330) : 2. (a) On a : S(0,0) = 1 S(?,?) = 0 et : (b) On a : S1 (x1) = S(x1,0) = exp(?x1) Comme U1 = S1 (X1), on en déduit l’expression suivante de la fonction copule : C avec u˜ = ?lnu. On obtient la copule EV appelée Gumbel II (TR-GDR, page 312). 3. (a) ? est une fonction de classe C2 qui vérifie ?(1) = 0, ?? (u) < 0 et ??? (u) > 0 (TR-GDR, page 293). (c) C’est la copule de Gumbel (TR-GDR, page 294) : C(u1,u2) = exp(?[(?lnu1)? + (?lnu1)?]1/?) On a (TR-GDR, page 308) : C(u1t,ut2) = exp(?[(?lnu1t)? +?(?lnut1)?]1/?) = exp(?(t? [(?lnu1) + (?lnu1)?])1/?) = exp(?t[(?lnu1)? + (?lnu1)?]1/?) = Ct (u1,u2) C’est donc une copule de valeurs extrêmes. (d) Soit f une fonction. On note y = f (x). On a dy = ?xf (x) dx, x = f?1 (y) et dx = ?yf?1 (y) dy. On en déduit que : On a donc : C Soient v1 et v2 deux nombres aléatoires uniformes. L’algorithme des distributions conditionnelles revient à poser : On en déduit que (TR-GDR, page 320) : On obtient l’algorithme de Genest et MacKay (1986). 4. (a) Soit X3 une variable gaussienne indépendante de X1 et X2. On a : X2 = ?X1 + ?1 ? ?2X3 On en déduit que (TR-GDR, page 296) : On obtient donc que : Pr (b) On a (TR-GDR, page 296) : C et : C On utilise la méthode des distributions conditionnelles. Soient v1 et v2 deux nombres aléatoires uniformes. On a : ou encore : 5. Une copule de valeurs extrêmes vérifie la propriété suivante pour tout t > 0 (TR-GDR, page 308) : C On a : C C et : C En effet, si on prend u1 = u2 = 3/2 et t = 2, on obtient : C et : C 2 De nouveaux exercices2.1 Calculs statistiquesSoit L la variable aléatoire représentant la perte d’un portefeuille. On note f et F les fonctions de densité et de répartition associées. 1. Exprimez F(x) à partir de f. 2. Donnez la définition de F?1 (?). 4. Même question si L ? LN (µ,?). 5. Même question si L ? E (?). 6. Calculez F?1 (95%) et F?1 (99%) dans le cas de la distribution discrète suivante :
{ }2.2 Contribution en risqueNous notons L la perte d’un portefeuille de n créances, et xi l’exposition au défaut de la i-ième créance. Nous avons : avec ei la perte unitaire de la i-ième créance. Nous notons F la fonction de distribution de L. 1. On suppose que e = (e1, ,en) ? N (0,?). Calculez la valeur en risque au seuil de confiance ?. 2. En déduire la valeur en risque marginale de la i-ième créance. Définissez alors la contribution en risque de la i-ième créance. 3. Vérifiez que la valeur en risque marginale est égale à : Interprétez ce résultat. 4. On se place dans le modèle de risque de crédit Bâle II. On a : ei = LGDi ×Di avec Di = 1{?i < Mi} l’indicatrice de défaut et Mi la maturité de la i-ième créance. Quelles sont les conditions à vérifier pour obtenir le résultat suivant : E[ei | L = F?1 (?)] = E[LGDi] × E[Di | L = F?1 (?)] 5. On suppose que le défaut intervient avant la maturité Mi si une variable latente Zi passe en dessous d’une certaine barrière Bi : ?i ? Mi ? Zi ? Bi On modélise Zi = ??X + ?1 ? ??i avec Zi, X et ?i trois variables aléatoires gaussiennes centrées réduites et indépendantes. X est le facteur (ou le risque systémique) et ?i est le risque individuel. Calculez la probabilité de défaut conditionnelle. 6. Montrez que, dans le modèle Bâle II, nous avons : E[ei | L = F?1 (?)] = E[LGDi] × E[Di | X = ??1 (1 ? ?)] 8. On suppose que le portefeuille est homogène, c’est-à-dire que les créances ont la même exposition au défaut, la même distribution de perte en cas de défaut et la même probabilité de défaut. En utilisant le résultat suivant : avec ?2 (x,y;?) la fonction de répartition de la distribution gaussienne bivariée de corrélation ? sur l’espace [??,x]×[??,y], calculez l’expected shortfall dans le cadre du modèle Bâle II. Commentez ce résultat. 2.3 Dérivés de crédit et corrélation de défaut1. On considère un CDS 3M sur une contrepartie X de maturité 3 ans et de notionel 1 million d’euros. Le spread actuel du CDS est égal à 200 pbs. Donnez le diagramme des flux du CDS en supposant que la jambe de protection est payée au moment du défaut et que le taux de recouvrement est fixe et égal à 60%. 2. On suppose que le temps de défaut de X est exponentiel de paramètre ?. On considère que le taux d’intérêt instantanné est constant et égal à r. Donnez l’expression théorique du spread s du CDS. 3. En utilisant un modèle de Merton, vous estimez que le paramètre ? est égal à 250 pbs. Quelle est la position d’arbitrage que vous pouvez mettre en place? 4. Définissez la copule Normale de dimension n. 5. On considère un panier de n crédits. Qu’appelle-t-on un CDS first-to-default (F2D), un CDS secondto-default (S2D) et un CDS last-to-default (L2D)? 6. Définissez la notion de corrélation de défaut. Quel est l’impact de la corrélation de défaut sur le spread des trois CDS précédents? 7. On suppose maintenant que n = 3. Montrez la relation suivante : L2D avec sCDSi le spread CDS du i-ième crédit. 2.4 Modélisation de la perte en cas de défaut1. Quelles différences faites vous entre le taux de recouvrement et la perte en cas de défaut? 2. On considère une banque qui octroie en moyenne 250000 crédits par an. Le montant moyen d’un crédit est égal à 50000 euros. On estime que la probabilité de défaut moyenne est égale à 1% et que le taux de recouvrement moyen est égal à 65%. Le coût total annuel du service contentieux est égal à 12,5 millions d’euros. Quelle est la perte en cas de défaut moyenne? 3. On rappelle que l’expression de la densité d’une distribution Bêta de paramètres a et b est : avec B. (a) Pourquoi la distribution Bêta est un bon candidat pour modéliser le paramètre LGD? Quel jeu de paramêtres (a,b) correspond à la distribution uniforme? (b) On considère un échantillon (x1, ,xn) de n pertes en cas de défaut. Ecrivez la fonction de log-vraisemblance. Déduisez-en les conditions de premier ordre de l’estimateur du maximum de vraisemblance. (c) On rappelle que les deux premiers moments de la distribution Bêta sont : Déduisez-en l’expression de l’estimateur des moments. 2.5 Produits exotiques et gestion des risques1. Définissez les notions de mark-to-market et mark-to-model. Pourquoi le backtesting de la VaR d’un portefeuille de dérivés exotiques pose-t-il un problème? 2. Qu’appelle-t-on le risque de modèle? 3. Quelles sont les différentes façons pour calculer la VaR d’un portefeuille de dérivés exotiques? 4. On considère la vente d’une option exotique sur un sous-jacent dont le prix actuel est 100. À cette date t, la valeur de l’option est égale à 6,78 euros. (b) A la date t + 2, la valeur du sous-jacent devient 100 alors que la volatilité implicite passe de 20% à 22%. Le trader constate un PnL négatif de ?2,37 euros. Pourquoi le PnL du trader est négatif? Pouvez-vous expliquer le PnL par les sensibilités sachant que le delta ?t+1 est égal à 43%, que le gamma ?t+1 vaut 2% et que le véga ?t+1 est estimé à 38%? (c) Quelles conclusions faites-vous en terme de risque de modèle? 2.6 Calcul des bornes de corrélation1. Donnez les définitions mathématiques des copules C?, C? et C+. 2. Définissez la copule Normale C(?) bivariée de corrélation ?. 3. Quelles sont les interprétations probabilistes des trois copules définies à la question 1? Déduisez-en que C(?=?1) = C?, C(?=0) = C? et C(?=1) = C+. 4. On considère le vecteur aléatoire (?,LGD) qui modélise la loi jointe du défaut ? et de la perte en cas de défaut LGD d’une contrepartie. On suppose que ? ? E? et LGD ? U[0,1]. (a) Montrez que la dépendance de (?,LGD) est maximale lorsque : LGD+e??? ? 1 = 0 (b) Montrez que la corrélation ?(?,LGD) vérifie l’inégalité suivante : |?(?,LGD) (c) Commentez ces résultats. 2.7 Le modèle exponentiel généralisé1. On note F et S les fonctions de répartition et de survie associée à la variable aléatoire ?. Définissez la fonction S(t) et déduisez-en l’expression de la fonction de densité f (t) associée. 2. Donnez la définition du taux de risque ?(t). Déduisez-en que le modèle exponentiel correspond au cas particulier ?(t) = ?. 3. Comment peut-on simuler la variable aléatoire ? dans le cas ?(t) = ?. 4. On suppose maintenant que : Donnez l’expression de la distribution de survie S(t). Déduisez-en l’expression de la fonction de densité f (t). Vérifiez que : constante S(t) avec S et f les fonctions de survie et de densité du temps de défaut ? associé au CDS, R le taux de recouvrement et T la maturité du CDS. Donnez l’égalité triangulaire lorsque ? ? E?. 6. On considère maintenant que ? est le modèle exponentiel généralisé de la question 4. On suppose que l’on dispose de trois spreads CDS de marché s1, s2 et s3 de maturité respectives 3 ans, 5 ans et 7 ans. Montrez que la calibration du modèle exponentiel généralisé revient à résoudre un système de 3 équations d’inconnues ?1, ?2 et ?3. Quel est le nom donné à cette méthode de calibration? 3 Correction des exercices supplémentaires3.1 Calculs statistiques1. Si on suppose que le support inférieur de L est ??, on a : F 2. On a : F?1 (?) = inf {x : F(x) ? ?} Dans le cas d’une fonction de distribution continue, on obtient : Pr{L ? F?1 (?)} = ? 3. On note ?(x) la fonction de répartition standardisée N (0,1). On a : Pr{L ? x} = ? Pr Pr On en déduit que : x ? µ ?1 (?) = ? ? On obtient finalement que : F?1 (?) = x = µ + ? · ??1 (?) Application : Si µ est égal à 10 millions d’euros et ? est égal à 3 millions d’euros, on a :
car ??1 (0,95) = 1,65 et ??1 (0,99) = 2,33. 4. Si L ? LN (µ,?), alors lnL ? N (µ,?). On en déduit que : = Pr{lnL ? lnx} d’où : On obtient : F?1 (?) = exp(µ + ? · ??1 (?)) 5. On a F(x) = 1 ? e??x. Comme 1 ? e??x = ?, on obtient : F?1 (?) = ?ln(1 ? ?) ? 6. On a :
{ ? }On en déduit que F?1 (95%) = 8 et F?1 (99%) = 10. 3.2 Contribution en risqueNous notons L la perte d’un portefeuille de n créances, et xi l’exposition au défaut de la i-ième créance. Nous avons : avec ei la perte unitaire de la i-ième créance. Nous notons F la fonction de distribution de L. 1. On a : On en déduit que : 2. On a : La valeur en risque marginale de la i-ième créance est donc (TR-GDR, page 497) : MR La contribution en risque de la i-ième créance est le produit de l’exposition et du risque marginal : RCi = xi × MRi 3. On considère le vecteur aléatoire Y = (e,L). Par construction, Y est un vecteur gaussien avec (TR-GDR, page 497) : La distribution conditionnelle de e lorsque L = ? est gaussienne et on a : On obtient donc : La VaR marginale de la créance i est donc égale à l’espérance conditionnelle de la perte unitaire ei de la créance i sachant que la perte du portefeuille est égale à la valeur en risque. 4. On fait les hypothèses suivantes (TR-GDR, page 179) : (i) la perte en cas de défaut est indépendante du temps de défaut; (ii) le portefeuille est infiniment granulaire. 5. On a : Pr car : PDi = Pr{?i ? Mi} = Pr{Zi ? Bi} = ?(Bi) 6. Dans le modèle Bâle II, on a : L = F?1 (?)
E[ei | L = F?1 (?)] = E[LGDi] × E[Di | X = ??1 (1 ? ?)] 7. On a (TR-GDR, page 182) : RCi = xi × MRi Figure 6 – Comparaison VaR / ES dans le modèle Bâle II 8. D’après la question 6, on a : L ? F?1 (?) ? X ? ??1 (1 ? ?) On note EAD l’exposition au défaut. Comme on a : L = n × EAD×E[LGD] on en déduit que : ES( = n × EAD = n × EAD = n × EAD×E[LGD] × (1 ? ?)?1 × C(1 ? ?,PD;??) avec C la copule gaussienne. L’expected shortfall est donc de la forme : ES(?) = n × EAD×E[LGD] × (1 ? ?)?1 × P avec P la probabilité correspondant au quadrant inférieur? [0,1 ? ??]×[0,PD] de la copule gaussienne de paramètre ?. On note aussi que l’on retrouve le paramètre ? des modèles Raroc (TR-GDR, page 524). On vérifie que l’on a pour ? = 0 : [LGD] × PD En fait, le rapport ? entre l’expected shortfall et la VaR est égale à : On peut montrer que ce rapport ? est une fonction croissante de ? et une fonction décroissante de PD. À titre d’illustration, on reporte sur le graphique 6 les valeurs de ES(?) et VaR(?) lorsque n = 1000, EAD = 1, E[LGD] = 50% et ? = 99%. 3.3 Dérivés de crédit et corrélation de défaut 1. On a : JF = = 5000 euros Au moment du défaut, la jambe de protection est égale à : JV = (1 ? 60%) × 1000000 = 400000 euros On en déduit le diagramme des flux suivant (TR-GDR, page 411) : ' Figure 7 – Diagramme des flux du point de vue de l’acheteur de la protection $ 2. On a (TR-GDR, page 410) : 3. A partir de l’égalité triangulaire, on en déduit que l’intensité de défaut valorisée par le marché est : pb On vend donc le CDS car on anticipe que le spread va baisser. 4. On a : C(u1, ,un;?) = ?n (??1 (u1), ,??1 (un);?) avec ?n la fonction de répartition de la distribution gaussienne à n dimensions. 5. Un CDS first-to-default (resp. second-to-default et last-to-default) est un CDS dont l’événement de crédit porte sur le premier (resp. second et dernier) défaut observé (TR-GDR, pages 413-416). sF2D (resp. sL2D) est une fonction décroissante (resp. croissante) de ? (TR-GDR, page 415). L’influence de ? sur le spread sS2D est plus complexe puisque la relation n’est pas monotone. 7. Acheter un CDS sur chaque crédit revient à acheter un CDS F2D, plus un CDS S2D et enfin un CDS LTD. On a donc: L2D 8. Considérons un panier de 3 crédits. On peut alors assimiler le CDS F2D à la tranche junior d’un CDO, le CDS S2D à la tranche mezzanine et enfin le CDS LTD à la tranche senior. Par construction, l’expected loss du panier ne dépend pas de la corrélation de défaut ?, ce qui revient à dire que CDS n’est pas une fonction de ?. On en déduit que sF2D + sS2D + sL2D ne dépend pas de ?. ? a une influence seulement sur la répartition de l’expected loss entre sF2D, sS2D et sL2D : sF2D (?1) + sS2D (?1) + sL2D (?1) = sF2D (?2) + sS2D (?2) + sL2D (?2) Une sous-estimation de la corrélation de défaut entraîne donc une sur-estimation de sF2D et une sous-estimation de sL2D. Lors de la crise des subprimes, on a assisté à une augmentation importante de la corrélation de défaut. Le vendeur de la tranche junior doit donc subir une perte, mais pas le vendeur de la tranche senior. Si toutes les tranches ont souffertes, on ne peut pas l’expliquer par une sous-estimation de la corrélation de défaut. La seule explication est la sous-estimation générale des probabilités de défaut et des pertes en cas de défaut. 3.4 Modélisation de la perte en cas de défaut1. On a (TR-GDR, page 447) : LGD = 1 ? R + g avec g la perte dûe à la gestion du recouvrement. Elle correspond donc à toutes les dépenses engagées (frais d’avocat, etc.) pour recouvrir la créance. Supposons par exemple une créance de 100 euros. On recouvre 70 euros, mais ce recouvrement a coûté 10 euros. On a alors R = 70% mais LGD = 60%. 2. Le montant de crédit est égal à : = 12,500 milliards d’euros La perte après recouvrement est égale à : L = N × (1 ? R) × PD = 43,75 millions d’euros Comme le coût total annuel du service contentieux C est égal à 12,5 millions d’euros, on en déduit que : LGD = = 45% Ce chiffre est supérieur au 35% obtenu si on ne tenait pas compte du coût du service contentieux. (a) Elle permet d’obenir l’ensemble des configuations possibles (U et U inversé). La distribution uniforme correspond à a = 1 et b = 1 : (b) On a : n n ?(a,b) = ?nlnB(a,b) + (a ? 1)?lnxi + (b ? 1)?ln(1 ? xi) i=1 i=1 Les conditions de premier ordre de l’estimateur du maximum de vraisemblance sont donc : et : (c) On a (TR-GDR, pages 452-453) : LGD LGD LGD 2. Le risque de modèle est l’ensemble des risques auxquels s’expose une institution financière dont la gestion des risques repose sur des modèles mathématiques (TR-GDR, page 98). Celui-ci comprend le risque opérationnel (erreurs d’implémentation, etc.), le risque de paramètre, le risque de spécification et le risque de couverture. 3. On peut calculer la VaR d’un portefeuille de dérivés exotiques en utilisant la méthode historique, les méthodes basées sur les sensibilités (VaR delta, VaR delta-gamma, etc.), la VaR Cornish-Fisher ou la méthode de Monte Carlo (TR-GDR, pages 93-95). 4. On a : PnL = C (St,?t) ? C (St+1,?t+1) avec St la valeur du sous-jacent et ?t la volatilité implicite à la date t. (a) La valeur du sous-jacent a baissé alors que la volatilité implicite est restée la même. La probabilité d’exercer l’option diminue, ce qui implique que : C (St+1,?t+1) ? C (St,?t) On en déduit que : PnL ? 0 On a (TR-GDR, page 94) : PnL = En utilisant les valeurs numériques de ?t, ?t et ?t, on obtient : PnL On arrive donc à expliquer le PnL par les sensibilités. (b) La valeur du sous-jacent a augmenté ainsi que la volatilité implicite. La probabilité d’exercer l’option augmente donc puisque ?t > 0 et ?t > 0, ce qui implique que : C (St+2,?t+2) ? C (St+1,?t+1) On en déduit que : PnL ? 0 On a : PnL On obtient le résultat numérique suivant : PnL En utilisant les sensibilités, on n’explique que seulement 90% du PnL. 3.6 Calcul des bornes de corrélation1. On a (TR-GDR, page 269) :
2. C’est la copule associée à la distribution bivariée gaussienne ?2 (x1,x2;?) (TR-GDR, page 296) : C 3. Soient 2 variables aléatoires X1 et X2. On a (TR-GDR, page 274) : (i) C?X1,X2? = C? si et seulement si il existe une fonction f décroissante telle que X2 = f (X1); (ii) C?X1,X2? = C? si et seulement si il X1 et X2 sont indépendantes; (iii) C?X1,X2? = C+ si et seulement si il existe une fonction f croissante telle que X2 = f (X1). Dans le cas de la distribution gaussienne ?2 (x1,x2;?), on sait que (TR-GDR, page 284) : (i) si ? = ?1, alors X2 = ?X1 ; (ii) si ? = 0, X1 et X2 sont indépendantes; (iii) si ? = 1, alors X2 = X1 ; Si ? = ?1 (resp. ? = 1), on a trouvé une fonction décroissante f (x) = ?x (resp. croissante f (x) = x) telle que X2 = f (X1). On en déduit donc que C(?=?1) = C?, C(?=0) = C? et C(?=1) = C+. 4. On note : U1 = 1 ? e??? U2 = LGD (a) La dépendance est maximale lorsque C??,LGD? = C+. Dans ce cas, on a U1 = U2 (TR-GDR, page 274) et on en déduit que : LGD+e??? ? 1 = 0 (b) Si C??,LGD? = C?, on a U1 = 1 ? U2 et donc LGD = e???. On sait que ?(?,LGD) ? [?min (?,LGD),?max (?,LGD)] avec ?min (?,LGD) la corrélation correspondant à la copule C? et ?max (?,LGD) la corrélation correspondant à la copule C+. Comme ? ? E? et LGD ? U[0,1], il vient que : ? (LGD) = Dans le cas où C = C?, on a LGD = e???. On a donc : E[? LGD] = On en déduit que : Dans le cas où C = C+, on a LGD = 1 ? e???. On a donc : E[? LGD] = On en déduit que : ?max (?,LGD) = On vérifie donc l’inégalité suivante : |?(?,LGD) On remarque donc que |?(?,LGD)| ne peut excéder 0,866. Voici donc un exemple pour lequel on ne peut pas atteindre les bornes ?1 et +1. 3.7 Le modèle exponentiel généralisé1. On a (TR-GDR, page 426) : S(t) = Pr{? ? t} = 1 ? F(t) et : f (t) = ?tF(t) = ?tS(t) 2. ?(t) est le taux de défaut instantané (TR-GDR, page 426) : Pr{t ? ? ? t + ?|} = lim ??0+ ?Pr{? ? t} 1 Pr{t ? ? ? t + ?|} = lim = S(t) Dans le cas du modèle exponentiel E (?), on obtient : f (t) S(t) 3. On note X = S(?). On remarque que X ? [0,1] et que : Pr{X ? x} = Pr{S(?) ?1x} = Pr({? 1? S?) (x)} = S S? (x) = x On en déduit que S(?) ? U[0,1] et ? ? S?1 (U[0,1]). Soit u un nombre aléatoire uniforme. On a donc : 4. On a (TR-GDR, page 427) : si t ? 3 Ssi 3 < t ? 5 si t > 5 et : si t ? 3 si 3 < t ? 5 si t > 5 On vérifie que : f ( S 5. Si ? ? E?, on obtient (TR-GDR, page 411) : dt 6. On a donc : ou encore : avec : A partir de la première équation, on estime le paramètre ?ˆ1. On peut alors résoudre la deuxième équation qui ne comporte plus qu’une seule inconnue ?2 : On estime alors le paramètre ?ˆ2. La troisième équation devient : C’est une équation à une seule inconnue ?3, on peut donc estime facilement le paramètre ?ˆ3. Cette méthode de calibration successive des paramètres est appelée la méthode du bootstrap. 4 Éléments de correction des questions de cours4.1 La réglementation Bâle II (TR-GDR, page 535)1. Les différences entre le ratio Cooke et Bâle II sont décrites aux pages 34-37 de TR-GDR. Les points sur lesquels il faut insiter sont : (a) le Nouvel Accord intègre le risque opérationnel; (b) le traitement du risque de marché ne change quasiment pas; (d) le rôle central de la surveillance prudentielle et la place du Pilier II dans le système; (e) le renforcement des directions des risques au détriment des directions financières. 2. Le ratio McDonough correspond en fait à 3 ratios. Ceux-ci sont définis précisément à la page 34 de TR-GDR. La différence entre les ratios tier one et tier two est aussi expliquée à cette même page. 3. Les trois Piliers du Nouvel Accord sont (TR-GDR, pages 35-36) : (a) l’exigence minimale de fonds propres; (b) la surveillance prudentielle; (c) la communication financière et la discipline de marché. Le premier pilier permet de définir les différents éléments nécessaires au calcul des fonds propres et du ratio McDonough. Attention, les fonds propres calculés avec le premier pilier correspondent à une exigence minimale. Le régulateur peut demander des fonds propres supplémentaires. C’est le rôle du second pilier. Le second pilier concerne donc la ligne de conduite que doit adopter le régulateur pour surveiller l’établissement financier. Il comprend deux grands axes. Le premier est le processus de validation du premier pilier. Le second concerne le calcul de fonds propres supplémentaires. Ceux-ci peuvent être exigés si l’établissement financier présente un risque particulier ou si les hypothèses de calcul utilisées dans le modèle Bâle II ne sont pas applicables à l’établissement financier. Les problèmes de portefeuilles non granulaires et de concentration de risque sont donc traités dans ce second pilier (et non dans le premier pilier). Le troisième pilier cherche à promouvoir une discipline de marché plus efficace. Par exemple, l’établissement financier doit communiquer sur les différents modèles internes utilisés pour gérer le risque, doit publier les pertes les plus significatives, etc. 4.2 Le risque de marché (TR-GDR, page 535)1. La distinction entre trading book et banking book est expliquée à la page 47 de TR-GDR. Il faut aussi reproduire le tableau page 47 de TR-GDR et illustrer les différences de traitement des risques Actions et Taux dans le banking book. Pour cela, on peut considérer le risque d’une participation majoritaire et celui lié aux plans d’épargne logement. 2. La valeur en risque est définie aux pages 27-28 de TR-GDR. 3. Donnez la formule des fonds propres qui se situe à la page 53 de TR-GDR. Les points à développer sont : (a) la période de détention (au minimum 10 jours); (b) le seuil de confiance de 99%; (c) la fenêtre de lissage de 60 jours; (d) le coefficient ?. 4. On a besoin de calculer une valeur en risque VaR1 pour l’exigence de fonds propres et une valeur en risque VaR2 pour le backtesting. La différence entre ces deux mesures concerne la période de détention, qui est de 10 jours pour VaR1 et de 1 jour pour VaR2. Expliquez aussi la notion de scaling (TR-GDR, page 74) et pourquoi il est plus pertinent de faire le backtesting avec la mesure VaR2 (TR-GDR, pages 56-57). 5. La réponse à cette question se trouve aux pages 90-103 de TR-GDR. Les points à développer sont : (a) le problème de l’identification des facteurs de risque (TR-GDR, page 91); (b) la différence entre un prix mark-to-market et un prix mark-to-model (TR-GDR, page 96); (c) la difficulté du backtesting (TR-GDR, pages 96-98); (d) le risque de modèle (TR-GDR, page 98). 6. La procédure de backtesting est présentée aux pages 56-58 de TR-GDR. Insistez sur le facteur complémentaire ? et son incidence sur l’exigence de fonds propres. 8. Le risque de contrepartie sur opérations de marché est défini aux pages 214-216 de TR-GDR. Précisez pourquoi il diffère du risque de crédit et pourquoi la perte en cas de défaut (EAD) est aléatoire (TR-GDR, page 216). Les trois méthodes pour l’exposition au défaut sont données à la page 222 de TR-GDR. Faites ensuite le lien avec l’exigence de fonds propres au titre du risque de crédit. 4.3 Le risque de crédit (TR-GDR, page 536)1. La définition du défaut est donnée à la page 162 de TR-GDR. 2. Reproduisez le tableau de la page 165 de TR-GDR. Insistez sur la notion de notation externe (TR-GDR, page 168) et sur les 16 catégories de risque. 3. Les 4 paramètres de la méthode IRB (EAD, LGD, LGD et M) sont définis à la page 175 de TR-GDR. 4. Les différences entre FIRB et AIRB sont expliquées aux pages 209-210 de TR-GDR. 5. Ces quatre mesures de risque sont définies à la page 29 de TR-GDR (voir aussi la page 491). 6. En 2001, lors de la sortie du CP2, le Comité de Bâle avait choisi la valeur en risque comme mesure du risque de crédit. À la suite des négociations avec l’industrie bancaire, le Comité de Bâle a finalement opté pour l’unexpected loss. L’argument de l’industrie bancaire était que la perte espérée faisait déjà l’objet d’un provisionnement et qu’elle impactait donc le compte de résultat. Pour celle-ci, il n’y avait donc pas de raison qu’elle soit comptabilisée au niveau des fonds propres, c’est-à-dire au niveau 7. La gestion du risque dans la banque de détail repose sur des systèmes de scoring et une segmentation de la clientèle. L’important est de bien distinguer les clients sains des clients qui peuvent faire défaut. La problématique de la tarification n’est donc pas centrale. Le taux d’un crédit est fixe et le même pour tous les clients (ou du moins pour ceux qui appartiennent au même segment). La question posée est alors la suivante : doit-on ou non accorder le crédit au client? Pour la banque commerciale, la problématique de la tarification est centrale. La question posée devient la suivante : à partir de quel taux peut-on accorder le crédit au client? Dans la banque de détail, un client qui a une probabilité de défaut élevée ne pourra donc pas avoir un crédit. Dans la banque commerciale, il obtiendra un crédit à condition qu’il paye un spread suffisamment élevé pour rémunérer ce risque. Ceci explique que les méthodes Raroc de tarification concernent principalement la banque de financement et d’investissement alors que la banque de détail utilise les outils de scoring (TR-GDR, pages 516523). 4.4 Le risque opérationnel (TR-GDR, page 536)1. La définition du risque opérationnel est donnée à la page 228 de TR-GDR. Les pages 229-230 présentent différents exemples de risque opérationnel, ainsi que la classification Bâle II en 7 types de risque. 2. La différence entre BIA et STA est expliquée aux pages 231-233 de TR-GDR. Insistez sur le découpage en 8 lignes métiers (TR-GDR, page 230), sur les coefficients ?i (TR-GDR, page 232) et sur la procédure de validation de la méthode STA. 4. Le formalisme mathématique de la méthode LDA est expliquée aux pages 236-237 de TR-GDR. 4.5 Les fonctions copules (TR-GDR, page 537)1. La définition d’une copule est page 260 de TR-GDR. Les bornes de Fréchet sont données à la page 269 de TR-GDR. Insistez sur l’interprétation probabiliste de ces bornes et sur les notions de contre-monotonie et co-monotonie (TR-GDR, page 274). 2. L’ordre stochastique de concordance est présenté aux pages 270-271 de TR-GDR. Expliquez son importance dans la définition des mesures de dépendance (TR-GDR, pages 280-282). 3. Un mesure de concordance est une mesure signée alors qu’une mesure de dépendance est comprise entre 0 et 1. Dans le cas gaussien, la mesure de corrélation est assimilée à une mesure de concordance, alors que c’est la valeur absolue de la corrélation qui s’apparente à une mesure de dépendance. Rappelez aussi les différentes propriétés de ces deux mesures (TR-GDR, pages 280-282). Notez aussi qu’une mesure de concordance nulle n’implique pas que les variables aléatoires soient indépendantes contrairement à une mesure de dépendance nulle. 4. Le tau de Kendall, le rhô de Spearman et le coefficient de Gini sont définis à la page 281 de TR-GDR. 5. L’expression classique de la copule gaussienne est donnée à la page 262 de TR-GDR. Il est préférable d’utiliser la formulation du théorème 8 de la page 296, ce qui permet d’en déduire la copule conditionnelle. 6. Les copules archimédiennes sont définies à la page 293 de TR-GDR, tandis que la représentation des copules de valeurs extrêmes correspond à la définition 12 de la page 308. 8. Les principaux résultats de la théorie des valeurs extrêmes multiples sont les suivants (TR-GDR, page 305) : (a) les marges vérifient le théorème de Fisher-Tippet (TR-GDR, page 139); (b) la dépendance est une copule de valeurs extrêmes (TR-GDR, page 308); (c) la dépendance entre les extrêmes est forcément positive (TR-GDR, page 310); (d) la dépendance de queue est conservée (TR-GDR, page 314); (e) les extrêmes sont indépendants si les variables aléatoires sous-jacentes n’ont pas de dépendance de queue (TR-GDR, page 314). 9. La copule empirique de Deheuvels est présentée aux pages 322-325 de TR-GDR. 10. On peut répondre à cette question en montrant que la VaR n’est pas forcément une mesure de risque sous-additive (TR-GDR, page 30). On peut utiliser les inégalités de Makarov comme exemple (TR- GDR, page 344). Une illustration se trouve à la page 12 de ce document (question 4 (b) de l’exercice 1.5). 4.6 Les dérivés de crédit (TR-GDR, page 538)1. Donnez la définition d’un CDS (TR-GDR, page 409). Précisez en quoi un CLN est proche d’un CDS. Expliquez le montage d’un CLN en insistant sur l’aspect de structuration et la notion de SPV (TR-GDR, page 406). 2. Les ABS sont définis à la page 406 de TR-GDR tandis que les différentes catégories d’ABS sont répertoriées à la page 407 de TR-GDR. 3. Le mécanisme d’un CDO est expliqué aux pages 416-421 de TR-GDR. Il faut développer les points suivants : (a) la technique du tranching (TR-GDR, page 416); (b) la notion de points d’attachement et de détachement; (c) la définition des tranches equity, mezzanine et senior (TR-GDR, page 419); (d) les points de ressemblance entre un CDO et les CDS F2D, S2D et L2D. 4. La différence entre les tranches equity, mezzanine et senior est expliquée à la page 419 de TR-GDR. 4.7 La gestion du risque de crédit (TR-GDR, page 538)1. Les modèles exponentiel généralisé et à intensité sont présentés respectivement aux pages 427 et 433 de TR-GDR. 2. La définition d’un générateur markovien se trouve à la page 441 de TR-GDR. 3. On peut utiliser un modèle économétrique, par exemple le modèle LossCalc (TR-GDR, page 449), une distribution discrète (TR-GDR, page 456) ou une distribution Bêta (TR-GDR, page 452). Donnez les avantages à utiliser ce dernier modèle. 4. Il y a principalement 3 types de modèles pour modéliser la dépendance des temps de défaut : (a) les modèles factoriels (TR-GDR, page 465); (b) les modèles de copules (TR-GDR, page 458); (c) les modèles de contagion (TR-GDR, page 470); Expliquez aussi pourquoi il y a une dualité entre les deux premiers modèles (TR-GDR, pages 465470). 5. Le risque incrémental est défini aux pages 493-497 de TR-GDR. La notion de contribution en risque est présentée aux pages 497-506 de TR-GDR. 6. Les modèles de gestion de portefeuille sont des modèles d’allocation (TR-GDR, pages 505-516) alors que les modèles Raroc sont des modèles de tarification (TR-GDR, pages 516-526). Ceci explique que la couverture du risque de crédit par le biais de CDS concerne la gestion de portefeuille (TRGDR, page 516). Présentez aussi les différentes utilisations d’un modèle Raroc en tant qu’outil transactionnel (TR-GDR, pages 522-523). . On suppose ici que les positions ne changent pas d’une période à l’autre. . car ??1(0,99) = ???1(0,01). . On prend l’opposé de G?1(?) car on considère ici les minima et non les maxima. . En fait, il y a une seconde racine : On note que e(t) peut prendre des valeurs négatives pour t autour de cette racine. On ignore ce problème pour calculer les différentes quantités. . car la volatilité journalière est égale à : [7] ,6% ? = ?250 = 2% . ? = 99%. . L’excès de kurtosis a très peu d’effet. . Celui-ci est mesuré en points de volatilité. . On a : x1 = = xn = EAD . On fait bien sûr l’hypothèse que les notionnels et les maturités sont les mêmes. [13] . Remarquez alors que la méthode des distributions conditionnelles est équivalente à la méthode de factorisation de Cholesky dans le cas d’une copule gaussienne. |