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

Épreuve de synthèse corrigée : probabilités et statistique (201-SN1-RE)

Une série de chapitre annonce son sujet ; une épreuve finale ne l'annonce pas, et c'est ce qu'il faut travailler avant l'examen.

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

Cette épreuve de synthèse couvre le cours 201-SN1-RE, Probabilités et statistique, compétence 0M01 du programme révisé de Sciences de la nature : statistique descriptive, probabilités (dénombrement, probabilités conditionnelles et formule de Bayes, variables aléatoires discrètes et continues, lois binomiale et normale), inférence (théorème central limite, intervalles de confiance et tests d'hypothèse) et lien entre deux variables (droite des moindres carrés, corrélation, test d'indépendance du khi-carré).

Onze exercices, 100 points, trois heures, dans les proportions du devis : la statistique descriptive (exercice 1), les probabilités, la plus grosse partie (exercices 2 à 6), l'inférence (exercices 7 à 9), la régression (exercice 10) et un problème de synthèse sur une étude environnementale, le sel de déglaçage et les grenouilles des bois (exercice 11). Les six cas du cours sont tous là : intervalle et test sur une moyenne avec un grand échantillon (exercices 11 et 7), avec la loi de Student (exercice 8), et sur une proportion (exercice 9).

À la fin, vous savez démontrer ce que le cours démontre (Bayes à l'exercice 3, E(aX+b)E(aX+b) et V(aX+b)V(aX+b) à l'exercice 4, l'espérance d'une densité par une intégrale impropre à l'exercice 5, V(X‾)=σ2nV(\overline{X})=\frac{\sigma^{2}}{n} à l'exercice 7), choisir la loi qui convient et rédiger une décision dans le contexte. Les pièges traités : l'écart-type divisé par nn (exercice 1), compter deux fois en dénombrant « au moins un » (exercice 2), confondre P(L∣C)P(L\mid C) et P(C∣L)P(C\mid L), incompatibles et indépendants (exercice 3), ajouter un coût fixe à un écart-type (exercice 4), oublier la correction de continuité (exercice 6), modéliser une maison par la loi normale (exercice 7), comparer à 1,961{,}96 avec 99 mesures ou choisir H1H_{1} après coup (exercice 8), mettre p^\hat{p} dans l'erreur type d'un test (exercice 9), extrapoler une droite (exercice 10) et lire une causalité dans un khi-carré (exercice 11).

Le fil de l'épreuve : décider sous incertitude, c'est d'abord nommer le modèle. Chaque calcul repose sur une hypothèse dite, indépendance, normalité, tirage avec ou sans remise, taille de l'échantillon, et la conclusion ne vaut que ce que vaut l'hypothèse. Faites l'épreuve d'une traite, calculatrice en main, avant d'ouvrir un seul corrigé ; chaque exercice nomme les séries du cours à reprendre.

  • 11 exercices corrigés
  • Durée : 3 h 00
  • Autocorrigée, 102 réponses
  • Calculatrice
Commencer l'exercice 1

Rappel de cours

  • • DESCRIPTIVE : s2=∑ni(ci−x‾)2n−1s^{2}=\dfrac{\sum n_{i}(c_{i}-\overline{x})^{2}}{n-1} pour un échantillon ; médiane et centiles par interpolation, Pk=L+k100n−Fnc×aP_{k}=L+\dfrac{\frac{k}{100}n-F}{n_{c}}\times a ; z=x−x‾sz=\dfrac{x-\overline{x}}{s} ; aberrante hors de [Q1−1,5 EI ;Q3+1,5 EI][Q_{1}-1{,}5\,EI\,;Q_{3}+1{,}5\,EI].
  • • DÉNOMBREMENT : A(n,k)=n!(n−k)!A(n,k)=\dfrac{n!}{(n-k)!} si l'ordre compte, (nk)=n!k! (n−k)!\dbinom{n}{k}=\dfrac{n!}{k!\,(n-k)!} sinon ; permutations d'objets identiques n!n1! n2!⋯\dfrac{n!}{n_{1}!\,n_{2}!\cdots} ; « au moins un » par le contraire.
  • • PROBABILITÉS : P(A∣B)=P(A∩B)P(B)P(A\mid B)=\dfrac{P(A\cap B)}{P(B)} ; probabilités totales sur une partition ; Bayes P(Bj∣A)=P(Bj)P(A∣Bj)∑iP(Bi)P(A∣Bi)P(B_{j}\mid A)=\dfrac{P(B_{j})P(A\mid B_{j})}{\sum_{i}P(B_{i})P(A\mid B_{i})}. Indépendants : P(A∩B)=P(A)P(B)P(A\cap B)=P(A)P(B) ; incompatibles : P(A∩B)=0P(A\cap B)=0.
  • • VARIABLES ALÉATOIRES : E(X)=∑xipiE(X)=\sum x_{i}p_{i} ou ∫xf(x) dx\int x f(x)\,dx ; V(X)=E(X2)−E(X)2V(X)=E(X^{2})-E(X)^{2} ; E(aX+b)=aE(X)+bE(aX+b)=aE(X)+b, V(aX+b)=a2V(X)V(aX+b)=a^{2}V(X) ; une densité est positive d'intégrale 11, F(x)=∫−∞xf(t) dtF(x)=\int_{-\infty}^{x}f(t)\,dt.
  • • LOIS : binomiale P(X=k)=(nk)pk(1−p)n−kP(X=k)=\binom{n}{k}p^{k}(1-p)^{n-k}, E=npE=np, V=np(1−p)V=np(1-p) ; normale z=x−μσz=\dfrac{x-\mu}{\sigma}, table de Φ\Phi avec zz au centième ; approximation si np≥5np\ge 5 et n(1−p)≥5n(1-p)\ge 5, avec correction de continuité.
  • • INFÉRENCE : x‾±zsn\overline{x}\pm z\dfrac{s}{\sqrt{n}} (n≥30n\ge 30) ; x‾±tα/2 ;n−1sn\overline{x}\pm t_{\alpha/2\,;n-1}\dfrac{s}{\sqrt{n}} (n<30n<30, population normale) ; p^±zp^(1−p^)n\hat{p}\pm z\sqrt{\dfrac{\hat{p}(1-\hat{p})}{n}}. Test : zz ou t=x‾−μ0s/nt=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}}, z=p^−p0p0(1−p0)/nz=\dfrac{\hat{p}-p_{0}}{\sqrt{p_{0}(1-p_{0})/n}} ; H1H_{1} choisie avant les données.
  • • DEUX VARIABLES : a=SxySxxa=\dfrac{S_{xy}}{S_{xx}}, b=y‾−ax‾b=\overline{y}-a\overline{x}, r=SxySxxSyyr=\dfrac{S_{xy}}{\sqrt{S_{xx}S_{yy}}} ; khi-carré ∑(O−E)2E\sum\dfrac{(O-E)^{2}}{E}, E=ligne×colonnen≥5E=\dfrac{\text{ligne}\times\text{colonne}}{n}\ge 5, (l−1)(c−1)(l-1)(c-1) degrés de liberté.

Partie A : statistique descriptive (/8)

Exercice 1 : Le diamètre des épinettes noires d'une placette : données groupées et position

9 réponses vérifiables

8 points Statistique descriptive : tendance centrale, dispersion et position

Lors d'un inventaire forestier, une équipe délimite au hasard une placette de 400400 m² dans une pessière de 300300 hectares et mesure le diamètre à hauteur de poitrine (DHP, pris à 1,31{,}3 m du sol) de chacune des 5050 épinettes noires qui s'y trouvent. Les données, fictives, sont groupées en classes de 44 cm dans l'histogramme ci-dessous ; on traite ces 5050 arbres comme un échantillon des épinettes noires de la pessière.

Conventions : diviseur n−1n-1 pour l'écart-type d'un échantillon ; médiane et centiles par interpolation dans leur classe, Pk=L+k100n−Fnc×aP_{k}=L+\dfrac{\frac{k}{100}n-F}{n_{c}}\times a, où LL est la borne inférieure de la classe, FF l'effectif cumulé avant elle, ncn_{c} son effectif et aa son amplitude.

0481216613159528121620242832DHP (cm)effectif
  • a) Identifiez la population, l'échantillon et la nature de la variable étudiée.
  • b) Calculez la moyenne et l'écart-type du DHP à partir des centres de classes.
  • c) Déterminez la classe modale, puis calculez la médiane par interpolation.
  • d) Comparez la moyenne et la médiane, puis décrivez la forme de la distribution à l'aide de l'histogramme.
  • e) Calculez le 90e centile P90P_{90} et interprétez-le.
  • f) Une épinette de la placette mesure 29,029{,}0 cm. Calculez sa cote zz, puis dites si elle est inhabituelle selon le repère ∣z∣>2|z|>2.

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

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

Réponses

  • a) Population : les épinettes noires de la pessière ; échantillon : les 5050 arbres ; variable quantitative continue
  • b) x‾=18,0\overline{x}=18{,}0 cm, s=5,24s=5{,}24 cm
  • c) Classe modale [16 ;20[[16\,;20[ ; Me=17,6Me=17{,}6 cm
  • d) x‾>Me\overline{x}>Me : légère asymétrie à droite
  • e) P90=25,6P_{90}=25{,}6 cm
  • f) z=2,10z=2{,}10 : inhabituelle selon le repère
a)

Population : les épinettes noires de la pessière, celles sur lesquelles porte l'inventaire.

Échantillon : les 5050 épinettes de la placette.

Variable : le DHP, quantitative CONTINUE : c'est une longueur, qui peut prendre toute valeur d'un intervalle ; le ruban qui arrondit au millimètre ne change pas sa nature.

Nuance : une placette est un échantillon en grappe ; des arbres voisins se ressemblent plus que des arbres tirés un par un dans toute la pessière, et l'estimation est moins précise que ne le laisserait croire n=50n=50.

b)

Centres de classes : 1010, 1414, 1818, 2222, 2626, 3030.

∑nici=6×10+13×14+15×18+9×22+5×26+2×30=900\displaystyle \displaystyle\sum n_{i}c_{i}=6\times 10+13\times 14+15\times 18+9\times 22+5\times 26+2\times 30=900

x‾=90050=18,0\overline{x}=\dfrac{900}{50}=18{,}0 cm

∑ni(ci−x‾)2=6×64+13×16+15×0+9×16+5×64+2×144=1344\displaystyle \displaystyle\sum n_{i}(c_{i}-\overline{x})^{2}=6\times 64+13\times 16+15\times 0+9\times 16+5\times 64+2\times 144=1344

s2=134449=27,43s^{2}=\dfrac{1344}{49}=27{,}43, donc s=5,24s=5{,}24 cm.

Piège : diviser par n=50n=50 donne 5,185{,}18 cm, l'écart-type d'une population ; ces 5050 arbres sont un échantillon.

c)

