Tag Archives: FDR

Fonction de répartition empirique

Ce matin, en cours, on avait utilisé une série d’observations pour décrire l’estimation nonparamétrique de la fonction de répartition

> Davis=read.table(
"http://socserv.socsci.mcmaster.ca/jfox/Books/Applied-Regression-2E/datasets/Davis.txt")

Comme il y une faute de saisie dans la base, on corrige, à la main,

> Davis[12,c(2,3)]=Davis[12,c(3,2)]

Notre variable d’intérêt sera la taille d’un individu,

> X=Davis$height
> X[1:10]
[1] 182 161 161 177 157 170 167 186 178 171

Pour tracer la fonction de répartition empirique,

https://latex.codecogs.com/gif.latex?\widehat{F}_n(x)=\frac{1}{n}\sum_{i=1}^n\boldsymbol{1}(\{x_i\leq%20x\})

on peut utiliser le code suivant

> n=length(X)
> x=seq(140,200)
> ggplot(data=data.frame(x=x), aes(x = x))+
+ stat_function(fun=Vectorize(function(x) sum(X<=x)/n),
+ geom="step", col = "red", size = 1) + 
+ ylab("Fonction de Repartition")

Statistiques, STT2700, cours 05/01

Petite précision suite au cours de ce matin: le cours du vendredi a lieu de 10.30 à 11.30 (il y avait une faute de frappe sur le plan de cours ici). Pour les billets liés au cours, il faut aller .

Côté travaux pratiques, ils commencent ce vendredi, et l’énoncé est en ligne ici (j’ai repris les exercices de la session d’automne).

Enfin, petit complément par rapport à une question posée à la fin du cours. Un élève m’avait demandé s’il ne fallait pas supposer la fonction de répartition dérivable pour définir une densité. Une variable aléatoire http://freakonometrics.blog.free.fr/public/maths/montone-deriv1.png est continue si sa fonction de répartition http://freakonometrics.blog.free.fr/public/maths/montone2.png est continue. http://freakonometrics.blog.free.fr/public/maths/montone2.png est par ailleurs croissante. Aussi, on peut montrer que http://freakonometrics.blog.free.fr/public/maths/montone2.png est dérivable presque partout. Comme le montrent Douglas S. Kurtz, Jaroslav Kurzweil et Charles Swartz dans Theories of integration: the integrals of Riemann, Lebesgue, Henstock,

ou plus particulièrement, on a le résultat suivant

On retrouve le même résultat, appelé théorème de Lebesgue dans l’ouvrage de Yeh, Real analysis: theory of measure and integration

Maintenant pour associer densité et dérivée, il faut faire un peu attention (cf le commentaire ici). Maintenant, en faisant quelques raccourcis (supposer la fonction à variations finies par exemple), on a le résultat suivant. Si http://freakonometrics.blog.free.fr/public/maths/montone2.png est supposée continue, la densité est une fonction http://freakonometrics.blog.free.fr/public/maths/montone4.png qui vérifie

http://freakonometrics.blog.free.fr/public/maths/montone3.png

En fait, comme le montrent les théorèmes ci-dessus, cette densité est unique, sauf aux points de non dérivabilité de http://freakonometrics.blog.free.fr/public/maths/montone2.png, qui sont de mesure nulle, autrement dit http://freakonometrics.blog.free.fr/public/maths/montone4.png est unique presque partout. Et, pour les points où http://freakonometrics.blog.free.fr/public/maths/montone2.png est dérivable,

http://freakonometrics.blog.free.fr/public/maths/montone5.png

et ailleurs… ? peu importe, en fait….
Considérons la fonction de répartition suivante, continue et dérivable partout, sauf en 0, et à variations finies.

La densité est alors, pour http://freakonometrics.blog.free.fr/public/maths/monontone6.png

http://freakonometrics.blog.free.fr/public/maths/monotone7.png

c’est à dire 1/2 sur http://freakonometrics.blog.free.fr/public/maths/monontone11.png et 1/4 sur http://freakonometrics.blog.free.fr/public/maths/monotone12.png,

et peu importe la valeur (finie) donnée en 0: on peut définir http://freakonometrics.blog.free.fr/public/maths/monotone13.png la dérivée à droite, et http://freakonometrics.blog.free.fr/public/maths/monotone15.png la dérivée à gauche, par exemple (si ces dérivées existent), de telle sorte que

