Probabilités et statistique 201-SN1-RE • Cégep, Sciences de la nature, à Montréal

Exercices corrigés : tests d'hypothèse (201-SN1-RE)

Voici dix exercices corrigés sur les tests d'hypothèse du cours Probabilités et statistique 201-SN1-RE du cégep, en Sciences de la nature, précédés d'une liste de calculs.

Ce que couvre la série, les pièges traités

Ils couvrent l'hypothèse nulle et l'hypothèse alternative, le test unilatéral et le test bilatéral, le seuil de signification, la statistique de test, la région critique, la valeur pp, les erreurs de première et de seconde espèce et la puissance, sur les trois cas du devis : le test sur une moyenne avec n≥30n\geq 30, le test sur une moyenne avec n<30n<30 par la loi de Student, et le test sur une proportion avec n≥30n\geq 30, ainsi que le lien entre un test et un intervalle de confiance.

L'exercice 0.1 aligne onze calculs directs : statistiques zz et tt, erreur type sous H0H_{0}, valeurs critiques, valeurs pp et décisions. Les situations viennent ensuite des sciences de la nature : le rendement d'un engrais (exercice 2), le phosphore d'un lac (exercice 3), les grenouilles malformées d'un étang puis de douze étangs (exercices 4 et 10), la mesure de gg au laboratoire (exercices 5 et 7), le plomb de l'eau potable face à une norme (exercice 6), un catalyseur (exercice 8) et la durée de floraison d'une plante, menée comme une étude complète (exercice 9).

À la fin de la série, vous savez poser H0H_{0} et H1H_{1} sur le paramètre (exercice 1), choisir entre la loi normale et la loi de Student, rédiger une décision dans le contexte, démontrer qu'un test bilatéral et un intervalle de confiance décident de la même façon (exercice 5), calculer le risque de seconde espèce et la puissance (exercices 6 et 10), et retrouver une valeur de la table de Student par une intégrale (exercice 7). Les pièges traités : écrire les hypothèses sur x‾\overline{x} (exercice 1), écarter une mesure gênante (exercice 3), mettre p^\hat{p} au lieu de p0p_{0} dans l'erreur type (exercice 4), comparer à 1,961{,}96 avec un petit échantillon (exercice 7) ou dans un test unilatéral (exercice 8), lire la valeur pp comme la probabilité de H0H_{0} (exercice 8), et multiplier les tests sans compter les fausses alarmes (exercice 10).

Le fil de la série : tout se calcule SOUS H0H_{0}, et rien ne prouve H0H_{0}. On suppose l'hypothèse nulle vraie pour calculer, avec μ0\mu_{0} au centre et p0p_{0} dans l'erreur type ; seule une surprise plus grande que le seuil fixé d'avance la fait tomber, et ne pas la rejeter ne la démontre jamais. Faites chaque exercice au complet avant d'ouvrir la correction.

  • 11 exercices corrigés
  • Durée : 3 h 00
  • Autocorrigée, 122 réponses
  • Calculatrice
Commencer l'exercice 0.1 La fiche du chapitre

Rappel de cours

  • • HYPOTHÈSES : sur un PARAMÈTRE de la population (μ\mu ou pp), jamais sur x‾\overline{x} ou p^\hat{p}. H0H_{0} porte l'égalité (μ=μ0\mu=\mu_{0}) ; H1H_{1} porte ce que l'étude veut établir : μ≠μ0\mu\neq\mu_{0} (bilatéral), μ<μ0\mu<\mu_{0} ou μ>μ0\mu>\mu_{0} (unilatéral). Le choix se fait AVANT de voir les données.
  • • STATISTIQUE DE TEST, calculée sous H0H_{0} : moyenne avec n≥30n\geq 30, z=x‾−μ0s/nz=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}} (loi normale) ; moyenne avec n<30n<30 et population normale, t=x‾−μ0s/nt=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}} (Student à n−1n-1 degrés de liberté) ; proportion avec np0≥5np_{0}\geq 5 et n(1−p0)≥5n(1-p_{0})\geq 5, z=p^−p0p0(1−p0)/nz=\dfrac{\hat{p}-p_{0}}{\sqrt{p_{0}(1-p_{0})/n}}.
  • • RÉGION CRITIQUE au seuil α\alpha : bilatéral, ∣z∣>zα/2|z|>z_{\alpha/2} (1,961{,}96 à 5 %5\,\%) ; unilatéral, z>zαz>z_{\alpha} ou z<−zαz<-z_{\alpha} (1,6451{,}645 à 5 %5\,\%, 2,3262{,}326 à 1 %1\,\%). Avec Student, la valeur critique se lit à n−1n-1 degrés de liberté.
  • • VALEUR pp : la probabilité, SI H0H_{0} est vraie, d'obtenir un résultat au moins aussi éloigné de H0H_{0} que celui observé, dans le sens de H1H_{1} (les deux queues en bilatéral). On rejette H0H_{0} si p≤αp\leq\alpha.
  • • ERREURS : première espèce, rejeter H0H_{0} vraie, de probabilité α\alpha ; seconde espèce, ne pas rejeter H0H_{0} fausse, de probabilité β\beta, qui dépend de la vraie valeur du paramètre. Puissance =1−β=1-\beta. À nn fixé, baisser α\alpha augmente β\beta ; augmenter nn réduit β\beta.
  • • CONCLUSION : « au seuil de 5 %5\,\%, les données permettent (ou ne permettent pas) de conclure que... », dans le vocabulaire du problème. On ne dit jamais « on accepte H0H_{0} ». Un test bilatéral au seuil α\alpha rejette H0H_{0} exactement quand μ0\mu_{0} sort de l'intervalle de confiance de niveau 1−α1-\alpha.

Calculs : statistiques de test, valeurs p et décisions

Exercice 0.1 : Statistiques de test, valeurs p et décisions : onze calculs directs

12 réponses vérifiables

Pour chaque ligne, calculez la quantité demandée à partir des données résumées, sans rédiger de conclusion. Valeurs critiques de la loi normale : z0,05=1,645z_{0{,}05}=1{,}645, z0,025=1,96z_{0{,}025}=1{,}96, z0,01=2,326z_{0{,}01}=2{,}326, z0,005=2,576z_{0{,}005}=2{,}576, où zαz_{\alpha} laisse une aire α\alpha à sa droite. Les probabilités de la loi normale se lisent dans la table, avec zz arrondi au centième. Statistiques au centième, valeurs pp à quatre décimales.

  • a) H0:μ=50H_{0}:\mu=50, x‾=52,4\overline{x}=52{,}4, s=6s=6, n=36n=36 : calculez la statistique zz.
  • b) H0:μ=12H_{0}:\mu=12, x‾=11,1\overline{x}=11{,}1, s=2,7s=2{,}7, n=81n=81 : calculez la statistique zz.
  • c) H0:μ=80H_{0}:\mu=80, x‾=83\overline{x}=83, s=5s=5, n=9n=9, population normale : calculez la statistique tt et son nombre de degrés de liberté.
  • d) H0:p=0,40H_{0}:p=0{,}40, 230230 succès sur n=500n=500 : calculez la statistique zz.
  • e) H0:p=0,2H_{0}:p=0{,}2, n=400n=400 : calculez l'erreur type de P^\hat{P} sous H0H_{0}.
  • f) Test unilatéral à gauche au seuil de 2,5 %2{,}5\,\% : donnez la valeur critique de zz.
  • g) Test bilatéral, z=−2,15z=-2{,}15 : calculez la valeur pp.
  • h) Test unilatéral à droite, z=1,32z=1{,}32 : calculez la valeur pp.
  • i) Test bilatéral au seuil de 5 %5\,\%, z=2,05z=2{,}05 : donnez la décision.
  • j) H0:μ=200H_{0}:\mu=200 contre H1:μ>200H_{1}:\mu>200, σ=15\sigma=15, n=25n=25, seuil de 5 %5\,\% : calculez la borne cc de la région critique x‾>c\overline{x}>c.
  • k) Test bilatéral au seuil de 5 %5\,\% sur une moyenne, population normale, n=15n=15, t=−2,40t=-2{,}40, et t0,025 ; 14=2,145t_{0{,}025\,;\,14}=2{,}145 : donnez la décision.

Tape tes réponses, la page te dit juste ou faux 0/12

a)
b)
c)
d)
e)
f)
g)
h)
i)
j)
k)
Voir la correction

Réponses

  • a) z=2,40z=2{,}40
  • b) z=−3,00z=-3{,}00
  • c) t=1,80t=1{,}80 à 88 degrés de liberté
  • d) z≈2,74z\approx 2{,}74
  • e) σP^=0,02\sigma_{\hat{P}}=0{,}02
  • f) z=−1,96z=-1{,}96
  • g) p≈0,0316p\approx 0{,}0316
  • h) p≈0,0934p\approx 0{,}0934
  • i) Rejeter H0H_{0}
  • j) c≈204,94c\approx 204{,}94
  • k) Rejeter H0H_{0}
a)

n≥30n\geq 30 : on remplace σ\sigma par ss et on utilise la loi normale.

z=x‾−μ0s/n=52,4−506/6=2,40\displaystyle z=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}}=\dfrac{52{,}4-50}{6/6}=2{,}40
b)
z=11,1−122,7/81=−0,90,3=−3,00\displaystyle z=\dfrac{11{,}1-12}{2{,}7/\sqrt{81}}=\dfrac{-0{,}9}{0{,}3}=-3{,}00

Le signe dit de quel côté de μ0\mu_{0} tombe l'échantillon ; on le garde.

c)

n<30n<30 et population normale : loi de Student à n−1=8n-1=8 degrés de liberté.

t=83−805/9=35/3=1,80\displaystyle t=\dfrac{83-80}{5/\sqrt{9}}=\dfrac{3}{5/3}=1{,}80
d)

p^=230500=0,46\hat{p}=\dfrac{230}{500}=0{,}46 et, sous H0H_{0}, σP^=0,40×0,60500≈0,02191\sigma_{\hat{P}}=\sqrt{\dfrac{0{,}40\times 0{,}60}{500}}\approx 0{,}02191.

z=0,46−0,400,02191≈2,74\displaystyle z=\dfrac{0{,}46-0{,}40}{0{,}02191}\approx 2{,}74
e)
σP^=0,2×0,8400=0,0004=0,02\displaystyle \sigma_{\hat{P}}=\sqrt{\dfrac{0{,}2\times 0{,}8}{400}}=\sqrt{0{,}0004}=0{,}02
f)

Toute l'aire α=0,025\alpha=0{,}025 est dans la queue gauche : la valeur critique laisse 2,5 %2{,}5\,\% à sa gauche, c'est −z0,025=−1,96-z_{0{,}025}=-1{,}96.

Piège : prendre −1,645-1{,}645 parce que le test est unilatéral. 1,6451{,}645 correspond à 5 %5\,\% dans une queue, pas à 2,5 %2{,}5\,\%.

g)

P(Z<−2,15)=1−0,9842=0,0158P(Z<-2{,}15)=1-0{,}9842=0{,}0158, et le test est bilatéral :

p=2×0,0158=0,0316\displaystyle p=2\times 0{,}0158=0{,}0316
h)
p=P(Z>1,32)=1−0,9066=0,0934\displaystyle p=P(Z>1{,}32)=1-0{,}9066=0{,}0934
i)

∣z∣=2,05>1,96|z|=2{,}05>1{,}96 : zz tombe dans la région critique, on rejette H0H_{0} au seuil de 5 %5\,\%.

j)

