Loi normale

En théorie des probabilités et en statistique, les lois normales sont parmi les lois de probabilité les plus utilisées pour modéliser des phénomènes naturels issus de plusieurs événements aléatoires. Elles sont en lien avec de nombreux objets mathématiques dont le mouvement brownien, le bruit blanc gaussien ou d'autres lois de probabilité. Elles sont également appelées lois gaussiennes, lois de Gauss ou lois de Laplace-Gauss des noms de Laplace (1749-1827) et Gauss (1777-1855), deux mathématiciens, astronomes et physiciens qui les ont étudiées.

Loi normale
Densité de probabilité La courbe rouge représente la fonction $\varphi$, densité de probabilité de la loi normale centrée réduite.
Fonction de répartition La courbe rouge représente la fonction $\Phi$, fonction de répartition de la loi normale centrée réduite.
Paramètres	$\mu$, espérance (nombre réel) $\sigma ^{2}>0$, variance (nombre réel)
Support	$\mathbb {R}$
Densité de probabilité	${\frac {1}{\sigma {\sqrt {2\pi }}}}\;\exp \left(-{\frac {\left(x-\mu \right)^{2}}{2\sigma ^{2}}}\right)\!~$
Fonction de répartition	${\frac {1}{2}}\left(1+\mathrm {erf} \left({\frac {x-\mu }{\sigma {\sqrt {2}}}}\right)\right)\!~$
Espérance	$\mu$
Médiane	$\mu$
Mode	$\mu$
Variance	$\sigma ^{2}$
Asymétrie	0
Kurtosis normalisé	0
Entropie	$\ln \left(\sigma {\sqrt {2\,\pi \,{\rm {e}}}}\right)\!~$
Fonction génératrice des moments	$\exp \left(\mu \,t+{\frac {\sigma ^{2}t^{2}}{2}}\right)$
Fonction caractéristique	$\exp \left(\mu \,{\rm {i}}\,t-{\frac {\sigma ^{2}t^{2}}{2}}\right)$

Plus formellement, une loi normale est une loi de probabilité absolument continue qui dépend de deux paramètres : son espérance, un nombre réel noté μ, et son écart type, un nombre réel positif noté σ. La densité de probabilité de la loi normale d'espérance μ et d'écart type σ est donnée par :

$f(x)={\frac {1}{\sigma {\sqrt {2\pi }}}}\operatorname {e} ^{-{\frac {1}{2}}\left({\frac {x-\mu }{\sigma }}\right)^{2}}$.

La courbe de cette densité est appelée courbe de Gauss ou courbe en cloche, entre autres. C'est la représentation la plus connue de ces lois. Lorsqu'une variable aléatoire X suit une loi normale, elle est dite gaussienne ou normale et il est habituel d'utiliser la notation avec la variance σ² :

$X\sim {\mathcal {N}}(\mu ,\sigma ^{2})$.

La loi normale de moyenne nulle et d'écart type unitaire, ${\mathcal {N}}(0,1)$, est appelée loi normale centrée réduite ou loi normale standard.

Parmi les lois de probabilité, les lois normales prennent une place particulière grâce au théorème central limite. En effet, elles correspondent au comportement, sous certaines conditions, d'une suite d'expériences aléatoires similaires et indépendantes lorsque le nombre d'expériences est très élevé. Grâce à cette propriété, une loi normale permet d'approcher d'autres lois et ainsi de modéliser de nombreuses études scientifiques comme des mesures d'erreurs ou des tests statistiques, en utilisant par exemple les tables de la loi normale centrée réduite.

Définition et explications informelles

Les diagrammes en bâtons représentent les lois discrètes de la somme de 1, 2, 3, 4 ou 5 dés. La courbe noire est la densité d'une loi normale vue comme limite des diagrammes en bâtons.

Les lois de probabilité permettent de décrire de manière théorique le caractère aléatoire d'une expérience qui est considérée comme aléatoire. Les lois normales en sont des cas particuliers. La manière historique de l'aborder est par approximation.

Lorsque le résultat de cette expérience aléatoire est à valeurs discrètes, par exemple la somme du lancer de deux dés vaut 2, 3… ou 12, une loi dite discrète modélise l'expérience. Les probabilités d'apparition de chaque valeur peuvent être représentées par des diagrammes en bâtons ou histogrammes (voir la figure ci-contre). Plusieurs scientifiques (voir Histoire de la loi normale) se sont intéressés à la réalisation d'un grand nombre d'expériences et au comportement de la loi de probabilité associée. Il apparaît que les fréquences d'apparition des valeurs possibles sont de plus en plus « lissées » (voir la figure ci-contre). Il existe une certaine répartition autour d'une valeur centrale ; ces probabilités peuvent être alors représentées par une courbe de Gauss ou courbe en cloche obtenue par calcul ou par expérience. Cette courbe est celle de la densité de probabilité d'une loi normale. Le rôle central de ces lois de probabilité vient du fait qu'elles sont la limite d'un grand nombre de lois de probabilité définies à partir de sommes, comme le montre le théorème central limite.

Une autre manière visuelle de voir apparaître cette courbe est réalisée par la planche de Galton. Des billes sont lâchées en haut de la planche ; à chaque étage, elles ont deux possibilités : aller à droite ou aller à gauche. Après plusieurs étages, elles ont donc eu plusieurs choix aléatoires. Lorsque le nombre de billes est grand, la répartition des billes suivant leur position est approximativement une loi normale.

Comme pour toute loi de probabilité, plusieurs définitions équivalentes des lois normales existent : par leur densité de probabilité (la courbe de Gauss), par leur fonction de répartition, par leur fonction caractéristique, etc. Une loi normale dépend de deux paramètres : le premier donne la moyenne, c'est-à-dire la valeur « centrale » (ou « médiane ») des valeurs possibles (par exemple, la moyenne de la somme de deux dés est 7) ; le deuxième paramètre renseigne sur la dispersion des valeurs autour de cette valeur centrale : plus ce paramètre est faible, plus les valeurs proches de la valeur centrale auront une forte probabilité d'apparaître. Beaucoup de grandeurs physiques peuvent être représentées par ces deux paramètres.

Lors de l'étude statistique d'une série d'observations d'une même grandeur, la moyenne des valeurs observées peut être considérée comme une variable aléatoire suivant une loi normale. La moyenne de cette loi normale est alors considérée comme la valeur « réelle » de la grandeur observée, et la dispersion de la loi renseigne sur l'« erreur » d'observation. C'est-à-dire qu'il est possible de calculer une valeur approchée de la probabilité qu'une variable suivant une loi normale soit dans un intervalle [μ – σ, μ + σ] autour de la moyenne μ. Il s'agit de pouvoir obtenir une approximation de la grandeur observée dans l'expérience en considérant les erreurs dues aux instruments de mesure ou autres.

Histoire

Le théorème central limite, Pierre-Simon de Laplace, Essai sur la philosophie des probabilités, 1840, 6^e éd., p. 90

Une des premières apparitions d'une loi normale est due à Abraham de Moivre en 1733 en approfondissant l'étude de la factorielle n! lors de l'étude d'un jeu de pile ou face. Il publie The Doctrine of Chances en 1756 dans lequel une loi normale apparaît comme limite d'une loi binomiale, ce qui sera à l'origine du théorème central limite. En 1777, Pierre-Simon de Laplace reprend ces travaux et obtient une bonne approximation de l'erreur entre cette loi normale et une loi binomiale grâce à la fonction gamma d'Euler. Dans son ouvrage publié en 1781^{[réf. souhaitée]}, Laplace donne une première table de cette loi. En 1809, Carl Friedrich Gauss assimile des erreurs d'observation en astronomie à la courbe, dite des erreurs, de la densité d'une loi normale.

Une loi normale est alors pleinement définie lorsque le premier théorème central limite, alors appelé théorème de Laplace, est énoncé par Laplace en 1812. Son nom « normale » est donné par Henri Poincaré à la fin du XIX^e siècle. Les lois normales portent également les noms de lois de Gauss ou lois de Laplace-Gauss en fonction de l'attribution de la paternité de la création de ces lois ; la dénomination de deuxième loi de Laplace est également utilisée occasionnellement.

Les études sur les lois normales se poursuivent durant le XIX^e siècle. Ainsi, de nouvelles tables numériques sont données en 1948 par Egon Sharpe Pearson, en 1952 par le National Bureau of Standards et en 1958 par Greenwood et Hartley.

Loi normale centrée réduite

Une loi normale est une loi de probabilité (c'est-à-dire une mesure N de masse totale unitaire) unidimensionnelle (c'est-à-dire à support réel $\mathbb {R}$). C'est une loi absolument continue, c'est-à-dire que la mesure est absolument continue par rapport à la mesure de Lebesgue. Autrement dit, il existe une densité de probabilité, souvent notée φ pour la loi normale centrée réduite, telle que : N(dx) = φ(x) dx. Elle est généralisée par la loi normale multidimensionnelle. La loi normale centrée réduite est appelée loi normale standard.

Définition par la fonction de densité

Fonction de densité de la loi normale centrée réduite (dite courbe de Gauss ou courbe en cloche).

La loi normale centrée réduite est la loi de probabilité absolument continue dont la densité de probabilité est donnée par la fonction $\varphi :\mathbb {R} \to \mathbb {R} _{+}$ définie par : $\varphi (t)={\frac {1}{\sqrt {2\pi }}}\operatorname {e} ^{-{\frac {1}{2}}t^{2}}$, pour tout $t\in \mathbb {R}$. Cette loi est dite centrée puisque son moment d'ordre 1 (espérance) vaut 0 et réduite puisque son moment d'ordre 2 (variance) vaut 1, tout comme son écart type. Le graphe de la densité φ est appelé fonction gaussienne, courbe de Gauss ou courbe en cloche. Cette loi est notée grâce à la première lettre de « normal », une variable aléatoire X qui suit la loi normale centrée réduite est notée : $X\sim {\mathcal {N}}(0,1)$.

Quelques remarques et propriétés immédiates (voir également les propriétés ci-dessous) :

• le calcul de l'intégrale de Gauss permet de démontrer que la fonction φ est une densité de probabilité par la formule : $\int _{-\infty }^{+\infty }\mathrm {exp} \left(-{\frac {t^{2}}{2}}\right)\,\mathrm {d} t={\sqrt {2\pi }}$ ;
• la densité φ est continue, uniformément bornée et paire ;
• cette parité fait que l'espérance et les moments d'ordres impairs sont nuls ;
• les moments d'ordres pairs sont donnés par $m_{2k}=(2k-1)\cdots 3\cdot 1={\frac {(2k)!}{2^{k}k!}}$(en particulier, Var(X) = m₂ = 1), d'après la relation de récurrence m_2k = (2k–1)m_2k–2 pour k ≥ 1, qui provient de l'intégration par parties suivante :$m_{2k}=\int _{-\infty }^{+\infty }t^{2k-1}t\varphi (t)\,\mathrm {d} t=-\int _{-\infty }^{+\infty }t^{2k-1}\varphi '(t)\,\mathrm {d} t=(2k-1)\int _{-\infty }^{+\infty }t^{2k-2}\varphi (t)\,\mathrm {d} t$.
• le maximum de la fonction φ est atteint en la moyenne 0 et vaut ${\frac {1}{\sqrt {2\pi }}}$ ;
• la fonction vérifie : $\lim _{x\to +\infty }\varphi (x)=\lim _{x\to -\infty }\varphi (x)=0$ ;
• la densité φ est infiniment dérivable ; un raisonnement par récurrence permet d'obtenir la formule : $\varphi ^{(n)}(x)=(-1)^{n}H_{n}(x)\varphi (x)$ où H_n est le n-ième polynôme d'Hermite ;
• la densité possède deux points d'inflexion en 1 et en –1. Ce sont les points en lesquels la dérivée seconde φ'' s'annule et change de signe. Les deux points se situent approximativement aux trois cinquièmes de la hauteur totale.

Définition par la fonction de répartition

Fonction de répartition de la loi normale centrée réduite.

Historiquement, une loi normale est apparue comme la loi limite dans le théorème central limite à l'aide de sa fonction de répartition. Il est alors utile de définir la loi par cette fonction. La loi normale est la loi de probabilité dont la fonction de répartition est donnée par la fonction $\Phi :\mathbb {R} \to \mathbb {R} _{+}$ définie par : $\Phi (x)={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{x}\operatorname {e} ^{-{\frac {1}{2}}t^{2}}\,\mathrm {d} t$, pour tout $x\in \mathbb {R}$. Elle donne la probabilité qu'une variable aléatoire de loi normale appartienne à un intervalle $[a,b]$ : $\mathbb {P} (X\in [a,b])=\Phi (b)-\Phi (a)$ (pour plus de détails de calcul, voir la section Tables numériques et calculs).