La classe modale est [16 ;20[[16\,;20[, la plus haute de l'histogramme (1515 arbres).

Effectifs cumulés : 66 ; 1919 ; 3434. La donnée de rang n2=25\frac{n}{2}=25 tombe dans [16 ;20[[16\,;20[, avec F=19F=19 et nc=15n_{c}=15.

Me=16+25−1915×4=17,6Me=16+\dfrac{25-19}{15}\times 4=17{,}6 cm

d)

x‾=18,0\overline{x}=18{,}0 cm dépasse un peu Me=17,6Me=17{,}6 cm.

L'histogramme s'étire davantage à droite : du centre de la classe modale (1818 cm), il y a 1414 cm jusqu'à la borne 3232 cm, contre 1010 cm jusqu'à la borne 88 cm, et les deux grands arbres de [28 ;32[[28\,;32[ tirent la moyenne vers le haut.

La distribution est légèrement asymétrique à droite. L'écart entre moyenne et médiane, 0,40{,}4 cm, reste petit devant s=5,24s=5{,}24 cm : l'asymétrie est faible.

e)

90100×50=45\frac{90}{100}\times 50=45 ; les effectifs cumulés valent 3434 avant [20 ;24[[20\,;24[ et 4343 avant [24 ;28[[24\,;28[ : le 45e arbre est dans [24 ;28[[24\,;28[, avec F=43F=43 et nc=5n_{c}=5.

P90=24+45−435×4=25,6P_{90}=24+\dfrac{45-43}{5}\times 4=25{,}6 cm

Environ 90 %90\ \% des épinettes de l'échantillon ont un DHP inférieur à 25,625{,}6 cm.

f)
z=29,0−18,05,24=2,10\displaystyle z=\dfrac{29{,}0-18{,}0}{5{,}24}=2{,}10

∣z∣>2|z|>2 : l'arbre est inhabituel selon ce repère, l'un des plus gros de la placette.

Nuance : le repère ∣z∣>2|z|>2 vient des distributions en cloche ; ici, par interpolation dans [28 ;32[[28\,;32[, environ 48,548{,}5 arbres sur 5050, soit 97 %97\ \%, ont un DHP inférieur à 29,029{,}0 cm, ce qui confirme une position extrême sans en faire une erreur de mesure.

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 B : probabilités (/43)

Exercice 2 : Des pièges photographiques : choisir les sites, compter les séquences

6 réponses vérifiables

8 points Dénombrement : permutations, arrangements, combinaisons, Pascal et binôme

Pour suivre le lynx du Canada, une équipe de biologistes dispose de 55 caméras à déclenchement automatique et a repéré 1212 sites candidats dans sa zone d'étude, numérotés de 11 à 1212 sur la carte ci-dessous. Quatre d'entre eux, les sites 77, 88, 1111 et 1212, sont en milieu humide. On installe au plus une caméra par site.

123456789101112milieu humidezone d'étude
  • a) Les cinq caméras sont identiques. Calculez le nombre de façons de choisir les sites équipés, en justifiant l'outil de dénombrement choisi.
  • b) Les caméras sont maintenant de cinq modèles différents, et l'on note quel modèle va sur quel site. Calculez le nombre d'installations possibles.
  • c) Avec des caméras identiques, calculez le nombre de choix de cinq sites qui comptent au moins un site en milieu humide.
  • d) Calculez le nombre de choix de cinq sites qui comptent exactement deux sites en milieu humide.
  • e) Montrez, en classant les choix selon leur nombre de sites humides, que ∑k=04(4k)(85−k)=(125)\displaystyle\sum_{k=0}^{4}\binom{4}{k}\binom{8}{5-k}=\binom{12}{5}.
  • f) Les cinq sites sont tirés au hasard parmi les douze, tous les choix ayant la même probabilité. Calculez la probabilité que deux sites exactement soient en milieu humide.
  • g) Une caméra enregistre six déclenchements dans une nuit : trois cerfs de Virginie, deux renards roux et un lynx. Calculez le nombre de séquences d'espèces possibles, dans l'ordre des déclenchements.

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

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

Réponses

  • a) (125)=792\binom{12}{5}=792 (combinaison)
  • b) A(12,5)=95 040A(12,5)=95\,040
  • c) 792−56=736792-56=736
  • d) (42)(83)=336\binom{4}{2}\binom{8}{3}=336
  • e) 56+280+336+112+8=79256+280+336+112+8=792 : partition selon le nombre kk de sites humides
  • f) P=1433≈0,4242P=\frac{14}{33}\approx 0{,}4242
  • g) 6!3! 2! 1!=60\frac{6!}{3!\,2!\,1!}=60
a)

On choisit un ENSEMBLE de 55 sites parmi 1212 : l'ordre ne compte pas (les caméras sont identiques) et il n'y a pas de répétition (une caméra au plus par site). C'est une combinaison.

(125)=12!5! 7!=12×11×10×9×8120=792\displaystyle \dbinom{12}{5}=\dfrac{12!}{5!\,7!}=\dfrac{12\times 11\times 10\times 9\times 8}{120}=792
b)

Le modèle 1 a 1212 sites possibles, le modèle 2 en a 1111, et ainsi de suite : l'ordre compte maintenant, c'est un arrangement.

A(12,5)=12×11×10×9×8=95 040=5!×792\displaystyle A(12,5)=12\times 11\times 10\times 9\times 8=95\,040=5!\times 792

Chaque ensemble de 55 sites donne 5!=1205!=120 façons de distribuer les modèles.

c)

On passe par le contraire : « aucun site humide » revient à choisir les 55 sites parmi les 88 autres.

(125)−(85)=792−56=736\displaystyle \dbinom{12}{5}-\dbinom{8}{5}=792-56=736

Piège : compter « un site humide, puis 44 sites quelconques », 4×(114)=13204\times\binom{11}{4}=1320, compte plusieurs fois un même choix qui a deux sites humides ou plus.

d)

Deux sites parmi les 44 humides ET trois parmi les 88 autres : principe de multiplication.

(42)(83)=6×56=336\displaystyle \dbinom{4}{2}\dbinom{8}{3}=6\times 56=336
e)

Tout choix de 55 sites contient exactement kk sites humides, pour un seul kk entre 00 et 44 : les cinq cas forment une partition des (125)\binom{12}{5} choix.

Pour un kk fixé, on choisit kk sites humides et 5−k5-k autres : (4k)(85−k)\binom{4}{k}\binom{8}{5-k} choix (principe de multiplication).

Principe d'addition : 56+280+336+112+8=792=(125)56+280+336+112+8=792=\dbinom{12}{5}.

f)

Les 792792 choix sont équiprobables, dont 336336 favorables.

P=336792=1433≈0,4242\displaystyle P=\dfrac{336}{792}=\dfrac{14}{33}\approx 0{,}4242
g)

Une séquence est une permutation de 66 objets dont certains sont identiques : 33 cerfs, 22 renards, 11 lynx.

6!3! 2! 1!=72012=60\displaystyle \dfrac{6!}{3!\,2!\,1!}=\dfrac{720}{12}=60

Piège : 6!=7206!=720 compterait comme différentes deux séquences qui ne diffèrent que par l'échange de deux cerfs, indiscernables sur la liste des espèces.

Exercice 3 : Un logiciel d'identification des félins : Bayes et indépendance conditionnelle

8 réponses vérifiables

9 points Probabilités et probabilités conditionnelles : axiomes, indépendance, Bayes

Les photos de félins prises par les caméras de l'exercice 2 sont triées par un logiciel. Sur une photo de félin prise au hasard, l'animal est un lynx roux (RR) avec la probabilité 0,550{,}55, un lynx du Canada (CC) avec la probabilité 0,300{,}30 et un chat domestique errant (DD) avec la probabilité 0,150{,}15. On note LL l'événement « le logiciel affiche l'étiquette lynx du Canada ». L'arbre donne P(L∣R)P(L\mid R), P(L∣C)P(L\mid C) et P(L∣D)P(L\mid D), estimées sur des photos déjà identifiées par des experts. Données fictives.

RLLCLLDLL0,550,080,920,300,900,100,150,020,98
  • a) Démontrez, à partir de la définition d'une probabilité conditionnelle, que P(C∣L)=P(C)P(L∣C)P(R)P(L∣R)+P(C)P(L∣C)+P(D)P(L∣D)P(C\mid L)=\dfrac{P(C)P(L\mid C)}{P(R)P(L\mid R)+P(C)P(L\mid C)+P(D)P(L\mid D)}.
  • b) Calculez P(L)P(L).
  • c) Calculez P(C∣L)P(C\mid L) et P(R∣L)P(R\mid L).
  • d) Les événements CC et RR sont-ils incompatibles ? Sont-ils indépendants ? Justifiez.
  • e) Le même animal est photographié deux fois de suite. On admet que, l'espèce étant fixée, les deux étiquettes sont indépendantes. Calculez la probabilité que l'animal soit un lynx du Canada sachant que les deux photos portent l'étiquette LL.
  • f) Sans connaître l'espèce, les événements L1L_{1} « la première photo est étiquetée LL » et L2L_{2} « la seconde photo est étiquetée LL » sont-ils indépendants ? Justifiez par un calcul.

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

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

Réponses

  • a) Définition, règle de multiplication et probabilités totales sur la partition RR, CC, DD
  • b) P(L)=0,317P(L)=0{,}317
  • c) P(C∣L)≈0,8517P(C\mid L)\approx 0{,}8517, P(R∣L)≈0,1388P(R\mid L)\approx 0{,}1388
  • d) Incompatibles, donc dépendants : P(C∩R)=0≠0,165P(C\cap R)=0\neq 0{,}165
  • e) P(C∣L1∩L2)≈0,9855P(C\mid L_{1}\cap L_{2})\approx 0{,}9855
  • f) P(L1∩L2)≈0,2466≠0,3172≈0,1005P(L_{1}\cap L_{2})\approx 0{,}2466\neq 0{,}317^{2}\approx 0{,}1005 : non indépendants
a)

Définition : P(C∣L)=P(C∩L)P(L)P(C\mid L)=\dfrac{P(C\cap L)}{P(L)}, et de même P(C∩L)=P(C)P(L∣C)P(C\cap L)=P(C)P(L\mid C) (règle de multiplication).

RR, CC, DD forment une partition de l'univers : ils sont deux à deux incompatibles et leur réunion est certaine. Donc L=(L∩R)∪(L∩C)∪(L∩D)L=(L\cap R)\cup(L\cap C)\cup(L\cap D), réunion d'événements incompatibles.

Probabilités totales : P(L)=P(R)P(L∣R)+P(C)P(L∣C)+P(D)P(L∣D)P(L)=P(R)P(L\mid R)+P(C)P(L\mid C)+P(D)P(L\mid D).

En remplaçant le numérateur et le dénominateur, on obtient la formule annoncée, qui suppose P(L)>0P(L)>0.

b)
P(L)=0,55×0,08+0,30×0,90+0,15×0,02\displaystyle P(L)=0{,}55\times 0{,}08+0{,}30\times 0{,}90+0{,}15\times 0{,}02
P(L)=0,044+0,270+0,003=0,317\displaystyle P(L)=0{,}044+0{,}270+0{,}003=0{,}317
c)

P(C∣L)=0,2700,317≈0,8517P(C\mid L)=\dfrac{0{,}270}{0{,}317}\approx 0{,}8517 et P(R∣L)=0,0440,317≈0,1388P(R\mid L)=\dfrac{0{,}044}{0{,}317}\approx 0{,}1388