http://freakonometrics.blog.free.fr/public/maths/monotone20.png

Bref, la densité est définie pour une loi continue (sans autre hypothèse), presque toujours. Pour aller plus loin, je renvoie ici par exemple….

 

Quantile de sommes, ou somme de quantiles ?

Je reprends ici un commentaire que j’ai longtemps entendu dans mon expérience d’actuaire dans une vie antérieure (et qu’on peut lire – entre les lignes le plus souvent – dans certaines réactions de risk managers). Les quantiles sont la base des mesures de risques, en particulier via la VaR en finance, mais quand j’étais actuaire, on calculait des 90/10 pour désigner des quantiles à 90% des pertes. Donc tout ça n’est pas spécialement nouveau. Je rappelle que la fonction quantile est simplement  où

Les quantiles sont calculés par type de risque (en gestion des risques), et forcément à la fin, on souhaite un quantile de la somme des risques, et il a longtemps été coutume de prendre la somme des quantiles. Je vais donc discuter un peu le résultat suivant,

  • L’égalité n’est pas vraie dans le cas indépendant

Suite à une espère de déformation professionnelle (à cause de la variance), beaucoup de monde pense que

mais ce résultat est faux…. Un exemple est d’ailleurs donné ci-dessous dans le cas de variables exponentielles.

  • L’égalité est vraie dans le cas comonotone…

En fait, le résultat précédant est vrai dans le cas comonotone. La comonotonie correspond au cas de dépendance parfaite positive, de corrélation maximale (si la corrélation existe), i.e.  avec  une fonction strictement croissante. Alors

  • … mais ce n’est pas un “worst case

Seconde bizarrerie, la comonotonie n’est pas le pire des cas. La comonotonie est souvent envisagée comme un worst case scenario mais il n’en est rien (cette idée se retrouve dans les QIS 3 par exemple, ici).
Formellement, si  (la classe de Fréchet), alors

et

En passant en dimension 2 (cas où la borne inférieure est effectivement une vraie fonction de répartition), si et  désignent respectivement des version anticomonotone et comonotone de , alors on peut se demander, en considérant une mesure de risque quelconque  si

En fait, Andre Tchen a montré en 1980 un résultat qui montre que la comonotonie est effectivement un worst case, mais seulement dans certains cas très limités.
Si est une fonction supermodulaire, i.e.

pour tout  et . Alors dans ce cas, pour tout ,

(la preuve se trouve ici). Je renvoie aux papiers de Michel Denuit, de Jan Dhaene ou de Marc Goovarts (ici ou ) pour des applications en tarification sur deux têtes, par exemple, ou sur les primes stop-loss en réassurance

  • Borner (numériquement) le quantile d’une somme

Commençons par un exemple, avec deux lois exponentielles, car c’est le plus simple, et surtout on peut calculer explicitement les bornes (dans le cas général, on se contentera de méthode numérique comme l’avait fait Khalil ici).
Si  et , où, pour être plus précis,

l’encadrement

est valide pour tout  quelle que soit la structure de dépendance entre et , où

Aussi, en prenant les inverses (car ces fonctions sont strictement croissante), on obtient une inégalité en terme de quantile, ou de Value-at-Risk,

pour tout .
Rappelons que dans le cas de variables indépendantes  alors que dans le cas comonotone . La figure ci-dessous montre les valeurs possible pour  la fonction de répartition de la somme, avec le cas indépendant en bleu, et le cas comotonone en rouge.

On peut aussi visualiser les valeurs possibles pour les quantiles de la somme. Dans un cas général, Ce problème a été traité par Makarov en 1981, et par Frank, Nelsen et Schweizer en 1997 (ici), tout en restant en dimension 2. En fait, ils ont travaillé sur les bornes de la fonction de répartition d’une somme mais c’est pareil, comme on l’a vu juste auparavant. Ils ont montré le résultat suivant: si  alors pour tout

désigne la copule anticomonotone (qui est une copule en dimension 2),

en posant , et

On peut retrouver des résultats plus fin en creusant dans la direction de l’arithmetic probability, par exemple avec la thèse de Williamson (ici, cette thèse est remarquable, et trop peu citée, malheureusement). Je ne parle que de la dimension 2, l’extension en dimension supérieure est délicate (je renvoie aux travaux de Paul Embrechts sur le sujet, ici ou , par exemple).