On rejette quand x‾−20015/25>1,645\dfrac{\overline{x}-200}{15/\sqrt{25}}>1{,}645, c'est-à-dire quand

x‾>200+1,645×155=204,935\displaystyle \overline{x}>200+1{,}645\times\dfrac{15}{5}=204{,}935

Donc c≈204,94c\approx 204{,}94.

k)

Loi de Student à 1414 degrés de liberté : on rejette quand ∣t∣>2,145|t|>2{,}145. Ici ∣t∣=2,40>2,145|t|=2{,}40>2{,}145 : on rejette H0H_{0}.

Coche ici les exercices faits ou à revoir : un compte gratuit, sans mot de passe, retient tes coches et tes réponses justes d'une visite à l'autre, te dit quel chapitre attaquer ensuite et te permet de demander l'exercice qui te manque. Crée ton espace, un courriel suffit.

Partie A : les bases (/50)

Exercice 1 : Poser H0 et H1 : sept études, sept paires d'hypothèses

23 réponses vérifiables

Un test d'hypothèse commence avant tout calcul : on nomme le paramètre de la POPULATION, on écrit l'hypothèse nulle H0H_{0}, qui porte l'égalité, puis l'hypothèse alternative H1H_{1}, celle que l'étude cherche à établir. Ces sept situations reviennent dans toute la série.

Pour chacune des questions a) à g), donnez le paramètre testé, la forme de H1H_{1} et la valeur de référence, puis écrivez H0H_{0} et H1H_{1}.

  • a) Engrais. Sans engrais, le rendement moyen d'une variété de blé est de 3,203{,}20 t/ha sur les parcelles de la station. Un agronome veut établir que son engrais AUGMENTE le rendement moyen.
  • b) Plomb. Une municipalité ne déclare l'eau d'un réseau conforme que si les données établissent que la concentration moyenne en plomb au robinet est INFÉRIEURE à une norme de 88 µg/L (norme fictive, fixée pour cette série).
  • c) Phosphore. L'association d'un lac a fixé un seuil de gestion de 2020 µg/L de phosphore total (valeur choisie pour l'exercice). Elle veut savoir si la concentration moyenne de l'été DÉPASSE ce seuil.
  • d) Mesure de gg. Une classe veut savoir si son montage de pendule donne une valeur de gg biaisée, dans un sens ou dans l'autre, par rapport à la valeur de référence 9,819{,}81 m/s².
  • e) Catalyseur. Sans catalyseur, une réaction dure en moyenne 120120 s dans les conditions du laboratoire. Une chimiste veut établir qu'un catalyseur RÉDUIT cette durée moyenne.
  • f) Grenouilles. Dans des étangs témoins, 5 %5\,\% des grenouilles présentent une malformation des membres (taux fictif). Un biologiste soupçonne un taux PLUS ÉLEVÉ dans un étang voisin de cultures.
  • g) Floraison. Des relevés anciens donnent une durée moyenne de floraison de 18,018{,}0 jours pour une plante de sous-bois. Une botaniste veut établir que cette durée a DIMINUÉ.
  • h) Pour l'engrais, un étudiant écrit « H0:x‾=3,20H_{0}:\overline{x}=3{,}20 et H1:x‾>3,20H_{1}:\overline{x}>3{,}20 ». Expliquez son erreur.
  • i) Pour le plomb, un conseiller municipal propose plutôt H0:μ=8H_{0}:\mu=8 contre H1:μ>8H_{1}:\mu>8, « pour ne déclarer l'eau contaminée que si on le prouve ». Expliquez ce que ce choix change quand l'échantillon est petit ou peu concluant.

Tape tes réponses, la page te dit juste ou faux 0/23

a)
b)
c)
d)
e)
f)
g)
h)
i)
Voir la correction

Réponses

  • a) H0:μ=3,20H_{0}:\mu=3{,}20 ; H1:μ>3,20H_{1}:\mu>3{,}20 (unilatéral à droite)
  • b) H0:μ=8H_{0}:\mu=8 ; H1:μ<8H_{1}:\mu<8 (unilatéral à gauche)
  • c) H0:μ=20H_{0}:\mu=20 ; H1:μ>20H_{1}:\mu>20 (unilatéral à droite)
  • d) H0:μ=9,81H_{0}:\mu=9{,}81 ; H1:μ≠9,81H_{1}:\mu\neq9{,}81 (bilatéral)
  • e) H0:μ=120H_{0}:\mu=120 ; H1:μ<120H_{1}:\mu<120 (unilatéral à gauche)
  • f) H0:p=0,05H_{0}:p=0{,}05 ; H1:p>0,05H_{1}:p>0{,}05 (unilatéral à droite)
  • g) H0:μ=18H_{0}:\mu=18 ; H1:μ<18H_{1}:\mu<18 (unilatéral à gauche)
  • h) Les hypothèses portent sur le paramètre μ\mu, jamais sur x‾\overline{x}, qui est connu
  • i) L'eau resterait en service par défaut, sans que sa conformité soit prouvée
a)

Paramètre μ\mu, rendement moyen de toutes les parcelles traitées.

H0:μ=3,20H_{0}:\mu=3{,}20 et H1:μ>3,20H_{1}:\mu>3{,}20.

Test unilatéral à droite, car seule une HAUSSE du rendement intéresse l'agronome ; une baisse ne ferait pas adopter l'engrais non plus, mais ce n'est pas ce qu'il cherche à établir.

b)

Paramètre μ\mu, concentration moyenne en plomb au robinet du réseau.

H0:μ=8H_{0}:\mu=8 et H1:μ<8H_{1}:\mu<8.

Test unilatéral à gauche, car c'est la conformité qui doit être PROUVÉE : on ne déclare l'eau sûre que si les données placent la moyenne sous la norme.

c)

Paramètre μ\mu, concentration moyenne de phosphore total du lac cet été.

H0:μ=20H_{0}:\mu=20 et H1:μ>20H_{1}:\mu>20.

Test unilatéral à droite, car l'association ne s'inquiète que d'un dépassement : une concentration plus basse que le seuil ne demande aucune action.

d)

Paramètre μ\mu, moyenne de toutes les mesures que donnerait le montage.

H0:μ=9,81H_{0}:\mu=9{,}81 et H1:μ≠9,81H_{1}:\mu\neq9{,}81.

Test bilatéral, car un biais dans un sens ou dans l'autre est un défaut du montage : on ne sait pas d'avance de quel côté il tombe.

e)

Paramètre μ\mu, durée moyenne de la réaction avec catalyseur.

H0:μ=120H_{0}:\mu=120 et H1:μ<120H_{1}:\mu<120.

Test unilatéral à gauche, car le catalyseur doit RÉDUIRE la durée ; une durée plus longue ne le rendrait pas utile.

f)

Paramètre pp, proportion de grenouilles malformées dans l'étang.

H0:p=0,05H_{0}:p=0{,}05 et H1:p>0,05H_{1}:p>0{,}05.

Test unilatéral à droite, car le soupçon porte sur un taux plus élevé que celui des étangs témoins.

g)

Paramètre μ\mu, durée moyenne de floraison actuelle.

H0:μ=18H_{0}:\mu=18 et H1:μ<18H_{1}:\mu<18.

Test unilatéral à gauche, car la question est une diminution ; une durée allongée serait une autre étude.

h)

Les hypothèses portent sur un PARAMÈTRE de la population, inconnu : ici μ\mu, le rendement moyen de toutes les parcelles qu'on traiterait avec l'engrais.

x‾\overline{x} est une statistique de l'échantillon : une fois la récolte pesée, sa valeur est connue, et une hypothèse sur elle n'a pas de sens.

Il faut écrire H0:μ=3,20H_{0}:\mu=3{,}20 et H1:μ>3,20H_{1}:\mu>3{,}20. On écrit l'égalité dans H0H_{0} même si l'on pense « μ≤3,20\mu\leq 3{,}20 » : c'est la valeur sous laquelle on fera tous les calculs, le cas limite entre les deux hypothèses.

Piège : écrire l'affirmation à établir dans H0H_{0}. Un test ne peut que rejeter H0H_{0} ou ne pas la rejeter ; ce qu'on veut établir va donc dans H1H_{1}.

i)

Un test ne rejette H0H_{0} que si les données sont nettement incompatibles avec elle. Avec H1:μ>8H_{1}:\mu>8, un petit échantillon peu concluant mène à ne pas rejeter μ=8\mu=8 : on ne prouve pas la contamination, et l'eau reste en service.

L'eau serait donc déclarée acceptable PAR DÉFAUT, simplement faute de données. Avec H1:μ<8H_{1}:\mu<8, c'est l'inverse : tant que la conformité n'est pas établie, l'eau n'est pas déclarée conforme.

Le choix de H1H_{1} décide qui porte la charge de la preuve. En santé publique, c'est à la conformité d'être prouvée, et l'exercice 6 chiffre les deux erreurs que ce choix répartit.

0−1,961,962,5 %2,5 %bilatéralH₁ : μ ≠ μ₀0−1,6455 %unilatéral à gaucheH₁ : μ < μ₀01,6455 %unilatéral à droiteH₁ : μ > μ₀

Exercice 2 : Test sur une moyenne avec n ≥ 30 : l'engrais augmente-t-il le rendement ?

9 réponses vérifiables

On reprend l'engrais de l'exercice 1 a). Sans engrais, le rendement moyen de la variété est de μ0=3,20\mu_{0}=3{,}20 t/ha. L'engrais est épandu sur n=40n=40 parcelles tirées au hasard dans la station. Sur les données brutes, le rendement moyen vaut x‾=3,38\overline{x}=3{,}38 t/ha et l'écart-type s=0,52s=0{,}52 t/ha. On travaille au seuil α=5 %\alpha=5\,\% ; on donne z0,05=1,645z_{0{,}05}=1{,}645 et z0,01=2,326z_{0{,}01}=2{,}326.

L'histogramme ci-dessous regroupe les 4040 rendements en classes.

24681012142,42,83,23,644,44,8rendement (t/ha)effectif
  • a) Formulez H0H_{0} et H1H_{1} en nommant le paramètre.
  • b) L'histogramme est nettement asymétrique. Expliquez pourquoi on peut tout de même utiliser la loi normale pour X‾\overline{X}, et pourquoi on a le droit de remplacer σ\sigma par ss.
  • c) Calculez l'erreur type de X‾\overline{X} et la statistique zz.
  • d) Déterminez la région critique au seuil de 5 %5\,\%, d'abord en zz, puis en rendement moyen x‾\overline{x}.
  • e) Prenez la décision et rédigez la conclusion dans le contexte.
  • f) Calculez la valeur pp du test et vérifiez qu'elle conduit à la même décision.
  • g) Reprenez la décision au seuil de 1 %1\,\%, en calculant la nouvelle borne de la région critique, et rédigez la conclusion.

Tape tes réponses, la page te dit juste ou faux 0/9

a)
b)
c)
d)
e)
f)
g)
Voir la correction

Réponses

  • a) H0:μ=3,20H_{0}:\mu=3{,}20 ; H1:μ>3,20H_{1}:\mu>3{,}20, où μ\mu est le rendement moyen avec engrais
  • b) Théorème central limite (n=40≥30n=40\geq 30) ; ss estime σ\sigma avec précision quand n≥30n\geq 30
  • c) σX‾≈0,0822\sigma_{\overline{X}}\approx 0{,}0822 t/ha ; z≈2,19z\approx 2{,}19
  • d) z>1,645z>1{,}645, soit x‾>3,335\overline{x}>3{,}335 t/ha
  • e) Rejeter H0H_{0} : au seuil de 5 %5\,\%, les données appuient une hausse du rendement moyen
  • f) p≈0,0143<0,05p\approx 0{,}0143<0{,}05 : même décision
  • g) c≈3,391c\approx 3{,}391 ; ne pas rejeter H0H_{0} au seuil de 1 %1\,\%
