Loi de Laplace (probabilités)

Modèle:Infobox Distribution statistiques En théorie des probabilités et en statistiques, la loi (distribution) de Laplace est une densité de probabilité continue, portant le nom de Pierre-Simon de Laplace. On la connaît aussi sous le nom de loi double exponentielle, car sa densité peut être vue comme l'association des densités de deux lois exponentielles, accolées dos à dos. La loi de Laplace s'obtient aussi comme résultat de la différence de deux variables exponentielles indépendantes.

Caractérisation

Densité de probabilité

Une variable aléatoire suit une loi de Laplace(μ, b) si sa densité de probabilité est

$f(x|\mu ,b)=\frac{1}{2b}\exp (-\frac{|x-\mu |}{b})$

$=\frac{1}{2b}\{\begin{array}{cc}\exp (-\frac{\mu -x}{b})& \text{si }x<\mu \\ \exp (-\frac{x-\mu }{b})& \text{si }x\ge \mu \end{array}$

Le réel μ est un paramètre de position et b > 0 un paramètre d'échelle. Si μ = 0 et b = 1, la loi de Laplace est dite standard et sa restriction à la demi-droite réelle positive est la loi exponentielle de paramètre 1/2.

La densité rappelle aussi celle de la loi normale ; toutefois, tandis que la loi normale est exprimée en termes de la différence au carré $(x-\mu {)}^2$, la loi de Laplace fait intervenir la différence absolue $|x-\mu |$. La loi de Laplace présente alors des queues plus épaisses que la loi normale.

Fonction de répartition

La densité de la loi de Laplace s'intègre aisément grâce à la présence de la valeur absolue. Sa fonction de répartition est :

$F(x)$	$=\underset{-\mathrm{\infty }}{\overset{x}{\int }}f(u)\mathrm{d}u$
	$=\{\begin{array}{cc}\frac{1}{2}\exp (-\frac{\mu -x}{b})& \text{si }x<\mu \\ 1-\frac{1}{2}\exp (-\frac{x-\mu }{b})& \text{si }x\ge \mu \end{array}$
	$=0,5[1+\mathrm{sgn}(x-\mu )(1-\exp (-|x-\mu |/b))].$

La réciproque de la fonction de répartition est

$F^{-1}(p)=\mu -b\mathrm{sgn}(p-0,5)\ln (1-2|p-0,5|).$

Tirer une variable selon la loi de Laplace

Étant donné une variable U, tirée selon une loi uniforme continue dans l'intervalle [-1/2, 1/2], la variable suivante

$X=\mu -b\mathrm{sgn}(U)\ln (1-2|U|)$

est distribuée selon la loi de Laplace de paramètres μ et b. Ce résultat provient de l'expression de l'inverse de la fonction de répartition et de la méthode de la transformée inverse.

Une variable Laplace(0, b) peut aussi se générer comme la différence de deux variables exponentielles, de paramètre 1/b, indépendantes. De même, une loi Laplace(0, 1) peut s'obtenir en considérant le logarithme du ratio de deux variables uniformes indépendantes.

Estimation des paramètres

Étant donné un échantillon de N variables iid xModèle:Ind, xModèle:Ind..., xModèle:Ind, un estimateur $\hat{\mu }$ de $\mu$ est la médiane empirique, et un estimateur par maximum de vraisemblance de b est

$\hat{b}=\frac{1}{N}\sum _{i=1}^N|x_i-\hat{\mu }|.$

Moments

Les moments centrés d'ordre r peuvent s'exprimer à partir des moments ordinaires :

${{\mu }_r}^{\prime }=\frac{1}{2}\sum _{k=0}^r(\genfrac{}{}{0ex}{}{r}{k})k!b^k{\mu }^{(r-k)}\{1+(-1{)}^k\}$

Lois associées

• Si $X\sim \mathrm{L}\mathrm{a}\mathrm{p}\mathrm{l}\mathrm{a}\mathrm{c}\mathrm{e}(0,b)$ alors $|X|\sim \mathrm{E}\mathrm{x}\mathrm{p}\mathrm{o}\mathrm{n}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{i}\mathrm{e}\mathrm{l}\mathrm{l}\mathrm{e}(b^{-1})$ est une loi exponentielle;
• Si $X\sim \mathrm{E}\mathrm{x}\mathrm{p}\mathrm{o}\mathrm{n}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{i}\mathrm{e}\mathrm{l}\mathrm{l}\mathrm{e}(\lambda )$ et Y indépendante de $X$ suit une loi de Rademacher, alors $XY\sim \mathrm{L}\mathrm{a}\mathrm{p}\mathrm{l}\mathrm{a}\mathrm{c}\mathrm{e}(0,{\lambda }^{-1})$;
• Si $X_1\sim \mathrm{E}\mathrm{x}\mathrm{p}\mathrm{o}\mathrm{n}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{i}\mathrm{e}\mathrm{l}\mathrm{l}\mathrm{e}({\lambda }_1)$ et $X_2\sim \mathrm{E}\mathrm{x}\mathrm{p}\mathrm{o}\mathrm{n}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{i}\mathrm{e}\mathrm{l}\mathrm{l}\mathrm{e}({\lambda }_2)$ indépendantes de $X_1$, alors ${\lambda }_1{X}_1-{\lambda }_2{X}_2\sim \mathrm{L}\mathrm{a}\mathrm{p}\mathrm{l}\mathrm{a}\mathrm{c}\mathrm{e}(0,1)$.