3 Association entre variables qualitatives

Lorsqu’on travaille avec des variables qualitatives, on utilise fréquemment le test du \(\chi^2\). Il en existe plusieurs versions suivant la question de recherche, mais la plus fréquente est le test du \(\chi^2\) d’indépendance, qui évalue si deux variables qualitatives sont associées. Pour l’utiliser, nous devons d’abord construire un tableau de contingence.

3.1 Tableau de contingence

ImportantDéfinition

Un tableau de contingence est un tableau qui croise deux variables qualitatives. Chaque case contient l’effectif des unités statistiques qui possèdent simultanément une modalité de la première variable et une modalité de la seconde.

Pour construire un tableau de contingence, il faut :

  1. Choisir quelle variable sera en ligne et quelle variable sera en colonne ;
  2. Mettre autant de lignes / colonnes qu’il y a de modalités (p. ex., croiser “sexe” (ligne) avec “tabagisme” (colonne), on aura deux lignes (“homme” et “femme”), et trois colonnes : “non fumeur”, “fumeur occasionnel” et “fumeur régulier”) ;
  3. Rajouter une dernière ligne / colonne “total” ;
  4. Remplir les cases avec les informations qu’on a à disposition.
NoteExemple

Dans l’étude HTA menée auprès de \(n=402\) individus, on dispose notamment de deux variables : avoir un IMC \(\ge30\) ou non, et avoir une HTA connue ou non. On sait que 116 personnes ont une HTA connue et que 48 personnes ont un IMC \(\ge30\), dont 24 présentent aussi une HTA. Grâce à ces informations, on peut commencer à remplir le tableau de contingence suivant :

Tableau 1: Tableau de contingence entre obésité et HTA connue dans l’enquête HTA, île Maurice, 19-21 avril 1999, 402 adultes de plus de 18 ans.
Statut pondéral HTA oui HTA non Total
IMC \(\ge30\) 24 ? 48
IMC \(<30\) ? ? ?
Total 116 ? 402

Pour le compléter, on peut utiliser les informations qu’on a à disposition. Par exemple, le total de la première ligne est 48, et on sait que 24 personnes ont une HTA connue. Il reste donc \(48-24=24\) personnes avec IMC \(\ge30\) mais sans HTA connue. On peut remplir les autres cases de la même façon, pour obtenir le tableau de contingence suivant :

Tableau 2: Tableau de contingence entre obésité et HTA connue dans l’enquête HTA, île Maurice, 19-21 avril 1999, 402 adultes de plus de 18 ans.
Statut pondéral HTA oui HTA non Total
IMC \(\ge30\) 24 24 48
IMC \(<30\) 92 262 354
Total 116 286 402

Pour ce qui suit, on distingue “deux types” de tableaux de contingence : le tableau de contingence observé (celui qu’on vient de construire en gros) et le tableau de contingence théorique.

3.2 Test du \(\chi^2\) d’indépendance

Le test du \(\chi^2\) d’indépendance évalue si deux variables qualitatives sont associées. Il s’applique à un tableau de contingence et compare les effectifs observés aux effectifs qui seraient attendus si les deux variables étaient indépendantes.

Le critère du \(\chi^2\) a été proposé par Karl Pearson en 1900 pour comparer des effectifs observés et attendus (Pearson 1900). Ronald Fisher en a ensuite étudié l’interprétation pour les tableaux de contingence en 1922 (Fisher 1922).

Dans l’exemple précédent, les hypothèses sont : \[ H_0:\text{les deux variables sont indépendantes} \] contre \[ H_1:\text{les deux variables ne sont pas indépendantes} \]

Équivalemment, ce test permet de comparer deux proportions. Ainsi, dans notre contexte, on peut écrire : \[ H_0:\text{la proportion de personnes avec HTA est la même chez les personnes obèses et non obèses} \] contre \[ H_1:\text{la proportion de personnes avec HTA est différente chez les personnes obèses et non obèses} \]

Sous \(H_0\), l’effectif attendu / théorique dans la case \((i,j)\) est : \[ E_{ij}=\frac{(\text{total ligne }i)\times(\text{total colonne }j)}{\text{total général}} \]

Pour que l’approximation par la loi du \(\chi^2\) soit raisonnable, on retiendra ici une règle simple et prudente : les effectifs théoriques doivent tous être au moins égaux à 5. Pour un tableau \(2\times2\) qui ne respecte pas cette règle, on préférera le test exact de Fisher.