Quelques remarques et propriétés immédiates :

• il n'existe pas d'expression analytique de la fonction de répartition Φ, c'est-à-dire qu'elle ne s'exprime pas à partir de fonctions usuelles mais devient elle-même une fonction usuelle ;
• elle s'exprime en fonction de la fonction d'erreur grâce aux deux formules équivalentes suivantes :
  1. $\Phi (x)={\frac {1}{2}}+{\frac {1}{2}}\operatorname {erf} \left({\frac {x}{\sqrt {2}}}\right)$,
  2. $\operatorname {erf} (x)=2\Phi \left(x{\sqrt {2}}\right)-1$ ;
• elle est dérivable une infinité de fois et vérifie Φ'(x) = φ(x). L'écriture équivalente dΦ(x) = φ(x)dx permet de définir l'intégrale de Lebesgue-Stieltjes par rapport à la loi normale ;
• elle est absolument continue et strictement croissante, c'est donc une bijection de $\mathbb {R}$ dans ]0 ; 1[. Sa réciproque Φ⁻¹ existe et s'appelle la fonction probit. Cette fonction est utilisée pour le modèle probit ;
• par parité de la loi, Φ(–x) = 1 – Φ(x) et ainsi Φ(0) = 1/2. Ceci montre que la médiane de la loi normale centrée réduite est 0 ;
• par définition de la fonction de répartition, $\Phi (x)=\mathbb {P} (X\leq x)$ lorsque la variable aléatoire X suit la loi normale centrée réduite, $X\sim {\mathcal {N}}(0,1)$. Pour obtenir les valeurs de cette probabilité, il faut approcher cette fonction par d'autres fonctions usuelles et il existe des tables de valeurs (voir la section Table de la loi normale ci-dessous).

Définition par la fonction caractéristique

Fonction caractéristique et fonction génératrice des moments de la loi normale centrée réduite.

La caractérisation d'une loi normale par sa fonction caractéristique présente un intérêt pour démontrer certaines propriétés, comme la stabilité par addition ou le théorème central limite. Cette fonction caractéristique $\phi :\mathbb {R} \to \mathbb {R} _{+}$, qui se calcule à partir de la densité de probabilité et caractérise la loi, est donnée par : $\phi (t)={\rm {e}}^{-{\frac {t^{2}}{2}}},\,\forall t\in \mathbb {R}$. Cette fonction caractéristique est égale, à une constante multiplicative près, à la densité de probabilité de la loi : on dit que la fonction caractéristique d'une gaussienne est gaussienne.

Si une variable aléatoire X suit la loi normale centrée réduite de fonction caractéristique ϕ définie ci-dessus, alors la transformation linéaire Y = aX+b admet pour fonction caractéristique : $\phi _{Y}(t)={\rm {e}}^{{\rm {i}}bt}\phi (at)$. C'est donc une variable aléatoire de loi normale de moyenne b et de variance a².

Définition par la fonction génératrice des moments

Une autre manière de définir une loi normale est par l'utilisation de sa fonction génératrice des moments $M:\mathbb {R} \to \mathbb {R} _{+},\;t\mapsto \sum _{n=0}^{\infty }m_{n}{\frac {t^{n}}{n!}}=\mathbb {E} \left(\operatorname {e} ^{tX}\right)$. Cette fonction, qui se calcule à partir de la fonction de densité et caractérise la loi, est donnée par : $M(t)={\rm {e}}^{\frac {t^{2}}{2}}$, pour tout $t\in \mathbb {R}$. On retrouve ainsi les valeurs des moments m_n (voir supra).

Loi normale générale

Définition

Plus généralement que la loi normale centrée réduite, une loi normale (non centrée et non réduite) est une loi de probabilité absolument continue dont l'un des quatre points suivants est vérifié :

• la densité de probabilité $\varphi :\mathbb {R} \to \mathbb {R} _{+}$ est donnée par :$\varphi (t)={\frac {1}{\sigma {\sqrt {2\pi }}}}\operatorname {e} ^{-{\frac {1}{2}}{\frac {(t-\mu )^{2}}{\sigma ^{2}}}}$, pour tout $t\in \mathbb {R}$ ;
• la fonction de répartition $F:\mathbb {R} \to \mathbb {R} _{+}$ est donnée par :$F(x)={\frac {1}{\sigma {\sqrt {2\pi }}}}\int _{-\infty }^{x}\mathrm {e} ^{-{\frac {1}{2}}{\frac {(t-\mu )^{2}}{\sigma ^{2}}}}\,\mathrm {d} t$, pour tout $x\in \mathbb {R}$ ;
• la fonction caractéristique $\phi :\mathbb {R} \to \mathbb {C}$ est donnée par :$\phi (t)={\rm {e}}^{\mu {\rm {i}}t-{\frac {1}{2}}\sigma ^{2}t^{2}}$, pour tout $t\in \mathbb {R}$ ;
• la fonction génératrice des moments $M:\mathbb {R} \to \mathbb {R} _{+}$ est donnée par :$M(t)={\rm {e}}^{\mu t+{\frac {1}{2}}\sigma ^{2}t^{2}}$, pour tout $t\in \mathbb {R}$,où $\mu \in \mathbb {R}$ et $\sigma \in \mathbb {R} _{+}^{\star }$.

Pour le cas où σ = 0, c'est une forme dégénérée de la loi normale, parfois appelée loi normale impropre. C'est alors la mesure de Dirac au point μ qui n'est pas absolument continue.

La valeur μ est la moyenne de la loi et σ est l'écart type alors que σ² en est la variance. Cette loi est notée grâce à la première lettre de « normal ». Une variable aléatoire X qui suit une loi normale est notée de deux manières différentes suivant les auteurs :$X\sim {\mathcal {N}}(\mu ,\sigma ){\text{ ou }}X\sim {\mathcal {N}}(\mu ,\sigma ^{2}).$ La deuxième notation a l'intérêt de pouvoir noter la stabilité par addition de manière simple ; elle sera utilisée dans cet article.

Remarques et propriétés immédiates

• Si la variable aléatoire X suit la loi normale centrée réduite ${\mathcal {N}}(0,1)$, alors la variable aléatoire σX + μ suit la loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$ de moyenne μ et de variance σ². Réciproquement, si Y suit la loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$, alors Y – μ/σ suit la loi normale centrée réduite. Dit autrement, toute loi normale peut s'obtenir par translation (shifting en anglais) et par dilatation (scaling en anglais) de la loi centrée réduite.Cette première propriété permet d'obtenir la formule très utile :$\mathbb {P} (Y\leq x)=\mathbb {P} \left({\frac {Y-\mu }{\sigma }}\leq {\frac {x-\mu }{\sigma }}\right)=\mathbb {P} \left(X\leq {\frac {x-\mu }{\sigma }}\right)$.Il est alors possible de déduire les propriétés d'une loi normale à partir de celles de la loi normale centrée réduite, et vice versa. La variable Y – μ/σ est parfois appelée la « standardisation » de Y ou « variable Y centrée réduite ».
• Plus généralement, si $X\sim {\mathcal {N}}(\mu ,\sigma ^{2})$ alors $aX+b\sim {\mathcal {N}}(a\mu +b,a^{2}\sigma ^{2})$ (on peut le lire sur les fonctions génératrices des moments, ou sur les densités).
• La densité f est symétrique par rapport à μ.
• Le maximum de la fonction f est atteint en μ et vaut ${\frac {1}{\sigma {\sqrt {2\pi }}}}$.
• La décroissance de la densité à droite et à gauche de μ est surexponentielle.
• Puisqu'une loi normale est une loi de probabilité absolument continue, l'événement [X = x] est négligeable, c'est-à-dire que presque sûrement une variable aléatoire de loi normale X n'est jamais égale à une valeur fixée x. Ceci se traduit mathématiquement par : $\mathbb {P} (X=x)=0$.
• La largeur à mi-hauteur permet de donner une valeur d'amplitude de la loi. C'est la largeur de la courbe à une hauteur qui vaut la moitié de la hauteur totale. Cette largeur à mi-hauteur de la loi normale est proportionnelle à l'écart type : H = 2 √2 ln (2) σ ≈ 2,3548 σ. Le facteur 2 est issu de la propriété de symétrie de la loi normale.
• La densité possède deux points d'inflexion en μ + σ et en μ – σ. Ce sont les points en lesquels la dérivée seconde f'' s'annule et change de signe. Les deux points se situent approximativement aux trois cinquièmes de la hauteur totale.
• Les lois normales sont des lois de la famille exponentielle, c'est-à-dire que leur densité s'écrit sous la forme :$f(x)=a(\theta )b(x)\mathrm {e} ^{-c(\theta )d(x)}$ ou, de manière équivalente, sous la forme$f(x)=\exp \left({\frac {x\theta _{1}-\beta (\theta _{1})}{\alpha (\theta _{2})}}\right),\ {\textrm {avec}}\ \theta _{1}=\mu ,\,\theta _{2}=\sigma ,\,\beta (\mu )=\mu ^{2}/2\ {\textrm {et}}\ \alpha (\sigma )=\sigma ^{2}$.

Propriétés

Autres caractérisations

En addition de la densité de probabilité, de la fonction de répartition, de la fonction caractéristique et de la fonction génératrice des moments, il existe d'autres caractérisations des lois normales.