Près de 14 %14\ \% des photos étiquetées « lynx du Canada » montrent en fait un lynx roux : le logiciel se trompe rarement sur un lynx roux (8 %8\ \%), mais les lynx roux sont presque deux fois plus nombreux.

Piège : confondre P(L∣C)=0,90P(L\mid C)=0{,}90 et P(C∣L)=0,85P(C\mid L)=0{,}85.

d)

Incompatibles : oui, un animal n'est pas des deux espèces à la fois, P(C∩R)=0P(C\cap R)=0.

Indépendants : non, car P(C)P(R)=0,30×0,55=0,165≠0=P(C∩R)P(C)P(R)=0{,}30\times 0{,}55=0{,}165\neq 0=P(C\cap R).

Deux événements incompatibles de probabilités non nulles sont toujours dépendants : savoir que c'est un lynx roux EXCLUT le lynx du Canada.

e)

Pour chaque espèce, la probabilité de deux étiquettes LL est le carré de P(L∣⋅)P(L\mid\cdot) (indépendance sachant l'espèce).

P(L1∩L2)=0,55×0,082+0,30×0,902+0,15×0,022=0,00352+0,243+0,00006=0,24658\displaystyle P(L_{1}\cap L_{2})=0{,}55\times 0{,}08^{2}+0{,}30\times 0{,}90^{2}+0{,}15\times 0{,}02^{2}=0{,}00352+0{,}243+0{,}00006=0{,}24658
P(C∣L1∩L2)=0,2430,24658≈0,9855\displaystyle P(C\mid L_{1}\cap L_{2})=\dfrac{0{,}243}{0{,}24658}\approx 0{,}9855

La seconde étiquette fait passer la probabilité de 0,850{,}85 à 0,990{,}99 : Bayes appliqué deux fois, la loi a posteriori de la première photo servant de loi a priori à la seconde.

Nuance : deux photos successives du même animal partagent l'angle, la lumière et la distance ; leurs erreurs sont sans doute liées, et 0,98550{,}9855 est une borne optimiste, valable dans le modèle admis.

f)

P(L1∩L2)=0,24658P(L_{1}\cap L_{2})=0{,}24658 (question e) et P(L1)P(L2)=0,3172≈0,1005P(L_{1})P(L_{2})=0{,}317^{2}\approx 0{,}1005.

Les deux nombres diffèrent : L1L_{1} et L2L_{2} ne sont PAS indépendants.

Une première étiquette LL rend plus probable un lynx du Canada, donc une seconde étiquette LL. L'indépendance sachant l'espèce n'entraîne pas l'indépendance tout court.

Exercice 4 : Analyser les écouvillons par lots : une variable aléatoire à deux valeurs

9 réponses vérifiables

8 points Variables aléatoires discrètes et continues : loi, densité, espérance et varianceProbabilités et probabilités conditionnelles : axiomes, indépendance, Bayes

Pour surveiller l'influenza aviaire chez les canards sauvages, un laboratoire reçoit des écouvillons prélevés sur des oiseaux. Pour économiser des analyses, il les regroupe par lots de k=5k=5 : il mélange les cinq écouvillons et fait une seule analyse du mélange. Si le mélange est négatif, les cinq écouvillons sont déclarés négatifs. S'il est positif, on analyse ensuite chacun des cinq écouvillons séparément, soit six analyses en tout.

Modèle adopté : chaque écouvillon est positif avec la probabilité p=0,05p=0{,}05, indépendamment des autres, et l'analyse du mélange est positive si et seulement si au moins un écouvillon l'est. On note XX le nombre d'analyses faites pour un lot.

  • a) Donnez la loi de probabilité de XX.
  • b) Calculez E(X)E(X), puis le nombre moyen d'analyses par écouvillon, et comparez à une analyse par écouvillon.
  • c) Calculez V(X)V(X) à l'aide de V(X)=E(X2)−E(X)2V(X)=E(X^{2})-E(X)^{2}, puis l'écart-type σ(X)\sigma(X).
  • d) Démontrez que, pour toute variable aléatoire discrète XX et tous réels aa et bb, E(aX+b)=aE(X)+bE(aX+b)=aE(X)+b et V(aX+b)=a2V(X)V(aX+b)=a^{2}V(X).
  • e) Chaque analyse coûte 3030 dollars et la préparation d'un lot 1212 dollars. Exprimez le coût CC d'un lot en fonction de XX, puis calculez E(C)E(C) et σ(C)\sigma(C).
  • f) Le laboratoire peut choisir la taille kk des lots. Montrez que le nombre moyen d'analyses par écouvillon vaut g(k)=1k+1−0,95kg(k)=\dfrac{1}{k}+1-0{,}95^{k}, puis déterminez l'entier kk entre 22 et 1010 qui le rend minimal.

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

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

Réponses

  • a) P(X=1)=0,955≈0,7738P(X=1)=0{,}95^{5}\approx 0{,}7738, P(X=6)≈0,2262P(X=6)\approx 0{,}2262
  • b) E(X)≈2,1311E(X)\approx 2{,}1311 ; 0,42620{,}4262 analyse par écouvillon
  • c) V(X)≈4,376V(X)\approx 4{,}376, σ(X)≈2,092\sigma(X)\approx 2{,}092
  • d) Linéarité de la somme et ∑pi=1\sum p_{i}=1 ; bb s'annule dans les écarts
  • e) C=30X+12C=30X+12 : E(C)≈75,93E(C)\approx 75{,}93 dollars, σ(C)≈62,76\sigma(C)\approx 62{,}76 dollars
  • f) k=5k=5, avec g(5)≈0,4262g(5)\approx 0{,}4262
a)

XX ne prend que deux valeurs : 11 (mélange négatif) ou 66 (mélange positif).

Le mélange est négatif si les cinq écouvillons le sont : par indépendance, P(X=1)=0,955≈0,7738P(X=1)=0{,}95^{5}\approx 0{,}7738.

P(X=6)=1−0,955≈0,2262P(X=6)=1-0{,}95^{5}\approx 0{,}2262 ; la somme vaut bien 11.

b)

E(X)=1×0,7738+6×0,2262=1+5×0,2262≈2,1311E(X)=1\times 0{,}7738+6\times 0{,}2262=1+5\times 0{,}2262\approx 2{,}1311 analyses par lot.

Par écouvillon : 2,13115≈0,4262\dfrac{2{,}1311}{5}\approx 0{,}4262 analyse, au lieu de 11.

Le regroupement économise en moyenne environ 57 %57\ \% des analyses, dans ce modèle.

c)
E(X2)=12×0,7738+62×0,2262≈8,9177\displaystyle E(X^{2})=1^{2}\times 0{,}7738+6^{2}\times 0{,}2262\approx 8{,}9177
V(X)=8,9177−2,13112≈4,3761\displaystyle V(X)=8{,}9177-2{,}1311^{2}\approx 4{,}3761

σ(X)=4,3761≈2,092\sigma(X)=\sqrt{4{,}3761}\approx 2{,}092 analyses.

Vérification : une variable à deux valeurs distantes de 55 a pour variance 52×0,7738×0,2262≈4,37615^{2}\times 0{,}7738\times 0{,}2262\approx 4{,}3761.

d)

Soit x1,…,xnx_{1},\dots,x_{n} les valeurs de XX et p1,…,pnp_{1},\dots,p_{n} leurs probabilités, avec ∑pi=1\sum p_{i}=1. aX+baX+b prend les valeurs axi+bax_{i}+b avec les mêmes probabilités.

E(aX+b)=∑(axi+b)pi=a∑xipi+b∑pi=aE(X)+b\begin{aligned}E(aX+b)&=\sum(ax_{i}+b)p_{i}\\&=a\sum x_{i}p_{i}+b\sum p_{i}\\&=aE(X)+b\end{aligned}
V(aX+b)=∑(axi+b−aE(X)−b)2pi=a2∑(xi−E(X))2pi=a2V(X)\begin{aligned}V(aX+b)&=\sum\left(ax_{i}+b-aE(X)-b\right)^{2}p_{i}\\&=a^{2}\sum(x_{i}-E(X))^{2}p_{i}\\&=a^{2}V(X)\end{aligned}

La constante bb décale toutes les valeurs sans les écarter : elle disparaît de la variance.

e)

C=30X+12C=30X+12, en dollars.

E(C)=30×2,1311+12≈75,93E(C)=30\times 2{,}1311+12\approx 75{,}93 dollars

σ(C)=∣30∣ σ(X)≈30×2,0919≈62,76\sigma(C)=|30|\,\sigma(X)\approx 30\times 2{,}0919\approx 62{,}76 dollars

Piège : ajouter 1212 à l'écart-type ; un coût fixe ne change pas la dispersion.

f)

Pour un lot de kk écouvillons : X=1X=1 avec la probabilité 0,95k0{,}95^{k}, X=k+1X=k+1 sinon.

E(X)=0,95k+(k+1)(1−0,95k)=1+k(1−0,95k)E(X)=0{,}95^{k}+(k+1)(1-0{,}95^{k})=1+k(1-0{,}95^{k}), et par écouvillon g(k)=E(X)k=1k+1−0,95kg(k)=\dfrac{E(X)}{k}=\dfrac{1}{k}+1-0{,}95^{k}.

Valeurs : g(4)≈0,4355g(4)\approx 0{,}4355 ; g(5)≈0,4262g(5)\approx 0{,}4262 ; g(6)≈0,4316g(6)\approx 0{,}4316 ; gg décroît de k=2k=2 à 55 puis croît jusqu'à 1010.

Le minimum est atteint pour k=5k=5 : de petits lots demandent beaucoup d'analyses de mélange, de grands lots sont trop souvent positifs.

Nuance : le modèle suppose les écouvillons indépendants et le mélange sans perte de sensibilité ; des oiseaux d'un même étang se contaminent entre eux, et la dilution d'un écouvillon faiblement positif peut le masquer.

Exercice 5 : La superficie des feux de foudre : une densité sur un intervalle infini

9 réponses vérifiables

9 points Variables aléatoires discrètes et continues : loi, densité, espérance et varianceLoi binomiale et loi normale : conditions, table, approximation

Dans une zone de protection intensive, les feux de forêt allumés par la foudre sont repérés et éteints vite, si bien que la plupart restent petits. Une équipe modélise la superficie brûlée XX d'un tel feu, en hectares, par une variable aléatoire continue de densité f(x)=k(1+x)4f(x)=\dfrac{k}{(1+x)^{4}} pour x≥0x\ge 0 et f(x)=0f(x)=0 pour x<0x<0, où kk est une constante. La figure représente ff pour la valeur de kk établie en a).

