Loi stable

La loi stable ou loi de Lévy tronquée, nommée d'après le mathématicien Paul Lévy, est une loi de probabilité utilisée en mathématiques, physique et analyse quantitative (finance de marché).

Variable aléatoire stable réelle

Définition

On dit qu'une variable aléatoire réelle $X$ est de loi stable si elle vérifie l'une des 3 propriétés équivalentes suivantes :

1. Pour tous réels strictement positifs $A$ et $B$, il existe un réel strictement positif $C$ et un réel $D$ tels que les variables aléatoires $A{X}_1+B{X}_2$ et $CX+D$ aient la même loi, où $X_1$ et $X_2$ sont des copies indépendantes de $X$.
2. Pour tout entier $n\ge 2$, il existe une constante strictement positive $C_n$ et un réel $D_n$ tels que les variables aléatoires $X_1+X_2+\dots +X_n$ et $C_{n}X+D_n$ aient la même loi, où $X_1,X_2,\dots ,X_n$ sont des copies indépendantes de $X$.
3. Il existe des réels $\alpha \in ]0,2]$, $\sigma \ge 0$, $\beta \in [-1,1]$ et $\mu \in ℝ$ telles que la fonction caractéristique de $X$ vérifie, pour tout $\theta \in ℝ$,

$𝔼\left[{\mathrm{e}}^{\mathrm{i}\theta X}\right]=\{\begin{array}{cc}\exp \{-{\sigma }^{\alpha }|\theta {|}^{\alpha }(1-\mathrm{i}\beta (\mathrm{sgn}\theta )\tan \frac{\pi \alpha }{2})+\mathrm{i}\mu \theta \}& \text{ si }\alpha \ne 1,\\ \\ \exp \{-\sigma |\theta |(1+\mathrm{i}\beta \frac{2}{\pi }(\mathrm{sgn}\theta )\ln |\sigma \theta |)+\mathrm{i}\mu \theta \}& \text{ si }\alpha =1,\end{array}$

où

$\mathrm{sgn}\theta =\{\begin{array}{cc}1& \text{ si }\theta >0,\\ 0& \text{ si }\theta =0,\\ -1& \text{ si }\theta <0.\end{array}$

Remarques :

• Les paramètres $\alpha \in ]0,2]$, $\sigma \ge 0$, $\beta \in [-1,1]$ et $\mu \in ℝ$ caractérisent la loi de $X$. On écrit alors $X\sim S_{\alpha }(\sigma ,\beta ,\mu )$.
• Le réel $\alpha$ dans $]0,2]$ est appelé paramètre de stabilité de $X$. Le réel positif $\sigma$ est appelé paramètre d'échelle de $X$.
• Les coefficients $A$, $B$ et $C$ sont liés par la relation $C^{\alpha }=A^{\alpha }+B^{\alpha }$.
• Pour tout $n\ge 2$, on a $C_n=n^{1/\alpha }$.

On dit qu'une variable aléatoire réelle $X$ est $\alpha$-stable si elle est stable et que son paramètre de stabilité est $\alpha$.

• Les lois 2-stables correspondent exactement aux lois normales. Pour ces lois, le paramètre $\beta$ n'a aucune influence. Plus précisément, la loi $S_2(\sigma ,\beta ,\mu )$ correspond à la loi normale $𝒩(\mu ,{\sigma }^2)$ de moyenne $\mu$ et de variance ${\sigma }^2$.

Propriétés des lois stables

• Si $X_1\sim S_{\alpha }({\sigma }_1,{\beta }_1,{\mu }_1)$ et $X_2\sim S_{\alpha }({\sigma }_2,{\beta }_2,{\mu }_2)$ sont indépendantes, alors $X_1+X_2\sim S_{\alpha }(\sigma ,\beta ,\mu )$ avec

$\sigma =({\sigma }_1^{\alpha }+{\sigma }_2^{\alpha }{)}^{1/\alpha },\beta =\frac{{\beta }_1{\sigma }_1^{\alpha }+{\beta }_2{\sigma }_2^{\alpha }}{{\sigma }_1^{\alpha }+{\sigma }_2^{\alpha }},\text{ et }\mu ={\mu }_1+{\mu }_2.$

