Loi Gamma

Loi Gamma
Densité de probabilité


Fonction de répartition

Paramètres	$k>0\,$ réel $\theta >0\,$ réel
Support	$x\in [0,+\infty [$
Densité de probabilité	${\frac {x^{k-1}\mathrm {e} ^{-{\frac {x}{\theta }}}}{\Gamma (k)\theta ^{k}}}$
Fonction de répartition	${\frac {\gamma (k,x/\theta )}{\Gamma (k)}}$
Espérance	$k\theta \,$
Médiane	pas d'expression formelle
Mode	$(k-1)\theta \,$ pour $k\geq 1\,$
Variance	$k\theta ^{2}\,$
Asymétrie	${\frac {2}{\sqrt {k}}}$
Kurtosis normalisé	${\frac {6}{k}}$
Entropie	$k\theta +(1-k)\ln(\theta )+\ln(\Gamma (k))\,$ $+(1-k)\psi (k)\,$
Fonction génératrice des moments	$(1-\theta \,t)^{-k}$ pour $t<1/\theta$
Fonction caractéristique	$\left(1-\mathrm {i} \theta t\right)^{-k}\,$
modifier

En théorie des probabilités et en statistiques, une distribution Gamma ou loi Gamma est un type de loi de probabilité de variables aléatoires réelles positives. La famille des distributions Gamma inclut, entre autres, la loi du χ² et les distributions exponentielles et la distribution d'Erlang. Une distribution Gamma est caractérisée par deux paramètres $k$ et $θ$ et qui affectent respectivement la forme et l'échelle de la représentation graphique de sa fonction de densité.

Les distributions Gamma sont utilisées pour modéliser une grande variété de phénomènes, et tout particulièrement les phénomènes se déroulant au cours du temps où par essence, le temps écoulé est une grandeur réelle positive ; c'est le cas par exemple dans l'analyse de survie.

Définition[modifier | modifier le code]

Paramétrage avec la forme k et l'échelle $θ$ [modifier | modifier le code]

Soient $k$ et $θ$ deux réels strictement positifs. Une variable aléatoire $X$ suit une loi Gamma de paramètres $k$ et $θ$ , ce que l'on note aussi $X\,\sim \Gamma (k,\theta )$ (où $Γ$ est la majuscule de la lettre grecque gamma) si sa fonction de densité de probabilité peut se mettre sous la forme :

$f(x;k,\theta )={\frac {x^{k-1}\mathrm {e} ^{-{\frac {x}{\theta }}}}{\Gamma (k)\theta ^{k}}}$ ,

pour tout $x > 0$ . Dans l'expression ci-dessus, $Γ$ désigne la fonction Gamma d'Euler. Le paramètre k s'appelle le paramètre de forme, et le paramètre $θ$ est un paramètre d'échelle.

Paramétrage avec la forme $α$ et l'intensité $β$ [modifier | modifier le code]

Alternativement, la distribution Gamma peut être paramétrée à l'aide d'un paramètre de forme $α = k$ et d'un paramètre d'intensité $\beta =1/\theta$ (rate parameter) :

$f(x;\alpha ,\beta )=x^{\alpha -1}{\frac {\beta ^{\alpha }\,\mathrm {e} ^{-\beta \,x}}{\Gamma (\alpha )}}\ \mathrm {pour} \ x>0$ .

Les deux paramétrages sont répandus, selon le contexte. On note la même notation $\Gamma (k,\theta )$ et $\Gamma (\alpha ,\beta )$ pour la loi pour les deux paramétrages. La notation est ambigüe, mais elle dépend du paramétrage choisi.

Propriétés[modifier | modifier le code]

Moyenne et variance[modifier | modifier le code]

La moyenne (espérance) d'une distribution gamma est le produit des paramètres de forme et d'échelle :

\mu =k\theta =\alpha /\beta

La variance est donnée par :

\sigma ^{2}=k\theta ^{2}=\alpha /\beta ^{2}

L'inverse de la racine carré du paramètre de forme donne le coefficient de variation :

\sigma /\mu =1/{\sqrt {k}}=1/{\sqrt {\alpha }}

.

Coefficient d'asymétrie[modifier | modifier le code]

Le coefficient d'asymétrie d'une distribution gamma ne dépend que du paramètre de forme et vaut $2/{\sqrt {k}}.$

Moments[modifier | modifier le code]

Pour tout n entier, le n-ième moment vaut :

\mathrm {E} [X^{n}]=\theta ^{n}{\frac {\Gamma (k+n)}{\Gamma (k)}}=\theta ^{n}\prod _{i=1}^{n}(k+i-1)

.

Somme[modifier | modifier le code]

Si chaque $X i$ suit la loi $Γ(k i, θ)$ pour i = 1, 2,..., N, et si les variables aléatoires $X i$ sont indépendantes, alors :

$\sum _{i=1}^{N}X_{i}\sim \Gamma \left(\sum _{i=1}^{N}k_{i},\theta \right)$ .

Changement d'échelle[modifier | modifier le code]

Soit $X$ une variable aléatoire qui suit une loi gamma de paramètres de forme k et d'échelle θ. Alors pour tout $t > 0$ , la variable $tX$ est distribuée selon une loi $\Gamma (k,t\theta )$ de paramètre de forme k et d'échelle tθ. Dit autrement pour le paramétrage (α, β), si $X$ suit une loi gamma de paramètres de forme α et d'intensité β, alors $tX$ est distribuée selon une loi de paramètre de forme $\alpha$ et d'intensité β/t, que l'on note également $\Gamma (\alpha ,\beta /t)$ .

Lien avec les autres distributions[modifier | modifier le code]

Contraintes sur les paramètres[modifier | modifier le code]

Si $X\sim {\Gamma }(k=1,\theta =1/\lambda )\,$ , alors $X$ a une distribution exponentielle de paramètre $λ$ . En effet, pour $k=1$ , l'expression de la densité $f(x;k,\theta ={\frac {1}{\lambda }})={\frac {x^{k-1}\mathrm {e} ^{-x/\theta }}{\theta ^{k}\Gamma (k)}}={\frac {x^{k-1}\lambda ^{k}\mathrm {e} ^{-\lambda x}}{\Gamma (k)}}$ devient $\lambda \mathrm {e} ^{-\lambda x}$ . Il s'agit bien de la densité de probabilité de la distribution exponentielle de paramètre $λ$ .
Si $X\sim {\Gamma }(k=\nu /2,\theta =2)\,$ , alors $X$ est identique à une variable qui suit $χ 2 (ν)$ , la loi du χ² avec $ν$ degrés de liberté.
Si $k$ est un entier, la loi Gamma est une distribution d'Erlang.
Si $X^{2}\sim {\Gamma }(3/2,2a^{2})\,$ , alors $X$ a une distribution de Maxwell-Boltzmann avec comme paramètre $a$ .^{[Information douteuse]}

Autres manipulations[modifier | modifier le code]

Si X a une distribution Γ(k, θ), alors 1/X a une distribution loi gamma inverse, de paramètres k et $θ -1$ .
Si X et Y sont distribuées indépendamment selon des lois Γ(α, θ) et Γ(β, θ) respectivement, alors X / (X + Y) a une distribution beta de paramètres α et β.
Si X_i sont distribuées selon des lois Γ(α_i, θ) respectivement, alors le vecteur (X₁ / S, ..., X_n / S), où S = X₁ + ... + X_n, suit une distribution de Dirichlet de paramètres α₁, ..., α_n.
Pour k grand, la distribution Gamma converge vers une loi normale, de moyenne $\mu =(k-1)\theta$ et de variance $\sigma ^{2}=(k-1)\theta ^{2}$ . De plus, quels que soient k et θ, en fixant de cette manière les constantes $\mu$ et $\sigma$ , les densités de probabilité de la distribution Gamma Γ(k, θ) et de la loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$ ont alors deux points d'inflexion aux mêmes abscisses, à savoir $\mu +\sigma =\theta (k-1)+\theta {\sqrt {k-1}}$ et $\mu -\sigma =\theta (k-1)-\theta {\sqrt {k-1}}$ .

Propriété de concentration[modifier | modifier le code]

Si $X\sim \chi ^{2}(k)$ , alors^[1] pour tout $t>0$ , $\mathbb {P} \left(X\geq k+2{\sqrt {kt}}+2t\right)\leq \mathrm {e} ^{-t}$ et $\mathbb {P} \left(X\leq k-2{\sqrt {kt}}\right)\leq \mathrm {e} ^{-t}$ .

Généralisation[modifier | modifier le code]

Une loi Gamma généralisée a été définie avec un troisième paramètre^[2]: $f(x;k,\theta ,b)={\frac {bx^{bk-1}\mathrm {e} ^{-\left({\frac {x}{\theta }}\right)^{b}}}{\Gamma (k)\theta ^{bk}}}$ , afin de réunir dans une même famille la loi Gamma, la loi de Weibull et la loi exponentielle.

Références[modifier | modifier le code]

↑ (en) Nicolas Verzelen et Elisabeth Gassiat, « Adaptative estimation of high-dimensional signal to noise ratios », arXiv,‎ 16 mars 2017, p. 41 (lire en ligne)
↑ (en) E. W. Stacy, « A Generalization of the Gamma Distribution », Ann. Math. Statist., vol. 33, n^o 3,‎ septembre 1962, p. 1187-1192 (DOI 10.1214/aoms/1177704481)

Portail des probabilités et de la statistique

[1] (en) Nicolas Verzelen et Elisabeth Gassiat, « Adaptative estimation of high-dimensional signal to noise ratios », arXiv,‎ 16 mars 2017, p. 41 (lire en ligne)

[2] (en) E. W. Stacy, « A Generalization of the Gamma Distribution », Ann. Math. Statist., vol. 33, n^o 3,‎ septembre 1962, p. 1187-1192 (DOI 10.1214/aoms/1177704481)

[1]

[2]