{"id":8105,"date":"2026-09-17T04:01:51","date_gmt":"2026-09-17T04:01:51","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/hypothese-egalite-des-variances-statistique\/"},"modified":"2026-09-17T04:01:55","modified_gmt":"2026-09-17T04:01:55","slug":"hypothese-egalite-des-variances-statistique","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=hypothese-egalite-des-variances-statistique","title":{"rendered":"Qu&rsquo;est-ce que l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances en statistique ?"},"content":{"rendered":"<p>L&rsquo;inf\u00e9rence statistique param\u00e9trique repose sur un ensemble de postulats th\u00e9oriques fondamentaux qui garantissent la validit\u00e9 math\u00e9matique des conclusions tir\u00e9es \u00e0 partir d&rsquo;\u00e9chantillons empiriques. Parmi ces conditions pr\u00e9alables, l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances \u2014 couramment d\u00e9sign\u00e9e sous les termes d&rsquo;homosc\u00e9dasticit\u00e9 ou d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances selon les cadres de mod\u00e9lisation \u2014 occupe une place \u00e0 la fois centrale et fr\u00e9quemment sous-estim\u00e9e. Trop souvent rel\u00e9gu\u00e9e au statut de simple case m\u00e9thodologique \u00e0 cocher lors de la v\u00e9rification routini\u00e8re des pr\u00e9requis, cette hypoth\u00e8se conditionne pourtant directement la distribution exacte des statistiques de test fondamentales, \u00e0 l&rsquo;instar des ratios de <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Student%27s_t-test\">Student<\/a> et de Fisher-Snedecor. Lorsque cette condition est satisfaite, les estimateurs de la dispersion r\u00e9siduelle convergent harmonieusement vers la variance populationnelle th\u00e9orique, permettant un contr\u00f4le rigoureux du taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce et garantissant une puissance statistique optimale.<\/p>\n<p>Cependant, dans la pratique de la recherche scientifique \u2014 tout particuli\u00e8rement au sein des sciences du comportement, de la psychologie quantitative, des neurosciences cognitives et de la m\u00e9decine clinique \u2014, les donn\u00e9es empiriques se conforment rarement aux canons de l&rsquo;homog\u00e9n\u00e9it\u00e9 parfaite. Les processus psychologiques, d\u00e9veloppementaux ou psychopathologiques induisent intrins\u00e8quement des asym\u00e9tries de variabilit\u00e9 : les groupes cliniques manifestent fr\u00e9quemment une dispersion accrue par rapport aux cohortes saines, les phases d&rsquo;apprentissage g\u00e9n\u00e8rent une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 temporelle de la performance, et les instruments psychom\u00e9triques r\u00e9v\u00e8lent des effets d&rsquo;\u00e9chelle qui modulent artificiellement la variabilit\u00e9 des scores. Dans ces contextes, postuler na\u00efvement que diff\u00e9rents groupes ou diff\u00e9rents niveaux de pr\u00e9dicteurs partagent une variance commune rel\u00e8ve d&rsquo;une simplification abusive susceptible de corrompre gravement l&rsquo;\u00e9difice inf\u00e9rentiel.<\/p>\n<p>L&rsquo;objet de cet article encyclop\u00e9dique est de diss\u00e9quer de mani\u00e8re exhaustive la nature, les fondements math\u00e9matiques, les implications pratiques et les solutions contemporaines relatives \u00e0 l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances. En explorant ses manifestations depuis le classique test <em>t<\/em> de comparaison de deux moyennes jusqu&rsquo;aux raffinements de la r\u00e9gression lin\u00e9aire multiple et des mod\u00e8les mixtes, nous mettrons en lumi\u00e8re les cons\u00e9quences concr\u00e8tes de sa violation sur l&rsquo;inflation du risque alpha et la d\u00e9gradation de la puissance d&rsquo;analyse. Enfin, nous proposerons un cadre d\u00e9cisionnel op\u00e9rationnel et rigoureux, guid\u00e9 par les recommandations contemporaines de la science ouverte et de la mod\u00e9lisation statistique robuste, afin de permettre aux chercheurs et analystes de traiter l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances non pas comme un obstacle paralysant, mais comme une information substantielle sur leurs objets d&rsquo;\u00e9tude.<\/p>\n<h2>1. Cadre conceptuel et d\u00e9finition de l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances<\/h2>\n<h3>1.1 D\u00e9finition formelle de l&rsquo;homosc\u00e9dasticit\u00e9 et de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9<\/h3>\n<p>Sur le plan s\u00e9mantique et m\u00e9thodologique, les termes \u00ab homog\u00e9n\u00e9it\u00e9 des variances \u00bb et \u00ab homosc\u00e9dasticit\u00e9 \u00bb d\u00e9signent fondamentalement le m\u00eame principe sous-jacent, mais leur emploi varie usuellement selon la nature du paradigme analytique mobilis\u00e9. Le vocable \u00ab homog\u00e9n\u00e9it\u00e9 des variances \u00bb est traditionnellement privil\u00e9gi\u00e9 dans le cadre des plans exp\u00e9rimentaux discrets, tels que les tests de comparaison de moyennes d&rsquo;\u00e9chantillons ind\u00e9pendants et les analyses de variance univari\u00e9es ou multivari\u00e9es (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Analysis_of_variance\">ANOVA<\/a>). \u00c0 l&rsquo;inverse, le terme \u00ab homosc\u00e9dasticit\u00e9 \u00bb \u2014 issu des racines grecques <em>homo<\/em> (m\u00eame) et <em>skedannumi<\/em> (dispersion) \u2014 trouve son terrain de pr\u00e9dilection dans la mod\u00e9lisation lin\u00e9aire g\u00e9n\u00e9rale, la r\u00e9gression par les moindres carr\u00e9s ordinaires (OLS) et l&rsquo;\u00e9conom\u00e9trie, o\u00f9 la dispersion des r\u00e9sidus est \u00e9valu\u00e9e le long d&rsquo;un continuum de valeurs pr\u00e9dictives continues.<\/p>\n<p>Sur le plan formel, dans une configuration impliquant <em>k<\/em> populations ind\u00e9pendantes, l&rsquo;hypoth\u00e8se nulle d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances postule l&rsquo;\u00e9quivalence stricte des variances populationnelles respectives. Math\u00e9matiquement, cela se formule de la mani\u00e8re suivante :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\sigma_1^2 = \\sigma_2^2 = dots = \\sigma_k^2 = \\sigma^2\">$$\\sigma_1^2 = \\sigma_2^2 = dots = \\sigma_k^2 = \\sigma^2$$<\/span><\/p>\n<p>Dans ce mod\u00e8le, chaque sous-population <em>j<\/em> (pour <em>j<\/em> allant de 1 \u00e0 <em>k<\/em>) poss\u00e8de une distribution th\u00e9orique dont le second moment centr\u00e9 est rigoureusement identique \u00e0 une constante commune <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span>, bien que leurs premiers moments respectifs (les moyennes populationnelles <span class=\"arabpsy-math-inline\" data-latex=\"\\mu_j\">$\\mu_j$<\/span>) puissent diverger sous l&rsquo;hypoth\u00e8se alternative de l&rsquo;investigateur.<\/p>\n<p>Dans le contexte de la r\u00e9gression lin\u00e9aire et de la mod\u00e9lisation conditionnelle, l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 se d\u00e9finit comme la non-constance de la variance conditionnelle des termes d&rsquo;erreur stochastiques <span class=\"arabpsy-math-inline\" data-latex=\"\\epsilon_i\">$\\epsilon_i$<\/span> sachant le vecteur de r\u00e9gresseurs <span class=\"arabpsy-math-inline\" data-latex=\"X_i\">$X_i$<\/span>. Formellement, la condition d&rsquo;homosc\u00e9dasticit\u00e9 s&rsquo;\u00e9crit :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Var}(\\epsilon_i mid X_i) = \\sigma^2 \\quad \\forall i in {1, dots, N}\">$$\\text{Var}(\\epsilon_i mid X_i) = \\sigma^2 \\quad \\forall i in {1, dots, N}$$<\/span><\/p>\n<p>D\u00e8s lors que cette variance conditionnelle n&rsquo;est plus une constante scalaire mais devient une fonction d\u00e9pendante des observations, des niveaux d&rsquo;un pr\u00e9dicteur ou de variables exog\u00e8nes latentes \u2014 de sorte que <span class=\"arabpsy-math-inline\" data-latex=\"\\text{Var}(\\epsilon_i mid X_i) = \\sigma_i^2\">$\\text{Var}(\\epsilon_i mid X_i) = \\sigma_i^2$<\/span> \u2014, le mod\u00e8le bascule dans l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9. Ce postulat constitue l&rsquo;une des pierres angulaires de l&rsquo;inf\u00e9rence param\u00e9trique classique, car son respect garantit que le bruit de mesure ou la variabilit\u00e9 non expliqu\u00e9e op\u00e8re avec une intensit\u00e9 homog\u00e8ne sur l&rsquo;ensemble du domaine d&rsquo;\u00e9tude, autorisant ainsi une estimation globale et unifi\u00e9e de l&rsquo;incertitude.<\/p>\n<h3>1.2 Origine et r\u00f4le de l&rsquo;hypoth\u00e8se dans la th\u00e9orie statistique classique<\/h3>\n<p>L&rsquo;origine de l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances est intrins\u00e8quement li\u00e9e \u00e0 la gen\u00e8se m\u00eame de la th\u00e9orie de l&rsquo;\u00e9chantillonnage et de l&rsquo;estimation optimale d\u00e9velopp\u00e9e au cours du XIXe et du d\u00e9but du XXe si\u00e8cle par des figures tut\u00e9laires telles que Carl Friedrich Gauss, William Sealy Gosset (sous le pseudonyme de Student) et Sir Ronald Aylmer Fisher. Dans la conceptualisation de la m\u00e9thode des moindres carr\u00e9s, Gauss a d\u00e9montr\u00e9 que la combinaison lin\u00e9aire optimale des observations n\u00e9cessitait l&rsquo;hypoth\u00e8se que chaque r\u00e9sidu partage\u00e2t la m\u00eame dispersion d&rsquo;erreur, condition sine qua non pour d\u00e9river la variance minimale des estimateurs.<\/p>\n<p>Lorsque Gosset formalise la distribution <em>t<\/em> en 1908 pour r\u00e9soudre le probl\u00e8me de l&rsquo;inf\u00e9rence sur de petits \u00e9chantillons dont la variance est inconnue, la mutualisation des variances devient un imp\u00e9ratif de calcul. Pour estimer l&rsquo;erreur-type de la diff\u00e9rence entre deux moyennes empiriques sans supposer connue la variance de la population, il est math\u00e9matiquement indispensable de supposer une variance d&rsquo;erreur unique. Cette \u00e9tape permet d&rsquo;agr\u00e9ger ou de regrouper (op\u00e9ration dite de <em>pooling<\/em>) les sommes de carr\u00e9s des \u00e9carts des deux \u00e9chantillons afin d&rsquo;obtenir un estimateur quadratique moyen commun non biais\u00e9, bas\u00e9 sur un nombre de degr\u00e9s de libert\u00e9 accru (<span class=\"arabpsy-math-inline\" data-latex=\"n_1 + n_2 - 2\">$n_1 + n_2 &#8211; 2$<\/span>).<\/p>\n<p>Fisher, en \u00e9tendant ces travaux \u00e0 la comparaison simultan\u00e9e de multiples groupes \u00e0 travers l&rsquo;invention de l&rsquo;analyse de variance, a institutionnalis\u00e9 cette n\u00e9cessit\u00e9. La d\u00e9rivation exacte de la distribution <em>F<\/em> de Fisher-Snedecor repose express\u00e9ment sur le ratio de deux variables al\u00e9atoires ind\u00e9pendantes suivant chacune une loi du Khi-carr\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"\\chi^2\">$\\chi^2$<\/span>), divis\u00e9es par leurs degr\u00e9s de libert\u00e9 respectifs. Cette d\u00e9rivation math\u00e9matique \u00e9l\u00e9gante n&rsquo;est rigoureusement exacte que si le num\u00e9rateur (variance intergroupe) et le d\u00e9nominateur (variance intragroupe ou carr\u00e9 moyen r\u00e9siduel) mesurent tous deux la m\u00eame quantit\u00e9 fondamentale <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span> sous l&rsquo;hypoth\u00e8se nulle. Si les variances sous-jacentes des groupes divergent, le d\u00e9nominateur ne suit plus une distribution du Khi-carr\u00e9 standard divis\u00e9e par ses degr\u00e9s de libert\u00e9, ce qui fausse les propri\u00e9t\u00e9s distributionnelles exactes de la statistique de test et invalide les tables de probabilit\u00e9 associ\u00e9es.<\/p>\n<h3>1.3 Pertinence \u00e9pist\u00e9mologique et contextuelle en psychologie quantitative<\/h3>\n<p>Dans le domaine des sciences psychologiques, cognitives et comportementales, l&rsquo;hypoth\u00e8se d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances acquiert une dimension \u00e9pist\u00e9mologique cruciale qui d\u00e9passe la simple m\u00e9canique math\u00e9matique. Les donn\u00e9es recueillies dans les recherches quantitatives en psychologie s&rsquo;\u00e9cartent tr\u00e8s fr\u00e9quemment des id\u00e9aux gaussiens th\u00e9oriques. Les \u00e9chelles d&rsquo;\u00e9valuation, qu&rsquo;il s&rsquo;agisse de questionnaires de personnalit\u00e9, de mesures de bien-\u00eatre subjectif ou de scores psychom\u00e9triques de d\u00e9pression, g\u00e9n\u00e8rent couramment des distributions marqu\u00e9es par une asym\u00e9trie substantielle et une variabilit\u00e9 h\u00e9t\u00e9rog\u00e8ne inh\u00e9rente aux ph\u00e9nom\u00e8nes mesur\u00e9s.<\/p>\n<p>Cette variabilit\u00e9 n&rsquo;est pas un simple artefact statistique, mais refl\u00e8te fr\u00e9quemment une r\u00e9alit\u00e9 substantielle : les populations cliniques (par exemple, des individus pr\u00e9sentant un trouble du spectre de l&rsquo;autisme, un \u00e9pisode d\u00e9pressif caract\u00e9ris\u00e9 ou une l\u00e9sion c\u00e9r\u00e9brale focale) manifestent typiquement une dispersion interindividuelle consid\u00e9rablement plus \u00e9lev\u00e9e que les populations t\u00e9moins neurotypiques ou asymptomatiques. Dans un groupe t\u00e9moin, les m\u00e9canismes de r\u00e9gulation ou les comp\u00e9tences cognitives de base tendent \u00e0 concentrer les scores autour d&rsquo;une moyenne normative, g\u00e9n\u00e9rant une variance restreinte. Dans le groupe clinique, l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des \u00e9tiologies, la s\u00e9v\u00e9rit\u00e9 variable des sympt\u00f4mes, les comorbidit\u00e9s et les m\u00e9canismes compensatoires individuels provoquent un \u00e9talement majeur de la distribution des performances, mat\u00e9rialis\u00e9 par une variance substantiellement accrue.<\/p>\n<p>Dans le cadre de la psychom\u00e9trie et de la validation des instruments de mesure, ignorer ces disparit\u00e9s de variance expose les chercheurs \u00e0 des erreurs diagnostiques majeures. Par exemple, lors de la normalisation d&rsquo;un test ou de l&rsquo;\u00e9valuation de l&rsquo;invariance de mesure \u00e0 travers diff\u00e9rents sous-groupes d\u00e9mographiques (tels que le genre, les tranches d&rsquo;\u00e2ge ou le niveau socio-\u00e9conomique), la pr\u00e9sence d&rsquo;une h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 non contr\u00f4l\u00e9e fausse les comparaisons de scores latents et conduit \u00e0 rejeter \u00e0 tort l&rsquo;\u00e9quivalence des construits mesur\u00e9s. De surcro\u00eet, les effets d&rsquo;intervention psychoth\u00e9rapeutique ou pharmacologique modifient souvent non seulement la moyenne du groupe trait\u00e9, mais \u00e9galement sa dispersion \u2014 certains patients r\u00e9pondant magnifiquement tandis que d&rsquo;autres stagnent \u2014, rompant ainsi le postulat d&rsquo;\u00e9galit\u00e9 indispensable aux analyses conventionnelles non ajust\u00e9es.<\/p>\n<h2>2. L&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances dans le test t de Student<\/h2>\n<h3>2.1 M\u00e9canisme math\u00e9matique de la variance combin\u00e9e<\/h3>\n<p>Le test <em>t<\/em> de Student pour deux \u00e9chantillons ind\u00e9pendants est con\u00e7u pour \u00e9valuer l&rsquo;hypoth\u00e8se nulle selon laquelle les moyennes de deux populations distinctes sont identiques (<span class=\"arabpsy-math-inline\" data-latex=\"H_0: \\mu_1 = \\mu_2\">$H_0: \\mu_1 = \\mu_2$<\/span>). Le calcul classique de la statistique <em>t<\/em> d\u00e9pend intimement du concept de variance combin\u00e9e ou pond\u00e9r\u00e9e, d\u00e9sign\u00e9e dans la litt\u00e9rature sous l&rsquo;anglicisme de <em>pooled variance<\/em> (<span class=\"arabpsy-math-inline\" data-latex=\"s_p^2\">$s_p^2$<\/span>). Ce m\u00e9canisme repose sur le postulat formel selon lequel les deux groupes sont issus de populations partageant une m\u00eame variance th\u00e9orique <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span>, les variances d&rsquo;\u00e9chantillonnage <span class=\"arabpsy-math-inline\" data-latex=\"s_1^2\">$s_1^2$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"s_2^2\">$s_2^2$<\/span> n&rsquo;\u00e9tant que des estimations fluctuantes de cette quantit\u00e9 unique.<\/p>\n<p>La formule math\u00e9matique de la variance combin\u00e9e s&rsquo;exprime comme une moyenne des variances d&rsquo;\u00e9chantillons, pond\u00e9r\u00e9e par leurs degr\u00e9s de libert\u00e9 respectifs :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"s_p^2 = \\frac{(n_1 - 1)s_1^2 + (n_2 - 1)s_2^2}{n_1 + n_2 - 2}\">$$s_p^2 = \\frac{(n_1 &#8211; 1)s_1^2 + (n_2 &#8211; 1)s_2^2}{n_1 + n_2 &#8211; 2}$$<\/span><\/p>\n<p>Dans cette \u00e9quation, <span class=\"arabpsy-math-inline\" data-latex=\"n_1\">$n_1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"n_2\">$n_2$<\/span> repr\u00e9sentent les effectifs des deux \u00e9chantillons, tandis que <span class=\"arabpsy-math-inline\" data-latex=\"n_1 + n_2 - 2\">$n_1 + n_2 &#8211; 2$<\/span> correspond aux degr\u00e9s de libert\u00e9 totaux allou\u00e9s \u00e0 l&rsquo;estimation de l&rsquo;erreur r\u00e9siduelle conjointe. L&rsquo;erreur-type de la diff\u00e9rence entre les deux moyennes empiriques (<span class=\"arabpsy-math-inline\" data-latex=\"\\bar{X}_1 - \\bar{X}_2\">$\\bar{X}_1 &#8211; \\bar{X}_2$<\/span>) est ensuite d\u00e9riv\u00e9e directement de cette variance combin\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"SE_{(\\bar{X}_1 - \\bar{X}_2)} = \\sqrt{s_p^2 \\left(\\frac{1}{n_1} + \\frac{1}{n_2}\\right)}\">$$SE_{(\\bar{X}_1 &#8211; \\bar{X}_2)} = \\sqrt{s_p^2 \\left(\\frac{1}{n_1} + \\frac{1}{n_2}\\right)}$$<\/span><\/p>\n<p>La statistique <em>t<\/em> finale prend alors la forme canonique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"t = \\frac{\\bar{X}_1 - \\bar{X}_2}{s_p \\sqrt{\\frac{1}{n_1} + \\frac{1}{n_2}}}\">$$t = \\frac{\\bar{X}_1 &#8211; \\bar{X}_2}{s_p \\sqrt{\\frac{1}{n_1} + \\frac{1}{n_2}}}$$<\/span><\/p>\n<p>La condition absolue de validit\u00e9 de cette statistique est que le ratio suive fid\u00e8lement une distribution th\u00e9orique de Student \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"n_1 + n_2 - 2\">$n_1 + n_2 &#8211; 2$<\/span> degr\u00e9s de libert\u00e9 sous <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span>. Cette convergence distributionnelle n&rsquo;est garantie que si <span class=\"arabpsy-math-inline\" data-latex=\"s_p^2\">$s_p^2$<\/span> est un estimateur sans biais et efficace de l&rsquo;erreur stochastique commune. Si la condition <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_1^2 = \\sigma_2^2\">$\\sigma_1^2 = \\sigma_2^2$<\/span> est compromise, <span class=\"arabpsy-math-inline\" data-latex=\"s_p^2\">$s_p^2$<\/span> cesse de repr\u00e9senter ad\u00e9quatement l&rsquo;incertitude sp\u00e9cifique associ\u00e9e \u00e0 chaque groupe, ce qui d\u00e9stabilise le d\u00e9nominateur de la statistique <em>t<\/em> et alt\u00e8re sa distribution de r\u00e9f\u00e9rence.<\/p>\n<h3>2.2 Impact des tailles d&rsquo;\u00e9chantillons sur la robustesse du test t<\/h3>\n<p>L&rsquo;impact de la violation de l&rsquo;homog\u00e9n\u00e9it\u00e9 des variances sur le test <em>t<\/em> classique d\u00e9pend de mani\u00e8re critique de l&rsquo;\u00e9quilibre des effectifs entre les deux groupes \u00e9tudi\u00e9s. Il s&rsquo;agit de l&rsquo;un des ph\u00e9nom\u00e8nes les plus remarquables de la th\u00e9orie des tests statistiques : l&rsquo;existence d&rsquo;une interaction majeure entre la disparit\u00e9 des variances et le rapport des tailles d&rsquo;\u00e9chantillons (<span class=\"arabpsy-math-inline\" data-latex=\"n_1 \/ n_2\">$n_1 \/ n_2$<\/span>).<\/p>\n<p>Lorsque le plan exp\u00e9rimental est rigoureusement \u00e9quilibr\u00e9, c&rsquo;est-\u00e0-dire lorsque <span class=\"arabpsy-math-inline\" data-latex=\"n_1 = n_2\">$n_1 = n_2$<\/span>, le test <em>t<\/em> de Student traditionnel d\u00e9montre une remarquable robustesse intrins\u00e8que face aux entorses \u00e0 l&rsquo;\u00e9galit\u00e9 des variances. Dans ce cas sp\u00e9cifique, la formule de la variance combin\u00e9e attribue un poids strictement \u00e9gal aux deux variances d&rsquo;\u00e9chantillon. M\u00eame si une population pr\u00e9sente une variance deux, trois ou quatre fois sup\u00e9rieure \u00e0 l&rsquo;autre, le taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce r\u00e9el (le taux <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha\">$\\alpha$<\/span> empirique) demeure remarquablement proche du seuil nominal conventionnel de 5 % (variant typiquement entre 0,045 et 0,065 selon l&rsquo;ampleur du d\u00e9s\u00e9quilibre de variance), \u00e0 condition que la distribution des donn\u00e9es ne pr\u00e9sente pas d&rsquo;asym\u00e9trie s\u00e9v\u00e8re.<\/p>\n<p>En revanche, d\u00e8s lors que le plan exp\u00e9rimental devient d\u00e9s\u00e9quilibr\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"n_1 \\neq n_2\">$n_1 \\neq n_2$<\/span>), la robustesse du test <em>t<\/em> classique s&rsquo;effondre litt\u00e9ralement, donnant lieu \u00e0 deux sc\u00e9narios hautement probl\u00e9matiques :<\/p>\n<ul>\n<li><strong>Le cas de l&rsquo;appariement positif :<\/strong> La taille d&rsquo;\u00e9chantillon la plus grande est associ\u00e9e \u00e0 la population pr\u00e9sentant la plus grande variance (<span class=\"arabpsy-math-inline\" data-latex=\"n_{\\text{grand}} \\rightarrow \\sigma_{\\text{grand}}^2\">$n_{\\text{grand}} \\rightarrow \\sigma_{\\text{grand}}^2$<\/span>). Dans cette situation, la formule de pooling surestime largement la variance r\u00e9siduelle globale en la tirant vers la valeur du grand groupe. L&rsquo;erreur-type conjointe est alors artificiellement gonfl\u00e9e, ce qui r\u00e9duit drastiquement la statistique <em>t<\/em> calcul\u00e9e. Le test devient hyper-conservateur : le risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce chute bien en de\u00e7\u00e0 du seuil nominal (souvent en dessous de 1 %), au prix d&rsquo;une perte catastrophique de puissance statistique (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Type_I_and_type_II_errors\">erreur de type II<\/a>).<\/li>\n<li><strong>Le cas de l&rsquo;appariement n\u00e9gatif :<\/strong> La taille d&rsquo;\u00e9chantillon la plus faible est associ\u00e9e \u00e0 la population pr\u00e9sentant la plus grande variance (<span class=\"arabpsy-math-inline\" data-latex=\"n_{\\text{petit}} \\rightarrow \\sigma_{\\text{grand}}^2\">$n_{\\text{petit}} \\rightarrow \\sigma_{\\text{grand}}^2$<\/span>). Il s&rsquo;agit du sc\u00e9nario le plus destructeur pour l&rsquo;int\u00e9grit\u00e9 scientifique. La variance combin\u00e9e est principalement d\u00e9termin\u00e9e par le grand \u00e9chantillon \u00e0 faible variance ; elle sous-estime tragiquement l&rsquo;incertitude r\u00e9elle pesant sur la diff\u00e9rence de moyennes. Le d\u00e9nominateur de la statistique <em>t<\/em> est anormalement contract\u00e9, catapultant la statistique <em>t<\/em> vers des valeurs extr\u00eames. Le taux de faux positifs s&rsquo;envole alors de mani\u00e8re dramatique, atteignant fr\u00e9quemment 15 %, 25 % voire 40 % pour un seuil nominal fix\u00e9 \u00e0 5 %. Le chercheur conclut ainsi \u00e0 une d\u00e9couverte statistiquement significative qui ne refl\u00e8te qu&rsquo;un artefact d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9.<\/li>\n<\/ul>\n<h3>2.3 Distinction entre test unilat\u00e9ral et bilat\u00e9ral face \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9<\/h3>\n<p>L&rsquo;asym\u00e9trie induite par des variances in\u00e9gales n&rsquo;exerce pas les m\u00eames ravages selon que le chercheur formule une hypoth\u00e8se directionnelle (test unilat\u00e9ral) ou non directionnelle (test bilat\u00e9ral). Dans un test bilat\u00e9ral, le risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce nominal de 5 % est r\u00e9parti \u00e9quitablement aux deux extr\u00e9mit\u00e9s de la distribution th\u00e9orique (2,5 % dans chaque queue). Bien que le taux d&rsquo;erreur global cumul\u00e9 soit d\u00e9form\u00e9 en cas d&rsquo;appariement n\u00e9gatif ou positif, la sym\u00e9trie de la r\u00e8gle de d\u00e9cision bilat\u00e9rale permet de temp\u00e9rer marginalement certaines distorsions unilat\u00e9rales extr\u00eames.<\/p>\n<p>\u00c0 l&rsquo;inverse, lors de la mise en \u0153uvre d&rsquo;un test unilat\u00e9ral \u2014 d\u00e9marche courante lorsqu&rsquo;un clinicien postule a priori qu&rsquo;un traitement exp\u00e9rimental ne peut que surpasser un placebo \u2014, la sensibilit\u00e9 \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 combin\u00e9e \u00e0 la forme des queues de distribution s&rsquo;accro\u00eet consid\u00e9rablement. Si la population affect\u00e9e par la plus forte variabilit\u00e9 correspond au groupe pour lequel on attend une moyenne plus \u00e9lev\u00e9e, la queue de distribution empirique situ\u00e9e du c\u00f4t\u00e9 de la r\u00e9gion de rejet devient anormalement \u00e9paisse. Il en r\u00e9sulte que le seuil critique unilat\u00e9ral est franchi avec une probabilit\u00e9 encore plus d\u00e9s\u00e9quilibr\u00e9e que dans le sc\u00e9nario bilat\u00e9ral.<\/p>\n<p>Les simulations statistiques men\u00e9es en m\u00e9thodologie psychologique d\u00e9montrent sans ambigu\u00eft\u00e9 que les tests directionnels unilat\u00e9raux appliqu\u00e9s \u00e0 des donn\u00e9es h\u00e9t\u00e9rosc\u00e9dastiques issues d&rsquo;\u00e9chantillons in\u00e9gaux g\u00e9n\u00e8rent des taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce totalement incontr\u00f4lables, pouvant d\u00e9passer les 30 % dans des contextes pourtant fr\u00e9quents en laboratoire (par exemple un groupe contr\u00f4le de 50 participants oppos\u00e9 \u00e0 un groupe exp\u00e9rimental clinique restreint de 20 participants). Par cons\u00e9quent, la recommandation m\u00e9thodologique moderne impose de proscrire rigoureusement le recours au test <em>t<\/em> de Student classique sous formulation unilat\u00e9rale d\u00e8s lors que l&rsquo;\u00e9galit\u00e9 absolue des variances ne peut \u00eatre d\u00e9montr\u00e9e formellement.<\/p>\n<h2>3. Le postulat d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances dans l&rsquo;ANOVA \u00e0 un ou plusieurs facteurs<\/h2>\n<h3>3.1 D\u00e9composition de la variance et calcul du ratio F<\/h3>\n<p>L&rsquo;analyse de variance \u00e0 un facteur (ANOVA <em>one-way<\/em>) g\u00e9n\u00e9ralise la logique du test <em>t<\/em> \u00e0 la comparaison simultan\u00e9e de trois groupes ou davantage. Son principe algorithmique repose sur le th\u00e9or\u00e8me de d\u00e9composition de la somme des carr\u00e9s : la variation totale observ\u00e9e dans les donn\u00e9es (<span class=\"arabpsy-math-inline\" data-latex=\"SS_{\\text{totale}}\">$SS_{\\text{totale}}$<\/span>) est scind\u00e9e de fa\u00e7on orthogonale en une variation attribuable aux diff\u00e9rences entre les moyennes des groupes (<span class=\"arabpsy-math-inline\" data-latex=\"SS_{\\text{inter}}\">$SS_{\\text{inter}}$<\/span> ou <em>Between-groups sum of squares<\/em>) et une variation attribuable aux fluctuations individuelles au sein de chaque groupe (<span class=\"arabpsy-math-inline\" data-latex=\"SS_{\\text{intra}}\">$SS_{\\text{intra}}$<\/span> ou <em>Within-groups sum of squares<\/em>, fr\u00e9quemment d\u00e9nomm\u00e9e somme des carr\u00e9s r\u00e9siduels, <span class=\"arabpsy-math-inline\" data-latex=\"SS_{\\text{r\u00e9siduelle}}\">$SS_{\\text{r\u00e9siduelle}}$<\/span>).<\/p>\n<p>Chacune de ces composantes est divis\u00e9e par ses degr\u00e9s de libert\u00e9 respectifs (<span class=\"arabpsy-math-inline\" data-latex=\"k - 1\">$k &#8211; 1$<\/span> pour la composante intergroupe, et <span class=\"arabpsy-math-inline\" data-latex=\"N - k\">$N &#8211; k$<\/span> pour la composante intragroupe) afin d&rsquo;obtenir des carr\u00e9s moyens. Le carr\u00e9 moyen r\u00e9siduel (<span class=\"arabpsy-math-inline\" data-latex=\"MS_{\\text{r\u00e9siduel}}\">$MS_{\\text{r\u00e9siduel}}$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"MS_{\\text{intra}}\">$MS_{\\text{intra}}$<\/span>) constitue la pi\u00e8ce ma\u00eetresse du dispositif : il repr\u00e9sente la moyenne arithm\u00e9tique pond\u00e9r\u00e9e des variances d&rsquo;\u00e9chantillon de l&rsquo;ensemble des <em>k<\/em> groupes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"MS_{\\text{r\u00e9siduel}} = \\frac{\\sum_{j=1}^k (n_j - 1)s_j^2}{N - k}\">$$MS_{\\text{r\u00e9siduel}} = \\frac{\\sum_{j=1}^k (n_j &#8211; 1)s_j^2}{N &#8211; k}$$<\/span><\/p>\n<p>La statistique inf\u00e9rentielle calcul\u00e9e est le ratio <em>F<\/em> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"F = \\frac{MS_{\\text{inter}}}{MS_{\\text{r\u00e9siduel}}}\">$$F = \\frac{MS_{\\text{inter}}}{MS_{\\text{r\u00e9siduel}}}$$<\/span><\/p>\n<p>Pour que cette statistique <em>F<\/em> suive exactement la distribution th\u00e9orique de Fisher-Snedecor avec les param\u00e8tres <span class=\"arabpsy-math-inline\" data-latex=\"k-1\">$k-1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"N-k\">$N-k$<\/span>, il est imp\u00e9ratif que les erreurs au sein de chaque niveau de traitement soient des variables al\u00e9atoires ind\u00e9pendantes et identiquement distribu\u00e9es suivant une loi normale d&rsquo;esp\u00e9rance nulle et de variance identique <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span>. En cas de violation de l&rsquo;homog\u00e9n\u00e9it\u00e9 des variances (<span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_1^2 \\neq \\sigma_2^2 \\neq dots \\neq \\sigma_k^2\">$\\sigma_1^2 \\neq \\sigma_2^2 \\neq dots \\neq \\sigma_k^2$<\/span>), le carr\u00e9 moyen r\u00e9siduel cesse de suivre une loi du Khi-carr\u00e9 standardis\u00e9e unifi\u00e9e. Le num\u00e9rateur et le d\u00e9nominateur cessent de satisfaire aux crit\u00e8res th\u00e9oriques d&rsquo;ind\u00e9pendance et d&rsquo;\u00e9chelles de dispersion conjointe. Il s&rsquo;ensuit une distorsion majeure de la fonction de densit\u00e9 de probabilit\u00e9 du ratio <em>F<\/em>, rendant les <em>p<\/em>-valeurs tabul\u00e9es hautement fallacieuses.<\/p>\n<h3>3.2 Plans factoriels et interactions complexes<\/h3>\n<p>La complexit\u00e9 de l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances s&rsquo;amplifie consid\u00e9rablement lorsque l&rsquo;on quitte le cadre univari\u00e9 \u00e0 un facteur pour explorer des plans factoriels crois\u00e9s (plans de type <span class=\"arabpsy-math-inline\" data-latex=\"2 \\times 2\">$2 \\times 2$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"3 \\times 2\">$3 \\times 2$<\/span>, ou <span class=\"arabpsy-math-inline\" data-latex=\"2 \\times 2 \\times 2\">$2 \\times 2 \\times 2$<\/span>). Dans ces dispositifs, l&rsquo;attention du chercheur ne se focalise pas uniquement sur les effets principaux des diff\u00e9rents facteurs, mais sonde de mani\u00e8re prioritaire l&rsquo;existence de termes d&rsquo;interaction, c&rsquo;est-\u00e0-dire l&rsquo;\u00e9valuation de la mani\u00e8re dont l&rsquo;effet d&rsquo;une variable ind\u00e9pendante se module en fonction des modalit\u00e9s d&rsquo;une autre.<\/p>\n<p>Dans les plans factoriels d\u00e9s\u00e9quilibr\u00e9s \u2014 o\u00f9 les cellules exp\u00e9rimentales comportent des effectifs in\u00e9gaux \u2014, l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 affecte de fa\u00e7on dissym\u00e9trique les effets principaux et les termes d&rsquo;interaction. Des simulations approfondies ont mis en exergue que les tests portant sur les interactions sont fr\u00e9quemment bien plus vuln\u00e9rables aux distorsions d&rsquo;\u00e9galit\u00e9 des variances que les tests d&rsquo;effets principaux. Une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 localis\u00e9e dans une seule cellule exp\u00e9rimentale (par exemple, un sous-groupe combinant un niveau d&rsquo;anxi\u00e9t\u00e9 \u00e9lev\u00e9 et une privation de sommeil s\u00e9v\u00e8re manifestant un \u00e9croulement de l&rsquo;homog\u00e9n\u00e9it\u00e9 de sa performance) est capable de contaminer l&rsquo;ensemble du terme d&rsquo;erreur r\u00e9siduel commun de l&rsquo;ANOVA factorielle.<\/p>\n<p>Cette contamination induit une cascade d&rsquo;erreurs m\u00e9thodologiques :<\/p>\n<ul>\n<li>Des interactions hautement significatives peuvent \u00e9merger de fa\u00e7on totalement fallacieuse par simple art\u00e9fact d&rsquo;une variance d\u00e9mesur\u00e9e dans une condition exp\u00e9rimentale particuli\u00e8re.<\/li>\n<li>Inversement, de v\u00e9ritables interactions th\u00e9oriques peuvent \u00eatre enti\u00e8rement masqu\u00e9es en raison d&rsquo;un gonflement artificiel du d\u00e9nominateur <span class=\"arabpsy-math-inline\" data-latex=\"MS_{\\text{r\u00e9siduel}}\">$MS_{\\text{r\u00e9siduel}}$<\/span>, induit par des cellules p\u00e9riph\u00e9riques instables.<\/li>\n<li>Le recours habituel aux d\u00e9compositions des sommes de carr\u00e9s de Type II ou de Type III ne pallie en rien cette d\u00e9ficience, car ces ajustements ne corrigent que les d\u00e9pendances structurelles li\u00e9es au d\u00e9s\u00e9quilibre des effectifs sans r\u00e9parer la rupture d&rsquo;homosc\u00e9dasticit\u00e9 sous-jacente.<\/li>\n<\/ul>\n<h3>3.3 Application sp\u00e9cifique aux mesures r\u00e9p\u00e9t\u00e9es et sph\u00e9ricit\u00e9<\/h3>\n<p>Lorsque les donn\u00e9es recueillies reposent sur des plans intra-participants ou longitudinaux, dans lesquels les m\u00eames sujets sont mesur\u00e9s \u00e0 travers diff\u00e9rentes conditions exp\u00e9rimentales ou points temporels, le concept d&rsquo;\u00e9galit\u00e9 des variances mute pour adopter une forme multidimensionnelle plus exigeante : la <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Mauchly%27s_sphericity_test\">sph\u00e9ricit\u00e9<\/a>. Dans l&rsquo;ANOVA \u00e0 mesures r\u00e9p\u00e9t\u00e9es, l&rsquo;hypoth\u00e8se d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances univari\u00e9e est insuffisante en raison de l&rsquo;interd\u00e9pendance structurelle des observations issues d&rsquo;un m\u00eame individu.<\/p>\n<p>La condition de sph\u00e9ricit\u00e9 \u2014 formalis\u00e9e math\u00e9matiquement par l&rsquo;hypoth\u00e8se de Mauchly \u2014 exige l&rsquo;\u00e9galit\u00e9 absolue des variances de toutes les diff\u00e9rences possibles entre les paires de conditions intra-sujets. Soient trois conditions r\u00e9p\u00e9t\u00e9es <span class=\"arabpsy-math-inline\" data-latex=\"A\">$A$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"C\">$C$<\/span>. L&rsquo;hypoth\u00e8se de sph\u00e9ricit\u00e9 impose la condition stricte suivante :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Var}(Y_A - Y_B) = \\text{Var}(Y_A - Y_C) = \\text{Var}(Y_B - Y_C) = \\text{constante}\">$$\\text{Var}(Y_A &#8211; Y_B) = \\text{Var}(Y_A &#8211; Y_C) = \\text{Var}(Y_B &#8211; Y_C) = \\text{constante}$$<\/span><\/p>\n<p>La sph\u00e9ricit\u00e9 repr\u00e9sente une condition moins restrictive que celle de sym\u00e9trie compos\u00e9e (qui requerrait \u00e0 la fois l&rsquo;\u00e9galit\u00e9 exacte de toutes les variances individuelles des conditions et l&rsquo;\u00e9galit\u00e9 parfaite de toutes les covariances entre paires de conditions). N\u00e9anmoins, dans la pratique empirique en psychologie et en neuro-imagerie, la sph\u00e9ricit\u00e9 est quasi syst\u00e9matiquement viol\u00e9e d\u00e8s lors que le nombre de niveaux du facteur r\u00e9p\u00e9t\u00e9 d\u00e9passe deux. Les points temporels plus rapproch\u00e9s manifestent presque toujours une corr\u00e9lation sup\u00e9rieure \u00e0 celle unissant des points \u00e9loign\u00e9s, ce qui fait in\u00e9vitablement fluctuer la variance des diff\u00e9rences.<\/p>\n<p>La violation de la sph\u00e9ricit\u00e9 a pour cons\u00e9quence in\u00e9luctable une inflation substantielle du risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce pour les tests <em>F<\/em> intra-sujets. L&rsquo;analyste se voit contraint d&rsquo;abandonner la statistique <em>F<\/em> usuelle au profit de versions p\u00e9nalis\u00e9es r\u00e9duisant math\u00e9matiquement les degr\u00e9s de libert\u00e9 du num\u00e9rateur et du d\u00e9nominateur via le facteur d&rsquo;ajustement epsilon (<span class=\"arabpsy-math-inline\" data-latex=\"epsilon\">$epsilon$<\/span>), tel que propos\u00e9 par Greenhouse et Geisser (proc\u00e9dure tr\u00e8s conservatrice) ou par Huynh et Feldt (proc\u00e9dure plus lib\u00e9rale adapt\u00e9e aux valeurs d&rsquo;epsilon proches de l&rsquo;unit\u00e9).<\/p>\n<h2>4. L&rsquo;hypoth\u00e8se d&rsquo;homosc\u00e9dasticit\u00e9 dans la r\u00e9gression lin\u00e9aire ordinaire (OLS)<\/h2>\n<h3>4.1 Le th\u00e9or\u00e8me de Gauss-Markov et les propri\u00e9t\u00e9s BLUE<\/h3>\n<p>Dans l&rsquo;architecture de la r\u00e9gression lin\u00e9aire univari\u00e9e ou multiple estim\u00e9e par la m\u00e9thode des moindres carr\u00e9s ordinaires (<em>Ordinary Least Squares<\/em> ou OLS), l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances des r\u00e9sidus constitue l&rsquo;un des piliers du v\u00e9n\u00e9rable <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Gauss%E2%80%93Markov_theorem\">th\u00e9or\u00e8me de Gauss-Markov<\/a>. Ce th\u00e9or\u00e8me fondamental \u00e9tablit les conditions sous lesquelles les estimateurs OLS des coefficients de r\u00e9gression <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}\">$\\hat{\\beta}$<\/span> sont qualifi\u00e9s de BLUE (<em>Best Linear Unbiased Estimator<\/em>), c&rsquo;est-\u00e0-dire les meilleurs estimateurs lin\u00e9aires sans biais.<\/p>\n<p>Pour qu&rsquo;un estimateur m\u00e9rite l&rsquo;attribut \u00ab <em>Best<\/em> \u00bb (meilleur), il doit pr\u00e9senter la variance d&rsquo;\u00e9chantillonnage minimale au sein de l&rsquo;ensemble de la classe des estimateurs lin\u00e9aires sans biais. Le th\u00e9or\u00e8me de Gauss-Markov prouve math\u00e9matiquement que cette efficacit\u00e9 optimale n&rsquo;est garantie que si la matrice de variance-covariance conditionnelle des erreurs stochastiques <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{\\epsilon}\">$\\mathbf{\\epsilon}$<\/span> est une matrice scalaire unitaire :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Var}(\\mathbf{\\epsilon} mid \\mathbf{X}) = \\sigma^2 \\mathbf{I}_N\">$$\\text{Var}(\\mathbf{\\epsilon} mid \\mathbf{X}) = \\sigma^2 \\mathbf{I}_N$$<\/span><\/p>\n<p>Cette condition implique deux r\u00e9quisits distincts : l&rsquo;absence d&rsquo;autocorr\u00e9lation entre les r\u00e9sidus (les \u00e9l\u00e9ments extradiagonaux de la matrice sont nuls) et l&rsquo;homosc\u00e9dasticit\u00e9 parfaite (tous les \u00e9l\u00e9ments de la diagonale principale sont rigoureusement \u00e9gaux \u00e0 une constante unique <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span>).<\/p>\n<p>Il est capital de comprendre une nuance th\u00e9orique majeure souvent mal comprise : la violation du postulat d&rsquo;homosc\u00e9dasticit\u00e9 n&rsquo;induit aucun biais dans l&rsquo;estimation ponctuelle des coefficients de r\u00e9gression. Les pentes estim\u00e9es <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}_1, \\hat{\\beta}_2, dots, \\hat{\\beta}_p\">$\\hat{\\beta}_1, \\hat{\\beta}_2, dots, \\hat{\\beta}_p$<\/span> demeurent math\u00e9matiquement non biais\u00e9es : leur esp\u00e9rance reste rigoureusement \u00e9gale aux vrais param\u00e8tres populationnels <span class=\"arabpsy-math-inline\" data-latex=\"\\beta\">$\\beta$<\/span>, sous r\u00e9serve que l&rsquo;hypoth\u00e8se d&rsquo;exog\u00e9n\u00e9it\u00e9 stricte (<span class=\"arabpsy-math-inline\" data-latex=\"E(\\mathbf{\\epsilon} mid \\mathbf{X}) = 0\">$E(\\mathbf{\\epsilon} mid \\mathbf{X}) = 0$<\/span>) soit satisfaite. Cependant, l&rsquo;estimateur OLS perd sa propri\u00e9t\u00e9 d&rsquo;efficacit\u00e9 statistique maximale : il cesse d&rsquo;\u00eatre BLUE. Il existe alors d&rsquo;autres familles d&rsquo;estimateurs lin\u00e9aires, notamment les moindres carr\u00e9s pond\u00e9r\u00e9s ou g\u00e9n\u00e9ralis\u00e9s, qui poss\u00e8dent une variance d&rsquo;\u00e9chantillonnage strictement inf\u00e9rieure.<\/p>\n<h3>4.2 Impact sur l&rsquo;estimation des erreurs-types des coefficients<\/h3>\n<p>Si la non-biaisitude des pentes est pr\u00e9serv\u00e9e en pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9, la m\u00e9canique inf\u00e9rentielle entourant ces coefficients est, quant \u00e0 elle, totalement compromise. Dans le cadre OLS classique, la formule analytique consacr\u00e9e pour estimer la matrice de variance-covariance des coefficients estim\u00e9s repose express\u00e9ment sur l&rsquo;hypoth\u00e8se d&rsquo;homog\u00e9n\u00e9it\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Var}(\\hat{\\mathbf{\\beta}}) = \\sigma^2 (\\mathbf{X}&#039;\\mathbf{X})^{-1}\">$$\\text{Var}(\\hat{\\mathbf{\\beta}}) = \\sigma^2 (\\mathbf{X}&rsquo;\\mathbf{X})^{-1}$$<\/span><\/p>\n<p>O\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span> est estim\u00e9 sans distinction par la variance r\u00e9siduelle globale <span class=\"arabpsy-math-inline\" data-latex=\"s_e^2 = \\frac{\\mathbf{e}&#039;\\mathbf{e}}{N - p - 1}\">$s_e^2 = \\frac{\\mathbf{e}&rsquo;\\mathbf{e}}{N &#8211; p &#8211; 1}$<\/span>. Or, si la v\u00e9ritable matrice de variance-covariance des perturbations est h\u00e9t\u00e9rosc\u00e9dastique et prend la forme <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{\\Omega} = \\text{diag}(\\sigma_1^2, \\sigma_2^2, dots, \\sigma_N^2)\">$\\mathbf{\\Omega} = \\text{diag}(\\sigma_1^2, \\sigma_2^2, dots, \\sigma_N^2)$<\/span>, la structure exacte de dispersion des estimateurs r\u00e9pond en r\u00e9alit\u00e9 \u00e0 la formulation matricielle dite \u00ab en sandwich \u00bb :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Var}(\\hat{\\mathbf{\\beta}})_{\\text{r\u00e9elle}} = (\\mathbf{X}&#039;\\mathbf{X})^{-1} (\\mathbf{X}&#039; \\mathbf{\\Omega} \\mathbf{X}) (\\mathbf{X}&#039;\\mathbf{X})^{-1}\">$$\\text{Var}(\\hat{\\mathbf{\\beta}})_{\\text{r\u00e9elle}} = (\\mathbf{X}&rsquo;\\mathbf{X})^{-1} (\\mathbf{X}&rsquo; \\mathbf{\\Omega} \\mathbf{X}) (\\mathbf{X}&rsquo;\\mathbf{X})^{-1}$$<\/span><\/p>\n<p>En employant aveugl\u00e9ment la formule standard au lieu de l&rsquo;expression exacte, le logiciel d&rsquo;analyse calcule des erreurs-types (<em>standard errors<\/em>, <span class=\"arabpsy-math-inline\" data-latex=\"SE(\\hat{\\beta}_j)\">$SE(\\hat{\\beta}_j)$<\/span>) qui sont biais\u00e9es de mani\u00e8re syst\u00e9matique. Dans la grande majorit\u00e9 des applications en sciences sociales, cette distorsion se traduit par une sous-estimation substantielle des v\u00e9ritables erreurs-types. En sous-estimant l&rsquo;erreur-type, le ratio de Wald classique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"t = \\frac{\\hat{\\beta}_j}{SE(\\hat{\\beta}_j)}\">$$t = \\frac{\\hat{\\beta}_j}{SE(\\hat{\\beta}_j)}$$<\/span><\/p>\n<p>est artificiellement gonfl\u00e9 vers le haut, \u00e9crasant les <em>p<\/em>-valeurs et produisant des r\u00e9sultats fallacieusement significatifs. De plus, les intervalles de confiance \u00e0 95 % calcul\u00e9s autour des pentes de r\u00e9gression deviennent excessivement \u00e9troits : ils \u00e9chouent \u00e0 capturer la vraie valeur du param\u00e8tre dans 95 % des r\u00e9plications th\u00e9oriques, affichant parfois des taux de couverture effectifs chutant sous la barre des 80 %.<\/p>\n<h3>4.3 Formes courantes d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 en mod\u00e9lisation psychologique<\/h3>\n<p>L&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 en mod\u00e9lisation psychologique et comportementale n&rsquo;appara\u00eet que rarement sous la forme d&rsquo;un bruit blanc al\u00e9atoire ; elle s&rsquo;organise le plus souvent selon des structures morphologiques conditionn\u00e9es par la nature des processus \u00e9tudi\u00e9s. On r\u00e9pertorie trois configurations embl\u00e9matiques particuli\u00e8rement r\u00e9pandues dans la discipline :<\/p>\n<ul>\n<li><strong>L&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 monotone croissante (ou d\u00e9croissante) li\u00e9e aux apprentissages ou \u00e0 la fatigue :<\/strong> Dans les protocoles chronom\u00e9triques d&rsquo;apprentissage moteur ou cognitif, la dispersion des temps de r\u00e9action est massive au d\u00e9but de l&rsquo;exp\u00e9rimentation lorsque les participants explorent des strat\u00e9gies divergentes. Au fur et \u00e0 mesure des blocs d&rsquo;essais, les performances s&rsquo;automatisent et convergent, provoquant une contraction progressive et continue de la variance r\u00e9siduelle \u00e0 mesure que les valeurs pr\u00e9dites diminuent. Inversement, lors d&rsquo;\u00e9preuves d&rsquo;endurance cognitive prolong\u00e9e, l&rsquo;apparition asynchrone de la fatigue chez certains sujets g\u00e9n\u00e8re un \u00e9clatement de la variance r\u00e9siduelle pour les valeurs \u00e9lev\u00e9es de dur\u00e9e.<\/li>\n<li><strong>L&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 structurelle li\u00e9e aux gradients d\u00e9veloppementaux ou d\u00e9mographiques :<\/strong> En psychologie du d\u00e9veloppement et en g\u00e9rontologie cognitive, la variabilit\u00e9 interindividuelle des capacit\u00e9s cognitives (m\u00e9moire de travail, vitesse d&rsquo;inhibition) augmente de fa\u00e7on structurelle avec l&rsquo;\u00e2ge chronologique. L&rsquo;inclusion de participants s&rsquo;\u00e9talant de l&rsquo;enfance \u00e0 l&rsquo;\u00e2ge adulte avanc\u00e9 induit syst\u00e9matiquement une variance r\u00e9siduelle conditionnelle qui cro\u00eet le long de l&rsquo;axe des \u00e2ges, dessinant une dispersion en \u00e9ventail classique.<\/li>\n<li><strong>Les effets de seuil (plancher et plafond) des instruments psychom\u00e9triques :<\/strong> Lorsqu&rsquo;une \u00e9chelle d&rsquo;auto-\u00e9valuation ou un test cognitif pr\u00e9sente une difficult\u00e9 inadapt\u00e9e \u00e0 la population \u00e9chantillonn\u00e9e, une fraction majeure des participants s&rsquo;agglutine au score minimal ou maximal possible. Vers ces extr\u00eames de l&rsquo;\u00e9chelle de mesure, la dispersion des scores est m\u00e9caniquement comprim\u00e9e \u00e0 z\u00e9ro par la borne de l&rsquo;instrument, alors qu&rsquo;elle s&rsquo;exprime librement au milieu de l&rsquo;\u00e9chelle. Ce ph\u00e9nom\u00e8ne engendre une h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 en forme de fuseau ou de cloche invers\u00e9e qui brise totalement l&rsquo;hypoth\u00e8se de variance constante.<\/li>\n<\/ul>\n<h2>5. \u00c9valuation visuelle et exploratoire de l&rsquo;\u00e9galit\u00e9 des variances<\/h2>\n<h3>5.1 Diagrammes en bo\u00eete (boxplots) et repr\u00e9sentations par dispersion<\/h3>\n<p>La premi\u00e8re ligne de d\u00e9fense de l&rsquo;analyste face aux \u00e9carts d&rsquo;homog\u00e9n\u00e9it\u00e9 r\u00e9side dans l&rsquo;analyse exploratoire et la visualisation graphique rigoureuse des donn\u00e9es brutes. Pour les configurations d&rsquo;\u00e9chantillons ind\u00e9pendants ou d&rsquo;ANOVA, le diagramme en bo\u00eete de Tukey (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Box_plot\">boxplot<\/a>) constitue l&rsquo;outil de criblage universel par excellence. L&rsquo;examen attentif d&rsquo;un boxplot offre une appr\u00e9ciation imm\u00e9diate et robuste de la dispersion sans d\u00e9pendre d&rsquo;hypoth\u00e8ses distributionnelles restrictives.<\/p>\n<p>L&rsquo;\u00e9valuation repose principalement sur la hauteur de la bo\u00eete centrale, qui mat\u00e9rialise l&rsquo;\u00e9cart interquartile (IQR, d\u00e9limit\u00e9 par le 25e et le 75e percentile et englobant 50 % des donn\u00e9es centrales). Lorsque les hauteurs des bo\u00eetes diff\u00e8rent visuellement de mani\u00e8re flagrante d&rsquo;un groupe \u00e0 l&rsquo;autre \u2014 par exemple si la bo\u00eete d&rsquo;un groupe clinique est deux ou trois fois plus haute que celle du groupe t\u00e9moin \u2014, l&rsquo;existence d&rsquo;une violation de l&rsquo;\u00e9galit\u00e9 des variances est hautement probable. Dans un second temps, l&rsquo;analyste examine l&rsquo;envergure totale des moustaches (qui s&rsquo;\u00e9tendent typiquement jusqu&rsquo;\u00e0 1,5 fois l&rsquo;IQR \u00e0 partir des charni\u00e8res) ainsi que le nombre et la distribution des valeurs consid\u00e9r\u00e9es comme aberrantes ou atypiques (<em>outliers<\/em>).<\/p>\n<p>Le boxplot permet \u00e9galement de dissocier simultan\u00e9ment la question de la dispersion et celle de la sym\u00e9trie de la distribution : la position relative de la ligne m\u00e9diane au sein de la bo\u00eete et l&rsquo;\u00e9galit\u00e9 de longueur des moustaches sup\u00e9rieure et inf\u00e9rieure renseignent sur l&rsquo;asym\u00e9trie (<em>skewness<\/em>). Pour un raffinement sup\u00e9rieur, la superposition d&rsquo;un nuage de points individualis\u00e9s avec instabilit\u00e9 horizontale al\u00e9atoire (<em>jittered stripchart<\/em>) ou le recours aux diagrammes en violon (<em>violin plots<\/em>) combinant l&rsquo;estimation de densit\u00e9 par noyau offre une transparence absolue sur l&rsquo;\u00e9talement des donn\u00e9es r\u00e9elles.<\/p>\n<h3>5.2 Graphiques des r\u00e9sidus versus valeurs pr\u00e9dites (Residuals vs Fitted)<\/h3>\n<p>Dans le cadre de la mod\u00e9lisation lin\u00e9aire et de la r\u00e9gression OLS, le graphique de diagnostic standard incontournable confronte les r\u00e9sidus ordinaires ou standardis\u00e9s (<span class=\"arabpsy-math-inline\" data-latex=\"e_i\">$e_i$<\/span>) port\u00e9s sur l&rsquo;axe des ordonn\u00e9es aux valeurs pr\u00e9dites ou ajust\u00e9es (<span class=\"arabpsy-math-inline\" data-latex=\"\\hat{Y}_i\">$\\hat{Y}_i$<\/span>) port\u00e9es sur l&rsquo;axe des abscisses. Lorsque l&rsquo;hypoth\u00e8se d&rsquo;homosc\u00e9dasticit\u00e9 est scrupuleusement respect\u00e9e dans la population \u00e9tudi\u00e9e, ce graphique affiche un nuage de points sans motif particulier, formant une bande horizontale r\u00e9guli\u00e8re d&rsquo;\u00e9paisseur constante tout au long de l&rsquo;axe des ordonn\u00e9es, r\u00e9partie de mani\u00e8re sym\u00e9trique et homog\u00e8ne de part et d&rsquo;autre de la ligne horizontale horizontale <span class=\"arabpsy-math-inline\" data-latex=\"y = 0\">$y = 0$<\/span>.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/boxplot-1.jpg\" alt=\"\" \/><\/figure>\n<p>L&rsquo;expression graphique embl\u00e9matique d&rsquo;une violation de l&rsquo;homosc\u00e9dasticit\u00e9 r\u00e9side dans l&rsquo;apparition d&rsquo;un motif caract\u00e9ristique en c\u00f4ne ou en entonnoir (commun\u00e9ment qualifi\u00e9 de <em>fan shape<\/em> ou <em>cone shape<\/em>). Dans cette configuration visuelle, les r\u00e9sidus sont tr\u00e8s regroup\u00e9s autour de z\u00e9ro pour les faibles valeurs de <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{Y}\">$\\hat{Y}$<\/span>, puis s&rsquo;\u00e9cartent de plus en plus vigoureusement \u00e0 mesure que <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{Y}\">$\\hat{Y}$<\/span> progresse vers la droite, trahissant une variance d&rsquo;erreur qui cro\u00eet en fonction directe de la pr\u00e9diction du mod\u00e8le.<\/p>\n<p>Ce graphique joue \u00e9galement un r\u00f4le capital de discrimination clinique entre deux violations distinctes : la non-lin\u00e9arit\u00e9 et l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9. Une tendance curviligne (en forme de parabole, de U ou de S) trac\u00e9e par les points ou par une courbe d&rsquo;ajustement non param\u00e9trique locale (LOESS) signale une d\u00e9faillance de la sp\u00e9cification de la fonction moyenne (non-lin\u00e9arit\u00e9), qui peut coexister avec l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des r\u00e9sidus ou en \u00eatre la cause directe. L&rsquo;isolation visuelle d&rsquo;une simple d\u00e9rive de l&rsquo;\u00e9paisseur du nuage de points sans courbure de la ligne centrale confirme quant \u00e0 elle une h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 pure.<\/p>\n<h3>5.3 Graphiques Scale-Location (Spread-Location)<\/h3>\n<p>Bien que le graphique des r\u00e9sidus bruts versus valeurs ajust\u00e9es soit tr\u00e8s populaire, il pr\u00e9sente une limite perceptuelle humaine bien document\u00e9e : l&rsquo;\u0153il \u00e9prouve des difficult\u00e9s \u00e0 comparer rigoureusement la largeur d&rsquo;un nuage de points lorsque la densit\u00e9 des observations varie le long de l&rsquo;axe ou lorsque des asym\u00e9tries unilat\u00e9rales sont pr\u00e9sentes. Pour surmonter cet \u00e9cueil cognitif, le graphique <em>Scale-Location<\/em> (ou <em>Spread-Location<\/em> plot) propose une transformation \u00e9l\u00e9gante de l&rsquo;axe vertical.<\/p>\n<p>L&rsquo;ordonn\u00e9e de ce graphique affiche la racine carr\u00e9e de la valeur absolue des r\u00e9sidus standardis\u00e9s :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"y_i = \\sqrt{|\\tilde{e}_i|}\">$$y_i = \\sqrt{|\\tilde{e}_i|}$$<\/span><\/p>\n<p>Cette transformation math\u00e9matique produit deux effets cruciaux. D&rsquo;une part, la prise de valeur absolue rabat l&rsquo;ensemble des r\u00e9sidus dans le quadrant positif, \u00e9liminant la composante directionnelle du signe pour ne conserver que la pure magnitude de l&rsquo;\u00e9cart. D&rsquo;autre part, l&rsquo;application de la racine carr\u00e9e temp\u00e8re la courbure naturelle induite par les distributions des carr\u00e9s et stabilise la variance de l&rsquo;indicateur lorsque les r\u00e9sidus d&rsquo;origine sont normaux. Sur ce graphique, l&rsquo;hypoth\u00e8se d&rsquo;homosc\u00e9dasticit\u00e9 se traduit par une ligne de lissage non param\u00e9trique (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Local_regression\">LOESS<\/a>) rigoureusement horizontale et plane.<\/p>\n<p>Si la ligne rouge LOESS d\u00e9vie de l&rsquo;horizontale en adoptant une trajectoire ascendante continue, l&rsquo;analyste dispose d&rsquo;une preuve visuelle \u00e9clatante d&rsquo;une augmentation de la dispersion des r\u00e9sidus corr\u00e9l\u00e9e aux valeurs ajust\u00e9es. Si la ligne plonge aux extr\u00e9mit\u00e9s ou forme une onde, cela r\u00e9v\u00e8le des zones localis\u00e9es de sur-dispersion ou de sous-dispersion qui guideront le choix d&rsquo;interventions cibl\u00e9es.<\/p>\n<h2>6. Tests d&rsquo;hypoth\u00e8se formels pour l&rsquo;\u00e9valuation de l&rsquo;\u00e9galit\u00e9 des variances<\/h2>\n<h3>6.1 Le test de Levene et ses variantes robustes (Brown-Forsythe)<\/h3>\n<p>Face aux subjectivit\u00e9s inh\u00e9rentes \u00e0 l&rsquo;interpr\u00e9tation purement visuelle des graphiques, la communaut\u00e9 statistique a \u00e9labor\u00e9 des tests d&rsquo;hypoth\u00e8se formels destin\u00e9s \u00e0 trancher statistiquement sur l&rsquo;\u00e9galit\u00e9 des variances. Le test historiquement le plus pr\u00e9conis\u00e9 dans les progiciels g\u00e9n\u00e9ralistes (SPSS, SAS) est le test de Levene (1960). L&rsquo;intuition fondamentale d&rsquo;Howard Levene a \u00e9t\u00e9 de transformer le probl\u00e8me de l&rsquo;\u00e9valuation de l&rsquo;\u00e9galit\u00e9 des variances en une classique analyse de variance \u00e0 un facteur appliqu\u00e9e non pas aux scores bruts, mais \u00e0 une m\u00e9trique d&rsquo;\u00e9cart individuel absolu.<\/p>\n<p>Dans sa version originelle, le test de Levene calcule pour chaque individu <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> appartenant au groupe <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span> la valeur absolue de l&rsquo;\u00e9cart \u00e0 la moyenne de son groupe :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"z_{ij} = |Y_{ij} - \\bar{Y}_{\\cdot j}|\">$$z_{ij} = |Y_{ij} &#8211; \\bar{Y}_{\\cdot j}|$$<\/span><\/p>\n<p>L&rsquo;analyse proc\u00e8de ensuite au calcul d&rsquo;une ANOVA standard sur ces scores d\u00e9riv\u00e9s <span class=\"arabpsy-math-inline\" data-latex=\"z_{ij}\">$z_{ij}$<\/span>. Si les variances des groupes d&rsquo;origine sont \u00e9gales, l&rsquo;esp\u00e9rance de ces distances absolues est identique entre les groupes, et le test <em>F<\/em> qui en r\u00e9sulte n&rsquo;est pas significatif. Cependant, la formulation initiale de Levene s&rsquo;est r\u00e9v\u00e9l\u00e9e hautement vuln\u00e9rable aux violations de l&rsquo;hypoth\u00e8se de normalit\u00e9 : si les donn\u00e9es brutes sont asym\u00e9triques ou pr\u00e9sentent des queues \u00e9paisses (kurtosis prononc\u00e9e), le test de Levene moyen g\u00e9n\u00e8re un exc\u00e8s dramatique d&rsquo;erreurs de premi\u00e8re esp\u00e8ce.<\/p>\n<p>Pour rem\u00e9dier \u00e0 cette tare, Brown et Forsythe (1974) ont introduit une modification cruciale : substituer la m\u00e9diane du groupe (<span class=\"arabpsy-math-inline\" data-latex=\"\\tilde{Y}_{\\cdot j}\">$\\tilde{Y}_{\\cdot j}$<\/span>) \u00e0 la moyenne arithm\u00e9tique lors du calcul des \u00e9carts absolus :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"z_{ij}^* = |Y_{ij} - \\tilde{Y}_{\\cdot j}|\">$$z_{ij}^* = |Y_{ij} &#8211; \\tilde{Y}_{\\cdot j}|$$<\/span><\/p>\n<p>Le test de Brown-Forsythe qui en d\u00e9coule conf\u00e8re une robustesse exceptionnelle face aux distributions asym\u00e9triques et aux contaminations par des valeurs extr\u00eames. L&rsquo;interpr\u00e9tation de la <em>p<\/em>-valeur suit la logique standard de l&rsquo;inf\u00e9rence : un r\u00e9sultat significatif (typiquement <span class=\"arabpsy-math-inline\" data-latex=\"p &lt; 0,05\">$p &lt; 0,05$<\/span>) conduit au rejet de l&#039;hypoth\u00e8se nulle d&#039;homog\u00e9n\u00e9it\u00e9 des variances, intimant \u00e0 l&#039;exp\u00e9rimentateur d&#039;adopter des m\u00e9thodes correctrices.<\/p>\n<h3>6.2 Le test de Bartlett et le test de Fligner-Killeen<\/h3>\n<p>Ant\u00e9rieurement aux travaux de Levene, Maurice Bartlett avait propos\u00e9 en 1937 un test d&rsquo;\u00e9galit\u00e9 des variances d\u00e9riv\u00e9 de la th\u00e9orie du rapport de vraisemblance (<em>likelihood ratio test<\/em>). Sous l&rsquo;hypoth\u00e8se formelle de normalit\u00e9 stricte des donn\u00e9es dans chaque population parente, la statistique de Bartlett s&rsquo;exprime comme suit :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"K = \\frac{(N - k) \\ln(s_p^2) - \\sum_{j=1}^k (n_j - 1) \\ln(s_j^2)}{1 + \\frac{1}{3(k - 1)} \\left( \\sum_{j=1}^k \\left( \\frac{1}{n_j - 1} \\right) - \\frac{1}{N - k} \\right)}\">$$K = \\frac{(N &#8211; k) \\ln(s_p^2) &#8211; \\sum_{j=1}^k (n_j &#8211; 1) \\ln(s_j^2)}{1 + \\frac{1}{3(k &#8211; 1)} \\left( \\sum_{j=1}^k \\left( \\frac{1}{n_j &#8211; 1} \\right) &#8211; \\frac{1}{N &#8211; k} \\right)}$$<\/span><\/p>\n<p>Sous <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span>, la statistique <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span> suit asymptotiquement une loi du Khi-carr\u00e9 \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"k - 1\">$k &#8211; 1$<\/span> degr\u00e9s de libert\u00e9. Lorsque la condition de normalit\u00e9 parfaite est satisfaite, le test de Bartlett est le plus puissant disponible sur le plan th\u00e9orique. Cependant, cette puissance exceptionnelle se paie au prix d&rsquo;une fragilit\u00e9 devenue proverbiale dans la litt\u00e9rature statistique. Le test de Bartlett est hypersensible aux moindres d\u00e9viations par rapport \u00e0 la loi normale univari\u00e9e : George Box a synth\u00e9tis\u00e9 cette faiblesse en \u00e9crivant avec ironie que \u00ab faire un test de Bartlett pour v\u00e9rifier l&rsquo;\u00e9galit\u00e9 des variances avant une ANOVA revient \u00e0 prendre la mer dans un canot de sauvetage pour v\u00e9rifier si les conditions m\u00e9t\u00e9orologiques sont assez s\u00fbres pour un paquebot \u00bb.<\/p>\n<p>\u00c0 l&rsquo;autre extr\u00e9mit\u00e9 du spectre se trouve le test de Fligner-Killeen (1976). Il s&rsquo;agit d&rsquo;une proc\u00e9dure non param\u00e9trique totalement affranchie de l&rsquo;hypoth\u00e8se de normalit\u00e9, bas\u00e9e sur le classement par rangs des scores centr\u00e9s par rapport \u00e0 leurs m\u00e9dianes de groupe, auxquels on applique ensuite une transformation bas\u00e9e sur les quantiles de la loi normale. Le test de Fligner-Killeen s&rsquo;av\u00e8re l&rsquo;un des outils les plus fiables et puissants lorsque les distributions sous-jacentes sont lourdement non-gaussiennes, bi-modales ou pollu\u00e9es par des valeurs aberrantes massives.<\/p>\n<h3>6.3 Limites inh\u00e9rentes \u00e0 l&rsquo;utilisation exclusive des tests formels<\/h3>\n<p>Bien que les tests statistiques formels d&rsquo;homog\u00e9n\u00e9it\u00e9 offrent un confort psychologique certain aux chercheurs en fournissant un verdict binaire rassurant (<span class=\"arabpsy-math-inline\" data-latex=\"p &gt; 0,05\">$p &gt; 0,05$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"p &lt; 0,05\">$p &lt; 0,05$<\/span>), la communaut\u00e9 m\u00e9thodologique contemporaine s&#039;accorde \u00e0 condamner leur usage exclusif et non critique. Deux faiblesses statistiques fondamentales expliquent cette prise de position unanime :<\/p>\n<ul>\n<li><strong>Le probl\u00e8me de la d\u00e9pendance \u00e0 la taille d&rsquo;\u00e9chantillon :<\/strong> Comme tout test de signification de l&rsquo;hypoth\u00e8se nulle (NHST), les tests de Levene, Bartlett ou Brown-Forsythe voient leur puissance statistique fluctuer drastiquement en fonction directe de la taille totale de l&rsquo;\u00e9chantillon (<span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span>). Dans les \u00e9tudes \u00e0 tr\u00e8s grands effectifs (par exemple en \u00e9pid\u00e9miologie ou dans les cohortes massives de psychologie en ligne avec <span class=\"arabpsy-math-inline\" data-latex=\"N &gt; 2000\">$N &gt; 2000$<\/span>), le test formel d\u00e9tectera des divergences microscopiques et totalement triviales de variance entre les groupes (qui n&rsquo;auraient strictement aucun impact d\u00e9l\u00e9t\u00e8re sur le test param\u00e9trique sous-jacent) et affichera des <em>p<\/em>-valeurs hyper-significatives (<span class=\"arabpsy-math-inline\" data-latex=\"p  0,05\">$p  0,05$<\/span>), ber\u00e7ant l&rsquo;exp\u00e9rimentateur dans une fausse illusion de s\u00e9curit\u00e9.<\/li>\n<li><strong>Le paradoxe du pr\u00e9-test (analyse en deux \u00e9tapes) :<\/strong> Conditionner l&rsquo;analyse principale (choisir entre le test <em>t<\/em> de Student et le test de Welch) sur le r\u00e9sultat formel d&rsquo;un pr\u00e9-test de Levene modifie et pollue le taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce global du protocole analytique. Cette strat\u00e9gie s\u00e9quentielle en cascade introduit une d\u00e9pendance statistique complexe qui alt\u00e8re les propri\u00e9t\u00e9s inf\u00e9rentielles de l&rsquo;ensemble de la d\u00e9marche.<\/li>\n<\/ul>\n<h2>7. Cons\u00e9quences statistiques de la violation du postulat<\/h2>\n<h3>7.1 Inflation ou d\u00e9flation du taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce (Type I)<\/h3>\n<p>L&rsquo;impact le plus redoutable de la violation de l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances r\u00e9side dans l&rsquo;alt\u00e9ration profonde du taux nominal d&rsquo;erreur de premi\u00e8re esp\u00e8ce (<span class=\"arabpsy-math-inline\" data-latex=\"\\alpha\">$\\alpha$<\/span>), traditionnellement calibr\u00e9 \u00e0 la valeur consensuelle de 5 % (<span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,05\">$\\alpha = 0,05$<\/span>). Dans un article s\u00e9minal publi\u00e9 en 1954, George E. P. Box a mod\u00e9lis\u00e9 math\u00e9matiquement les distorsions exactes subies par les tests inf\u00e9rentiels param\u00e9triques lorsque les postulats sous-jacents sont bris\u00e9s. Ses d\u00e9rivations analytiques, corrobor\u00e9es depuis par des millions de simulations de Monte-Carlo, ont formalis\u00e9 la dichotomie fondamentale entre inflation et d\u00e9flation du risque d&rsquo;erreur selon la structure de la covariance des donn\u00e9es.<\/p>\n<p>Comme \u00e9voqu\u00e9 pr\u00e9c\u00e9demment, la distorsion prend des proportions alarmantes dans les plans d\u00e9s\u00e9quilibr\u00e9s. Consid\u00e9rons un cas classique en recherche exp\u00e9rimentale : un chercheur compare un groupe t\u00e9moin volumineux (<span class=\"arabpsy-math-inline\" data-latex=\"n_1 = 80\">$n_1 = 80$<\/span>) et un groupe clinique rare (<span class=\"arabpsy-math-inline\" data-latex=\"n_2 = 20\">$n_2 = 20$<\/span>). Si le groupe clinique manifeste une variabilit\u00e9 quatre fois sup\u00e9rieure \u00e0 celle du groupe t\u00e9moin (<span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_2^2 = 4\\sigma_1^2\">$\\sigma_2^2 = 4\\sigma_1^2$<\/span>, soit un ratio d&rsquo;\u00e9carts-types de 2:1), et que le chercheur utilise le test <em>t<\/em> de Student conventionnel ou l&rsquo;ANOVA classique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"s_p^2 \\approx \\frac{79 s_1^2 + 19 (4s_1^2)}{98} = \\frac{155}{98} s_1^2 \\approx 1,58 s_1^2\">$$s_p^2 \\approx \\frac{79 s_1^2 + 19 (4s_1^2)}{98} = \\frac{155}{98} s_1^2 \\approx 1,58 s_1^2$$<\/span><\/p>\n<p>Le calcul montre que la variance combin\u00e9e est principalement r\u00e9gie par la variance du grand \u00e9chantillon. Or, l&rsquo;erreur-type r\u00e9elle de la diff\u00e9rence fait peser un poids massif sur le terme <span class=\"arabpsy-math-inline\" data-latex=\"s_2^2\/n_2 = 4s_1^2\/20 = 0,20 s_1^2\">$s_2^2\/n_2 = 4s_1^2\/20 = 0,20 s_1^2$<\/span>, alors que la formule combin\u00e9e lui attribue artificiellement <span class=\"arabpsy-math-inline\" data-latex=\"s_p^2\/n_2 \\approx 1,58 s_1^2 \/ 20 \\approx 0,079 s_1^2\">$s_p^2\/n_2 \\approx 1,58 s_1^2 \/ 20 \\approx 0,079 s_1^2$<\/span>. L&rsquo;erreur-type est sous-estim\u00e9e de plus de la moiti\u00e9. Le taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce r\u00e9el explose pour atteindre approximativement 18 % \u00e0 22 % au lieu des 5 % suppos\u00e9s.<\/p>\n<p>Dans un tel contexte, pr\u00e8s d&rsquo;une d\u00e9cision de rejet de l&rsquo;hypoth\u00e8se nulle sur cinq constitue une illusion statistique (un faux positif pur). Ce ph\u00e9nom\u00e8ne a constitu\u00e9 l&rsquo;un des moteurs silencieux mais destructeurs de la crise de la reproductibilit\u00e9 qui a frapp\u00e9 la psychologie quantitative au cours de la derni\u00e8re d\u00e9cennie.<\/p>\n<h3>7.2 Perte de puissance statistique et erreur de seconde esp\u00e8ce (Type II)<\/h3>\n<p>\u00c0 l&rsquo;autre extr\u00e9mit\u00e9 du spectre d&rsquo;interaction, le sc\u00e9nario de l&rsquo;appariement positif (<span class=\"arabpsy-math-inline\" data-latex=\"n_{\\text{grand}} \\rightarrow \\sigma_{\\text{grand}}^2\">$n_{\\text{grand}} \\rightarrow \\sigma_{\\text{grand}}^2$<\/span>) engendre une pathologie m\u00e9thodologique tout aussi pernicieuse, bien que moins m\u00e9diatis\u00e9e que l&rsquo;inflation des faux positifs : l&rsquo;effondrement de la puissance statistique (<span class=\"arabpsy-math-inline\" data-latex=\"1 - \\beta\">$1 &#8211; \\beta$<\/span>) et l&rsquo;explosion corr\u00e9lative de l&rsquo;<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Type_I_and_type_II_errors#Type_II_error\">erreur de type II<\/a> (taux <span class=\"arabpsy-math-inline\" data-latex=\"\\beta\">$\\beta$<\/span>).<\/p>\n<p>Dans cette configuration, l&rsquo;attribution d&rsquo;un poids pr\u00e9pond\u00e9rant \u00e0 la forte variance du grand \u00e9chantillon au sein de la variance combin\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"s_p^2\">$s_p^2$<\/span> entra\u00eene un gonflement massif du d\u00e9nominateur du test <em>t<\/em> ou du carr\u00e9 moyen r\u00e9siduel de l&rsquo;ANOVA. En cons\u00e9quence directe, le test statistique devient excessivement conservateur, son taux <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha\">$\\alpha$<\/span> effectif s&rsquo;effondrant parfois sous la barre de 0,5 %. D\u00e8s lors, m\u00eame en pr\u00e9sence d&rsquo;un effet exp\u00e9rimental r\u00e9el, cliniquement ou th\u00e9oriquement substantiel dans la population cible, le ratio empirique \u00e9choue syst\u00e9matiquement \u00e0 d\u00e9passer la valeur critique tabul\u00e9e.<\/p>\n<p>Ce dysfonctionnement alt\u00e8re fondamentalement l&rsquo;estimation de la taille d&rsquo;effet standardis\u00e9e. Le <span class=\"arabpsy-math-inline\" data-latex=\"d\">$d$<\/span> de Cohen classique, calcul\u00e9 \u00e0 l&rsquo;aide de la variance combin\u00e9e d\u00e9form\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"d = \\frac{\\bar{X}_1 - \\bar{X}_2}{s_p}\">$$d = \\frac{\\bar{X}_1 &#8211; \\bar{X}_2}{s_p}$$<\/span><\/p>\n<p>se trouve substantiellement compress\u00e9 vers z\u00e9ro lorsque la variance combin\u00e9e est surestim\u00e9e par rapport \u00e0 la dispersion propre de la population de r\u00e9f\u00e9rence. Il en r\u00e9sulte que les revues de litt\u00e9rature m\u00e9ta-analytiques int\u00e9grant ces estimations biais\u00e9es sous-\u00e9valuent syst\u00e9matiquement la magnitude r\u00e9elle des ph\u00e9nom\u00e8nes psychologiques \u00e9tudi\u00e9s, et que les calculs de puissance a priori r\u00e9alis\u00e9s pour dimensionner les \u00e9tudes ult\u00e9rieures sont fauss\u00e9s d\u00e8s leur base.<\/p>\n<h3>7.3 Invalidation des comparaisons post-hoc multiples<\/h3>\n<p>Lorsqu&rsquo;une analyse de variance \u00e0 trois groupes ou plus r\u00e9v\u00e8le un effet omnibus significatif, la d\u00e9marche exp\u00e9rimentale usuelle consiste \u00e0 d\u00e9ployer une batterie de comparaisons par paires (tests post-hoc) afin de localiser pr\u00e9cis\u00e9ment l&rsquo;origine des contrastes significatifs tout en jugulant l&rsquo;inflation globale du risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce de la famille de tests (<em>Family-Wise Error Rate<\/em>, FWER). La grande majorit\u00e9 des proc\u00e9dures post-hoc canoniques, au premier rang desquelles tr\u00f4ne le test de la plus petite diff\u00e9rence significative honn\u00eate de Tukey (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Tukey%27s_range_test\">Tukey HSD<\/a>), s&rsquo;appuient directement sur le carr\u00e9 moyen r\u00e9siduel unifi\u00e9 de l&rsquo;ANOVA :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"q = \\frac{\\bar{Y}_i - \\bar{Y}_j}{\\sqrt{\\frac{MS_{\\text{r\u00e9siduel}}}{2} \\left( \\frac{1}{n_i} + \\frac{1}{n_j} \\right)}}\">$$q = \\frac{\\bar{Y}_i &#8211; \\bar{Y}_j}{\\sqrt{\\frac{MS_{\\text{r\u00e9siduel}}}{2} \\left( \\frac{1}{n_i} + \\frac{1}{n_j} \\right)}}$$<\/span><\/p>\n<p>Cette architecture repose sur le postulat non n\u00e9gociable que chaque paire de sous-groupes partage la m\u00eame dispersion d&rsquo;erreur, valablement r\u00e9sum\u00e9e par le terme global <span class=\"arabpsy-math-inline\" data-latex=\"MS_{\\text{r\u00e9siduel}}\">$MS_{\\text{r\u00e9siduel}}$<\/span>. En cas de rupture de l&rsquo;homog\u00e9n\u00e9it\u00e9 des variances, cette m\u00e9canique s&rsquo;effondre int\u00e9gralement :<\/p>\n<ul>\n<li>Pour les paires de groupes dont les variances r\u00e9elles sont intrins\u00e8quement faibles par rapport aux autres conditions de l&rsquo;\u00e9tude, l&rsquo;utilisation du <span class=\"arabpsy-math-inline\" data-latex=\"MS_{\\text{r\u00e9siduel}}\">$MS_{\\text{r\u00e9siduel}}$<\/span> (tir\u00e9 vers le haut par les groupes tr\u00e8s dispers\u00e9s) produit des tests post-hoc excessivement conservateurs, incapables de d\u00e9tecter des contrastes r\u00e9els manifestes.<\/li>\n<li>Inversement, pour les paires de groupes dont les variances empiriques sont les plus colossales de l&rsquo;exp\u00e9rience, le d\u00e9nominateur de Tukey sous-estime l&rsquo;erreur propre \u00e0 ce contraste sp\u00e9cifique, d\u00e9clenchant une prolif\u00e9ration anarchique de fausses significations post-hoc.<\/li>\n<li>L&rsquo;assurance math\u00e9matique d&rsquo;un contr\u00f4le strict du FWER \u00e0 exactement 5 % sur l&rsquo;ensemble de la famille de tests n&rsquo;est plus du tout garantie, ce qui invalide l&rsquo;assise probabiliste des conclusions de l&rsquo;article.<\/li>\n<\/ul>\n<h2>8. Solutions param\u00e9triques directes : corrections et tests adapt\u00e9s<\/h2>\n<h3>8.1 Le test t de Welch pour deux \u00e9chantillons ind\u00e9pendants<\/h3>\n<p>Face \u00e0 la faillibilit\u00e9 prouv\u00e9e du test <em>t<\/em> de Student classique en pr\u00e9sence de variances h\u00e9t\u00e9rog\u00e8nes, Bernard Lewis Welch a d\u00e9velopp\u00e9 en 1938 et 1947 une alternative d&rsquo;une \u00e9l\u00e9gance statistique remarquable : le test <em>t<\/em> de Welch (parfois baptis\u00e9 test de Welch-Aspin ou test <em>t<\/em> \u00e0 variances in\u00e9gales). Le g\u00e9nie de la mod\u00e9lisation de Welch a \u00e9t\u00e9 d&rsquo;abandonner purement et simplement le principe du regroupement unifi\u00e9 des variances (<span class=\"arabpsy-math-inline\" data-latex=\"s_p^2\">$s_p^2$<\/span>) au profit du calcul direct de l&rsquo;erreur-type sans mise en commun :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"SE_{(\\bar{X}_1 - \\bar{X}_2)} = \\sqrt{\\frac{s_1^2}{n_1} + \\frac{s_2^2}{n_2}}\">$$SE_{(\\bar{X}_1 &#8211; \\bar{X}_2)} = \\sqrt{\\frac{s_1^2}{n_1} + \\frac{s_2^2}{n_2}}$$<\/span><\/p>\n<p>La statistique <span class=\"arabpsy-math-inline\" data-latex=\"t^*\">$t^*$<\/span> de Welch prend ainsi la forme :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"t^* = \\frac{\\bar{X}_1 - \\bar{X}_2}{\\sqrt{\\frac{s_1^2}{n_1} + \\frac{s_2^2}{n_2}}}\">$$t^* = \\frac{\\bar{X}_1 &#8211; \\bar{X}_2}{\\sqrt{\\frac{s_1^2}{n_1} + \\frac{s_2^2}{n_2}}}$$<\/span><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/het2-2.jpg\" alt=\"\" \/><\/figure>\n<p>Puisque les variances d&rsquo;\u00e9chantillonnage ne sont plus mutualis\u00e9es, la distribution exacte de cette statistique ne correspond plus rigoureusement \u00e0 une loi de Student standard bas\u00e9e sur <span class=\"arabpsy-math-inline\" data-latex=\"n_1 + n_2 - 2\">$n_1 + n_2 &#8211; 2$<\/span> degr\u00e9s de libert\u00e9. Pour surmonter cet obstacle th\u00e9orique (connu sous le nom d&rsquo;impasse de Behrens-Fisher), Welch a mobilis\u00e9 l&rsquo;approximation des degr\u00e9s de libert\u00e9 de Satterthwaite pour d\u00e9river une valeur ajust\u00e9e fractionnaire des degr\u00e9s de libert\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"\\nu\">$\\nu$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"df_{\\text{Welch}}\">$df_{\\text{Welch}}$<\/span>) :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\nu = \\frac{\\left(\\frac{s_1^2}{n_1} + \\frac{s_2^2}{n_2}\\right)^2}{\\frac{\\left(\\frac{s_1^2}{n_1}\\right)^2}{n_1 - 1} + \\frac{\\left(\\frac{s_2^2}{n_2}\\right)^2}{n_2 - 1}}\">$$\\nu = \\frac{\\left(\\frac{s_1^2}{n_1} + \\frac{s_2^2}{n_2}\\right)^2}{\\frac{\\left(\\frac{s_1^2}{n_1}\\right)^2}{n_1 &#8211; 1} + \\frac{\\left(\\frac{s_2^2}{n_2}\\right)^2}{n_2 &#8211; 1}}$$<\/span><\/p>\n<p>Cette \u00e9quation ing\u00e9nieuse ajuste dynamiquement les degr\u00e9s de libert\u00e9 en fonction du poids relatif des variances de chaque groupe. Si les variances sont rigoureusement identiques et les effectifs \u00e9gaux, <span class=\"arabpsy-math-inline\" data-latex=\"\\nu\">$\\nu$<\/span> converge pr\u00e9cis\u00e9ment vers <span class=\"arabpsy-math-inline\" data-latex=\"n_1 + n_2 - 2\">$n_1 + n_2 &#8211; 2$<\/span>. Si une variance domine de mani\u00e8re disproportionn\u00e9e un petit effectif, <span class=\"arabpsy-math-inline\" data-latex=\"\\nu\">$\\nu$<\/span> s&rsquo;effondre pour approcher les degr\u00e9s de libert\u00e9 du seul groupe instable, p\u00e9nalisant ad\u00e9quatement la valeur critique du test et maintenant le taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce rigoureusement riv\u00e9 \u00e0 5 %.<\/p>\n<p>Une pl\u00e9thore de travaux m\u00e9thodologiques majeurs (Delacre, Lakens, &amp; Leys, 2017 ; Ruxton, 2006) ont d\u00e9montr\u00e9 de mani\u00e8re irr\u00e9futable que le test <em>t<\/em> de Welch pr\u00e9sente des performances strictement \u00e9quivalentes au test de Student lorsque les variances sont \u00e9gales, tout en assurant une robustesse absolue lorsqu&rsquo;elles divergent. Ces auteurs recommandent instamment d&rsquo;adopter le test <em>t<\/em> de Welch comme <strong>option par d\u00e9faut syst\u00e9matique<\/strong> pour l&rsquo;ensemble des comparaisons de deux moyennes ind\u00e9pendantes, rendant obsol\u00e8te la routine du pr\u00e9-test de Levene.<\/p>\n<h3>8.2 L&rsquo;ANOVA de Welch pour trois groupes et plus<\/h3>\n<p>De m\u00eame que le test <em>t<\/em> de Welch surpasse son homologue classique pour deux groupes, l&rsquo;<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Welch%27s_t-test\">ANOVA de Welch<\/a> (1951) fournit la solution de premier plan pour tester les diff\u00e9rences de moyennes entre trois groupes ou plus sous h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 av\u00e9r\u00e9e. L&rsquo;ANOVA conventionnelle calcule une moyenne non pond\u00e9r\u00e9e des moyennes de groupes pour d\u00e9river le carr\u00e9 moyen intergroupe, et utilise une variance d&rsquo;erreur globale commune. L&rsquo;ANOVA de Welch brise cette sym\u00e9trie en pond\u00e9rant explicitement chaque moyenne de groupe par l&rsquo;inverse de sa variance d&rsquo;\u00e9chantillonnage sp\u00e9cifique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"w_j = \\frac{n_j}{s_j^2}\">$$w_j = \\frac{n_j}{s_j^2}$$<\/span><\/p>\n<p>Cette formulation octroie un poids statistique pr\u00e9pond\u00e9rant aux groupes dont la variabilit\u00e9 est la plus restreinte (ceux qui offrent l&rsquo;information la plus stable et pr\u00e9cise) et temp\u00e8re proportionnellement l&rsquo;influence des groupes dont la dispersion est erratique. La moyenne g\u00e9n\u00e9rale ajust\u00e9e pond\u00e9r\u00e9e s&rsquo;\u00e9crit alors :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\bar{X}^* = \\frac{\\sum_{j=1}^k w_j \\bar{X}_j}{\\sum_{j=1}^k w_j}\">$$\\bar{X}^* = \\frac{\\sum_{j=1}^k w_j \\bar{X}_j}{\\sum_{j=1}^k w_j}$$<\/span><\/p>\n<p>La statistique <span class=\"arabpsy-math-inline\" data-latex=\"F^*\">$F^*$<\/span> de Welch est calcul\u00e9e \u00e0 partir de ces pond\u00e9rations, et est \u00e9valu\u00e9e face \u00e0 une distribution de Fisher dont les degr\u00e9s de libert\u00e9 r\u00e9siduels sont ajust\u00e9s selon une formule d&rsquo;approximation complexe analogue \u00e0 celle de Satterthwaite. En compl\u00e9ment direct de l&rsquo;ANOVA de Welch, la proc\u00e9dure de comparaisons multiples de <strong>Games-Howell<\/strong> doit \u00eatre syst\u00e9matiquement substitu\u00e9e au Tukey HSD. Le test de Games-Howell n&rsquo;utilise aucun carr\u00e9 moyen commun : il applique la formule de Welch et l&rsquo;approximation de Satterthwaite \u00e0 chaque paire sp\u00e9cifique de groupes compar\u00e9s, tout en int\u00e9grant un contr\u00f4le rigoureux du taux d&rsquo;erreur global de la famille de contrastes bas\u00e9 sur la loi de l&rsquo;\u00e9tendue studentis\u00e9e.<\/p>\n<h3>8.3 Le test de Brown-Forsythe pour les comparaisons multiples<\/h3>\n<p>Outre l&rsquo;approche de Welch, une m\u00e9thode param\u00e9trique alternative remarquable pour rem\u00e9dier \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances au sein des analyses multigroupes a \u00e9t\u00e9 propos\u00e9e par Morton B. Brown et Alan B. Forsythe en 1974 : l&rsquo;ANOVA de Brown-Forsythe modifi\u00e9e. Contrairement \u00e0 l&rsquo;ANOVA de Welch qui module la moyenne g\u00e9n\u00e9rale par des pond\u00e9rations inverses de la variance, la statistique <span class=\"arabpsy-math-inline\" data-latex=\"F^*\">$F^*$<\/span> de Brown-Forsythe conserve le calcul conventionnel de la somme des carr\u00e9s intergroupes au num\u00e9rateur, mais modifie profond\u00e9ment le d\u00e9nominateur.<\/p>\n<p>Le d\u00e9nominateur de Brown-Forsythe n&rsquo;est pas le <span class=\"arabpsy-math-inline\" data-latex=\"MS_{\\text{r\u00e9siduel}}\">$MS_{\\text{r\u00e9siduel}}$<\/span> unifi\u00e9, mais une somme pond\u00e9r\u00e9e des variances d&rsquo;\u00e9chantillon individuelles :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{D\u00e9nominateur}_{BF} = \\sum_{j=1}^k \\left( 1 - \\frac{n_j}{N} \\right) s_j^2\">$$\\text{D\u00e9nominateur}_{BF} = \\sum_{j=1}^k \\left( 1 &#8211; \\frac{n_j}{N} \\right) s_j^2$$<\/span><\/p>\n<p>Les degr\u00e9s de libert\u00e9 du d\u00e9nominateur sont ensuite corrig\u00e9s par une d\u00e9rivation d&rsquo;approximation de type Satterthwaite afin d&rsquo;assurer que le ratio suive convenablement une loi <em>F<\/em> sous <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span>. Les \u00e9tudes de simulation comparatives r\u00e9v\u00e8lent que l&rsquo;ANOVA de Welch et l&rsquo;ANOVA de Brown-Forsythe sont toutes deux immens\u00e9ment sup\u00e9rieures \u00e0 l&rsquo;ANOVA classique d\u00e8s que l&rsquo;homog\u00e9n\u00e9it\u00e9 des variances fait d\u00e9faut. Cependant, de subtiles nuances op\u00e9rationnelles les distinguent :<\/p>\n<ul>\n<li>L&rsquo;ANOVA de Welch conserve une puissance statistique sup\u00e9rieure dans la vaste majorit\u00e9 des configurations exp\u00e9rimentales, en particulier lorsqu&rsquo;une ou deux moyennes extr\u00eames co\u00efncident avec des groupes \u00e0 faible variance.<\/li>\n<li>L&rsquo;ANOVA de Brown-Forsythe se r\u00e9v\u00e8le l\u00e9g\u00e8rement plus performante et robuste que celle de Welch dans un cas pr\u00e9cis : lorsque le design exp\u00e9rimental comporte des groupes d&rsquo;effectifs extr\u00eamement in\u00e9gaux combin\u00e9s \u00e0 des distributions mod\u00e9r\u00e9ment asym\u00e9triques ou leptokurtiques, ou lorsqu&rsquo;un groupe d\u00e9viant cumule \u00e0 la fois une taille tr\u00e8s r\u00e9duite et une variance extr\u00eamement faible.<\/li>\n<\/ul>\n<h2>9. Transformations de variables pour la stabilisation de la variance<\/h2>\n<h3>9.1 Transformations math\u00e9matiques classiques et leurs cibles<\/h3>\n<p>Historiquement, avant la d\u00e9mocratisation de la puissance de calcul informatique et la g\u00e9n\u00e9ralisation logicielle des approximations de Welch et des estimateurs robustes, la r\u00e9ponse m\u00e9thodologique standard face \u00e0 des variances h\u00e9t\u00e9rog\u00e8nes consistait \u00e0 appliquer une transformation math\u00e9matique non lin\u00e9aire \u00e0 l&rsquo;ensemble de la variable d\u00e9pendante (<span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>). L&rsquo;objectif math\u00e9matique d&rsquo;une transformation stabilisatrice de variance est de modifier l&rsquo;\u00e9chelle de mesure de telle sorte que la dispersion des donn\u00e9es cesse de corr\u00e9ler avec la grandeur de la moyenne.<\/p>\n<p>Le choix de la transformation fonctionnelle optimale d\u00e9pend intimement de la relation observ\u00e9e entre la moyenne du groupe (<span class=\"arabpsy-math-inline\" data-latex=\"\\mu\">$\\mu$<\/span>) et sa variance (<span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span>) :<\/p>\n<ul>\n<li><strong>La transformation logarithmique (<span class=\"arabpsy-math-inline\" data-latex=\"ln(Y)\">$ln(Y)$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"\\log_{10}(Y + c)\">$\\log_{10}(Y + c)$<\/span>) :<\/strong> Elle constitue la solution canonique lorsque l&rsquo;\u00e9cart-type est proportionnel \u00e0 la moyenne (<span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_j propto \\mu_j\">$\\sigma_j propto \\mu_j$<\/span>), ce qui correspond \u00e0 un coefficient de variation constant \u00e0 travers les niveaux de facteurs. En psychologie cognitive, les donn\u00e9es chronom\u00e9triques (temps de r\u00e9action) pr\u00e9sentent typiquement cette caract\u00e9ristique : des r\u00e9ponses lentes g\u00e9n\u00e8rent une dispersion proportionnellement dilat\u00e9e. Le logarithme comprime drastiquement les valeurs \u00e9lev\u00e9es tout en \u00e9tirant les scores faibles, r\u00e9gularisant simultan\u00e9ment l&rsquo;homog\u00e9n\u00e9it\u00e9 des variances et r\u00e9duisant l&rsquo;asym\u00e9trie positive.<\/li>\n<li><strong>La transformation racine carr\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"\\sqrt{Y}\">$\\sqrt{Y}$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"\\sqrt{Y + 0,5}\">$\\sqrt{Y + 0,5}$<\/span>) :<\/strong> Elle s&rsquo;applique de mani\u00e8re privil\u00e9gi\u00e9e aux donn\u00e9es de comptage d&rsquo;occurrences discr\u00e8tes, lesquelles suivent classiquement une distribution th\u00e9orique de Poisson o\u00f9 la variance est intrins\u00e8quement \u00e9gale \u00e0 la moyenne (<span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_j^2 = \\mu_j\">$\\sigma_j^2 = \\mu_j$<\/span>). C&rsquo;est fr\u00e9quemment le cas du nombre d&rsquo;erreurs commises dans une t\u00e2che attentionnelle ou du nombre d&rsquo;intrusions mn\u00e9siques.<\/li>\n<li><strong>La transformation arc sinus de la racine carr\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"\\arcsin(\\sqrt{p})\">$\\arcsin(\\sqrt{p})$<\/span>) :<\/strong> Con\u00e7ue pour stabiliser la variance des proportions, ratios ou pourcentages issus de processus binomiaux, o\u00f9 la variance d\u00e9pend directement de la probabilit\u00e9 de succ\u00e8s : <span class=\"arabpsy-math-inline\" data-latex=\"\\text{Var}(p) = \\frac{\\pi(1-\\pi)}{n}\">$\\text{Var}(p) = \\frac{\\pi(1-\\pi)}{n}$<\/span>. La dispersion s&rsquo;effondre naturellement vers les extr\u00e9mit\u00e9s (0 et 1) et culmine \u00e0 0,5 ; la fonction angulaire \u00e9tire les queues de distribution pour r\u00e9tablir une variance conditionnelle plane.<\/li>\n<\/ul>\n<h3>9.2 La famille de transformations param\u00e9triques Box-Cox<\/h3>\n<p>Plut\u00f4t que d&rsquo;appliquer de mani\u00e8re empirique et discr\u00e9tionnaire l&rsquo;une des fonctions classiques \u00e9voqu\u00e9es ci-dessus, George Box et David Cox ont d\u00e9velopp\u00e9 en 1964 une m\u00e9thode formalis\u00e9e permettant d&rsquo;identifier algorithmiquement la transformation math\u00e9matique optimale au sein d&rsquo;une famille continue de fonctions puissances. La <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Power_transform#Box%E2%80%93Cox_transformation\">transformation de Box-Cox<\/a> est param\u00e9tr\u00e9e par un exposant scalaire continu, g\u00e9n\u00e9ralement not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"lambda\">$lambda$<\/span> (lambda) :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Y^{(\\lambda)} = \\begin{\\cases} \\frac{Y^\\lambda - 1}{\\lambda} &amp;a\\mp; \\text{si } \\lambda \\neq 0  \\ln(Y) &amp;a\\mp; \\text{si } \\lambda = 0 \\end{\\cases}\">$$Y^{(\\lambda)} = \\begin{\\cases} \\frac{Y^\\lambda &#8211; 1}{\\lambda} &#038;a\\mp; \\text{si } \\lambda \\neq 0  \\ln(Y) &#038;a\\mp; \\text{si } \\lambda = 0 \\end{\\cases}$$<\/span><\/p>\n<p>La valeur num\u00e9rique du param\u00e8tre <span class=\"arabpsy-math-inline\" data-latex=\"lambda\">$lambda$<\/span> est estim\u00e9e par la m\u00e9thode du maximum de vraisemblance (ou par maximisation du profil de log-vraisemblance) sur les donn\u00e9es observ\u00e9es. Cette famille englobe un continuum de transformations usuelles : <span class=\"arabpsy-math-inline\" data-latex=\"lambda = 1\">$lambda = 1$<\/span> correspond \u00e0 l&rsquo;absence de transformation (identit\u00e9 lin\u00e9aire) ; <span class=\"arabpsy-math-inline\" data-latex=\"lambda = 0,5\">$lambda = 0,5$<\/span> \u00e9quivaut \u00e0 la racine carr\u00e9e ; <span class=\"arabpsy-math-inline\" data-latex=\"lambda = 0\">$lambda = 0$<\/span> correspond au logarithme n\u00e9p\u00e9rien ; <span class=\"arabpsy-math-inline\" data-latex=\"lambda = -1\">$lambda = -1$<\/span> repr\u00e9sente la transformation inverse (<span class=\"arabpsy-math-inline\" data-latex=\"1\/Y\">$1\/Y$<\/span>).<\/p>\n<p>Une contrainte computationnelle absolue de la formulation classique de Box-Cox est qu&rsquo;elle exige imp\u00e9rativement des valeurs brutes strictement positives (<span class=\"arabpsy-math-inline\" data-latex=\"Y_i &gt; 0\">$Y_i &gt; 0$<\/span>). Pour \u00e9tendre cette approche \u00e0 des variables comportementales pouvant int\u00e9grer des z\u00e9ros stricts ou des scores n\u00e9gatifs (comme des indices d&rsquo;amor\u00e7age ou des variations de temp\u00e9rature corporelle), In-Kwon Yeo et Richard A. Johnson ont d\u00e9velopp\u00e9 en 2000 une extension g\u00e9n\u00e9ralis\u00e9e (la transformation de Yeo-Johnson), qui pr\u00e9serve les propri\u00e9t\u00e9s remarquables de stabilisation de variance de Box-Cox tout en autorisant le traitement sans rupture des valeurs n\u00e9gatives ou nulles.<\/p>\n<h3>9.3 D\u00e9fis d&rsquo;interpr\u00e9tation clinique et psychologique post-transformation<\/h3>\n<p>Bien que les transformations math\u00e9matiques r\u00e9solvent avec efficacit\u00e9 les violations conjointes de normalit\u00e9 et d&rsquo;homosc\u00e9dasticit\u00e9, leur impl\u00e9mentation en sciences humaines et m\u00e9dicales soul\u00e8ve des difficult\u00e9s interpr\u00e9tatives et communicationnelles majeures que les analystes ne sauraient ignorer. Transformer une variable modifie fondamentalement l&rsquo;unit\u00e9 de mesure et la nature substantielle de l&rsquo;\u00e9chelle psychologique sous-jacente.<\/p>\n<p>Par exemple, si un chercheur applique une transformation logarithmique sur des scores d&rsquo;anxi\u00e9t\u00e9 mesur\u00e9s par le questionnaire STAI, l&rsquo;analyse inf\u00e9rentielle teste d\u00e9sormais l&rsquo;\u00e9galit\u00e9 des moyennes des logarithmes des scores (<span class=\"arabpsy-math-inline\" data-latex=\"\\overline{\\ln(Y)}\">$\\overline{\\ln(Y)}$<\/span>), et non plus les moyennes des scores originaux. Les conclusions tir\u00e9es portent d\u00e8s lors sur une m\u00e9trique compress\u00e9e difficilement transposable au diagnostic clinique usuel. Pour restituer une intelligibilit\u00e9 descriptive, le chercheur proc\u00e8de souvent \u00e0 une transformation inverse (dite <em>back-transformation<\/em>) par exponentiation des moyennes transform\u00e9es. Or, en vertu de l&rsquo;in\u00e9galit\u00e9 de Jensen, l&rsquo;exponentielle de la moyenne d&rsquo;une s\u00e9rie de logarithmes ne restitue pas la moyenne arithm\u00e9tique originale, mais correspond rigoureusement \u00e0 la <strong>moyenne g\u00e9om\u00e9trique<\/strong> des donn\u00e9es brutes, qui est syst\u00e9matiquement inf\u00e9rieure \u00e0 la moyenne arithm\u00e9tique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\exp(E[\\ln(Y)]) = \\text{Moyenne G\u00e9om\u00e9trique} le E[Y]\">$$\\exp(E[\\ln(Y)]) = \\text{Moyenne G\u00e9om\u00e9trique} le E[Y]$$<\/span><\/p>\n<p>De surcro\u00eet, les contrastes, pentes de r\u00e9gression et intervalles de confiance transform\u00e9s de nouveau dans l&rsquo;\u00e9chelle d&rsquo;origine cessent d&rsquo;\u00eatre sym\u00e9triques, complexifiant la pr\u00e9sentation formelle des r\u00e9sultats. Dans le cadre de plans factoriels, appliquer une transformation non lin\u00e9aire alt\u00e8re profond\u00e9ment la nature math\u00e9matique des interactions : une interaction multiplicative dans l&rsquo;\u00e9chelle brute peut \u00eatre int\u00e9gralement r\u00e9sorb\u00e9e dans l&rsquo;\u00e9chelle logarithmique, ou inversement, une interaction totalement factice peut \u00e9merger par distorsion d&rsquo;\u00e9chelle. C&rsquo;est pourquoi, face \u00e0 ces \u00e9cueils interpr\u00e9tatifs, la communaut\u00e9 statistique contemporaine privil\u00e9gie d\u00e9sormais massivement les approches directes sans transformation d&rsquo;\u00e9chelle, fond\u00e9es sur les tests de Welch ou les mod\u00e8les lin\u00e9aires robustes.<\/p>\n<h2>10. Approches avanc\u00e9es en r\u00e9gression : h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 et mod\u00e9lisation robuste<\/h2>\n<h3>10.1 Erreurs-types robustes \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 (estimateurs sandwich)<\/h3>\n<p>Dans la mod\u00e9lisation par r\u00e9gression lin\u00e9aire multiple OLS, la strat\u00e9gie contemporaine dominante pour neutraliser d\u00e9finitivement l&rsquo;impact de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 sans avoir recours \u00e0 des transformations d&rsquo;\u00e9chelles d\u00e9stabilisantes repose sur l&rsquo;impl\u00e9mentation d&rsquo;erreurs-types robustes \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 (<em>Heteroscedasticity-Consistent Standard Errors<\/em>, ou estimateurs HC). Cette m\u00e9thodologie a \u00e9t\u00e9 initialement th\u00e9oris\u00e9e par Halbert White en 1980, qui a d\u00e9riv\u00e9 l&rsquo;estimateur canonique baptis\u00e9 <strong>HC0<\/strong>.<\/p>\n<p>L&rsquo;estimateur sandwich de White conserve int\u00e9gralement les coefficients de r\u00e9gression originaux <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\mathbf{\\beta}}_{\\text{OLS}}\">$\\hat{\\mathbf{\\beta}}_{\\text{OLS}}$<\/span>, garantissant ainsi l&rsquo;intelligibilit\u00e9 directe des pentes. Il substitue en revanche \u00e0 la matrice de variance conventionnelle l&rsquo;estimation empirique suivante :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Var}_{\\text{HC0}}(\\hat{\\mathbf{\\beta}}) = (\\mathbf{X}&#039;\\mathbf{X})^{-1} \\left( \\sum_{i=1}^N e_i^2 \\mathbf{x}_i \\mathbf{x}_i&#039; \\right) (\\mathbf{X}&#039;\\mathbf{X})^{-1}\">$$\\text{Var}_{\\text{HC0}}(\\hat{\\mathbf{\\beta}}) = (\\mathbf{X}&rsquo;\\mathbf{X})^{-1} \\left( \\sum_{i=1}^N e_i^2 \\mathbf{x}_i \\mathbf{x}_i&rsquo; \\right) (\\mathbf{X}&rsquo;\\mathbf{X})^{-1}$$<\/span><\/p>\n<p>Dans cette structure, les carr\u00e9s des r\u00e9sidus individuels observ\u00e9s (<span class=\"arabpsy-math-inline\" data-latex=\"e_i^2\">$e_i^2$<\/span>) servent d&rsquo;estimateurs directs et non contraints des variances sp\u00e9cifiques conditionnelles <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_i^2\">$\\sigma_i^2$<\/span>. Bien que l&rsquo;estimateur HC0 soit asymptotiquement parfait (convergeant rigoureusement vers la vraie variance lorsque <span class=\"arabpsy-math-inline\" data-latex=\"N to \\infty\">$N to \\infty$<\/span>), des recherches ult\u00e9rieures ont d\u00e9montr\u00e9 qu&rsquo;il manifeste un biais s\u00e9v\u00e8re sur des \u00e9chantillons de taille modeste (<span class=\"arabpsy-math-inline\" data-latex=\"N &lt; 250\">$N &lt; 250$<\/span>), sous-estimant les erreurs-types et maintenant une inflation d&#039;erreur de Type I.<\/p>\n<p>Pour r\u00e9soudre cette d\u00e9faillance sur petits \u00e9chantillons, une s\u00e9rie de raffinements successifs a \u00e9t\u00e9 mod\u00e9lis\u00e9e :<\/p>\n<ul>\n<li><strong>HC1 (MacKinnon &amp; White, 1985) :<\/strong> Introduit une correction par les degr\u00e9s de libert\u00e9 en multipliant l&rsquo;ensemble par le facteur scalaire <span class=\"arabpsy-math-inline\" data-latex=\"\\frac{N}{N - p - 1}\">$\\frac{N}{N &#8211; p &#8211; 1}$<\/span>, temp\u00e9rant l\u00e9g\u00e8rement le biais sur petits effectifs.<\/li>\n<li><strong>HC2 (MacKinnon &amp; White, 1985) :<\/strong> Divise le carr\u00e9 de chaque r\u00e9sidu par la quantit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"(1 - h_{ii})\">$(1 &#8211; h_{ii})$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> repr\u00e9sente le levier statistique (<em>leverage<\/em>) de la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me observation issu de la matrice de projection chapeau (<span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{H}\">$\\mathbf{H}$<\/span>). Les points dot\u00e9s d&rsquo;un levier \u00e9lev\u00e9 voient ainsi leur influence r\u00e9siduelle correctement r\u00e9\u00e9valu\u00e9e.<\/li>\n<li><strong>HC3 (Davidson &amp; MacKinnon, 1993) :<\/strong> P\u00e9nalise plus vigoureusement les points \u00e0 fort levier en divisant chaque r\u00e9sidu au carr\u00e9 par <span class=\"arabpsy-math-inline\" data-latex=\"(1 - h_{ii})^2\">$(1 &#8211; h_{ii})^2$<\/span>. Les vastes simulations conduites par Long et Ervin (2000) ont consacr\u00e9 l&rsquo;estimateur <strong>HC3<\/strong> comme le standard d&rsquo;excellence pour la recherche quantitative : il offre un contr\u00f4le irr\u00e9prochable du taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce, m\u00eame avec de faibles effectifs (<span class=\"arabpsy-math-inline\" data-latex=\"N &lt; 50\">$N &lt; 50$<\/span>), ind\u00e9pendamment de la structure ou de la gravit\u00e9 de l&#039;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 pr\u00e9sente.<\/li>\n<li><strong>HC4 (Cribari-Neto, 2004) :<\/strong> Con\u00e7u sp\u00e9cifiquement pour juguler les effets d\u00e9l\u00e9t\u00e8res d&rsquo;observations combinant \u00e0 la fois un fort levier et une forte variance r\u00e9siduelle, en exploitant un exposant modul\u00e9 par le ratio de levier.<\/li>\n<\/ul>\n<h3>10.2 La r\u00e9gression par les moindres carr\u00e9s pond\u00e9r\u00e9s (WLS)<\/h3>\n<p>Si les erreurs-types robustes HC corrigent parfaitement l&rsquo;inf\u00e9rence (tests <em>t<\/em> et intervalles de confiance) sans alt\u00e9rer les coefficients OLS, elles ne restaurent pas l&rsquo;efficacit\u00e9 statistique th\u00e9orique : les estimateurs OLS demeurent sous-optimaux sur le plan de la variance d&rsquo;\u00e9chantillonnage minimale. Lorsque la forme pr\u00e9cise de la variance conditionnelle des erreurs est connue ou peut \u00eatre estim\u00e9e de mani\u00e8re cr\u00e9dible, la m\u00e9thode des <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Weighted_least_squares\">moindres carr\u00e9s pond\u00e9r\u00e9s<\/a> (<em>Weighted Least Squares<\/em>, WLS) permet de recouvrer pleinement le statut de meilleur estimateur lin\u00e9aire sans biais (BLUE).<\/p>\n<p>Le principe fondamental de la r\u00e9gression WLS consiste \u00e0 assigner \u00e0 chaque observation <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> un poids statistique (<span class=\"arabpsy-math-inline\" data-latex=\"w_i\">$w_i$<\/span>) inversement proportionnel \u00e0 sa variance d&rsquo;erreur attendue :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"w_i = \\frac{1}{\\sigma_i^2}\">$$w_i = \\frac{1}{\\sigma_i^2}$$<\/span><\/p>\n<p>La minimisation math\u00e9matique ne porte plus sur la somme ordinaire des carr\u00e9s des r\u00e9sidus, mais sur la somme pond\u00e9r\u00e9e de ces \u00e9carts :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Minimiser } \\sum_{i=1}^N w_i (Y_i - \\hat{Y}_i)^2\">$$\\text{Minimiser } \\sum_{i=1}^N w_i (Y_i &#8211; \\hat{Y}_i)^2$$<\/span><\/p>\n<p>Sur le plan g\u00e9om\u00e9trique et conceptuel, la m\u00e9thode WLS r\u00e9\u00e9chelonne l&rsquo;espace des observations : les points appartenant \u00e0 des segments o\u00f9 la dispersion r\u00e9siduelle est colossale voient leur poids d&rsquo;entra\u00eenement fortement att\u00e9nu\u00e9, \u00e9vitant qu&rsquo;ils ne tirent erratiquement la ligne de r\u00e9gression vers eux. \u00c0 l&rsquo;inverse, les observations collect\u00e9es avec une pr\u00e9cision microscopique (faible dispersion) re\u00e7oivent une influence pr\u00e9pond\u00e9rante dans la d\u00e9termination des pentes. L&rsquo;estimateur analytique WLS s&rsquo;exprime sous forme matricielle directe :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{\\mathbf{\\beta}}_{\\text{WLS}} = (\\mathbf{X}&#039;\\mathbf{W}\\mathbf{X})^{-1} \\mathbf{X}&#039;\\mathbf{W}\\mathbf{Y}\">$$\\hat{\\mathbf{\\beta}}_{\\text{WLS}} = (\\mathbf{X}&rsquo;\\mathbf{W}\\mathbf{X})^{-1} \\mathbf{X}&rsquo;\\mathbf{W}\\mathbf{Y}$$<\/span><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{W}\">$\\mathbf{W}$<\/span> est une matrice diagonale contenant les poids <span class=\"arabpsy-math-inline\" data-latex=\"w_i\">$w_i$<\/span>. Le d\u00e9fi majeur de l&rsquo;impl\u00e9mentation WLS en pratique r\u00e9side dans le fait que les v\u00e9ritables variances <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_i^2\">$\\sigma_i^2$<\/span> sont g\u00e9n\u00e9ralement inconnues de l&rsquo;investigateur. Elles doivent \u00eatre estim\u00e9es dans un premier temps en r\u00e9gressant le logarithme des r\u00e9sidus OLS au carr\u00e9 sur les variables explicatives suspect\u00e9es de piloter l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9, d\u00e9marche connue sous le nom de moindres carr\u00e9s g\u00e9n\u00e9ralis\u00e9s r\u00e9alisables (<em>Feasible Generalized Least Squares<\/em>, FGLS). Si cette mod\u00e9lisation de la variance d&rsquo;erreur est mal sp\u00e9cifi\u00e9e, l&rsquo;estimateur WLS peut paradoxalement s&rsquo;av\u00e9rer moins pr\u00e9cis que le mod\u00e8le OLS standard.<\/p>\n<h3>10.3 Mod\u00e8les lin\u00e9aires mixtes et mod\u00e8les de localisation-\u00e9chelle<\/h3>\n<p>L&rsquo;essor de la mod\u00e9lisation statistique avanc\u00e9e a permis de transcender l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 en cessant de la consid\u00e9rer comme une simple nuisance \u00e0 corriger, pour en faire une composante th\u00e9orique noble du mod\u00e8le. Les mod\u00e8les lin\u00e9aires mixtes (<em>Linear Mixed Models<\/em>, LMM), amplement exploit\u00e9s dans l&rsquo;analyse de donn\u00e9es hi\u00e9rarchiques ou longitudinales, permettent de param\u00e9trer explicitement et directement la structure de covariance des r\u00e9sidus intra-individuels ou inter-groupes.<\/p>\n<p>Au sein des progiciels statistiques contemporains (comme le package <code>nlme<\/code> en langage R), des fonctions de variance explicites (telles que <code>varIdent<\/code>) permettent de d\u00e9clarer que chaque sous-groupe d&rsquo;un plan exp\u00e9rimental poss\u00e8de sa propre variance d&rsquo;erreur r\u00e9siduelle stochastique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Var}(\\mathbf{\\epsilon}_j) = \\sigma^2 \\mathbf{\\Lambda}_j\">$$\\text{Var}(\\mathbf{\\epsilon}_j) = \\sigma^2 \\mathbf{\\Lambda}_j$$<\/span><\/p>\n<p>Cette sp\u00e9cification permet d&rsquo;estimer simultan\u00e9ment les param\u00e8tres de moyennes \u00e0 effets fixes et les ratios d&rsquo;\u00e9cart-type propres \u00e0 chaque modalit\u00e9 sans forcer un pooling artificiel de la variance. L&rsquo;inf\u00e9rence sur les effets fixes est alors calcul\u00e9e en tenant compte de ces variances h\u00e9t\u00e9rog\u00e8nes estim\u00e9es par le maximum de vraisemblance restreint (REML).<\/p>\n<p>L&rsquo;aboutissement contemporain de cette d\u00e9marche conceptuelle s&rsquo;incarne dans les mod\u00e8les de localisation-\u00e9chelle (<em>Location-Scale Models<\/em> ou <em>Double Generalized Linear Models<\/em>, DGLM), popularis\u00e9s notamment par Donald Hedeker en psychologie quantitative. Ces architectures sophistiqu\u00e9es estiment deux \u00e9quations de r\u00e9gression simultan\u00e9es et interconnect\u00e9es \u00e0 partir d&rsquo;un m\u00eame jeu de donn\u00e9es :<\/p>\n<ul>\n<li><strong>L&rsquo;\u00e9quation de localisation (Location) :<\/strong> Elle mod\u00e9lise la moyenne conditionnelle de la variable d\u00e9pendante (<span class=\"arabpsy-math-inline\" data-latex=\"\\mu_i = E[Y_i mid X_i]\">$\\mu_i = E[Y_i mid X_i]$<\/span>) \u00e0 travers une combinaison lin\u00e9aire standard de r\u00e9gresseurs.<\/li>\n<li><strong>L&rsquo;\u00e9quation d&rsquo;\u00e9chelle (Scale) :<\/strong> Elle mod\u00e9lise explicitement le logarithme de la variance r\u00e9siduelle (<span class=\"arabpsy-math-inline\" data-latex=\"\\ln(\\sigma_i^2)\">$\\ln(\\sigma_i^2)$<\/span>) en fonction d&rsquo;un ensemble de covariables sp\u00e9cifiques.<\/li>\n<\/ul>\n<p>Cette m\u00e9thodologie ouvre des perspectives scientifiques r\u00e9volutionnaires en psychologie diff\u00e9rentielle : elle autorise l&rsquo;exp\u00e9rimentateur \u00e0 tester formellement si une intervention th\u00e9rapeutique modifie l&rsquo;instabilit\u00e9 ou la variabilit\u00e9 intra-individuelle des patients (l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9), ind\u00e9pendamment de son effet sur leur niveau moyen de performance cognitive ou affective.<\/p>\n<h2>11. Alternatives non param\u00e9triques et approches par r\u00e9\u00e9chantillonnage<\/h2>\n<h3>11.1 Tests non param\u00e9triques bas\u00e9s sur les rangs et leurs limites<\/h3>\n<p>Lorsqu&rsquo;un chercheur observe que ses donn\u00e9es violent simultan\u00e9ment les postulats de normalit\u00e9 univari\u00e9e et d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances, un r\u00e9flexe m\u00e9thodologique extr\u00eamement r\u00e9pandu consiste \u00e0 abandonner les tests param\u00e9triques pour se r\u00e9fugier dans l&rsquo;arsenal des tests non param\u00e9triques bas\u00e9s sur les rangs. Le <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Mann%E2%80%93Whitney_U_test\">test U de Mann-Whitney-Wilcoxon<\/a> (pour deux \u00e9chantillons) et le test de Kruskal-Wallis (pour trois groupes ou plus) sont ainsi couramment mobilis\u00e9s sous l&rsquo;assertion rassurante qu&rsquo;ils sont exempts de toute contrainte distributionnelle.<\/p>\n<p>Or, cette croyance rel\u00e8ve d&rsquo;une <strong>m\u00e9conception th\u00e9orique majeure<\/strong> d\u00e9nonc\u00e9e avec insistance par les math\u00e9maticiens statisticiens. Si les tests de Mann-Whitney et de Kruskal-Wallis ne requi\u00e8rent effectivement aucune distribution gaussienne sous-jacente, leur utilisation en tant que tests stricts d&rsquo;\u00e9galit\u00e9 des m\u00e9dianes ou des param\u00e8tres de position exige express\u00e9ment l&rsquo;hypoth\u00e8se de formes distributionnelles identiques entre les populations compar\u00e9es, ce qui inclut express\u00e9ment <strong>l&rsquo;\u00e9galit\u00e9 des variances<\/strong>.<\/p>\n<p>Le test de Wilcoxon \u00e9value en r\u00e9alit\u00e9 l&rsquo;hypoth\u00e8se nulle de dominance stochastique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"H_0: P(X_1 &gt; X_2) = 0,5\">$$H_0: P(X_1 &gt; X_2) = 0,5$$<\/span><\/p>\n<p>Si deux populations partagent une m\u00eame m\u00e9diane th\u00e9orique mais poss\u00e8dent des variances radicalement in\u00e9gales, la distribution conjointe des rangs s&rsquo;en trouve profond\u00e9ment d\u00e9s\u00e9quilibr\u00e9e. Les travaux de simulation ont d\u00e9montr\u00e9 que le test de Mann-Whitney-Wilcoxon appliqu\u00e9 \u00e0 des groupes d&rsquo;effectifs in\u00e9gaux et de variances h\u00e9t\u00e9rog\u00e8nes subit une inflation de l&rsquo;erreur de Type I tout aussi cataclysmique que le test <em>t<\/em> de Student classique (le taux <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha\">$\\alpha$<\/span> d\u00e9passant r\u00e9guli\u00e8rement les 25 % en cas d&rsquo;appariement n\u00e9gatif). Par cons\u00e9quent, substituer un test sur les rangs classique au test <em>t<\/em> pour r\u00e9soudre un probl\u00e8me d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 constitue une erreur de jugement statistique.<\/p>\n<h3>11.2 M\u00e9thodes de r\u00e9\u00e9chantillonnage (Bootstrapping)<\/h3>\n<p>Face aux d\u00e9faillances des m\u00e9thodes analytiques fond\u00e9es sur des hypoth\u00e8ses rigides d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances, les techniques de calcul intensif par r\u00e9\u00e9chantillonnage informatique, au premier plan desquelles s&rsquo;impose le <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Bootstrapping_(statistics)\">bootstrapping<\/a> non param\u00e9trique th\u00e9oris\u00e9 par Bradley Efron (1979), fournissent une voie de r\u00e9solution empirique particuli\u00e8rement puissante.<\/p>\n<p>Le bootstrap s&rsquo;affranchit de la n\u00e9cessit\u00e9 de d\u00e9river une loi de distribution th\u00e9orique ferm\u00e9e sous des postulats d&rsquo;homog\u00e9n\u00e9it\u00e9 stricte. Son principe op\u00e9ratoire consiste \u00e0 traiter l&rsquo;\u00e9chantillon empirique initial comme une pseudo-population repr\u00e9sentative, au sein de laquelle l&rsquo;ordinateur tire par tirage avec remise des milliers de r\u00e9-\u00e9chantillons (typiquement <span class=\"arabpsy-math-inline\" data-latex=\"B = 2000\">$B = 2000$<\/span> \u00e0 $10,000$ it\u00e9rations) de taille <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span>. Pour chaque r\u00e9plication bootstrap, la statistique d&rsquo;int\u00e9r\u00eat \u2014 par exemple la diff\u00e9rence brute de moyennes (<span class=\"arabpsy-math-inline\" data-latex=\"\\bar{X}_1^* - \\bar{X}_2^*\">$\\bar{X}_1^* &#8211; \\bar{X}_2^*$<\/span>) ou une pente de r\u00e9gression <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}_j^*\">$\\hat{\\beta}_j^*$<\/span> \u2014 est calcul\u00e9e et consign\u00e9e.<\/p>\n<p>La compilation ordonn\u00e9e de ces milliers d&rsquo;estimations empiriques trace la distribution d&rsquo;\u00e9chantillonnage directe du param\u00e8tre. L&rsquo;inf\u00e9rence s&rsquo;effectue ensuite par la construction d&rsquo;intervalles de confiance empiriques :<\/p>\n<ul>\n<li><strong>L&rsquo;intervalle de confiance percentile bootstrap :<\/strong> Les bornes inf\u00e9rieure et sup\u00e9rieure \u00e0 95 % sont simplement d\u00e9finies par le 2,5e et le 97,5e percentile de la distribution empirique g\u00e9n\u00e9r\u00e9e. Cet intervalle capture la variance intrins\u00e8que sp\u00e9cifique de chaque \u00e9chantillon sans rien postuler sur leur \u00e9quivalence pr\u00e9alable.<\/li>\n<li><strong>L&rsquo;intervalle de confiance BCa (Bias-Corrected and Accelerated) :<\/strong> Il affine la m\u00e9thode des percentiles en introduisant une correction math\u00e9matique formalis\u00e9e pour l&rsquo;asym\u00e9trie empirique (param\u00e8tre d&rsquo;acc\u00e9l\u00e9ration <span class=\"arabpsy-math-inline\" data-latex=\"a\">$a$<\/span>) et le biais m\u00e9dian (param\u00e8tre <span class=\"arabpsy-math-inline\" data-latex=\"z_0\">$z_0$<\/span>). L&rsquo;intervalle BCa d\u00e9montre une pr\u00e9cision de couverture remarquable face \u00e0 des structures de variances h\u00e9t\u00e9rog\u00e8nes conjugu\u00e9es \u00e0 des distributions fortement asym\u00e9triques.<\/li>\n<\/ul>\n<p>En psychologie moderne, le bootstrap constitue l&rsquo;\u00e9talon-or incontest\u00e9 pour tester les effets indirects dans les mod\u00e8les de m\u00e9diation (Preacher &amp; Hayes, 2008), o\u00f9 les distributions des produits de coefficients violent de mani\u00e8re structurelle les conditions canoniques d&rsquo;homog\u00e9n\u00e9it\u00e9 et de normalit\u00e9.<\/p>\n<h3>11.3 Tests de permutation exacte et asymptotique<\/h3>\n<p>Les tests de randomisation ou de permutation, initialement conceptualis\u00e9s par Ronald Fisher dans les ann\u00e9es 1930, offrent une alternative non param\u00e9trique rigoureuse pour l&rsquo;\u00e9valuation de l&rsquo;hypoth\u00e8se nulle sans d\u00e9pendance aux lois de probabilit\u00e9 classiques. Leur fondement logique repose sur le principe d&rsquo;interchangeabilit\u00e9 (<em>exchangeability<\/em>) : sous l&rsquo;hypoth\u00e8se nulle d&rsquo;absence totale d&rsquo;effet ou de diff\u00e9rence entre les traitements, les \u00e9tiquettes assignant chaque observation \u00e0 un groupe exp\u00e9rimental donn\u00e9 sont purement arbitraires.<\/p>\n<p>L&rsquo;algorithme de permutation r\u00e9assigne al\u00e9atoirement les scores observ\u00e9s entre les diff\u00e9rents groupes exp\u00e9rimentaux en pr\u00e9servant les tailles d&rsquo;\u00e9chantillons originelles, puis recalcule la statistique de divergence choisie \u00e0 chaque it\u00e9ration. En \u00e9valuant la proportion de permutations r\u00e9assign\u00e9es qui affichent une statistique au moins aussi extr\u00eame que celle mesur\u00e9e sur les donn\u00e9es r\u00e9elles, l&rsquo;exp\u00e9rimentateur d\u00e9rive une <em>p<\/em>-valeur exacte ou asymptotique exempte de toute approximation param\u00e9trique.<\/p>\n<p>Cependant, une mise en garde \u00e9pist\u00e9mologique cruciale s&rsquo;impose : le principe fondamental d&rsquo;interchangeabilit\u00e9 sous l&rsquo;hypoth\u00e8se nulle n&rsquo;est rigoureusement satisfait que si les distributions des populations d&rsquo;origine sont parfaitement identiques en tout point (m\u00eame forme, m\u00eame m\u00e9diane, et <strong>m\u00eame variance<\/strong>). Si les groupes manifestent des variances intrins\u00e8quement distinctes, les observations ne sont plus th\u00e9oriquement interchangeables sous <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span>. Un test de permutation standard utilisant la diff\u00e9rence brute de moyennes <span class=\"arabpsy-math-inline\" data-latex=\"(\\bar{X}_1 - \\bar{X}_2)\">$(\\bar{X}_1 &#8211; \\bar{X}_2)$<\/span> comme statistique de pivot verra son taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce s\u00e9v\u00e8rement biais\u00e9 en cas d&rsquo;effectifs in\u00e9gaux et de variances h\u00e9t\u00e9rog\u00e8nes.<\/p>\n<p>Pour immuniser les tests de permutation face \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances, les statisticiens contemporains pr\u00e9conisent d&rsquo;utiliser non pas une statistique brute, mais une statistique <strong>\u00e9tudiante de Welch<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"t^*\">$t^*$<\/span>) comme m\u00e9trique de permutation (proc\u00e9dure de permutation dite studentis\u00e9e ou test de Welch perm\u00e9). Cette hybridation algorithmique combine l&rsquo;immunit\u00e9 structurelle de la correction de Welch face aux variances in\u00e9gales et la flexibilit\u00e9 distributionnelle totale du r\u00e9\u00e9chantillonnage exact.<\/p>\n<h2>12. Protocole d\u00e9cisionnel et bonnes pratiques de rapportage en recherche psychologique<\/h2>\n<h3>12.1 Arbre d\u00e9cisionnel pour l&rsquo;analyse des donn\u00e9es quantitatives<\/h3>\n<p>Afin de traduire l&rsquo;ensemble de ces imp\u00e9ratifs th\u00e9oriques en une pratique de laboratoire op\u00e9rationnelle, reproductible et robuste, la formalisation d&rsquo;un arbre d\u00e9cisionnel standardis\u00e9 s&rsquo;av\u00e8re indispensable pour tout chercheur en psychologie quantitative. Ce protocole syst\u00e9matique s&rsquo;articule en trois \u00e9tapes s\u00e9quentielles int\u00e9gr\u00e9es :<\/p>\n<ol>\n<li><strong>\u00c9tape 1 : Exploration visuelle et descriptive pr\u00e9alable :<\/strong>\n<ul>\n<li>G\u00e9n\u00e9rer syst\u00e9matiquement des diagrammes en bo\u00eete enrichis (boxplots superpos\u00e9s aux points r\u00e9els individuels) pour chaque modalit\u00e9 et condition. Examiner l&rsquo;homog\u00e9n\u00e9it\u00e9 de l&rsquo;\u00e9cart interquartile (hauteur des bo\u00eetes) et la sym\u00e9trie des moustaches.<\/li>\n<li>En r\u00e9gression, inspecter rigoureusement le graphique des r\u00e9sidus versus valeurs pr\u00e9dites (recherche active d&rsquo;un motif en entonnoir) ainsi que le graphique <em>Scale-Location<\/em> (horizontalit\u00e9 de la courbe LOESS).<\/li>\n<li>Calculer et rapporter le ratio des variances maximales et minimales (<span class=\"arabpsy-math-inline\" data-latex=\"s_{\\text{\\max}}^2 \/ s_{\\text{\\min}}^2\">$s_{\\text{\\max}}^2 \/ s_{\\text{\\min}}^2$<\/span>). Si ce ratio exc\u00e8de 1,5 avec des effectifs in\u00e9gaux, ou s&rsquo;il exc\u00e8de 4 avec des effectifs strictement \u00e9gaux, une alerte m\u00e9thodologique majeure doit \u00eatre consign\u00e9e.<\/li>\n<\/ul>\n<\/li>\n<li><strong>\u00c9tape 2 : \u00c9valuation conjointe de la normalit\u00e9 et de l&rsquo;\u00e9quilibre des effectifs :<\/strong>\n<ul>\n<li>V\u00e9rifier le ratio des tailles d&rsquo;\u00e9chantillons (<span class=\"arabpsy-math-inline\" data-latex=\"n_{\\text{grand}} \/ n_{\\text{petit}}\">$n_{\\text{grand}} \/ n_{\\text{petit}}$<\/span>). D\u00e8s lors que ce ratio s&rsquo;\u00e9carte de 1:1 (plans d\u00e9s\u00e9quilibr\u00e9s), la tol\u00e9rance de l&rsquo;analyse classique face \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 devient strictement nulle.<\/li>\n<li>En cas de doute sur la s\u00e9v\u00e9rit\u00e9 de la violation, calculer le test de Brown-Forsythe (Levene bas\u00e9 sur la m\u00e9diane) ou de Fligner-Killeen en guise d&rsquo;information diagnostique contextuelle, sans jamais faire d\u00e9pendre l&rsquo;analyse principale exclusivement de son franchissement binaire du seuil 0,05.<\/li>\n<\/ul>\n<\/li>\n<li><strong>\u00c9tape 3 : S\u00e9lection raisonn\u00e9e et automatis\u00e9e du mod\u00e8le inf\u00e9rentiel :<\/strong>\n<ul>\n<li><em>Comparaison de deux groupes :<\/em> Appliquer de mani\u00e8re syst\u00e9matique et inconditionnelle le <strong>test <em>t<\/em> de Welch<\/strong> (avec degr\u00e9s de libert\u00e9 de Satterthwaite), ind\u00e9pendamment du r\u00e9sultat des pr\u00e9-tests. Cette recommandation contemporaine garantit une protection absolue sans co\u00fbt statistique r\u00e9el.<\/li>\n<li><em>Comparaison de trois groupes ou plus :<\/em> Privil\u00e9gier l&rsquo;<strong>ANOVA de Welch<\/strong> associ\u00e9e aux tests post-hoc de <strong>Games-Howell<\/strong>. En pr\u00e9sence de violations conjointes de la normalit\u00e9, mobiliser le bootstrap BCa ou des tests de permutation studentis\u00e9s.<\/li>\n<li><em>R\u00e9gression lin\u00e9aire multiple :<\/em> Impl\u00e9menter de routine les <strong>erreurs-types robustes de type HC3<\/strong> (ou HC4 en pr\u00e9sence d&rsquo;observations \u00e0 fort levier) ou calibrer un mod\u00e8le lin\u00e9aire mixte incluant une fonction de variance explicite (<code>varIdent<\/code>).<\/li>\n<\/ul>\n<\/li>\n<\/ol>\n<h3>12.2 Normes de rapportage selon l&rsquo;APA (American Psychological Association, 7e \u00e9d.)<\/h3>\n<p>L&rsquo;int\u00e9grit\u00e9 de la communication scientifique impose une transparence absolue lors de la restitution des analyses statistiques dans les manuscrits soumis \u00e0 publication. Selon les directives de l&rsquo;<em>American Psychological Association<\/em> (APA, 7e \u00e9dition, Section 6.44), le chercheur ne doit jamais masquer la pr\u00e9sence d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 ni omettre d&rsquo;expliciter les ajustements statistiques appliqu\u00e9s pour y r\u00e9pondre.<\/p>\n<p>La restitution formelle de ces r\u00e9sultats r\u00e9pond \u00e0 des conventions d&rsquo;\u00e9criture typographiques et math\u00e9matiques standardis\u00e9es tr\u00e8s pr\u00e9cises :<\/p>\n<ul>\n<li><strong>Rapportage d&rsquo;un test t de Welch :<\/strong> Il est imp\u00e9ratif de mentionner explicitement que la correction de Welch pour variances in\u00e9gales a \u00e9t\u00e9 appliqu\u00e9e. La statistique <em>t<\/em> doit \u00eatre rapport\u00e9e avec les <strong>degr\u00e9s de libert\u00e9 fractionnaires<\/strong> r\u00e9els issus de la formule d&rsquo;approximation de Satterthwaite, sans jamais les arrondir arbitrairement \u00e0 des entiers :<br \/>\n    <br \/>\n    <em>\u00ab L&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances \u00e9tant compromise, un test t de Welch pour \u00e9chantillons ind\u00e9pendants a \u00e9t\u00e9 conduit. La performance mn\u00e9sique moyenne du groupe exp\u00e9rimental (M = 14,32, SD = 4,87) s&rsquo;est av\u00e9r\u00e9e significativement sup\u00e9rieure \u00e0 celle du groupe t\u00e9moin (M = 11,05, SD = 2,12), t(34,28) = 2,84, p = 0,007, d = 0,86, IC \u00e0 95 % [0,94, 5,60]. \u00bb<\/em><\/li>\n<li><strong>Rapportage d&rsquo;une ANOVA de Welch et contrastes associ\u00e9s :<\/strong><br \/>\n    <br \/>\n    <em>\u00ab En raison d&rsquo;une disparit\u00e9 substantielle de variabilit\u00e9 entre les trois cohortes, une ANOVA robuste de Welch a \u00e9t\u00e9 mobilis\u00e9e. L&rsquo;effet omnibus du niveau d&rsquo;anxi\u00e9t\u00e9 s&rsquo;est r\u00e9v\u00e9l\u00e9 significatif, F_Welch(2, 28,45) = 8,92, p &lt; 0,001, omega^2 = 0,14. Les comparaisons par paires a posteriori conduites \u00e0 l&#039;aide de la proc\u00e9dure de Games-Howell ont indiqu\u00e9 que le groupe \u00e0 anxi\u00e9t\u00e9 s\u00e9v\u00e8re pr\u00e9sentait des temps de r\u00e9ponse significativement plus \u00e9lev\u00e9s que le groupe asymptomatique (p = 0,002) et le groupe mod\u00e9r\u00e9 (p = 0,028). \u00bb<\/em><\/li>\n<li><strong>Rapportage en r\u00e9gression lin\u00e9aire robuste :<\/strong><br \/>\n    <br \/>\n    <em>\u00ab Face \u00e0 la pr\u00e9sence av\u00e9r\u00e9e d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 diagnostiqu\u00e9e par l&rsquo;inspection des r\u00e9sidus vs fitted, les param\u00e8tres de r\u00e9gression ont \u00e9t\u00e9 estim\u00e9s par la m\u00e9thode OLS conjointe \u00e0 l&rsquo;estimateur d&rsquo;erreurs-types robustes HC3 de Davidson-MacKinnon. Le score d&rsquo;inhibition cognitive pr\u00e9dit significativement la fluidit\u00e9 verbale, B = 0,42, SE_HC3 = 0,11, t(118) = 3,82, p &lt; 0,001, IC \u00e0 95 % [0,20, 0,64]. \u00bb<\/em><\/li>\n<\/ul>\n<h3>12.3 Transparence, science ouverte et reproductibilit\u00e9 computationnelle<\/h3>\n<p>Enfin, le traitement ad\u00e9quat de l&rsquo;\u00e9galit\u00e9 des variances s&rsquo;inscrit au c\u0153ur du mouvement plan\u00e9taire pour la science ouverte (<em>Open Science<\/em>) et la reproductibilit\u00e9 computationnelle. La flexibilit\u00e9 dont dispose l&rsquo;analyste face au choix de ses m\u00e9thodes de v\u00e9rification de postulats et d&rsquo;ajustements statistiques constitue une source majeure de degr\u00e9s de libert\u00e9 non document\u00e9s (<em>researcher degrees of freedom<\/em>), susceptibles d&rsquo;alimenter des pratiques de recherche douteuses (<em>QRP<\/em>) telles que le <em>p<\/em>-hacking.<\/p>\n<p>Naviguer entre le test <em>t<\/em> classique, le test de Welch, une transformation logarithmique ou un test non param\u00e9trique jusqu&rsquo;\u00e0 ce que la <em>p<\/em>-valeur bascule opportun\u00e9ment sous le seuil symbolique de 0,05 corrompt l&rsquo;\u00e9difice scientifique. Pour \u00e9radiquer ces biais syst\u00e9miques, plusieurs pratiques vertueuses doivent \u00eatre syst\u00e9matiquement adopt\u00e9es :<\/p>\n<ul>\n<li><strong>Le pr\u00e9-enregistrement des analyses (Preregistration) :<\/strong> Les chercheurs doivent sp\u00e9cifier pr\u00e9alablement sur des plateformes d\u00e9di\u00e9es (telles que le <em>Center for Open Science<\/em> &#8211; OSF) leur plan analytique d&rsquo;investigation. Ce document doit consigner formellement la d\u00e9cision d&#8217;employer par d\u00e9faut les proc\u00e9dures robustes (comme le test de Welch ou les erreurs-types HC3), coupant court \u00e0 toute tentation d&rsquo;arbitrage post-hoc.<\/li>\n<li><strong>Le partage int\u00e9gral des scripts computationnels :<\/strong> Les scripts d&rsquo;analyse r\u00e9dig\u00e9s sous des environnements de programmation statistique ouverts (R, Python) ou des logiciels libres pilot\u00e9s par menus transparents (JASP, Jamovi) doivent \u00eatre d\u00e9pos\u00e9s en acc\u00e8s public dans des d\u00e9p\u00f4ts p\u00e9rennes. Le code source doit consigner sans fard l&rsquo;int\u00e9gralit\u00e9 du pipeline m\u00e9thodologique : scripts de diagnostic visuel, calculs d&rsquo;erreurs-types robustes et mod\u00e9lisations alternatives d&rsquo;analyse de sensibilit\u00e9.<\/li>\n<li><strong>La publication d&rsquo;analyses de sensibilit\u00e9 multivers :<\/strong> Lorsque la fronti\u00e8re entre homosc\u00e9dasticit\u00e9 et h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 est ambigu\u00eb, rapporter en annexes ou en mat\u00e9riel suppl\u00e9mentaire la matrice compl\u00e8te des r\u00e9sultats obtenus sous diff\u00e9rents prismes analytiques (test classique, test de Welch, r\u00e9gression robuste, bootstrap) conf\u00e8re aux conclusions de l&rsquo;article une robustesse et une autorit\u00e9 empirique in\u00e9branlables.<\/li>\n<\/ul>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<p>Box, G. E. P. (1954). Some theorems on quadratic forms applied in the study of analysis of variance problems, I. Effect of inequality of variance in the one-way classification. <em>The Annals of Mathematical Statistics<\/em>, 25(2), 290\u2013302. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aoms\/1177728786\">https:\/\/doi.org\/10.1214\/aoms\/1177728786<\/a><\/p>\n<p>Box, G. E. P., &amp; Cox, D. R. (1964). An analysis of transformations. <em>Journal of the Royal Statistical Society: Series B (Methodological)<\/em>, 26(2), 211\u2013243. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x\">https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x<\/a><\/p>\n<p>Brown, M. B., &amp; Forsythe, A. B. (1974). Robust tests for the equality of variances. <em>Journal of the American Statistical Association<\/em>, 69(346), 364\u2013367. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1974.10482955\">https:\/\/doi.org\/10.1080\/01621459.1974.10482955<\/a><\/p>\n<p>Cribari-Neto, F. (2004). Asymptotic inference under heteroskedasticity of unknown form. <em>Computational Statistics &amp; Data Analysis<\/em>, 45(2), 215\u2013233. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/S0167-9473(02)00366-3\">https:\/\/doi.org\/10.1016\/S0167-9473(02)00366-3<\/a><\/p>\n<p>Davidson, R., &amp; MacKinnon, J. G. (1993). <em>Estimation and inference in econometrics<\/em>. Oxford University Press.<\/p>\n<p>Delacre, M., Lakens, D., &amp; Leys, C. (2017). Why psychologists should by default use Welch\u2019s t-test instead of Student\u2019s t-test. <em>International Review of Social Psychology<\/em>, 30(1), 92\u2013101. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.5334\/irsp.82\">https:\/\/doi.org\/10.5334\/irsp.82<\/a><\/p>\n<p>Efron, B. (1979). Bootstrap methods: Another look at the jackknife. <em>The Annals of Statistics<\/em>, 7(1), 1\u201326. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aos\/1176344552\">https:\/\/doi.org\/10.1214\/aos\/1176344552<\/a><\/p>\n<p>Fisher, R. A. (1925). <em>Statistical methods for research workers<\/em>. Oliver and Boyd.<\/p>\n<p>Fligner, M. A., &amp; Killeen, T. J. (1976). Distribution-free two-sample tests for scale. <em>Journal of the American Statistical Association<\/em>, 71(353), 210\u2013213. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1976.10481517\">https:\/\/doi.org\/10.1080\/01621459.1976.10481517<\/a><\/p>\n<p>Games, P. A., &amp; Howell, J. F. (1976). Pairwise multiple comparison procedures with unequal n&rsquo;s and\/or variances: A Monte Carlo study. <em>Journal of Educational Statistics<\/em>, 1(2), 113\u2013125. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1164979\">https:\/\/doi.org\/10.2307\/1164979<\/a><\/p>\n<p>Hedeker, D., Mermelstein, R. J., &amp; Demirtas, H. (2008). An application of a mixed-effects location scale model for analysis of ecological momentary assessment (EMA) data. <em>Biometrics<\/em>, 64(2), 627\u2013634. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.1541-0420.2007.00898.x\">https:\/\/doi.org\/10.1111\/j.1541-0420.2007.00898.x<\/a><\/p>\n<p>Levene, H. (1960). Robust tests for equality of variances. In I. Olkin, S. G. Ghurye, W. Hoeffding, W. G. Madow, &amp; H. B. Mann (Eds.), <em>Contributions to probability and statistics: Essays in honor of Harold Hotelling<\/em> (pp. 278\u2013292). Stanford University Press.<\/p>\n<p>Long, J. S., &amp; Ervin, L. H. (2000). Using heteroscedasticity consistent standard errors in the linear regression model. <em>The American Statistician<\/em>, 54(3), 217\u2013224. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00031305.2000.10474549\">https:\/\/doi.org\/10.1080\/00031305.2000.10474549<\/a><\/p>\n<p>MacKinnon, J. G., &amp; White, H. (1985). Some heteroskedasticity-consistent covariance matrix estimators with improved finite sample properties. <em>Journal of Econometrics<\/em>, 29(3), 305\u2013325. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/0304-4076(85)90158-7\">https:\/\/doi.org\/10.1016\/0304-4076(85)90158-7<\/a><\/p>\n<p>Mauchly, J. W. (1940). Significance test for sphericity of a normal n-variate distribution. <em>The Annals of Mathematical Statistics<\/em>, 11(2), 204\u2013209. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aoms\/1177731909\">https:\/\/doi.org\/10.1214\/aoms\/1177731909<\/a><\/p>\n<p>Preacher, K. J., &amp; Hayes, A. F. (2008). Asymptotic and resampling strategies for assessing and comparing indirect effects in multiple mediator models. <em>Behavior Research Methods<\/em>, 40(3), 879\u2013891. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.3758\/BRM.40.3.879\">https:\/\/doi.org\/10.3758\/BRM.40.3.879<\/a><\/p>\n<p>Ruxton, G. D. (2006). The unequal variance t-test is underused. <em>Behavioral Ecology<\/em>, 17(4), 688\u2013690. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/beheco\/ark016\">https:\/\/doi.org\/10.1093\/beheco\/ark016<\/a><\/p>\n<p>Student. (1908). The probable error of a mean. <em>Biometrika<\/em>, 6(1), 1\u201325. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/6.1.1\">https:\/\/doi.org\/10.1093\/biomet\/6.1.1<\/a><\/p>\n<p>Welch, B. L. (1938). The significance of the difference between two means when the population variances are unequal. <em>Biometrika<\/em>, 29(3\/4), 350\u2013362. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/29.3-4.350\">https:\/\/doi.org\/10.1093\/biomet\/29.3-4.350<\/a><\/p>\n<p>Welch, B. L. (1947). The generalization of &lsquo;Student&rsquo;s&rsquo; problem when several different population variances are involved. <em>Biometrika<\/em>, 34(1\/2), 28\u201335. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/34.1-2.28\">https:\/\/doi.org\/10.1093\/biomet\/34.1-2.28<\/a><\/p>\n<p>Welch, B. L. (1951). On the comparison of several mean values: An alternative approach. <em>Biometrika<\/em>, 38(3\/4), 330\u2013336. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/38.3-4.330\">https:\/\/doi.org\/10.1093\/biomet\/38.3-4.330<\/a><\/p>\n<p>White, H. (1980). A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity. <em>Econometrica<\/em>, 48(4), 817\u2013838. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1912934\">https:\/\/doi.org\/10.2307\/1912934<\/a><\/p>\n<p>Yeo, I. K., &amp; Johnson, R. A. (2000). A new family of power transformations to improve normality or symmetry. <em>Biometrika<\/em>, 87(4), 954\u2013959. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/87.4.954\">https:\/\/doi.org\/10.1093\/biomet\/87.4.954<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Comprendre l&rsquo;hypoth\u00e8se d&rsquo;\u00e9galit\u00e9 des variances en statistique : d\u00e9finition, tests associ\u00e9s (ANOVA, r\u00e9gression), m\u00e9thodes de diagnostic et solutions pratiques.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1907,1858],"tags":[5167,2461,1942,2460,3321],"class_list":["post-8105","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-et-statistiques","category-psychometrie-et-recherche-quantitative","tag-anova-de-welch","tag-homoscedasticite","tag-statistiques-en-psychologie","tag-test-de-levene","tag-test-t-de-student"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8105","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat"}],"about":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/types\/psych_stat"}],"author":[{"embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=8105"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8105\/revisions"}],"predecessor-version":[{"id":8108,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8105\/revisions\/8108"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=8105"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=8105"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=8105"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}