• Si $X\sim S_{\alpha }(\sigma ,\beta ,\mu )$ et $a\in ℝ$, alors $X+a\sim S_{\alpha }(\sigma ,\beta ,\mu +a)$.
• Si $X\sim S_{\alpha }(\sigma ,\beta ,\mu )$ avec $\alpha \in ]0,2[$, alors

$\{\begin{array}{c}\underset{\lambda \to +\mathrm{\infty }}{\lim }{\lambda }^{\alpha }\mathbb{P}(X>\lambda )=C_{\alpha }\frac{1+\beta }{2}{\sigma }^{\alpha },\\ \underset{\lambda \to -\mathrm{\infty }}{\lim }{\lambda }^{\alpha }\mathbb{P}(X<\lambda )=C_{\alpha }\frac{1-\beta }{2}{\sigma }^{\alpha },\end{array}$

où

$C_{\alpha }={(\underset{0}{\overset{+\mathrm{\infty }}{\int }}{x}^{-\alpha }\sin x\mathrm{d}x)}^{-1}=\frac{2\mathrm{\Gamma }(\alpha )\sin (\pi \alpha /2)}{\pi }$

.

• Si $X\sim S_{\alpha }(\sigma ,\beta ,\mu )$ avec $\alpha \in ]0,2[$, alors

$\{\begin{array}{cc}𝔼[|X{|}^p]<+\mathrm{\infty }& \text{ si }p\in ]0,\alpha [,\\ 𝔼[|X{|}^p]=+\mathrm{\infty }& \text{ si }p\ge \alpha .\end{array}$

Cas symétrique

On dit que $X$ est de loi symétrique $\alpha$-stable si $X$ est $\alpha$-stable et que les variables aléatoires $X$ et $-X$ sont identiquement distribuées.

• $X$ est de loi symétrique $\alpha$-stable si, et seulement si, $X\sim S_{\alpha }(\sigma ,0,0)$. On note simplement dans ce cas $X\sim S_{\alpha }S(\sigma )$.
• $X$ est de loi symétrique $\alpha$-stable si, et seulement si, sa fonction caractéristique vérifie pour tout $\theta \in ℝ$ l'égalité $𝔼\left[{\mathrm{e}}^{\mathrm{i}\theta X}\right]={\mathrm{e}}^{-{\sigma }^{\alpha }|\theta {|}^{\alpha }}$, où $\sigma$ est le paramètre d'échelle de $X$.

Domaine d'attraction d'une loi stable

L'étude des lois stables vient, en fait, de l'étude de la convergence de sommes de variables aléatoires indépendantes et identiquement distribuées (i.i.d) normalisées de manière affine. Les lois stables sont alors les seules lois limites possibles.

Domaine d'attraction

Modèle:Théorème Une loi de probabilité $\nu$ sur $ℝ$ est dite dégénérée si elle et tout son poids sur un point, autrement dit, si il existe $x\in ℝ$ tel que $\nu (\{x\})=1$. Les lois dégénérées sont toutes stables. Une loi stable est dégénérée si et seulement si son paramètre $\sigma$ est nul.

Deux lois de probabilités $\nu$ et ${\nu }^{\prime }$ sur $ℝ$ sont dites de même type si on peut passer de l'une à l'autre par une transformation affine, autrement dit, si il existe $a\in ℝ$ et $b>0$ tels que pour tout $x\in ℝ$, $\nu (]-\mathrm{\infty },x])={\nu }^{\prime }(]-\mathrm{\infty },bx+a])$. Cela revient encore à dire que, pour toute variable aléatoire $X\sim \nu$ on a que $(bX+a)\sim {\nu }^{\prime }$. Avoir le même type définit une relation d'équivalence et si deux lois ont le même type alors la stabilité de l'une implique celle de l'autre. Plus précisément si $X\sim S_{\alpha }(\sigma ,\beta ,\mu )$ et si $Y=bX+a$ avec $a\in ℝ$ et $b>0$, alors $Y\sim S_{\alpha }(\sigma b,\beta ,b\mu +a)$. Deux lois stables non dégénérées ont le même type si et seulement si elles partagent les mêmes paramètres $\alpha <2$ et $\beta$ ou si elles sont toutes les deux 2-stables (normales). Toutes les lois dégénérées ont le même type. Modèle:Théorème Modèle:Théorème À noter que si $X$ a une variance finie, alors la condition du théorème est automatiquement vérifiée et on retrouve presque la conclusion du théorème central limite. Ce qui empêche de retrouver exactement le théorème central limite est que le théorème ci-dessus ne donne pas de suites $(A_n{)}_{n\ge 1}$ et $(B_n{)}_{n\ge 1}$ pour lesquels la convergence peut s'obtenir ni les paramètres $\mu$ et $\sigma$ de la loi normale limite associée. Le théorème central limite informe qu'en prenant $A_n=0$ et $B_n=\sqrt{n}$ il y a convergence en loi vers la loi normale $𝒩(\mu ,{\sigma }^2)$ où $\mu :=𝔼[X]$ et ${\sigma }^2:=\mathrm{V}\mathrm{a}\mathrm{r}(X)$.