NoteExemple

On a trouvé avant que :

Tableau 3: Tableau de contingence entre obésité et HTA connue dans l’enquête HTA, île Maurice, 19-21 avril 1999, 402 adultes de plus de 18 ans.
Statut pondéral HTA oui HTA non Total
IMC \(\ge30\) 24 24 48
IMC \(<30\) 92 262 354
Total 116 286 402

L’effectif attendu / théorique de personnes obèses avec HTA sous indépendance est : \[ E_{11}=\frac{48\times116}{402}\approx13,85 \]

En faisant cela pour toutes les cases, on obtient le tableau de contingence théorique suivant :

Tableau 4: Tableau de contingence entre obésité et HTA connue dans l’enquête HTA, île Maurice, 19-21 avril 1999, 402 adultes de plus de 18 ans.
Statut pondéral HTA oui HTA non Total
IMC \(\ge30\) \(\frac{48\times116}{402} \approx 13,85\) \(\frac{48\times286}{402} \approx 34,15\) 48
IMC \(<30\) \(\frac{354\times116}{402} \approx 102,15\) \(\frac{354\times286}{402} \approx 251,85\) 354
Total 116 286 402

Comme tous les effectifs attendus sont supérieurs à 5, le test du \(\chi^2\) d’indépendance est applicable.

Un effectif attendu n’est pas un effectif réellement observé : il peut donc être décimal. Il représente ce que l’on attendrait en moyenne si les deux variables étaient indépendantes, en conservant les totaux de lignes et de colonnes.

Ici, on en observe 24. Cette case contribue donc fortement à l’écart entre le tableau observé et le tableau attendu sous indépendance.

La statistique de test est : \[ \chi^2=\sum_{i,j}\frac{(O_{ij}-E_{ij})^2}{E_{ij}} \underset{H_0}{\overset{\mathrm{approx.}}{\sim}} \chi^2_{(r-1)(c-1)} \] où \(O_{ij}\) est l’effectif observé, \(E_{ij}\) l’effectif attendu, \(r\) le nombre de lignes et \(c\) le nombre de colonnes du tableau.

Au seuil \(\alpha\), la région critique est : \[ \chi^2>\chi^2_{1-\alpha,(r-1)(c-1)} \] Avec \(\alpha=5\%\) dans un tableau 2x2, la valeur critique est \(\chi^2_{0,95,1}\approx3,84\).

NoteExemple

En combinant en un seul tableau les effectifs observés (en bleu) et attendus (en orange), on obtient le tableau suivant :

Tableau 5: Tableau de contingence observé et attendu entre obésité et HTA connue dans l’enquête HTA, île Maurice, 19-21 avril 1999, 402 adultes de plus de 18 ans.
Statut pondéral HTA oui HTA non Total
IMC \(\ge30\) 24 / 13,85 24 / 34,15 48
IMC \(<30\) 92 / 102,15 262 / 251,85 354
Total 116 286 402

Le \(\chi^2\) d’indépendance vaut alors : \[ \begin{aligned} \chi^2 &= \overbrace{\frac{(24-13,85)^2}{13,85}}^{\text{"IMC $\ge30$ et HTA oui"}} + \overbrace{\frac{(24-34,15)^2}{34,15}}^{\text{"IMC $\ge30$ et HTA non"}} + \overbrace{\frac{(92-102,15)^2}{102,15}}^{\text{"IMC $<30$ et HTA oui"}} + \overbrace{\frac{(262-251,85)^2}{251,85}}^{\text{"IMC $<30$ et HTA non"}} \\ &\approx 11,87 \end{aligned} \]

Avec 2 lignes et 2 colonnes, on a \((r-1)(c-1)=(2-1)(2-1)=1\) degré de liberté. La valeur critique d’une loi du \(\chi^2\) à 1 degré de liberté au seuil de 5 % est environ 3,84. Comme \(11,87>3,84\), on rejette l’indépendance au seuil de 5 %. Il existe donc une association entre obésité et HTA connue dans cet échantillon.

3.3 Correction de Yates

Pour un tableau \(2\times2\) (UNIQUEMENT !), on rencontre parfois le test du \(\chi^2\) avec correction de continuité de Yates. Il est adéquat d’utiliser cette correction quand un effectif théorique est compris entre 2,5 et 5. Cependant, ce n’est pas pour autant magique ; et si les effectifs attendus sont trop faibles, le test exact de Fisher reste préférable.