00,511,522,53012345superficie x (ha)f(x) = 3/(1 + x)⁴
  • a) Montrez que k=3k=3.
  • b) Déterminez la fonction de répartition F(x)=P(X≤x)F(x)=P(X\le x) pour x≥0x\ge 0, puis calculez P(X>1)P(X>1) et P(0,5≤X≤2)P(0{,}5\le X\le 2).
  • c) Calculez la médiane de XX.
  • d) Calculez l'espérance E(X)E(X) par une intégration par parties.
  • e) Calculez E(X2)E(X^{2}), puis la variance et l'écart-type de XX. Comparez l'espérance et la médiane, et interprétez.
  • f) Un feu a déjà brûlé plus de 11 hectare quand une équipe arrive. Calculez la probabilité qu'il dépasse 22 hectares.
  • g) Pendant une saison, 1212 feux de foudre éclatent dans la zone, et l'on admet leurs superficies indépendantes. Calculez la probabilité qu'au moins 33 d'entre eux dépassent 11 hectare.

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) ∫0+∞k(1+x)4dx=k3=1\int_{0}^{+\infty}\frac{k}{(1+x)^{4}}dx=\frac{k}{3}=1, donc k=3k=3
  • b) F(x)=1−(1+x)−3F(x)=1-(1+x)^{-3} ; P(X>1)=0,125P(X>1)=0{,}125 ; P(0,5≤X≤2)≈0,2593P(0{,}5\le X\le 2)\approx 0{,}2593
  • c) m=23−1≈0,260m=\sqrt[3]{2}-1\approx 0{,}260 ha
  • d) E(X)=0,5E(X)=0{,}5 ha
  • e) E(X2)=1E(X^{2})=1, V(X)=0,75V(X)=0{,}75, σ≈0,866\sigma\approx 0{,}866 ha ; E(X)≈2 mE(X)\approx 2\,m : forte asymétrie à droite
  • f) 827≈0,2963\frac{8}{27}\approx 0{,}2963
  • g) Y∼B(12 ;0,125)Y\sim\mathcal{B}(12\,;0{,}125), P(Y≥3)≈0,1820P(Y\ge 3)\approx 0{,}1820
a)

Une densité est positive et d'intégrale totale 11 : f≥0f\ge 0 exige k≥0k\ge 0, puis

∫0+∞k(1+x)4 dx=lim⁡b→+∞[−k3(1+x)3]0b=lim⁡b→+∞(k3−k3(1+b)3)=k3\displaystyle \displaystyle\int_{0}^{+\infty}\frac{k}{(1+x)^{4}}\,dx=\lim_{b\to+\infty}\left[-\frac{k}{3(1+x)^{3}}\right]_{0}^{b}=\lim_{b\to+\infty}\left(\frac{k}{3}-\frac{k}{3(1+b)^{3}}\right)=\frac{k}{3}

k3=1\dfrac{k}{3}=1 donne k=3k=3.

b)

F(x)=∫0x3(1+t)4 dt=[−1(1+t)3]0x=1−1(1+x)3F(x)=\displaystyle\int_{0}^{x}\frac{3}{(1+t)^{4}}\,dt=\left[-\frac{1}{(1+t)^{3}}\right]_{0}^{x}=1-\frac{1}{(1+x)^{3}} pour x≥0x\ge 0 (F(x)=0F(x)=0 pour x<0x<0).

P(X>1)=1−F(1)=123=0,125\displaystyle P(X>1)=1-F(1)=\dfrac{1}{2^{3}}=0{,}125
P(0,5≤X≤2)=F(2)−F(0,5)=11,53−133≈0,2963−0,0370=0,2593\displaystyle P(0{,}5\le X\le 2)=F(2)-F(0{,}5)=\dfrac{1}{1{,}5^{3}}-\dfrac{1}{3^{3}}\approx 0{,}2963-0{,}0370=0{,}2593

Pour une variable continue, P(X=0,5)=0P(X=0{,}5)=0 : inclure ou non les bornes ne change rien.

c)

La médiane mm vérifie F(m)=12F(m)=\frac{1}{2} : (1+m)3=2(1+m)^{3}=2, donc m=23−1≈0,260m=\sqrt[3]{2}-1\approx 0{,}260 hectare.

La moitié des feux brûlent moins de 0,260{,}26 ha, soit environ 2 6002\,600 m².

d)

E(X)=∫0+∞x 3(1+x)4 dxE(X)=\displaystyle\int_{0}^{+\infty}x\,\frac{3}{(1+x)^{4}}\,dx. Par parties, u=xu=x et v′=3(1+x)4v'=\dfrac{3}{(1+x)^{4}}, donc u′=1u'=1 et v=−1(1+x)3v=-\dfrac{1}{(1+x)^{3}} :

E(X)=lim⁡b→+∞[−x(1+x)3]0b+∫0+∞dx(1+x)3=0+[−12(1+x)2]0+∞=12\displaystyle \displaystyle E(X)=\lim_{b\to+\infty}\left[-\frac{x}{(1+x)^{3}}\right]_{0}^{b}+\int_{0}^{+\infty}\frac{dx}{(1+x)^{3}}=0+\left[-\frac{1}{2(1+x)^{2}}\right]_{0}^{+\infty}=\frac{1}{2}

Le crochet tend vers 00 car b(1+b)3→0\dfrac{b}{(1+b)^{3}}\to 0. E(X)=0,5E(X)=0{,}5 hectare.

e)

Par parties encore, u=x2u=x^{2} et v=−1(1+x)3v=-\dfrac{1}{(1+x)^{3}} ; le crochet tend vers 00 :

E(X2)=∫0+∞2x(1+x)3 dx\displaystyle E(X^{2})=\int_{0}^{+\infty}\frac{2x}{(1+x)^{3}}\,dx, et avec w=1+xw=1+x : 2∫1+∞w−1w3 dw=2[−1w+12w2]1+∞=2×12=1\displaystyle 2\int_{1}^{+\infty}\frac{w-1}{w^{3}}\,dw=2\left[-\frac{1}{w}+\frac{1}{2w^{2}}\right]_{1}^{+\infty}=2\times\frac{1}{2}=1

V(X)=E(X2)−E(X)2=1−0,25=0,75V(X)=E(X^{2})-E(X)^{2}=1-0{,}25=0{,}75 et σ(X)=0,75≈0,866\sigma(X)=\sqrt{0{,}75}\approx 0{,}866 hectare.

E(X)=0,5E(X)=0{,}5 ha vaut presque le double de la médiane 0,260{,}26 ha : la distribution est très asymétrique à droite. Quelques grands feux pèsent lourd dans la moyenne, alors que le feu typique reste petit.

Nuance : la densité donne une probabilité non nulle à des feux immenses ; dans la réalité, la superficie est bornée par la zone, et le modèle ne vaut que pour l'ordre de grandeur des feux observés.

f)
P(X>2∣X>1)=P(X>2)P(X>1)=1/271/8=827≈0,2963\displaystyle P(X>2\mid X>1)=\dfrac{P(X>2)}{P(X>1)}=\dfrac{1/27}{1/8}=\dfrac{8}{27}\approx 0{,}2963

L'événement X>2X>2 est inclus dans X>1X>1, d'où l'intersection égale à X>2X>2.

Sans condition, P(X>2)≈0,0370P(X>2)\approx 0{,}0370 : un feu qui a déjà dépassé 11 ha a huit fois plus de chances de dépasser 22 ha qu'un feu qui commence.

g)

Chaque feu dépasse 11 ha avec la probabilité q=0,125q=0{,}125 (question b), indépendamment des autres : le nombre YY de tels feux suit la loi binomiale B(12 ;0,125)\mathcal{B}(12\,;0{,}125).

P(Y=0)=0,87512≈0,2014P(Y=0)=0{,}875^{12}\approx 0{,}2014 ; P(Y=1)=12×0,125×0,87511≈0,3453P(Y=1)=12\times 0{,}125\times 0{,}875^{11}\approx 0{,}3453 ; P(Y=2)=(122)×0,1252×0,87510≈0,2713P(Y=2)=\dbinom{12}{2}\times 0{,}125^{2}\times 0{,}875^{10}\approx 0{,}2713

P(Y≥3)=1−(0,2014+0,3453+0,2713)≈0,1820\displaystyle P(Y\ge 3)=1-(0{,}2014+0{,}3453+0{,}2713)\approx 0{,}1820

Nuance : des feux d'un même orage, par temps sec et venteux, ont des superficies liées ; l'indépendance est une hypothèse du modèle.

00,250,50,751012345médiane : 0,26 haF(x) = 1 − (1 + x)⁻³superficie x (ha)
La fonction de répartition F(x)=1−(1+x)−3F(x)=1-(1+x)^{-3} : la médiane se lit où FF vaut 0,50{,}5, vers 0,260{,}26 ha.

Exercice 6 : Des saumoneaux marqués franchissent un barrage : de la loi normale à la binomiale

8 réponses vérifiables

9 points Loi binomiale et loi normale : conditions, table, approximation

Sur une rivière à saumon équipée d'un barrage hydroélectrique, des biologistes suivent la dévalaison des saumoneaux, les jeunes saumons atlantiques qui descendent vers la mer. On modélise la longueur LL d'un saumoneau par une loi normale de moyenne μ=14,6\mu=14{,}6 cm et d'écart-type σ=1,1\sigma=1{,}1 cm (valeurs fictives). Un émetteur acoustique ne peut être implanté que dans un poisson d'au moins 1313 cm.

Pour la loi normale, arrondissez zz au centième et lisez la table de Φ(z)=P(Z≤z)\Phi(z)=P(Z\le z). Critère de ce cours pour l'approximation normale d'une binomiale : np≥5np\ge 5 et n(1−p)≥5n(1-p)\ge 5.

111213141516171813 cmμ = 14,6longueur L (cm)
  • a) Calculez la proportion des saumoneaux qui peuvent recevoir un émetteur.
  • b) Les biologistes écartent les 10 %10\ \% de saumoneaux les plus petits. Calculez la longueur minimale retenue, sachant que la table donne Φ(−1,28)≈0,10\Phi(-1{,}28)\approx 0{,}10.
  • c) On admet que chaque saumoneau marqué franchit le barrage vivant avec la probabilité 0,930{,}93, indépendamment des autres. On marque 1010 saumoneaux et l'on note NN le nombre de survivants. Justifiez l'emploi d'une loi binomiale pour NN, puis calculez la probabilité qu'ils survivent tous, et celle qu'au moins 99 survivent.
  • d) On marque 400400 saumoneaux et l'on note XX le nombre de survivants. Calculez E(X)E(X) et σ(X)\sigma(X), puis vérifiez les conditions de l'approximation normale.
  • e) Calculez P(X≤360)P(X\le 360) par l'approximation normale, avec la correction de continuité.
  • f) Cette année, 360360 des 400400 saumoneaux marqués ont franchi le barrage. Ce résultat rend-il douteux le taux de survie de 0,930{,}93 ? Justifiez.

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) z≈−1,45z\approx -1{,}45 ; P(L≥13)≈0,9265P(L\ge 13)\approx 0{,}9265
  • b) ℓ≈13,19\ell\approx 13{,}19 cm
  • c) P(N=10)≈0,4840P(N=10)\approx 0{,}4840, P(N≥9)≈0,8483P(N\ge 9)\approx 0{,}8483
  • d) E(X)=372E(X)=372, σ(X)≈5,10\sigma(X)\approx 5{,}10 ; np=372np=372, n(1−p)=28n(1-p)=28
  • e) P(X≤360)≈Φ(−2,25)=0,0122P(X\le 360)\approx\Phi(-2{,}25)=0{,}0122
  • f) Oui : environ 1,2 %1{,}2\ \% sous le taux annoncé (exact : 1,55 %1{,}55\ \%), le taux devient douteux
