Loi hypergéométrique

Un article de Wikipédia, l'encyclopédie libre.

Une loi hypergéométrique de paramètres n, p et A correspond au modèle suivant:

On tire simultanément n boules dans une urne contenant pA boules gagnantes et qA boules perdantes (avec q = 1 - p). On compte alors le nombre de boules gagnantes extraites et on appelle X la variable aléatoire donnant le nombre de boules gagnantes.

L'univers X(Ω) est l'ensemble des entiers de 0 à n.

La variable aléatoire suit une loi de probabilité définie par

p(k)=\frac{C_{pA}^kC_{qA}^{n-k}}{C_A^n}

Cette loi de probabilité s'appelle la loi hypergéométrique de paramètres (n ; p ; A). Il est nécessaire que p soit un réel compris entre 0 et 1, que pA soit entier et que n ≤ A

Sommaire

[modifier] Calcul de p(k)

Il s'agit d'un tirage simultané (c'est-à-dire non ordonné et sans remise) de n éléments parmi A. Tirage que l'on considère comme équiprobable.

La combinatoire permet de dire que le cardinal de l'univers est C_A^n.

L'événement « X=k » signifie que l'on a tiré k boules gagnantes parmi pA et n - k boules perdantes parmi qA. Le cardinal de cet événement est donc C_{pA}^kC_{qA}^{n-k}.

La probabilité de l'événement est donc p(k)=\frac{C_{pA}^kC_{qA}^{n-k}}{C_A^n}

[modifier] Espérance, variance et écart type

L'espérance d'une variable aléatoire suivant une loi hypergéométrique est la même que dans le cas binômiale. X\, suit une loi hypergéométrique de paramètres n, N_A, N\,, alors son espérance est \mathbb{E}(X)=\frac{nN_A}{N}=np\,. (On a posé p=\frac{N_A}{N}\,)


Démonstration

On se donne: X \sim \mathcal{H}(n,N_A,N) où N=N_A+N_B \,

(si on se rapporte à un modèle d'urnes à tirage simultané (c'est-à-dire non ordonné et sans remise). On a donc N_A\,: le nombre de boules de type "réussite" et N_B\,: le nombre de boules de type "echec".)

\mathbb{P}(X=k)=\frac{C_{N_A}^kC_{N_B}^{n-k}}{C_N^n}

Numérotons de 1 à N_A\, les boules de type "réussite" et définissons pour tout k\, compris entre 1 et N_A\, l'événement:

E_k=\,{ on a tiré parmi les n boules la boule de type "réussite" k }. Comme le nombre total X\, de boules de type "réussite" tirées est X=\sum_{k=1}^{N_A}1_{E_k}\,

(où 1 est la Fonction indicatrice de E_k\,), par linéarité de l'espérance, \mathbb{E}(X)=N_A\mathbb{P}(E_1)\,.

Evaluons maintenant \mathbb{P}(E_1)\,. En passant au complémentaire,

\mathbb{P}\bar{(E_1)}=\frac{C_{N-1}^n}{C_N^n}=\frac{(N-1)!}{n!(N-1-n)!}=\frac{n!(N-n)!}{N!}=\frac{N-n}{N}\, qui est la probabilité de ne jamais tirer une boule donnée.

Donc \mathbb{P}(E_1)=1-\frac{N-n}{N}=\frac{n}{N}=\frac{n}{N_A+N_B}\,

On en conclut donc que \mathbb{E}(X)=\frac{nN_A}{N_A+N_B}=\frac{nN_A}{N}\,

Si on pose \frac{N_A}{N}=p\, qui est exactement la probabilité d'avoir un succès, on a bien \mathbb{E}(X)=np\,.

La variance d'une variable aléatoire suivant une loi hypergéométrique de paramètres n, p, A est npq\frac{A - n}{A - 1}

L'écart type est alors \sqrt{npq}\sqrt{\frac{A - n}{A - 1}}.

[modifier] Convergence

Pour n petit devant A, la loi hypergéométrique est équivalente à une loi binomiale de paramètres n et p. En fait, on considère que, pour A grand, tirer simultanément n boules revient à effectuer n fois une épreuve de Bernoulli dont la probabilité de succès serait p (p est la proportion de boules gagnantes dans l'ensemble des boules), car il est très peu probable de retomber sur la même boule, même si on la replace dans l'urne.

En pratique, on remplace la loi hypergéométrique de paramètres (n ; p ; A) par une loi binomiale de paramètres (n ; p) dès que 10n < A.

Un exemple très classique de ce remplacement concerne le sondage. On considère fréquemment le sondage de n personnes comme n sondages indépendants alors qu'en réalité le sondage est exhaustif (on n'interroge jamais deux fois la même personne). Comme n(nombre de personnes interrogées) < A(population sondée)/10, cette approximation est légitime.

[modifier] Voir aussi