a)

Le paramètre est μ\mu, le rendement moyen de TOUTES les parcelles de la station qu'on traiterait avec l'engrais.

H0:μ=3,20H_{0}:\mu=3{,}20 t/ha (l'engrais ne change rien) et H1:μ>3,20H_{1}:\mu>3{,}20 t/ha (il augmente le rendement moyen). Test unilatéral à droite.

b)

Le théorème central limite porte sur la MOYENNE : pour n≥30n\geq 30 tirages indépendants, X‾\overline{X} suit approximativement une loi normale de moyenne μ\mu et d'écart-type σn\dfrac{\sigma}{\sqrt{n}}, quelle que soit la forme de la distribution des rendements, pourvu que σ\sigma soit fini.

L'asymétrie de l'histogramme décrit les DONNÉES ; elle ne dit rien d'interdit sur la loi de X‾\overline{X}. Avec n=40n=40, l'approximation normale est raisonnable, d'autant que l'asymétrie reste modérée.

Pour la même raison, avec n≥30n\geq 30, ss estime σ\sigma avec assez de précision pour qu'on garde la loi normale : c'est la convention du cours pour une moyenne avec n≥30n\geq 30.

Piège : conclure qu'un test est impossible parce que les données ne sont pas normales. Cette condition ne sert que pour n<30n<30 (exercice 3).

c)

σX‾≈sn=0,5240≈0,0822\sigma_{\overline{X}}\approx\dfrac{s}{\sqrt{n}}=\dfrac{0{,}52}{\sqrt{40}}\approx 0{,}0822 t/ha

z=x‾−μ0s/n=3,38−3,200,0822≈2,19\displaystyle z=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}}=\dfrac{3{,}38-3{,}20}{0{,}0822}\approx 2{,}19

Le calcul se fait SOUS H0H_{0} : la loi de X‾\overline{X} est centrée sur μ0=3,20\mu_{0}=3{,}20, et c'est l'écart à cette valeur qu'on mesure en erreurs types.

d)

Test unilatéral à droite au seuil de 5 %5\,\% : on rejette H0H_{0} quand z>1,645z>1{,}645.

En rendement : x‾−3,200,0822>1,645  ⟺  x‾>3,20+1,645×0,0822≈3,335\dfrac{\overline{x}-3{,}20}{0{,}0822}>1{,}645 \iff \overline{x}>3{,}20+1{,}645\times 0{,}0822\approx 3{,}335 t/ha.

La région critique se fixe AVANT de regarder x‾\overline{x} : elle ne dépend que de μ0\mu_{0}, de l'erreur type et du seuil.

e)

z≈2,19>1,645z\approx 2{,}19>1{,}645, ou encore x‾=3,38>3,335\overline{x}=3{,}38>3{,}335 : on rejette H0H_{0}.

Conclusion : au seuil de 5 %5\,\%, les données des 4040 parcelles permettent de conclure que l'engrais augmente le rendement moyen de cette variété dans cette station.

La conclusion parle du rendement moyen et de la station étudiée ; elle ne dit rien d'une autre variété ni d'un autre sol.

f)
p=P(Z>2,19)=1−0,9857=0,0143\displaystyle p=P(Z>2{,}19)=1-0{,}9857=0{,}0143

Si l'engrais était sans effet, un rendement moyen d'au moins 3,383{,}38 t/ha sur 4040 parcelles n'apparaîtrait que dans environ 1,4 %1{,}4\,\% des échantillons.

p=0,0143<0,05p=0{,}0143<0{,}05 : on rejette H0H_{0}, comme en e). Les deux règles, région critique et valeur pp, donnent toujours la même décision.

g)

Au seuil de 1 %1\,\%, on rejette quand z>2,326z>2{,}326, soit x‾>3,20+2,326×0,0822≈3,391\overline{x}>3{,}20+2{,}326\times 0{,}0822\approx 3{,}391 t/ha.

x‾=3,38<3,391\overline{x}=3{,}38<3{,}391 (ou p=0,0143>0,01p=0{,}0143>0{,}01) : on ne rejette pas H0H_{0}.

Conclusion : au seuil de 1 %1\,\%, les données ne permettent pas de conclure à une hausse du rendement moyen.

Piège : écrire « au seuil de 1 %1\,\%, l'engrais est sans effet ». Ne pas rejeter H0H_{0} ne la démontre pas : les mêmes données montrent un effet probable, simplement pas assez net pour un seuil aussi exigeant.

3,203,335x̄ = 3,38non-rejetrejet : 5 %X̄ sous H₀

Exercice 3 : Test de Student, n < 30 : le phosphore du lac dépasse-t-il le seuil ?

8 réponses vérifiables

On reprend le lac de l'exercice 1 c). L'association a prélevé, cet été, n=10n=10 échantillons d'eau en des points et à des dates tirés au hasard. Concentrations de phosphore total, en µg/L : 22,4;19,6;23,4;21,1;17,2;24,1;20,3;22,7;18,8;21,422{,}4 ; 19{,}6 ; 23{,}4 ; 21{,}1 ; 17{,}2 ; 24{,}1 ; 20{,}3 ; 22{,}7 ; 18{,}8 ; 21{,}4.

On suppose que les concentrations de l'été suivent une loi normale. On teste H0:μ=20H_{0}:\mu=20 contre H1:μ>20H_{1}:\mu>20 au seuil de 5 %5\,\%. Valeurs de la loi de Student, où tα ; νt_{\alpha\,;\,\nu} laisse une aire α\alpha à sa droite : t0,10 ; 9=1,383t_{0{,}10\,;\,9}=1{,}383, t0,05 ; 9=1,833t_{0{,}05\,;\,9}=1{,}833, t0,025 ; 9=2,262t_{0{,}025\,;\,9}=2{,}262 et t0,05 ; 8=1,860t_{0{,}05\,;\,8}=1{,}860.

Le diagramme en points ci-dessous place les dix mesures par rapport au seuil.

1617181920212223242526seuil : 20phosphore total (µg/L)
  • a) Calculez la moyenne x‾\overline{x} et l'écart-type ss de l'échantillon.
  • b) Expliquez pourquoi on utilise la loi de Student et non la loi normale, quelle condition le permet, et comment le diagramme la rend plausible. Donnez le nombre de degrés de liberté.
  • c) Calculez la statistique tt.
  • d) Prenez la décision au seuil de 5 %5\,\% et rédigez la conclusion.
  • e) Encadrez la valeur pp du test à l'aide des valeurs de la table.
  • f) Un membre de l'association propose d'écarter la plus faible mesure, 17,217{,}2 µg/L, « sûrement une erreur ». Recalculez tt sans elle, puis expliquez pourquoi ce second calcul n'a pas le droit de remplacer le premier.

Tape tes réponses, la page te dit juste ou faux 0/8

a)
b)
c)
d)
e)
f)
Voir la correction

Réponses

  • a) x‾=21,10\overline{x}=21{,}10 µg/L ; s≈2,16s\approx 2{,}16 µg/L
  • b) n<30n<30 et σ\sigma inconnu ; population normale supposée ; 99 degrés de liberté
  • c) t≈1,61t\approx 1{,}61
  • d) Ne pas rejeter H0H_{0} : les données ne permettent pas de conclure que la moyenne dépasse 2020 µg/L
  • e) 0,05<p<0,100{,}05<p<0{,}10
  • f) t≈2,60t\approx 2{,}60, mais retirer une mesure sans raison expérimentale fausse le test
a)

x‾=211,010=21,10\overline{x}=\dfrac{211{,}0}{10}=21{,}10 µg/L.

s=∑(xi−x‾)2n−1=42,029≈2,16s=\sqrt{\dfrac{\sum (x_{i}-\overline{x})^{2}}{n-1}}=\sqrt{\dfrac{42{,}02}{9}}\approx 2{,}16 µg/L

À la calculatrice : mode statistique, liste des dix valeurs, et lire sxs_{x} (et non σx\sigma_{x}, qui divise par nn).

b)

Avec n=10<30n=10<30, le théorème central limite ne garantit plus que X‾\overline{X} soit normale, et ss, calculé sur dix valeurs, est une estimation incertaine de σ\sigma.

Si la population des concentrations est NORMALE, la statistique T=X‾−μ0S/nT=\dfrac{\overline{X}-\mu_{0}}{S/\sqrt{n}} suit exactement une loi de Student à n−1=9n-1=9 degrés de liberté : c'est la condition qui remplace n≥30n\geq 30.

Le diagramme ne prouve pas la normalité, mais il la rend plausible : une seule bosse, pas de forte asymétrie, aucune valeur isolée loin des autres.

Piège : utiliser z0,05=1,645z_{0{,}05}=1{,}645. Les queues de la loi de Student sont plus lourdes ; sa valeur critique, 1,8331{,}833, est plus grande.

c)
t=x‾−μ0s/n=21,10−202,161/10=1,100,683≈1,61\displaystyle t=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}}=\dfrac{21{,}10-20}{2{,}161/\sqrt{10}}=\dfrac{1{,}10}{0{,}683}\approx 1{,}61
d)

Test unilatéral à droite, 99 degrés de liberté : on rejette H0H_{0} si t>1,833t>1{,}833.

t≈1,61<1,833t\approx 1{,}61<1{,}833 : on ne rejette pas H0H_{0}.

Conclusion : au seuil de 5 %5\,\%, ces dix mesures ne permettent pas de conclure que la concentration moyenne de phosphore de l'été dépasse 2020 µg/L.

Piège : écrire « le lac respecte le seuil ». La moyenne observée, 21,1021{,}10, est au-dessus de 2020 ; le test dit seulement que cet écart peut venir du hasard de dix prélèvements.

e)

La valeur pp est l'aire à droite de t=1,61t=1{,}61 sous la loi de Student à 99 degrés de liberté.

1,383<1,61<1,8331{,}383<1{,}61<1{,}833, donc 0,05<p<0,100{,}05<p<0{,}10.

Une table de Student ne donne qu'un encadrement ; un logiciel donnerait p≈0,071p\approx 0{,}071, cohérent avec l'encadrement.

f)

Sans 17,217{,}2 : x‾≈21,53\overline{x}\approx 21{,}53, s≈1,77s\approx 1{,}77, n=9n=9, donc 88 degrés de liberté et

t=21,53−201,772/9≈2,60>1,860\displaystyle t=\dfrac{21{,}53-20}{1{,}772/\sqrt{9}}\approx 2{,}60>1{,}860

Ce second test rejetterait H0H_{0}, mais il n'a aucune valeur : la mesure a été retirée APRÈS avoir vu qu'elle gênait la conclusion souhaitée. On a choisi les données pour obtenir le résultat.

Une mesure ne s'écarte que pour une raison expérimentale documentée (flacon contaminé, erreur de protocole), décidée sans regarder l'effet sur le test. Sinon, la décision reste celle de d).

Exercice 4 : Test sur une proportion : les grenouilles malformées de l'étang

11 réponses vérifiables

