Aller au contenu

L'Établi des Maths Cours

Outils pour apprendre et pratiquer

PASS / L.AS · Statistique descriptive · PA-11

Cours et exercices corrigés de PASS / L.AS : Corrélation et régression linéaire

Séries bivariées et nuage de points, covariance, coefficient de corrélation linéaire de Pearson, droite des moindres carrés, coefficient de détermination, prédiction, corrélation et causalité.

Objectifs du chapitre

Le cours

1. Nuage de points et point moyen

Une série bivariée associe à chaque individu deux mesures (xi;yi)(x_i ; y_i), i=1,…,ni = 1, \dots, n (par exemple l'âge et la pression systolique d'un patient). On la représente par un nuage de points.
Le point moyen est G(xˉ;yˉ)G(\bar{x} ; \bar{y}), où xˉ\bar{x} et yˉ\bar{y} sont les moyennes des deux séries.
Avant tout calcul, on regarde la forme du nuage : allongé selon une droite (liaison linéaire), courbe (liaison non linéaire), ou sans structure.

2. Covariance et coefficient de corrélation

Covariance : cov(x,y)=1n∑i=1n(xi−xˉ)(yi−yˉ)=xy‾−xˉ yˉ.\mathrm{cov}(x, y) = \dfrac{1}{n}\sum_{i=1}^{n}(x_i - \bar{x})(y_i - \bar{y}) = \overline{xy} - \bar{x}\,\bar{y}. Elle est positive quand xx et yy varient dans le même sens, négative sinon ; elle dépend des unités.
Coefficient de corrélation linéaire de Pearson : r=cov(x,y)sx sy,r = \dfrac{\mathrm{cov}(x, y)}{s_x\, s_y}, où sxs_x et sys_y sont les écarts types des deux séries (calculés avec le même diviseur que la covariance).
rr est sans unité, toujours compris entre −1-1 et 11 ; ∣r∣=1|r| = 1 si et seulement si les points sont exactement alignés. rr ne change pas si l'on change d'unité (transformation affine de coefficient positif).

3. Droite des moindres carrés

La droite de régression de yy en xx (ou des moindres carrés) est la droite y=ax+by = ax + b qui rend minimale la somme des carrés des écarts verticaux ∑(yi−axi−b)2\sum (y_i - ax_i - b)^2. On a a=cov(x,y)sx2,b=yˉ−axˉ.a = \dfrac{\mathrm{cov}(x, y)}{s_x^2}, \qquad b = \bar{y} - a\bar{x}.
La droite passe par le point moyen GG. La pente aa a le signe de rr et s'exprime en « unité de yy par unité de xx » : c'est la variation moyenne de yy quand xx augmente d'une unité.
On a aussi a=r sysxa = r\,\dfrac{s_y}{s_x}. Attention : la droite de régression de xx en yy est en général différente.

4. Coefficient de détermination, prédiction et causalité