a)
z=13−14,61,1≈−1,45\displaystyle z=\dfrac{13-14{,}6}{1{,}1}\approx -1{,}45
P(L≥13)=1−Φ(−1,45)=Φ(1,45)=0,9265\displaystyle P(L\ge 13)=1-\Phi(-1{,}45)=\Phi(1{,}45)=0{,}9265

Environ 92,7 %92{,}7\ \% des saumoneaux mesurent au moins 1313 cm.

b)

On cherche ℓ\ell tel que P(L<ℓ)=0,10P(L<\ell)=0{,}10, soit z=−1,28z=-1{,}28.

ℓ=μ+zσ=14,6−1,28×1,1≈13,19\ell=\mu+z\sigma=14{,}6-1{,}28\times 1{,}1\approx 13{,}19 cm

Piège : prendre z=+1,28z=+1{,}28, qui donne le 90e centile, 16,0116{,}01 cm.

c)

NN suit une loi binomiale B(10 ;0,93)\mathcal{B}(10\,;0{,}93) : n=10n=10 épreuves fixées, deux issues (survit ou non), même probabilité p=0,93p=0{,}93, épreuves indépendantes par hypothèse.

P(N=10)=0,9310≈0,4840\displaystyle P(N=10)=0{,}93^{10}\approx 0{,}4840
P(N≥9)=0,9310+10×0,07×0,939≈0,4840+0,3643=0,8483\displaystyle P(N\ge 9)=0{,}93^{10}+10\times 0{,}07\times 0{,}93^{9}\approx 0{,}4840+0{,}3643=0{,}8483

Nuance : les saumoneaux dévalent en groupes ; un prédateur posté à la sortie des turbines en frappe plusieurs à la fois, et l'indépendance reste une hypothèse de modèle.

d)
E(X)=np=400×0,93=372\displaystyle E(X)=np=400\times 0{,}93=372
σ(X)=np(1−p)=26,04≈5,10\displaystyle \sigma(X)=\sqrt{np(1-p)}=\sqrt{26{,}04}\approx 5{,}10

np=372≥5np=372\ge 5 et n(1−p)=28≥5n(1-p)=28\ge 5 : l'approximation par la loi normale N(372 ;5,102)N(372\,;5{,}10^{2}) est acceptée selon le critère du cours.

e)

XX est entière : X≤360X\le 360 correspond à Y≤360,5Y\le 360{,}5 pour la variable normale YY.

z=360,5−3725,10≈−2,25\displaystyle z=\dfrac{360{,}5-372}{5{,}10}\approx -2{,}25
P(X≤360)≈Φ(−2,25)=0,0122\displaystyle P(X\le 360)\approx\Phi(-2{,}25)=0{,}0122

Piège : sans correction, z=−2,35z=-2{,}35 et l'on trouve 0,00940{,}0094, qui néglige la moitié de la barre de 360360.

f)

Si le taux de survie valait 0,930{,}93, on n'observerait 360360 survivants ou moins qu'environ une année sur 8080.

Le résultat est donc peu compatible avec 0,930{,}93 : il rend ce taux douteux, sans dire pourquoi (année exceptionnelle, prédation, débit des turbines).

Remarque : le calcul exact par la loi binomiale donne P(X≤360)≈0,0155P(X\le 360)\approx 0{,}0155. Avec pp proche de 11, la loi est un peu asymétrique et l'approximation normale sous-estime cette queue ; la conclusion ne change pas.

C'est déjà le raisonnement d'un test d'hypothèse, que les exercices 7 à 9 rendent formel.

Partie C : inférence (/26)

Exercice 7 : Le radon des sous-sols : théorème central limite et test sur une moyenne

10 réponses vérifiables

9 points Intervalles de confiance : moyenne, Student et proportionTests d'hypothèse : moyenne, Student et proportionLoi binomiale et loi normale : conditions, table, approximation

Le radon, un gaz radioactif issu de l'uranium des sols et des roches, s'accumule dans les sous-sols. Une campagne régionale (fictive) a mesuré la concentration moyenne annuelle de radon dans 600600 sous-sols ; l'histogramme ci-dessous les regroupe en classes de 5050 Bq/m³. On admet que, pour l'ensemble des maisons de la région, la moyenne vaut μ=95\mu=95 Bq/m³ et l'écart-type σ=80\sigma=80 Bq/m³. Au Canada, la ligne directrice recommande d'agir dans une maison au-delà de 200200 Bq/m³.

Valeurs utiles : z0,05=1,645z_{0{,}05}=1{,}645 et z0,01=2,326z_{0{,}01}=2{,}326, qui laissent une aire de 0,050{,}05 et 0,010{,}01 à leur droite ; table de Φ\Phi avec zz arrondi au centième.

050100150200050100150200250300350400450radon (Bq/m³)maisons
  • a) Peut-on modéliser la concentration XX d'une maison par la loi normale de moyenne 9595 et d'écart-type 8080 ? Justifiez à l'aide de l'histogramme et du calcul de P(X<0)P(X<0) dans ce modèle.
  • b) Soit X1,…,XnX_{1},\dots,X_{n} les concentrations de nn maisons tirées au hasard avec remise, et X‾\overline{X} leur moyenne. Démontrez que E(X‾)=μE(\overline{X})=\mu et V(X‾)=σ2nV(\overline{X})=\dfrac{\sigma^{2}}{n}.
  • c) On tire 5050 maisons au hasard. Énoncez ce que le théorème central limite permet d'affirmer sur X‾\overline{X}, puis calculez P(X‾>120)P(\overline{X}>120).
  • d) Une municipalité bâtie sur un socle granitique veut savoir si la concentration moyenne de ses maisons dépasse la moyenne régionale. Sur 4040 maisons tirées au hasard, elle obtient x‾=118\overline{x}=118 Bq/m³ et s=86s=86 Bq/m³. Posez les hypothèses, puis calculez la statistique de test.
  • e) Concluez au seuil de 5 %5\ \%, puis au seuil de 1 %1\ \%, et calculez la valeur pp du test.
  • f) Le maire annonce : « Le test prouve que nos maisons dépassent la ligne directrice. » Choisissez la conclusion correcte.

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

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

Réponses

  • a) Non : distribution asymétrique, et P(X<0)≈0,1170P(X<0)\approx 0{,}1170 serait absurde
  • b) Linéarité de EE ; indépendance : V(∑Xi)=nσ2V(\sum X_{i})=n\sigma^{2}, d'où σ2n\frac{\sigma^{2}}{n}
  • c) X‾\overline{X} approximativement normale, σX‾≈11,31\sigma_{\overline{X}}\approx 11{,}31 ; P(X‾>120)≈0,0136P(\overline{X}>120)\approx 0{,}0136
  • d) H0:μ=95H_{0}:\mu=95, H1:μ>95H_{1}:\mu>95 ; z≈1,69z\approx 1{,}69
  • e) Rejet à 5 %5\ \%, non-rejet à 1 %1\ \% ; valeur p ≈0,0455\approx 0{,}0455
  • f) Le test compare la moyenne à 9595, il ne dit rien d'une maison ni du seuil de 200200
a)

Non. L'histogramme est fortement asymétrique à droite : un mode près de 00, une longue queue jusqu'à 450450 Bq/m³.

Dans le modèle normal : z=0−9580≈−1,19z=\dfrac{0-95}{80}\approx -1{,}19 et P(X<0)=Φ(−1,19)=0,1170P(X<0)=\Phi(-1{,}19)=0{,}1170.

Ce modèle donnerait à 11,7 %11{,}7\ \% des maisons une concentration NÉGATIVE, ce qui est impossible : la loi normale ne décrit pas une maison.

b)

Tirage avec remise : les XiX_{i} sont indépendantes et ont toutes la loi de la population, d'espérance μ\mu et de variance σ2\sigma^{2}.

E(X‾)=1n(E(X1)+⋯+E(Xn))=1n×nμ=μE(\overline{X})=\dfrac{1}{n}\left(E(X_{1})+\dots+E(X_{n})\right)=\dfrac{1}{n}\times n\mu=\mu (linéarité de l'espérance).

V(X‾)=1n2(V(X1)+⋯+V(Xn))=nσ2n2=σ2nV(\overline{X})=\dfrac{1}{n^{2}}\left(V(X_{1})+\dots+V(X_{n})\right)=\dfrac{n\sigma^{2}}{n^{2}}=\dfrac{\sigma^{2}}{n} : V(aY)=a2V(Y)V(aY)=a^{2}V(Y), et la variance d'une somme de variables INDÉPENDANTES est la somme des variances.

c)

Les 5050 mesures sont indépendantes et de même loi, et n=50≥30n=50\ge 30 : par le théorème central limite, X‾\overline{X} suit approximativement une loi normale, même si XX ne la suit pas.

Son espérance est 9595 et son écart-type, l'erreur type, vaut 8050≈11,31\dfrac{80}{\sqrt{50}}\approx 11{,}31 Bq/m³.

z=120−9511,31≈2,21z=\dfrac{120-95}{11{,}31}\approx 2{,}21 et P(X‾>120)=1−Φ(2,21)=1−0,9864=0,0136P(\overline{X}>120)=1-\Phi(2{,}21)=1-0{,}9864=0{,}0136

Piège : diviser par σ=80\sigma=80 au lieu de l'erreur type donnerait z≈0,31z\approx 0{,}31 ; c'est une question sur la MOYENNE de 5050 maisons, pas sur une maison.

d)

Paramètre : μ\mu, la concentration moyenne des maisons de la municipalité.

H0:μ=95H_{0}:\mu=95 contre H1:μ>95H_{1}:\mu>95 (unilatéral à droite : l'étude cherche un dépassement).

n=40≥30n=40\ge 30 : par le théorème central limite, on utilise la loi normale, avec ss à la place de σ\sigma.

z=x‾−μ0s/n=118−9586/40=2313,60≈1,69\displaystyle z=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}}=\dfrac{118-95}{86/\sqrt{40}}=\dfrac{23}{13{,}60}\approx 1{,}69
e)

Au seuil de 5 %5\ \% : z=1,69>1,645z=1{,}69>1{,}645, on rejette H0H_{0}.

Au seuil de 1 %1\ \% : z=1,69<2,326z=1{,}69<2{,}326, on ne rejette pas H0H_{0}.

Valeur pp : P(Z≥1,69)=1−Φ(1,69)=1−0,9545=0,0455P(Z\ge 1{,}69)=1-\Phi(1{,}69)=1-0{,}9545=0{,}0455, comprise entre 0,010{,}01 et 0,050{,}05, ce qui explique les deux décisions.

Au seuil de 5 %5\ \%, les données permettent de conclure que la concentration moyenne de la municipalité dépasse la moyenne régionale ; au seuil de 1 %1\ \%, elles ne le permettent pas.

f)

Le test porte sur la MOYENNE municipale et la compare à 9595 Bq/m³, pas à 200200 Bq/m³.

Il ne dit rien de la proportion de maisons au-delà de 200200 Bq/m³ : avec une distribution aussi asymétrique, une moyenne de 118118 s'accorde avec une majorité de maisons sous la ligne directrice et quelques-unes très au-dessus.