• Caractérisation due à Georges Darmois (1951) et Sergueï Bernstein (1954) : si deux variables aléatoires X₁ et X₂ sont indépendantes et de même loi et si les deux variables aléatoires X₁ + X₂ et X₁ – X₂ sont également indépendantes, alors la loi commune X₁ et X₂ est une loi normale.
• Caractérisation due à Charles Stein (1972) : les lois normales sont les seules lois de probabilité (mesures de probabilité) $\mathbb {P}$ telles que, pour toute fonction g de classe C¹ (c'est-à-dire dérivable et de dérivée continue) :$\int _{\mathbb {R} }g'(x)\mathrm {d} \mathbb {P} (x)=\int _{\mathbb {R} }xg(x)\mathrm {d} \mathbb {P} (x)$.

Moments

Le moment d'ordre 1 est appelé la moyenne (μ) et est donné en paramètre d'une loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$. Le deuxième paramètre est son écart type (σ), c'est-à-dire la racine carrée de la variance qui est par définition la moyenne des carrés des écarts à la moyenne. Il est alors également intéressant d'obtenir les moments centrés d'une loi normale, ils sont donnés par :

${\begin{cases}\mu _{2k}=\mathbb {E} [(X-\mu )^{2k}]={\frac {(2\,k)!}{2^{k}k!}}\sigma ^{2k}\\\mu _{2k+1}=\mathbb {E} [(X-\mu )^{2k+1}]=0\end{cases}}$ pour $k\geq 0$ et X une variable aléatoire de loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$.

Le moment ordinaire m_n peut s'obtenir à partir des moments d'ordre inférieur à n – 1 et du moment centré d'ordre n, en utilisant la formule qui exprime μ_n en fonction de m₀, m₁, ..., m_n. Les premiers moments d'une loi normale sont alors : $m_{1}=\mathbb {E} [X]=\mu ,\quad m_{2}=\mathbb {E} [X^{2}]=\sigma ^{2}+\mu ^{2},\quad m_{3}=\mathbb {E} [X^{3}]=3\mu \sigma ^{2}+\mu ^{3},\quad m_{4}=\mathbb {E} [X^{4}]=3\sigma ^{4}+6\sigma ^{2}\mu ^{2}+\mu ^{4}$.

Calcul direct

Grâce à la symétrie autour de μ de la fonction de densité d'une loi normale, les moments centrés d'ordre impair sont tous nuls.

Des moments d'ordre pairs de la loi normale centrée réduite (voir supra), on déduit la formule des moments centrés : $\mu _{2k}={\frac {(2\,k)!}{2^{k}k!}}\sigma ^{2k}$.

Par la fonction génératrice des moments

Les moments centrés (μ_n, n ≥ 0) d'une loi peuvent s'obtenir à partir de la fonction génératrice des moments centrés. Le cas particulier μ = 0 de la fonction génératrice des moments (voir supra) donne : $M_{\text{centré}}(t)={\rm {e}}^{\frac {\sigma ^{2}t^{2}}{2}}=\sum _{k=0}^{\infty }{\frac {1}{k!}}\left({\frac {\sigma ^{2}t^{2}}{2}}\right)^{k}$.

Comme par ailleurs on a (pour toute loi) $M_{\text{centré}}(t)=\sum _{n=0}^{\infty }{\frac {1}{n!}}\mu _{n}t^{n}$, on en déduit, par identification des coefficients des deux séries, les moments centrés d'une loi normale (voir supra).

Quant aux moments ordinaires, leur fonction génératrice permet d'établir la relation de récurrence :

$m_{n+1}=\mu m_{n}+n\sigma ^{2}m_{n-1}\quad (n\geq 1)$.

Densités de probabilités de lois de kurtosis différents. Loi de Laplace en rouge, loi sécante hyperbolique en orange, loi logistique en vert, loi normale en noir, loi du cosinus surélevé en cyan, loi du demi-cercle en bleu et loi uniforme en violet.

Asymétrie et aplatissement

L'asymétrie γ₁, le kurtosis β₂ et le kurtosis normalisé γ₂ s'obtiennent à partir des formules des moments :

$\gamma _{1}={\frac {\mu _{3}}{\sigma ^{3}}}=0$ ;
$\beta _{2}={\frac {\mu _{4}}{\sigma ^{4}}}=3$ ;
$\gamma _{2}=\beta _{2}-3=0$.

Les lois normales servent de point de référence pour la comparaison des épaisseurs de traîne : si une loi possède un kurtosis normalisé γ₂ > 0, alors la loi possède une traîne plus épaisse qu'une loi normale et est dite leptokurtique ; à l'inverse si γ₂ < 0, la loi possède une traîne moins épaisse qu'une loi normale et est appelée platikurtique ; les lois de kurtosis normalisé nul possèdent une traîne comparable à la loi normale et sont dites mésokurtiques.

Cumulants

La fonction caractéristique permet d'obtenir la fonction génératrice des cumulants par la formule $\ln(\phi (t))=\sum _{n=1}^{+\infty }K_{n}{\frac {({\rm {i}}t)^{n}}{n!}}$ et permet d'obtenir les cumulants : K₁ = μ, K₂ = σ² et K_n = 0 pour n ≥ 3.

Théorèmes de convergence

Lorsque le nombre n de variables augmente, la densité de probabilité de la variable S_n (centrée réduite) se rapproche de la courbe en cloche d'une loi normale.

La première version du théorème central limite, appelé alors théorème de Moivre-Laplace, a été énoncée dans le cas de variables aléatoires de loi de Bernoulli. De manière plus générale, si X₁, X₂, ..., X_n sont des variables indépendantes et identiquement distribuées de variance finie et si la somme est notée S_n = X₁ + X₂ + ... + X_n, alors pour tout a < b $\lim _{n\to +\infty }\mathbb {P} \left(a\leq {\frac {S_{n}-\mathbb {E} [S_{n}]}{\sqrt {\operatorname {Var} (S_{n})}}}\leq b\right)=\int _{a}^{b}\varphi (x)\,\mathrm {d} x$ où φ est la densité de probabilité de la loi normale centrée réduite.

Ce théorème signifie que tout ce qui peut être considéré comme étant la somme d'une grande quantité de petites valeurs aléatoires indépendantes et identiquement distribuées est approximativement de loi normale. Ceci montre le caractère central des lois normales en théorie des probabilités. Un énoncé physique de ce théorème peut être formulé : Si une grandeur physique subit l'influence additive d'un nombre important de facteurs indépendants et si l'influence de chaque facteur pris séparément est petite, alors la distribution de cette grandeur est une distribution gaussienne.

Ce théorème central limite est valide pour toute loi de probabilité initiale des variables iid X₁, X₂, ..., X_n ayant un écart type fini, il permet d'obtenir de bonnes approximations de la somme S_n, par exemple :

• si les variables X_i sont de loi de Bernoulli ${\mathcal {B}}(p)$, alors S_n suit approximativement une loi normale ${\mathcal {N}}(np,np(1-p))$. Cette approximation est satisfaisante dans le cas où np(1–p) > 10 ;
• si les variables X_i sont de loi du χ² : χ²(1), alors S_n suit approximativement une loi normale ${\mathcal {N}}(n,4n^{2})$ ;
• si les variables X_i sont de loi exponentielle : ${\mathcal {E}}(\lambda )$, alors S_n suit approximativement une loi normale ${\mathcal {N}}\left({\frac {n}{\lambda }},{\frac {n}{\lambda ^{2}}}\right)$.

Il existe des versions plus générales de ce théorème, par exemple en considérant des variables aléatoires indépendantes, pas de même loi mais ayant des variances petites comparées à celle de leur moyenne. Un théorème de Gnedenko et Kolmogorov (1954) stipule qu'une variable aléatoire normale est la somme d'un grand nombre de variables aléatoires indépendantes petites dont aucune n'est prépondérante :

Théorème — Considérons une suite de variables aléatoires (X_n, n ≥ 1) dont chacune est la somme d'un nombre fini de variables aléatoires $X_{n,1},\dots ,X_{n,k_{n}}$ avec $k_{n}\to +\infty$.

Pour tout ε > 0, introduisons la variable aléatoire tronquée :

$X^{\varepsilon }={\begin{cases}X&{\text{ si }}|X|\leq \varepsilon {\text{ ;}}\\0&{\text{ sinon}}\end{cases}}$

et supposons :

• $\sum _{1\leq k\leq n}|X_{nk}|{\underset {n\to \infty }{\longrightarrow }}0$ (en probabilité) ;
• pour tout ε > 0, $\sum _{1\leq k\leq n}\mathbb {E} [X_{nk}^{\varepsilon }]{\underset {n\to \infty }{\longrightarrow }}\mu$ et $\sum _{1\leq k\leq n}{\text{Var}}[X_{nk}^{\varepsilon }]{\underset {n\to \infty }{\longrightarrow }}\sigma ^{2}$.

Alors la loi de X_n converge vers la loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$.

Stabilités et famille normale

Stabilité par additivité (propriété de conservation)

Les lois normales sont stables par additivité, c'est-à-dire que la somme de deux variables aléatoires indépendantes de lois normales est elle-même une variable aléatoire de loi normale. Plus explicitement : si $X_{1}\sim {\mathcal {N}}(\mu _{1},\sigma _{1}^{2})$, $X_{2}\sim {\mathcal {N}}(\mu _{2},\sigma _{2}^{2})$ et X₁ et X₂ sont indépendantes, alors la variable aléatoire $X_{1}+X_{2}$ suit la loi normale ${\mathcal {N}}(\mu _{1}+\mu _{2},\sigma _{1}^{2}+\sigma _{2}^{2})$.

Cette propriété se généralise pour n variables, c'est-à-dire que si pour tout $i\in \{1,2,\dots ,n\}$, les variables aléatoires X_i suivent une loi normale ${\mathcal {N}}(\mu _{i},\sigma _{i}^{2})$ et sont indépendantes, alors la somme X₁ + X₂ + ... + X_n suit une loi normale ${\mathcal {N}}(\mu _{1}+\mu _{2}+\dots +\mu _{n},\sigma _{1}^{2}+\sigma _{2}^{2}+\dots +\sigma _{n}^{2})$.

Cette propriété se démontre directement au moyen des fonctions caractéristiques. La densité de probabilité de la somme de deux variables indépendantes de loi normale est donnée par la convolution des deux densités. Ceci se traduit par les formules de convolution de fonctions ou de convolution de mesures normales que l'on note ${\mathcal {N}}_{\mu _{1},\sigma _{1}^{2}}$ : $\varphi \left({\frac {x-\mu _{1}}{\sigma _{1}}}\right)\ast \varphi \left({\frac {x-\mu _{2}}{\sigma _{2}}}\right)=\varphi \left({\frac {x-(\mu _{1}+\mu _{2})}{\sqrt {\sigma _{1}^{2}+\sigma _{2}^{2}}}}\right)$ et ${\mathcal {N}}_{\mu _{1},\sigma _{1}^{2}}\ast {\mathcal {N}}_{\mu _{2},\sigma _{2}^{2}}={\mathcal {N}}_{\mu _{1}+\mu _{2},\sigma _{1}^{2}+\sigma _{2}^{2}}$.

Il ne faut pas confondre avec la loi dont la densité est la somme de densités de lois normales (voir la section Constructions à partir de la loi normale ci-dessous).

Différentes densités de probabilité de lois stables dont la loi normale est un cas particulier : la courbe noire est la courbe en cloche.

Famille normale

L'ensemble de fonctions $\{\varphi ({\frac {x-\mu }{\sigma }})\mid \mu \in \mathbb {R} ,\sigma >0\}$ forme la famille dite famille normale. La famille normale est également le nom de l'ensemble des lois normales $\{{\mathcal {N}}_{\mu ,\sigma ^{2}}\mid \mu \in \mathbb {R} ,\sigma >0\}$. La famille de fonctions est fermée pour la convolution au sens où : la fonction φ est la seule qui engendre la famille ; si la convolution de deux densités est dans la famille alors les deux fonctions sont dans la famille ; et toute densité convolée un nombre suffisamment grand de fois et convenablement renormalisée est proche d'une fonction de la famille normale. Les trois théorèmes suivants donnent plus de précisions mathématiques.

1. Théorème : si pour une fonction de densité f de moyenne 0 et d'écart type 1, il existe $\mu \in \mathbb {R}$ et $\sigma \in \mathbb {R} _{+}^{*}$ satisfaisant :$f\left({\frac {x-\mu _{1}}{\sigma _{1}}}\right)\ast f\left({\frac {x-\mu _{2}}{\sigma _{2}}}\right)=f\left({\frac {x-\mu }{\sigma }}\right)$,alors $f\equiv \varphi$ est la densité de la loi normale centrée réduite.
2. Théorème de Lévy-Cramér (1936) (conjecturé par Paul Lévy en 1935) : si deux fonctions de densité, f₁ et f₂, vérifient :$f_{1}(x)\ast f_{2}(x)=\varphi \left({\frac {x-\mu }{\sigma }}\right)$,alors$f_{1}(x)=\varphi \left({\frac {x-\mu _{1}}{\sigma _{1}}}\right)$ et $f_{2}(x)=\varphi \left({\frac {x-\mu _{2}}{\sigma _{2}}}\right)$ avec μ₁ + μ₂ = μ et σ₁ + σ₂ = σ. Autrement dit, si la somme de deux variables aléatoires indépendantes est normale, alors les deux variables sont de lois normales.
3. Théorème : si f est la densité commune de n variables aléatoires indépendantes de moyenne 0 et d'écart type 1, alors la convolée n fois de f converge uniformément en x : $\left(f(x/{\sqrt {n}})\right)^{\ast n}\to \varphi (x)$ (ce théorème est équivalent au théorème central limite). Il ne faut pas confondre cette famille normale avec la famille normale de fonctions holomorphes.

Stabilité par linéarité

Les lois normales sont stables par linéarité : si α ≥ 0 et β sont deux réels et $X\sim {\mathcal {N}}(\mu ,\sigma ^{2})$, alors la variable aléatoire αX + β suit la loi normale ${\mathcal {N}}(\alpha \mu +\beta ,\alpha ^{2}\sigma ^{2})$.

Grâce aux stabilités par addition et par linéarité, une loi normale est un cas particulier de loi stable avec pour paramètre de stabilité α = 2. Parmi les lois stables, les lois normales, la loi de Lévy (α = 1/2) et la loi de Cauchy (α = 1) sont les seules à posséder une expression analytique de leur fonction de densité.

Stabilité par moyenne

Les lois normales sont stables par moyennisation, c'est-à-dire si X₁, X₂, ..., X_n sont des variables aléatoires indépendantes suivant respectivement les lois normales ${\mathcal {N}}(\mu _{1},\sigma _{1}^{2}),{\mathcal {N}}(\mu _{2},\sigma _{2}^{2}),\dots ,{\mathcal {N}}(\mu _{n},\sigma _{n}^{2})$, alors la moyenne 1/n(X₁ + X₂ + ... + X_n) suit la loi ${\mathcal {N}}\left({\tfrac {\mu _{1}+\mu _{2}+\dots +\mu _{n}}{n}},{\tfrac {\sigma _{1}^{2}+\sigma _{2}^{2}+\dots +\sigma _{n}^{2}}{n^{2}}}\right)$.

Convexité

Les lois normales ne sont pas convexe, c'est-à-dire que l'inégalité $\lambda \mathbb {P} (A)+(1-\lambda )\mathbb {P} (B)\leq \mathbb {P} (\lambda A+(1-\lambda )B)$ pour tous boréliens A et B n'est pas vérifiée lorsque la mesure $\mathbb {P}$ est normale. Cependant, lorsque l'on normalise cette inégalité avec l'inverse de la fonction de répartition de la loi normale centrée réduite, on obtient le théorème suivant, analogue à l'inégalité de Brunn-Minkowski-Lusternik pour la mesure de Lebesgue dans $\mathbb {R} ^{n}$ :

Inégalité de Ehrhard — Pour la mesure normale standard ${\mathcal {N}}_{0,1}$, pour tous ensembles boréliens A et B et pour tout λ ∈ ]0 ; 1[,

$\lambda \Phi ^{-1}\left({\mathcal {N}}_{0,1}(A)\right)+(1-\lambda )\Phi ^{-1}\left({\mathcal {N}}_{0,1}(B)\right)\leq \Phi ^{-1}\left({\mathcal {N}}_{0,1}(\lambda A+(1-\lambda )B)\right)$.

En fait, les lois normales font partie de la famille des distributions de mesures log-concaves, c'est-à-dire vérifiant pour tous boréliens A et B et tout $\lambda \in ]0,1[$,

$\mathbb {P} (\lambda A+(1-\lambda )B)\geq \mathbb {P} (A)^{\lambda }\mathbb {P} (B)^{1-\lambda }.$

Entropie et quantité d'information

Entropie de Shannon

L'entropie de Shannon d'une loi de probabilité absolument continue de densité donnée par f permet de mesurer une quantité d'information et est définie par : $H=-\int _{-\infty }^{+\infty }f(x)\ln f(x)\,\mathrm {d} x$. Dans l'ensemble des lois absolument continues de variance σ² fixée, les lois normales ${\mathcal {N}}(\cdot ,\sigma ^{2})$ sont d'entropie maximum. L'entropie maximum, pour une loi normale donc, est donnée par : H = ln (σ √2πe). Ainsi la théorie de maximisation de l'entropie dit que, même si elle n'est pas la meilleure loi adaptée aux valeurs, une loi normale ajustée aux valeurs est adéquate pour prendre une décision.

Il y a également une connexion entre la convergence de suites de lois de probabilité vers une loi normale et la croissance de l'entropie, ce qui en fait un outil majeur dans la théorie de l'information.

La quantité d'information de Fisher

L'information de Fisher d'une loi à densité de probabilité est une autre notion de quantité d'information. Pour une densité f, elle est donnée par : $I=\int _{-\infty }^{+\infty }\left({\frac {f'(x)}{f(x)}}\right)^{2}f(x)\,\mathrm {d} x$. Pour toute densité suffisamment régulière d'une loi centrée réduite, cette information vérifie I ≥ 1. Les lois normales se distinguent des autres densités puisque l'inégalité précédente est une égalité si et seulement si la densité est celle de la loi normale centrée réduite.