Frank Yates a présenté cette correction en 1934 dans son étude des tableaux à petits effectifs (Yates 1934).

La statistique corrigée remplace l’écart \(O_{ij}-E_{ij}\) par un écart réduit de 0,5 en valeur absolue : \[ \chi^2_{\mathrm{Yates}} =\sum_{i,j}\frac{\max(0\, ; \,|O_{ij}-E_{ij}|-0,5)^2}{E_{ij}} \underset{H_0}{\overset{\mathrm{approx.}}{\sim}} \chi^2_{1} \]

Au seuil \(\alpha\), la région critique reste : \[ \chi^2_{\mathrm{Yates}}>\chi^2_{1-\alpha,(r-1)(c-1)} \]

NoteExemple

On étudie dans une petite étude cas-témoins si le nombre d’accouchements est associé au cancer du sein. On observe le tableau suivant :

Tableau 6: Tableau observé d’une petite étude cas-témoins sur cancer du sein et nombre d’accouchements.
Nombre d’accouchements Cancer du sein Pas de cancer du sein Total
\(\ge2\) 1 5 6
\(<2\) 11 8 19
Total 12 13 25

Le tableau de contingence attendu sous indépendance est :

Tableau 7: Tableau attendu d’une petite étude cas-témoins sur cancer du sein et nombre d’accouchements.
Nombre d’accouchements Cancer du sein Pas de cancer du sein Total
\(\ge2\) \(\frac{6\times12}{25}\approx2,88\) \(\frac{6\times13}{25}\approx3,12\) 6
\(<2\) \(\frac{19\times12}{25}\approx9,12\) \(\frac{19\times13}{25}\approx9,88\) 19
Total 12 13 25

Deux effectifs attendus sont inférieurs à 5. Pour illustrer le calcul de la correction de Yates, on obtient : \[ \chi^2_{\mathrm{Yates}} = \frac{(|1-2,88|-0,5)^2}{2,88} + \frac{(|5-3,12|-0,5)^2}{3,12} + \frac{(|11-9,12|-0,5)^2}{9,12} + \frac{(|8-9,88|-0,5)^2}{9,88} \approx 1,6732 \] Comme avant, la valeur critique d’une loi du \(\chi^2\) à 1 degré de liberté au seuil de 5 % est environ 3,84. Comme \(1,6732<3,84\), on ne rejette pas l’indépendance au seuil de 5 %. Les données ne mettent donc pas en évidence d’association entre le nombre d’accouchements et le cancer du sein ; elles ne prouvent pas pour autant l’absence d’association.

Note

La correction de Yates peut aussi être trop conservatrice (cela signifie qu’on peut ne pas rejeter \(H_0\) alors qu’il y a une association réelle). Si les effectifs sont vraiment faibles dans un tableau \(2\times2\), le test exact de Fisher est souvent plus approprié.

3.4 Test exact de Fisher

Le test exact de Fisher est une alternative au \(\chi^2\) pour les tableaux \(2\times2\) lorsque les effectifs sont petits. Ce test se fait difficilement à la main, mais il est très simple à réaliser avec un logiciel statistique. Il est dit “exact” car sa p-valeur ne repose pas sur l’approximation asymptotique par une loi du \(\chi^2\).

Fisher a illustré un raisonnement exact en 1935 avec son célèbre exemple de dégustation de thé (Fisher 1935). Yates avait déjà décrit un calcul exact pour de petits tableaux en 1934 (Yates 1934).

Les références

Fisher, Ronald A. 1922. « On the Interpretation of χ² from Contingency Tables, and the Calculation of P ». Journal of the Royal Statistical Society 85 (1): 87‑94. https://doi.org/10.1111/j.2397-2335.1922.tb00768.x.
Fisher, Ronald A. 1935. The Design of Experiments. Oliver; Boyd. https://search.worldcat.org/title/The-design-of-experiments/oclc/2417943.
Pearson, Karl. 1900. « On the Criterion that a Given System of Deviations from the Probable in the Case of a Correlated System of Variables Is Such that It Can Be Reasonably Supposed to Have Arisen from Random Sampling ». The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 50 (302): 157‑75. https://doi.org/10.1080/14786440009463897.
Yates, Frank. 1934. « Contingency Tables Involving Small Numbers and the χ² Test ». Supplement to the Journal of the Royal Statistical Society 1 (2): 217‑35. https://doi.org/10.2307/2983604.