Seule une mesure dans chaque maison dit si elle doit agir.

Piège : déclarer le test invalide parce que les concentrations ne sont pas normales ; avec n=40n=40, le théorème central limite rend la statistique approximativement normale.

0−2−11,6452,326z = 1,691 %rejet : 5 %Z sous H₀
La statistique sous H0H_{0} : z≈1,69z\approx 1{,}69 tombe dans la région de rejet à 5 %5\ \%, mais pas au-delà de 2,3262{,}326, la valeur critique à 1 %1\ \%.

Exercice 8 : Le calcium d'une source : intervalle et test de Student

9 réponses vérifiables

9 points Intervalles de confiance : moyenne, Student et proportionTests d'hypothèse : moyenne, Student et proportion

Une source karstique alimente le réseau d'eau d'un village. Un rapport hydrogéologique ancien donnait une concentration moyenne de calcium de 48,048{,}0 mg/L. Cette année, un laboratoire de cégep fait 99 prélèvements à des dates tirées au hasard au cours de la saison et dose le calcium par titrage. Résultats, en mg/L (données fictives) : 46,246{,}2 ; 48,148{,}1 ; 48,448{,}4 ; 45,845{,}8 ; 47,547{,}5 ; 45,745{,}7 ; 48,648{,}6 ; 47,347{,}3 ; 47,247{,}2.

On suppose que les concentrations de la saison suivent une loi normale. Valeurs de la loi de Student à 88 degrés de liberté : t0,025 ;8=2,306t_{0{,}025\,;8}=2{,}306 et t0,05 ;8=1,860t_{0{,}05\,;8}=1{,}860, qui laissent une aire de 0,0250{,}025 et 0,050{,}05 à leur droite.

45,045,546,046,547,047,548,048,549,0rapport ancien : 48,0calcium (mg/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 plutôt que la loi normale, et donnez son nombre de degrés de liberté.
  • c) Construisez l'intervalle de confiance à 95 %95\ \% de la concentration moyenne μ\mu de la saison.
  • d) Testez H0:μ=48,0H_{0}:\mu=48{,}0 contre H1:μ≠48,0H_{1}:\mu\neq 48{,}0 au seuil de 5 %5\ \% : statistique de test, valeur critique et décision rédigée.
  • e) Un étudiant compare la statistique à 1,961{,}96. Sa décision serait-elle différente ? Expliquez son erreur.
  • f) Après avoir vu les données, une technicienne propose de tester plutôt H1:μ<48,0H_{1}:\mu<48{,}0, qui donnerait un rejet. Expliquez pourquoi cette démarche est refusée.
  • g) Vérifiez que la décision du d) concorde avec l'intervalle du c), et dites pourquoi.

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

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

Réponses

  • a) x‾=47,20\overline{x}=47{,}20 mg/L, s≈1,091s\approx 1{,}091 mg/L
  • b) n<30n<30 et σ\sigma inconnu, population supposée normale : Student à 88 ddl
  • c) IC95 %=[46,36 ;48,04]IC_{95\,\%}=[46{,}36\,;48{,}04] mg/L
  • d) t≈−2,20t\approx -2{,}20, ∣t∣<2,306|t|<2{,}306 : on ne rejette pas H0H_{0}
  • e) Oui, 1,961{,}96 ferait rejeter : c'est la valeur de la loi normale, pas de Student à 88 ddl
  • f) Le sens de H1H_{1} se fixe avant les données ; choisi après, le risque réel double
  • g) 48,0∈[46,36 ;48,04]48{,}0\in[46{,}36\,;48{,}04] et non-rejet : test bilatéral et intervalle concordent toujours
a)

∑xi=424,8\sum x_{i}=424{,}8, donc x‾=424,89=47,20\overline{x}=\dfrac{424{,}8}{9}=47{,}20 mg/L.

Écarts à la moyenne : −1,0-1{,}0 ; 0,90{,}9 ; 1,21{,}2 ; −1,4-1{,}4 ; 0,30{,}3 ; −1,5-1{,}5 ; 1,41{,}4 ; 0,10{,}1 ; 00, de carrés de somme 9,529{,}52.

s=9,528=1,19≈1,091s=\sqrt{\dfrac{9{,}52}{8}}=\sqrt{1{,}19}\approx 1{,}091 mg/L

b)

L'échantillon est petit (n=9<30n=9<30) et l'écart-type de la population est inconnu : on l'estime par ss, ce qui ajoute de l'incertitude.

Si la population est normale, T=X‾−μS/nT=\dfrac{\overline{X}-\mu}{S/\sqrt{n}} suit la loi de Student à n−1=8n-1=8 degrés de liberté, plus étalée que la loi normale.

Le diagramme des 99 valeurs ne montre ni asymétrie forte ni valeur isolée : l'hypothèse de normalité est plausible, sans être prouvée par si peu de données.

c)

Erreur type : s9=1,0913≈0,3636\dfrac{s}{\sqrt{9}}=\dfrac{1{,}091}{3}\approx 0{,}3636 mg/L.

Marge : 2,306×0,3636≈0,8392{,}306\times 0{,}3636\approx 0{,}839 mg/L.

IC95 %=[47,20−0,84 ;47,20+0,84]=[46,36 ;48,04]IC_{95\,\%}=[47{,}20-0{,}84\,;47{,}20+0{,}84]=[46{,}36\,;48{,}04] mg/L

d)
t=x‾−μ0s/n=47,20−48,00,3636≈−2,20\displaystyle t=\dfrac{\overline{x}-\mu_{0}}{s/\sqrt{n}}=\dfrac{47{,}20-48{,}0}{0{,}3636}\approx -2{,}20

Test bilatéral à 5 %5\ \% : on rejette si ∣t∣>t0,025 ;8=2,306|t|>t_{0{,}025\,;8}=2{,}306.

∣t∣=2,20<2,306|t|=2{,}20<2{,}306 : on ne rejette pas H0H_{0}.

Au seuil de 5 %5\ \%, les données ne permettent pas de conclure que la concentration moyenne de calcium a changé depuis le rapport. Elles ne prouvent pas non plus qu'elle est restée à 48,048{,}0 mg/L : avec 99 mesures, un écart réel de près de 11 mg/L passe facilement inaperçu.

e)

Oui : avec 1,961{,}96, ∣t∣=2,20>1,96|t|=2{,}20>1{,}96 conduirait à rejeter H0H_{0}.

Piège : 1,961{,}96 est la valeur critique de la loi normale ; à 88 degrés de liberté, la loi de Student a des queues plus épaisses, et la valeur critique monte à 2,3062{,}306. Prendre 1,961{,}96 fait commettre l'erreur de première espèce plus souvent que 5 %5\ \% des fois.

f)

Le sens de H1H_{1} se fixe AVANT de voir les données, d'après la question de l'étude, ici « la concentration a-t-elle changé ? », qui est bilatérale.

Choisir le côté après coup, celui où tombe x‾\overline{x}, revient à tester dans les deux sens en ne gardant que le meilleur : le vrai risque de première espèce de cette démarche vaut 2×5 %=10 %2\times 5\ \%=10\ \%, pas 5 %5\ \%.

Le rejet obtenu ainsi (t=−2,20<−1,860t=-2{,}20<-1{,}860) n'a donc pas la valeur annoncée.

g)

μ0=48,0\mu_{0}=48{,}0 appartient à [46,36 ;48,04][46{,}36\,;48{,}04], de justesse, et le test ne rejette pas H0H_{0}, de justesse aussi.

C'est attendu : ∣t∣≤2,306|t|\le 2{,}306 équivaut à ∣x‾−μ0∣≤2,306 sn|\overline{x}-\mu_{0}|\le 2{,}306\,\dfrac{s}{\sqrt{n}}, c'est-à-dire à μ0∈[x‾−2,306 sn ; x‾+2,306 sn]\mu_{0}\in\left[\overline{x}-2{,}306\,\dfrac{s}{\sqrt{n}}\,;\,\overline{x}+2{,}306\,\dfrac{s}{\sqrt{n}}\right]. Un test bilatéral au seuil α\alpha et l'intervalle de niveau 1−α1-\alpha décident toujours de la même façon.

0−2,3062,306−1,961,96t = −2,20loi normaleT sous H₀, 8 ddl
Sous H0H_{0}, TT suit la loi de Student à 88 ddl (violet), plus étalée que la loi normale (pointillé gris). t≈−2,20t\approx -2{,}20 dépasse −1,96-1{,}96 mais pas −2,306-2{,}306 : pas de rejet.

Exercice 9 : Les tiques et la maladie de Lyme : intervalle et test sur une proportion

9 réponses vérifiables

8 points Intervalles de confiance : moyenne, Student et proportionTests d'hypothèse : moyenne, Student et proportion

Dans une région où la maladie de Lyme progresse, une équipe de santé publique récolte des tiques à pattes noires en traînant un drap blanc dans la végétation de sites tirés au hasard. Elle analyse 300300 tiques adultes : 6363 portent la bactérie responsable de la maladie, Borrelia burgdorferi. L'an dernier, la même campagne avait estimé cette proportion à 15 %15\ \%. Données fictives.

Valeurs utiles : z0,025=1,96z_{0{,}025}=1{,}96 et z0,01=2,326z_{0{,}01}=2{,}326 ; table de Φ\Phi avec zz arrondi au centième.

  • a) Calculez la proportion p^\hat{p} de tiques porteuses de l'échantillon, puis l'intervalle de confiance à 95 %95\ \% de la proportion pp dans la région.
  • b) Choisissez l'interprétation correcte du niveau de confiance de 95 %95\ \%.
  • c) Testez au seuil de 1 %1\ \% si la proportion de tiques porteuses a augmenté depuis l'an dernier : hypothèses, conditions, statistique de test et décision.
  • d) Calculez la valeur p du test et rédigez la conclusion.
  • e) L'équipe veut une marge d'erreur d'au plus 33 points de pourcentage à 95 %95\ \% l'an prochain. Calculez la taille d'échantillon nécessaire en prenant p^\hat{p} comme estimation préalable, puis sans aucune estimation préalable.

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

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

Réponses

  • a) p^=0,21\hat{p}=0{,}21 ; IC95 %=[0,164 ;0,256]IC_{95\,\%}=[0{,}164\,;0{,}256]
  • b) 95 %95\ \% des échantillons donnent un intervalle qui contient pp
  • c) H1:p>0,15H_{1}:p>0{,}15 ; z≈2,91>2,326z\approx 2{,}91>2{,}326 : rejet de H0H_{0}
  • d) Valeur p ≈0,0018\approx 0{,}0018 : la proportion a augmenté (seuil de 1 %1\ \%)
  • e) n=709n=709 ; sans estimation, n=1068n=1068
a)
p^=63300=0,21\displaystyle \hat{p}=\dfrac{63}{300}=0{,}21

Conditions : n=300≥30n=300\ge 30, np^=63≥5n\hat{p}=63\ge 5 et n(1−p^)=237≥5n(1-\hat{p})=237\ge 5.

Erreur type estimée : 0,21×0,79300≈0,0235\sqrt{\dfrac{0{,}21\times 0{,}79}{300}}\approx 0{,}0235 ; marge : 1,96×0,0235≈0,0461{,}96\times 0{,}0235\approx 0{,}046.