Distance entre lois

La divergence de Kullback-Leibler entre deux lois permet de mesurer une distance entre les deux lois, ou une perte d'information entre les deux lois. La divergence de Kullback-Leibler entre les deux lois normales ${\mathcal {N}}(\mu _{1},\sigma _{1}^{2})$ et ${\mathcal {N}}(\mu _{2},\sigma _{2}^{2})$ est : $D_{KL}({\mathcal {N}}(\mu _{1},\sigma _{1}^{2})\|{\mathcal {N}}(\mu _{2},\sigma _{2}^{2}))=\log \left({\frac {\sigma _{2}}{\sigma _{1}}}\right)+{\frac {1}{2}}\left({\frac {\sigma _{1}^{2}}{\sigma _{2}^{2}}}+{\frac {(\mu _{2}-\mu _{1})^{2}}{\sigma _{2}^{2}}}-1\right)$. Cette divergence est nulle pour μ₁ = μ₂ et σ₁ = σ₂ ; de plus, elle croît lorsque $|\mu _{1}-\mu _{2}|$ croît.

Approximation de la fonction de répartition

Il n'existe pas d'expression analytique pour la fonction de répartition Φ de la loi normale centrée réduite, c'est-à-dire qu'il n'existe pas de formule simple entre la fonction de répartition et les fonctions classiques telles que les fonctions polynomiales, exponentielle, logarithmique, trigonométriques, etc. Cependant la fonction de répartition apparaît dans plusieurs résultats à vocation à être appliqués, il est donc important de mieux cerner cette fonction. Différentes écritures sous forme de séries ou de fractions continues généralisées sont possibles.

Pour les valeurs de $0<x\ll 1$, la fonction de répartition de la loi normale centrée réduite s'écrit sous la forme : $\Phi (x)={\frac {1}{2}}+{\frac {1}{\sqrt {2\pi }}}\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{n!2^{n}(2n+1)}}x^{2n+1}={\frac {1}{2}}+{\frac {1}{\sqrt {2\pi }}}\left(x-{\frac {x^{3}}{6}}+{\frac {x^{5}}{40}}+\dots \right)$, ou sous la forme : $\Phi (x)={\frac {1}{2}}+\varphi (x)\sum _{n=0}^{\infty }{\frac {1}{1\cdot 3\cdot 5\dots (2n+1)}}x^{2n+1}={\frac {1}{2}}+\varphi (x)\left(x+{\frac {x^{3}}{3}}+{\frac {x^{5}}{15}}+\dots \right)$.

Pour $1\ll x$, la fonction de répartition de la loi normale centrée réduite s'écrit sous la forme : $\Phi (x)=1-{\frac {\varphi (x)}{x}}\left(1-{\frac {1}{x^{2}}}+{\frac {1\cdot 3}{x^{4}}}-{\frac {1\cdot 3\cdot 5}{x^{6}}}+\dots +{\frac {1\cdot 3\dots (2n-1)}{x^{2n}}}\right)+R_{n}$ avec $R_{n}=(-1)^{n+1}1\cdot 3\dots (2n+1)\int _{x}^{\infty }{\frac {\varphi (y)}{y^{2n+2}}}\,\mathrm {d} y$.

De manière plus numérique et facilement calculable, les approximations suivantes donnent des valeurs de la fonction de répartition Φ de la loi normale centrée réduite avec :

• une erreur de l'ordre de 10⁻⁵ : pour x > 0, $\Phi (x)=1-{\frac {{\rm {e}}^{-{\frac {x^{2}}{2}}}}{\sqrt {2\pi }}}\left({\frac {0{,}4361836}{1+0{,}33267\,x}}+{\frac {-0{,}1201676}{(1+0{,}33267\,x)^{2}}}+{\frac {0{,}9372980}{(1+0{,}33267\,x)^{3}}}\right)+\epsilon (x)$ où $|\epsilon (x)|<10^{-5}$ ;
• une erreur de l'ordre de $2{,}5\,.\,10^{-4}$ : pour $x>0$ : $\Phi (x)\approx 1-{\frac {1}{2\left(1+0{,}196854\,x+0{,}115194\,x^{2}+0{,}000344\,x^{3}+0{,}019527\,x^{4}\right)^{4}}}$ ;
• une erreur de l'ordre de $10^{-2}$ : $\Phi (x)={\begin{cases}0{,}1x(4{,}4-x)+0,5&{\text{ pour }}0\leq x\leq 2{,}2\\0{,}99&{\text{ pour }}2{,}2\leq x\leq 2{,}6\\1&{\text{ pour }}x\geq 2{,}6.\end{cases}}$

Voici un exemple d'algorithme pour le langage C :

double Phi(double x){
   long double s=x,t=0,b=x,q=x*x,i=1;

   while(s!=t)
      s = (t=s) + (b*=q/(i+=2));
   return 0.5 + s*exp(-0.5*q - 0.91893853320467274178L);
}

Une autre écriture de la fonction de répartition de la loi normale centrée réduite utilise une fraction continue généralisée : $\Phi (x{\sqrt {2}})={\frac {1}{2}}-{\cfrac {1}{\sqrt {\pi }}}{\cfrac {{\cfrac {1}{2}}{\rm {e}}^{-x^{2}}}{x+{\cfrac {1}{2x+{\cfrac {2}{x+{\cfrac {3}{2x+{\cfrac {4}{x+\dots }}}}}}}}}}$.

Tables numériques et calculs

Fonction quantile d'une loi normale centrée réduite.

Comme mentionné dans la section précédente, il est utile de bien connaître la fonction de répartition Φ pour les applications numériques. Des tables de valeurs ont alors été calculées pour la fonction de répartition, mais également pour son inverse, ce qui permet d'obtenir les quantiles et les intervalles de confiance pour un seuil de tolérance fixé.^{[réf. nécessaire]}

Table de valeurs de la fonction de répartition

La table suivante donne les valeurs de la fonction de répartition $\Phi (x)=\mathbb {P} [X\leq x]$, lorsque X suit la loi normale centrée réduite ${\mathcal {N}}(0,1)$.

Les valeurs en début de lignes donnent la première partie de la variable, les valeurs en début de colonnes donnent la deuxième partie. Ainsi la case de la deuxième ligne et troisième colonne donne : Φ(0,12) = 0,54776.

La courbe en cloche est la fonction de densité. La droite verticale est la valeur x. La surface de la partie colorée sous la courbe est la valeur de $\mathbb {P} [X\leq x]=\Phi (x)$.

La courbe en cloche est la fonction de densité. Les droites verticales sont les valeurs x₁ et x₂. La surface de la partie colorée sous la courbe est la valeur de $\mathbb {P} [x_{1}\leq X\leq x_{2}]=\Phi (x_{2})-\Phi (x_{1})$.

La courbe en cloche est la fonction de densité. La droite verticale est la valeur x. La surface de la partie colorée sous la courbe est la valeur de $\mathbb {P} [X\geq x]=1-\Phi (x)$.

Table de valeur de la fonction de répartition de la loi normale

$\Phi (x)$	0,00	0,01	0,02	0,03	0,04	0,05	0,06	0,07	0,08	0,09
0,0	0,50000	0,50399	0,50798	0,51197	0,51595	0,51994	0,52392	0,52790	0,53188	0,53586
0,1	0,53983	0,54380	0,54776	0,55172	0,55567	0,55962	0,56356	0,56749	0,57142	0,57535
0,2	0,57926	0,58317	0,58706	0,59095	0,59483	0,59871	0,60257	0,60642	0,61026	0,61409
0,3	0,61791	0,62172	0,62552	0,62930	0,63307	0,63683	0,64058	0,64431	0,64803	0,65173
0,4	0,65542	0,65910	0,66276	0,66640	0,67003	0,67364	0,67724	0,68082	0,68439	0,68793
0,5	0,69146	0,69497	0,69847	0,70194	0,70540	0,70884	0,71226	0,71566	0,71904	0,72240
0,6	0,72575	0,72907	0,73237	0,73565	0,73891	0,74215	0,74537	0,74857	0,75175	0,75490
0,7	0,75804	0,76115	0,76424	0,76730	0,77035	0,77337	0,77637	0,77935	0,78230	0,78524
0,8	0,78814	0,79103	0,79389	0,79673	0,79955	0,80234	0,80511	0,80785	0,81057	0,81327
0,9	0,81594	0,81859	0,82121	0,82381	0,82639	0,82894	0,83147	0,83398	0,83646	0,83891
1,0	0,84134	0,84375	0,84614	0,84849	0,85083	0,85314	0,85543	0,85769	0,85993	0,86214
1,1	0,86433	0,86650	0,86864	0,87076	0,87286	0,87493	0,87698	0,87900	0,88100	0,88298
1,2	0,88493	0,88686	0,88877	0,89065	0,89251	0,89435	0,89617	0,89796	0,89973	0,90147
1,3	0,90320	0,90490	0,90658	0,90824	0,90988	0,91149	0,91309	0,91466	0,91621	0,91774
1,4	0,91924	0,92073	0,92220	0,92364	0,92507	0,92647	0,92785	0,92922	0,93056	0,93189
1,5	0,93319	0,93448	0,93574	0,93699	0,93822	0,93943	0,94062	0,94179	0,94295	0,94408
1,6	0,94520	0,94630	0,94738	0,94845	0,94950	0,95053	0,95154	0,95254	0,95352	0,95449
1,7	0,95543	0,95637	0,95728	0,95818	0,95907	0,95994	0,96080	0,96164	0,96246	0,96327
1,8	0,96407	0,96485	0,96562	0,96638	0,96712	0,96784	0,96856	0,96926	0,96995	0,97062
1,9	0,97128	0,97193	0,97257	0,97320	0,97381	0,97441	0,97500	0,97558	0,97615	0,97670
2,0	0,97725	0,97778	0,97831	0,97882	0,97932	0,97982	0,98030	0,98077	0,98124	0,98169
2,1	0,98214	0,98257	0,98300	0,98341	0,98382	0,98422	0,98461	0,98500	0,98537	0,98574
2,2	0,98610	0,98645	0,98679	0,98713	0,98745	0,98778	0,98809	0,98840	0,98870	0,98899
2,3	0,98928	0,98956	0,98983	0,99010	0,99036	0,99061	0,99086	0,99111	0,99134	0,99158
2,4	0,99180	0,99202	0,99224	0,99245	0,99266	0,99286	0,99305	0,99324	0,99343	0,99361
2,5	0,99379	0,99396	0,99413	0,99430	0,99446	0,99461	0,99477	0,99492	0,99506	0,99520
2,6	0,99534	0,99547	0,99560	0,99573	0,99585	0,99598	0,99609	0,99621	0,99632	0,99643
2,7	0,99653	0,99664	0,99674	0,99683	0,99693	0,99702	0,99711	0,99720	0,99728	0,99736
2,8	0,99744	0,99752	0,99760	0,99767	0,99774	0,99781	0,99788	0,99795	0,99801	0,99807
2,9	0,99813	0,99819	0,99825	0,99831	0,99836	0,99841	0,99846	0,99851	0,99856	0,99861
3,0	0,99865	0,99869	0,99874	0,99878	0,99882	0,99886	0,99889	0,99893	0,99896	0,99900
3,1	0,99903	0,99906	0,99910	0,99913	0,99916	0,99918	0,99921	0,99924	0,99926	0,99929
3,2	0,99931	0,99934	0,99936	0,99938	0,99940	0,99942	0,99944	0,99946	0,99948	0,99950
3,3	0,99952	0,99953	0,99955	0,99957	0,99958	0,99960	0,99961	0,99962	0,99964	0,99965
3,4	0,99966	0,99968	0,99969	0,99970	0,99971	0,99972	0,99973	0,99974	0,99975	0,99976
3,5	0,99977	0,99978	0,99978	0,99979	0,99980	0,99981	0,99981	0,99982	0,99983	0,99983
3,6	0,99984	0,99985	0,99985	0,99986	0,99986	0,99987	0,99987	0,99988	0,99988	0,99989
3,7	0,99989	0,99990	0,99990	0,99990	0,99991	0,99992	0,99992	0,99992	0,99992	0,99992
3,8	0,99993	0,99993	0,99993	0,99994	0,99994	0,99994	0,99994	0,99995	0,99995	0,99995
3,9	0,99995	0,99995	0,99996	0,99996	0,99996	0,99996	0,99996	0,99996	0,99997	0,99997

 