Ce théorème peut s'appliquer à des lois de variance infinie. Par exemple si $X$ suit la loi de Pareto suivante : $\mathbb{P}(X>x)=x^{-3}$ pour tout $x\ge 1$, alors la condition du théorème est satisfaite, pourtant $X$ n'a pas une variance finie. Dans cet exemple précis, il est possible de trouver des suites $(A_n{)}_{n\ge 1}$ et $(B_n{)}_{n\ge 1}$ convenables pour que la convergence se fasse grâce au théorème suivant.Modèle:ThéorèmeA noter que si $a>2$ dans le théorème précédent, alors la variance de $X_1$ est finie et on retrouve exactement le théorème central limite. L'intérêt du théorème précédent par rapport au théorème central limite est donc l'étude du cas $a=2$.Modèle:Théorème

Vecteur aléatoire stable et variable aléatoire stable complexe

Vecteur aléatoire stable

On dit qu'un vecteur aléatoire $X=(X_1,\dots ,X_d)$ de ${ℝ}^d$ est de loi stable s'il vérifie une des 2 propriétés équivalentes suivantes :

1. Pour tous réels strictement positifs $A$ et $B$, il existe un réel strictement positif $C$ et un vecteur $D$ de ${ℝ}^d$ tels que les vecteurs aléatoires $A{X}^{(1)}+B{X}^{(2)}$ et $CX+D$ aient la même loi, où $X^{(1)}$ et $X^{(2)}$ sont des copies indépendantes de $X$.
2. Il existe une mesure finie $\mathrm{\Gamma }$ sur la sphère $S_d$ de ${ℝ}^d$ et un vecteur ${\mu }^0\in {ℝ}^d$ tels que la fonction caractéristique de $X$ vérifie, pour tout $({\theta }_1,\dots ,{\theta }_d)\in {ℝ}^d$,

$𝔼[\exp \left(\mathrm{i}\sum _{l=1}^d{\theta }_l{X}_l\right)]=\{\begin{array}{cc}\exp \{-\underset{S_d}{\int }|⟨\theta ,s⟩{|}^{\alpha }(1-\mathrm{i}\mathrm{sgn}(⟨\theta ,s⟩)\tan \frac{\pi \alpha }{2}))\mathrm{\Gamma }(\mathrm{d}s)+\mathrm{i}(\theta ,{\mu }^0)\}& \text{ si }\alpha \ne 1,\\ \\ \exp \{-\underset{S_d}{\int }|⟨\theta ,s⟩|(1+\mathrm{i}\frac{2}{\pi }\mathrm{sgn}(⟨\theta ,s⟩)\ln |⟨\theta ,s⟩|)\mathrm{\Gamma }(\mathrm{d}s)+\mathrm{i}(\theta ,{\mu }^0)\}& \text{ si }\alpha =1,\end{array}$

où

$⟨\cdot ,\cdot ⟩$

est le produit scalaire classique sur

${ℝ}^d$

.

Remarques :

• La paire $(\mathrm{\Gamma },{\mu }^0)$ est unique.
• Le réel $\alpha$ est appelé paramètre de stabilité de $X$.
• Les coefficients $A$, $B$ et $C$ sont liés par la relation $C^{\alpha }=A^{\alpha }+B^{\alpha }$.
• On dit que $X$ est de loi symétrique $\alpha$-stable si $X$ est $\alpha$-stable et que les variables aléatoires $X$ et $-X$ sont identiquement distribuées. Dans ce cas, sa fonction caractéristique est donnée, pour tout $({\theta }_1,\dots ,{\theta }_d)\in {ℝ}^d$, par $𝔼[\exp \left(\mathrm{i}\sum _{l=1}^d{\theta }_l{X}_l\right)]=\exp \{-\underset{S_d}{\int }|⟨\theta ,s⟩{|}^{\alpha }\mathrm{\Gamma }(\mathrm{d}s)\}$.

Propriétés des vecteurs aléatoires stables

• Si $X=(X_1,\dots ,X_d)$ est un vecteur $\alpha$-stable, alors, pour tous réels $b_1,\dots ,b_d$, la variable aléatoire réelle $\sum _{l=1}^d{b}_l{X}_l$ est $\alpha$-stable.
• Si $\alpha \in [1,2]$ et, pour tous réels $b_1,\dots ,b_d$, la variable aléatoire réelle $\sum _{l=1}^d{b}_l{X}_l$ est $\alpha$-stable, alors le vecteur $X=(X_1,\dots ,X_d)$ est $\alpha$-stable.
• Si, pour tous réels $b_1,\dots ,b_d$, la variable aléatoire réelle $\sum _{l=1}^d{b}_l{X}_l$ est symétrique $\alpha$-stable, alors le vecteur $X=(X_1,\dots ,X_d)$ est symétrique $\alpha$-stable.

Variable aléatoire stable complexe

On dit qu'une variable aléatoire complexe $Z=X+\mathrm{i}Y$ est de loi $\alpha$-stable, si le vecteur $(X,Y)$ de ${ℝ}^2$ est $\alpha$-stable.