On reprend l'étang de l'exercice 1 f). Dans les étangs témoins de la région, la proportion de grenouilles présentant une malformation des membres est de p0=0,05p_{0}=0{,}05 (taux fictif). Dans un étang voisin de cultures, une équipe capture au hasard n=220n=220 grenouilles et en trouve 1919 malformées. On teste au seuil de 1 %1\,\%, avec z0,01=2,326z_{0{,}01}=2{,}326.

  • a) Formulez H0H_{0} et H1H_{1}, puis calculez la proportion observée p^\hat{p}.
  • b) Vérifiez les conditions d'application de l'approximation normale.
  • c) Calculez l'erreur type de P^\hat{P} sous H0H_{0}, puis la statistique zz.
  • d) Prenez la décision au seuil de 1 %1\,\% et rédigez la conclusion.
  • e) Calculez la valeur pp du test.
  • f) Un étudiant calcule l'erreur type avec p^\hat{p} au lieu de p0p_{0}. Calculez la statistique qu'il obtient et la décision qu'il prend au seuil de 1 %1\,\%. Expliquez pourquoi c'est p0p_{0} qui doit figurer dans l'erreur type.
  • g) Dans son rapport, le biologiste écrit : « les pesticides des cultures voisines causent les malformations ». Dites ce que le test établit, et ce qu'il n'établit pas.

Tape tes réponses, la page te dit juste ou faux 0/11

a)
b)
c)
d)
e)
f)
g)
Voir la correction

Réponses

  • a) H0:p=0,05H_{0}:p=0{,}05 ; H1:p>0,05H_{1}:p>0{,}05 ; p^=19220≈0,0864\hat{p}=\frac{19}{220}\approx 0{,}0864
  • b) np0=11≥5np_{0}=11\geq 5 et n(1−p0)=209≥5n(1-p_{0})=209\geq 5 ; échantillon tiré au hasard
  • c) σP^≈0,0147\sigma_{\hat{P}}\approx 0{,}0147 ; z≈2,47z\approx 2{,}47
  • d) Rejeter H0H_{0} : au seuil de 1 %1\,\%, la proportion de malformées dépasse 5 %5\,\% dans cet étang
  • e) p≈0,0068p\approx 0{,}0068
  • f) z≈1,92z\approx 1{,}92 : il ne rejetterait pas H0H_{0} ; or le calcul se fait sous H0H_{0}, donc avec p0p_{0}
  • g) Une proportion supérieure à 5 %5\,\% dans cet étang ; pas la cause, l'étude est observationnelle
a)

Paramètre : pp, proportion de grenouilles malformées dans tout l'étang.

H0:p=0,05H_{0}:p=0{,}05 et H1:p>0,05H_{1}:p>0{,}05 : unilatéral à droite.

p^=19220≈0,0864\displaystyle \hat{p}=\dfrac{19}{220}\approx 0{,}0864
b)

np0=220×0,05=11≥5np_{0}=220\times 0{,}05=11\geq 5 et n(1−p0)=220×0,95=209≥5n(1-p_{0})=220\times 0{,}95=209\geq 5 : sous H0H_{0}, P^\hat{P} suit approximativement une loi normale.

Les conditions se vérifient avec p0p_{0}, la valeur sous laquelle on calcule. Il faut aussi des captures au hasard et indépendantes : si l'équipe n'avait capturé que les grenouilles lentes, plus souvent malformées, la proportion serait biaisée, et aucun calcul ne le corrigerait.

c)
σP^=p0(1−p0)n=0,05×0,95220≈0,0147\displaystyle \sigma_{\hat{P}}=\sqrt{\dfrac{p_{0}(1-p_{0})}{n}}=\sqrt{\dfrac{0{,}05\times 0{,}95}{220}}\approx 0{,}0147
z=p^−p0σP^=0,0864−0,050,0147≈2,47\displaystyle z=\dfrac{\hat{p}-p_{0}}{\sigma_{\hat{P}}}=\dfrac{0{,}0864-0{,}05}{0{,}0147}\approx 2{,}47
d)

On rejette H0H_{0} si z>2,326z>2{,}326. Ici z≈2,47>2,326z\approx 2{,}47>2{,}326 : on rejette H0H_{0}.

Conclusion : au seuil de 1 %1\,\%, les données permettent de conclure que la proportion de grenouilles malformées dans cet étang dépasse le taux de 5 %5\,\% des étangs témoins.

e)
p=P(Z>2,47)=1−0,9932=0,0068\displaystyle p=P(Z>2{,}47)=1-0{,}9932=0{,}0068

Si le taux de l'étang était de 5 %5\,\%, on ne trouverait 1919 malformées ou plus sur 220220 que dans environ 0,7 %0{,}7\,\% des échantillons. p<0,01p<0{,}01 : même décision qu'en d).

f)

Avec p^\hat{p} : 0,0864×0,9136220≈0,0189\sqrt{\dfrac{0{,}0864\times 0{,}9136}{220}}\approx 0{,}0189 et z≈0,03640,0189≈1,92<2,326z\approx\dfrac{0{,}0364}{0{,}0189}\approx 1{,}92<2{,}326 : il ne rejetterait pas H0H_{0}.

Or tout le test raisonne SOUS H0H_{0} : on calcule la loi de P^\hat{P} en supposant p=0,05p=0{,}05, et dans ce monde l'erreur type vaut 0,05×0,95220\sqrt{\frac{0{,}05\times 0{,}95}{220}}. Mettre p^\hat{p}, c'est mélanger deux hypothèses dans un même calcul.

Piège : la formule de l'intervalle de confiance, qui utilise p^\hat{p} parce qu'il n'y a pas de valeur supposée. Ici l'erreur coûte toute la conclusion, comme le montre la figure.

Plus p^\hat{p} s'éloigne de p0p_{0}, plus les deux erreurs types diffèrent, et plus ce choix pèse sur la décision.

g)

Le test établit une seule chose : dans cet étang, la proportion de malformées dépasse 5 %5\,\%, avec un risque de première espèce de 1 %1\,\%.

Il ne dit rien de la CAUSE : l'étude est observationnelle. L'étang diffère des étangs témoins par les cultures voisines, mais peut-être aussi par sa température, ses parasites (certains trématodes provoquent des malformations des membres) ou sa faune.

Pour attribuer la hausse aux pesticides, il faudrait mesurer leur concentration dans plusieurs étangs, ou une expérience contrôlée en laboratoire. Le rapport doit écrire « est associée à la proximité des cultures, sans que la cause soit établie ».

012,3264avec p̂ : 1,92avec p₀ : 2,47rejet : 1 %

Exercice 5 : Test bilatéral et intervalle de confiance : la mesure de g est-elle biaisée ?

11 réponses vérifiables

On reprend la classe de l'exercice 1 d). Chaque équipe mesure la période d'un pendule et en déduit une valeur de gg par g=4π2LT2g=\dfrac{4\pi^{2}L}{T^{2}}. Sur n=50n=50 mesures, la moyenne vaut x‾=9,774\overline{x}=9{,}774 m/s² et l'écart-type s=0,092s=0{,}092 m/s². La valeur de référence au laboratoire est g0=9,81g_{0}=9{,}81 m/s².

On donne z0,025=1,96z_{0{,}025}=1{,}96 et z0,005=2,576z_{0{,}005}=2{,}576.

  • a) Formulez H0H_{0} et H1H_{1}, et justifiez le choix d'un test bilatéral.
  • b) Calculez l'erreur type et la statistique zz, puis prenez la décision au seuil de 5 %5\,\%.
  • c) Calculez la valeur pp du test.
  • d) Construisez l'intervalle de confiance à 95 %95\,\% pour μ\mu, au millième. Contient-il 9,819{,}81 ?
  • e) Démontrez que, pour un test bilatéral au seuil α\alpha, on rejette H0:μ=μ0H_{0}:\mu=\mu_{0} si et seulement si μ0\mu_{0} n'appartient pas à l'intervalle de confiance de niveau 1−α1-\alpha construit avec la même erreur type.
  • f) Sans refaire le test, déduisez la décision au seuil de 1 %1\,\% de l'intervalle de confiance à 99 %99\,\%.
  • g) Rédigez la conclusion pour la classe, et choisissez une cause physique compatible avec le sens de l'écart.

Tape tes réponses, la page te dit juste ou faux 0/11

a)
b)
c)
d)
f)
g)
Voir la correction

Réponses

  • a) H0:μ=9,81H_{0}:\mu=9{,}81 ; H1:μ≠9,81H_{1}:\mu\neq 9{,}81 (un biais dans un sens ou dans l'autre est un défaut)
  • b) σX‾≈0,0130\sigma_{\overline{X}}\approx 0{,}0130 m/s² ; z≈−2,77z\approx -2{,}77 ; rejeter H0H_{0}
  • c) p≈0,0056p\approx 0{,}0056
  • d) [9,748 ; 9,800][9{,}748\,;\,9{,}800] m/s² ; il ne contient pas 9,819{,}81
  • e) ∣z∣>zα/2  ⟺  ∣x‾−μ0∣>zα/2 σX‾  ⟺  μ0∉[x‾±zα/2 σX‾]|z|>z_{\alpha/2} \iff |\overline{x}-\mu_{0}|>z_{\alpha/2}\,\sigma_{\overline{X}} \iff \mu_{0}\notin\left[\overline{x}\pm z_{\alpha/2}\,\sigma_{\overline{X}}\right]
  • f) [9,740 ; 9,808][9{,}740\,;\,9{,}808] ne contient pas 9,819{,}81 : rejeter H0H_{0} au seuil de 1 %1\,\%
  • g) Le montage donne des valeurs de gg trop faibles en moyenne ; grande amplitude des oscillations
a)

Paramètre : μ\mu, la moyenne de toutes les valeurs de gg que donnerait ce montage. H0:μ=9,81H_{0}:\mu=9{,}81 m/s² (pas de biais) et H1:μ≠9,81H_{1}:\mu\neq 9{,}81 m/s².

Bilatéral : on cherche un défaut du montage, et l'on ne sait pas d'avance s'il surestime ou sous-estime gg. Le choix se fait AVANT de voir que x‾<9,81\overline{x}<9{,}81.

b)

n=50≥30n=50\geq 30 : σX‾≈0,09250≈0,0130\sigma_{\overline{X}}\approx\dfrac{0{,}092}{\sqrt{50}}\approx 0{,}0130 m/s².

z=9,774−9,810,0130≈−2,77\displaystyle z=\dfrac{9{,}774-9{,}81}{0{,}0130}\approx -2{,}77

Bilatéral au seuil de 5 %5\,\% : on rejette si ∣z∣>1,96|z|>1{,}96. Ici ∣z∣≈2,77>1,96|z|\approx 2{,}77>1{,}96 : on rejette H0H_{0}.

c)
p=2×P(Z<−2,77)=2×(1−0,9972)=2×0,0028=0,0056\displaystyle p=2\times P(Z<-2{,}77)=2\times(1-0{,}9972)=2\times 0{,}0028=0{,}0056

Si le montage n'avait aucun biais, une moyenne de 5050 mesures au moins aussi éloignée de 9,819{,}81, d'un côté ou de l'autre, n'apparaîtrait que dans environ 0,56 %0{,}56\,\% des cas.

d)
x‾±1,96 σX‾=9,774±1,96×0,01301=9,774±0,0255\displaystyle \overline{x}\pm 1{,}96\,\sigma_{\overline{X}}=9{,}774\pm 1{,}96\times 0{,}01301=9{,}774\pm 0{,}0255

L'intervalle est [9,748 ; 9,800][9{,}748\,;\,9{,}800] m/s² : il ne contient pas 9,819{,}81, ce qui donne la même décision qu'en b).

e)

Notons e=σX‾e=\sigma_{\overline{X}} l'erreur type, la même dans le test et dans l'intervalle. Le test bilatéral rejette H0H_{0} si et seulement si