Tables de valeurs des quantiles

Les deux tables suivantes donnent les valeurs du quantile $q_{p}$ de la loi normale centrée réduite ${\mathcal {N}}(0,1)$ défini par $q_{p}=\Phi ^{-1}(p)$.

Les valeurs en début de ligne donne la première partie de la variable, les valeurs en début de colonne donne la deuxième partie. Ainsi la case de la deuxième ligne et troisième colonne donne : $q_{0{,}62}=\Phi ^{-1}(0{,}62)=0{,}3055$.

Table de valeurs des quantiles $q_{p}$

$q_{p}$	0,00	0,01	0,02	0,03	0,04	0,05	0,06	0,07	0,08	0,09
0,50	0,0000	0,0251	0,0502	0,0753	0,1004	0,1257	0,1510	0,1764	0,2019	0,2275
0,60	0,2533	0,2793	0,3055	0,3319	0,3585	0,3853	0,4125	0,4399	0,4677	0,4959
0,70	0,5244	0,5534	0,5828	0,6128	0,6433	0,6745	0,7063	0,7388	0,7722	0,8064
0,80	0,8416	0,8779	0,9154	0,9542	0,9945	1,036	1,080	1,126	1,175	1,227
0,90	1,282	1,341	1,405	1,476	1,555	1,645	1,751	1,881	2,054	2,326

Cette table donne les valeurs des quantiles pour p grand.

Table de valeur des quantiles $q_{p}$

p	0,975	0,995	0,999	0,9995	0,9999	0,99995	0,99999	0,999995
$q_{p}$	1,9600	2,5758	3,0902	3,2905	3,7190	3,8906	4,2649	4,4172

 

Les tables sont données pour les valeurs positives de la loi normale centrée réduite. Grâce aux formules de la fonction de répartition, il est possible d'obtenir d'autres valeurs.

Les valeurs négatives de la fonction de répartition sont données par la formule Φ(–x) = 1 – Φ(x). Par exemple : $\Phi (-1{,}07)=\mathbb {P} [X\leq -1{,}07]\approx 1-0{,}85769=0{,}14231$ pour $X\sim {\mathcal {N}}(0,1)$.

Les valeurs de la fonction de répartition de la loi générale s'obtiennent par la formule $F(y)=\Phi ({\frac {y-\mu }{\sigma }})$. Par exemple : $F(12{,}14)=\mathbb {P} [Y\leq 12{,}14]=\mathbb {P} \left[{\frac {Y-10}{2}}\leq {\frac {12{,}14-10}{2}}\right]=\mathbb {P} [X\leq 1{,}07]=\Phi (1{,}07)\approx 0{,}85769$, pour $Y\sim {\mathcal {N}}(10{,}2^{2})$.

La table de valeurs permet également d'obtenir la probabilité qu'une variable aléatoire de loi normale $X\sim {\mathcal {N}}(0,1)$ appartienne à un intervalle donné [a , b] par la formule : $\mathbb {P} \left[X\in [a,b]\right]=\mathbb {P} [X\leq b]-\mathbb {P} [X<a]=\Phi (b)-\Phi (a)$. Par exemple :

• $\mathbb {P} [X\geq 1{,}07]=1-\mathbb {P} [X<1{,}07]=1-\mathbb {P} [X\leq 1{,}07]\approx 0{,}14231$ pour $X\sim {\mathcal {N}}(0,1)$ ;
• $\mathbb {P} [0\leq X\leq 1{,}07]=\Phi (1{,}07)-\Phi (0)=\Phi (1{,}07)-0{,}5\approx 0{,}85769-0{,}5=0{,}35769$ pour $X\sim {\mathcal {N}}(0,1)$.

Plages de normalité, intervalles de confiance

Un des intérêts de calculer des probabilités sur des intervalles est l'utilisation des intervalles de confiance pour les tests statistiques. Une loi normale est définie par deux valeurs : sa moyenne μ et son écart type σ. Ainsi il est utile de s'intéresser aux intervalles du type [μ – rσ, μ + rσ]. $\mathbb {P} [\mu -r\sigma \leq Y\leq \mu +r\sigma ]=\Phi (r)-(1-\Phi (r))=2\Phi (r)-1$ pour $Y\sim {\mathcal {N}}(\mu ,\sigma ^{2})$.

Table de valeurs des intervalles de confiance

La courbe en cloche est la densité de probabilité. Les surfaces des zones colorées sous la courbe correspondent aux probabilités des intervalles [μ – rσ, μ + rσ].

La table suivante s'obtient grâce aux tables précédentes et donne les probabilités : $\mathbb {P} _{r}=\mathbb {P} [\mu -r\sigma \leq Y\leq \mu +r\sigma ]=2\Phi (r)-1{\text{ pour }}Y\sim {\mathcal {N}}(\mu ,\sigma ^{2})$.

Table de valeurs des intervalles de confiance

r	0,0	0,5	1,0	1,5	2,0	2,5	3,0	3,5
$\mathbb {P} _{r}$	0,00	0,3829	0,6827	0,8664	0,9545	0,9876	0,9973	0,9995

 

Représentation d'une boîte à moustaches et le lien avec les quantiles d'une loi normale.

Cette table de valeurs des intervalles de confiance permet d'obtenir les plages de normalité pour un niveau de confiance donné. Pour $Y\sim {\mathcal {N}}(\mu ,\sigma ^{2})$, le tableau donne :

• $\mathbb {P} (\mu -\sigma \leq Y\leq \mu +\sigma )\approx 0{,}6827$.L'intervalle [μ – σ, μ + σ] est la plage de normalité au niveau de confiance 68 % ;
• $\mathbb {P} (\mu -0{,}5H\leq Y\leq \mu +0{,}5H)\approx 0{,}76$.L'intervalle [μ – 0,5H, μ + 0,5H], H étant la largeur à mi-hauteur, est la plage de normalité au niveau de confiance 76 % ;
• $\mathbb {P} (\mu -2\sigma \leq Y\leq \mu +2\sigma )\approx 0{,}9545$

L'intervalle [μ – 2σ, μ + 2σ] est la plage de normalité au niveau de confiance 95.4 % ;

• $\mathbb {P} (\mu -1.96\sigma \leq Y\leq \mu +1.96\sigma )\approx 0{,}95$

L'intervalle [μ – 1.96σ, μ + 1.96σ] est la plage de normalité au niveau de confiance 95 % (voir 97,5ème centile)

• $\mathbb {P} (\mu -3\sigma \leq Y\leq \mu +3\sigma )\approx 0{,}9973$L'intervalle [μ – 3σ, μ + 3σ] est la plage de normalité au niveau de confiance 99 %.

Inversement, lorsque la valeur de la probabilité α ∈ [0 ; 1] est fixée, il existe une unique valeur $r>0$ telle que : $\mathbb {P} (\mu -r\sigma \leq Y\leq \mu +r\sigma )=2\Phi (r)-1=\alpha$. L'intervalle [μ – rσ, μ + rσ] est appelé plage de normalité ou intervalle de confiance au niveau de confiance α. Pour une loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$ et le seuil α donnés, la méthode pour retrouver cette valeur r consiste à utiliser le tableau de valeur des quantiles (ci-dessus) pour trouver la valeur r telle que Φ(r) = α + 1/2 ; l'intervalle de confiance est alors [μ – rσ, μ + rσ].

Par exemple, la plage de normalité au niveau de confiance 95 % d'une loi normale ${\mathcal {N}}(10{,}2^{2})$ est l'intervalle [10 – 2r, 10 + 2r] où r vérifie Φ(r) = 0,95 + 1/2 = 0,975, soit r = q_0,975 ≈ 1,96, l'intervalle est donc : [6,08 ; 13,92] aux arrondis près.

Liens avec d'autres lois

Grâce à son rôle central parmi les lois de probabilité et dans les applications, les lois normales possèdent beaucoup de liens avec les autres lois. Certaines lois sont même construites à partir d'une loi normale pour mieux correspondre aux applications.

Lois usuelles

Différentes lois du $\chi$ et $\chi ^{2}$

Lois	en fonction de variables de loi normale
loi du χ²	$\sum _{i=1}^{k}\left({\frac {X_{i}-\mu _{i}}{\sigma _{i}}}\right)^{2}$
loi du χ² non centrée	$\sum _{i=1}^{k}\left({\frac {X_{i}}{\sigma _{i}}}\right)^{2}$
loi du χ	${\sqrt {\sum _{i=1}^{k}\left({\frac {X_{i}-\mu _{i}}{\sigma _{i}}}\right)^{2}}}$
loi du χ non centrée	${\sqrt {\sum _{i=1}^{k}\left({\frac {X_{i}}{\sigma _{i}}}\right)^{2}}}$

Lois unidimensionnelles