On dit de plus que la loi de $Z$ est isotrope si, pour tout $\varphi \in [0,2\pi [$, les variables aléatoires ${\mathrm{e}}^{\mathrm{i}\varphi }Z$ et $Z$ sont identiquement distribuées. Dans ce cas, sa fonction caractéristique vérifie, pour tous complexes $\theta ={\theta }_1+\mathrm{i}{\theta }_2$, $𝔼\left[{\mathrm{e}}^{\mathrm{i}({\theta }_1{X}_1+{\theta }_1{X}_1)}\right]={\mathrm{e}}^{-{\sigma }^{\alpha }|\theta {|}^{\alpha }}$, où $\sigma$ est un réel positif appelé paramètre d'échelle de $Z$.

Représentation en série de LePage

Cas symétrique réel

Soit $\alpha \in ]0,2[$. On pose $a(\alpha )={(\underset{0}{\overset{+\mathrm{\infty }}{\int }}{x}^{-\alpha }\sin (x)\mathrm{d}x)}^{-1/\alpha }$. Soit $\{{\mathrm{\Gamma }}_m:m\in \mathbb{N}\}$ et $\{Z_m:m\in \mathbb{N}\}$ deux processus mutuellement indépendants de variables aléatoires définis sur le même espace de probabilité $(\mathrm{\Omega },𝒢,\mathbb{P})$ satisfaisant aux propriétés suivantes :

1. Les ${\mathrm{\Gamma }}_m$, $m\in \mathbb{N}$, sont les temps d'arrivée d'un processus de Poisson d'intensité 1 ; c'est-à-dire, pour tous $m\in \mathbb{N}$, on a ${\mathrm{\Gamma }}_m=\sum _{n=1}^m{\nu }_n$, où $({\nu }_n{)}_{n\in \mathbb{N}}$ est une suite de variables aléatoires exponentielles de paramètre 1 indépendantes.
2. Les $Z_m$, $m\in \mathbb{N}$ sont des variables aléatoires réelles, symétriques, indépendantes, identiquement distribuées et vérifiant $𝔼[|Z_m{|}^{\alpha }]<+\mathrm{\infty }$.

Alors la série $\sum _{m=1}^{+\mathrm{\infty }}{Z}_m{\mathrm{\Gamma }}_m^{-1/\alpha }$ converge presque sûrement. De plus, elle est de loi symétrique $\alpha$-stable et son paramètre d'échelle $\sigma$ vérifie $\sigma =a(\alpha {)}^{-1}{(𝔼[|Z_1{|}^{\alpha }])}^{1/\alpha }$.

Cas isotrope complexe

Soit $\alpha \in ]0,2[$. On pose $a(\alpha )={(\underset{0}{\overset{+\mathrm{\infty }}{\int }}{x}^{-\alpha }\sin (x)\mathrm{d}x)}^{-1/\alpha }$. Soit $\{{\mathrm{\Gamma }}_m:m\in \mathbb{N}\}$ et $\{Z_m:m\in \mathbb{N}\}$ deux processus mutuellement indépendants de variables aléatoires définis sur le même espace de probabilité $(\mathrm{\Omega },𝒢,\mathbb{P})$ satisfaisant aux propriétés suivantes :

1. Les ${\mathrm{\Gamma }}_m$, $m\in \mathbb{N}$, sont les temps d'arrivée d'un processus de Poisson d'intensité 1.
2. Les $Z_m$, $m\in \mathbb{N}$, sont des variables aléatoires complexes, isotropes, indépendantes, identiquement distribuées et vérifiant $𝔼[|\text{Re}(Z_m){|}^{\alpha }]<+\mathrm{\infty }$, où $\text{Re}(Z_m)$ désigne la partie réelle de $Z_m$.

Alors la série $\sum _{m=1}^{+\mathrm{\infty }}{Z}_m{\mathrm{\Gamma }}_m^{-1/\alpha }$ converge presque sûrement. De plus, elle est de loi isotrope $\alpha$-stable et son paramètre d'échelle $\sigma$ vérifie $\sigma =a(\alpha {)}^{-1}{(𝔼[|{\text{Re}(Z_1)|}^{\alpha }])}^{1/\alpha }$.

Liens avec d'autres lois

Elle a pour cas particuliers :

• La loi de Lévy (paramètres α=1/2 et beta=1), définie par une formule analytique explicite.
• La loi normale (paramètre α=2), définie par une formule analytique explicite.
• La loi de Cauchy (paramètre α=1), définie par une formule analytique explicite.

Gnedenko et Kolmogorov ont établi une généralisation du théorème central limite selon laquelle la somme de variables aléatoires ayant des queues de loi décroissantes selon 1/|x|Modèle:Exp avec 0 < α < 2 (ayant donc une variance infinie) tend vers une loi stable de paramètre α.

Références

Modèle:Palette Modèle:Portail