∣x‾−μ0e∣>zα/2  ⟺  ∣x‾−μ0∣>zα/2 e\displaystyle \left|\dfrac{\overline{x}-\mu_{0}}{e}\right|>z_{\alpha/2} \iff |\overline{x}-\mu_{0}|>z_{\alpha/2}\,e
  ⟺  μ0<x‾−zα/2 e ou μ0>x‾+zα/2 e\displaystyle \iff \mu_{0}<\overline{x}-z_{\alpha/2}\,e \ \text{ou}\ \mu_{0}>\overline{x}+z_{\alpha/2}\,e
  ⟺  μ0∉[x‾−zα/2 e ; x‾+zα/2 e]\displaystyle \iff \mu_{0}\notin\left[\overline{x}-z_{\alpha/2}\,e\,;\,\overline{x}+z_{\alpha/2}\,e\right]

C'est exactement l'intervalle de confiance de niveau 1−α1-\alpha. Chaque étape est une équivalence, donc les deux règles donnent toujours la même décision.

Deux conditions tiennent cette équivalence : un test BILATÉRAL, et la même erreur type des deux côtés. Pour une proportion, le test utilise p0p_{0} et l'intervalle p^\hat{p} (exercice 4) : l'équivalence n'est alors qu'approximative.

f)

x‾±2,576×0,01301=9,774±0,0335\overline{x}\pm 2{,}576\times 0{,}01301=9{,}774\pm 0{,}0335, soit [9,740 ; 9,808][9{,}740\,;\,9{,}808] m/s².

9,81>9,8089{,}81>9{,}808 : la valeur de référence est hors de l'intervalle à 99 %99\,\%, donc on rejette H0H_{0} au seuil de 1 %1\,\%, sans autre calcul. C'est cohérent avec p=0,0056<0,01p=0{,}0056<0{,}01.

La marge est mince : 9,819{,}81 n'est qu'à 0,0020{,}002 m/s² de la borne. Arrondir trop tôt l'erreur type à 0,0130{,}013 suffirait à brouiller la lecture.

g)

Conclusion : au seuil de 1 %1\,\%, les données des 5050 mesures permettent de conclure que le montage donne en moyenne une valeur de gg différente de 9,819{,}81 m/s², et plus faible.

g=4π2LT2g=\frac{4\pi^{2}L}{T^{2}} sous-estime gg si TT est trop grand. La formule n'est valable que pour de petites oscillations : à grande amplitude, la période réelle s'allonge, ce qui baisse gg à chaque mesure, dans le même sens.

Une erreur aléatoire de chronométrage disperse les mesures sans les décaler toutes du même côté : elle augmente ss, pas l'écart moyen. Mesurer LL jusqu'au centre de la masse est la bonne pratique.

Piège : conclure que gg vaut 9,7749{,}774 m/s² au laboratoire. Le test juge le MONTAGE, pas la valeur de gg.

0−1,961,96z = −2,77Z sous H₀9,729,749,769,789,809,829,84IC à 95 %g₀ = 9,81g (m/s²)

Partie B : problèmes et raisonnement (/50)

Exercice 6 : Les deux erreurs d'un test : le plomb dans l'eau potable

10 réponses vérifiables

On reprend le réseau d'eau de l'exercice 1 b), avec la norme fictive de 88 µg/L. Le laboratoire prélève n=40n=40 échantillons au robinet, tirés au hasard dans le réseau. Ses années de contrôle permettent de considérer l'écart-type des concentrations comme connu : σ=2,4\sigma=2{,}4 µg/L. La municipalité teste H0:μ=8H_{0}:\mu=8 contre H1:μ<8H_{1}:\mu<8 au seuil de 5 %5\,\%, et ne déclare l'eau conforme que si elle rejette H0H_{0}. On donne z0,05=1,645z_{0{,}05}=1{,}645 et z0,01=2,326z_{0{,}01}=2{,}326.

La figure ci-dessous montre la loi de X‾\overline{X} si μ=8\mu=8 (hypothèse nulle) et si μ=7\mu=7 (une eau nettement conforme).

6789μ = 8 (H₀)μ = 7x̄ (µg/L)
  • a) Décrivez, dans ce contexte, l'erreur de première espèce et l'erreur de seconde espèce. Dites laquelle ce choix de H1H_{1} rend rare, et pourquoi.
  • b) Déterminez la région critique en concentration moyenne : x‾<c\overline{x}<c.
  • c) Démontrez que, si μ=8\mu=8, la probabilité de rejeter H0H_{0} vaut 0,050{,}05, dans le modèle normal de X‾\overline{X}.
  • d) La vraie moyenne du réseau est en réalité μ=7\mu=7 µg/L. Calculez le risque de seconde espèce β\beta et la puissance du test.
  • e) Le laboratoire double l'effort : n=80n=80 échantillons, même seuil. Calculez la nouvelle borne cc et le nouveau β\beta pour μ=7\mu=7.
  • f) On revient à n=40n=40, mais au seuil de 1 %1\,\%. Calculez β\beta pour μ=7\mu=7, puis dites ce qu'il faut faire pour réduire les deux risques à la fois.
  • g) Les 4040 échantillons donnent x‾=7,45\overline{x}=7{,}45 µg/L. Calculez zz, prenez la décision et rédigez la conclusion de la municipalité.

Tape tes réponses, la page te dit juste ou faux 0/10

a)
b)
d)
e)
f)
g)
Voir la correction

Réponses

  • a) Première espèce : déclarer conforme une eau qui ne l'est pas (risque sanitaire, contrôlé à 5 %5\,\%) ; seconde espèce : déclarer non conforme une eau conforme
  • b) x‾<7,376\overline{x}<7{,}376 µg/L
  • c) P(X‾<c)=P(Z<−1,645)=0,05P(\overline{X}<c)=P\left(Z<-1{,}645\right)=0{,}05 quand μ=8\mu=8
  • d) β≈0,161\beta\approx 0{,}161 ; puissance ≈0,839\approx 0{,}839
  • e) c≈7,559c\approx 7{,}559 ; β≈0,019\beta\approx 0{,}019
  • f) β≈0,378\beta\approx 0{,}378 ; augmenter nn
  • g) z≈−1,45z\approx -1{,}45 ; ne pas rejeter H0H_{0} : l'eau ne peut pas être déclarée conforme
a)

Erreur de PREMIÈRE espèce : rejeter H0H_{0} alors qu'elle est vraie, c'est-à-dire déclarer l'eau conforme alors que sa concentration moyenne atteint la norme. C'est un risque pour la santé.

Erreur de SECONDE espèce : ne pas rejeter H0H_{0} alors que H1H_{1} est vraie, c'est-à-dire refuser de déclarer conforme une eau qui l'est. Elle coûte de nouvelles analyses ou des travaux inutiles.

Le seuil α=5 %\alpha=5\,\% borne la première, la plus grave : avec H1:μ<8H_{1}:\mu<8, c'est la conformité qui doit être prouvée, et l'erreur que l'on contrôle est celle qui exposerait la population.

b)

Erreur type : σX‾=2,440≈0,3795\sigma_{\overline{X}}=\dfrac{2{,}4}{\sqrt{40}}\approx 0{,}3795 µg/L.

Unilatéral à gauche : on rejette si x‾−80,3795<−1,645\dfrac{\overline{x}-8}{0{,}3795}<-1{,}645, soit

x‾<8−1,645×0,3795≈7,376\overline{x}<8-1{,}645\times 0{,}3795\approx 7{,}376 µg/L

c)

Si μ=8\mu=8, X‾\overline{X} suit approximativement la loi normale de moyenne 88 et d'écart-type σX‾\sigma_{\overline{X}} (n=40≥30n=40\geq 30). Alors Z=X‾−8σX‾Z=\dfrac{\overline{X}-8}{\sigma_{\overline{X}}} suit la loi normale centrée réduite.

P(rejet)=P(X‾<8−1,645 σX‾)=P(X‾−8σX‾<−1,645)=P(Z<−1,645)=0,05\displaystyle P(\text{rejet})=P(\overline{X}<8-1{,}645\,\sigma_{\overline{X}})=P\left(\dfrac{\overline{X}-8}{\sigma_{\overline{X}}}<-1{,}645\right)=P(Z<-1{,}645)=0{,}05

La borne cc a été construite pour cela : le seuil α\alpha EST la probabilité de rejeter H0H_{0} quand elle est vraie, aux arrondis de la table près, et quelle que soit la valeur de σ\sigma et de nn.

d)

Si μ=7\mu=7, X‾\overline{X} est centrée sur 77 avec la même erreur type. On ne rejette pas H0H_{0} quand x‾≥7,376\overline{x}\geq 7{,}376 :

β=P(X‾≥7,376)=P(Z≥7,376−70,3795)=P(Z≥0,99)=1−0,8389=0,1611\displaystyle \beta=P(\overline{X}\geq 7{,}376)=P\left(Z\geq\dfrac{7{,}376-7}{0{,}3795}\right)=P(Z\geq 0{,}99)=1-0{,}8389=0{,}1611

La puissance vaut 1−β≈0,8391-\beta\approx 0{,}839 : une eau à 77 µg/L en moyenne n'est déclarée conforme que dans environ 84 %84\,\% des campagnes de 4040 échantillons. C'est l'aire bleue de la figure du corrigé.

Piège : croire que β=1−α=0,95\beta=1-\alpha=0{,}95. β\beta dépend de la vraie valeur de μ\mu, que le test ne connaît pas ; α\alpha et β\beta ne se calculent pas sous la même loi.

e)

σX‾=2,480≈0,2683\sigma_{\overline{X}}=\dfrac{2{,}4}{\sqrt{80}}\approx 0{,}2683 et c=8−1,645×0,2683≈7,559c=8-1{,}645\times 0{,}2683\approx 7{,}559 µg/L.

β=P(Z≥7,559−70,2683)=P(Z≥2,08)=1−0,9812=0,0188\displaystyle \beta=P\left(Z\geq\dfrac{7{,}559-7}{0{,}2683}\right)=P(Z\geq 2{,}08)=1-0{,}9812=0{,}0188

Doubler nn divise l'erreur type par 2\sqrt{2} : les deux courbes se resserrent, se recouvrent beaucoup moins, et β\beta passe d'environ 16 %16\,\% à environ 2 %2\,\%, à α\alpha inchangé.

f)

c=8−2,326×0,3795≈7,117c=8-2{,}326\times 0{,}3795\approx 7{,}117 µg/L, et

β=P(Z≥7,117−70,3795)=P(Z≥0,31)=1−0,6217=0,3783\displaystyle \beta=P\left(Z\geq\dfrac{7{,}117-7}{0{,}3795}\right)=P(Z\geq 0{,}31)=1-0{,}6217=0{,}3783

À nn fixé, baisser α\alpha recule la borne cc vers la gauche : on rejette moins souvent, donc on rate plus souvent une eau conforme. β\beta passe de 0,1610{,}161 à 0,3780{,}378.

Seule une augmentation de nn réduit les deux risques à la fois, comme en e). Un test bilatéral, lui, rendrait le test MOINS puissant contre μ<8\mu<8.

g)

z=7,45−80,3795≈−1,45>−1,645z=\dfrac{7{,}45-8}{0{,}3795}\approx -1{,}45>-1{,}645 : zz n'est pas dans la région critique, on ne rejette pas H0H_{0}.

Conclusion : au seuil de 5 %5\,\%, les données ne permettent pas d'établir que la concentration moyenne est sous la norme ; l'eau ne peut donc pas être déclarée conforme sur la base de cette campagne.