• Si une variable aléatoire $X$ suit la loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$, alors la variable aléatoire $\exp(X)$ suit la loi log-normale.
• Si U et V sont deux variables aléatoires indépendantes de loi uniforme sur [0, 1], alors les deux variables aléatoires $X={\sqrt {-2\ln(U)}}\,\cos(2\pi V)$ et $Y={\sqrt {-2\ln(U)}}\,\sin(2\pi V)$ sont de loi normale centrée réduite. De plus X et Y sont indépendantes. Ces deux formules sont utilisées pour simuler la loi normale.
• Si les variables $X_{1},X_{2},\dots ,X_{n}$ sont indépendantes et de loi commune ${\mathcal {N}}(0,1)$, alors la somme de leur carré : $\sum _{k=1}^{n}X_{k}^{2}$ suit une loi du χ² à n degrés de liberté : $\chi ^{2}(n)$. La formule s'étend pour des variables normales non centrées et non réduites. De plus, le même type de lien existe avec la loi du χ² non centrée, la loi du χ et la loi du χ non centrée (voir le tableau ci-contre).
• Si la variable U suit la loi normale centrée réduite : ${\mathcal {N}}(0,1)$, si V suit une loi du χ² à n degrés de liberté $\chi ^{2}(n)$ et si U et V sont indépendantes, alors la variable ${\frac {U}{\sqrt {\frac {V}{n}}}}$ suit une loi de Student à n degrés de liberté : $t(n)$.
• Si $X$ est une variable aléatoire de loi normale centrée réduite et $U$ de loi uniforme sur [0, 1], alors ${\frac {X}{U}}$ est de loi dite de Slash.
• Pour une variable aléatoire $X$ de loi normale centrée réduite ${\mathcal {N}}(0,1)$, la variable $\mathrm {signe} (X)|X|^{p}$ est de loi normale puissance p. Pour $p=1$, cette variable est de loi normale centrée réduite.
• Si $Z_{1}$ et $Z_{2}$ sont deux variables aléatoires indépendantes de loi normale centrée réduite, alors le quotient ${\frac {Z_{1}}{Z_{2}}}$ suit la loi de Cauchy de paramètre 0 et 1 : $Cau(0,1)$. Dans le cas où $Z_{1}$ et $Z_{2}$ sont deux gaussiennes quelconques (non centrées, non réduites), le quotient ${\frac {Z_{1}}{Z_{2}}}$ suit une loi complexe dont la densité s'exprime en fonction de polynômes d'Hermite (l'expression exacte est donnée par Pham-Gia en 2006).

Lois multidimensionnelles

• Il existe une version multidimensionnelle d'une loi normale, appelée loi normale multidimensionnelle, loi multinormale ou loi de Gauss à plusieurs variables. Lorsque $X_{1},X_{2},\dots ,X_{n}$ sont des variables aléatoires de lois normales, alors la loi de probabilité du vecteur aléatoire $(X_{1},X_{2},\dots ,X_{n})$ est de loi normale multidimensionnelle. Sa densité de probabilité prend la même forme que la densité d'une loi normale mais avec une écriture matricielle. Si le vecteur aléatoire $(X_{1},X_{2})$ est de loi normale multidimensionnelle ${\mathcal {N}}(\mu ,\mathbf {\Sigma } )$ où μ est le vecteur des moyennes et $\mathbf {\Sigma }$ est la matrice de variance-covariance, alors la loi conditionnelle $(X_{1}|X_{2}=x)$ de $X_{1}$ sachant que $X_{2}=x$ est la loi normale ${\mathcal {N}}(\mu _{1|x},\sigma _{1|x})$ :Si $(X_{1},X_{2})\sim {\mathcal {N}}\left(\left({\begin{matrix}\mu _{1}\\\mu _{2}\end{matrix}}\right),\left({\begin{matrix}\sigma _{11}&\sigma _{12}\\\sigma _{21}&\sigma _{22}\end{matrix}}\right)\right)$, alors $(X_{1}|X_{2}=x)\sim {\mathcal {N}}(\mu _{1|x},\sigma _{1|x})$ avec $\mu _{1|x}=\mu _{1}+{\frac {\sigma _{12}}{\sigma _{22}}}\left(x-\mu _{2}\right)$ et $\sigma _{1|x}=\sigma _{11}-{\frac {\sigma _{12}\sigma _{21}}{\sigma _{22}}}$.
• La loi de la norme euclidienne d'un vecteur dont les coordonnées sont indépendantes et de lois normales centrées, de même variance, est une loi de Rayleigh.

Il est à noter que la loi inverse-gaussienne et loi inverse-gaussienne généralisée n'ont pas de lien avec une formule simple créée à partir de variables de loi normale, mais ont une relation avec le mouvement brownien.

Lois normales généralisées

Plusieurs généralisations de la loi normale ont été introduites afin de changer sa forme, son asymétrie, son support, etc.

Un nouveau paramètre $\beta >0$ dit de forme peut être introduit dans une loi normale pour obtenir une loi normale généralisée. Cette famille de lois contient les lois normales, c'est le cas pour $\beta =2$, mais également la loi de Laplace pour $\beta =1$. La nouvelle densité de probabilité est donnée par $f(x)={\frac {\beta }{2\alpha \Gamma (1/\beta )}}\operatorname {e} ^{-\left({\frac {|x-\mu |}{\sigma }}\right)^{\beta }}$.

Il existe une manière de changer l'asymétrie d'une loi normale afin d'obtenir une loi dite loi normale asymétrique (skew normal distribution en anglais). L'introduction d'un paramètre $\lambda \in \mathbb {R}$ permet d'obtenir une loi normale lorsque $\lambda =0$, une asymétrie vers la droite lorsque $\lambda >0$ et une asymétrie vers la gauche lorsque $\lambda <0$. La densité de cette loi est donnée par : $f(x)=2\varphi (x)\Phi (\lambda x)$.

Afin de changer le support d'une loi normale et notamment de le rendre borné, une modification possible de cette loi est de la tronquer. Elle est alors changée d'échelle pour que les parties coupées se répartissent sur l'ensemble des valeurs gardées (à la différence de la loi repliée, voir ci-dessous). La loi normale centrée réduite tronquée en –T et en T a pour support l'intervalle $[-T,T]$ et sa fonction de densité se définit par : $f(x)={\begin{cases}{\frac {\varphi (x)}{2\Phi (T)-1}}&{\text{ si }}x\in [-T,T]\\0&{\text{ sinon }}.\end{cases}}$

Il est également possible de tronquer une loi normale d'un seul côté. Elle est alors appelée « loi normale rectifiée ». Si une variable aléatoire $X$ suit la loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$, alors $\max(X,0)$ suit la loi normale rectifiée.

Une autre manière de changer le support de la loi normale est de « replier » la densité à partir d'une valeur, la loi obtenue est une loi normale repliée. Les valeurs retirées, par exemple $\left]-\infty ,0\right[$, sont alors réparties proche de la valeur charnière, 0 ici (à la différence de la loi tronquée, voir ci-dessus). La densité de probabilité de la loi normale repliée en 0 est donnée par : $f(x)={\begin{cases}{\frac {1}{\sigma {\sqrt {2\pi }}}}\,\exp \left(-{\frac {(x+\mu )^{2}}{2\sigma ^{2}}}\right)+{\frac {1}{\sigma {\sqrt {2\pi }}}}\,\exp \left(-{\frac {(x-\mu )^{2}}{2\sigma ^{2}}}\right)&{\text{ si }}x\geq 0\\0&{\text{ sinon.}}\end{cases}}$

Une version généralisée de la loi log-normale permet d'obtenir une famille de lois comprenant les lois normales comme cas particulier. La famille est définie à partir de trois paramètres : un paramètre de position μ, un paramètre d'échelle σ et un paramètre de forme $\kappa \in \mathbb {R}$. Lorsque $\kappa =0$, cette loi log-normale généralisée est la loi normale. La densité est donnée par : $f(x)={\frac {\varphi (y)}{\alpha -\kappa (x-\xi )}}$, où $y={\begin{cases}-{\frac {1}{\kappa }}\log \left[1-{\frac {\kappa (x-\xi )}{\alpha }}\right]&{\text{si }}\kappa \neq 0\\{\frac {x-\xi }{\alpha }}&{\text{si }}\kappa =0.\end{cases}}$

Différentes formes pour la densité de la loi normale généralisée.	Différentes formes pour la densité de la loi normale asymétrique.
Loi normale centrée réduite tronquée en 1,5 pour la courbe rouge et en 2,5 pour la courbe bleue.	En vert, la densité de la loi normale repliée en 0.
Différentes formes pour la densité de la loi log-normale.

Constructions à partir d'une loi normale

Mélange de lois

En bleu : densité d'une combinaison linéaire de deux densités normales.

Un mélange gaussien est une loi de probabilité dont la densité est définie par une combinaison linéaire de deux densités de lois normales. Si l'on note $f_{1}$ la densité de ${\mathcal {N}}(\mu _{1},\sigma _{1}^{2})$ et $f_{2}$ la densité de ${\mathcal {N}}(\mu _{2},\sigma _{2}^{2})$, alors $\lambda f_{1}+(1-\lambda )f_{2}$ est la densité d'une loi de probabilité dite mélange gaussien.

Il ne faut pas confondre la combinaison linéaire de deux variables aléatoires indépendantes de loi normale, qui reste une variable gaussienne, et la combinaison linéaire de leurs deux densités, qui permet d'obtenir une loi qui n'est pas une loi normale.

Les modes des deux lois normales sont donnés par μ₁ et μ₂, le mélange gaussien est alors une loi bimodale. Ses maxima locaux sont proches de mais non égaux aux valeurs μ₁ et μ₂.

Généralités

Il est possible de construire d'autres densités de probabilité grâce à la densité $\varphi$ de la loi normale centrée réduite. Harald Cramér énonce en 1926 un résultat général : si une densité de probabilité $g$ est deux fois dérivable, si l'intégrale $\int (g''(x))^{2}{\rm {e}}^{x^{2}/2}\,\mathrm {d} x$ converge et si $\lim _{+\infty }g=\lim _{-\infty }g=0$, alors la fonction $g$ peut être développée en une série absolument et uniformément convergente en fonction des dérivées de la densité de la loi normale centrée réduite et des polynômes d'Hermite $H_{k}$ : $g(x)=\sum _{k=0}^{\infty }{\frac {1}{k!}}\varphi ^{(k)}(x)\int g(y)H_{k}(y)~{\rm {d}}y$.

Utilisations

Historiquement, les lois normales sont introduites lors d'études d'objets célestes ou de jeux de hasard. Elles sont ensuite étudiées et généralisée mathématiquement puis elles sont utilisées dans de nombreuses autres applications : en mathématiques, dans d'autres sciences exactes, dans des sciences plus appliquées ou des sciences humaines et sociales. Voici une sélection d'exemples.

Balistique

Au XIX^e siècle, pour améliorer les précisions des tirs de l'artillerie, de nombreux tirs de canons sont réalisés. Il est observé que la direction et la portée sont assimilables à des lois normales. Cette compréhension permet de mieux entraîner les servants pour régler les tirs. Ces lois normales proviennent de différents facteurs comme les conditions climatiques, mais également de l'usure du matériel militaire. La dispersion des points d'impact, et donc de la loi, renseigne sur l'état du matériel et sur le nombre éventuel de tirs anormaux. L'ajustement à une loi normale est alors effectué par le test de Lhoste sur une série de 200 tirs. Le mathématicien Jules Haag applique la méthode pour 2 680 tirs de différentes portées et de différentes directions.

Quotient intellectuel

Le quotient intellectuel (QI) a pour objectif de donner une valeur numérique à l'intelligence humaine. En 1939, David Wechsler donne une définition à ce quotient de manière statistique. Une note de 100 est donnée à la moyenne des valeurs obtenues dans une population de même âge et 15 points sont retranchés pour un écart égal à l'écart type obtenu à partir des valeurs de la population testée. Pour cette raison, en pratique, la courbe de répartition du QI est modélisée par la courbe en cloche de la loi normale centrée en 100 et d'écart type 15 : ${\mathcal {N}}(100,\ 15^{2})$. Cependant cette modélisation est remise en cause par certains scientifiques. En effet, les résultats des tests dépendraient des classes sociales ; la population ne serait donc plus homogène, c'est-à-dire que la propriété d'indépendance des individus ne serait pas vérifiée. Le quotient intellectuel ne serait alors qu'une approximation de mesure de l'intelligence humaine dont on ne connaît pas l'erreur.^{[réf. nécessaire]}

Anatomie humaine

Exemple de courbe de croissance du poids.

Un caractère observable et mesurable dans une population d'individus comparables a souvent une fréquence modélisée par une loi normale. C'est le cas par exemple de la taille humaine pour un âge donné (en séparant les hommes et les femmes), de la taille des becs dans une population d'oiseaux comme les pinsons de Darwin étudiés par Darwin. Plus précisément, un caractère mesurable dans une population peut être modélisé à l'aide d'une loi normale s'il est codé génétiquement par de nombreux allèles ou par de nombreux loci ou si le caractère dépend d'un grand nombre d'effets environnementaux.

Les courbes de croissance données par l'OMS, et présentes par exemple dans les carnets de santé, sont issues de modélisations grâce à une loi normale. Grâce à une étude détaillée des centiles mesurés dans une population d'âge fixé et grâce à des tests statistiques d'adéquation, les répartitions du poids et de la taille par tranche d'âge ont été modélisées par des lois de probabilité. Parmi ces lois on retrouve les lois normales, la loi normale de Box-Cox (en) (généralisation de la loi normale), la loi Student de Box-Cox (généralisation de la loi normale de Box-Cox) ou encore la loi exponentielle-puissance de Box-Cox. Graphiquement, pour chaque âge, c'est-à-dire pour chaque axe vertical, la médiane m est représentée (elle donne la courbe centrale) et les deux valeurs de m + σ et m – σ où σ est l'écart type, donnent les deux courbes et ainsi représentent l'évolution d'un intervalle de fluctuation.^{[réf. nécessaire]}

Traitement du signal et mesures physiques

Un filtre gaussien a été appliqué à l'image du haut, issue d'un journal, pour obtenir l'image du bas, plus lisse, moins granuleuse.

Lorsqu'un signal est transmis, une perte d'information apparaît à cause du moyen de transmission ou du décodage du signal. Lorsqu'une mesure physique est effectuée, une incertitude sur le résultat peut provenir d'une imprécision de l'appareil de mesure ou d'une impossibilité à obtenir la valeur théorique. Une méthode pour modéliser de tels phénomènes est de considérer un modèle déterministe (non aléatoire) pour le signal ou la mesure et d'y ajouter ou multiplier un terme aléatoire qui représente la perturbation aléatoire, parfois appelée erreur ou bruit. Dans beaucoup de cas cette erreur additive est supposée de loi normale, de loi log-normale dans le cas multiplicatif. C'est le cas, par exemple, pour la transmission d'un signal à travers un câble électrique. Lorsque le processus dépend du temps, le signal ou la mesure est alors modélisé grâce à un bruit blanc (voir ci-dessus).

En traitement d'images, une loi normale est utilisée pour améliorer les images et notamment diminuer le bruit, c'est-à-dire les imperfections de l'image. Un lissage grâce à un filtre gaussien est alors utilisé.^{[réf. nécessaire]}

Économie

Les prix de certaines denrées sont données par une bourse, c'est le cas du cours du blé, du coton brut ou de l'or. Au temps $t$, le prix $Z(t)$ évolue jusqu'au temps $t+T$ par l'accroissement $Z(t+T)-Z(t)$. En 1900, Louis Bachelier postule que cet accroissement suit une loi normale de moyenne nulle et dont la variance dépend de $t$ et $T$. Cependant ce modèle satisfait peu l'observation faite des marchés financiers. D'autres mathématiciens proposent alors d'améliorer ce modèle en supposant que c'est l'accroissement $\ln Z(t+T)-\ln Z(t)$ qui suit une loi normale, c'est-à-dire que l'accroissement du prix suit une loi log-normale. Cette hypothèse est à la base du modèle et de la formule de Black-Scholes utilisés massivement par l'industrie financière.

Ce modèle est encore amélioré, par Benoît Mandelbrot notamment, en supposant que l'accroissement suit une loi stable (la loi normale est un cas particulier de loi stable). Il apparaît alors le mouvement brownien dont l'accroissement est de loi normale et le processus de Lévy (stable) dont l'accroissement stable pour modéliser les courbes des marchés.

Mathématiques

Bruit blanc gaussien unidimensionnel.

Les lois normales sont utilisées dans plusieurs domaines des mathématiques. Le bruit blanc gaussien est un processus stochastique tel qu'en tout point, le processus est une variable aléatoire de loi normale indépendante du processus aux autres points. Le mouvement brownien $(B(t),t\geq 0)$ est un processus stochastique dont les accroissements sont indépendants, stationnaires et de loi normale. Notamment pour une valeur $t>0$ fixée, la variable aléatoire $B(t)$ suit la loi normale ${\mathcal {N}}(0,t)$. Ce processus aléatoire possède de nombreuses applications, il fait un lien entre l'équation de la chaleur et la loi normale. Lorsque l'extrémité d'une tige métallique est chauffée pendant un court instant, la chaleur se propage le long de la tige sous la forme d'une courbe en cloche.

Les lois normales ont également des applications dans des domaines mathématiques non aléatoires comme la théorie des nombres. Tout nombre entier n peut s'écrire comme un produit de puissances de nombres premiers. Notons $\omega (n)$ le nombre de nombres premiers différents dans cette décomposition. Par exemple, puisque $60=2^{2}\times 3\times 5$, $\omega (60)=3$. Le théorème d'Erdős-Kac assure que cette fonction $n\mapsto \omega (n)$ pour $n\leq N$ est apparentée à la densité d'une loi normale ${\mathcal {N}}\left(\ln \ln(N),{\sqrt {\ln \ln(N)}}\right)$. C'est-à-dire que pour un grand nombre de l'ordre de $1000000000=10^{9}$, il y a une forte probabilité pour que le nombre de diviseurs premiers soit 3, puisque $\ln \ln(10^{9})\approx 3{,}03$.

Tests et estimations

Critères de normalité

Quatre valeurs ainsi que la droite de Henry représentées sur un papier gausso-arithmétique.

Il est important de savoir si des valeurs sont distribuées suivant une loi normale. Quelques critères peuvent être étudiés avant de réaliser un test statistique (voir la section Tests de normalité ci-dessous).

Le premier critère, le plus simple, consiste à tracer le diagramme en bâtons de la distribution et à vérifier visuellement si le diagramme est en forme de « cloche ». Ce critère, subjectif, permet cependant d'éliminer une partie des distributions jugées alors non gaussiennes.

[réf. nécessaire]

De manière plus précise, l'utilisation des plages de normalité permet de comparer avec les fréquences observées facilement calculables. Le critère consiste à utiliser les plages de normalité ou intervalles de confiance. Lorsque des valeurs suivent la loi normale :

• 68 % d'entre elles sont dans l'intervalle $[{\overline {x}}-\sigma ,{\overline {x}}+\sigma ]$ ;
• 95 % d'entre elles sont dans l'intervalle $[{\overline {x}}-2\sigma ,{\overline {x}}+2\sigma ]$ ;
• 99,7 % d'entre elles sont dans l'intervalle $[{\overline {x}}-3\sigma ,{\overline {x}}+3\sigma ]$.

Si ce n'est pas le cas, le choix de modéliser la loi des valeurs observées par la loi normale n'est pas conseillé.

[réf. nécessaire]

La droite de Henry permet de faire un ajustement des valeurs observées avec une loi normale. C'est-à-dire qu'en représentant la droite de Henry, il est possible de porter un diagnostic sur la nature normale ou non de la distribution et, dans le cas où celle-ci a des chances d'être normale, elle permet d'en déterminer la moyenne et l'écart type. Les valeurs $(x_{i},i\leq n)$ sont observées et représentées par leur fonction de répartition empirique $F_{n}$. Elles sont gaussiennes si les points $(x_{i},F_{n}(x_{i}))$ représentés sur un papier gausso-arithmétique sont alignés suivant une droite dite de Henri. Un papier gausso-arithmétique est gradué avec une échelle arithmétique en abscisse et graduée suivant l'inverse de la fonction de répartition de la loi normale centrée réduite $\Phi ^{-1}$ en ordonnée.

Ces critères sont nécessaires mais non suffisants. Cependant, il ne suffit pas de remplir les critères pour affirmer que les valeurs suivent la loi normale.^{[réf. nécessaire]}

Tests de normalité

Grâce à son rôle dans le théorème central limite, les lois normales se retrouvent dans de nombreux tests statistiques dits gaussiens ou asymptotiquement gaussiens. L'hypothèse dite de normalité est faite sur une loi a priori dans un test d'adéquation pour indiquer que cette loi suit, approximativement, une loi normale. Il existe plusieurs tests de normalité.

• Un test du χ² d'adéquation à une loi normale est possible pour tester si une série de k valeurs observées suit une loi normale. Dans ce type de test, l'hypothèse nulle est : la distribution observée peut être approchée par une loi normale. Après avoir regroupé les k valeurs observées en classes, il faut calculer les probabilités qu'une variable aléatoire de loi normale appartienne à chaque classe en estimant les paramètres de la loi grâce aux valeurs observées. Ces probabilités peuvent être obtenues avec les tables numériques d'une loi normale. Si l'hypothèse nulle est vraie, la statistique du χ² calculée à partir des valeurs observées et des probabilités précédentes suit une loi du χ². Le nombre de degré de liberté est k – 1 si la moyenne et l'écart type sont connus, k – 2 si l'un des deux paramètres est inconnu, ou k – 3 si les deux paramètres sont inconnus. L'hypothèse nulle est rejetée si la statistique du χ² est supérieure à la valeur obtenue grâce à la table de la loi du χ² au seuil α.
• Le test de Lilliefors est basé sur la comparaison entre la fonction de répartition d'une loi normale et la fonction de répartition empirique, c'est une adaptation du test de Kolmogorov-Smirnov. Les avis sont partagés sur la puissance de ce test, il est performant autour de la moyenne mais l'est moins pour la comparaison des queues de distribution. Les valeurs observées $(x_{i},i\leq n)$ sont rangées par ordre croissant $(x_{(i)},i\leq n)$, les valeurs $F_{i}=\Phi \left((x_{(i)}-{\overline {x}})/s\right)$ sont les fréquences théoriques de la loi normale centrée réduite associées aux valeurs standardisées. Si la statistique :$D=\max _{i=1,\dots ,n}\left(F_{i}-{\frac {i-1}{n}};{\frac {i}{n}}-F_{i}\right)$est supérieure à une valeur critique calculée grâce au seuil α et à la taille de l'échantillon, alors l'hypothèse de normalité est rejetée au seuil α.
• Le test d'Anderson-Darling est une autre version du test de Kolmogorov-Smirnov mieux adaptée à l'étude des queues de distribution. En reprenant les mêmes notations que le test de Lilliefors, si la statistique :$A=-n-{\frac {1}{n}}\sum _{i=1}^{n}(2i-1)\left(\ln(F_{i})+\ln(1-F_{n-i+1})\right)$est supérieure à une valeur critique calculée grâce au seuil α et à la taille de l'échantillon, alors l'hypothèse de normalité est rejetée au seuil α.
• Le test de D'Agostino est basé sur les coefficients de symétrie et d'aplatissement. Il est particulièrement efficace à partir de $n\geq 20$ valeurs observées. Même si l'idée de ce test est simple, les formules sont plus compliquées à écrire. L'idée est de construire des modifications des coefficients de symétrie et d'aplatissement pour obtenir des variables $z_{1}$ et $z_{2}$ de loi normale centrée réduite. Il faut alors effectuer un test du χ² avec la statistique $z_{1}^{2}+z_{2}^{2}$.
• Le test de Jarque-Bera est également basé sur les coefficients de symétrie et d'aplatissement. Ce test n'est intéressant que pour un nombre élevé de valeurs observées. En considérant les deux estimateurs :$b_{1}={\frac {{\frac {1}{n}}\sum _{i=1}^{n}(x_{i}-{\overline {x}})^{3}}{\left({\frac {1}{n}}\sum _{i=1}^{n}(x_{i}-{\overline {x}})^{2}\right)^{\frac {3}{2}}}}$ et $b_{2}={\frac {{\frac {1}{n}}\sum _{i=1}^{n}(x_{i}-{\overline {x}})^{4}}{\left({\frac {1}{n}}\sum _{i=1}^{n}(x_{i}-{\overline {x}})^{2}\right)^{2}}},$comme précédemment, il faut effectuer un test du χ² avec la statistique $T=n\left(b_{1}^{2}/6+(b_{2}-3)^{2}/24\right)$.
• Le test de Shapiro-Wilk (proposé en 1965) est efficace pour les petits échantillons de moins de 50 valeurs. Les valeurs observées $(x_{i},i\leq n)$ sont rangées par ordre croissant $(x_{(i)},i\leq n)$ et des coefficients $a_{i}$ sont calculés à partir de quantiles, moyenne, variance et covariance d'une loi normale. Si la statistique$W={\frac {\left(\sum _{i=1}^{[n/2]}a_{i}\left(x_{(n-i+1)}-x_{(i)}\right)\right)^{2}}{\sum _{i=1}^{n}(x_{i}-{\overline {x}})^{2}}}$est inférieure à une valeur critique calculée grâce au seuil α et à la taille de l'échantillon, alors l'hypothèse de normalité est rejetée au seuil α.

Estimations des paramètres

Lorsqu'un phénomène aléatoire est observé et qu'il est considéré comme pouvant être modélisé par une loi normale, une des questions que l'on peut se poser est : que valent les paramètres μ et σ de la loi normale ${\mathcal {N}}(\mu ,\sigma ^{2})$ ? Une estimation est alors à effectuer. Les observations récupérées lors de l'observation du phénomène sont notées par des variables aléatoires $X_{1},X_{2},\dots ,X_{n}$, les notations de la moyenne arithmétique et de la moyenne des carrés sont également utiles : $S_{n}={\frac {1}{n}}(X_{1}+X_{2}+\dots +X_{n})$ et $T_{n-1}^{2}={\frac {1}{n-1}}\sum _{k=1}^{n}(X_{k}-S_{n})^{2}$. Ces deux valeurs sont respectivement des estimateurs de la moyenne et de la variance qui se calculent à partir des valeurs observées. Puisque les variables $X_{1},X_{2},\dots ,X_{n}$ sont de loi normale, alors $S_{n}$ est de loi ${\mathcal {N}}(\mu ,{\frac {\sigma ^{2}}{n}})$ et ${\frac {n-1}{\sigma ^{2}}}T_{n-1}^{2}$ est indépendante de S_n et de loi du χ² à n – 1 degrés de liberté.

Estimation de la moyenne μ (lorsque l'écart type σ est connu)

Une méthode est de chercher un intervalle de confiance à un seuil α autour de la moyenne théorique μ. En utilisant les quantiles d'ordre ${\frac {\alpha }{2}}$ et $1-{\frac {\alpha }{2}}$, la formule définissant les quantiles permet d'obtenir : $\mathbb {P} \left(S_{n}-{\frac {\sigma }{\sqrt {n}}}q_{\alpha /2}\leq \mu \leq S_{n}+{\frac {\sigma }{\sqrt {n}}}q_{\alpha /2}\right)\geq 1-\alpha$. Grâce aux valeurs observées et aux tables numériques de la loi normale centrée réduite (voir la table), il est alors possible de donner les valeurs numériques de l'intervalle $\left[S_{n}-{\frac {\sigma }{\sqrt {n}}}q_{\alpha /2},S_{n}+{\frac {\sigma }{\sqrt {n}}}q_{\alpha /2}\right]$, intervalle de confiance pour μ au seuil α.

Estimation de la moyenne μ (lorsque l'écart type σ est inconnu)

Une méthode est d'utiliser une variable intermédiaire qui peut s'écrire à l'aide de nouvelles variables aléatoires $U={\frac {S_{n}-\mu }{\sigma /{\sqrt {n}}}}$ de loi ${\mathcal {N}}(0,1)$ et $V={\frac {n-1}{\sigma ^{2}}}T_{n-1}^{2}$ de loi $\chi ^{2}(n-1)$ : ${\sqrt {n}}{\frac {S_{n}-\mu }{T_{n-1}}}={\frac {U{\sqrt {n-1}}}{\sqrt {V}}}$ est de loi de Student $t(n-1)$. En utilisant les quantiles d'ordre ${\frac {\alpha }{2}}$ et $1-{\frac {\alpha }{2}}$, la formule définissant les quantiles permet d'obtenir : $\mathbb {P} \left(S_{n}+{\frac {T_{n-1}}{\sqrt {n}}}q_{\alpha /2}\leq \mu \leq S_{n}-{\frac {T_{n-1}}{\sqrt {n}}}q_{\alpha /2}\right)\geq 1-\alpha$. Grâce aux valeurs observées et aux tables numériques des lois de Student, il est alors possible de donner les valeurs numériques de l'intervalle $\left[S_{n}+{\frac {T_{n-1}}{\sqrt {n}}}q_{\alpha /2},S_{n}-{\frac {T_{n-1}}{\sqrt {n}}}q_{\alpha /2}\right]$, intervalle de confiance pour μ au seuil α.

Estimation de l'écart type σ (lorsque la moyenne μ est inconnue)

La méthode est la même que la précédente. L'introduction de la variable aléatoire $T_{n-1}^{2}{\frac {n-1}{\sigma ^{2}}}$ de loi du χ² à n – 1 degrés de liberté permet d'obtenir : $\mathbb {P} \left({\sqrt {T_{n-1}^{2}{\frac {n-1}{q_{1-\alpha /2}}}}}\leq \sigma \leq {\sqrt {T_{n-1}^{2}{\frac {n-1}{q_{\alpha /2}}}}}\right)\geq 1-\alpha$ où $q_{1-\alpha /2}$ et $q_{\alpha /2}$ sont les quantiles de la loi du χ² à n – 1 degrés de liberté que l'on peut obtenir à partir de la table numérique du χ². L'intervalle $\left[{\sqrt {T_{n-1}^{2}{\frac {n-1}{q_{1-\alpha /2}}}}},{\sqrt {T_{n-1}^{2}{\frac {n-1}{q_{\alpha /2}}}}}\right]$ est l'intervalle de confiance au seuil α.

Simulation

Pour étudier un phénomène aléatoire dans lequel intervient une variable normale dont les paramètres sont connus ou estimés, une approche analytique est souvent trop complexe à développer. Dans un tel cas, il est possible de recourir à une méthode de simulation, en particulier à la méthode de Monte-Carlo qui consiste à générer un échantillon artificiel de valeurs indépendantes de la variable, ceci à l'aide d'un ordinateur. Les logiciels ou les langages de programmation possèdent en général un générateur de nombres pseudo-aléatoires ayant une distribution uniforme sur ]0, 1[. Il s'agit alors de transformer cette variable de loi $U(]0,1[)$ en une variable ${\mathcal {N}}(0,1)$ (l'adaptation à d'autres valeurs des paramètres ne pose pas de problème).

Approches à éviter

• De manière générale, on peut exploiter la fonction réciproque de la fonction de répartition : en l'occurrence, la variable aléatoire $\Phi ^{-1}(U)$ suit la loi normale centrée réduite ; cette méthode est cependant malcommode, faute d'expressions simples des fonctions $\Phi$ et $\Phi ^{-1}$ ; par ailleurs, les résultats sont numériquement insatisfaisants.
• Si $U_{1},U_{2},\dots ,U_{12}$ sont douze variables indépendantes de loi uniforme sur [0, 1], alors la variable $\sum _{k=1}^{12}U_{k}-6$ est de moyenne nulle et d'écart type unitaire. Ainsi, grâce au théorème central limite, cette variable suit approximativement la loi normale centrée réduite. C'est une manière simple de générer une loi normale, cependant l'approximation reste imprécise.

Approches efficientes

• Un meilleur algorithme est la méthode de Box-Muller qui utilise une représentation polaire de deux coordonnées uniformes donnée par les formules : Si ${\begin{cases}U\sim {\mathcal {U}}(0,1)\\V\sim {\mathcal {U}}(0,1)\end{cases}}$ alors ${\begin{cases}{\sqrt {-2\ln(U)}}\cos(2\pi V)\sim {\mathcal {N}}(0,1)\\{\sqrt {-2\ln(U)}}\sin(2\pi V)\sim {\mathcal {N}}(0,1)\end{cases}}$et les deux variables obtenues sont indépendantes. Cet algorithme est simple à réaliser, mais le calcul d'un logarithme, d'une racine carrée et d'une fonction trigonométrique ralentit le traitement.
• Une amélioration a été proposée par Marsaglia et Bray en 1964, en remplaçant les cosinus et sinus par les variables $V_{1}/{\sqrt {W}}$ et $V_{2}/{\sqrt {W}}$ où $V_{1}$ et $V_{2}$ sont indépendantes de loi ${\mathcal {U}}(-1,1)$ et $W=V_{1}^{2}+V_{2}^{2}$, lorsque $W<1$ (on rejette les couples qui ne vérifient pas cette dernière condition). Ainsi :${\begin{cases}V_{1}{\sqrt {-2{\dfrac {\ln W}{W}}}}\sim {\mathcal {N}}(0,1)\\V_{2}{\sqrt {-2{\dfrac {\ln W}{W}}}}\sim {\mathcal {N}}(0,1).\end{cases}}$Cet algorithme n'est pas plus lourd à mettre en œuvre et la simulation gagne en vitesse.
• Pour un grand nombre de tirages aléatoires, la méthode Ziggourat est encore plus rapide, mais sa mise en œuvre est plus complexe.

Hommages

Par son utilisation généralisée dans les sciences, une loi normale, souvent par l'utilisation de la courbe en cloche, est mise en avant dans différents contextes et est utilisée pour représenter l'universalité d'une répartition statistique, entre autres.

Francis Galton parle des lois normales dans son œuvre Natural Inheritence de 1889 en ces termes élogieux :

« Je ne connais rien d'autre si propre à frapper l'imagination que cette merveilleuse forme d'ordre cosmique donnée par la Loi de Fréquence des Erreurs... Elle règne avec sérénité et en toute abnégation au milieu de la confusion sauvage. »

— Francis Galton

En 1989, un hommage est rendu à Carl Friedrich Gauss en imprimant un billet à son effigie, la courbe en cloche est également présente sur le billet. Des pierres tombales portent le signe de la courbe en cloche, c'est le cas pour certains mathématiciens.

Le statisticien William J. Youden écrit en 1962 une explication du but et de la position des lois normales dans les sciences. Il la présente en calligramme sous forme de courbe en cloche :

THE NORMAL LAW OF ERROR STANDS OUT IN THE EXPERIENCE OF MANKIND AS ONE OF THE BROADEST GENERALIZATIONS OF NATURAL PHILOSOPHY ♦ IT SERVES AS THE GUIDING INSTRUMENT IN RESEARCHES IN THE PHYSICAL AND SOCIAL SCIENCES AND IN MEDICINE AGRICULTURE AND ENGINEERING ♦ IT IS AN INDISPENSABLE TOOL FOR THE ANALYSIS AND THE INTERPRETATION OF THE BASIC DATA OBTAINED BY OBSERVATION AND EXPERIMENT

« La loi normale des erreurs se distingue dans l'expérience de l'humanité comme une des plus larges généralisations de la philosophie naturelle ♦ Elle sert de guide dans la recherche en sciences physique et sociale, en médecine, en agriculture et en ingénierie ♦ C'est un outil indispensable pour l'analyse et l'interprétation des données de base obtenues par l'observation et l'expérience. »

Notes et références

Voir aussi

Bibliographie

 : document utilisé comme source pour la rédaction de cet article.

• (en) Milton Abramowitz et Irene Stegun, Handbook of Mathematical Functions, New York, Dover, 1972, 9^e éd., 1047 p. (ISBN 0-486-61272-4, lire en ligne), chap. 26 (« Probability Functions »), p. 927-996
• Patrick Bogaert, Probabilités pour scientifiques et ingénieurs, Paris, De Boeck, 2006, 387 p. (ISBN 2-8041-4794-0, lire en ligne)
• (en) Harald Cramér, Random Variables and Probability Distributions, Cambridge University Press, 1970, 3^e éd., 123 p. (ISBN 0-521-60486-9, lire en ligne)
• Yadolah Dodge (en), Statistique - dictionnaire encyclopédique, Springer Verlag, 2004, 637 p. (ISBN 2-287-21325-2, lire en ligne)
• Jean-Jacques Droesbeke, Michel Lejeune et Gilbert Saporta, Modèles statistiques pour données qualitatives, Technip, 2005, 295 p. (ISBN 2-7108-0855-2, lire en ligne)
• (en) Joseph Arthur Greenwood et H. O. Hartley, Guide to Tables in Mathematical Statistics, Princeton University Press, 1962, 1014 p.
• (en) Charles Miller Grinstead et James Laurie Snell, Introduction to Probability, AMS, 1997, 2^e éd., 519 p. (ISBN 0-8218-0749-8, lire en ligne)
• (en) J. R. M. Hosking et James R. Wallis, Regional Frequency Analysis: An Approach Based on L-Moments, Cambridge University Press, 1997, 224 p. (ISBN 978-0-521-43045-6, lire en ligne)
• (en) M. A. Lifshits, Gaussian Random Functions, Kluwer Academic Publishers, 1995, 339 p. (ISBN 0-7923-3385-3, lire en ligne)
• (en) NBS, A Guide to Tables of the Normal Probability Integral, U.S. Govt. Print. Off., 1952, 16 p.
• Konstantin Protassov, Analyse statistique des données expérimentales, EDP Sciences, 2002, 148 p. (ISBN 978-2-75980113-8, lire en ligne)
• Martine Quinio Benamo, Probabilités et Statistique aujourd'hui, L'Harmattan, 2005, 277 p. (ISBN 2-7475-9799-7, lire en ligne)
• (en) Mark Ridley (en), Evolution, Blackwell, 2004, 3^e éd., 751 p. (ISBN 1-4051-0345-0)
• Sheldon M. Ross, Initiation aux probabilités, PPUR, 2007, 592 p. (ISBN 978-2-88074-738-1, lire en ligne)
• (en) Stephen Stigler, Statistics on the Table, Harvard University Press, 1999, 499 p. (ISBN 0-674-83601-4, lire en ligne)
• Philippe Tassi et Sylvia Legait, Théorie des probabilités en vue des applications statistiques, Technip, 1990, 367 p. (ISBN 2-7108-0582-0, lire en ligne)
• Alain Yger et Jacques-Arthur Weil, Mathématiques appliquées, Pearson Education, 2009, 890 p. (ISBN 978-2-7440-7352-6, lire en ligne)

Articles connexes

• Loi normale multidimensionnelle
• Fonction d'erreur
• Théorème central limite

Liens externes

• Ressources relatives à la santé :
  • Medical Subject Headings
  • WikiSkripta
• Notices dans des dictionnaires ou encyclopédies généralistes :

  • Britannica
  • Den Store Danske Encyklopædi
  • Larousse
  • Store norske leksikon
• Notices d'autorité :

  • BnF (données)
  • LCCN
  • GND
  • Israël
  • Tchéquie
• (en) Generating Gaussian Random Numbers
• Dix vidéos (sur YouTube) de moins de dix minutes chacune dans lesquelles Saïd Chermak explique les propriétés de la loi normale.
• Animation en ligne de la planche de Galton.