IC95 %=[0,21−0,046 ;0,21+0,046]=[0,164 ;0,256]\displaystyle IC_{95\,\%}=[0{,}21-0{,}046\,;0{,}21+0{,}046]=[0{,}164\,;0{,}256]
b)

pp est un nombre fixe, inconnu ; c'est l'intervalle qui varie d'un échantillon à l'autre.

95 %95\ \% est la proportion des échantillons possibles pour lesquels la MÉTHODE produit un intervalle qui contient pp. L'intervalle [0,164 ;0,256][0{,}164\,;0{,}256], lui, contient pp ou ne le contient pas.

c)

H0:p=0,15H_{0}:p=0{,}15 contre H1:p>0,15H_{1}:p>0{,}15 (unilatéral à droite : l'étude cherche une hausse).

Conditions sous H0H_{0} : np0=45≥5np_{0}=45\ge 5 et n(1−p0)=255≥5n(1-p_{0})=255\ge 5.

L'erreur type se calcule SOUS H0H_{0}, avec p0p_{0} : 0,15×0,85300≈0,02062\sqrt{\dfrac{0{,}15\times 0{,}85}{300}}\approx 0{,}02062.

z=0,21−0,150,02062≈2,91>2,326z=\dfrac{0{,}21-0{,}15}{0{,}02062}\approx 2{,}91>2{,}326 : on rejette H0H_{0} au seuil de 1 %1\ \%.

Piège : mettre p^\hat{p} dans l'erreur type donne z≈2,55z\approx 2{,}55 ; la décision est la même ici, mais la statistique n'est plus celle du test.

Nuance : des tiques récoltées sur un même site ne sont pas tout à fait indépendantes ; tirer de nombreux sites au hasard limite cet effet de grappe.

d)

Valeur pp : P(Z≥2,91)=1−Φ(2,91)=1−0,9982=0,0018P(Z\ge 2{,}91)=1-\Phi(2{,}91)=1-0{,}9982=0{,}0018, bien sous 0,010{,}01.

Au seuil de 1 %1\ \%, les données permettent de conclure que la proportion de tiques porteuses de Borrelia burgdorferi dans la région a augmenté depuis l'an dernier.

Le test ne dit pas de combien : l'intervalle du a) situe la proportion actuelle entre environ 16 %16\ \% et 26 %26\ \%.

e)

On veut 1,96p(1−p)n≤0,031{,}96\sqrt{\dfrac{p(1-p)}{n}}\le 0{,}03, soit n≥1,962 p(1−p)0,032n\ge\dfrac{1{,}96^{2}\,p(1-p)}{0{,}03^{2}}.

Avec p=0,21p=0{,}21 : n≥3,8416×0,16590,0009≈708,1n\ge\dfrac{3{,}8416\times 0{,}1659}{0{,}0009}\approx 708{,}1, donc n=709n=709 tiques (arrondi VERS LE HAUT).

Sans estimation préalable, on prend p=0,5p=0{,}5, qui rend p(1−p)p(1-p) maximal : n≥3,8416×0,250,0009≈1067,1n\ge\dfrac{3{,}8416\times 0{,}25}{0{,}0009}\approx 1067{,}1, donc n=1068n=1068.

Piège : arrondir 708,1708{,}1 à 708708 donnerait une marge un peu supérieure à 33 points.

0,100,120,140,160,180,200,220,240,260,280,30IC à 95 % : [0,164 ; 0,256]p₀ = 0,15proportion p
L'intervalle de confiance à 95 %95\ \% ne contient pas p0=0,15p_{0}=0{,}15, la proportion de l'an dernier : cohérent avec le rejet du test.

Partie D : lien entre deux variables (/8)

Exercice 10 : Oxygène dissous et température : la droite des moindres carrés et ses limites

10 réponses vérifiables

8 points Régression linéaire et test d'indépendance du khi-carré

À sept dates d'une même année, une équipe mesure vers midi, au même point d'une rivière bien brassée, la température xx de l'eau (en °C) et la concentration yy d'oxygène dissous (en mg/L). Données fictives mais plausibles :

xx : 44 ; 88 ; 1212 ; 1515 ; 1919 ; 2222 ; 2525 • yy : 12,612{,}6 ; 11,511{,}5 ; 10,610{,}6 ; 9,99{,}9 ; 9,29{,}2 ; 8,68{,}6 ; 8,28{,}2.

On donne ∑xi=105\sum x_{i}=105, ∑yi=70,6\sum y_{i}=70{,}6, ∑xi2=1919\sum x_{i}^{2}=1919, ∑yi2=727,22\sum y_{i}^{2}=727{,}22 et ∑xiyi=987,1\sum x_{i}y_{i}=987{,}1. Notations : Sxx=∑xi2−(∑xi)2nS_{xx}=\sum x_{i}^{2}-\frac{(\sum x_{i})^{2}}{n}, Sxy=∑xiyi−∑xi∑yinS_{xy}=\sum x_{i}y_{i}-\frac{\sum x_{i}\sum y_{i}}{n}, SyyS_{yy} de même.

78910111213051015202530T (°C)O₂ dissous (mg/L)
  • a) Calculez SxxS_{xx} et SxyS_{xy}, puis la pente aa et l'ordonnée à l'origine bb de la droite des moindres carrés y^=ax+b\hat{y}=ax+b.
  • b) Démontrez que la droite des moindres carrés passe par le point moyen (x‾ ;y‾)(\overline{x}\,;\overline{y}) et que la somme des résidus ei=yi−y^ie_{i}=y_{i}-\hat{y}_{i} est nulle, sachant que b=y‾−ax‾b=\overline{y}-a\overline{x}.
  • c) Calculez le coefficient de corrélation rr et le coefficient de détermination r2r^{2}, puis interprétez r2r^{2}.
  • d) Estimez la concentration d'oxygène dissous quand l'eau est à 1717 °C.
  • e) Les résidus valent, dans l'ordre des températures : 0,220{,}22 ; −0,05-0{,}05 ; −0,11-0{,}11 ; −0,19-0{,}19 ; −0,05-0{,}05 ; −0,02-0{,}02 ; 0,200{,}20. Décrivez leur allure et ce qu'elle indique.
  • f) Un étudiant prédit l'oxygène dissous d'une canicule où l'eau atteint 3232 °C. Calculez sa prédiction, puis critiquez-la.
  • g) Choisissez l'affirmation correcte sur la causalité.

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

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

Réponses

  • a) Sxx=344S_{xx}=344, Sxy=−71,9S_{xy}=-71{,}9 ; a≈−0,2090a\approx -0{,}2090, b≈13,22b\approx 13{,}22
  • b) ax‾+b=y‾a\overline{x}+b=\overline{y} et ∑ei=n(y‾−ax‾−b)=0\sum e_{i}=n(\overline{y}-a\overline{x}-b)=0
  • c) r≈−0,9954r\approx -0{,}9954, r2≈0,991r^{2}\approx 0{,}991 : 99 %99\ \% de la variation expliquée
  • d) y^(17)≈9,67\hat{y}(17)\approx 9{,}67 mg/L (interpolation)
  • e) Arc : positifs aux extrémités, négatifs au milieu ; légère courbure
  • f) y^(32)≈6,53\hat{y}(32)\approx 6{,}53 mg/L : extrapolation, probablement sous-estimée
  • g) Association seule, sans preuve de causalité ; le mécanisme physique l'appuie
a)
Sxx=1919−10527=1919−1575=344\displaystyle S_{xx}=1919-\dfrac{105^{2}}{7}=1919-1575=344
Sxy=987,1−105×70,67=987,1−1059=−71,9\displaystyle S_{xy}=987{,}1-\dfrac{105\times 70{,}6}{7}=987{,}1-1059=-71{,}9

a=SxySxx=−71,9344≈−0,2090a=\dfrac{S_{xy}}{S_{xx}}=\dfrac{-71{,}9}{344}\approx -0{,}2090 mg/L par °C

x‾=15\overline{x}=15, y‾=70,67≈10,086\overline{y}=\dfrac{70{,}6}{7}\approx 10{,}086, donc b=y‾−ax‾≈10,086+0,2090×15≈13,22b=\overline{y}-a\overline{x}\approx 10{,}086+0{,}2090\times 15\approx 13{,}22 mg/L

y^=13,22−0,209 x\hat{y}=13{,}22-0{,}209\,x : chaque degré de plus s'accompagne en moyenne de 0,210{,}21 mg/L d'oxygène en moins.

b)

En x=x‾x=\overline{x} : y^=ax‾+b=ax‾+y‾−ax‾=y‾\hat{y}=a\overline{x}+b=a\overline{x}+\overline{y}-a\overline{x}=\overline{y}, donc la droite passe par (x‾ ;y‾)(\overline{x}\,;\overline{y}).

∑ei=∑(yi−axi−b)=ny‾−anx‾−nb=n(y‾−ax‾−b)=0\begin{aligned}\sum e_{i}&=\sum\left(y_{i}-ax_{i}-b\right)\\&=n\overline{y}-an\overline{x}-nb\\&=n\left(\overline{y}-a\overline{x}-b\right)=0\end{aligned}

Les points au-dessus et au-dessous de la droite se compensent exactement.

c)
Syy=727,22−70,627≈15,169\displaystyle S_{yy}=727{,}22-\dfrac{70{,}6^{2}}{7}\approx 15{,}169

r=SxySxxSyy=−71,9344×15,169≈−0,9954r=\dfrac{S_{xy}}{\sqrt{S_{xx}S_{yy}}}=\dfrac{-71{,}9}{\sqrt{344\times 15{,}169}}\approx -0{,}9954 et r2≈0,991r^{2}\approx 0{,}991

Liaison linéaire négative très forte : environ 99 %99\ \% de la variation de l'oxygène dissous entre ces sept dates est expliquée par la droite en fonction de la température.

d)

1717 °C est dans la plage observée (44 à 2525 °C) : c'est une interpolation.

y^=13,22−0,2090×17≈9,67\hat{y}=13{,}22-0{,}2090\times 17\approx 9{,}67 mg/L

e)

Les résidus sont positifs aux deux extrémités et négatifs au milieu : ils dessinent un arc, pas un nuage sans forme.

La relation est légèrement COURBE (convexe) : la droite est un bon résumé dans la plage, malgré r2=0,99r^{2}=0{,}99, mais elle n'est pas la forme exacte de la relation.

Un rr proche de −1-1 ne garantit pas qu'une droite soit le bon modèle : il faut regarder les résidus.

f)

y^=13,22−0,2090×32≈6,53\hat{y}=13{,}22-0{,}2090\times 32\approx 6{,}53 mg/L

C'est une EXTRAPOLATION : 3232 °C sort de la plage observée, où rien ne garantit que la relation reste la même.

Les résidus positifs aux extrémités montrent que les données remontent au-dessus de la droite : prolongée, elle sous-estime probablement l'oxygène dissous à 3232 °C, dans ce modèle d'eau bien brassée.

Piège : une prédiction hors de la plage des xx observés n'a pas la fiabilité que suggère r2r^{2}.