Piège : écrire « l'eau est contaminée ». La moyenne observée est sous 88 ; elle n'en est simplement pas assez loin pour exclure le hasard. Une campagne de 8080 échantillons, plus puissante, pourrait trancher.

6897,376μ = 8 (H₀)μ = 7αβα = 0,05 ; β ≈ 0,16x̄ (µg/L)

Exercice 7 : Student à 2 degrés de liberté : la valeur de la table démontrée par l'intégrale

9 réponses vérifiables

Une équipe mesure gg en filmant la chute libre d'une bille, et n'a le temps de faire que trois essais : 9,529{,}52 ; 9,779{,}77 ; 9,649{,}64 m/s². On suppose les mesures issues d'une loi normale, et l'on teste H0:μ=9,81H_{0}:\mu=9{,}81 contre H1:μ≠9,81H_{1}:\mu\neq 9{,}81 au seuil de 5 %5\,\%.

Avec n=3n=3, la statistique T=X‾−μ0S/3T=\dfrac{\overline{X}-\mu_{0}}{S/\sqrt{3}} suit, sous H0H_{0}, une loi de Student à 22 degrés de liberté, dont la densité est f(t)=(2+t2)−3/2f(t)=\left(2+t^{2}\right)^{-3/2} pour tout réel tt. La figure compare cette densité à celle de la loi normale centrée réduite.

−6−4−2024623456loi normaleStudent, 2 ddlqueue droite agrandie (hauteurs × 4)
  • a) Montrez que F(t)=12+t22+t2F(t)=\dfrac{1}{2}+\dfrac{t}{2\sqrt{2+t^{2}}} vérifie F′(t)=f(t)F'(t)=f(t), puis calculez ses limites en −∞-\infty et en +∞+\infty. Concluez que ff est une densité de probabilité et que FF est sa fonction de répartition.
  • b) Déduisez-en la valeur critique t0,025 ; 2t_{0{,}025\,;\,2}, telle que P(T>t0,025 ; 2)=0,025P(T>t_{0{,}025\,;\,2})=0{,}025, en résolvant F(t)=0,975F(t)=0{,}975. Calculez de même t0,05 ; 2t_{0{,}05\,;\,2}.
  • c) Calculez x‾\overline{x}, ss et la statistique tt.
  • d) Prenez la décision au seuil de 5 %5\,\% et rédigez la conclusion.
  • e) Un camarade compare ∣t∣|t| à 1,961{,}96, « comme d'habitude », et rejette H0H_{0}. Calculez le vrai risque de première espèce de sa règle, P(∣T∣>1,96)P(|T|>1{,}96).
  • f) Montrez que la valeur pp du test vaut 1−∣t∣2+t21-\dfrac{|t|}{\sqrt{2+t^{2}}}, puis calculez-la.
  • g) À l'exercice 5, cinquante mesures détectaient un écart de 0,0360{,}036 m/s². Ici l'écart observé est de 0,170{,}17 m/s², et le test ne conclut pas. Expliquez.

Tape tes réponses, la page te dit juste ou faux 0/9

b)
c)
d)
e)
f)
g)
Voir la correction

Réponses

  • a) F′=fF'=f, lim⁡−∞F=0\lim_{-\infty}F=0, lim⁡+∞F=1\lim_{+\infty}F=1 : f≥0f\geq 0 et ∫−∞+∞f=1\int_{-\infty}^{+\infty}f=1
  • b) t0,025 ; 2=2×0,9521−0,952≈4,303t_{0{,}025\,;\,2}=\sqrt{\frac{2\times 0{,}95^{2}}{1-0{,}95^{2}}}\approx 4{,}303 ; t0,05 ; 2≈2,920t_{0{,}05\,;\,2}\approx 2{,}920
  • c) x‾≈9,643\overline{x}\approx 9{,}643 m/s² ; s≈0,125s\approx 0{,}125 m/s² ; t≈−2,31t\approx -2{,}31
  • d) Ne pas rejeter H0H_{0} : trois mesures ne permettent pas de conclure à un biais
  • e) P(∣T∣>1,96)≈0,189P(|T|>1{,}96)\approx 0{,}189, près de quatre fois le seuil annoncé
  • f) p=1−2,3092+2,3092≈0,147p=1-\frac{2{,}309}{\sqrt{2+2{,}309^{2}}}\approx 0{,}147
  • g) Avec n=3n=3, l'erreur type et la valeur critique sont grandes : le test a peu de puissance
a)

On dérive le quotient u(t)=t2+t2u(t)=\dfrac{t}{\sqrt{2+t^{2}}} :

u′(t)=2+t2−t⋅t2+t22+t2=(2+t2)−t2(2+t2)3/2=2(2+t2)3/2\displaystyle u'(t)=\dfrac{\sqrt{2+t^{2}}-t\cdot\frac{t}{\sqrt{2+t^{2}}}}{2+t^{2}}=\dfrac{(2+t^{2})-t^{2}}{(2+t^{2})^{3/2}}=\dfrac{2}{(2+t^{2})^{3/2}}

Donc F′(t)=12 u′(t)=(2+t2)−3/2=f(t)F'(t)=\dfrac{1}{2}\,u'(t)=(2+t^{2})^{-3/2}=f(t).

Pour ∣t∣|t| grand, 2+t2=∣t∣1+2/t2\sqrt{2+t^{2}}=|t|\sqrt{1+2/t^{2}}, donc u(t)→−1u(t)\to -1 en −∞-\infty et u(t)→1u(t)\to 1 en +∞+\infty. Ainsi lim⁡−∞F=0\lim_{-\infty}F=0 et lim⁡+∞F=1\lim_{+\infty}F=1.

f(t)>0f(t)>0 et ∫−∞+∞f(t) dt=lim⁡+∞F−lim⁡−∞F=1\int_{-\infty}^{+\infty}f(t)\,dt=\lim_{+\infty}F-\lim_{-\infty}F=1 : ff est une densité. Comme F′=fF'=f et F→0F\to 0 en −∞-\infty, F(x)=∫−∞xf(t) dt=P(T≤x)F(x)=\int_{-\infty}^{x}f(t)\,dt=P(T\leq x) : c'est la fonction de répartition.

b)

F(t)=0,975  ⟺  t2+t2=0,95F(t)=0{,}975 \iff \dfrac{t}{\sqrt{2+t^{2}}}=0{,}95, avec t>0t>0. On élève au carré :

t2=0,9025 (2+t2)  ⟺  t2=2×0,90250,0975≈18,513\displaystyle t^{2}=0{,}9025\,(2+t^{2}) \iff t^{2}=\dfrac{2\times 0{,}9025}{0{,}0975}\approx 18{,}513

t0,025 ; 2≈4,303t_{0{,}025\,;\,2}\approx 4{,}303, exactement la valeur de la table. De même, F(t)=0,95F(t)=0{,}95 donne t2+t2=0,9\dfrac{t}{\sqrt{2+t^{2}}}=0{,}9, t2=2×0,810,19t^{2}=\dfrac{2\times 0{,}81}{0{,}19} et t0,05 ; 2≈2,920t_{0{,}05\,;\,2}\approx 2{,}920.

La table de Student n'est donc pas une boîte noire : chacune de ses valeurs est une borne d'intégrale de la densité.

c)

x‾=9,52+9,77+9,643≈9,643\overline{x}=\dfrac{9{,}52+9{,}77+9{,}64}{3}\approx 9{,}643 m/s²

s=(−0,1233)2+0,12672+(−0,0033)22≈0,125s=\sqrt{\dfrac{(-0{,}1233)^{2}+0{,}1267^{2}+(-0{,}0033)^{2}}{2}}\approx 0{,}125 m/s²

t=9,643−9,810,125/3≈−0,16670,0722≈−2,31\displaystyle t=\dfrac{9{,}643-9{,}81}{0{,}125/\sqrt{3}}\approx\dfrac{-0{,}1667}{0{,}0722}\approx -2{,}31
d)

Bilatéral à 22 degrés de liberté : on rejette si ∣t∣>4,303|t|>4{,}303. Ici ∣t∣≈2,31<4,303|t|\approx 2{,}31<4{,}303 : on ne rejette pas H0H_{0}.

Conclusion : au seuil de 5 %5\,\%, ces trois essais ne permettent pas de conclure que la chute filmée donne une valeur de gg biaisée.

e)

Par symétrie de ff, P(∣T∣>1,96)=2(1−F(1,96))=1−1,962+1,962≈1−0,8109=0,189P(|T|>1{,}96)=2\left(1-F(1{,}96)\right)=1-\dfrac{1{,}96}{\sqrt{2+1{,}96^{2}}}\approx 1-0{,}8109=0{,}189

La règle « ∣t∣>1,96|t|>1{,}96 » annonce un risque de 5 %5\,\% mais en prend près de 19 %19\,\% : avec deux degrés de liberté, les queues de Student sont beaucoup plus lourdes que celles de la loi normale (aire rouge de la figure du corrigé).

Piège : utiliser la loi normale avec un petit échantillon. Le rejet du camarade n'est pas une découverte, c'est l'effet de la mauvaise table.

f)
p=P(∣T∣≥∣t∣)=2(1−F(∣t∣))=2(12−∣t∣22+t2)=1−∣t∣2+t2\displaystyle p=P(|T|\geq |t|)=2\left(1-F(|t|)\right)=2\left(\dfrac{1}{2}-\dfrac{|t|}{2\sqrt{2+t^{2}}}\right)=1-\dfrac{|t|}{\sqrt{2+t^{2}}}

Avec ∣t∣≈2,309|t|\approx 2{,}309 : p≈1−2,3097,331≈1−0,853=0,147>0,05p\approx 1-\dfrac{2{,}309}{\sqrt{7{,}331}}\approx 1-0{,}853=0{,}147>0{,}05, même décision qu'en d).

g)

À l'exercice 5, l'erreur type valait 0,0130{,}013 m/s² et la valeur critique 1,961{,}96. Ici, s/3≈0,072s/\sqrt{3}\approx 0{,}072 m/s² et la valeur critique 4,3034{,}303 : l'écart doit dépasser environ 4,303×0,072≈0,314{,}303\times 0{,}072\approx 0{,}31 m/s² pour être détecté.

Un biais réel de 0,170{,}17 m/s² passerait donc le plus souvent inaperçu : le risque de seconde espèce est grand, la puissance faible.

Ne pas rejeter H0H_{0} ne prouve pas que la chute filmée est juste ; cela dit que trois essais ne suffisent pas à trancher. Pour conclure, il faut plus d'essais, ce que l'exercice 6 chiffre avec β\beta.

0−4,3034,303−1,961,961,9634,30356loi normaleStudent, 2 ddlqueue droite agrandie (hauteurs × 4)

Exercice 8 : Cinq affirmations à corriger

9 réponses vérifiables

On reprend le catalyseur de l'exercice 1 e). Sans catalyseur, la réaction dure en moyenne μ0=120\mu_{0}=120 s. La chimiste teste H0:μ=120H_{0}:\mu=120 contre H1:μ<120H_{1}:\mu<120 au seuil de 5 %5\,\%, sur n=36n=36 essais avec catalyseur ; l'écart-type des durées est d'environ s=18s=18 s, donc l'erreur type vaut 33 s.

