La série A de Gram-Charlier (nommée en l'honneur de Jørgen Pedersen Gram et Carl Charlier) et la série d'Edgeworth (nommée en l'honneur de Francis Ysidro Edgeworth) sont des séries qui se rapprochent d'une distribution de probabilité exprimée à partir de ses cumulants [1]. Les séries sont identiques, mais l'arrangement des termes (et donc la précision de la troncature de la série) diffère[2]. Le principe de ces développements est d'écrire la fonction caractéristique de la distribution dont la fonction de densité de probabilité f doit être approchée en fonction de la fonction caractéristique d'une distribution avec des propriétés connues et appropriées, et de récupérer f par une transformée de Fourier inverse.

Série A de Gram-Charlier

modifier

On considère une variable aléatoire continue. On note   la fonction caractéristique de sa distribution dont la fonction de densité est f, et   ses cumulants. On développe en termes de distribution connue avec la fonction de densité de probabilité ψ, sa fonction caractéristique  , et ses cumulants   . La densité ψ est généralement choisie comme étant celle de la distribution normale, mais d'autres choix sont également possibles. Par la définition des cumulants, on a (voir Wallace, 1958) [3]:

  et
 

qui donne l'identité formelle suivante :

 

Par les propriétés de la transformée de Fourier,   est la transformée de Fourier de  , où D est l'opérateur différentiel par rapport à x . Ainsi, après avoir changé   avec   des deux côtés de l'équation, on trouve pour f le développement formel

 

Si ψ est choisi comme la densité normale

 

avec la moyenne et la variance données par f, c'est-à-dire comme moyenne   et comme variance  , alors le développement devient

 

puisque   pour tout r > 2, car les cumulants d'ordres supérieurs de la distribution normale sont nuls. En développant les termes exponentiels et réunissant les termes selon l'ordre des dérivées, on arrive à la série A de Gram-Charlier. Un tel développement peut être écrit de manière compacte avec les polynômes de Bell ainsi :

 

Puisque la dérivée n-ième de la fonction gaussienne   est donnée en termes de polynôme d'Hermite par

 

on obtient l'expression finale de la série A de Gram-Charlier comme

 

L'intégration de la série donne la fonction de répartition :

 

  est la fonction de répartition de la loi normale.

Si on n'inclut que les deux premiers termes de correction à la loi normale, il vient :

 

avec   et   .

On notera que cette expression n'est pas assurée d'être positive et ne définit donc pas une loi de probabilité valide. La série A de Gram – Charlier diverge dans de nombreux cas classiques - elle ne converge que si   décroit plus vite que   à l'infini (Cramér 1957). Lorsqu'elle ne converge pas, la série n'est pas non plus un véritable développement asymptotique, car il n'est pas possible d'estimer l'erreur du développement. Pour cette raison, la série d'Edgeworth est généralement préférée à la série A de Gram-Charlier.

Série d'Edgeworth

modifier

Edgeworth a établi un développement similaire en tant qu'amélioration du théorème central limite[4]. L'avantage de la série d'Edgeworth est que l'erreur est contrôlée, de sorte qu'il s'agit d'un véritable développement asymptotique.

Soient   une suite de variables aléatoires indépendantes et identiquement distribuées de moyenne   et variance  , et soient   leurs sommes standardisées :

 

On note   les fonctions de répartition des variables  . Alors par le théorème central limite,

 

pour chaque  , tant que la moyenne et la variance sont finies.

La standardisation de   assure que les deux premiers cumulants de   valent   et   On suppose maintenant que, en plus d'avoir une moyenne   et une variance   finies, les variables aléatoires iid   ont des cumulants plus élevés   bien définis. A partir des propriétés d'additivité et d'homogénéité des cumulants, les cumulants de   en termes de cumulants de   sont pour   ,

 

Si on développe l'expression formelle de la fonction caractéristique   de   en termes de distribution normale standard, c'est-à-dire si l'on pose

 

alors les différences entre cumulants dans le développement sont

 
 
 

La série A de Gram-Charlier pour la fonction de densité de   est maintenant

 

La série d'Edgeworth est développée de manière similaire à la série A de Gram-Charlier, sauf que maintenant les termes sont regroupés en fonction des puissances de n. Les coefficients du terme n-m/2 peuvent être obtenus en rassemblant les monômes des polynômes de Bell correspondant aux partitions entières de m. Ainsi, on obtient la fonction caractéristique :

 

  est un polynôme de degré  . Ainsi, après transformée de Fourier inverse, la fonction de densité   s'exprime comme suit :

 

De même, en intégrant la série, on obtient la fonction de répartition :

 

On peut écrire explicitement le polynôme   comme

 

où la sommation est sur toutes les partitions entières de m telles que   et   et  

Par exemple, si m = 3, alors il y a trois façons de partitionner ce nombre : 1 + 1 + 1 = 2 + 1 = 3. Il faut donc examiner trois cas :

  • 1 + 1 + 1 = 1 · k1, ce qui donne k1 = 3, l1 = 3 et s = 9.
  • 1 + 2 = 1 · k1 + 2 · k2, ce qui donne k1 = 1, k2 = 1, l1 = 3, l2 = 4 et s = 7.
  • 3 = 3 · k3, ce qui donne k3 = 1, l3 = 5 et s = 5.

Ainsi, le polynôme recherché est

 

Les cinq premiers termes du développement sont

 

Ici, ϕ(j)(x) est la dérivée j -ième de ϕ(·) au point x. En rappelant que les dérivées de la densité de la distribution normale sont liées à la densité normale par  , (où   est le polynôme d'Hermite d'ordre n), ceci explique les représentations alternatives en termes de fonction de densité. Blinnikov et Moessner (1998) ont donné un algorithme simple pour calculer les termes d'ordre supérieur du développement.

Il faut noter que dans le cas d'une distribution de réseau (qui a des valeurs discrètes), le développement d'Edgeworth doit être ajusté pour tenir compte des sauts discontinus entre les points du réseau[5].

Illustration : densité de la moyenne d'échantillon de trois distributions χ2

modifier
 
Densité de la moyenne de l'échantillon de trois variables χ2. Le graphique compare la densité réelle, l'approximation normale et deux développements d'Edgeworth.

On prend l'exemple de trois variables iid suivant toutes la loi du χ2 à deux degrés de liberté :

  et la moyenne de l'échantillon   .

On peut utiliser plusieurs distributions pour approcher   :

  • la distribution exacte, qui suit une distribution gamma :   .
  • la distribution normale asymptotique :   .
  • deux développements d'Edgeworth, de degrés 2 et 3.

Discussion des résultats

modifier
  • Pour les échantillons finis, une série d'Edgeworth n'est pas garantie d'être une loi de probabilité appropriée car les valeurs de la fonction de répartition à certains points peuvent aller au-delà de l'intervalle   .
  • Ils garantissent (asymptotiquement) des erreurs absolues, mais les erreurs relatives peuvent être facilement évaluées en comparant le terme principal d'Edgeworth dans le reste avec le terme principal global [6].

Voir aussi

modifier

Références

modifier
  1. (en) A. Stuart et M.G. Kendall, The advanced theory of statistics, Hafner Publishing Company, .
  2. (en) J.E. Kolassa, Series approximation methods in statistics, vol. 88, Springer Science & Business Media, .
  3. D.L. Wallace, « Asymptotic Approximations to Distributions », Annals of Mathematical Statistics, vol. 29, no 3,‎ , p. 635–654 (DOI 10.1214/aoms/1177706528, JSTOR 2237255)
  4. Hall, P. (2013). The bootstrap and Edgeworth expansion. Springer Science & Business Media.
  5. John E. Kolassa et Peter McCullagh, « Edgeworth series for lattice distributions », Annals of Statistics, vol. 18, no 2,‎ , p. 981–985 (DOI 10.1214/aos/1176347637, JSTOR 2242145)
  6. John E. Kolassa, Series approximation methods in statistics, 3rd, (ISBN 0387322272)

Bibliographie

modifier
  • (en) Harald Cramér, Mathematical Methods of Statistics, Princeton, Princeton University Press, .
  • (en) D.L. Wallace, « Asymptotic approximations to distributions », Annals of Mathematical Statistics, vol. 29, no 3,‎ , p. 635–654 (DOI 10.1214/aoms/1177706528  )
  • (en) M. Kendall et A. Stuart, The advanced theory of statistics, vol. 1: Distribution theory, New York, Macmillan, .
  • (en) Peter McCullagh, Tensor Methods in Statistics, London, Chapman and Hall, .
  • (en) David Cox et Ole Barndorff-Nielsen, Asymptotic Techniques for Use in Statistics, London, Chapman and Hall, .
  • (en) P. Hall, The Bootstrap and Edgeworth Expansion, New York, Springer, .
  • (en) « Edgeworth series », dans Michiel Hazewinkel, Encyclopædia of Mathematics, Springer, (ISBN 978-1556080104, lire en ligne)
  • (en) S. Blinnikov et R. Moessner, « Expansions for nearly Gaussian distributions », Astronomy and Astrophysics Supplement Series, vol. 130,‎ , p. 193–205 (DOI 10.1051/aas:1998221, Bibcode 1998A&AS..130..193B, arXiv astro-ph/9711239, lire en ligne)
  • (en) Douglas Martin et Rohit Arora, « Inefficiency and bias of modified value-at-risk and expected shortfall », Journal of Risk, vol. 19, no 6,‎ , p. 59–84 (DOI 10.21314/JOR.2017.365)
  • (en) J.E. Kolassa, Series Approximation Methods in Statistics, New York, Springer, .