g)

La régression mesure une association ; seule, elle ne prouve pas que la température cause la baisse.

Ici, un mécanisme physique connu, la solubilité d'un gaz qui diminue quand l'eau se réchauffe, appuie une causalité ; d'autres facteurs varient pourtant avec la saison (photosynthèse, décomposition, débit) et contribuent aussi.

78910111213051015202530T (°C)O₂ dissous (mg/L)(x̄ ; ȳ)?ŷ = 13,22 − 0,209 x
La droite des moindres carrés passe par le point moyen (vert). Les résidus, trop petits pour se lire à cette échelle, sont donnés en e) ; le prolongement au-delà de 2525 °C (rouge) n'est pas garanti.

Partie E : problème de synthèse (/15)

Exercice 11 : Problème de synthèse : le sel de déglaçage et les grenouilles des bois

15 réponses vérifiables

15 points Statistique descriptive : tendance centrale, dispersion et positionIntervalles de confiance : moyenne, Student et proportionDénombrement : permutations, arrangements, combinaisons, Pascal et binômeRégression linéaire et test d'indépendance du khi-carré

Le sel épandu sur les routes en hiver ruisselle vers les milieux humides voisins. Pour en mesurer l'effet, une équipe d'écologie tire au hasard 4040 étangs parmi ceux d'une région situés à moins de 200200 m d'une route. Données fictives.

Document 1. Concentration de chlorures mesurée au printemps dans les 4040 étangs, en mg/L : minimum 1818, Q1=46Q_{1}=46, médiane 7878, Q3=152Q_{3}=152, maximum 610610 ; les deux plus fortes valeurs sont 340340 et 610610 ; moyenne x‾=112\overline{x}=112 et écart-type s=96s=96. Quatorze étangs dépassent 120120 mg/L, le repère d'exposition prolongée retenu par l'étude. La boîte à moustaches ci-dessous résume ces données.

Document 2. Au même printemps, l'équipe note dans chaque étang la présence ou l'absence de masses d'œufs de grenouille des bois, selon la classe de chlorures : faible (moins de 5050 mg/L), moyenne (de 5050 à 120120 mg/L), élevée (plus de 120120 mg/L).

Valeurs utiles : z0,025=1,96z_{0{,}025}=1{,}96 ; loi du khi-carré à 22 degrés de liberté, χ0,05 ;22=5,991\chi^{2}_{0{,}05\,;2}=5{,}991 et χ0,01 ;22=9,210\chi^{2}_{0{,}01\,;2}=9{,}210, valeurs qui laissent une aire de 0,050{,}05 et 0,010{,}01 à leur droite.

Document 2FaibleMoyenneÉlevéeTotal
Pontes présentes108321
Pontes absentes261119
Total12141440
0100200300400500600120 mg/Lchlorures (mg/L)
  • a) À l'aide du document 1, calculez l'écart interquartile et les bornes au-delà desquelles une concentration est dite aberrante, puis dites combien d'étangs sont aberrants.
  • b) Décrivez la forme de la distribution, puis dites quelle mesure, de la moyenne ou de la médiane, décrit le mieux un étang typique.
  • c) Construisez l'intervalle de confiance à 95 %95\ \% de la concentration moyenne μ\mu de chlorures des étangs de la région, en justifiant la méthode malgré l'asymétrie.
  • d) Une technicienne retourne dans 33 des 4040 étangs, tirés au hasard sans remise. Calculez la probabilité qu'au moins un d'entre eux dépasse 120120 mg/L.
  • e) Pour le document 2, calculez les effectifs théoriques sous l'hypothèse d'indépendance entre la présence de pontes et la classe de chlorures, puis vérifiez la condition d'application du test du khi-carré.
  • f) Effectuez le test d'indépendance du khi-carré : hypothèses, degrés de liberté, statistique et décision aux seuils de 5 % et de 1 %.
  • g) Choisissez la conclusion que l'étude permet de rédiger, et justifiez à l'aide des proportions de pontes par classe.

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

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

Réponses

  • a) EI=106EI=106 mg/L ; bornes −113-113 et 311311 ; 22 étangs aberrants (340340 et 610610)
  • b) Fortement asymétrique à droite (x‾=112>Me=78\overline{x}=112>Me=78) : la médiane
  • c) IC95 %≈[82,25 ;141,75]IC_{95\,\%}\approx[82{,}25\,;141{,}75] mg/L
  • d) 1−26009880≈0,73681-\frac{2600}{9880}\approx 0{,}7368
  • e) 6,36{,}3 ; 7,357{,}35 ; 7,357{,}35 / 5,75{,}7 ; 6,656{,}65 ; 6,656{,}65 : tous ≥5\ge 5
  • f) 22 ddl ; χ2≈10,12>9,210\chi^{2}\approx 10{,}12>9{,}210 : rejet de l'indépendance à 5 %5\ \% et à 1 %1\ \%
  • g) Pontes : 83 %83\ \%, 57 %57\ \%, 21 %21\ \% ; lien établi, causalité non prouvée (étude d'observation)
a)

EI=Q3−Q1=152−46=106EI=Q_{3}-Q_{1}=152-46=106 mg/L

Bornes : Q1−1,5 EI=46−159=−113Q_{1}-1{,}5\,EI=46-159=-113 et Q3+1,5 EI=152+159=311Q_{3}+1{,}5\,EI=152+159=311 mg/L.

Aucune concentration n'est sous −113-113 ; au-dessus de 311311, il y a 340340 et 610610 : 22 étangs aberrants.

Aberrant ne veut pas dire faux : ces deux étangs sont peut-être au pied d'un dépôt de sel ou d'une forte pente ; on les garde, et on les vérifie.

b)

La médiane (7878) est à 3232 mg/L de Q1Q_{1} et à 7474 mg/L de Q3Q_{3}, donc bien plus près de Q1Q_{1} ; la moustache droite s'étire jusqu'à 610610, et x‾=112>Me=78\overline{x}=112>Me=78 : la distribution est fortement asymétrique à droite.

La médiane décrit mieux un étang typique : la moyenne est tirée vers le haut par quelques étangs très salés.

c)

n=40≥30n=40\ge 30 et les étangs sont tirés au hasard : par le théorème central limite, X‾\overline{X} est approximativement normale malgré l'asymétrie des concentrations ; avec nn grand, on remplace σ\sigma par ss.

Erreur type : 9640≈15,18\dfrac{96}{\sqrt{40}}\approx 15{,}18 mg/L ; marge : 1,96×15,18≈29,751{,}96\times 15{,}18\approx 29{,}75 mg/L.

IC95 %=[112−29,75 ;112+29,75]≈[82,25 ;141,75]IC_{95\,\%}=[112-29{,}75\,;112+29{,}75]\approx[82{,}25\,;141{,}75] mg/L

Nuance : avec une distribution aussi étirée, n=40n=40 reste juste ; les deux valeurs aberrantes pèsent lourd sur x‾\overline{x} et sur ss.

d)

Contraire : aucun des 33 étangs ne dépasse 120120 mg/L, soit 33 étangs choisis parmi les 2626 autres.

P=1−(263)(403)=1−26009880≈1−0,2632=0,7368\displaystyle P=1-\dfrac{\binom{26}{3}}{\binom{40}{3}}=1-\dfrac{2600}{9880}\approx 1-0{,}2632=0{,}7368

Piège : une binomiale 1−(2640)3≈0,72541-\left(\frac{26}{40}\right)^{3}\approx 0{,}7254 supposerait un tirage AVEC remise ; sans remise sur 4040 étangs, on dénombre.

e)

Effectif théorique : total de la ligne×total de la colonne40\dfrac{\text{total de la ligne}\times\text{total de la colonne}}{40}.

Pontes présentes : 21×1240=6,3\dfrac{21\times 12}{40}=6{,}3 ; 21×1440=7,35\dfrac{21\times 14}{40}=7{,}35 ; 7,357{,}35.

Pontes absentes : 19×1240=5,7\dfrac{19\times 12}{40}=5{,}7 ; 6,656{,}65 ; 6,656{,}65.

Le plus petit vaut 5,7≥55{,}7\ge 5 : la condition des effectifs théoriques est satisfaite, de justesse.

f)

H0H_{0} : la présence de pontes et la classe de chlorures sont indépendantes ; H1H_{1} : elles sont liées.

Degrés de liberté : (2−1)(3−1)=2(2-1)(3-1)=2.

χ2=(10−6,3)26,3+(8−7,35)27,35+(3−7,35)27,35+(2−5,7)25,7+(6−6,65)26,65+(11−6,65)26,65\displaystyle \chi^{2}=\dfrac{(10-6{,}3)^{2}}{6{,}3}+\dfrac{(8-7{,}35)^{2}}{7{,}35}+\dfrac{(3-7{,}35)^{2}}{7{,}35}+\dfrac{(2-5{,}7)^{2}}{5{,}7}+\dfrac{(6-6{,}65)^{2}}{6{,}65}+\dfrac{(11-6{,}65)^{2}}{6{,}65}
χ2≈2,173+0,057+2,574+2,402+0,064+2,845≈10,12\displaystyle \chi^{2}\approx 2{,}173+0{,}057+2{,}574+2{,}402+0{,}064+2{,}845\approx 10{,}12

10,12>5,99110{,}12>5{,}991 et 10,12>9,21010{,}12>9{,}210 : on rejette H0H_{0} aux deux seuils.

Les plus grandes contributions viennent des classes extrêmes : trop de pontes en eau peu salée, trop peu en eau très salée, par rapport à l'indépendance.

g)

Proportions d'étangs avec pontes : 1012≈83 %\frac{10}{12}\approx 83\ \% en classe faible, 814≈57 %\frac{8}{14}\approx 57\ \% en classe moyenne, 314≈21 %\frac{3}{14}\approx 21\ \% en classe élevée.

Conclusion : au seuil de 1 %1\ \%, les données montrent un lien entre la salinité des étangs et la présence de pontes de grenouille des bois ; les pontes sont plus rares dans les étangs les plus salés.

L'étude est une étude d'OBSERVATION : elle ne prouve pas que le sel CAUSE l'absence de pontes. Les étangs très salés sont sans doute aussi plus proches de routes plus passantes (bruit, mortalité routière, sédiments) ; seule une expérience contrôlée, ou la prise en compte de ces facteurs, départagerait.

Elle vaut pour les étangs proches des routes de cette région, la population d'où viennent les 4040 étangs.

02468101214khi-carré, 2 ddl5,99189,211214χ² = 10,121 %rejet : 5 %queue droite agrandie (hauteurs × 8)
Loi du khi-carré à 22 degrés de liberté. En bas, la queue droite agrandie : χ2≈10,12\chi^{2}\approx 10{,}12 dépasse 5,9915{,}991 (5 %5\ \%) et 9,219{,}21 (1 %1\ \%).
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é.

© Ahmed Squalli Houssaini. Série publiée sur www.letuteurscientifique.ca/exercices/sn1-epreuve-synthese. 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. On prépare l'épreuve finale sur des sujets complets, en travaillant d'abord ce qui rapporte le plus de points : nommer le modèle, choisir la bonne loi et rédiger la conclusion dans le contexte.

Site par Studio Squalli