Chacune des cinq affirmations suivantes, entendues au laboratoire, est fausse. Pour chacune, choisissez l'énoncé correct, faites le calcul demandé et expliquez l'erreur.

  • a) « Le test donne une valeur pp de 0,030{,}03 : il y a donc 3 %3\,\% de chances que le catalyseur soit sans effet. »
  • b) « Avec x‾=116\overline{x}=116 s, le test ne rejette pas H0H_{0} : le catalyseur n'a donc aucun effet. » Calculez zz.
  • c) « Le test est au seuil de 5 %5\,\%, donc on rejette H0H_{0} dès que ∣z∣>1,96|z|>1{,}96. » Calculez zz pour x‾=114,6\overline{x}=114{,}6 s et prenez la bonne décision.
  • d) « Au seuil de 5 %5\,\%, quand on rejette H0H_{0}, on a 5 %5\,\% de chances de se tromper. » Le laboratoire estime que 10 %10\,\% des catalyseurs qu'il essaie sont efficaces, et que le test détecte un catalyseur efficace avec une puissance de 0,800{,}80. Calculez la probabilité qu'un catalyseur déclaré efficace ne le soit pas.
  • e) « Augmenter le nombre d'essais diminue le risque de première espèce. » Pour un catalyseur dont la vraie moyenne est 112112 s, calculez β\beta avec 3636 essais, puis avec 100100 essais (erreur type 1,81{,}8 s).

Tape tes réponses, la page te dit juste ou faux 0/9

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) La valeur pp se calcule SI H0H_{0} est vraie ; ce n'est pas la probabilité de H0H_{0}
  • b) z≈−1,33z\approx -1{,}33 : pas de preuve d'effet au seuil de 5 %5\,\%, ce qui n'est pas une preuve d'absence d'effet
  • c) z=−1,80<−1,645z=-1{,}80<-1{,}645 : rejeter H0H_{0} (la règle ∣z∣>1,96|z|>1{,}96 est celle du test bilatéral)
  • d) P(sans effet∣rejet)=0,36P(\text{sans effet}\mid\text{rejet})=0{,}36, et non 0,050{,}05
  • e) β≈0,154\beta\approx 0{,}154 puis β≈0,003\beta\approx 0{,}003 : nn fait baisser β\beta, α\alpha reste à 5 %5\,\%
a)

La valeur pp se calcule en SUPPOSANT H0H_{0} vraie : c'est P(un eˊcart au moins aussi grand∣H0)P(\text{un écart au moins aussi grand}\mid H_{0}).

Correction : « si le catalyseur était sans effet, une durée moyenne aussi courte ou plus courte n'apparaîtrait que dans 3 %3\,\% des séries de 3636 essais ».

La probabilité que H0H_{0} soit vraie sachant les données est une autre quantité, P(H0∣donneˊes)P(H_{0}\mid\text{données}), que le test ne calcule pas ; la question d) montre à quel point elle peut différer.

b)

z=116−1203≈−1,33>−1,645z=\dfrac{116-120}{3}\approx -1{,}33>-1{,}645 : on ne rejette pas H0H_{0}.

Correction : « au seuil de 5 %5\,\%, les données ne permettent pas de conclure que le catalyseur réduit la durée moyenne ».

Une durée moyenne de 116116 s est compatible avec un catalyseur sans effet, mais aussi avec un catalyseur qui gagne quelques secondes : l'absence de preuve n'est pas la preuve d'une absence.

c)
z=114,6−1203=−1,80\displaystyle z=\dfrac{114{,}6-120}{3}=-1{,}80

Le test est UNILATÉRAL à gauche : toute l'aire α=0,05\alpha=0{,}05 est dans la queue gauche, et l'on rejette si z<−1,645z<-1{,}645. Ici −1,80<−1,645-1{,}80<-1{,}645 : on rejette H0H_{0}.

La règle ∣z∣>1,96|z|>1{,}96 est celle du test bilatéral, qui partage α\alpha entre deux queues. Ici elle aurait conduit à ne pas rejeter : l'erreur change la conclusion.

d)

Notons EE « le catalyseur est efficace » et RR « le test rejette H0H_{0} ». P(E)=0,10P(E)=0{,}10, P(R∣E)=0,80P(R\mid E)=0{,}80 (puissance), P(R∣E‾)=0,05P(R\mid\overline{E})=0{,}05 (seuil).

P(E‾∣R)=P(E‾) P(R∣E‾)P(E‾) P(R∣E‾)+P(E) P(R∣E)=0,9×0,050,045+0,08=0,0450,125=0,36\displaystyle P(\overline{E}\mid R)=\dfrac{P(\overline{E})\,P(R\mid\overline{E})}{P(\overline{E})\,P(R\mid\overline{E})+P(E)\,P(R\mid E)}=\dfrac{0{,}9\times 0{,}05}{0{,}045+0{,}08}=\dfrac{0{,}045}{0{,}125}=0{,}36

Plus d'un catalyseur « efficace » sur trois ne l'est pas. Le seuil α=P(R∣E‾)\alpha=P(R\mid\overline{E}) n'est pas P(E‾∣R)P(\overline{E}\mid R) : on retrouve la formule de Bayes, et la confusion entre une probabilité conditionnelle et sa réciproque.

Correction : « quand H0H_{0} est vraie, on la rejette dans 5 %5\,\% des cas ».

e)

Avec 3636 essais, on rejette si x‾<120−1,645×3=115,065\overline{x}<120-1{,}645\times 3=115{,}065. Si μ=112\mu=112 :

β=P(Z≥115,065−1123)=P(Z≥1,02)=1−0,8461≈0,154\displaystyle \beta=P\left(Z\geq\dfrac{115{,}065-112}{3}\right)=P(Z\geq 1{,}02)=1-0{,}8461\approx 0{,}154

Avec 100100 essais, c=120−1,645×1,8=117,039c=120-1{,}645\times 1{,}8=117{,}039 et β=P(Z≥117,039−1121,8)=P(Z≥2,80)=1−0,9974≈0,003\beta=P\left(Z\geq\dfrac{117{,}039-112}{1{,}8}\right)=P(Z\geq 2{,}80)=1-0{,}9974\approx 0{,}003.

Correction : « augmenter nn diminue le risque de seconde espèce β\beta ; le risque de première espèce reste égal au seuil choisi, 5 %5\,\% ». α\alpha est fixé par l'expérimentateur, pas par la taille de l'échantillon.

Exercice 9 : Problème : une étude complète, la floraison a-t-elle raccourci ?

9 réponses vérifiables

On reprend la plante de sous-bois de l'exercice 1 g), une trille d'une érablière. Des relevés anciens (valeurs fictives) donnent une durée moyenne de floraison de 18,018{,}0 jours par plant. Ce printemps, une botaniste suit n=14n=14 plants choisis au hasard dans la population et note la durée de floraison de chacun, en jours : 15;18;14;19;17;20;13;17;16;15;18;16;17;1915 ; 18 ; 14 ; 19 ; 17 ; 20 ; 13 ; 17 ; 16 ; 15 ; 18 ; 16 ; 17 ; 19.

Elle veut savoir si la durée moyenne de floraison a diminué, au seuil de 5 %5\,\%. Valeurs de Student à 1313 degrés de liberté : t0,05 ; 13=1,771t_{0{,}05\,;\,13}=1{,}771, t0,025 ; 13=2,160t_{0{,}025\,;\,13}=2{,}160, t0,01 ; 13=2,650t_{0{,}01\,;\,13}=2{,}650. Le diagramme en points ci-dessous représente les données.

12131415161718192021référence : 18 jdurée de floraison (jours)
  • a) Calculez la moyenne x‾\overline{x} et l'écart-type ss de l'échantillon.
  • b) Choisissez la loi de la statistique de test. Justifiez ce choix et dites ce que le diagramme permet de vérifier.
  • c) Formulez H0H_{0} et H1H_{1}, puis calculez la statistique tt.
  • d) Prenez la décision au seuil de 5 %5\,\%.
  • e) Encadrez la valeur pp à l'aide des valeurs de la table.
  • f) Calculez la borne supérieure de confiance à 95 %95\,\%, x‾+t0,05 ; 13 sn\overline{x}+t_{0{,}05\,;\,13}\,\dfrac{s}{\sqrt{n}}, et montrez qu'elle conduit à la même décision.
  • g) Rédigez la conclusion de l'étude en trois phrases : le résultat, sa portée, ses limites. Le test permet-il d'attribuer la baisse au réchauffement du climat ?

Tape tes réponses, la page te dit juste ou faux 0/9

a)
b)
c)
d)
e)
f)
g)
Voir la correction

Réponses

  • a) x‾≈16,71\overline{x}\approx 16{,}71 j ; s≈2,02s\approx 2{,}02 j
  • b) Student à 1313 degrés de liberté : n<30n<30, σ\sigma inconnu, population supposée normale ; le diagramme ne montre ni forte asymétrie ni valeur aberrante
  • c) H0:μ=18H_{0}:\mu=18 ; H1:μ<18H_{1}:\mu<18 ; t≈−2,39t\approx -2{,}39
  • d) Rejeter H0H_{0} : t<−1,771t<-1{,}771
  • e) 0,01<p<0,0250{,}01<p<0{,}025
  • f) 17,67<1817{,}67<18 : la valeur de référence est exclue, même décision
  • g) Durée moyenne plus courte que 1818 j dans cette population ce printemps ; une seule année, cause non établie : non
a)

x‾=23414≈16,71\overline{x}=\dfrac{234}{14}\approx 16{,}71 jours.

s=∑(xi−x‾)213=52,85713≈2,02s=\sqrt{\dfrac{\sum (x_{i}-\overline{x})^{2}}{13}}=\sqrt{\dfrac{52{,}857}{13}}\approx 2{,}02 jours

b)

n=14<30n=14<30 et σ\sigma est inconnu : si les durées de floraison suivent une loi normale dans la population, T=X‾−μ0S/nT=\dfrac{\overline{X}-\mu_{0}}{S/\sqrt{n}} suit une loi de Student à 1313 degrés de liberté.

La normalité ne se prouve pas avec 1414 valeurs, mais le diagramme la rend plausible : une seule bosse centrée vers 1616 ou 1717 jours, une étendue de 1313 à 2020 jours, pas d'asymétrie marquée, aucune valeur isolée.

Une durée entière en jours est une variable discrète ; la loi normale en est un modèle, raisonnable quand les valeurs prises sont nombreuses.

c)

μ\mu : durée moyenne de floraison actuelle de la population. H0:μ=18H_{0}:\mu=18 et H1:μ<18H_{1}:\mu<18, unilatéral à gauche.

t=16,714−182,016/14=−1,2860,539≈−2,39\displaystyle t=\dfrac{16{,}714-18}{2{,}016/\sqrt{14}}=\dfrac{-1{,}286}{0{,}539}\approx -2{,}39
d)

On rejette H0H_{0} si t<−1,771t<-1{,}771. Ici t≈−2,39<−1,771t\approx -2{,}39<-1{,}771 : on rejette H0H_{0} au seuil de 5 %5\,\%.

e)

La valeur pp est l'aire à gauche de −2,39-2{,}39, égale par symétrie à l'aire à droite de 2,392{,}39.

2,160<2,39<2,6502{,}160<2{,}39<2{,}650, donc 0,01<p<0,0250{,}01<p<0{,}025 : on rejetterait aussi au seuil de 2,5 %2{,}5\,\%, mais pas au seuil de 1 %1\,\%.

f)

16,714+1,771×2,01614≈16,714+0,954≈17,6716{,}714+1{,}771\times\dfrac{2{,}016}{\sqrt{14}}\approx 16{,}714+0{,}954\approx 17{,}67 jours.

Avec une confiance de 95 %95\,\%, la durée moyenne actuelle est inférieure à 17,6717{,}67 jours. La valeur 1818 est au-dessus de cette borne, donc exclue : on rejette H0H_{0}, comme en d).

