La loi binomiale compte les succès dans un nombre fixé d’épreuves indépendantes. Deux lois voisines répondent à d’autres questions, tout aussi fréquentes. La loi géométrique renverse la question : on répète une épreuve de Bernoulli de probabilité p jusqu’au premier succès, et l’on compte le nombre d’épreuves nécessaires. La loi hypergéométrique change les conditions du tirage : on prélève n objets sans remise dans une population finie de N objets dont K sont « marqués », et l’on compte les objets marqués obtenus. Les tirages ne sont alors plus indépendants : chaque objet sorti modifie la composition de ce qui reste.
Loi géométrique G(p)
Loi hypergéométrique H(N, K, n)
Question
combien d’épreuves jusqu’au premier succès ?
combien d’objets marqués parmi n tirés sans remise ?
Valeurs
1, 2, 3, … (sans limite)
de max(0, n − (N − K)) à min(n, K)
P(X = k)
(1 − p)ᵏ⁻¹ · p
C(K, k) · C(N − K, n − k) / C(N, n)
Espérance
1 / p
n · K / N
Variance
(1 − p) / p²
n · (K/N) · (1 − K/N) · (N − n)/(N − 1)
Propriété clé
absence de mémoire
proche de la binomiale quand N est grand devant n
Remarque de vocabulaire : certains auteurs appellent « loi géométrique » le nombre d’échecs avant le premier succès, qui commence à 0. Cette page compte les épreuves, succès compris : X commence à 1.
1
Attendre le premier succès
Pour que le premier succès arrive à la k-ième épreuve, il faut d’abord k − 1 échecs, puis un succès : par indépendance, P(X = k) = qᵏ⁻¹·p avec q = 1 − p. Ces probabilités forment une suite géométrique de raison q — d’où le nom. Elles décroissent toujours : la valeur la plus probable est 1, quelle que soit p. Pourtant l’attente moyenne est 1/p : six lancers en moyenne pour obtenir un 6.
La partie en cours : les épreuves s’enchaînent jusqu’au premier succès (en rouge).
Barres bleues : loi exacte P(X = k). Points orange : fréquences observées sur les parties jouées. Trait vert : l’espérance 1/p ; trait violet : la médiane.
Espérance 1/p—
Contrôle : Σ P(X ≥ k) pour k ≥ 1—autre formule de l’espérance, sommée sur la queue
Médiane : plus petit k avec P(X ≤ k) ≥ ½—
Écart-type √q / p—
Parties jouées : attente moyenne—jouez des parties
La loi géométrique est très dissymétrique : la valeur la plus probable est 1, la médiane est nettement plus petite que la moyenne, et de longues attentes restent possibles. Pour le 6, il faut au moins 4 lancers une fois sur deux, mais plus de 12 lancers environ une fois sur neuf, et plus de 25 lancers environ une fois sur cent. La probabilité d’attendre plus de k épreuves est simplement qᵏ : il faut que les k premières soient toutes des échecs.
2
Le hasard n’a pas de mémoire
Un joueur a lancé dix fois le dé sans voir de 6. Le 6 est-il « en retard », donc plus probable au prochain lancer ? Non : le dé ignore son passé. Mathématiquement, sachant que m échecs se sont déjà produits, l’attente restante suit exactement la même loi géométrique qu’au départ : P(X > m + n | X > m) = qᵐ⁺ⁿ / qᵐ = qⁿ = P(X > n). La loi géométrique est la seule loi à valeurs entières qui ait cette propriété. Comparez avec une situation sans remise : chercher l’as de pique en retournant les cartes d’un paquet de 52, où chaque carte retournée rapproche bel et bien du but.
Barres bleues : loi de l’attente totale depuis le début. Barres rouges : loi de l’attente restante, sachant que m échecs se sont déjà produits (loi conditionnelle, calculée exactement). Points orange : attentes restantes observées sur les parties simulées qui ont dépassé m.
P(succès dans les n prochaines épreuves | m échecs)—
P(succès dans les n premières épreuves)—
Plus grand écart entre les deux lois—
Attente restante moyenne : exacte et observée—
Avec le dé, les barres rouges recouvrent exactement les bleues, quel que soit m : l’écart calculé est nul. C’est l’« erreur du joueur » : croire qu’un numéro absent depuis longtemps du Loto ou de la roulette va sortir. Avec le paquet de cartes, au contraire, chaque carte retournée sans succès réduit le nombre de cartes restantes : après 40 cartes, l’as de pique est forcément parmi les 12 dernières, et l’attente restante n’est plus que de 6,5 cartes en moyenne. La mémoire vient du tirage sans remise — celui de la loi hypergéométrique.
3
Tirer sans remise : la loi hypergéométrique
Une urne contient N boules, dont K rouges. On en tire n d’un coup, ou l’une après l’autre sans les remettre. Toutes les poignées de n boules ont la même probabilité, et il y en a C(N, n). Celles qui contiennent exactement k rouges se forment en choisissant k rouges parmi K et n − k blanches parmi N − K : d’où P(X = k) = C(K, k)·C(N − K, n − k) / C(N, n). Si l’on remettait chaque boule, ce serait la loi binomiale B(n, K/N). Les deux lois ont la même moyenne, mais le tirage sans remise est moins dispersé : chaque rouge tirée rend la suivante un peu moins probable.
L’urne et la poignée en cours : les boules tirées quittent l’urne une à une.
Barres bleues : loi hypergéométrique exacte. Points rouges : loi binomiale de même moyenne (tirage avec remise). Points orange : fréquences observées sur les poignées tirées.
Écart à la binomiale (distance en variation totale)—
Poignées tirées—tirez des poignées
Le facteur (N − n)/(N − 1), appelé facteur de correction pour population finie, mesure ce que le tirage sans remise fait gagner en précision. Si l’on tire toute l’urne (n = N), il vaut 0 : on connaît K exactement. Si l’urne est immense devant la poignée, il vaut presque 1 et les deux lois se confondent. C’est pourquoi un sondage de 1 000 personnes a la même précision dans une ville de 100 000 habitants que dans un pays de 60 millions.
4
Trois applications : le Loto, les petits lots, la dame qui goûte le thé
La loi hypergéométrique s’impose chaque fois qu’on tire dans un ensemble fini sans remettre : les boules d’un jeu de loterie, les pièces d’un petit lot contrôlé en entier, les tasses d’une expérience. Dans chacun des trois cas, les probabilités sont des fractions exactes, calculées ici en nombres entiers.
Dans un loto à N numéros dont on en tire T, vous cochez T numéros. Le nombre de bons numéros de votre grille suit la loi H(N, T, T) : les numéros tirés sont les boules « rouges », votre grille est la poignée.
À gauche, la grille et le dernier tirage (cases bleues : votre grille ; cercles rouges : numéros tirés). À droite, la loi exacte du nombre de bons numéros, en échelle logarithmique ; points orange : tirages simulés.
Contrôle : somme des fractions exactes—
Tirages simulés—
Un lot de N pièces en contient D défectueuses. On en prélève n sans remise et l’on accepte le lot si l’échantillon contient au plus c défauts. Pour un petit lot, l’échantillon représente une part importante du total et la loi binomiale, qui suppose un tirage avec remise, se trompe.
Probabilité d’accepter le lot selon son nombre D de pièces défectueuses : loi hypergéométrique exacte (bleu) et approximation binomiale (rouge pointillé).
P(accepter) exacte (hypergéométrique)—
Approximation binomiale—
Contrôle : 20 000 lots simulés—
Vers 1920, à Rothamsted, la biologiste Muriel Bristol affirme reconnaître si le lait a été versé avant ou après le thé. Ronald Fisher lui présente 2m tasses, dont m « lait d’abord », dans un ordre aléatoire ; elle doit désigner les m tasses « lait d’abord ». Si elle n’a aucun talent, son choix est une poignée de m tasses tirée au hasard, et son nombre de bonnes réponses suit H(2m, m, m). C’est le premier test exact de l’histoire, décrit par Fisher en 1935.
Loi du nombre de tasses « lait d’abord » bien reconnues par quelqu’un qui répond au hasard ; en rouge, les résultats au moins aussi bons que celui observé.
p-valeur : P(au moins autant de bonnes réponses par hasard)—
Nombre de choix possibles C(2m, m)—
Conclusion au seuil de 5 %—
Au Loto à 6 numéros sur 49, une grille a une chance sur 13 983 816 d’avoir les six bons numéros, et près de 44 % de chances de n’en avoir aucun. Pour un petit lot, l’approximation binomiale surestime la probabilité d’accepter un mauvais lot : en réalité, tirer sans remise fouille mieux le lot. Quant au thé, avec 4 tasses de chaque sorte, tout reconnaître par hasard n’arrive qu’une fois sur 70 : c’est assez rare pour croire au talent ; 3 bonnes réponses sur 4, en revanche, arrivent une fois sur quatre environ par simple hasard.
5
Signature — compter les poissons d’un lac sans les pêcher tous
Combien y a-t-il de poissons dans ce lac ? On en pêche K, on les marque et on les relâche. Quelques jours plus tard, ils se sont mélangés aux autres ; on pêche à nouveau n poissons et l’on compte k poissons marqués. Si le lac contient N poissons, k suit la loi hypergéométrique H(N, K, n). L’idée est d’égaler les proportions : la part de marqués dans la seconde pêche, k/n, doit ressembler à la part de marqués dans le lac, K/N. D’où l’estimation de Lincoln-Petersen N̂ = K·n / k. La même méthode sert à estimer des populations d’oiseaux, de baleines, de sans-abri ou d’erreurs dans un texte.
Le lac (la vraie taille N est cachée au pêcheur) : en rouge les poissons marqués lors de la première pêche ; le filet entoure la seconde pêche. À droite, la vraisemblance L(N) = P(k marqués | N) en fonction de la taille supposée N : son maximum est l’estimation de Lincoln-Petersen.
Marqués dans la seconde pêche—
Estimation de Lincoln-Petersen K·n/k—
Contrôle : maximum de la vraisemblance, cherché sur tous les N—doit valoir la partie entière de K·n/k
Mais une estimation se juge sur de nombreuses répétitions. Répétons l’expérience des milliers de fois sur le même lac : la formule de Lincoln-Petersen surestime N en moyenne, et elle échoue complètement quand k = 0. En 1951, Douglas Chapman propose une correction presque sans biais : N̂ = (K + 1)(n + 1)/(k + 1) − 1.
Distribution des estimations sur les expériences répétées : Lincoln-Petersen (rouge) et Chapman (bleu) ; le trait vert est la vraie taille N.
Lincoln-Petersen : moyenne (si k ≥ 1) exacte et simulée—
Chapman : moyenne exacte et simulée—
P(k = 0) : aucun marqué repêché—Lincoln-Petersen donne alors l’infini
Intervalle où tombent 95 % des estimations de Chapman—
Les moyennes « exactes » sont calculées en sommant l’estimation sur toutes les valeurs possibles de k, pondérées par la loi hypergéométrique : la simulation les retrouve. Le biais de Lincoln-Petersen vient de la division par k : quand k est petit, l’estimation explose, et ces cas tirent la moyenne vers le haut. L’estimateur de Chapman est exactement sans biais dès que K + n ≥ N, et reste presque sans biais au-delà. La précision dépend surtout du nombre de marqués repêchés : avec une dizaine seulement, l’estimation varie du simple au double.
Repères historiques
1657Christiaan Huygens pose, parmi les problèmes de son traité, celui de tirer des jetons d’une urne sans les remettre : un des premiers calculs hypergéométriques.
1713Jacques Bernoulli, dans l’Ars Conjectandi, étudie les tirages avec et sans remise ; Abraham de Moivre calcule en 1718 des temps d’attente qui relèvent de la loi géométrique.
1786Pierre-Simon Laplace estime la population de la France en rapportant le nombre de naissances de tout le pays à celui de communes recensées : l’ancêtre de la capture-recapture.
1896Le biologiste danois Carl Petersen marque des plies pour estimer les populations de poissons ; Frederick Lincoln applique l’idée aux canards en 1930.
1935Ronald Fisher décrit l’expérience de la dame qui goûte le thé dans The Design of Experiments et fonde le test exact sur la loi hypergéométrique.
1951Douglas Chapman propose l’estimateur corrigé qui porte son nom.
Ce qu’il faut retenir
Géométrique : attendre
P(X = k) = qᵏ⁻¹p, espérance 1/p, P(X > k) = qᵏ. La loi n’a pas de mémoire : un succès « en retard » n’est pas plus probable.
Hypergéométrique : tirer sans remettre
P(X = k) = C(K, k)C(N − K, n − k)/C(N, n), espérance nK/N. Elle est moins dispersée que la binomiale, d’un facteur (N − n)/(N − 1).
Des applications concrètes
Jeux de loterie, contrôle de petits lots, test exact de Fisher, capture-recapture : chaque fois qu’on tire dans un ensemble fini, la loi hypergéométrique donne les probabilités exactes.
Pour aller plus loin
Le pendant continu. La loi exponentielle, densité λe^(−λt), est la seule loi continue sans mémoire : c’est la limite de la loi géométrique quand les épreuves deviennent très rapprochées et très peu probables. Elle décrit la durée de vie d’un atome radioactif ou l’attente entre deux appels.
Attendre r succès. Le nombre d’épreuves nécessaires pour obtenir r succès est une somme de r variables géométriques indépendantes : il suit la loi binomiale négative, d’espérance r/p. Sans remise, le même problème donne la loi hypergéométrique négative, celle de l’as de pique de la planche 2.
De l’hypergéométrique à la binomiale. Quand N et K tendent vers l’infini avec K/N = p fixé, H(N, K, n) tend vers B(n, p). L’écart en variation totale est de l’ordre de n/N, ce que mesure la planche 3.