Lemme de Neyman-Pearson

lemme en statistiques

En statistique, selon le lemme de Neyman-Pearson, lorsque l'on veut effectuer un test d'hypothèse entre deux hypothèses H0 : θ = θ0 et H1 : θ = θ1, pour un échantillon 𝐱=(X1,…,Xn), alors le test du rapport de vraisemblance, qui rejette H0 en faveur de H1 lorsque ℒ(𝐱,θ0)ℒ(𝐱,θ1)≤kα, où kα est tel que

Lemme de Neyman-Pearson
Type
Inventeurs
Nommé en référence à
Formule
Λ(x)=L(x∣θ0)L(x∣θ1)≤ηVoir et modifier les données sur Wikidata
P(ℒ(𝐱,θ0)ℒ(𝐱,θ1)≤kα|H0)=α, est le test le plus puissant de niveau α.

Ce lemme est nommé d'après Jerzy Neyman et Egon Sharpe Pearson dans un article publié en 1933[1].

En pratique, la plupart du temps, le rapport de vraisemblance lui-même n'est pas explicitement utilisé dans le test. En effet, le test du rapport de vraisemblance ci-dessus est souvent équivalent à un test de la forme T≤tα pour une statistique T plus simple, et le test est effectué sous cette forme-ci.

Démonstration

modifier

Théorème : La région de rejet R0 optimale est définie par l'ensemble des points 𝐱=(x1,…,xn)∈ℝn tels que

ℒ(𝐱,θ0)ℒ(𝐱,θ1)≤kα

où la constante kα est telle que P(𝐱∈R0|θ0)=α. À noter qu'on a les relations suivantes :

P(𝐃n∈R0|θ0)=α=∫R0ℒ(𝐱;θ0) d𝐱
P(𝐃n∈R0|θ1)=1−β=∫R0ℒ(𝐱;θ1) d𝐱


où Dn=(x'1,…,x'n) est l'échantillon.

Démonstration :

Montrons tout d'abord que lorsque f𝒳(.;θ) est une densité bornée, il existe toujours une constante k telle que

P(ℒ(𝐱,θ0)ℒ(𝐱,θ1)>k|H0)=α.
En effet, lorsque k=0, cette probabilité vaut 1. D'autre part, cette probabilité décroit monotonément et continument vers zéro, lorsque k→∞. Par conséquent, il doit exister une valeur finie de k, appelée kα, qui satisfait l'égalité, ∀α∈]0;1[.
Désignons alors par R0, le sous-ensemble de ℝn suivant,
R0≜{𝐱∈ℝn|ℒ(𝐱,θ0)ℒ(𝐱,θ1)≤kα},
et soit R une autre partie de ℝn, telle que P(𝐱∈R|θ0)≤α.

Montrons que P(𝐱∈R0|θ1)>P(𝐱∈R|θ1):

P(𝐱∈R0|θ1)−P(𝐱∈R|θ1)=∫R0ℒ(𝐱;θ1) d𝐱−∫Rℒ(𝐱;θ1) d𝐱=∫R0∖Rℒ(𝐱;θ1) d𝐱−∫R∖R0ℒ(𝐱;θ1) d𝐱

Or ℒ(𝐱;θ1)≥kαℒ(𝐱;θ0) sur R0 et ℒ(𝐱;θ1)<kαℒ(𝐱;θ0) en dehors

P(𝐱∈R0|θ1)−P(𝐱∈R|θ1)≥kα(∫R0∖Rℒ(𝐱;θ0) d𝐱−∫R∖R0ℒ(𝐱;θ0) d𝐱)≥kα(∫R0ℒ(𝐱;θ0) d𝐱−∫Rℒ(𝐱;θ0) d𝐱)

La première intégrale vaut α par construction, la deuxième est majorée par α, on obtient:

P(𝐱∈R0|θ1)−P(𝐱∈R|θ1)≥0 ce qui conclut.

Notes et références

modifier
  1. ↑ (en) J. Neyman et E. S. Pearson, « IX. On the problem of the most efficient tests of statistical hypotheses », Phil. Trans. R. Soc. Lond. A, vol. 231, nos 694-706,‎ , p. 289–337 (ISSN 0264-3952, DOI 10.1098/rsta.1933.0009, lire en ligne)

Liens externes

modifier

Voir aussi

modifier