C'est la même équivalence qu'à l'exercice 5, en version unilatérale : un test unilatéral au seuil de 5 %5\,\% correspond à une borne de confiance UNILATÉRALE à 95 %95\,\%, et non à l'intervalle bilatéral à 95 %95\,\%.

g)

Résultat : au seuil de 5 %5\,\%, les données permettent de conclure que la durée moyenne de floraison de cette population, ce printemps, est inférieure aux 1818 jours des relevés anciens ; elle est estimée à 16,716{,}7 jours (0,01<p<0,0250{,}01<p<0{,}025).

Portée : la conclusion vaut pour cette érablière et ce printemps, si les 1414 plants ont bien été tirés au hasard.

Limites : une seule saison a été observée, et la référence ancienne a pu être mesurée autrement. Le test compare deux moyennes, il ne désigne aucune cause : chaleur, sécheresse, ombrage ou broutage peuvent raccourcir une floraison.

Piège : écrire « le réchauffement climatique a raccourci la floraison ». Attribuer une cause demande des séries sur plusieurs années, avec la température mesurée, et une étude conçue pour cela.

0−1,77124t = −2,39rejet : 5 %T sous H₀, 13 ddl

Exercice 10 : Problème : surveiller douze étangs, fausses alarmes et puissance

11 réponses vérifiables

Après l'étude de l'exercice 4, un programme régional surveille chaque année 1212 étangs. Dans chacun, on examine n=150n=150 grenouilles capturées au hasard, et l'on teste H0:p=0,05H_{0}:p=0{,}05 contre H1:p>0,05H_{1}:p>0{,}05, comme à l'exercice 4, avec l'approximation normale (np0=7,5≥5np_{0}=7{,}5\geq 5) sans correction de continuité.

Le diagramme ci-dessous donne le nombre de grenouilles malformées trouvées cette année dans chaque étang. On donne z0,05=1,645z_{0{,}05}=1{,}645 et z0,10=1,282z_{0{,}10}=1{,}282.

4812161629344115768716859101071112128pointillé : 7,5 attendus sous H₀étangmalforméessur 150
  • a) Au seuil de 5 %5\,\% par étang, déterminez le plus petit nombre kk de grenouilles malformées sur 150150 qui fait rejeter H0H_{0}. Combien d'étangs sont signalés cette année ?
  • b) Supposons que les 1212 étangs soient tous sains (p=0,05p=0{,}05) et que les tests soient indépendants, chacun avec un risque de première espèce de 5 %5\,\%. Soit NN le nombre de fausses alarmes. Donnez la loi de NN, puis calculez P(N≥1)P(N\geq 1) et E(N)E(N).
  • c) Quel seuil aa faudrait-il prendre dans chaque étang pour que la probabilité d'au moins une fausse alarme sur les douze soit de 5 %5\,\% ? Ce seuil correspond à za≈2,63z_{a}\approx 2{,}63 : calculez le nouveau nombre minimal k′k' et dites combien d'étangs restent signalés.
  • d) Au seuil de 5 %5\,\% par étang, calculez la puissance du test contre un étang où p=0,10p=0{,}10.
  • e) Montrez que, pour détecter p1=0,10p_{1}=0{,}10 avec une puissance de 0,900{,}90 au seuil de 5 %5\,\%, il faut n≥z0,05p0(1−p0)+z0,10p1(1−p1)p1−p0\sqrt{n}\geq\dfrac{z_{0{,}05}\sqrt{p_{0}(1-p_{0})}+z_{0{,}10}\sqrt{p_{1}(1-p_{1})}}{p_{1}-p_{0}}, puis calculez le plus petit nn.
  • f) Rédigez la recommandation du programme pour l'étang 7.

Tape tes réponses, la page te dit juste ou faux 0/11

a)
b)
c)
d)
e)
f)
Voir la correction

Réponses

  • a) k=12k=12 ; deux étangs signalés, le 7 et le 11
  • b) N∼B(12 ; 0,05)N\sim\mathcal{B}(12\,;\,0{,}05) ; P(N≥1)=1−0,9512≈0,460P(N\geq 1)=1-0{,}95^{12}\approx 0{,}460 ; E(N)=0,6E(N)=0{,}6
  • c) a=1−0,951/12≈0,00427a=1-0{,}95^{1/12}\approx 0{,}00427 ; k′=15k'=15 ; seul l'étang 7 reste signalé
  • d) Puissance ≈0,802\approx 0{,}802
  • e) n≥14,86\sqrt{n}\geq 14{,}86, soit n≥221n\geq 221
  • f) Étang 7 : 1616 malformées sur 150150, signal maintenu même avec le seuil corrigé ; enquête recommandée
a)

σP^=0,05×0,95150≈0,01780\sigma_{\hat{P}}=\sqrt{\dfrac{0{,}05\times 0{,}95}{150}}\approx 0{,}01780. On rejette si p^>0,05+1,645×0,01780≈0,07927\hat{p}>0{,}05+1{,}645\times 0{,}01780\approx 0{,}07927.

En nombre : 150×0,07927≈11,89150\times 0{,}07927\approx 11{,}89, donc il faut au moins k=12k=12 grenouilles malformées.

Sur le diagramme, deux étangs atteignent 1212 : l'étang 7 (1616) et l'étang 11 (1212). L'étang 4, avec 1111, n'est pas signalé.

b)

Chaque étang sain déclenche une fausse alarme avec probabilité 0,050{,}05, indépendamment des autres : NN suit la loi binomiale de paramètres 1212 et 0,050{,}05.

P(N≥1)=1−P(N=0)=1−0,9512≈1−0,540=0,460\displaystyle P(N\geq 1)=1-P(N=0)=1-0{,}95^{12}\approx 1-0{,}540=0{,}460
E(N)=12×0,05=0,6\displaystyle E(N)=12\times 0{,}05=0{,}6

Même si tous les étangs étaient sains, le programme signalerait au moins un étang presque une année sur deux. Chaque test a un risque de 5 %5\,\% ; une famille de douze tests en a un bien plus grand.

c)

On veut 1−(1−a)12=0,051-(1-a)^{12}=0{,}05, soit (1−a)12=0,95(1-a)^{12}=0{,}95 :

a=1−0,951/12≈0,00427\displaystyle a=1-0{,}95^{1/12}\approx 0{,}00427

Avec za≈2,63z_{a}\approx 2{,}63 : p^>0,05+2,63×0,01780≈0,0968\hat{p}>0{,}05+2{,}63\times 0{,}01780\approx 0{,}0968, soit plus de 150×0,0968≈14,52150\times 0{,}0968\approx 14{,}52 malformées : k′=15k'=15.

Seul l'étang 7 (16≥1516\geq 15) reste signalé. Le signal de l'étang 11 est compatible avec une fausse alarme parmi douze tests.

Ce partage du risque entre plusieurs tests dépasse le minimum du devis ; il ne demande que la loi binomiale et l'indépendance.

d)

Si p=0,10p=0{,}10, P^\hat{P} suit approximativement la loi normale de moyenne 0,100{,}10 et d'écart-type 0,10×0,90150≈0,02449\sqrt{\dfrac{0{,}10\times 0{,}90}{150}}\approx 0{,}02449.

Puissance =P(P^>0,07927)=P(Z>0,07927−0,100,02449)=P(Z>−0,85)=0,8023=P(\hat{P}>0{,}07927)=P\left(Z>\dfrac{0{,}07927-0{,}10}{0{,}02449}\right)=P(Z>-0{,}85)=0{,}8023

Un étang dont le taux a doublé échappe au programme environ une année sur cinq : β≈0,20\beta\approx 0{,}20.

Piège : calculer la puissance avec l'erreur type de H0H_{0}. Sous l'alternative, P^\hat{P} se disperse autour de p1p_{1}, avec p1(1−p1)/n\sqrt{p_{1}(1-p_{1})/n}.

e)

On rejette si p^>c\hat{p}>c, avec c=p0+z0,05p0(1−p0)nc=p_{0}+z_{0{,}05}\sqrt{\dfrac{p_{0}(1-p_{0})}{n}} (le seuil fixe α\alpha).

La puissance vaut 0,900{,}90 ou plus si P(P^>c)≥0,90P(\hat{P}>c)\geq 0{,}90 quand p=p1p=p_{1}, c'est-à-dire si c≤p1−z0,10p1(1−p1)nc\leq p_{1}-z_{0{,}10}\sqrt{\dfrac{p_{1}(1-p_{1})}{n}}.

On remplace cc et l'on multiplie par n>0\sqrt{n}>0 : n (p1−p0)≥z0,05p0(1−p0)+z0,10p1(1−p1)\sqrt{n}\,(p_{1}-p_{0})\geq z_{0{,}05}\sqrt{p_{0}(1-p_{0})}+z_{0{,}10}\sqrt{p_{1}(1-p_{1})}, d'où l'inégalité annoncée.

n≥1,645×0,2179+1,282×0,30,05≈0,74310,05≈14,86\displaystyle \sqrt{n}\geq\dfrac{1{,}645\times 0{,}2179+1{,}282\times 0{,}3}{0{,}05}\approx\dfrac{0{,}7431}{0{,}05}\approx 14{,}86

n≥220,9n\geq 220{,}9 : il faut examiner au moins 221221 grenouilles par étang.

f)

L'étang 7 compte 1616 malformées sur 150150, soit p^≈0,107\hat{p}\approx 0{,}107 et z=0,1067−0,050,0178≈3,18z=\dfrac{0{,}1067-0{,}05}{0{,}0178}\approx 3{,}18 : il dépasse le seuil corrigé 2,632{,}63.

Recommandation : son signal ne s'explique pas raisonnablement par le nombre de tests ; une enquête sur cet étang est justifiée (eau, sédiments, parasites), avec 221221 captures l'an prochain pour un test assez puissant.

L'étang 11 sera simplement suivi : un signal isolé parmi douze tests est attendu une année sur deux. Comme à l'exercice 4, aucun de ces tests ne désigne la cause des malformations.

Il manque quelque chose dans cette série ? Un type d'exercice que ton prof donne, un énoncé qui te bloque, une erreur repérée : dis-le-moi, je lis chaque demande.

L'envoi passe par ton espace gratuit : c'est ce qui me permet de te répondre et de te prévenir quand l'exercice est ajouté. Un courriel suffit, pas de mot de passe. Crée ton espace. Pas envie de créer un compte ? Écris-moi directement.

De quoi s'agit-il ?
Je lis chaque demande moi-même, personne d'autre ne la voit. Inutile d'écrire ton nom ou celui de ton école.
Quand ta demande est traitée, tu es prévenu par courriel et dans ton espace.

En envoyant, tu acceptes que je lise ta demande pour y répondre. Elle est gardée avec ton compte, effacée avec lui, et jamais publiée : politique de confidentialité.

Ce chapitre resservira dans

Les chapitres qui le réclament en amont, plus tard dans l'année ou dans les années suivantes.

© Ahmed Squalli Houssaini. Série publiée sur www.letuteurscientifique.ca/exercices/sn1-tests-hypotheses. Libre pour l'usage personnel et en classe ; sa republication ailleurs demande une autorisation écrite (mentions légales).

Vous cherchez un tuteur en probabilités et statistique 201-SN1-RE à Montréal ?

Contactez-moi pour une première séance. Un test d'hypothèse se corrige d'abord sur la démarche : les hypothèses sur le bon paramètre, la bonne loi, la décision et une conclusion rédigée dans le contexte. C'est là que se jouent la plupart des points de l'épreuve.

Site par Studio Squalli