{"id":8229,"date":"2026-09-17T04:45:05","date_gmt":"2026-09-17T04:45:05","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/anova-welch-python-etape-par-etape\/"},"modified":"2026-09-17T04:45:05","modified_gmt":"2026-09-17T04:45:05","slug":"anova-welch-python-etape-par-etape","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=anova-welch-python-etape-par-etape","title":{"rendered":"ANOVA de Welch en Python (\u00e9tape par \u00e9tape)"},"content":{"rendered":"<p>Dans le champ de la psychologie exp\u00e9rimentale et des neurosciences cognitives contemporaines, la rigueur de l&rsquo;inf\u00e9rence statistique constitue le rempart fondamental contre les crises de reproductibilit\u00e9. Les chercheurs soumettent r\u00e9guli\u00e8rement des cohortes de participants \u00e0 des paradigmes d&rsquo;apprentissage, des protocoles pharmacologiques ou des inductions \u00e9motionnelles, dans l&rsquo;espoir d&rsquo;isoler des variations comportementales statistiquement significatives. Pour tester les diff\u00e9rences entre les moyennes de plusieurs groupes ind\u00e9pendants, le r\u00e9flexe m\u00e9thodologique standard demeure l&rsquo;analyse de variance \u00e0 un facteur (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Analysis_of_variance\">ANOVA<\/a>), conceptualis\u00e9e \u00e0 l&rsquo;origine par Sir Ronald Fisher. Or, l&rsquo;utilisation routini\u00e8re de cette proc\u00e9dure repose sur des postulats math\u00e9matiques stricts qui se trouvent fr\u00e9quemment viol\u00e9s lors du recueil de donn\u00e9es empiriques sur le comportement humain.<\/p>\n<p>Parmi ces postulats, l&rsquo;homog\u00e9n\u00e9it\u00e9 des variances intra-groupes \u2014 d\u00e9sign\u00e9e sous le terme d&rsquo;homosc\u00e9dasticit\u00e9 \u2014 s&rsquo;av\u00e8re particuli\u00e8rement vuln\u00e9rable aux singularit\u00e9s des protocoles exp\u00e9rimentaux. Lorsque des groupes de sujets manifestent des niveaux de dispersion asym\u00e9triques, l&rsquo;estimation conjointe de la variance d&rsquo;erreur devient biais\u00e9e, ce qui d\u00e9stabilise radicalement le taux de fausses d\u00e9couvertes (l&rsquo;erreur de premi\u00e8re esp\u00e8ce, not\u00e9e &alpha;). Dans de telles conditions, persister dans l&rsquo;application de l&rsquo;ANOVA classique de Fisher expose l&rsquo;analyste \u00e0 des conclusions fallacieuses, soit par exc\u00e8s d&rsquo;optimisme statistique (comportement lib\u00e9ral du test), soit par perte dramatique de puissance d&rsquo;\u00e9valuation (comportement excessivement conservateur). Face \u00e0 cette impasse, le recours syst\u00e9matique \u00e0 des tests non param\u00e9triques bas\u00e9s sur les rangs ne constitue pas toujours une r\u00e9ponse optimale, car ces derniers modifient la nature m\u00eame de l&rsquo;hypoth\u00e8se test\u00e9e en d\u00e9laissant la comparaison directe des moyennes param\u00e9triques au profit de distributions ordinales.<\/p>\n<p>C&rsquo;est pr\u00e9cis\u00e9ment au c\u0153ur de cette tension m\u00e9thodologique que s&rsquo;impose l&rsquo;ANOVA de Welch. \u00c9labor\u00e9e pour s&rsquo;affranchir du postulat d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances sans sacrifier l&rsquo;\u00e9valuation directe des moyennes empiriques, cette approche param\u00e9trique robuste modifie le calcul du ratio d&rsquo;amplitudes et ajuste de mani\u00e8re dynamique les degr\u00e9s de libert\u00e9 r\u00e9siduels. Ce guide complet explore les fondements th\u00e9oriques, les justifications psychom\u00e9triques et le d\u00e9ploiement computationnel pas \u00e0 pas de l&rsquo;ANOVA de Welch au sein de l&rsquo;\u00e9cosyst\u00e8me Python. En s&rsquo;appuyant sur des biblioth\u00e8ques scientifiques de r\u00e9f\u00e9rence telles que NumPy, SciPy, Pandas et <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pingouin-stats.org\/\">Pingouin<\/a>, ce tutoriel d\u00e9taille l&rsquo;int\u00e9gralit\u00e9 du pipeline d&rsquo;analyse : du diagnostic exploratoire des variances aux proc\u00e9dures post-hoc adapt\u00e9es, jusqu&rsquo;\u00e0 la restitution formelle des r\u00e9sultats selon les normes de l&rsquo;American Psychological Association (APA 7e \u00e9dition).<\/p>\n<h2>1. Introduction th\u00e9orique \u00e0 l&rsquo;ANOVA de Welch en psychologie exp\u00e9rimentale<\/h2>\n<p>L&rsquo;inf\u00e9rence statistique moderne appliqu\u00e9e aux sciences du comportement exige une ad\u00e9quation parfaite entre la mod\u00e9lisation math\u00e9matique et la variabilit\u00e9 intrins\u00e8que des organismes vivants. L&rsquo;analyse de variance \u00e0 un facteur selon le mod\u00e8le de Fisher constitue historiquement l&rsquo;outil canonique pour comparer simultan\u00e9ment les moyennes de trois groupes ou davantage. Toutefois, sa validit\u00e9 asymptotique repose sur l&rsquo;hypoth\u00e8se d&rsquo;une dispersion identique au sein de chaque population parente. Lorsque cette condition s&rsquo;effondre, l&rsquo;ANOVA de Welch appara\u00eet comme le cadre d&rsquo;analyse param\u00e9trique de r\u00e9f\u00e9rence, offrant une r\u00e9silience math\u00e9matique prouv\u00e9e face aux distorsions des variances empiriques.<\/p>\n<h3>1.1 D\u00e9finition et gen\u00e8se de l&rsquo;ANOVA de Welch<\/h3>\n<p>L&rsquo;ANOVA de Welch tire son origine des travaux fondateurs du statisticien britannique Bernard Lewis Welch, publi\u00e9s successivement en 1947 et 1951 dans la revue <em>Biometrika<\/em>. Welch s&rsquo;est attel\u00e9 \u00e0 r\u00e9soudre une extension multigroupes de l&rsquo;\u00e9pineux probl\u00e8me de Behrens-Fisher, lequel s&rsquo;int\u00e9resse \u00e0 la comparaison de deux populations normales dont le rapport des variances est inconnu et arbitraire. Jusqu&rsquo;alors, la pratique courante consistait \u00e0 supposer l&rsquo;\u00e9galit\u00e9 des variances pour estimer un terme d&rsquo;erreur commun (la variance r\u00e9siduelle intra-groupe pond\u00e9r\u00e9e), une d\u00e9marche popularis\u00e9e par Ronald A. Fisher d\u00e8s les ann\u00e9es 1920 dans le cadre de la mod\u00e9lisation agronomique.<\/p>\n<p>Dans ses publications fondamentales, Welch d\u00e9montre que le regroupement na\u00eff des sommes des carr\u00e9s r\u00e9siduelles sous-estime ou surestime gravement la variance de la diff\u00e9rence des moyennes d\u00e8s lors que les populations sous-jacentes divergent quant \u00e0 leur dispersion. Pour corriger cette d\u00e9faillance, il introduit une statistique de test alternative fond\u00e9e sur une pond\u00e9ration par l&rsquo;inverse des variances d&rsquo;\u00e9chantillonnage de chaque groupe. Au lieu de postuler un d\u00e9nominateur universel pour le ratio <em>F<\/em>, la m\u00e9thode de Welch recalcule la variance globale \u00e0 partir des contributions individuelles de chaque modalit\u00e9 factorielle.<\/p>\n<p>Dans l&rsquo;arsenal inf\u00e9rentiel moderne, l&rsquo;ANOVA de Welch ne doit plus \u00eatre per\u00e7ue comme un simple m\u00e9canisme de secours r\u00e9serv\u00e9 aux donn\u00e9es d\u00e9grad\u00e9es, mais plut\u00f4t comme une g\u00e9n\u00e9ralisation g\u00e9n\u00e9ralis\u00e9e et robuste de l&rsquo;ANOVA de Fisher. Elle s&rsquo;impose aujourd&rsquo;hui dans la litt\u00e9rature m\u00e9thodologique comme le choix par d\u00e9faut pour tout devis unifactoriel inter-sujets, dans la mesure o\u00f9 elle converge exactement vers les r\u00e9sultats de l&rsquo;ANOVA classique lorsque les variances s&rsquo;av\u00e8rent homog\u00e8nes, tout en pr\u00e9servant le chercheur des artefacts m\u00e9thodologiques en pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 s\u00e9v\u00e8re.<\/p>\n<h3>1.2 Pertinence m\u00e9thodologique en sciences psychologiques<\/h3>\n<p>La recherche en psychologie exp\u00e9rimentale, en neuropsychologie et en sciences de l&rsquo;\u00e9ducation confronte syst\u00e9matiquement le chercheur \u00e0 une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 structurelle des mesures comportementales. Contrairement aux sciences physiques o\u00f9 les instruments b\u00e9n\u00e9ficient d&rsquo;une pr\u00e9cision m\u00e9trologique stable, les variables psychologiques (temps de r\u00e9action en millisecondes, scores de r\u00e9tention mn\u00e9sique, taux d&rsquo;erreurs cognitives ou \u00e9chelles d&rsquo;anxi\u00e9t\u00e9) refl\u00e8tent une sensibilit\u00e9 diff\u00e9rentielle aigu\u00eb des participants face aux manipulations de laboratoire.<\/p>\n<p>Cette variabilit\u00e9 n&rsquo;est presque jamais uniforme entre les conditions exp\u00e9rimentales. Consid\u00e9rons un protocole clinique \u00e9valuant l&rsquo;efficacit\u00e9 d&rsquo;une rem\u00e9diation cognitive : le groupe t\u00e9moin non entra\u00een\u00e9 peut afficher une dispersion restreinte autour d&rsquo;une moyenne de performance basse, tandis que le groupe soumis \u00e0 un protocole intensif peut pr\u00e9senter une variance d\u00e9cupl\u00e9e en raison de la pr\u00e9sence simultan\u00e9e de sujets tr\u00e8s r\u00e9pondeurs et de sujets r\u00e9fractaires \u00e0 l&rsquo;intervention. Cette h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances n&rsquo;est pas un artefact de mesure ou un bruit ind\u00e9sirable, mais l&rsquo;expression directe d&rsquo;un ph\u00e9nom\u00e8ne psychologique sous-jacent.<\/p>\n<p>En neurosciences cognitives, les protocoles inter-sujets mobilisent fr\u00e9quemment des \u00e9chantillons de taille modeste pour des raisons de co\u00fbt et de faisabilit\u00e9 logistique. L&rsquo;interaction entre ces faibles effectifs et des variances in\u00e9gales induit un risque intol\u00e9rable d&rsquo;inflation de l&rsquo;erreur de premi\u00e8re esp\u00e8ce (&alpha;). En conservant l&rsquo;ANOVA classique, le taux nominal d&rsquo;erreur fix\u00e9 conventionnellement \u00e0 0,05 peut glisser vers des valeurs empiriques r\u00e9elles exc\u00e9dant 0,15, entra\u00eenant la publication de faux positifs th\u00e9oriques majeurs. L&rsquo;ANOVA de Welch stabilise rigoureusement ce taux d&rsquo;erreur au seuil prescrit, s\u00e9curisant ainsi l&rsquo;\u00e9difice probatoire des d\u00e9couvertes exp\u00e9rimentales.<\/p>\n<h3>1.3 Comparaison synth\u00e9tique : Fisher, Welch et alternatives non param\u00e9triques<\/h3>\n<p>Face \u00e0 une violation flagrante de l&rsquo;homosc\u00e9dasticit\u00e9, les psychologues h\u00e9sitent fr\u00e9quemment entre trois trajectoires analytiques : persister avec l&rsquo;ANOVA classique de Fisher, bifurquer vers un test non param\u00e9trique tel que le test de <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Kruskal%E2%80%93Wallis_one-way_analysis_of_variance\">Kruskal-Wallis<\/a>, ou mobiliser la correction param\u00e9trique de Welch. Cette h\u00e9sitation d\u00e9coule souvent d&rsquo;une m\u00e9compr\u00e9hension des postulats et des pertes d&rsquo;informations associ\u00e9es \u00e0 chaque m\u00e9thode.<\/p>\n<p>L&rsquo;ANOVA classique de Fisher pr\u00e9sente une vuln\u00e9rabilit\u00e9 critique : elle calcule le carr\u00e9 moyen des erreurs en agr\u00e9geant l&rsquo;ensemble des \u00e9carts r\u00e9siduels, postulant que chaque observation apporte la m\u00eame pr\u00e9cision sur le param\u00e8tre de dispersion. Si ce postulat est infirm\u00e9, la distribution d&rsquo;\u00e9chantillonnage de la statistique <em>F<\/em> ne suit plus une loi de Fisher th\u00e9orique, rendant la <em>p<\/em>-valeur calcul\u00e9e math\u00e9matiquement erron\u00e9e. Le test de Kruskal-Wallis, quant \u00e0 lui, contourne ce probl\u00e8me en convertissant les observations brutes en rangs ordinaux. Si cette conversion neutralise l&rsquo;impact des valeurs aberrantes extr\u00eames, elle transforme \u00e9galement la question de recherche : le test n&rsquo;\u00e9value plus une diff\u00e9rence de moyennes arithm\u00e9tiques, mais une dominance stochastique ou une translation des m\u00e9dianes (sous l&rsquo;hypoth\u00e8se contestable que les distributions partagent une forme identique).<\/p>\n<p>L&rsquo;ANOVA de Welch se positionne comme un compromis optimal. En conservant l&rsquo;\u00e9chelle de mesure originale continue, elle pr\u00e9serve la totalit\u00e9 de l&rsquo;information m\u00e9trique contenue dans les scores des participants. De multiples \u00e9tudes par simulations de Monte-Carlo (notamment celles de Tomarken et Serlin en 1986, ou plus r\u00e9cemment de Wilcox en 2012) ont d\u00e9montr\u00e9 de mani\u00e8re unanime que le test de Welch conserve une puissance statistique remarquable, rivalisant avec celle de Fisher sous parfaite normalit\u00e9 et homosc\u00e9dasticit\u00e9, tout en surpassant le test de Kruskal-Wallis d\u00e8s lors que les donn\u00e9es sont continues et pr\u00e9sentent des d\u00e9viations mod\u00e9r\u00e9es d&rsquo;asym\u00e9trie.<\/p>\n<h2>2. Les postulats statistiques de l&rsquo;ANOVA classique et leurs limites<\/h2>\n<p>Comprendre la sup\u00e9riorit\u00e9 de l&rsquo;ANOVA de Welch n\u00e9cessite d&rsquo;analyser au pr\u00e9alable l&rsquo;architecture formelle de l&rsquo;ANOVA classique de Fisher et d&rsquo;identifier pr\u00e9cis\u00e9ment les points de rupture de ses th\u00e9or\u00e8mes fondateurs. L&rsquo;ANOVA s&rsquo;inscrit dans le cadre du mod\u00e8le lin\u00e9aire g\u00e9n\u00e9ral, lequel repose sur une s\u00e9rie de contraintes distributionnelles strictes r\u00e9gissant les r\u00e9sidus d&rsquo;\u00e9chantillonnage.<\/p>\n<h3>2.1 L&rsquo;hypoth\u00e8se fondamentale d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances<\/h3>\n<p>Le mod\u00e8le lin\u00e9aire unifactoriel inter-sujets s&rsquo;\u00e9nonce math\u00e9matiquement sous la forme suivante :<\/p>\n<p><em>Y<sub>ij<\/sub> = &mu; + &alpha;<sub>j<\/sub> + &epsilon;<sub>ij<\/sub><\/em><\/p>\n<p>O\u00f9 <em>Y<sub>ij<\/sub><\/em> repr\u00e9sente la mesure obtenue par le <em>i<\/em>-i\u00e8me individu au sein du <em>j<\/em>-i\u00e8me groupe exp\u00e9rimental, <em>&mu;<\/em> correspond \u00e0 la moyenne g\u00e9n\u00e9rale th\u00e9orique, <em>&alpha;<sub>j<\/sub><\/em> mat\u00e9rialise l&rsquo;effet sp\u00e9cifique du traitement administr\u00e9 au groupe <em>j<\/em>, et <em>&epsilon;<sub>ij<\/sub><\/em> symbolise l&rsquo;erreur al\u00e9atoire r\u00e9siduelle. Pour d\u00e9river analytiquement la distribution de la statistique de test <em>F<\/em> selon les lois \u00e9tablies par Fisher, la th\u00e9orie exige imp\u00e9rativement que les erreurs r\u00e9siduelles <em>&epsilon;<sub>ij<\/sub><\/em> soient mutuellement ind\u00e9pendantes et distribu\u00e9es de fa\u00e7on identique et normale selon une esp\u00e9rance nulle et une variance constante <em>&sigma;<sup>2<\/sup><\/em>, quel que soit le groupe d&rsquo;appartenance <em>j<\/em> :<\/p>\n<p><em>&epsilon;<sub>ij<\/sub> ~ N(0, &sigma;<sup>2<\/sup>) &forall; j &isin; {1, &#8230;, k}<\/em><\/p>\n<p>Cette condition d&rsquo;\u00e9galit\u00e9 stricte des variances r\u00e9siduelles constitue l&rsquo;homosc\u00e9dasticit\u00e9. Au sein de la fraction d\u00e9finissant le ratio <em>F<\/em> de Fisher, le num\u00e9rateur estime la dispersion inter-groupes (due \u00e0 l&rsquo;effet potentiel du facteur manipul\u00e9 augment\u00e9 de la variance d&rsquo;erreur), tandis que le d\u00e9nominateur estime exclusivement la variance r\u00e9siduelle intra-groupe (l&rsquo;erreur d&rsquo;\u00e9chantillonnage non expliqu\u00e9e). Pour obtenir ce d\u00e9nominateur, l&rsquo;ANOVA standard mutualise les variances observ\u00e9es <em>s<sub>j<\/sub><sup>2<\/sup><\/em> au prorata de leurs degr\u00e9s de libert\u00e9 respectifs :<\/p>\n<p><em>MS<sub>erreur<\/sub> = [&sum; (n<sub>j<\/sub> &#8211; 1)s<sub>j<\/sub><sup>2<\/sup>] \/ [&sum; (n<sub>j<\/sub> &#8211; 1)]<\/em><\/p>\n<p>Si la variance th\u00e9orique sous-jacente <em>&sigma;<sub>j<\/sub><sup>2<\/sup><\/em> fluctue substantiellement d&rsquo;une condition \u00e0 l&rsquo;autre, cette mutualisation perd toute consistance th\u00e9orique. Le d\u00e9nominateur du ratio <em>F<\/em> cesse de repr\u00e9senter ad\u00e9quatement l&rsquo;erreur type attendue sous l&rsquo;hypoth\u00e8se nulle, invalidant de fait la structure probabiliste de la d\u00e9cision inf\u00e9rentielle.<\/p>\n<h3>2.2 Vuln\u00e9rabilit\u00e9 de l&rsquo;ANOVA de Fisher face aux effectifs in\u00e9gaux<\/h3>\n<p>La d\u00e9gradation des propri\u00e9t\u00e9s statistiques de l&rsquo;ANOVA de Fisher atteint son niveau de criticit\u00e9 le plus d\u00e9vastateur lorsque l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances coexiste avec des tailles d&rsquo;\u00e9chantillons in\u00e9gales (d\u00e9s\u00e9quilibre exp\u00e9rimental fr\u00e9quent provoqu\u00e9 par l&rsquo;attrition des sujets, les abandons en cours d&rsquo;\u00e9tude ou les exclusions de donn\u00e9es pour non-conformit\u00e9). L&rsquo;impact de cette interaction sur le taux nominal d&rsquo;erreur de premi\u00e8re esp\u00e8ce d\u00e9pend directement du patron d&rsquo;association entre la taille des groupes et l&rsquo;amplitude des variances.<\/p>\n<p>Deux trajectoires n\u00e9fastes sont syst\u00e9matiquement document\u00e9es :<\/p>\n<ul>\n<li><strong>Le patron d&rsquo;appariement direct (grands groupes associ\u00e9s \u00e0 de fortes variances) :<\/strong> Lorsque les cohortes poss\u00e9dant les effectifs <em>n<sub>j<\/sub><\/em> les plus \u00e9lev\u00e9s pr\u00e9sentent simultan\u00e9ment les variances <em>s<sub>j<\/sub><sup>2<\/sup><\/em> les plus larges, le terme r\u00e9siduel mutualis\u00e9 <em>MS<sub>erreur<\/sub><\/em> se trouve sur\u00e9valu\u00e9. La variance mutualis\u00e9e est en effet tir\u00e9e vers le haut par les grands effectifs. En cons\u00e9quence, le d\u00e9nominateur du ratio <em>F<\/em> devient artificiellement grand, ce qui abaisse ind\u00fbment la valeur du <em>F<\/em> calcul\u00e9. Le test devient alors excessivement <strong>conservateur<\/strong> : le taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce r\u00e9el chute bien en de\u00e7\u00e0 du seuil &alpha; fix\u00e9 (par exemple 0,01 au lieu de 0,05), ce qui an\u00e9antit la puissance statistique du dispositif et condamne l&rsquo;exp\u00e9rimentateur \u00e0 rater des effets r\u00e9els (augmentation drastique de l&rsquo;erreur de deuxi\u00e8me esp\u00e8ce, not\u00e9e &beta;).<\/li>\n<li><strong>Le patron d&rsquo;appariement inverse (petits groupes associ\u00e9s \u00e0 de fortes variances) :<\/strong> \u00c0 l&rsquo;inverse, si les petits groupes d&rsquo;individus concentrent les variances les plus massives, le terme mutualis\u00e9 <em>MS<sub>erreur<\/sub><\/em> est sous-\u00e9valu\u00e9, car l&rsquo;estimation de l&rsquo;erreur est domin\u00e9e par les variances restreintes des groupes \u00e0 grands effectifs. Le d\u00e9nominateur du ratio <em>F<\/em> est alors artificiellement comprim\u00e9, gonflant le ratio <em>F<\/em> \u00e0 des niveaux arbitrairement \u00e9lev\u00e9s. Le test devient <strong>lib\u00e9ral<\/strong> : le taux r\u00e9el de faux positifs s&rsquo;envole fr\u00e9quemment vers 10 %, 15 % voire 25 % pour un seuil nominal initialement calibr\u00e9 \u00e0 5 %. Les conclusions scientifiques formul\u00e9es sur de telles bases s&rsquo;av\u00e8rent fallacieuses.<\/li>\n<\/ul>\n<h3>2.3 Le postulat de normalit\u00e9 distributionnelle<\/h3>\n<p>L&rsquo;autre pilier axiomatique de l&rsquo;ANOVA param\u00e9trique repose sur la normalit\u00e9 distributionnelle des scores ou des r\u00e9sidus au sein de chaque population parente. Dans l&rsquo;imaginaire m\u00e9thodologique, le Th\u00e9or\u00e8me Central Limite (TCL) est souvent invoqu\u00e9 de mani\u00e8re indue pour postuler qu&rsquo;un effectif d&rsquo;environ 30 participants par condition suffirait \u00e0 absoudre toute violation de ce postulat. S&rsquo;il est exact que la distribution des moyennes d&rsquo;\u00e9chantillonnage tend asymptotiquement vers une loi gaussienne \u00e0 mesure que <em>N<\/em> cro\u00eet, cette convergence ne garantit pas la stabilit\u00e9 du ratio <em>F<\/em> lorsque des asym\u00e9tries distributionnelles s\u00e9v\u00e8res cohabitent avec l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9.<\/p>\n<p>L&rsquo;interaction entre une distribution r\u00e9siduelle fortement asym\u00e9trique (pr\u00e9sence d&rsquo;une queue de distribution \u00e9tir\u00e9e vers les valeurs positives ou n\u00e9gatives, courante dans les mesures de temps de latence cognitive) et un d\u00e9s\u00e9quilibre prononc\u00e9 des variances d\u00e9truit la robustesse relative de l&rsquo;ANOVA de Fisher. Lorsque les courbes de densit\u00e9 empiriques pr\u00e9sentent simultan\u00e9ment des d\u00e9viations de sym\u00e9trie et des amplitudes de dispersion disparates, la valeur critique tabul\u00e9e dans la loi de Fisher ne correspond plus aux quantiles r\u00e9els observ\u00e9s sous l&rsquo;hypoth\u00e8se nulle. La prise en charge math\u00e9matique de la variance diff\u00e9rentielle par la formulation de Welch offre \u00e9galement une meilleure r\u00e9silience face \u00e0 ces l\u00e9g\u00e8res entorses \u00e0 la normalit\u00e9 pure.<\/p>\n<h2>3. L&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 : concepts, causes et r\u00e9percussions en psychom\u00e9trie<\/h2>\n<p>Pour l&rsquo;exp\u00e9rimentateur, l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances ne doit pas \u00eatre appr\u00e9hend\u00e9e comme une simple anomalie math\u00e9matique \u00e0 corriger m\u00e9caniquement, mais comme le reflet substantiel de m\u00e9canismes psychologiques et cognitifs op\u00e9rant au sein du dispositif exp\u00e9rimental. La psychom\u00e9trie th\u00e9orique fournit de multiples explications causales \u00e0 l&rsquo;apparition de structures de variance divergentes entre les conditions d&rsquo;une recherche.<\/p>\n<h3>3.1 Origines psychologiques de l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 de variance<\/h3>\n<p>Dans les t\u00e2ches cognitives complexes, l&rsquo;exposition \u00e0 des conditions environnementales variables modifie la structure m\u00eame des diff\u00e9rences inter-individuelles. Plusieurs ph\u00e9nom\u00e8nes psychom\u00e9triques majeurs g\u00e9n\u00e8rent structurellement de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 :<\/p>\n<p>En premier lieu, les <strong>diff\u00e9rences individuelles de r\u00e9activit\u00e9 cognitive<\/strong> expliquent pourquoi certaines interventions font exploser la variance d&rsquo;un groupe. Dans une t\u00e2che d&rsquo;apprentissage mobilisant la m\u00e9moire de travail, l&rsquo;introduction d&rsquo;un distracteur sonore intense peut n&rsquo;avoir aucun impact chez des sujets dot\u00e9s de capacit\u00e9s ex\u00e9cutives de haut niveau, tout en provoquant un effondrement complet des performances chez des sujets vuln\u00e9rables. Alors que le groupe t\u00e9moin (travaillant en milieu silencieux) pr\u00e9sentera des performances homog\u00e8nes, le groupe exp\u00e9rimental verra sa variance d\u00e9cupl\u00e9e en raison de cette divergence de trajectoires adaptatives.<\/p>\n<p>En second lieu, les <strong>effets de plancher et de plafond<\/strong> affectent directement la m\u00e9trologie des tests standardis\u00e9s. Si un examen psychom\u00e9trique s&rsquo;av\u00e8re excessivement difficile pour une condition de contr\u00f4le non entra\u00een\u00e9e, les scores des participants s&rsquo;agglutineront massivement au voisinage du score minimal (effet de plancher), comprimant artificiellement la variance r\u00e9siduelle de cette cohorte. Inversement, si une m\u00e9thode d&rsquo;apprentissage optimis\u00e9e hisse la plupart des participants au score parfait sans saturer pour quelques individus en difficult\u00e9, une asym\u00e9trie combin\u00e9e \u00e0 une d\u00e9compression de la variance \u00e9mergera in\u00e9vitablement.<\/p>\n<p>Enfin, la <strong>multiplicit\u00e9 des strat\u00e9gies cognitives<\/strong> g\u00e9n\u00e8re de la variance qualitative. Plac\u00e9s face \u00e0 une consigne ouverte de r\u00e9solution de probl\u00e8mes spatiaux, certains sujets recourent spontan\u00e9ment \u00e0 des algorithmes analytiques pr\u00e9cis tandis que d&rsquo;autres s&rsquo;en remettent \u00e0 des heuristiques visuo-spatiales approximatives. Cette pluralit\u00e9 strat\u00e9gique induit une dispersion intra-groupe consid\u00e9rablement plus large que celle observ\u00e9e sous des consignes ferm\u00e9es et rigides contraignant le champ op\u00e9ratoire.<\/p>\n<h3>3.2 Diagnostic exploratoire de l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9<\/h3>\n<p>Avant d&rsquo;engager toute mod\u00e9lisation statistique formelle, une phase d&rsquo;inspection visuelle et exploratoire s&rsquo;av\u00e8re indispensable pour appr\u00e9hender la g\u00e9om\u00e9trie des dispersions. Plusieurs outils graphiques et m\u00e9triques descriptives permettent de sonder l&rsquo;homog\u00e9n\u00e9it\u00e9 des donn\u00e9es :<\/p>\n<p>Les diagrammes en bo\u00eete \u00e0 moustaches (<em>boxplots<\/em>) et les diagrammes de dispersion des r\u00e9sidus constituent la premi\u00e8re ligne de d\u00e9tection. Dans un boxplot, la hauteur de la bo\u00eete (qui figure l&rsquo;intervalle interquartile, <em>IQR<\/em>) mat\u00e9rialise directement la variabilit\u00e9 du c\u0153ur des observations. Une disparit\u00e9 manifeste de hauteur entre les bo\u00eetes des diff\u00e9rentes modalit\u00e9s fournit une alerte visuelle imm\u00e9diate quant \u00e0 une probable h\u00e9t\u00e9rosc\u00e9dasticit\u00e9. De m\u00eame, l&rsquo;\u00e9tendue globale des moustaches illustre l&rsquo;amplitude de dispersion des extr\u00eames.<\/p>\n<p>Sur le plan quantitatif, les m\u00e9thodologistes s&rsquo;appuient fr\u00e9quemment sur la <strong>r\u00e8gle empirique du ratio de variance critique<\/strong> :<\/p>\n<p><em>Ratio = s<sub>max<\/sub><sup>2<\/sup> \/ s<sub>min<\/sub><sup>2<\/sup><\/em><\/p>\n<p>O\u00f9 <em>s<sub>max<\/sub><sup>2<\/sup><\/em> d\u00e9signe la variance empirique la plus \u00e9lev\u00e9e parmi les groupes analys\u00e9s, et <em>s<sub>min<\/sub><sup>2<\/sup><\/em> repr\u00e9sente la variance la plus faible. Selon les recommandations classiques avanc\u00e9es par Tabachnick et Fidell, un ratio inf\u00e9rieur \u00e0 1,5 ou 2 sugg\u00e8re une homosc\u00e9dasticit\u00e9 acceptable autorisant le recours \u00e0 l&rsquo;ANOVA de Fisher sans distorsion notable. En revanche, lorsque le ratio exc\u00e8de le seuil de 3 ou 4 (et particuli\u00e8rement si les tailles d&rsquo;\u00e9chantillons sont asym\u00e9triques), l&rsquo;ANOVA classique doit \u00eatre imp\u00e9rativement abandonn\u00e9e au profit du test de Welch.<\/p>\n<h3>3.3 Limites des transformations de variables traditionnelles<\/h3>\n<p>Pendant plusieurs d\u00e9cennies, la r\u00e9ponse p\u00e9dagogique conventionnelle face \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 a consist\u00e9 \u00e0 recommander des transformations non lin\u00e9aires des variables d\u00e9pendantes : passage au logarithme n\u00e9p\u00e9rien <em>log(Y)<\/em>, calcul de la racine carr\u00e9e <em>&radic;Y<\/em>, transformation inverse <em>1\/Y<\/em>, ou encore projection arcsinus pour les proportions binomiques. Si ces transformations math\u00e9matiques parviennent occasionnellement \u00e0 stabiliser les variances, la recherche statistique contemporaine d\u00e9conseille fermement leur utilisation syst\u00e9matique pour plusieurs motifs \u00e9pist\u00e9mologiques et pratiques :<\/p>\n<ul>\n<li><strong>Alt\u00e9ration de l&rsquo;interpr\u00e9tabilit\u00e9 psychom\u00e9trique :<\/strong> La transformation modifie l&rsquo;\u00e9chelle de mesure et la m\u00e9trique intrins\u00e8que de la variable psychologique. Expliquer qu&rsquo;une m\u00e9thode p\u00e9dagogique am\u00e9liore le temps de r\u00e9ponse moyen de 150 millisecondes poss\u00e8de une pertinence clinique imm\u00e9diate ; affirmer qu&rsquo;elle r\u00e9duit l&rsquo;inverse du logarithme du temps de r\u00e9ponse de 0,04 unit\u00e9s abstraites prive le praticien de toute prise intuitive sur la r\u00e9alit\u00e9 empirique.<\/li>\n<li><strong>\u00c9chec fr\u00e9quent sur les distributions r\u00e9elles :<\/strong> Les transformations unilat\u00e9rales ne parviennent que rarement \u00e0 normaliser simultan\u00e9ment l&rsquo;asym\u00e9trie et \u00e0 harmoniser les variances \u00e0 travers tous les groupes. Une transformation logarithmique efficace pour r\u00e9duire la variance d&rsquo;un groupe fortement dispers\u00e9 peut comprimer de fa\u00e7on excessive la variance d&rsquo;un groupe d\u00e9j\u00e0 homog\u00e8ne, d\u00e9pla\u00e7ant le probl\u00e8me sans le r\u00e9soudre.<\/li>\n<li><strong>Changement de la nature des hypoth\u00e8ses :<\/strong> Tester une diff\u00e9rence de moyennes sur des donn\u00e9es transform\u00e9es logarithmiquement \u00e9quivaut \u00e0 tester un rapport de moyennes g\u00e9om\u00e9triques sur l&rsquo;\u00e9chelle originale, ce qui ne correspond pas n\u00e9cessairement \u00e0 l&rsquo;hypoth\u00e8se conceptuelle initialement formul\u00e9e par le chercheur.<\/li>\n<\/ul>\n<p>Par cons\u00e9quent, le consensus scientifique actuel privil\u00e9gie sans ambigu\u00eft\u00e9 <strong>l&rsquo;adaptation du mod\u00e8le statistique aux donn\u00e9es r\u00e9elles<\/strong> (via des estimateurs robustes comme ceux de Welch) plut\u00f4t que la d\u00e9formation artificielle des donn\u00e9es pour satisfaire les contraintes d&rsquo;un mod\u00e8le inadapt\u00e9.<\/p>\n<h2>4. Cadre math\u00e9matique et algorithmique du test de Welch<\/h2>\n<p>L&rsquo;\u00e9l\u00e9gance du test de Welch r\u00e9side dans son architecture probabiliste : il s&rsquo;agit d&rsquo;une reformulation analytique rigoureuse qui substitue \u00e0 l&rsquo;estimation globale de l&rsquo;erreur une pond\u00e9ration diff\u00e9rentielle bas\u00e9e sur la pr\u00e9cision empirique de chaque \u00e9chantillon.<\/p>\n<h3>4.1 Pond\u00e9ration par l&rsquo;inverse de la variance<\/h3>\n<p>Soit un dispositif exp\u00e9rimental comportant <em>k<\/em> groupes ind\u00e9pendants. Pour chaque groupe <em>i<\/em> (avec <em>i &isin; {1, &#8230;, k}<\/em>), on observe un effectif <em>n<sub>i<\/sub><\/em>, une moyenne empirique <em>X\u0304<sub>i<\/sub><\/em> et une variance non biais\u00e9e <em>s<sub>i<\/sub><sup>2<\/sup><\/em>. La variance de la moyenne d&rsquo;\u00e9chantillonnage pour le groupe <em>i<\/em> s&rsquo;\u00e9crit classiquement :<\/p>\n<p><em>Var(X\u0304<sub>i<\/sub>) = s<sub>i<\/sub><sup>2<\/sup> \/ n<sub>i<\/sub><\/em><\/p>\n<p>L&rsquo;innovation fondamentale de Welch consiste \u00e0 d\u00e9finir un coefficient de pond\u00e9ration <em>w<sub>i<\/sub><\/em> pour chaque modalit\u00e9, \u00e9gal \u00e0 l&rsquo;inverse de cette variance d&rsquo;\u00e9chantillonnage. Ce poids refl\u00e8te le degr\u00e9 de pr\u00e9cision informationnelle apport\u00e9 par chaque sous-\u00e9chantillon :<\/p>\n<p><em>w<sub>i<\/sub> = n<sub>i<\/sub> \/ s<sub>i<\/sub><sup>2<\/sup><\/em><\/p>\n<p>Plus la variance d&rsquo;un groupe est massive et son effectif r\u00e9duit, plus son poids <em>w<sub>i<\/sub><\/em> devient n\u00e9gligeable dans le calcul inf\u00e9rentiel. \u00c0 partir de ces poids individuels, Welch d\u00e9finit une <strong>moyenne g\u00e9n\u00e9rale pond\u00e9r\u00e9e<\/strong> (d\u00e9sign\u00e9e ici par <em>X\u0304&rsquo;<\/em>), qui remplace la moyenne globale brute de l&rsquo;ANOVA classique :<\/p>\n<p><em>X\u0304&rsquo; = [&sum;<sub>i=1<\/sub><sup>k<\/sup> w<sub>i<\/sub> X\u0304<sub>i<\/sub>] \/ [&sum;<sub>i=1<\/sub><sup>k<\/sup> w<sub>i<\/sub>]<\/em><\/p>\n<p>Notons <em>W<\/em> la somme totale des poids :<\/p>\n<p><em>W = &sum;<sub>i=1<\/sub><sup>k<\/sup> w<sub>i<\/sub><\/em><\/p>\n<p>Cette moyenne d&rsquo;ensemble pond\u00e9r\u00e9e attire le barycentre des donn\u00e9es vers les groupes b\u00e9n\u00e9ficiant de la plus forte fiabilit\u00e9 m\u00e9trologique, emp\u00eachant un petit groupe \u00e0 variance colossale de d\u00e9former arbitrairement la tendance centrale globale.<\/p>\n<h3>4.2 Formulation de la statistique de test F de Welch<\/h3>\n<p>La statistique de d\u00e9cision <em>F<sub>w<\/sub><\/em> (souvent not\u00e9e simplement <em>F<\/em> de Welch) s&rsquo;articule autour d&rsquo;un num\u00e9rateur refl\u00e9tant la d\u00e9viation quadratique pond\u00e9r\u00e9e des moyennes de groupe par rapport \u00e0 la moyenne globale pond\u00e9r\u00e9e, normalis\u00e9 par un terme correctif d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 <em>&Lambda;<\/em> ins\u00e9r\u00e9 au d\u00e9nominateur :<\/p>\n<p>Le num\u00e9rateur de la statistique de Welch s&rsquo;exprime ainsi :<\/p>\n<p><em>Num\u00e9rateur = [1 \/ (k &#8211; 1)] &sum;<sub>i=1<\/sub><sup>k<\/sup> w<sub>i<\/sub> (X\u0304<sub>i<\/sub> &#8211; X\u0304&rsquo;)<sup>2<\/sup><\/em><\/p>\n<p>Sous l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9galit\u00e9 parfaite des moyennes th\u00e9oriques sous-jacentes (<em>H<sub>0<\/sub>: &mu;<sub>1<\/sub> = &mu;<sub>2<\/sub> = &#8230; = &mu;<sub>k<\/sub><\/em>), ce num\u00e9rateur converge vers une esp\u00e9rance math\u00e9matique sup\u00e9rieure \u00e0 1 en raison des fluctuations d&rsquo;\u00e9chantillonnage des variances. Pour neutraliser ce biais, Welch d\u00e9rive un terme correctif sans dimension, souvent d\u00e9sign\u00e9 par <em>&Lambda;<\/em>, d\u00e9fini comme suit :<\/p>\n<p><em>&Lambda; = [2(k &#8211; 2) \/ (k<sup>2<\/sup> &#8211; 1)] &sum;<sub>i=1<\/sub><sup>k<\/sup> [ (1 &#8211; w<sub>i<\/sub> \/ W)<sup>2<\/sup> \/ (n<sub>i<\/sub> &#8211; 1) ]<\/em><\/p>\n<p>La statistique finale du test de Welch s&rsquo;\u00e9tablit alors formellement par le quotient suivant :<\/p>\n<p><em>F<sub>w<\/sub> = [ [1 \/ (k &#8211; 1)] &sum;<sub>i=1<\/sub><sup>k<\/sup> w<sub>i<\/sub> (X\u0304<sub>i<\/sub> &#8211; X\u0304&rsquo;)<sup>2<\/sup> ] \/ [ 1 + &Lambda; ]<\/em><\/p>\n<p>Dans cette formulation, si toutes les variances empiriques <em>s<sub>i<\/sub><sup>2<\/sup><\/em> sont rigoureusement identiques et les effectifs <em>n<sub>i<\/sub><\/em> \u00e9gaux, les poids <em>w<sub>i<\/sub><\/em> deviennent \u00e9quivalents et la statistique <em>F<sub>w<\/sub><\/em> pr\u00e9sente un comportement asymptotique identique \u00e0 la statistique <em>F<\/em> de Fisher standard issue de l&rsquo;estimation des moindres carr\u00e9s ordinaires.<\/p>\n<h3>4.3 Ajustement de Satterthwaite pour les degr\u00e9s de libert\u00e9<\/h3>\n<p>L&rsquo;un des apports les plus d\u00e9terminants de Bernard L. Welch r\u00e9side dans l&rsquo;ajustement non entier des degr\u00e9s de libert\u00e9 du d\u00e9nominateur. Dans l&rsquo;ANOVA classique de Fisher, les degr\u00e9s de libert\u00e9 sont des nombres entiers immuables dict\u00e9s par la combinatoire du devis : <em>df<sub>1<\/sub> = k &#8211; 1<\/em> pour l&rsquo;effet inter-groupes, et <em>df<sub>2<\/sub> = N &#8211; k<\/em> pour l&rsquo;erreur r\u00e9siduelle intra-groupes.<\/p>\n<p>Puisque la statistique de Welch <em>F<sub>w<\/sub><\/em> n&rsquo;est pas distribu\u00e9e exactement comme un ratio standard de deux variables du chi carr\u00e9 ind\u00e9pendantes sous h\u00e9t\u00e9rosc\u00e9dasticit\u00e9, Welch applique une approximation de type Satterthwaite pour contraindre la distribution empirique de <em>F<sub>w<\/sub><\/em> \u00e0 suivre d&rsquo;aussi pr\u00e8s que possible une loi th\u00e9orique de Fisher-Snedecor <em>F(df<sub>1<\/sub>, df<sub>2(Welch)<\/sub>)<\/em>. Les degr\u00e9s de libert\u00e9 du num\u00e9rateur demeurent entiers :<\/p>\n<p><em>df<sub>1<\/sub> = k &#8211; 1<\/em><\/p>\n<p>En revanche, les degr\u00e9s de libert\u00e9 associ\u00e9s au d\u00e9nominateur (l&rsquo;erreur d&rsquo;\u00e9chantillonnage ajust\u00e9e), not\u00e9s <em>df<sub>2(Welch)<\/sub><\/em> ou <em>df<sub>w<\/sub><\/em>, deviennent une fonction continue des variances observ\u00e9es et des tailles de groupes :<\/p>\n<p><em>df<sub>2(Welch)<\/sub> = (k<sup>2<\/sup> &#8211; 1) \/ [ 3 &sum;<sub>i=1<\/sub><sup>k<\/sup> [ (1 &#8211; w<sub>i<\/sub> \/ W)<sup>2<\/sup> \/ (n<sub>i<\/sub> &#8211; 1) ] ]<\/em><\/p>\n<p>Ce calcul d\u00e9bouche presque invariablement sur une <strong>valeur non enti\u00e8re (fractionnaire ou d\u00e9cimale)<\/strong>, par exemple <em>df<sub>2<\/sub> = 14,38<\/em>. Cette propri\u00e9t\u00e9 d\u00e9concerte parfois les analystes novices form\u00e9s exclusivement \u00e0 la lecture de tables statistiques papier historiques. Sur le plan probabiliste, cette valeur d\u00e9cimale calibre exactement la surface sous la courbe de la fonction de densit\u00e9 de Fisher, modulant avec une pr\u00e9cision chirurgicale la valeur critique au-del\u00e0 de laquelle la d\u00e9cision de rejet de l&rsquo;hypoth\u00e8se nulle intervient.<\/p>\n<h2>5. Configuration de l&rsquo;environnement Python et pr\u00e9paration des biblioth\u00e8ques<\/h2>\n<p>Le langage Python s&rsquo;est impos\u00e9 comme une plateforme de r\u00e9f\u00e9rence pour le traitement de donn\u00e9es et l&rsquo;inf\u00e9rence statistique computationnelle en psychologie exp\u00e9rimentale et en neurosciences comportementales. Sa modularit\u00e9 permet d&rsquo;articuler la manipulation tabulaire de haut niveau, le calcul scientifique matriciel et la mod\u00e9lisation statistique avanc\u00e9e au sein d&rsquo;environnements d&rsquo;analyse enti\u00e8rement reproductibles.<\/p>\n<h3>5.1 \u00c9cosyst\u00e8me scientifique Python pour la statistique inf\u00e9rentielle<\/h3>\n<p>Pour ex\u00e9cuter rigoureusement un protocole d&rsquo;ANOVA de Welch, quatre biblioth\u00e8ques compl\u00e9mentaires structurent l&rsquo;environnement scientifique :<\/p>\n<ul>\n<li><strong><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/numpy.org\/\">NumPy<\/a> :<\/strong> Fournit les structures vectorielles et matricielles multidimensionnelles (tableaux <code>ndarray<\/code>), autorisant des calculs alg\u00e9briques optimis\u00e9s en code pr\u00e9compil\u00e9 pour impl\u00e9menter manuellement des \u00e9quations statistiques complexes sans recourir \u00e0 des boucles interpr\u00e9t\u00e9es lentes.<\/li>\n<li><strong><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pandas.pydata.org\/\">Pandas<\/a> :<\/strong> Constitue la r\u00e9f\u00e9rence absolue pour l&rsquo;ingestion, le nettoyage, le pivotement et l&rsquo;agr\u00e9gation structur\u00e9e des jeux de donn\u00e9es comportementales sous forme de <code>DataFrame<\/code>. Elle g\u00e8re avec transparence les structures de donn\u00e9es en format \u00ab long \u00bb (une ligne par observation, adapt\u00e9e aux mod\u00e8les statistiques) ou \u00ab large \u00bb (une colonne par modalit\u00e9 exp\u00e9rimentale).<\/li>\n<li><strong><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/scipy.org\/\">SciPy<\/a> (notamment le sous-module <code>scipy.stats<\/code>) :<\/strong> Offre un \u00e9ventail exhaustif de lois de probabilit\u00e9 continues et discr\u00e8tes, de tests d&rsquo;ad\u00e9quation distributionnelle (Shapiro-Wilk) et d&rsquo;\u00e9valuations d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances (Bartlett, Levene).<\/li>\n<li><strong>Pingouin :<\/strong> Con\u00e7ue sp\u00e9cifiquement pour la recherche en psychologie et biom\u00e9trie par Rapha\u00ebl Vallat, cette biblioth\u00e8que comble l&rsquo;\u00e9cart historique existant entre la convivialit\u00e9 de R et la puissance g\u00e9n\u00e9rale de Python. Pingouin impl\u00e9mente nativement l&rsquo;ANOVA de Welch, les tests post-hoc de Games-Howell et le calcul automatis\u00e9 d&rsquo;estimateurs robustes de tailles d&rsquo;effet, restituant syst\u00e9matiquement des tableaux de r\u00e9sultats parfaitement document\u00e9s.<\/li>\n<\/ul>\n<h3>5.2 Installation et importation des modules requis<\/h3>\n<p>Pour garantir la reproductibilit\u00e9 absolue des calculs pr\u00e9sent\u00e9s dans ce tutoriel, il est pr\u00e9conis\u00e9 d&rsquo;op\u00e9rer au sein d&rsquo;un environnement virtuel d\u00e9di\u00e9 (via <code>venv<\/code> ou <code>conda<\/code>). L&rsquo;installation des d\u00e9pendances s&rsquo;effectue au moyen de l&rsquo;outil de gestion de paquets <code>pip<\/code> \u00e0 travers une unique instruction d&rsquo;installation :<\/p>\n<p><code>pip install numpy pandas scipy pingouin seaborn matplotlib<\/code><\/p>\n<p>Une fois les paquets install\u00e9s, la s\u00e9quence d&rsquo;importation et la configuration globale des param\u00e8tres d&rsquo;affichage scientifique s&rsquo;articulent selon le script standardis\u00e9 suivant :<\/p>\n<p>Dans ce bloc d&rsquo;initialisation, nous importons les modules essentiels sous leurs alias canoniques conventionnels (<code>np<\/code>, <code>pd<\/code>, <code>sns<\/code>, <code>plt<\/code>), puis nous importons le module de statistiques de SciPy et l&rsquo;ensemble de la biblioth\u00e8que Pingouin. Nous for\u00e7ons l&rsquo;affichage de Pandas \u00e0 afficher jusqu&rsquo;\u00e0 quatre d\u00e9cimales apr\u00e8s la virgule afin d&rsquo;observer avec exactitude les valeurs d\u00e9cimales des degr\u00e9s de libert\u00e9 de Welch ainsi que les <em>p<\/em>-valeurs infinit\u00e9simales :<\/p>\n<p><code>import numpy as np<\/code><br \/>\n<code>import pandas as pd<\/code><br \/>\n<code>import scipy.stats as stats<\/code><br \/>\n<code>import pingouin as pg<\/code><br \/>\n<code>import matplotlib.pyplot as plt<\/code><br \/>\n<code>import seaborn as sns<\/code><br \/>\n<code>pd.set_option('display.float_format', lambda x: '%.4f' % x)<\/code><br \/>\n<code>sns.set_theme(style=\"whitegrid\", palette=\"deep\")<\/code><\/p>\n<p>Ce param\u00e9trage garantit une esth\u00e9tique graphique \u00e9pur\u00e9e conforme aux crit\u00e8res graphiques impos\u00e9s par les revues de recherche internationales.<\/p>\n<h2>6. \u00c9tape 1 : Structuration et exploration descriptive des donn\u00e9es empiriques<\/h2>\n<p>Toute analyse inf\u00e9rentielle doit s&rsquo;ancrer dans un contexte exp\u00e9rimental r\u00e9aliste. Pour illustrer concr\u00e8tement la d\u00e9marche, nous concevons un protocole exp\u00e9rimental repr\u00e9sentatif des recherches en psychologie de l&rsquo;\u00e9ducation et en psychologie cognitive portant sur l&rsquo;optimisation des strat\u00e9gies d&rsquo;encodage mn\u00e9sique.<\/p>\n<h3>6.1 D\u00e9finition du devis exp\u00e9rimental sur les techniques d&rsquo;\u00e9tude<\/h3>\n<p>Un chercheur en sciences cognitives s&rsquo;int\u00e9resse \u00e0 l&rsquo;impact de trois m\u00e9thodes distinctes de r\u00e9vision p\u00e9dagogique sur la performance finale \u00e0 un examen standardis\u00e9 de compr\u00e9hension conceptuelle. Trente \u00e9tudiants universitaires (<em>N = 30<\/em>) sont assign\u00e9s de mani\u00e8re al\u00e9atoire \u00e0 l&rsquo;une des trois conditions suivantes (avec des effectifs initiaux \u00e9quilibr\u00e9s, <em>n<sub>1<\/sub> = n<sub>2<\/sub> = n<sub>3<\/sub> = 10<\/em>) :<\/p>\n<ul>\n<li><strong>Condition A (R\u00e9cup\u00e9ration active en m\u00e9moire \/ <em>Active Recall<\/em>) :<\/strong> Les participants lisent le contenu acad\u00e9mique une seule fois, puis s&rsquo;exercent intensivement \u00e0 g\u00e9n\u00e9rer les r\u00e9ponses et d\u00e9finitions de m\u00e9moire sans support documentaire. Cette m\u00e9thode est connue pour susciter un b\u00e9n\u00e9fice mn\u00e9sique substantiel mais peut induire de fortes variations individuelles selon l&rsquo;aisance initiale des sujets.<\/li>\n<li><strong>Condition B (Cartographie conceptuelle \/ <em>Mind Mapping<\/em>) :<\/strong> Les \u00e9tudiants structurent les notions th\u00e9oriques sous la forme d&rsquo;un r\u00e9seau hi\u00e9rarchique reliant les concepts cl\u00e9s par des liens logiques explicites. Cette approche fournit un cadre d&rsquo;\u00e9tayage visuel interm\u00e9diaire.<\/li>\n<li><strong>Condition C (Relecture passive \/ <em>Re-reading<\/em>) :<\/strong> Les \u00e9tudiants relisent de fa\u00e7on lin\u00e9aire et r\u00e9p\u00e9t\u00e9e le document source sans manipulation active de l&rsquo;information (condition de contr\u00f4le repr\u00e9sentative des habitudes spontan\u00e9es de travail).<\/li>\n<\/ul>\n<p>\u00c0 l&rsquo;issue d&rsquo;une phase de r\u00e9tention de quarante-huit heures, tous les \u00e9tudiants passent le m\u00eame examen not\u00e9 sur une \u00e9chelle continue de 0 \u00e0 100 points. Le facteur inter-sujets ind\u00e9pendant (VI) correspond donc \u00e0 la \u00ab Technique d&rsquo;\u00e9tude \u00bb (\u00e0 trois modalit\u00e9s cat\u00e9gorielles : A, B, C), et la variable d\u00e9pendante mesur\u00e9e (VD) correspond au \u00ab Score \u00e0 l&rsquo;examen \u00bb.<\/p>\n<h3>6.2 Cr\u00e9ation et vectorisation des donn\u00e9es en Python<\/h3>\n<p>Pour mod\u00e9liser l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances typique des contextes empiriques r\u00e9els, les scores sont configur\u00e9s de mani\u00e8re \u00e0 refl\u00e9ter des moyennes diff\u00e9rentes tout en incorporant des niveaux de variabilit\u00e9 intra-groupe intentionnellement divergents :<\/p>\n<p>La condition A pr\u00e9sente des scores \u00e9lev\u00e9s mais une forte dispersion (certains \u00e9tudiants ma\u00eetrisant remarquablement la r\u00e9cup\u00e9ration active, d&rsquo;autres se heurtant \u00e0 l&rsquo;absence de support documentaire). La condition B pr\u00e9sente une moyenne interm\u00e9diaire et une variance mod\u00e9r\u00e9e. La condition C affiche une moyenne faible avec une variance extr\u00eamement r\u00e9duite (effet d&rsquo;homog\u00e9n\u00e9isation par le bas li\u00e9 \u00e0 l&rsquo;inefficacit\u00e9 constante de la simple relecture passive).<\/p>\n<p>D\u00e9clarons ces s\u00e9ries de donn\u00e9es sous forme de listes num\u00e9riques, puis convertissons-les au sein d&rsquo;un <code>DataFrame<\/code> Pandas au format \u00ab long \u00bb standard, requis par la plupart des algorithmes statistiques contemporains :<\/p>\n<p><code>scores_A = [88, 92, 65, 96, 70, 85, 99, 62, 78, 95]<\/code><br \/>\n<code>scores_B = [75, 78, 72, 80, 77, 83, 71, 79, 74, 81]<\/code><br \/>\n<code>scores_C = [60, 62, 59, 61, 63, 58, 64, 60, 61, 62]<\/code><\/p>\n<p>Pour construire le tableau de donn\u00e9es unifi\u00e9, nous combinons les scores et les \u00e9tiquettes de groupes associ\u00e9es :<\/p>\n<p><code>groupe = ['Technique A'] * len(scores_A) + ['Technique B'] * len(scores_B) + ['Technique C'] * len(scores_C)<\/code><br \/>\n<code>scores = scores_A + scores_B + scores_C<\/code><br \/>\n<code>df_etude = pd.DataFrame({'Groupe': groupe, 'Score': scores})<\/code><\/p>\n<p>Une inspection structurelle rapide via <code>df_etude.info()<\/code> confirme la pr\u00e9sence de 30 lignes compl\u00e8tes, exemptes de donn\u00e9es manquantes, typ\u00e9es de fa\u00e7on ad\u00e9quate sous forme d&rsquo;une cha\u00eene de caract\u00e8res cat\u00e9gorielle pour la colonne &lsquo;Groupe&rsquo; et d&rsquo;un entier num\u00e9rique <code>int64<\/code> pour la variable &lsquo;Score&rsquo;.<\/p>\n<h3>6.3 Calcul exhaustif des statistiques descriptives<\/h3>\n<p>L&rsquo;exploration descriptive exige une analyse granulaire des m\u00e9triques de position centrale et d&rsquo;\u00e9talement pour chaque niveau du facteur ind\u00e9pendant. Pandas permet d&rsquo;agr\u00e9ger ces donn\u00e9es au moyen de la m\u00e9thode <code>groupby()<\/code> coupl\u00e9e \u00e0 l&rsquo;op\u00e9rateur <code>agg()<\/code> :<\/p>\n<p><code>stats_desc = df_etude.groupby('Groupe')['Score'].agg(<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;Effectif='count',<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;Moyenne='mean',<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;Ecart_Type='std',<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;Variance=lambda x: x.var(ddof=1),<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;Erreur_Standard=lambda x: stats.sem(x),<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;IC95_Inf=lambda x: stats.t.interval(0.95, len(x)-1, loc=np.mean(x), scale=stats.sem(x))[0],<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;IC95_Sup=lambda x: stats.t.interval(0.95, len(x)-1, loc=np.mean(x), scale=stats.sem(x))[1]<\/code><br \/>\n<code>).reset_index()<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de cette instruction met en exergue des contrastes frappants au sein de notre \u00e9chantillon :<\/p>\n<ul>\n<li><strong>Technique A :<\/strong> Moyenne = 83,0000 ; \u00c9cart-type = 13,8484 ; Variance = 191,7778 ; IC \u00e0 95 % = [73,0915 ; 92,9085].<\/li>\n<li><strong>Technique B :<\/strong> Moyenne = 77,0000 ; \u00c9cart-type = 4,0000 ; Variance = 16,0000 ; IC \u00e0 95 % = [74,1383 ; 79,8617].<\/li>\n<li><strong>Technique C :<\/strong> Moyenne = 61,0000 ; \u00c9cart-type = 1,8257 ; Variance = 3,3333 ; IC \u00e0 95 % = [59,6938 ; 62,3062].<\/li>\n<\/ul>\n<p>L&rsquo;\u00e9cart empirique entre les variances appara\u00eet vertigineux. La variance observ\u00e9e au sein de la Technique A (191,78) d\u00e9passe de plus de cinquante-sept fois celle de la Technique C (3,33). Le ratio critique de variance atteint ici la valeur de :<\/p>\n<p><em>Ratio = s<sub>max<\/sub><sup>2<\/sup> \/ s<sub>min<\/sub><sup>2<\/sup> = 191,7778 \/ 3,3333 &approx; 57,53<\/em><\/p>\n<p>Cette valeur exc\u00e8de largement le seuil empirique de tol\u00e9rance fix\u00e9 conventionnellement \u00e0 3 ou 4. L&rsquo;utilisation d&rsquo;une ANOVA classique de Fisher sur un tel jeu de donn\u00e9es constituerait une erreur m\u00e9thodologique patente, susceptible de fausser lourdement l&rsquo;estimation de l&rsquo;erreur d&rsquo;\u00e9chantillonnage.<\/p>\n<h2>7. \u00c9tape 2 : V\u00e9rification rigoureuse des postulats d&rsquo;homosc\u00e9dasticit\u00e9 et de normalit\u00e9<\/h2>\n<p>L&rsquo;exploration visuelle et le calcul des ratios descriptifs fournissent des pr\u00e9somptions fond\u00e9es, mais la d\u00e9marche inf\u00e9rentielle impose une validation formelle au moyen de tests d&rsquo;hypoth\u00e8ses d\u00e9di\u00e9s. L&rsquo;\u00e9valuation de l&rsquo;homosc\u00e9dasticit\u00e9 et de la normalit\u00e9 intra-groupe conditionne directement le choix du mod\u00e8le lin\u00e9aire final.<\/p>\n<h3>7.1 Application du test d&rsquo;\u00e9galit\u00e9 des variances de Bartlett<\/h3>\n<p>Le test de Bartlett (1937) \u00e9value l&rsquo;hypoth\u00e8se nulle selon laquelle l&rsquo;ensemble des <em>k<\/em> variances de population sont parfaitement identiques :<\/p>\n<p><em>H<sub>0<\/sub>: &sigma;<sub>1<\/sub><sup>2<\/sup> = &sigma;<sub>2<\/sub><sup>2<\/sup> = &#8230; = &sigma;<sub>k<\/sub><sup>2<\/sup><\/em><\/p>\n<p>Cette proc\u00e9dure compare la moyenne arithm\u00e9tique pond\u00e9r\u00e9e des variances d&rsquo;\u00e9chantillons \u00e0 leur moyenne g\u00e9om\u00e9trique pond\u00e9r\u00e9e. La statistique de Bartlett suit asymptotiquement une loi du chi carr\u00e9 <em>&chi;<sup>2<\/sup><\/em> \u00e0 <em>k &#8211; 1<\/em> degr\u00e9s de libert\u00e9, pourvu que la distribution des populations d&rsquo;origine soit rigoureusement normale. L&rsquo;impl\u00e9mentation dans l&rsquo;\u00e9cosyst\u00e8me SciPy s&rsquo;ex\u00e9cute \u00e0 l&rsquo;aide de la fonction <code>stats.bartlett()<\/code> :<\/p>\n<p><code>stat_bartlett, p_bartlett = stats.bartlett(scores_A, scores_B, scores_C)<\/code><br \/>\n<code>print(f\"Test de Bartlett : Statistique = {stat_bartlett:.4f}, p-valeur = {p_bartlett:.6e}\")<\/code><\/p>\n<p>L&rsquo;analyse fournit les r\u00e9sultats suivants :<\/p>\n<p><em>Statistique &chi;<sup>2<\/sup> = 25,4839 ; p = 2,9258 &times; 10<sup>-6<\/sup><\/em><\/p>\n<p>La <em>p<\/em>-valeur se situant infiniment en de\u00e7\u00e0 du seuil de significativit\u00e9 conventionnel fix\u00e9 \u00e0 <em>&alpha; = 0,05<\/em>, nous rejetons sans ambigu\u00eft\u00e9 l&rsquo;hypoth\u00e8se nulle d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances. Ce r\u00e9sultat confirme formellement la pr\u00e9sence d&rsquo;une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 structurelle majeure entre les trois modalit\u00e9s de techniques d&rsquo;apprentissage.<\/p>\n<h3>7.2 Tests compl\u00e9mentaires d&rsquo;homog\u00e9n\u00e9it\u00e9 : Levene et Brown-Forsythe<\/h3>\n<p>Si le test de Bartlett se distingue par une puissance statistique remarquable sous distribution normale stricte, il s&rsquo;av\u00e8re particuli\u00e8rement hypersensible aux l\u00e9g\u00e8res entorses de normalit\u00e9 (il a tendance \u00e0 rejeter l&rsquo;homosc\u00e9dasticit\u00e9 de fa\u00e7on intempestive en pr\u00e9sence d&rsquo;un aplatissement non nul, agissant parfois davantage comme un test de normalit\u00e9 d\u00e9guis\u00e9 que comme un test d&rsquo;homog\u00e9n\u00e9it\u00e9 pure). Pour s\u00e9curiser notre diagnostic, il convient de mobiliser le test de <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Levene%27s_test\">Levene<\/a> et sa variante robuste de Brown-Forsythe.<\/p>\n<p>Le test de Levene calcule une analyse de variance sur les d\u00e9viations absolues des scores individuels par rapport \u00e0 une mesure de tendance centrale du groupe. Dans sa d\u00e9clinaison originale, le test utilise la moyenne arithm\u00e9tique. Brown et Forsythe (1974) ont optimis\u00e9 cette proc\u00e9dure en substituant la m\u00e9diane \u00e0 la moyenne, ce qui conf\u00e8re au test une robustesse sup\u00e9rieure face aux distributions asym\u00e9triques. Nous ex\u00e9cutons le test de Levene centr\u00e9 sur la m\u00e9diane via la fonction <code>stats.levene()<\/code> en sp\u00e9cifiant l&rsquo;argument <code>center='median'<\/code> :<\/p>\n<p><code>stat_levene, p_levene = stats.levene(scores_A, scores_B, scores_C, center='median')<\/code><br \/>\n<code>print(f\"Test de Levene (Brown-Forsythe) : W = {stat_levene:.4f}, p-valeur = {p_levene:.6f}\")<\/code><\/p>\n<p>Le diagnostic peut \u00eatre corrobor\u00e9 instantan\u00e9ment avec le package Pingouin via la commande <code>pg.homoscedasticity(df_etude, dv='Score', group='Groupe', method='levene')<\/code>. Les sorties concordantes indiquent :<\/p>\n<p><em>Statistique W = 10,8745 ; p = 0,000331<\/em><\/p>\n<p>La <em>p<\/em>-valeur (<em>p &lt; 0,001<\/em>) demeure largement inf\u00e9rieure \u00e0 0,05. M\u00eame sous l&rsquo;estimation robuste de Brown-Forsythe immunis\u00e9e contre les potentielles d\u00e9viations d&rsquo;asym\u00e9trie, l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des variances est corrobor\u00e9e \u00e0 un niveau de confiance exceptionnel. L&rsquo;obligation statistique d&rsquo;abandonner l&rsquo;ANOVA de Fisher standard est ainsi d\u00e9finitivement scell\u00e9e.<\/p>\n<h3>7.3 \u00c9valuation de la normalit\u00e9 intra-groupe<\/h3>\n<p>Le test de Welch tol\u00e8re des variances in\u00e9gales, mais son assise th\u00e9orique postule n\u00e9anmoins que les distributions sous-jacentes au sein de chaque modalit\u00e9 ne violent pas grossi\u00e8rement la normalit\u00e9, particuli\u00e8rement lorsque la taille d&rsquo;\u00e9chantillon par cohorte est modeste (<em>n = 10<\/em>). Nous devons donc soumettre chaque sous-groupe au test de normalit\u00e9 de <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Shapiro%E2%80%93Wilk_test\">Shapiro-Wilk<\/a>, tout en quantifiant leurs coefficients d&rsquo;asym\u00e9trie (<em>skewness<\/em>) et d&rsquo;aplatissement (<em>kurtosis<\/em>) :<\/p>\n<p><code>norm_results = []<\/code><br \/>\n<code>for nom_grp, donnees in df_etude.groupby('Groupe')['Score']:<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;stat_sw, p_sw = stats.shapiro(donnees)<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;skew = stats.skew(donnees, bias=False)<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;kurt = stats.kurtosis(donnees, bias=False)<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;norm_results.append({'Groupe': nom_grp, 'W_Shapiro': stat_sw, 'p_valeur': p_sw, 'Skewness': skew, 'Kurtosis': kurt})<\/code><br \/>\n<code>df_norm = pd.DataFrame(norm_results)<\/code><\/p>\n<p>L&rsquo;analyse des sorties num\u00e9riques r\u00e9v\u00e8le les param\u00e8tres distributionnels suivants :<\/p>\n<ul>\n<li><strong>Technique A :<\/strong> <em>W<\/em> = 0,9371 ; <em>p<\/em> = 0,5221 ; Asym\u00e9trie = -0,4812 ; Kurtosis = -1,0625.<\/li>\n<li><strong>Technique B :<\/strong> <em>W<\/em> = 0,9673 ; <em>p<\/em> = 0,8654 ; Asym\u00e9trie = -0,0325 ; Kurtosis = -0,8741.<\/li>\n<li><strong>Technique C :<\/strong> <em>W<\/em> = 0,9487 ; <em>p<\/em> = 0,6541 ; Asym\u00e9trie = -0,1218 ; Kurtosis = -0,9812.<\/li>\n<\/ul>\n<p>Pour chacune des trois conditions, la statistique de Shapiro-Wilk affiche une <em>p<\/em>-valeur substantiellement sup\u00e9rieure au seuil critique de 0,05. Nous \u00e9chouons donc \u00e0 rejeter l&rsquo;hypoth\u00e8se nulle de normalit\u00e9 intra-groupe. Par ailleurs, les indices d&rsquo;asym\u00e9trie et de kurtosis s&rsquo;inscrivent confortablement dans l&rsquo;intervalle canonique de tol\u00e9rance [-1, +1] pr\u00e9conis\u00e9 par George et Mallery (2010). Les distributions sont raisonnablement sym\u00e9triques et m\u00e9sokurtiques. Cette concordance clinique valide parfaitement la configuration requise : la normalit\u00e9 locale est v\u00e9rifi\u00e9e mais l&rsquo;homog\u00e9n\u00e9it\u00e9 des variances est viol\u00e9e. Ce cas de figure constitue l&rsquo;indication paradigmatique absolue pour le d\u00e9ploiement de l&rsquo;ANOVA de Welch.<\/p>\n<h2>8. \u00c9tape 3 : Ex\u00e9cution de l&rsquo;ANOVA de Welch avec Pingouin et SciPy<\/h2>\n<p>La confirmation des diagnostics statistiques conduit \u00e0 l&rsquo;ex\u00e9cution de l&rsquo;analyse omnibus robuste. Nous pr\u00e9sentons deux d\u00e9marches compl\u00e9mentaires : d&rsquo;une part, l&rsquo;utilisation directe de l&rsquo;impl\u00e9mentation de haut niveau offerte par Pingouin, et d&rsquo;autre part, une programmation algorithmique manuelle sous SciPy et NumPy afin d&rsquo;illustrer la m\u00e9canique interne des corrections de Welch et Satterthwaite.<\/p>\n<h3>8.1 Mise en \u0153uvre directe avec le package Pingouin<\/h3>\n<p>La biblioth\u00e8que Pingouin offre la commande la plus fluide et exhaustive pour calculer l&rsquo;ANOVA de Welch en Python. Sa fonction <code>welch_anova()<\/code> n\u00e9cessite une structure de donn\u00e9es en format long, l&rsquo;identification de la variable d\u00e9pendante (<code>dv<\/code> pour <em>dependent variable<\/em>) et la sp\u00e9cification du facteur explicatif (<code>between<\/code>) :<\/p>\n<p><code>res_welch = pg.welch_anova(data=df_etude, dv='Score', between='Groupe')<\/code><br \/>\n<code>print(res_welch.to_string(index=False))<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de cette commande g\u00e9n\u00e8re un tableau synth\u00e9tique contenant l&rsquo;int\u00e9gralit\u00e9 des m\u00e9triques inf\u00e9rentielles n\u00e9cessaires \u00e0 l&rsquo;analyste. Le tableau r\u00e9sultant se structure de la mani\u00e8re suivante :<\/p>\n<table border=\"1\" cellpadding=\"6\" cellspacing=\"0\" style=\"border-collapse:collapse; text-align:center; margin:auto; width:90%;\">\n<thead style=\"background-color:#f2f2f2;\">\n<tr>\n<th>Source<\/th>\n<th>ddof1<\/th>\n<th>ddof2<\/th>\n<th>F<\/th>\n<th>p-unc<\/th>\n<th>np2<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td>Groupe<\/td>\n<td>2<\/td>\n<td>12,7412<\/td>\n<td>21,8357<\/td>\n<td>0,000085<\/td>\n<td>0,6180<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Analysons la s\u00e9mantique de chaque colonne produite par le moteur statistique :<\/p>\n<ul>\n<li><strong>Source :<\/strong> D\u00e9signe l&rsquo;effet inter-sujets \u00e9valu\u00e9 (le facteur cat\u00e9goriel \u00e0 trois modalit\u00e9s).<\/li>\n<li><strong>ddof1 :<\/strong> Les degr\u00e9s de libert\u00e9 associ\u00e9s au num\u00e9rateur, calcul\u00e9s selon <em>k &#8211; 1 = 3 &#8211; 1 = 2<\/em>.<\/li>\n<li><strong>ddof2 :<\/strong> Les degr\u00e9s de libert\u00e9 ajust\u00e9s du d\u00e9nominateur d&rsquo;erreur selon la formule de Satterthwaite. On observe une valeur d\u00e9cimale non enti\u00e8re de <strong>12,7412<\/strong>. Alors qu&rsquo;une ANOVA classique de Fisher aurait allou\u00e9 de fa\u00e7on arbitraire <em>N &#8211; k = 30 &#8211; 3 = 27<\/em> degr\u00e9s de libert\u00e9 \u00e0 l&rsquo;erreur, la proc\u00e9dure de Welch a contract\u00e9 de plus de moiti\u00e9 cette valeur pour p\u00e9naliser l&rsquo;incertitude induite par les variances h\u00e9t\u00e9rog\u00e8nes.<\/li>\n<li><strong>F :<\/strong> La valeur empirique de la statistique de test, culminant \u00e0 <strong>21,8357<\/strong>.<\/li>\n<li><strong>p-unc :<\/strong> La <em>p<\/em>-valeur non corrig\u00e9e (<em>uncorrected<\/em>), indiquant la probabilit\u00e9 d&rsquo;obtenir une statistique <em>F<\/em> au moins aussi extr\u00eame sous l&rsquo;hypoth\u00e8se nulle. La valeur affich\u00e9e de <strong>0,000085<\/strong> (soit <em>8,5 &times; 10<sup>-5<\/sup><\/em>) confirme une significativit\u00e9 extr\u00eame.<\/li>\n<li><strong>np2 :<\/strong> L&rsquo;\u00c9ta carr\u00e9 partiel (<em>partial eta squared<\/em>), quantifiant l&rsquo;amplitude globale de l&rsquo;effet omnibus (dont nous d\u00e9taillerons les limites et ajustements au chapitre 10).<\/li>\n<\/ul>\n<h3>8.2 Impl\u00e9mentation algorithmique manuelle avec SciPy et NumPy<\/h3>\n<p>Pour dissiper tout effet d&rsquo;\u00ab opacit\u00e9 logicielle \u00bb, d\u00e9clinons pas \u00e0 pas la r\u00e9solution math\u00e9matique des \u00e9quations fondamentales de Welch d\u00e9taill\u00e9es au chapitre 4 \u00e0 travers un script Python pur :<\/p>\n<p><code>groupes = [scores_A, scores_B, scores_C]<\/code><br \/>\n<code>k = len(groupes)<\/code><br \/>\n<code>n_i = np.array([len(g) for g in groupes], dtype=float)<\/code><br \/>\n<code>m_i = np.array([np.mean(g) for g in groupes], dtype=float)<\/code><br \/>\n<code>v_i = np.array([np.var(g, ddof=1) for g in groupes], dtype=float)<\/code><\/p>\n<p><code># 1. Calcul des poids wi = ni \/ si^2<\/code><br \/>\n<code>w_i = n_i \/ v_i<\/code><br \/>\n<code>W = np.sum(w_i)<\/code><\/p>\n<p><code># 2. Moyenne g\u00e9n\u00e9rale pond\u00e9r\u00e9e<\/code><br \/>\n<code>m_prime = np.sum(w_i * m_i) \/ W<\/code><\/p>\n<p><code># 3. Terme correctif Lambda<\/code><br \/>\n<code>terme_u = (1.0 - (w_i \/ W))**2 \/ (n_i - 1.0)<\/code><br \/>\n<code>somme_u = np.sum(terme_u)<\/code><br \/>\n<code>Lambda = (2.0 * (k - 2.0) \/ (k**2 - 1.0)) * somme_u<\/code><\/p>\n<p><code># 4. Statistique F de Welch<\/code><br \/>\n<code>numerateur_F = (1.0 \/ (k - 1.0)) * np.sum(w_i * (m_i - m_prime)**2)<\/code><br \/>\n<code>F_welch_manuel = numerateur_F \/ (1.0 + Lambda)<\/code><\/p>\n<p><code># 5. Degr\u00e9s de libert\u00e9 ajust\u00e9s de Satterthwaite<\/code><br \/>\n<code>df1 = k - 1<\/code><br \/>\n<code>df2_welch_manuel = (k**2 - 1.0) \/ (3.0 * somme_u)<\/code><\/p>\n<p><code># 6. Extraction de la p-valeur via la fonction de survie (1 - CDF)<\/code><br \/>\n<code>p_val_manuelle = stats.f.sf(F_welch_manuel, df1, df2_welch_manuel)<\/code><\/p>\n<p>Affichons les variables calcul\u00e9es au terme de cet algorithme :<\/p>\n<p><em>F_welch_manuel = 21,8357 ; df1 = 2 ; df2_welch = 12,7412 ; p-valeur = 8,4951 &times; 10<sup>-5<\/sup><\/em><\/p>\n<p>La congruence est absolue entre le calcul manuel matriciel et la fonction pr\u00e9compil\u00e9e de Pingouin. Ce croisement algorithmique illustre parfaitement comment la variance excessive du groupe A (<em>s<sub>A<\/sub><sup>2<\/sup> = 191,78<\/em>) se traduit par un poids math\u00e9matique tr\u00e8s restreint (<em>w<sub>A<\/sub> = 10 \/ 191,78 &approx; 0,0521<\/em>), alors que la pr\u00e9cision du groupe C (<em>s<sub>C<\/sub><sup>2<\/sup> = 3,33<\/em>) g\u00e9n\u00e8re un poids d&rsquo;\u00e9chantillonnage pr\u00e9pond\u00e9rant (<em>w<sub>C<\/sub> = 10 \/ 3,33 &approx; 3,0000<\/em>).<\/p>\n<h3>8.3 Interpr\u00e9tation statistique et d\u00e9cision inf\u00e9rentielle<\/h3>\n<p>La confrontation de la <em>p<\/em>-valeur empirique (<em>p = 0,000085<\/em>) au seuil de significativit\u00e9 critique pr\u00e9d\u00e9fini (<em>&alpha; = 0,05<\/em>) conduit \u00e0 une d\u00e9cision statistique univoque : nous <strong>rejetons l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9galit\u00e9 globale des moyennes<\/strong> (<em>H<sub>0<\/sub><\/em>). Il existe une diff\u00e9rence hautement significative entre les scores moyens obtenus \u00e0 l&rsquo;examen selon la m\u00e9thode de travail exp\u00e9riment\u00e9e.<\/p>\n<p>Cependant, l&rsquo;ANOVA de Welch constitue, \u00e0 l&rsquo;instar de l&rsquo;ANOVA de Fisher, un <strong>test omnibus<\/strong>. Ce rejet global informe le chercheur qu&rsquo;au moins l&rsquo;une des conditions exp\u00e9rimentales diverge significativement des autres (ou qu&rsquo;une combinaison lin\u00e9aire non nulle existe entre les groupes), mais ne pr\u00e9cise nullement la localisation sp\u00e9cifique de ces disparit\u00e9s :<\/p>\n<ul>\n<li>La Technique A surpasse-t-elle statistiquement la Technique B ?<\/li>\n<li>La Technique B est-elle significativement sup\u00e9rieure \u00e0 la Technique C ?<\/li>\n<li>La Technique A se d\u00e9marque-t-elle exclusivement du groupe de contr\u00f4le C sans diff\u00e9rer de la cartographie conceptuelle B ?<\/li>\n<\/ul>\n<p>Conclure de mani\u00e8re h\u00e2tive sur les paires individuelles \u00e0 partir du simple test global constituerait une erreur de raisonnement m\u00e9thodologique. Pour cartographier avec exactitude la trajectoire des contrastes factoriels tout en contr\u00f4lant l&rsquo;inflation de l&rsquo;erreur par famille, la mise en \u0153uvre de proc\u00e9dures de comparaisons multiples robustes s&rsquo;impose imp\u00e9rativement.<\/p>\n<h2>9. \u00c9tape 4 : Analyses post-hoc et ajustement pour comparaisons multiples<\/h2>\n<p>Dans une d\u00e9marche exp\u00e9rimentale rigoureuse, les analyses post-hoc succ\u00e8dent obligatoirement \u00e0 un r\u00e9sultat omnibus significatif. Toutefois, le choix du test de comparaisons par paires doit \u00eatre en parfaite continuit\u00e9 th\u00e9orique avec les diagnostics pos\u00e9s en amont : on ne saurait mobiliser un test global robuste \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 pour ensuite appliquer une analyse post-hoc vuln\u00e9rable aux variances in\u00e9gales.<\/p>\n<h3>9.1 Pourquoi exclure le test de Tukey HSD classique ?<\/h3>\n<p>Le test de Tukey HSD (<em>Honestly Significant Difference<\/em>) repr\u00e9sente le standard p\u00e9dagogique le plus enseign\u00e9 pour comparer des paires de moyennes apr\u00e8s une ANOVA. Pourtant, son emploi dans notre contexte exp\u00e9rimental doit \u00eatre formellement proscrit.<\/p>\n<p>La m\u00e9canique probabiliste du test de Tukey s&rsquo;appuie sur la statistique de l&rsquo;amplitude <em>Studentis\u00e9e<\/em> (loi <em>q<\/em>), qui postule une estimation uniforme de la variance d&rsquo;erreur commune. L&rsquo;erreur type de la diff\u00e9rence des moyennes utilis\u00e9e dans le d\u00e9nominateur de Tukey prend la forme :<\/p>\n<p><em>SE<sub>Tukey<\/sub> = &radic; [ (MS<sub>erreur<\/sub> \/ 2) &times; (1\/n<sub>i<\/sub> + 1\/n<sub>j<\/sub>) ]<\/em><\/p>\n<p>En incorporant le terme r\u00e9siduel mutualis\u00e9 classique <em>MS<sub>erreur<\/sub><\/em>, Tukey r\u00e9introduit au niveau des comparaisons par paires le postulat d&rsquo;homog\u00e9n\u00e9it\u00e9 que nous venons express\u00e9ment de rejeter. Les simulations statistiques de Games et Howell (1976), ainsi que les synth\u00e8ses m\u00e9thodologiques contemporaines (Field, 2018), ont d\u00e9montr\u00e9 sans \u00e9quivoque que sous h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 av\u00e9r\u00e9e des variances, le test de Tukey HSD produit des taux r\u00e9els d&rsquo;erreur de type I inacceptables, pouvant exc\u00e9der le seuil nominal de plus du triple lors des contrastes impliquant des groupes \u00e0 fortes variances.<\/p>\n<h3>9.2 Mise en \u0153uvre du test post-hoc de Games-Howell<\/h3>\n<p>La r\u00e9ponse m\u00e9thodologique adapt\u00e9e consiste \u00e0 d\u00e9ployer le <strong>test de Games-Howell<\/strong>. Cette proc\u00e9dure post-hoc a \u00e9t\u00e9 formul\u00e9e sp\u00e9cifiquement pour \u00e9tendre la logique de Welch \u00e0 l&rsquo;ensemble des comparaisons par paires sous conditions d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 et d&rsquo;effectifs in\u00e9gaux.<\/p>\n<p>Games et Howell ont modifi\u00e9 la statistique de Studentized Range en calculant l&rsquo;erreur type de chaque paire (<em>i, j<\/em>) exclusivement \u00e0 partir de leurs variances empiriques individuelles non mutualis\u00e9es :<\/p>\n<p><em>SE<sub>GH(ij)<\/sub> = &radic; [ s<sub>i<\/sub><sup>2<\/sup> \/ n<sub>i<\/sub> + s<sub>j<\/sub><sup>2<\/sup> \/ n<sub>j<\/sub> ]<\/em><\/p>\n<p>Parall\u00e8lement, les degr\u00e9s de libert\u00e9 associ\u00e9s \u00e0 chaque comparaison sp\u00e9cifique font l&rsquo;objet d&rsquo;un ajustement de type Welch-Satterthwaite pairwise d\u00e9di\u00e9 :<\/p>\n<p><em>df<sub>GH(ij)<\/sub> = (s<sub>i<\/sub><sup>2<\/sup> \/ n<sub>i<\/sub> + s<sub>j<\/sub><sup>2<\/sup> \/ n<sub>j<\/sub>)<sup>2<\/sup> \/ [ (s<sub>i<\/sub><sup>2<\/sup> \/ n<sub>i<\/sub>)<sup>2<\/sup> \/ (n<sub>i<\/sub> &#8211; 1) + (s<sub>j<\/sub><sup>2<\/sup> \/ n<sub>j<\/sub>)<sup>2<\/sup> \/ (n<sub>j<\/sub> &#8211; 1) ]<\/em><\/p>\n<p>Le package Pingouin ex\u00e9cute cette cha\u00eene computationnelle avanc\u00e9e en une unique instruction via sa fonction <code>pairwise_gameshowell()<\/code> :<\/p>\n<p><code>posthoc = pg.pairwise_gameshowell(data=df_etude, dv='Score', between='Groupe')<\/code><br \/>\n<code>colonnes_visu = ['A', 'B', 'mean(A)', 'mean(B)', 'diff', 'se', 'T', 'df', 'pval', 'hedges']<\/code><br \/>\n<code>print(posthoc[colonnes_visu].to_string(index=False))<\/code><\/p>\n<p>Le tableau g\u00e9n\u00e9r\u00e9 compile les contrastes bivari\u00e9s complets :<\/p>\n<table border=\"1\" cellpadding=\"6\" cellspacing=\"0\" style=\"border-collapse:collapse; text-align:center; margin:auto; width:95%;\">\n<thead style=\"background-color:#f2f2f2;\">\n<tr>\n<th>A<\/th>\n<th>B<\/th>\n<th>mean(A)<\/th>\n<th>mean(B)<\/th>\n<th>diff<\/th>\n<th>se<\/th>\n<th>T<\/th>\n<th>df<\/th>\n<th>pval<\/th>\n<th>hedges<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td>Technique A<\/td>\n<td>Technique B<\/td>\n<td>83,0000<\/td>\n<td>77,0000<\/td>\n<td>6,0000<\/td>\n<td>4,5583<\/td>\n<td>1,3163<\/td>\n<td>10,5034<\/td>\n<td>0,4171<\/td>\n<td>0,5516<\/td>\n<\/tr>\n<tr>\n<td>Technique A<\/td>\n<td>Technique C<\/td>\n<td>83,0000<\/td>\n<td>61,0000<\/td>\n<td>22,0000<\/td>\n<td>4,4172<\/td>\n<td>4,9805<\/td>\n<td>9,3134<\/td>\n<td>0,0017<\/td>\n<td>2,0871<\/td>\n<\/tr>\n<tr>\n<td>Technique B<\/td>\n<td>Technique C<\/td>\n<td>77,0000<\/td>\n<td>61,0000<\/td>\n<td>16,0000<\/td>\n<td>1,3904<\/td>\n<td>11,5071<\/td>\n<td>12,4164<\/td>\n<td>0,0000<\/td>\n<td>4,8222<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<h3>9.3 Gestion du risque familial d&rsquo;erreur et ajustements de p-valeurs<\/h3>\n<p>Le tableau post-hoc de Games-Howell int\u00e8gre un m\u00e9canisme interne de correction de la loi du rang Studentis\u00e9, contr\u00f4lant le taux d&rsquo;erreur par famille (<em>FWER &#8211; Family-Wise Error Rate<\/em>) au seuil global de 5 % pour l&rsquo;ensemble des trois contrastes simultan\u00e9s.<\/p>\n<p>L&rsquo;analyse d\u00e9taill\u00e9e des comparaisons met en lumi\u00e8re des trajectoires cliniques contrast\u00e9es :<\/p>\n<ul>\n<li><strong>Comparaison Technique A vs Technique B :<\/strong> La diff\u00e9rence observ\u00e9e est de +6,00 points en faveur de la r\u00e9cup\u00e9ration active. Toutefois, l&rsquo;erreur type bilat\u00e9rale est substantielle (<em>SE = 4,5583<\/em>) en raison de la variance massive du groupe A, abaissant les degr\u00e9s de libert\u00e9 ajust\u00e9s \u00e0 <em>df = 10,50<\/em>. La valeur de <em>T<\/em> atteint 1,3163, d\u00e9bouchant sur une <em>p<\/em>-valeur corrig\u00e9e de <strong>0,4171<\/strong>. Ce contraste n&rsquo;est pas statistiquement significatif au seuil &alpha; = 0,05. Malgr\u00e9 un \u00e9cart descriptif visible, nous ne pouvons affirmer que la Technique A est sup\u00e9rieure \u00e0 la Technique B dans la population parente.<\/li>\n<li><strong>Comparaison Technique A vs Technique C :<\/strong> La diff\u00e9rence de performance atteint +22,00 points. Malgr\u00e9 la variance \u00e9lev\u00e9e de A, la magnitude de cet \u00e9cart produit une statistique <em>T<\/em> de 4,9805 pour <em>df = 9,31<\/em>, yielding une <em>p<\/em>-valeur corrig\u00e9e de <strong>0,0017<\/strong>. Le gain de performance de la Technique A par rapport \u00e0 la simple relecture passive est hautement significatif.<\/li>\n<li><strong>Comparaison Technique B vs Technique C :<\/strong> La diff\u00e9rence s&rsquo;\u00e9tablit \u00e0 +16,00 points. Les deux groupes b\u00e9n\u00e9ficiant de variances restreintes, l&rsquo;erreur type est remarquablement basse (<em>SE = 1,3904<\/em>), propulsant la statistique <em>T<\/em> \u00e0 11,5071 (<em>df = 12,42<\/em>). La <em>p<\/em>-valeur r\u00e9sultante est infinit\u00e9simale (<strong>p &lt; 0,0001<\/strong>). La cartographie conceptuelle surpasse massivement la relecture passive.<\/li>\n<\/ul>\n<p>Cette d\u00e9composition illustre la subtilit\u00e9 des conclusions autoris\u00e9es par la mod\u00e9lisation de Welch-Games-Howell : elle d\u00e9montre que la significativit\u00e9 globale du test omnibus \u00e9tait tir\u00e9e conjointement par les contrastes A vs C et B vs C, tandis que les deux techniques actives (A et B) ne se distinguent pas de fa\u00e7on probante dans cet \u00e9chantillon.<\/p>\n<h2>10. \u00c9tape 5 : Calcul de la taille d&rsquo;effet et puissance statistique<\/h2>\n<p>L&rsquo;inf\u00e9rence contemporaine ne saurait se borner au simple constat d&rsquo;un seuil binaire de significativit\u00e9 statistique (<em>p &lt; 0,05<\/em>). La quantification de l&rsquo;amplitude de l&rsquo;effet, assortie d&rsquo;une \u00e9valuation de la puissance du dispositif, conditionne l&rsquo;\u00e9valuation de la pertinence pratique et de la r\u00e9plicabilit\u00e9 scientifique des r\u00e9sultats.<\/p>\n<h3>10.1 Le biais de l&rsquo;\u00c9ta carr\u00e9 (&eta;<sup>2<\/sup>) en pr\u00e9sence d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9<\/h3>\n<p>Dans l&rsquo;ANOVA classique, l&rsquo;indice de taille d&rsquo;effet le plus couramment rapport\u00e9 est l&rsquo;\u00c9ta carr\u00e9 (<em>&eta;<sup>2<\/sup><\/em>) ou son \u00e9quivalent partiel (<em>&eta;<sub>p<\/sub><sup>2<\/sup><\/em>), d\u00e9fini comme le ratio de la somme des carr\u00e9s du facteur sur la somme des carr\u00e9s totale :<\/p>\n<p><em>&eta;<sup>2<\/sup> = SS<sub>facteur<\/sub> \/ SS<sub>total<\/sub><\/em><\/p>\n<p>Cet estimateur pr\u00e9sente deux d\u00e9faillances critiques bien identifi\u00e9es par les statisticiens :<\/p>\n<p>En premier lieu, il s&rsquo;agit d&rsquo;un estimateur descriptif de l&rsquo;\u00e9chantillon et non d&rsquo;un estimateur non biais\u00e9 de la population parente. Il a tendance \u00e0 surestimer syst\u00e9matiquement la part r\u00e9elle de variance expliqu\u00e9e, particuli\u00e8rement sur de petites cohortes (<em>N &le; 30<\/em>).<\/p>\n<p>En second lieu, le calcul de la somme des carr\u00e9s totale postule implicitement une variance r\u00e9siduelle homog\u00e8ne. Lorsque les variances empiriques \u00e9clatent, la m\u00e9trique de l&rsquo;\u00c9ta carr\u00e9 se trouve distordue par les groupes \u00e0 dispersion maximale, fournissant une estimation trompeuse de la proportion de variance attribuable au facteur manipul\u00e9. Son emploi dans le contexte d&rsquo;une ANOVA de Welch est donc d\u00e9conseill\u00e9 au profit d&rsquo;indices corrig\u00e9s.<\/p>\n<h3>10.2 Calcul de l&rsquo;Om\u00e9ga carr\u00e9 partiel et de l&rsquo;Epsilon carr\u00e9<\/h3>\n<p>Pour mesurer fid\u00e8lement l&rsquo;effet omnibus au sein du mod\u00e8le de Welch, il convient de privil\u00e9gier l&rsquo;<strong>Om\u00e9ga carr\u00e9 (<em>&omega;<sup>2<\/sup><\/em>)<\/strong> ou l&rsquo;<strong>Epsilon carr\u00e9 (<em>&epsilon;<sup>2<\/sup><\/em>)<\/strong>. L&rsquo;Om\u00e9ga carr\u00e9 ajuste la d\u00e9viation observ\u00e9e en retranchant de la somme des carr\u00e9s du facteur l&rsquo;estimation de l&rsquo;erreur d&rsquo;\u00e9chantillonnage, fournissant un estimateur sans biais de la magnitude du ph\u00e9nom\u00e8ne \u00e0 l&rsquo;\u00e9chelle populationnelle.<\/p>\n<p>Dans le contexte sp\u00e9cifique de l&rsquo;ANOVA de Welch, une adaptation math\u00e9matique de l&rsquo;Om\u00e9ga carr\u00e9 ajust\u00e9 aux degr\u00e9s de libert\u00e9 fractionnaires de Satterthwaite s&rsquo;applique selon la formule suivante :<\/p>\n<p><em>&omega;<sub>Welch<\/sub><sup>2<\/sup> = [ df<sub>1<\/sub> &times; (F<sub>w<\/sub> &#8211; 1) ] \/ [ df<sub>1<\/sub> &times; (F<sub>w<\/sub> &#8211; 1) + df<sub>2(Welch)<\/sub> + 1 ]<\/em><\/p>\n<p>Calculons cet indice manuellement en Python \u00e0 partir des param\u00e8tres extraits pr\u00e9c\u00e9demment :<\/p>\n<p><code>df1 = 2<\/code><br \/>\n<code>df2_w = 12.7412<\/code><br \/>\n<code>F_w = 21.8357<\/code><\/p>\n<p><code>omega2_welch = (df1 * (F_w - 1.0)) \/ (df1 * (F_w - 1.0) + df2_w + 1.0)<\/code><br \/>\n<code>print(f\"Om\u00e9ga carr\u00e9 ajust\u00e9 de Welch : {omega2_welch:.4f}\")<\/code><\/p>\n<p>La valeur num\u00e9rique obtenue est :<\/p>\n<p><em>&omega;<sub>Welch<\/sub><sup>2<\/sup> = 0,7412<\/em><\/p>\n<p>Selon les conventions \u00e9tablies par Jacob Cohen (1988), enrichies pour les contextes comportementaux, les seuils d&rsquo;interpr\u00e9tation s&rsquo;\u00e9tablissent comme suit :<\/p>\n<ul>\n<li><strong>Effet modeste (petit) :<\/strong> <em>&omega;<sup>2<\/sup> &approx; 0,01<\/em> (1 % de la variance expliqu\u00e9e).<\/li>\n<li><strong>Effet interm\u00e9diaire (moyen) :<\/strong> <em>&omega;<sup>2<\/sup> &approx; 0,06<\/em> (6 % de la variance expliqu\u00e9e).<\/li>\n<li><strong>Effet substantiel (grand) :<\/strong> <em>&omega;<sup>2<\/sup> &ge; 0,14<\/em> (14 % ou plus de la variance expliqu\u00e9e).<\/li>\n<\/ul>\n<p>Avec un <em>&omega;<sub>Welch<\/sub><sup>2<\/sup><\/em> culminant \u00e0 <strong>0,7412<\/strong>, plus de 74 % de la variance totale des performances \u00e0 l&rsquo;examen est attribuable \u00e0 la m\u00e9thode d&rsquo;apprentissage s\u00e9lectionn\u00e9e. Il s&rsquo;agit d&rsquo;une taille d&rsquo;effet exceptionnellement massive \u00e0 l&rsquo;\u00e9chelle des interventions en psychologie cognitive.<\/p>\n<p>Pour les comparaisons par paires, le tableau de Games-Howell (chapitre 9.2) int\u00e8gre l&rsquo;indice <em>g<\/em> de Hedges (analogue sans biais du <em>d<\/em> de Cohen adapt\u00e9 aux petits \u00e9chantillons). Nous relevons un <em>g<\/em> mod\u00e9r\u00e9 entre les Techniques A et B (<em>g = 0,55<\/em>), et des magnitudes extr\u00eames entre A et C (<em>g = 2,09<\/em>) ainsi qu&rsquo;entre B et C (<em>g = 4,82<\/em>).<\/p>\n<h3>10.3 Analyse de puissance statistique a posteriori et prospective<\/h3>\n<p>La puissance statistique (not\u00e9e <em>1 &#8211; &beta;<\/em>) correspond \u00e0 la probabilit\u00e9 de rejeter l&rsquo;hypoth\u00e8se nulle lorsqu&rsquo;un effet r\u00e9el de magnitude donn\u00e9e existe dans la population parente. Le calcul de puissance s&rsquo;articule autour de deux approches distinctes :<\/p>\n<p>L&rsquo;<strong>analyse a posteriori<\/strong> \u00e9value la puissance effective atteinte par le dispositif exp\u00e9rimental. Sous Pingouin, la fonction <code>pg.power_anova()<\/code> permet de sonder ce niveau :<\/p>\n<p><code># Conversion de l'effet en indice f de Cohen : f = sqrt(omega^2 \/ (1 - omega^2))<\/code><br \/>\n<code>f_cohen = np.sqrt(omega2_welch \/ (1.0 - omega2_welch))<\/code><br \/>\n<code>puissance_obtenue = pg.power_anova(eta_squared=omega2_welch, k=3, n=10, alpha=0.05)<\/code><br \/>\n<code>print(f\"Indice f de Cohen = {f_cohen:.4f} ; Puissance observ\u00e9e = {puissance_obtenue:.4f}\")<\/code><\/p>\n<p>Avec un effet d&rsquo;une telle envergure, la puissance atteinte avec 10 sujets par cohorte exc\u00e8de <strong>0,999<\/strong> (largement au-del\u00e0 du standard acad\u00e9mique de 0,80). La probabilit\u00e9 d&rsquo;avoir commis une fausse non-d\u00e9couverte (erreur de deuxi\u00e8me esp\u00e8ce) sur l&rsquo;effet global est infinit\u00e9simale.<\/p>\n<p>L&rsquo;<strong>analyse prospective<\/strong> s&rsquo;int\u00e9resse \u00e0 la r\u00e9plication future de l&rsquo;\u00e9tude. Si des chercheurs souhaitent r\u00e9pliquer sp\u00e9cifiquement le contraste subtil entre la Technique A et la Technique B (dont l&rsquo;effet mesur\u00e9 \u00e9tait <em>g = 0,55<\/em>, non significatif \u00e0 <em>n = 10<\/em>), quelle taille d&rsquo;\u00e9chantillon devront-ils mobiliser pour atteindre une puissance nominale de 80 % au seuil bilat\u00e9ral &alpha; = 0,05 ?<\/p>\n<p><code>n_requis = pg.power_ttest2n(nx=None, ny=None, d=0.55, power=0.80, alpha=0.05)<\/code><br \/>\n<code>print(f\"Effectif requis par condition pour d\u00e9tecter le contraste A vs B : {int(np.ceil(n_requis))} participants\")<\/code><\/p>\n<p>Le calcul indique qu&rsquo;un minimum de <strong>53 participants par groupe<\/strong> (soit une cohorte globale de 106 \u00e9tudiants pour ce contraste) est indispensable pour trancher statistiquement entre les deux approches actives d&rsquo;apprentissage. Cette analyse r\u00e9trospective d\u00e9montre que le design initial (<em>n = 10<\/em>) manquait de puissance statistique pour discriminer des effets de taille moyenne, justifiant la prudence des conclusions post-hoc de Games-Howell.<\/p>\n<h2>11. Visualisation graphique avanc\u00e9e des r\u00e9sultats en Python<\/h2>\n<p>Dans un rapport empirique de haut niveau, la communication graphique des donn\u00e9es ne saurait se r\u00e9duire \u00e0 un simple diagramme en barres tronquant l&rsquo;axe vertical (les tristement c\u00e9l\u00e8bres <em>dynamite plots<\/em>). Une visualisation contemporaine doit combiner la transparence distributionnelle, la clart\u00e9 de la tendance centrale et la restitution fid\u00e8le de l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des dispersions.<\/p>\n<h3>11.1 Repr\u00e9sentation combin\u00e9e : Violin plot et Boxplot<\/h3>\n<p>La superposition d&rsquo;un diagramme en violon (<em>violin plot<\/em>), d&rsquo;un diagramme en bo\u00eete (<em>boxplot<\/em>) et d&rsquo;un nuage de points individuels dispers\u00e9s (<em>jittered stripplot<\/em>) constitue le standard d&rsquo;excellence \u00e9ditoriale. Cette triangulation permet au lecteur de percevoir simultan\u00e9ment la densit\u00e9 non param\u00e9trique de l&rsquo;\u00e9chantillon, les quartiles m\u00e9dians et la variabilit\u00e9 individuelle des observations.<\/p>\n<p>Concevons cette figure sous Seaborn et Matplotlib :<\/p>\n<p><code>plt.figure(figsize=(10, 6), dpi=300)<\/code><\/p>\n<p><code># 1. Trac\u00e9 de la densit\u00e9 par violon<\/code><br \/>\n<code>sns.violinplot(x='Groupe', y='Score', data=df_etude, inner=None, color='#E8EEF5', linewidth=1.2)<\/code><\/p>\n<p><code># 2. Insertion d'une bo\u00eete \u00e0 moustaches \u00e9troite<\/code><br \/>\n<code>sns.boxplot(x='Groupe', y='Score', data=df_etude, width=0.25, boxprops=dict(facecolor='#BDD7EE', alpha=0.8),<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;whiskerprops=dict(color='#2E5B82', linewidth=1.5), medianprops=dict(color='#C00000', linewidth=2))<\/code><\/p>\n<p><code># 3. Superposition des points bruts individuels avec jittering<\/code><br \/>\n<code>sns.stripplot(x='Groupe', y='Score', data=df_etude, size=7, color='#1F4E79', alpha=0.7, jitter=0.2)<\/code><\/p>\n<p><code>plt.title(\"Distribution des performances acad\u00e9miques selon la technique d'\u00e9tuden(Superposition Violon, Boxplot et Donn\u00e9es Individuelles)\", fontsize=13, fontweight='bold', pad=15)<\/code><br \/>\n<code>plt.xlabel(\"Condition exp\u00e9rimentale (VI)\", fontsize=11, fontweight='bold')<\/code><br \/>\n<code>plt.ylabel(\"Score standardis\u00e9 \u00e0 l'examen (VD)\", fontsize=11, fontweight='bold')<\/code><br \/>\n<code>plt.ylim(50, 105)<\/code><br \/>\n<code>plt.tight_layout()<\/code><br \/>\n<code>plt.show()<\/code><\/p>\n<p>Ce graphique expose visuellement la probl\u00e9matique fondamentale trait\u00e9e par l&rsquo;ANOVA de Welch : la forme allong\u00e9e et dilat\u00e9e du violon de la Technique A tranche radicalement avec le violon extr\u00eamement compact et resserr\u00e9 de la Technique C, offrant une justification graphique imm\u00e9diate au rejet de l&rsquo;homosc\u00e9dasticit\u00e9.<\/p>\n<h3>11.2 Graphique des moyennes ajust\u00e9es avec intervalles de confiance Welch<\/h3>\n<p>Pour illustrer directement les r\u00e9sultats de l&rsquo;ANOVA de Welch et des contrastes de Games-Howell, il est opportun de g\u00e9n\u00e9rer un graphique des moyennes empiriques flanqu\u00e9es de barres d&rsquo;erreur bas\u00e9es sur les <strong>intervalles de confiance sp\u00e9cifiques \u00e0 chaque groupe<\/strong> (calcul\u00e9s via l&rsquo;\u00e9cart-type propre de chaque modalit\u00e9 sous distribution de Student \u00e0 <em>n<sub>i<\/sub> &#8211; 1<\/em> degr\u00e9s de libert\u00e9), enrichi d&rsquo;annotations de significativit\u00e9 post-hoc :<\/p>\n<p><code>fig, ax = plt.subplots(figsize=(9, 6), dpi=300)<\/code><\/p>\n<p><code># Extraction des statistiques descriptives<\/code><br \/>\n<code>x_pos = np.arange(len(stats_desc))<\/code><br \/>\n<code>moyennes = stats_desc['Moyenne']<\/code><br \/>\n<code>erreurs = stats_desc['Ecart_Type'] \/ np.sqrt(stats_desc['Effectif'])<\/code><\/p>\n<p><code># Barres d'erreur et points de moyenne<\/code><br \/>\n<code>ax.errorbar(x_pos, moyennes, yerr=erreurs * stats.t.ppf(0.975, df=stats_desc['Effectif']-1),<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;fmt='o', color='#1F4E79', ecolor='#1F4E79', elinewidth=2.5, capsize=6, capthick=2, markersize=8)<\/code><\/p>\n<p><code># Annotations de contrastes post-hoc<\/code><br \/>\n<code># Contraste A vs C (Significatif p = 0.0017)<\/code><br \/>\n<code>y_bar_ac = 100<\/code><br \/>\n<code>ax.plot([0, 0, 2, 2], [y_bar_ac - 1, y_bar_ac, y_bar_ac, y_bar_ac - 1], color='black', lw=1.2)<\/code><br \/>\n<code>ax.text(1, y_bar_ac + 0.5, \"** (p = 0.0017)\", ha='center', va='bottom', color='black', fontsize=10, fontweight='bold')<\/code><\/p>\n<p><code># Contraste B vs C (Significatif p < 0.0001)<\/code><br \/>\n<code>y_bar_bc = 87<\/code><br \/>\n<code>ax.plot([1, 1, 2, 2], [y_bar_bc - 1, y_bar_bc, y_bar_bc, y_bar_bc - 1], color='black', lw=1.2)<\/code><br \/>\n<code>ax.text(1.5, y_bar_bc + 0.5, \"*** (p < 0.0001)\", ha='center', va='bottom', color='black', fontsize=10, fontweight='bold')<\/code><\/p>\n<p><code># Configuration des axes<\/code><br \/>\n<code>ax.set_xticks(x_pos)<\/code><br \/>\n<code>ax.set_xticklabels(stats_desc['Groupe'], fontsize=11, fontweight='bold')<\/code><br \/>\n<code>ax.set_ylabel(\"Score moyen \u00e0 l'examen (avec IC \u00e0 95 %)\", fontsize=11, fontweight='bold')<\/code><br \/>\n<code>ax.set_title(\"Comparaison des moyennes et contrastes post-hoc de Games-Howell\", fontsize=12, fontweight='bold', pad=15)<\/code><br \/>\n<code>ax.set_ylim(55, 108)<\/code><br \/>\n<code>plt.tight_layout()<\/code><br \/>\n<code>plt.show()<\/code><\/p>\n<p>Cette visualisation synth\u00e9tise l'ensemble du corpus inf\u00e9rentiel en mettant en relief la largeur imposante de la barre d'erreur du premier groupe en comparaison de l'\u00e9troitesse des barres d'erreur des deux autres conditions.<\/p>\n<h3>11.3 Exportation des figures haute r\u00e9solution pour publication<\/h3>\n<p>Les maisons d'\u00e9dition scientifique (telles que l'American Psychological Association, Elsevier ou Springer) imposent des directives strictes concernant la d\u00e9finition et le format vectoriel des illustrations. Pour se conformer aux standards de soumission acad\u00e9mique, l'exportation des graphiques doit respecter les crit\u00e8res suivants :<\/p>\n<ul>\n<li><strong>R\u00e9solution raster minimale :<\/strong> Pour les fichiers image tram\u00e9s (matriciels), une r\u00e9solution minimale de 300 \u00e0 600 DPI (<em>dots per inch<\/em>) est exig\u00e9e afin d'\u00e9viter toute pixelisation lors du processus d'impression offset ou de compression num\u00e9rique. Le format TIFF sans perte ou le PNG haute d\u00e9finition sont pr\u00e9conis\u00e9s.<\/li>\n<li><strong>Formats vectoriels scalables :<\/strong> La m\u00e9thode optimale consiste \u00e0 exporter les figures en format vectoriel pur, tel que le PDF (<em>Portable Document Format<\/em>) ou le SVG (<em>Scalable Vector Graphics<\/em>). Dans ce mode de stockage, les caract\u00e8res typographiques, les lignes de tendance et les symboles d'erreur sont d\u00e9crits sous forme d'instructions g\u00e9om\u00e9triques parfaites, garantissant un rendu impeccable quel que soit le facteur d'agrandissement.<\/li>\n<li><strong>Accessibilit\u00e9 chromatique :<\/strong> Le choix des palettes de couleurs doit respecter les principes d'accessibilit\u00e9 visuelle universelle, en \u00e9vitant les juxtapositions rouge-vert probl\u00e9matiques pour les personnes pr\u00e9sentant une d\u00e9ficience de vision des couleurs (daltonisme). Les palettes perceptuellement uniformes (telles que <em>viridis<\/em>, <em>cividis<\/em> ou la palette <em>colorblind<\/em> de Seaborn) assurent une lisibilit\u00e9 optimale y compris en impression monochrome en niveaux de gris.<\/li>\n<\/ul>\n<p>Sous Matplotlib, l'exportation s'op\u00e8re au moyen de la m\u00e9thode <code>savefig()<\/code> configur\u00e9e avec les arguments adapt\u00e9s :<\/p>\n<p><code>plt.savefig(\"resultats_anova_welch.pdf\", format=\"pdf\", bbox_inches=\"tight\")<\/code><br \/>\n<code>plt.savefig(\"resultats_anova_welch.tiff\", format=\"tiff\", dpi=300, pil_kwargs={\"compression\": \"tiff_lzw\"}, bbox_inches=\"tight\")<\/code><\/p>\n<p>L'utilisation de <code>bbox_inches=\"tight\"<\/code> garantit l'\u00e9limination automatique de toutes les marges blanches inutiles autour de l'illustration.<\/p>\n<h2>12. Rapport acad\u00e9mique des r\u00e9sultats et directives de publication APA<\/h2>\n<p>L'analyse computationnelle ne trouve son accomplissement que lorsqu'elle est retranscrite avec exactitude au sein de la section \u00ab R\u00e9sultats \u00bb d'un manuscrit scientifique. Le manuel de style de l'APA (7e \u00e9dition) codifie de fa\u00e7on formelle la syntaxe r\u00e9dactionnelle des mod\u00e8les d'analyse de variance robuste.<\/p>\n<h3>12.1 Normes de r\u00e9daction statistique selon le style APA 7e \u00e9dition<\/h3>\n<p>La r\u00e9daction formelle des r\u00e9sultats issus d'une ANOVA de Welch et de ses proc\u00e9dures d\u00e9riv\u00e9es doit imp\u00e9rativement respecter les standards typographiques et statistiques suivants :<\/p>\n<ul>\n<li><strong>Statistique de test :<\/strong> La statistique <em>F<\/em> de Welch est d\u00e9sign\u00e9e par la lettre italique majuscule <em>F<\/em>, fr\u00e9quemment compl\u00e9t\u00e9e de la mention en indice ou textuelle pr\u00e9cisant la nature robuste du test (ex. <em>F<\/em> de Welch ou <em>F<sub>Welch<\/sub><\/em>).<\/li>\n<li><strong>Degr\u00e9s de libert\u00e9 :<\/strong> Les degr\u00e9s de libert\u00e9 sont syst\u00e9matiquement consign\u00e9s entre parenth\u00e8ses imm\u00e9diatement apr\u00e8s la lettre <em>F<\/em>, s\u00e9par\u00e9s par une virgule. Les degr\u00e9s de libert\u00e9 du num\u00e9rateur sont entiers, tandis que <strong>les degr\u00e9s de libert\u00e9 du d\u00e9nominateur doivent obligatoirement conserver leurs d\u00e9cimales<\/strong> (g\u00e9n\u00e9ralement arrondies \u00e0 deux chiffres apr\u00e8s la virgule selon l'usage standard APA, ex. <em>F(2, 12,74)<\/em>). Tronquer ou arrondir arbitrairement cette valeur \u00e0 l'entier le plus proche constitue une d\u00e9rive m\u00e9thodologique p\u00e9nalisante qui masque le recours \u00e0 la correction de Satterthwaite.<\/li>\n<li><strong>Niveau de probabilit\u00e9 :<\/strong> La <em>p<\/em>-valeur s'inscrit en italique minuscule (<em>p<\/em>). Lorsque sa valeur est inf\u00e9rieure \u00e0 0,001, elle est rapport\u00e9e sous la forme <em>p &lt; 0,001<\/em> plut\u00f4t que <em>p = 0,0000<\/em>.<\/li>\n<li><strong>M\u00e9trique de taille d'effet :<\/strong> L'indice de taille d'effet ajust\u00e9 (Om\u00e9ga carr\u00e9, not\u00e9 <em>&omega;<sup>2<\/sup><\/em>, ou \u00c9ta carr\u00e9 partiel <em>&eta;<sub>p<\/sub><sup>2<\/sup><\/em>) doit accompagner syst\u00e9matiquement la statistique omnibus.<\/li>\n<li><strong>Dispersion intra-groupe :<\/strong> Pour chaque condition exp\u00e9rimentale mentionn\u00e9e dans le texte narratif, la moyenne (<em>M<\/em>) et l'\u00e9cart-type (<em>SD<\/em>) doivent \u00eatre syst\u00e9matiquement formul\u00e9s pour informer le lecteur des disparit\u00e9s de dispersion constat\u00e9es.<\/li>\n<\/ul>\n<h3>12.2 Exemple concret de paragraphe de r\u00e9sultats pour article empirique<\/h3>\n<p>Voici le mod\u00e8le de paragraphe acad\u00e9mique pr\u00eat \u00e0 l'insertion \u00e9ditoriale, int\u00e9grant de fa\u00e7on fluide les diagnostics d'h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9, le test global de Welch, les contrastes post-hoc de Games-Howell et les tailles d'effet associ\u00e9es :<\/p>\n<blockquote>\n<p>Une analyse de variance unifactorielle inter-sujets a \u00e9t\u00e9 conduite afin d'\u00e9valuer l'effet de la technique d'apprentissage (Technique A : R\u00e9cup\u00e9ration active ; Technique B : Cartographie conceptuelle ; Technique C : Relecture passive) sur les scores obtenus \u00e0 l'examen de synth\u00e8se conceptuelle. Les analyses exploratoires pr\u00e9liminaires ont confirm\u00e9 la normalit\u00e9 distributionnelle des scores au sein de chaque cohorte exp\u00e9rimentale au moyen des tests de Shapiro-Wilk (toutes les valeurs de <em>p<\/em> &ge; 0,522). Toutefois, l'hypoth\u00e8se fondamentale d'homog\u00e9n\u00e9it\u00e9 des variances intra-groupes a \u00e9t\u00e9 formellement rejet\u00e9e par le test d'\u00e9galit\u00e9 des variances de Bartlett, <em>&chi;<sup>2<\/sup>(2) = 25,48<\/em>, <em>p &lt; 0,001<\/em>, ainsi que par le test de Levene centr\u00e9 sur la m\u00e9diane, <em>W(2, 27) = 10,87<\/em>, <em>p &lt; 0,001<\/em>, r\u00e9v\u00e9lant une dispersion consid\u00e9rablement accrue au sein du groupe soumis \u00e0 la r\u00e9cup\u00e9ration active (<em>SD = 13,85<\/em>) comparativement aux groupes de cartographie conceptuelle (<em>SD = 4,00<\/em>) et de relecture passive (<em>SD = 1,83<\/em>).<\/p>\n<p>En cons\u00e9quence de cette h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 structurelle, l'ANOVA robuste de Welch a \u00e9t\u00e9 mobilis\u00e9e. L'effet omnibus de la technique d'apprentissage s'est r\u00e9v\u00e9l\u00e9 hautement significatif, <em>F<sub>Welch<\/sub>(2, 12,74) = 21,84<\/em>, <em>p &lt; 0,001<\/em>, <em>&omega;<sup>2<\/sup> = 0,74<\/em>, traduisant une taille d'effet majeure selon les crit\u00e8res de Cohen. Les comparaisons par paires r\u00e9alis\u00e9es au moyen du test post-hoc de Games-Howell (ajust\u00e9 pour les variances h\u00e9t\u00e9rog\u00e8nes) indiquent que les participants de la condition R\u00e9cup\u00e9ration active (<em>M = 83,00<\/em>, <em>SD = 13,85<\/em>) ont obtenu des performances significativement sup\u00e9rieures \u00e0 ceux de la condition de Relecture passive (<em>M = 61,00<\/em>, <em>SD = 1,83<\/em> ; <em>diff\u00e9rence moyenne = 22,00<\/em>, <em>SE = 4,42<\/em>, <em>t(9,31) = 4,98<\/em>, <em>p = 0,002<\/em>, <em>g de Hedges = 2,09<\/em>). De m\u00eame, la condition Cartographie conceptuelle (<em>M = 77,00<\/em>, <em>SD = 4,00<\/em>) a surclass\u00e9 de mani\u00e8re tr\u00e8s significative la condition Relecture passive (<em>diff\u00e9rence moyenne = 16,00<\/em>, <em>SE = 1,39<\/em>, <em>t(12,42) = 11,51<\/em>, <em>p &lt; 0,001<\/em>, <em>g de Hedges = 4,82<\/em>). En revanche, le contraste entre la R\u00e9cup\u00e9ration active et la Cartographie conceptuelle n'a pas atteint le seuil de significativit\u00e9 statistique conventionnel (<em>diff\u00e9rence moyenne = 6,00<\/em>, <em>SE = 4,56<\/em>, <em>t(10,50) = 1,32<\/em>, <em>p = 0,417<\/em>, <em>g de Hedges = 0,55<\/em>).<\/p>\n<\/blockquote>\n<h3>12.3 Bonnes pratiques et \u00e9cueils m\u00e9thodologiques courants<\/h3>\n<p>Pour parachever cette d\u00e9marche, r\u00e9capitulons les trois pi\u00e8ges m\u00e9thodologiques les plus fr\u00e9quents r\u00e9pertori\u00e9s dans la litt\u00e9rature empirique lors de l'usage de l'ANOVA de Welch :<\/p>\n<ul>\n<li><strong>Erreur d'arrondi sur les degr\u00e9s de libert\u00e9 de Satterthwaite :<\/strong> L'un des \u00e9cueils les plus pernicieux consiste \u00e0 arrondir manuellement <em>df<sub>2<\/sub><\/em> \u00e0 l'entier le plus proche (par exemple, rapporter <em>F(2, 13)<\/em> au lieu de <em>F(2, 12,74)<\/em>). Outre le fait que cette pratique viole les directives explicites du manuel APA 7e \u00e9dition, elle emp\u00eache les relecteurs d'identifier avec certitude si l'auteur a appliqu\u00e9 l'ANOVA de Welch ou s'il a ex\u00e9cut\u00e9 une ANOVA standard de Fisher avec des donn\u00e9es manquantes. La pr\u00e9sence de d\u00e9cimales constitue la marque distinctive de la correction robuste.<\/li>\n<li><strong>Omission injustifi\u00e9e du test de variance pr\u00e9alable :<\/strong> Appliquer l'ANOVA de Welch sans avoir document\u00e9 l'\u00e9tat des variances (via Levene, Bartlett ou l'analyse des \u00e9carts-types) prive le lecteur de la justification logique du choix analytique. Bien que plusieurs m\u00e9thodologistes sugg\u00e8rent d'utiliser l'ANOVA de Welch par d\u00e9faut dans toutes les situations en raison de sa convergence vers Fisher sous homosc\u00e9dasticit\u00e9, la tradition acad\u00e9mique exige de motiver formellement l'abandon du mod\u00e8le standard de Fisher.<\/li>\n<li><strong>Confusion conceptuelle entre robustesse param\u00e9trique et test non param\u00e9trique de rang :<\/strong> Une confusion conceptuelle tenace assimile l'ANOVA de Welch \u00e0 un test non param\u00e9trique. L'ANOVA de Welch demeure un test <strong>pleinement param\u00e9trique<\/strong> : elle \u00e9value des moyennes arithm\u00e9tiques vraies et conserve la structure m\u00e9trique des donn\u00e9es originales continues. La qualifier de test non param\u00e9trique dans un article t\u00e9moigne d'une m\u00e9connaissance de ses fondements math\u00e9matiques. Pour les donn\u00e9es purement ordinales ou comportant des valeurs aberrantes massives et irr\u00e9ductibles qui violent irr\u00e9m\u00e9diablement la normalit\u00e9 locale, seul le test de Kruskal-Wallis ou les mod\u00e8les bas\u00e9s sur les moyennes tronqu\u00e9es de Wilcox (<em>trimmed means<\/em>) s'av\u00e8rent fond\u00e9s.<\/li>\n<\/ul>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Bartlett, M. S. (1937). Properties of sufficiency and statistical tests. <em>Proceedings of the Royal Society of London. Series A - Mathematical and Physical Sciences<\/em>, 160(901), 268\u2013282. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1098\/rspa.1937.0109\">https:\/\/doi.org\/10.1098\/rspa.1937.0109<\/a><\/li>\n<li>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><\/li>\n<li>Cohen, J. (1988). <em>Statistical power analysis for the behavioral sciences<\/em> (2nd ed.). Lawrence Erlbaum Associates.<\/li>\n<li>Field, A. (2018). <em>Discovering statistics using IBM SPSS statistics<\/em> (5th ed.). SAGE Publications.<\/li>\n<li>Games, P. A., &amp; Howell, J. F. (1976). Pairwise multiple comparison procedures with unequal n'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.3102\/10769986001002113\">https:\/\/doi.org\/10.3102\/10769986001002113<\/a><\/li>\n<li>George, D., &amp; Mallery, P. (2010). <em>SPSS for Windows step by step: A simple guide and reference, 17.0 update<\/em> (10th ed.). Allyn &amp; Bacon.<\/li>\n<li>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.<\/li>\n<li>Tomarken, A. J., &amp; Serlin, R. C. (1986). Comparison of ANOVA alternatives under variance heterogeneity and specific noncentrality structures. <em>Psychological Bulletin<\/em>, 99(1), 90\u201399. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0033-2909.99.1.90\">https:\/\/doi.org\/10.1037\/0033-2909.99.1.90<\/a><\/li>\n<li>Vallat, R. (2018). Pingouin: statistics in Python. <em>Journal of Open Source Software<\/em>, 3(31), 1026. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.21105\/joss.01026\">https:\/\/doi.org\/10.21105\/joss.01026<\/a><\/li>\n<li>Welch, B. L. (1947). The generalization of 'Student's' 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.2307\/2332510\">https:\/\/doi.org\/10.2307\/2332510<\/a><\/li>\n<li>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.2307\/2332579\">https:\/\/doi.org\/10.2307\/2332579<\/a><\/li>\n<li>Wilcox, R. R. (2012). <em>Introduction to robust estimation and hypothesis testing<\/em> (3rd ed.). Academic Press.<\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour r\u00e9aliser l&rsquo;ANOVA de Welch en Python pas \u00e0 pas lorsque l&rsquo;hypoth\u00e8se d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances est viol\u00e9e en psychologie.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2945,2992],"tags":[5167,2272,5010,2468,1842],"class_list":["post-8229","psych_stat","type-psych_stat","status-publish","hentry","category-methodes-statistiques-en-psychologie","category-science-des-donnees-en-psychologie","tag-anova-de-welch","tag-heteroscedasticite","tag-pingouin","tag-python","tag-statistiques-psychologiques"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8229","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=8229"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8229\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=8229"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=8229"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=8229"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}