• 12 heures
  • Moyenne

Ce cours est visible gratuitement en ligne.

course.header.alt.is_video

course.header.alt.is_certifying

J'ai tout compris !

Mis à jour le 30/03/2020

Testez une moyenne ou une variance

Connectez-vous ou inscrivez-vous gratuitement pour bénéficier de toutes les fonctionnalités de ce cours !

Tester la moyenne (théorique) d'un échantillon

Le coin méthodologique

Supposons que nous ayons à disposition un échantillon i.i.d de loi N(μ,σ2), ou un grand échantillon i.i.d quelconque (en pratique de taille supérieure à 30).
On considère les tests suivants (c'est le test 3 qui est considéré dans le second exemple introductif) :

  1. {H0:μ=μ0H1:μ>μ0

  2. {H0:μ=μ0H1:μ<μ0

  3. {H0:μ=μ0H1:μμ0

On pourrait utiliser X¯μ0 (ou X¯) comme statistique de test, malheureusement sa loi dépend du paramètre σ inconnu, c'est pourquoi on utilise :

T=nX¯μ0S

On a sous H0 :

TT(n1)

Les régions critiques au niveau α sont :

  1.   W={T>tn1,1α}.

  2.   W={T<tn1,1α} .

  3.   W={|T|>tn1,1α2}.

On constate bien ici que la forme de la région critique est donnée par l'hypothèse alternative H1.

Les p-valeurs sont :

  1.  P(T(n1)nx¯μ0s).

  2.  P(T(n1)nx¯μ0s).

  3.  P(|T(n1)|nx¯μ0s)=2P(T(n1)nx¯μ0s).

Le coin R : exemple de la consommation d'essence

On teste :

{H0:μ=μ0H1:μμ0

avec  μ0=31

Pour α=5 , on a :

  •  tn1,1α2=t127,0.9751.97

  •  |n x¯μ0s||128 31.45312.16|2.35 .

On rejette H0 au niveau de test 5% car :

|n x¯μ0s|=2.35>1.97tn1,1α2 .

La p-valeur vaut :
 p-valeur=P(|T(n1)|n x¯μ0s)=2(1P(T(n1)n x¯μ0s))2(1P(T(127)2.35))0.02 .

On constate qu’on rejette bien l’hypothèse nulle au niveau de test 5% ( 0.02<0.05 ) mais pas au niveau de test 1% ( 0.020.05 )

En pratique, le data analyst peut utiliser la commande  t.test  pour résoudre ce test :

alpha <- 0.05
t.test(essence$conso,mu=31,alternative="two.sided")
## 
## Results of Hypothesis Test
## --------------------------
## 
## Null Hypothesis: mean = 31
## 
## Alternative Hypothesis: True mean is not equal to 31
## 
## Test Name: One Sample t-test
## 
## Estimated Parameter(s): mean of x = 31.44945
## 
## Data: essence$conso
## 
## Test Statistic: t = 2.354358
## 
## Test Statistic Parameter: df = 127
## 
## P-value: 0.02008833
## 
## 95% Confidence Interval: LCL = 31.07169
## UCL = 31.82722

Tester la variance (théorique) d'un échantillon

Supposons que nous ayons à disposition un échantillon i.i.d de loi N(μ,σ2), ou un grand échantillon i.i.d quelconque (en pratique de taille supérieure à 30).
On considère les tests suivants :

  1. {H0:σ2=σ02H1:σ2>σ02

  2. {H0:σ2=σ02H1:σ2<σ02

  3. {H0:σ2=σ02H1:σ2σ02

On utilise comme statistique de test :

K=(n1)S2σ02

On a sous H0 :

  Kχ2(n1)

Les régions critiques au niveau α sont :

  1.  W={K>χn1,1α2} .

  2.  W={K<χn1,α2} .

  3.  W={K<χn1,α22 ou K>χn1,1α22} .

De là on tire les p-valeurs...

Le coin R : exemple de la consommation d'essence

Considérons par exemple le test suivant :

{H0:σ2=σ02H1:σ2>σ02

avec σ02=4.5 .

En pratique, le data analyst pourra de nouveau utiliser la commande  varTest  issue du package  EnvStats  pour obtenir cet intervalle de confiance :

library(EnvStats)
alpha <- 0.05
varTest(essence$conso,sigma.squared=4.5,alternative="greater")

## 
## Results of Hypothesis Test
## --------------------------
## 
## Null Hypothesis: variance = 4.5
## 
## Alternative Hypothesis: True variance is greater than 4.5
## 
## Test Name: Chi-Squared Test on Variance
## 
## Estimated Parameter(s): variance = 4.66481
## 
## Data: essence$conso
## 
## Test Statistic: Chi-Squared = 131.6513
## 
## Test Statistic Parameter: df = 127
## 
## P-value: 0.3706697
## 
## 95% Confidence Interval: LCL = 3.839436
## UCL = Inf

On ne rejette donc pas l’hypothèse nulle au niveau de test 5%, la p-valeur vaut en effet environ 0.37 (elle n’est évidemment pas inférieure à 0.05).

Exemple de certificat de réussite
Exemple de certificat de réussite