Coefficient de détermination : r2r^2, entre 00 et 11, est la part de la variance de yy « expliquée » par la relation linéaire avec xx. Avec r=0,8r = 0{,}8, r2=0,64r^2 = 0{,}64 : 64 %64\,\% de la variabilité de yy est expliquée par la droite.
Prédiction : pour une valeur x0x_0, on prédit y^=ax0+b\hat{y} = ax_0 + b. Elle n'est raisonnable que dans le domaine des xx observés : au-delà, c'est une extrapolation, risquée car rien ne garantit que la relation reste linéaire.
Corrélation n'est pas causalité : une forte corrélation peut venir d'un facteur de confusion lié aux deux variables (l'âge, la saison…), du hasard, ou d'une causalité inverse. Seul un essai contrôlé randomisé permet de conclure à un effet causal.
Test de nullité de rr (aperçu) : pour tester H0H_0 « pas de corrélation linéaire dans la population », on calcule t=rn−21−r2t = \dfrac{r\sqrt{n - 2}}{\sqrt{1 - r^2}}, que l'on compare au seuil de la loi de Student à n−2n - 2 degrés de liberté (proche de 22 au risque 5 %5\,\% dès que nn dépasse quelques dizaines). Avec un grand effectif, un rr faible peut être significatif sans être cliniquement important.

5. Aux partiels : pièges fréquents

  • Diviser la covariance par sxsys_x s_y pour la pente : c'est le coefficient rr ; la pente est cov(x,y)sx2\dfrac{\mathrm{cov}(x, y)}{s_x^2}.
  • Oublier de vérifier que la droite passe par G(xˉ;yˉ)G(\bar{x} ; \bar{y}) : c'est un contrôle rapide de bb.
  • Croire que r=0r = 0 signifie « aucune liaison » : il n'y a pas de liaison linéaire, mais une liaison courbe (en U) reste possible.
  • Confondre rr et r2r^2 : r=0,5r = 0{,}5 n'explique que 25 %25\,\% de la variance.
  • Interpréter une corrélation comme un effet causal, ou prédire loin hors du domaine observé.

Méthodes et exemples rédigés

Méthode : Calculer rr et la droite des moindres carrés

Exemple

Dose xx (en mg) : 11 ; 22 ; 33 ; 44 ; 55. Baisse de pression yy (en mmHg) : 22 ; 44 ; 55 ; 44 ; 55.

  1. Étape 1 (Point moyen) : xˉ=3\bar{x} = 3, yˉ=205=4\bar{y} = \dfrac{20}{5} = 4, donc G(3;4)G(3 ; 4).

  2. Étape 2 (Écarts et produits) : écarts xi−xˉx_i - \bar{x} : −2,−1,0,1,2-2, -1, 0, 1, 2 ; écarts yi−yˉy_i - \bar{y} : −2,0,1,0,1-2, 0, 1, 0, 1 ; somme des produits : 4+0+0+0+2=64 + 0 + 0 + 0 + 2 = 6.

  3. Étape 3 (Covariance et variances) : cov=65=1,2\mathrm{cov} = \dfrac{6}{5} = 1{,}2 ; sx2=105=2s_x^2 = \dfrac{10}{5} = 2 ; sy2=65=1,2s_y^2 = \dfrac{6}{5} = 1{,}2.

  4. Étape 4 (Droite) : a=1,22=0,6a = \dfrac{1{,}2}{2} = 0{,}6 et b=4−0,6×3=2,2b = 4 - 0{,}6 \times 3 = 2{,}2 : y=0,6x+2,2y = 0{,}6x + 2{,}2.

  5. Étape 5 (Corrélation) : r2=1,222×1,2=0,6r^2 = \dfrac{1{,}2^2}{2 \times 1{,}2} = 0{,}6, donc r=0,6≈0,77r = \sqrt{0{,}6} \approx 0{,}77 : 60 %60\,\% de la variance de yy est expliquée par la dose.

Méthode : Prédire et critiquer une prédiction

Exemple

Chez des enfants de 22 à 1010 ans, la droite de régression de la taille yy (cm) en fonction de l'âge xx (ans) est y=6x+77y = 6x + 77.

  1. Étape 1 (Domaine) : x0=6x_0 = 6 ans est dans le domaine observé ([2;10][2 ; 10]) : la prédiction est légitime.

  2. Étape 2 (Calcul) : y^=6×6+77=113\hat{y} = 6 \times 6 + 77 = 113 cm (valeur moyenne attendue, pas la taille exacte d'un enfant).

  3. Étape 3 (Critique) : pour x0=30x_0 = 30 ans, la formule donnerait 257257 cm : l'extrapolation est absurde car la croissance s'arrête.

Pièges fréquents

Questions de révision

Essaie de répondre de tête avant d'afficher la réponse.

Formule de la covariance (König-Huygens) ?

cov(x,y)=xy‾−xˉ yˉ\mathrm{cov}(x, y) = \overline{xy} - \bar{x}\,\bar{y}.

Coefficient de corrélation de Pearson ?

r=cov(x,y)sxsyr = \dfrac{\mathrm{cov}(x, y)}{s_x s_y}, compris entre −1-1 et 11.

Pente et ordonnée à l'origine de la droite des moindres carrés ?

a=cov(x,y)sx2a = \dfrac{\mathrm{cov}(x, y)}{s_x^2}, b=yˉ−axˉb = \bar{y} - a\bar{x} : la droite passe par G(xˉ;yˉ)G(\bar{x} ; \bar{y}).

Que représente r2r^2 ?

La part de la variance de yy expliquée par la régression linéaire sur xx.

Pourquoi une corrélation ne prouve-t-elle pas la causalité ?

Un facteur de confusion, le hasard ou une causalité inverse peuvent l'expliquer.

Exercices corrigés

Exercice 1 Coordonnée du point moyen

Palier 1 · découverte

On relève l'âge xx (en années) et la pression systolique yy (en mmHg) de 55 patients :

  • âges : 3030 ; 4040 ; 5050 ; 6060 ; 7070 ;
  • pressions : 120120 ; 125125 ; 130130 ; 140140 ; 135135.

Le point moyen du nuage est G(50;yˉ)G(50 ; \bar{y}). Calculez yˉ\bar{y}, en mmHg.

Un indice

L'ordonnée du point moyen est la moyenne de la seconde série, sans lien avec la première.

Voir la correction

Réponse attendue : 130

yˉ=120+125+130+140+1355=6505=\bar{y} = \dfrac{120 + 125 + 130 + 140 + 135}{5} = \dfrac{650}{5} = 130130 mmHg, donc G(50;130)G(50 ; 130).
Toute droite des moindres carrés de ce nuage passera par ce point.

Exercice 2 Équation de la droite des moindres carrés

Palier 2 · application

On étudie la dose xx d'un diurétique (en mg) et la diurèse yy (en mL/h) chez des patients. On obtient xˉ=10\bar{x} = 10, yˉ=50\bar{y} = 50, cov(x,y)=18\mathrm{cov}(x, y) = 18 et sx2=9s_x^2 = 9.
Quelle est l'équation de la droite de régression de yy en xx ?

  1. y=2x+30y = 2x + 30
  2. y=0,5x+45y = 0{,}5x + 45
  3. y=2x+50y = 2x + 50
  4. y=18x−130y = 18x - 130
Un indice

Calculez d'abord la pente, puis utilisez le fait que la droite passe par le point moyen.

Voir la correction

Bonne réponse : A. y=2x+30y = 2x + 30

a=cov(x,y)sx2=189=2a = \dfrac{\mathrm{cov}(x, y)}{s_x^2} = \dfrac{18}{9} = 2 et b=yˉ−axˉ=50−20=30b = \bar{y} - a\bar{x} = 50 - 20 = 30 : y=2x+30y = 2x + 30.
y=2x+50y = 2x + 50 oublie de retrancher axˉa\bar{x} ; y=18x−130y = 18x - 130 prend la covariance pour la pente ; y=0,5x+45y = 0{,}5x + 45 inverse le rapport.

Exercice 3 Pente à partir des données

Palier 3 · approfondissement

On étudie l'effet de la supplémentation en vitamine D. Pour 44 patients, on note la dose quotidienne xx (en centaines d'UI) et l'augmentation de la concentration sérique yy (en ng/mL) après 33 mois :
(2;3)(2 ; 3), (4;7)(4 ; 7), (6;9)(6 ; 9), (8;13)(8 ; 13).
Calculez la pente aa de la droite de régression de yy en xx (réponse décimale exacte).

Un indice

Calculez le point moyen, puis les écarts à la moyenne de chaque série.

Voir la correction

Réponse attendue : 1.6

xˉ=5\bar{x} = 5 et yˉ=8\bar{y} = 8. Écarts : xi−xˉ=−3,−1,1,3x_i - \bar{x} = -3, -1, 1, 3 et yi−yˉ=−5,−1,1,5y_i - \bar{y} = -5, -1, 1, 5.
cov(x,y)=15+1+1+154=8\mathrm{cov}(x, y) = \dfrac{15 + 1 + 1 + 15}{4} = 8 et sx2=9+1+1+94=5s_x^2 = \dfrac{9 + 1 + 1 + 9}{4} = 5.
a=85=a = \dfrac{8}{5} = 1,61{,}6 ng/mL par centaine d'UI ; comme b=8−1,6×5=0b = 8 - 1{,}6 \times 5 = 0, la droite est y=1,6xy = 1{,}6x.

Exercice 4 Score de douleur et rééducation

Palier 4 · défi

On suit le score de douleur yy (sur 1010) d'un patient en rééducation au fil des semaines xx :
(1;10)(1 ; 10), (2;8)(2 ; 8), (3;8)(3 ; 8), (4;5)(4 ; 5), (5;4)(5 ; 4).
Déterminez la droite de régression y=ax+by = ax + b de yy en xx, puis donnez la valeur de bb (réponse décimale exacte).

Un indice

Calculez le point moyen, la covariance et la variance de xx avant d'obtenir la pente, puis bb.

Voir la correction

Réponse attendue : 11.5

xˉ=3\bar{x} = 3, yˉ=355=7\bar{y} = \dfrac{35}{5} = 7. Écarts : xi−xˉ=−2,−1,0,1,2x_i - \bar{x} = -2, -1, 0, 1, 2 ; yi−yˉ=3,1,1,−2,−3y_i - \bar{y} = 3, 1, 1, -2, -3.
cov(x,y)=−6−1+0−2−65=−3\mathrm{cov}(x, y) = \dfrac{-6 - 1 + 0 - 2 - 6}{5} = -3 et sx2=4+1+0+1+45=2s_x^2 = \dfrac{4 + 1 + 0 + 1 + 4}{5} = 2, donc a=−1,5a = -1{,}5.
b=yˉ−axˉ=7+4,5=b = \bar{y} - a\bar{x} = 7 + 4{,}5 = 11,511{,}5 : y=−1,5x+11,5y = -1{,}5x + 11{,}5. Le score baisse en moyenne de 1,51{,}5 point par semaine.