{"id":2764,"date":"2026-09-06T20:24:09","date_gmt":"2026-09-06T20:24:09","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-effectuer-anova-deux-facteurs-stata\/"},"modified":"2026-09-06T20:24:19","modified_gmt":"2026-09-06T20:24:19","slug":"comment-effectuer-anova-deux-facteurs-stata","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-effectuer-anova-deux-facteurs-stata","title":{"rendered":"Comment effectuer une ANOVA \u00e0 deux facteurs dans Stata"},"content":{"rendered":"<p>L&rsquo;analyse de variance factorielle, couramment d\u00e9sign\u00e9e sous l&rsquo;acronyme d&rsquo;ANOVA \u00e0 deux facteurs (ou <em>two-way ANOVA<\/em>), constitue l&rsquo;un des piliers m\u00e9thodologiques les plus robustes et les plus r\u00e9pandus au sein de la recherche contemporaine en sciences comportementales, en psychologie exp\u00e9rimentale, en neurosciences et en m\u00e9decine clinique. Alors que l&rsquo;analyse univari\u00e9e simple limite l&rsquo;investigation \u00e0 l&rsquo;influence isol\u00e9e d&rsquo;un facteur unique sur une m\u00e9trique quantitative continue, le protocole factoriel permet d&rsquo;examiner simultan\u00e9ment l&rsquo;action conjointe de deux variables cat\u00e9gorielles ind\u00e9pendantes. Cette architecture analytique d\u00e9passe la simple juxtaposition de mod\u00e8les ind\u00e9pendants en ouvrant la voie \u00e0 la d\u00e9tection, \u00e0 l&rsquo;estimation et \u00e0 l&rsquo;interpr\u00e9tation des effets d&rsquo;interaction. Ces derniers permettent de d\u00e9terminer si l&rsquo;impact d&rsquo;un premier facteur se trouve modalis\u00e9, att\u00e9nu\u00e9 ou exacerb\u00e9 selon les modalit\u00e9s prises par le second facteur, refl\u00e9tant ainsi la complexit\u00e9 intrins\u00e8que des ph\u00e9nom\u00e8nes biologiques et psychologiques.<\/p>\n<p>Dans cet \u00e9cosyst\u00e8me m\u00e9thodologique, le logiciel statistique <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.stata.com\">Stata<\/a> s&rsquo;impose comme une r\u00e9f\u00e9rence incontournable au sein de la communaut\u00e9 acad\u00e9mique internationale. Gr\u00e2ce \u00e0 son moteur computationnel rigoureux, sa syntaxe concise et sa panoplie d&rsquo;outils d&rsquo;inf\u00e9rence post-estimation, Stata offre aux chercheurs une flexibilit\u00e9 d&rsquo;analyse in\u00e9gal\u00e9e. Le logiciel facilite non seulement l&rsquo;estimation des mod\u00e8les lin\u00e9aires sous-jacents, mais permet \u00e9galement d&rsquo;aborder avec une grande pr\u00e9cision le calcul des moyennes marginales pr\u00e9dictives, le trac\u00e9 graphique des interactions complexes, ainsi que la correction minutieuse des erreurs de type I inh\u00e9rentes aux comparaisons multiples. N\u00e9anmoins, l&rsquo;ex\u00e9cution aveugle de lignes de commande sans une compr\u00e9hension profonde des principes math\u00e9matiques et des postulats de distribution expose l&rsquo;analyste \u00e0 des conclusions erron\u00e9es ou \u00e0 des biais d&rsquo;interpr\u00e9tation majeurs.<\/p>\n<p>Ce guide exhaustif a \u00e9t\u00e9 con\u00e7u pour accompagner les chercheurs, biostatisticiens et doctorants dans la mise en \u0153uvre rigoureuse d&rsquo;une ANOVA factorielle \u00e0 deux facteurs sous l&rsquo;environnement Stata. De la th\u00e9orisation du mod\u00e8le lin\u00e9aire g\u00e9n\u00e9ral aux exigences formelles de publication selon les normes de l&rsquo;American Psychological Association (APA 7e \u00e9dition), chaque \u00e9tape du pipeline analytique est minutieusement diss\u00e9qu\u00e9e. Nous explorerons les fondements th\u00e9oriques des plans factoriels, la v\u00e9rification formelle des postulats d&rsquo;homosc\u00e9dasticit\u00e9 et de normalit\u00e9 r\u00e9siduelle, l&rsquo;ajustement des sommes des carr\u00e9s dans les plans d\u00e9s\u00e9quilibr\u00e9s, l&rsquo;analyse fine des effets simples, la quantification normalis\u00e9e des tailles d&rsquo;effet, ainsi que les strat\u00e9gies de mod\u00e9lisation alternative en pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 persistante.<\/p>\n<h2>1. Fondements th\u00e9oriques de l&rsquo;ANOVA \u00e0 deux facteurs en sciences comportementales<\/h2>\n<h3>1.1 D\u00e9finition et objectifs du plan factoriel<\/h3>\n<p>Le plan d&rsquo;exp\u00e9rience factoriel \u00e0 deux variables ind\u00e9pendantes discr\u00e8tes repr\u00e9sente une modalit\u00e9 exp\u00e9rimentale o\u00f9 chaque niveau d&rsquo;un premier facteur cat\u00e9goriel est crois\u00e9 avec chaque niveau d&rsquo;un second facteur cat\u00e9goriel. Contrairement aux approches r\u00e9ductionnistes qui \u00e9tudient un facteur \u00e0 la fois tout en maintenant artificiellement les autres variables constantes, le plan factoriel capture l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 naturelle des contextes d&rsquo;investigation. Dans un protocole factoriel complet, dit plan factoriel compl\u00e8tement crois\u00e9, l&rsquo;ensemble des cellules exp\u00e9rimentales issues de la combinatoire des modalit\u00e9s des deux facteurs est renseign\u00e9 par des observations empiriques, ce qui permet d&rsquo;isoler \u00e0 la fois les dynamiques propres \u00e0 chaque variable et leurs interd\u00e9pendances structurelles.<\/p>\n<p>D&rsquo;un point de vue m\u00e9thodologique, l&#8217;emploi d&rsquo;un plan factoriel pr\u00e9sente des avantages d\u00e9terminants par rapport \u00e0 la multiplication s\u00e9quentielle de tests t de Student ou d&rsquo;ANOVA univari\u00e9es \u00e0 un seul facteur. En premier lieu, il permet de pr\u00e9server l&rsquo;\u00e9conomie g\u00e9n\u00e9rale du protocole exp\u00e9rimental : un seul \u00e9chantillon global de participants permet de tester simultan\u00e9ment plusieurs hypoth\u00e8ses de recherche, r\u00e9duisant ainsi drastiquement les co\u00fbts humains et logistiques de l&rsquo;\u00e9chantillonnage. En deuxi\u00e8me lieu, l&rsquo;analyse conjointe permet de contr\u00f4ler plus efficacement la variance r\u00e9siduelle non expliqu\u00e9e. En extrayant syst\u00e9matiquement de l&rsquo;erreur exp\u00e9rimentale la part de variance attribuable au second facteur ainsi qu&rsquo;\u00e0 l&rsquo;interaction, la pr\u00e9cision statistique du test se trouve accrue, ce qui augmente substantiellement la puissance statistique globale du dispositif.<\/p>\n<p>Les applications pratiques de ce paradigme abondent en psychologie clinique, en sciences cognitives et en pharmacologie exp\u00e9rimentale. Par exemple, un chercheur en neuropsychopharmacologie peut vouloir \u00e9valuer l&rsquo;efficacit\u00e9 diff\u00e9rentielle de plusieurs mol\u00e9cules th\u00e9rapeutiques (Facteur A : Mol\u00e9cule active, Placebo, Traitement de r\u00e9f\u00e9rence) selon le g\u00e9notype sp\u00e9cifique des patients (Facteur B : All\u00e8le \u00e0 risque versus All\u00e8le protecteur). De m\u00eame, en ergonomie cognitive, on mesurera le temps de r\u00e9action face \u00e0 diff\u00e9rentes interfaces num\u00e9riques selon le niveau d&rsquo;expertise des utilisateurs (D\u00e9butants versus Experts). Ces paradigmes d\u00e9montrent que l&rsquo;effet d&rsquo;une intervention est rarement universel : il est quasi-syst\u00e9matiquement conditionn\u00e9 par les caract\u00e9ristiques individuelles ou contextuelles des sujets observ\u00e9s.<\/p>\n<p>Il est enfin fondamental d&rsquo;op\u00e9rer une distinction conceptuelle nette entre les facteurs \u00e0 effets fixes, les facteurs \u00e0 effets al\u00e9atoires et les mod\u00e8les dits mixtes. Un facteur est qualifi\u00e9 d&rsquo;effet fixe lorsque l&rsquo;ensemble exhaustif des modalit\u00e9s qui int\u00e9ressent l&rsquo;investigateur est inclus dans le design d&rsquo;\u00e9tude (par exemple, le sexe biologique ou un nombre d\u00e9termin\u00e9 de dosages pharmacologiques). Les inf\u00e9rences statistiques tir\u00e9es du mod\u00e8le se limitent alors strictement \u00e0 ces niveaux pr\u00e9cis. Inversement, un facteur est dit \u00e0 effets al\u00e9atoires lorsque les modalit\u00e9s incluses constituent un \u00e9chantillon repr\u00e9sentatif tir\u00e9 au hasard d&rsquo;une population plus vaste de modalit\u00e9s potentielles (par exemple, des classes d&rsquo;\u00e9coles, des sites cliniques ou des th\u00e9rapeutes). Enfin, les plans mixtes combinent simultan\u00e9ment au moins un facteur \u00e0 effets fixes et un facteur \u00e0 effets al\u00e9atoires, n\u00e9cessitant des ajustements computationnels sp\u00e9cifiques au niveau des carr\u00e9s moyens r\u00e9siduels.<\/p>\n<h3>1.2 Effets principaux et effet d&rsquo;interaction<\/h3>\n<p>L&rsquo;architecture logique de l&rsquo;ANOVA \u00e0 deux facteurs repose sur la partition tripartite de la variance expliqu\u00e9e entre les effets principaux des deux facteurs et leur terme d&rsquo;interaction. L&rsquo;effet principal d&rsquo;un premier facteur repr\u00e9sente l&rsquo;impact moyen global de cette variable ind\u00e9pendante sur la variable d\u00e9pendante continue, calcul\u00e9 en agr\u00e9geant l&rsquo;ensemble des modalit\u00e9s du second facteur. Math\u00e9matiquement, cela revient \u00e0 comparer les moyennes marginales de chaque modalit\u00e9 du facteur consid\u00e9r\u00e9 comme si le second facteur n&rsquo;exer\u00e7ait aucune influence modulatrice ou \u00e9tait parfaitement \u00e9quilibr\u00e9. L&rsquo;analyste examine ici la question fondamentale suivante : ind\u00e9pendamment de toute autre variable consid\u00e9r\u00e9e dans l&rsquo;\u00e9tude, existe-t-il une divergence statistiquement significative entre les niveaux de ce facteur ?<\/p>\n<p>De mani\u00e8re rigoureusement analogue, l&rsquo;effet principal du second facteur appr\u00e9hende l&rsquo;action moyenne isol\u00e9e de cette seconde variable sur la r\u00e9ponse observ\u00e9e, apr\u00e8s avoir neutralis\u00e9 les fluctuations induites par les diff\u00e9rents paliers du premier pr\u00e9dicteur. Bien que ces effets principaux offrent une vision globale du ph\u00e9nom\u00e8ne, leur interpr\u00e9tation isol\u00e9e peut s&rsquo;av\u00e9rer trompeuse en pr\u00e9sence d&rsquo;une interd\u00e9pendance fonctionnelle entre les variables explicatives. C&rsquo;est pr\u00e9cis\u00e9ment l\u00e0 qu&rsquo;intervient le concept math\u00e9matique et conceptuel de l&rsquo;interaction bidirectionnelle, formalis\u00e9 par le produit crois\u00e9 des deux facteurs dans l&rsquo;\u00e9quation lin\u00e9aire.<\/p>\n<p>L&rsquo;interaction statistique quantifie la mesure dans laquelle l&rsquo;effet d&rsquo;un facteur varie en amplitude ou en direction en fonction du niveau auquel se situe le second facteur. Sur le plan th\u00e9orique, l&rsquo;existence d&rsquo;une interaction indique que les deux facteurs ne se contentent pas d&rsquo;additionner leurs influences respectives de mani\u00e8re ind\u00e9pendante, mais agissent de concert pour produire un r\u00e9sultat singulier. Deux grandes typologies d&rsquo;interactions m\u00e9ritent une attention analytique rigoureuse : les interactions ordinales et les interactions crois\u00e9es (ou non-ordinales).<\/p>\n<p>Dans une interaction ordinale, l&rsquo;amplitude de l&rsquo;effet d&rsquo;un facteur varie selon les niveaux de l&rsquo;autre, mais la hi\u00e9rarchie relative des moyennes demeure constante : les courbes de profil repr\u00e9sent\u00e9es graphiquement ne se coupent pas, conservant le m\u00eame sens de pente bien que l&rsquo;\u00e9cart entre elles fluctue. En revanche, dans une interaction crois\u00e9e, les droites ou courbes de r\u00e9ponse se croisent de mani\u00e8re flagrante : l&rsquo;effet d&rsquo;un facteur peut \u00eatre hautement positif sous une condition donn\u00e9e, devenir totalement nul sous une autre, voire s&rsquo;inverser compl\u00e8tement sous une troisi\u00e8me condition. La pr\u00e9sence d&rsquo;une interaction crois\u00e9e invalide l&rsquo;interpr\u00e9tation substantielle des effets principaux pris isol\u00e9ment, car affirmer qu&rsquo;un facteur produit un effet global moyen occulte le fait empirique que son action est diam\u00e9tralement oppos\u00e9e selon le sous-groupe consid\u00e9r\u00e9.<\/p>\n<h3>1.3 Le mod\u00e8le lin\u00e9aire sous-jacent<\/h3>\n<p>Sur le plan formel, l&rsquo;ANOVA factorielle \u00e0 deux facteurs s&rsquo;enracine pleinement dans le cadre unifi\u00e9 du <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/online.stat.psu.edu\/stat500\/lesson\/10\">mod\u00e8le lin\u00e9aire g\u00e9n\u00e9ral<\/a>. Pour une observation individuelle donn\u00e9e, identifi\u00e9e par l&rsquo;indice <em>k<\/em>, appartenant au niveau <em>i<\/em> du premier facteur (comportant <em>I<\/em> modalit\u00e9s) et au niveau <em>j<\/em> du second facteur (comportant <em>J<\/em> modalit\u00e9s), l&rsquo;\u00e9quation fondamentale du mod\u00e8le s&rsquo;\u00e9nonce sous sa forme param\u00e9trique classique :<\/p>\n<p><em>Y<sub>ijk<\/sub> = \u03bc + \u03b1<sub>i<\/sub> + \u03b2<sub>j<\/sub> + (\u03b1\u03b2)<sub>ij<\/sub> + \u03b5<sub>ijk<\/sub><\/em><\/p>\n<p>Dans cette formulation, <em>Y<sub>ijk<\/sub><\/em> repr\u00e9sente la valeur observ\u00e9e de la variable continue d\u00e9pendante pour le <em>k<\/em>-i\u00e8me sujet au sein de la cellule (<em>i, j<\/em>). Le param\u00e8tre <em>\u03bc<\/em> d\u00e9signe la moyenne g\u00e9n\u00e9rale th\u00e9orique (ou grand mean) de l&rsquo;ensemble de la population statistique. Le terme <em>\u03b1<sub>i<\/sub><\/em> correspond \u00e0 l&rsquo;effet sp\u00e9cifique du niveau <em>i<\/em> du facteur A, mesur\u00e9 comme l&rsquo;\u00e9cart diff\u00e9rentiel par rapport \u00e0 la moyenne globale sous la contrainte d&rsquo;identifiabilit\u00e9 usuelle voulant que la somme des <em>\u03b1<sub>i<\/sub><\/em> soit \u00e9gale \u00e0 z\u00e9ro. De m\u00eame, <em>\u03b2<sub>j<\/sub><\/em> mat\u00e9rialise l&rsquo;effet propre au niveau <em>j<\/em> du facteur B, sous la contrainte \u00e9quivalente que la somme des <em>\u03b2<sub>j<\/sub><\/em> s&rsquo;annule. Le terme <em>(\u03b1\u03b2)<sub>ij<\/sub><\/em> repr\u00e9sente la composante d&rsquo;interaction non additive sp\u00e9cifique \u00e0 la cellule (<em>i, j<\/em>), mesurant la d\u00e9viation de la moyenne de cellule par rapport aux pr\u00e9dictions purement additives des effets principaux. Enfin, <em>\u03b5<sub>ijk<\/sub><\/em> incarne le terme d&rsquo;erreur r\u00e9siduelle al\u00e9atoire, postul\u00e9e ind\u00e9pendante et identiquement distribu\u00e9e selon une loi normale d&rsquo;esp\u00e9rance nulle et de variance constante.<\/p>\n<p>L&rsquo;estimation des param\u00e8tres repose sur la d\u00e9composition arithm\u00e9tique rigoureuse de la variance totale observ\u00e9e, exprim\u00e9e sous la forme de la somme des carr\u00e9s totale (Total Sum of Squares, SST). Cette variabilit\u00e9 globale est scind\u00e9e en deux macro-composantes orthogonales : la somme des carr\u00e9s inter-groupes ou du mod\u00e8le (Model Sum of Squares, SSM) et la somme des carr\u00e9s intra-groupes ou r\u00e9siduelle (Error Sum of Squares, SSE). La composante du mod\u00e8le est elle-m\u00eame subdivis\u00e9e pour isoler la contribution unique de chaque source de variation :<\/p>\n<p><em>SST = SSA + SSB + SSAB + SSE<\/em><\/p>\n<p>Chacune de ces sommes de carr\u00e9s est indissociable d&rsquo;un nombre sp\u00e9cifique de degr\u00e9s de libert\u00e9 (ddl), refl\u00e9tant le nombre de comparaisons ind\u00e9pendantes possibles au sein de la matrice factorielle. Si l&rsquo;on note <em>N<\/em> l&rsquo;effectif total d&rsquo;observations, la distribution des degr\u00e9s de libert\u00e9 s&rsquo;\u00e9tablit comme suit :<\/p>\n<ul>\n<li>Facteur A : <em>df<sub>A<\/sub> = I &#8211; 1<\/em><\/li>\n<li>Facteur B : <em>df<sub>B<\/sub> = J &#8211; 1<\/em><\/li>\n<li>Interaction AB : <em>df<sub>AB<\/sub> = (I &#8211; 1) \u00d7 (J &#8211; 1)<\/em><\/li>\n<li>Erreur r\u00e9siduelle : <em>df<sub>E<\/sub> = N &#8211; (I \u00d7 J)<\/em><\/li>\n<li>Total : <em>df<sub>Total<\/sub> = N &#8211; 1<\/em><\/li>\n<\/ul>\n<p>Le ratio de chaque somme des carr\u00e9s par son degr\u00e9 de libert\u00e9 respectif fournit les carr\u00e9s moyens (Mean Squares, MS). La statistique d&rsquo;inf\u00e9rence de Fisher-Snedecor (statistique <em>F<\/em>) est ensuite g\u00e9n\u00e9r\u00e9e en calculant le rapport entre le carr\u00e9 moyen d&rsquo;un effet donn\u00e9 et le carr\u00e9 moyen de l&rsquo;erreur r\u00e9siduelle, constituant le fondement du rejet probabiliste de l&rsquo;hypoth\u00e8se nulle.<\/p>\n<h2>2. Formulation des hypoth\u00e8ses statistiques et sp\u00e9cification du mod\u00e8le<\/h2>\n<h3>2.1 Hypoth\u00e8ses nulles et alternatives<\/h3>\n<p>L&rsquo;impl\u00e9mentation op\u00e9rationnelle de l&rsquo;ANOVA factorielle \u00e0 deux facteurs n\u00e9cessite l&rsquo;\u00e9nonciation probabiliste formelle de trois couples distincts d&rsquo;hypoth\u00e8ses statistiques. Chaque couple met en concurrence une hypoth\u00e8se nulle (<em>H<sub>0<\/sub><\/em>), stipulant l&rsquo;absence de tout effet syst\u00e9matique au sein de la population parente, et une hypoth\u00e8se alternative (<em>H<sub>1<\/sub><\/em>), postulant l&rsquo;existence d&rsquo;une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 significative des param\u00e8tres mod\u00e9lis\u00e9s. Contrairement aux tests unilat\u00e9raux orient\u00e9s, l&rsquo;ANOVA op\u00e8re un test omnibus global qui \u00e9value l&rsquo;\u00e9galit\u00e9 simultan\u00e9e d&rsquo;un ensemble de param\u00e8tres.<\/p>\n<p>Le premier ensemble d&rsquo;hypoth\u00e8ses cible l&rsquo;effet principal du facteur A. L&rsquo;hypoth\u00e8se nulle postule que l&rsquo;ensemble des effets diff\u00e9rentiels <em>\u03b1<sub>i<\/sub><\/em> est rigoureusement nul, ou de fa\u00e7on \u00e9quivalente, que les moyennes marginales associ\u00e9es \u00e0 chaque niveau de A sont parfaitement identiques entre elles :<\/p>\n<p><em>H<sub>0(A)<\/sub> : \u03bc<sub>1.<\/sub> = \u03bc<sub>2.<\/sub> = &#8230; = \u03bc<sub>I.<\/sub><\/em> versus <em>H<sub>1(A)<\/sub> : \u2203 (i, i&rsquo;) tel que \u03bc<sub>i.<\/sub> \u2260 \u03bc<sub>i&rsquo;.<\/sub><\/em><\/p>\n<p>Le deuxi\u00e8me ensemble concerne l&rsquo;effet principal du facteur B. L&rsquo;hypoth\u00e8se nulle formalise l&rsquo;invariance des moyennes marginales du facteur B \u00e0 travers l&rsquo;ensemble de ses <em>J<\/em> modalit\u00e9s :<\/p>\n<p><em>H<sub>0(B)<\/sub> : \u03bc<sub>.1<\/sub> = \u03bc<sub>.2<\/sub> = &#8230; = \u03bc<sub>.J<\/sub><\/em> versus <em>H<sub>1(B)<\/sub> : \u2203 (j, j&rsquo;) tel que \u03bc<sub>.j<\/sub> \u2260 \u03bc<sub>.j&rsquo;<\/sub><\/em><\/p>\n<p>Le troisi\u00e8me ensemble, souvent le plus d\u00e9terminant sur le plan de la validit\u00e9 th\u00e9orique, interroge la significativit\u00e9 du terme d&rsquo;interaction bidirectionnelle. L&rsquo;hypoth\u00e8se nulle d&rsquo;absence d&rsquo;interaction stipule que la diff\u00e9rence de moyennes entre deux niveaux quelconques du facteur A demeure rigoureusement invariante quel que soit le niveau consid\u00e9r\u00e9 du facteur B, ce qui se traduit par l&rsquo;annulation conjointe de tous les coefficients <em>(\u03b1\u03b2)<sub>ij<\/sub><\/em> :<\/p>\n<p><em>H<sub>0(AB)<\/sub> : (\u03b1\u03b2)<sub>ij<\/sub> = 0 pour tout couple (i, j)<\/em> versus <em>H<sub>1(AB)<\/sub> : \u2203 (i, j) tel que (\u03b1\u03b2)<sub>ij<\/sub> \u2260 0<\/em><\/p>\n<p>Les cons\u00e9quences inf\u00e9rentielles d\u00e9coulant du rejet de l&rsquo;une ou plusieurs de ces hypoth\u00e8ses conditionnent la totalit\u00e9 de la trajectoire analytique subs\u00e9quente. Si l&rsquo;hypoth\u00e8se nulle de l&rsquo;interaction <em>H<sub>0(AB)<\/sub><\/em> est maintenue (p &ge; 0,05), l&rsquo;analyste est fond\u00e9 \u00e0 interpr\u00e9ter directement les effets principaux <em>H<sub>0(A)<\/sub><\/em> et <em>H<sub>0(B)<\/sub><\/em>. En revanche, le rejet formel de <em>H<sub>0(AB)<\/sub><\/em> signifie que le mod\u00e8le additif est inad\u00e9quat pour d\u00e9crire la complexit\u00e9 des donn\u00e9es : l&rsquo;interpr\u00e9tation des effets principaux globaux devient m\u00e9thodologiquement subsidiaire, voire trompeuse, imposant une d\u00e9composition d\u00e9taill\u00e9e des effets simples au sein de chaque sous-groupe factoriel.<\/p>\n<h3>2.2 \u00c9quilibrage des plans exp\u00e9rimentaux<\/h3>\n<p>L&rsquo;\u00e9quilibre structurel d&rsquo;un plan factoriel d\u00e9pend de la r\u00e9partition num\u00e9rique des observations au sein des diff\u00e9rentes cellules exp\u00e9rimentales form\u00e9es par le croisement des facteurs. Un plan d&rsquo;exp\u00e9rience est qualifi\u00e9 d&rsquo;\u00e9quilibr\u00e9 (ou balanc\u00e9) lorsque l&rsquo;effectif d&rsquo;\u00e9chantillonnage est strictement identique dans chacune des <em>I \u00d7 J<\/em> cellules du dispositif (<em>n<sub>ij<\/sub> = n<\/em> pour tout <em>i<\/em> et tout <em>j<\/em>). Dans cette configuration id\u00e9ale, les facteurs A, B et leur terme d&rsquo;interaction b\u00e9n\u00e9ficient d&rsquo;une propri\u00e9t\u00e9 math\u00e9matique fondamentale : l&rsquo;orthogonalit\u00e9 parfaite. L&rsquo;orthogonalit\u00e9 garantit que les sommes des carr\u00e9s attribuables \u00e0 chaque source de variation sont totalement ind\u00e9pendantes les unes des autres. La variance expliqu\u00e9e par un facteur ne chevauche en aucune mani\u00e8re celle expliqu\u00e9e par un autre, assurant une d\u00e9composition univoque de la variance globale sans ambigu\u00eft\u00e9 de calcul.<\/p>\n<p>Dans la pratique de la recherche empirique, particuli\u00e8rement dans les \u00e9tudes cliniques, sociologiques ou \u00e9pid\u00e9miologiques, les plans d&rsquo;exp\u00e9rience deviennent fr\u00e9quemment d\u00e9s\u00e9quilibr\u00e9s. Cette distorsion r\u00e9sulte d&rsquo;abandons en cours d&rsquo;exp\u00e9rimentation (attrition des participants), de donn\u00e9es manquantes non al\u00e9atoires, ou de contraintes d&rsquo;\u00e9chantillonnage inh\u00e9rentes aux populations rares. D\u00e8s lors que les effectifs de cellule deviennent in\u00e9gaux, l&rsquo;orthogonalit\u00e9 est irr\u00e9versiblement rompue : les facteurs se trouvent corr\u00e9l\u00e9s entre eux. Cette colin\u00e9arit\u00e9 artificielle induit un ph\u00e9nom\u00e8ne de redondance o\u00f9 certaines fractions de la variance expliqu\u00e9e peuvent \u00eatre attribu\u00e9es indistinctement au facteur A, au facteur B, ou \u00e0 leur interaction conjointe.<\/p>\n<p>Face \u00e0 cette complexit\u00e9, la m\u00e9thode d&rsquo;estimation d\u00e9pend crucialement du type de d\u00e9composition des sommes des carr\u00e9s retenu dans le moteur d&rsquo;estimation statistique. La th\u00e9orie statistique distingue classiquement trois approches fondamentales, d\u00e9sign\u00e9es sous les vocables de sommes des carr\u00e9s de Type I, Type II et Type III :<\/p>\n<ul>\n<li><strong>Type I (Sommes des carr\u00e9s s\u00e9quentielles ou hi\u00e9rarchiques) :<\/strong> Dans cette modalit\u00e9, les facteurs sont introduits successivement dans le mod\u00e8le selon un ordre sp\u00e9cifi\u00e9 par l&rsquo;utilisateur. Le premier facteur incorpor\u00e9 capte toute la variance qu&rsquo;il peut expliquer ; le second facteur ne capte que la variance incr\u00e9mentale non expliqu\u00e9e par le premier, et le terme d&rsquo;interaction ne re\u00e7oit que la variance r\u00e9siduelle finale. Cette approche d\u00e9pend lourdement de l&rsquo;ordre d&rsquo;entr\u00e9e des variables, ce qui la rend inappropri\u00e9e pour tester des hypoth\u00e8ses factorielles sans pr\u00e9s\u00e9ance th\u00e9orique absolue.<\/li>\n<li><strong>Type II (Sommes des carr\u00e9s hi\u00e9rarchiques ajust\u00e9es) :<\/strong> Chaque effet principal est estim\u00e9 en \u00e9tant ajust\u00e9 pour l&rsquo;autre effet principal, mais en ignorant d\u00e9lib\u00e9r\u00e9ment le terme d&rsquo;interaction d&rsquo;ordre sup\u00e9rieur. Cette approche est hautement recommand\u00e9e par de nombreux statisticiens lorsque l&rsquo;interaction est pr\u00e9sum\u00e9e absente ou scientifiquement n\u00e9gligeable, car elle offre une puissance statistique accrue pour l&rsquo;\u00e9valuation des effets principaux.<\/li>\n<li><strong>Type III (Sommes des carr\u00e9s marginales ou partielles) :<\/strong> Dans cette m\u00e9thode, la somme des carr\u00e9s associ\u00e9e \u00e0 chaque effet est calcul\u00e9e en contr\u00f4lant rigoureusement pour l&rsquo;ensemble des autres termes pr\u00e9sents dans le mod\u00e8le, y compris les interactions d&rsquo;ordre sup\u00e9rieur. Chaque effet refl\u00e8te ainsi la variance unique et non partag\u00e9e qu&rsquo;il apporte \u00e0 la pr\u00e9diction de la r\u00e9ponse. C&rsquo;est l&rsquo;approche standard adopt\u00e9e par d\u00e9faut dans la grande majorit\u00e9 des progiciels statistiques contemporains pour les plans d\u00e9s\u00e9quilibr\u00e9s, car elle garantit que les tests des effets principaux restent invariants face aux variations d&rsquo;effectifs marginaux.<\/li>\n<\/ul>\n<p>Dans l&rsquo;environnement Stata, la commande native d&rsquo;analyse factorielle ex\u00e9cute par d\u00e9faut les sommes des carr\u00e9s partielles de Type III, garantissant une coh\u00e9rence optimale dans l&rsquo;\u00e9valuation des contrastes factoriels au sein de sous-groupes asym\u00e9triques.<\/p>\n<h2>3. V\u00e9rification pr\u00e9alable des postulats de l&rsquo;ANOVA sous Stata<\/h2>\n<h3>3.1 Ind\u00e9pendance des observations<\/h3>\n<p>L&rsquo;ind\u00e9pendance statistique des observations repr\u00e9sente le postulat le plus critique et le moins n\u00e9gociable de l&rsquo;analyse de variance factorielle. Cette condition fondamentale exige que la valeur prise par la variable d\u00e9pendante pour un individu donn\u00e9 ne soit en aucun cas influenc\u00e9e, conditionn\u00e9e ou li\u00e9e \u00e0 la valeur mesur\u00e9e chez un quelconque autre participant de l&rsquo;\u00e9chantillon. Math\u00e9matiquement, cela implique que la covariance entre les erreurs r\u00e9siduelles de deux observations distinctes soit strictement nulle : <em>Cov(\u03b5<sub>k<\/sub>, \u03b5<sub>k&rsquo;<\/sub>) = 0<\/em> pour tout <em>k \u2260 k&rsquo;<\/em>.<\/p>\n<p>Contrairement aux postulats de distributionnelle normale ou d&rsquo;\u00e9galit\u00e9 des variances, l&rsquo;ind\u00e9pendance ne peut \u00eatre \u00e9valu\u00e9e ou certifi\u00e9e par un test statistique formel post-hoc. Elle rel\u00e8ve exclusivement de la validit\u00e9 interne du dispositif exp\u00e9rimental et de la rigueur de l&rsquo;\u00e9chantillonnage. Le recours \u00e0 un \u00e9chantillonnage probabiliste al\u00e9atoire et \u00e0 une assignation randomis\u00e9e en double aveugle des sujets au sein des diff\u00e9rentes conditions exp\u00e9rimentales constitue la seule garantie m\u00e9thodologique assurant la rupture de tout lien de d\u00e9pendance syst\u00e9matique entre les unit\u00e9s d&rsquo;observation.<\/p>\n<p>La violation de cette hypoth\u00e8se survient fr\u00e9quemment lorsqu&rsquo;il existe des structures de regroupement naturel, telles que des patients trait\u00e9s par un m\u00eame clinicien, des \u00e9l\u00e8ves appartenant \u00e0 une m\u00eame classe scolaire, ou des mesures collect\u00e9es s\u00e9quentiellement dans le temps sans neutralisation de l&rsquo;autocorr\u00e9lation. Cette corr\u00e9lation intra-classe engendre une sous-estimation dramatique de la variance d&rsquo;\u00e9chantillonnage de l&rsquo;erreur r\u00e9siduelle, ce qui gonfle artificiellement la statistique de test <em>F<\/em> et d\u00e9multiplie le taux d&rsquo;erreur de type I (rejet intempestif de l&rsquo;hypoth\u00e8se nulle en l&rsquo;absence d&rsquo;effet r\u00e9el). Il convient en outre de distinguer scrupuleusement les facteurs inter-sujets (o\u00f9 chaque participant n&rsquo;est assign\u00e9 qu&rsquo;\u00e0 une seule et unique cellule) des plans \u00e0 mesures r\u00e9p\u00e9t\u00e9es ou intra-sujets, ces derniers exigeant des sp\u00e9cifications d&rsquo;ANOVA pour mesures r\u00e9p\u00e9t\u00e9es ou des mod\u00e8les lin\u00e9aires mixtes hi\u00e9rarchiques.<\/p>\n<h3>3.2 Normalit\u00e9 de la distribution des r\u00e9sidus<\/h3>\n<p>L&rsquo;inf\u00e9rence param\u00e9trique de l&rsquo;ANOVA repose sur l&rsquo;hypoth\u00e8se formelle selon laquelle les termes d&rsquo;erreur r\u00e9siduelle <em>\u03b5<sub>ijk<\/sub><\/em> suivent une loi normale au sein de chaque sous-population associ\u00e9e \u00e0 une cellule exp\u00e9rimentale. Il importe de souligner une m\u00e9prise conceptuelle r\u00e9currente dans la litt\u00e9rature appliqu\u00e9e : ce n&rsquo;est pas la variable d\u00e9pendante brute qui doit pr\u00e9senter une allure gaussienne \u00e0 l&rsquo;\u00e9chelle de l&rsquo;\u00e9chantillon agr\u00e9g\u00e9, mais bien la composante r\u00e9siduelle non expliqu\u00e9e par le mod\u00e8le apr\u00e8s soustraction des moyennes de cellules respectives. \u00c9valuer la normalit\u00e9 de la variable d\u00e9pendante sans la scinder par groupe constitue une erreur m\u00e9thodologique fr\u00e9quente, car la juxtaposition de plusieurs sous-groupes pr\u00e9sentant des moyennes distinctes engendre naturellement une distribution bimodale ou multimodale, parfaitement compatible avec des r\u00e9sidus locaux rigoureusement normaux.<\/p>\n<p>Dans l&rsquo;\u00e9cosyst\u00e8me Stata, l&rsquo;\u00e9valuation de cette normalit\u00e9 r\u00e9siduelle s&rsquo;op\u00e8re par une combinaison articul\u00e9e d&rsquo;outils visuels diagnostiques et de proc\u00e9dures de test formelles. L&rsquo;approche graphique repose principalement sur la construction de diagrammes quantile-quantile (Q-Q plots). Apr\u00e8s avoir ajust\u00e9 le mod\u00e8le factoriel et g\u00e9n\u00e9r\u00e9 les r\u00e9sidus standardis\u00e9s \u00e0 l&rsquo;aide de la commande <em>predict res, rstandard<\/em>, la commande Stata <em>qnorm res<\/em> permet de tracer les quantiles de la distribution empirique des r\u00e9sidus en regard des quantiles th\u00e9oriques issus d&rsquo;une loi normale standardis\u00e9e. Tout \u00e9cart syst\u00e9matique par rapport \u00e0 la diagonale de r\u00e9f\u00e9rence mat\u00e9rialise une distorsion morphologique de la distribution, qu&rsquo;il s&rsquo;agisse d&rsquo;une asym\u00e9trie marqu\u00e9e ou d&rsquo;un aplatissement anomal (queues lourdes ou l\u00e9g\u00e8res).<\/p>\n<p>Sur le plan des tests inf\u00e9rentiels quantitatifs, Stata offre plusieurs commandes sp\u00e9cialis\u00e9es. Le test de D&rsquo;Agostino, d&rsquo;asym\u00e9trie et de kurtosis, ex\u00e9cut\u00e9 via la commande <em>sktest res<\/em>, \u00e9value simultan\u00e9ment si le coefficient d&rsquo;asym\u00e9trie (skewness) et le coefficient d&rsquo;aplatissement (kurtosis) d\u00e9vient de mani\u00e8re significative des valeurs attendues pour une loi normale (respectivement 0 et 3). Conjointement, le test de Shapiro-Wilk, sollicit\u00e9 au moyen de la syntaxe <em>swilk res<\/em>, teste la conformit\u00e9 globale de la distribution r\u00e9siduelle pour des \u00e9chantillons de taille mod\u00e9r\u00e9e (N &lt; 2000). Une valeur p inf\u00e9rieure au seuil critique conventionnel de 0,05 conduit au rejet de l&rsquo;hypoth\u00e8se de normalit\u00e9.<\/p>\n<p>Il convient n\u00e9anmoins de temp\u00e9rer la s\u00e9v\u00e9rit\u00e9 de ces diagnostics \u00e0 la lumi\u00e8re du th\u00e9or\u00e8me central limite. L&rsquo;ANOVA factorielle d\u00e9montre une robustesse reconnue face \u00e0 des d\u00e9viations mod\u00e9r\u00e9es de la normalit\u00e9, \u00e0 condition que les distributions ne pr\u00e9sentent pas d&rsquo;asym\u00e9trie extr\u00eame et que les effectifs par cellule soient raisonnablement dimensionn\u00e9s (<em>n<\/em> &gt; 20 \u00e0 30 observations par groupe). En pr\u00e9sence d&rsquo;\u00e9chantillons volumineux, les tests formels comme celui de Shapiro-Wilk deviennent par ailleurs hyper-sensibles \u00e0 des d\u00e9viations minimes et sans cons\u00e9quence pratique sur le taux d&rsquo;erreur nominal, conf\u00e9rant au jugement graphique par diagramme quantile-quantile une primaut\u00e9 d\u00e9cisionnelle.<\/p>\n<h3>3.3 Homosc\u00e9dasticit\u00e9 des variances<\/h3>\n<p>Le postulat d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances, ou homosc\u00e9dasticit\u00e9, stipule que la variance de la variable d\u00e9pendante doit demeurer invariante \u00e0 travers l&rsquo;ensemble des cellules exp\u00e9rimentales d\u00e9finies par le croisement des facteurs : <em>\u03c3\u00b2<sub>11<\/sub> = \u03c3\u00b2<sub>12<\/sub> = &#8230; = \u03c3\u00b2<sub>IJ<\/sub> = \u03c3\u00b2<\/em>. La pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 compromet l&rsquo;exactitude de l&rsquo;estimation de l&rsquo;erreur r\u00e9siduelle globale (le d\u00e9nominateur du ratio <em>F<\/em>), introduisant des biais d&rsquo;inf\u00e9rence substantiels. Si des cellules \u00e0 faibles effectifs pr\u00e9sentent de fortes variances tandis que des cellules \u00e0 grands effectifs pr\u00e9sentent de faibles variances, le test <em>F<\/em> devient ind\u00fbment lib\u00e9ral, augmentant massivement les faux positifs. Dans la configuration sym\u00e9trique inverse, le test devient excessivement conservateur, \u00e9touffant la d\u00e9couverte d&rsquo;effets scientifiquement authentiques.<\/p>\n<p>L&rsquo;\u00e9valuation statistique de l&rsquo;homosc\u00e9dasticit\u00e9 s&rsquo;effectue sous Stata principalement par le biais du test de Levene, compl\u00e9t\u00e9 par ses variantes robustes propos\u00e9es par Brown et Forsythe. La commande Stata <em>robvar<\/em> constitue l&rsquo;outil de r\u00e9f\u00e9rence pour cette v\u00e9rification. Appliqu\u00e9e \u00e0 une variable de regroupement crois\u00e9e synth\u00e9tisant l&rsquo;appartenance factorielle, la commande calcule trois m\u00e9triques compl\u00e9mentaires : la statistique de Levene bas\u00e9e sur l&rsquo;\u00e9cart moyen par rapport \u00e0 la moyenne du groupe, la statistique de Brown-Forsythe bas\u00e9e sur la m\u00e9diane (hautement r\u00e9sistante aux valeurs aberrantes), et une variante fond\u00e9e sur la moyenne tronqu\u00e9e \u00e0 10 %. L&rsquo;obtention d&rsquo;une valeur de p sup\u00e9rieure \u00e0 0,05 valide le postulat d&rsquo;homog\u00e9n\u00e9it\u00e9.<\/p>\n<p>Historiquement, le test de Bartlett (disponible sous Stata via l&rsquo;option <em>oneway depvar group, bartlett<\/em>) a \u00e9t\u00e9 largement employ\u00e9 pour tester l&rsquo;\u00e9galit\u00e9 des variances. Toutefois, la litt\u00e9rature m\u00e9thodologique d\u00e9conseille son usage routinier en raison de son hypersensibilit\u00e9 extr\u00eame aux d\u00e9rogations \u00e0 la normalit\u00e9 : si les donn\u00e9es d\u00e9vient tant soit peu de la distribution gaussienne, le test de Bartlett rejettera erron\u00e9ment l&rsquo;hypoth\u00e8se d&rsquo;homosc\u00e9dasticit\u00e9, confondant asym\u00e9trie distributionnelle et disparit\u00e9 de variance.<\/p>\n<p>Lorsque l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 est formellement constat\u00e9e au sein du plan factoriel, plusieurs strat\u00e9gies de rem\u00e9diation s&rsquo;offrent \u00e0 l&rsquo;analyste. La premi\u00e8re consiste \u00e0 appliquer des transformations math\u00e9matiques stabilisatrices de variance sur la variable d\u00e9pendante, telles que la transformation logarithmique (particuli\u00e8rement adapt\u00e9e aux asym\u00e9tries positives), la transformation racine carr\u00e9e (pour les donn\u00e9es de comptage) ou la transformation r\u00e9ciproque. La seconde approche, plus moderne et respectueuse de l&rsquo;\u00e9chelle d&rsquo;origine des mesures, consiste \u00e0 abandonner l&rsquo;hypoth\u00e8se homosc\u00e9dastique au profit d&rsquo;une estimation lin\u00e9aire robuste des variances en utilisant des erreurs-types h\u00e9t\u00e9rosc\u00e9dastiques-consistantes (mod\u00e8le de Huber-White), ou \u00e0 recourir aux mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/stats.oarc.ucla.edu\/stata\/output\/glm\/\">GLM<\/a>).<\/p>\n<h2>4. Pr\u00e9paration de l&rsquo;environnement de travail et du jeu de donn\u00e9es dans Stata<\/h2>\n<h3>4.1 Chargement du jeu de donn\u00e9es d&rsquo;exemple<\/h3>\n<p>Afin d&rsquo;illustrer de mani\u00e8re concr\u00e8te et reproductible chaque \u00e9tape du processus computationnel, nous utiliserons un jeu de donn\u00e9es canonique issu directement des serveurs officiels de StataCorp : l&rsquo;archive <em>systolic<\/em>. Cette ressource p\u00e9dagogique et biom\u00e9dicale permet d&rsquo;\u00e9tudier l&rsquo;\u00e9volution de la pression art\u00e9rielle systolique sous l&rsquo;influence combin\u00e9e de diff\u00e9rents protocoles th\u00e9rapeutiques et du profil physiopathologique des sujets.<\/p>\n<p>Pour charger directement ce jeu de donn\u00e9es dans la m\u00e9moire vive de la session Stata active, il convient d&rsquo;ex\u00e9cuter la commande webuse dans la console ou au sein d&rsquo;un fichier script :<\/p>\n<p><em>webuse systolic, clear<\/em><\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/oneWayStata1-1.png\" alt=\"Systolic dataset in Stata example\" \/><figcaption>Systolic dataset in Stata example<\/figcaption><\/figure>\n<p>Le jeu de donn\u00e9es se compose de trois variables principales, totalisant 58 observations individuelles distribu\u00e9es au sein d&rsquo;une structure factorielle d\u00e9s\u00e9quilibr\u00e9e :<\/p>\n<ul>\n<li><strong>systolic :<\/strong> Variable continue d\u00e9pendante mesurant la variation relative (en millim\u00e8tres de mercure, mmHg) de la pression art\u00e9rielle systolique observ\u00e9e chez le patient apr\u00e8s administration de l&rsquo;intervention clinique.<\/li>\n<li><strong>drug :<\/strong> Premier facteur cat\u00e9goriel ind\u00e9pendant repr\u00e9sentant le type de traitement pharmacologique administr\u00e9. Ce facteur comporte 4 modalit\u00e9s distinctes (cod\u00e9es num\u00e9riquement de 1 \u00e0 4), o\u00f9 le niveau 1 correspond typiquement \u00e0 un groupe t\u00e9moin recevant un placebo, et les niveaux 2, 3 et 4 d\u00e9signent diff\u00e9rentes mol\u00e9cules hypotensives actives \u00e0 l&rsquo;\u00e9tude.<\/li>\n<li><strong>disease :<\/strong> Second facteur cat\u00e9goriel ind\u00e9pendant caract\u00e9risant la pathologie clinique sous-jacente du patient. Cette variable comporte 3 modalit\u00e9s distinctes (cod\u00e9es de 1 \u00e0 3), refl\u00e9tant diff\u00e9rents degr\u00e9s ou formes \u00e9tiologiques d&rsquo;hypertension art\u00e9rielle.<\/li>\n<\/ul>\n<p>Le croisement de ces deux facteurs g\u00e9n\u00e8re un plan factoriel complet 4 \u00d7 3, comprenant 12 cellules exp\u00e9rimentales distinctes. La particularit\u00e9 m\u00e9thodologique fondamentale de cet \u00e9chantillon de 58 observations r\u00e9side dans son asym\u00e9trie num\u00e9rique : le nombre de patients par cellule n&rsquo;est pas strictement identique (variant entre 4 et 6 sujets selon les croisements), ce qui en fait un cas d&rsquo;\u00e9cole id\u00e9al pour appr\u00e9hender le traitement rigoureux des plans factoriels d\u00e9s\u00e9quilibr\u00e9s et l&#8217;emploi des sommes des carr\u00e9s partielles de Type III.<\/p>\n<h3>4.2 Inspection des donn\u00e9es brutes<\/h3>\n<p>D\u00e8s le chargement de la base d&rsquo;investigation, une d\u00e9marche rigoureuse impose un audit minutieux des variables afin d&rsquo;identifier d&rsquo;\u00e9ventuelles anomalies de saisie, des valeurs manquantes non document\u00e9es ou des distorsions d&rsquo;\u00e9chantillonnage. Dans Stata, cette exploration pr\u00e9liminaire s&rsquo;initie par l&rsquo;inspection visuelle directe de la matrice de donn\u00e9es au moyen de la commande <em>browse<\/em>, laquelle ouvre le visualiseur sans risque d&rsquo;alt\u00e9ration accidentelle des enregistrements (contrairement \u00e0 la commande <em>edit<\/em>, qui autorise la modification directe des cellules et doit \u00eatre \u00e9vit\u00e9e en phase d&rsquo;audit pur).<\/p>\n<p>L&rsquo;interrogation synth\u00e9tique de l&rsquo;architecture du fichier s&rsquo;ex\u00e9cute \u00e0 l&rsquo;aide des commandes fondamentales <em>describe<\/em> et <em>codebook<\/em>. La commande <em>describe<\/em> renseigne sur le nombre total d&rsquo;observations, le nombre de colonnes, les types de stockage machine (byte, int, float, double) et l&rsquo;existence d&rsquo;\u00e9tiquettes pr\u00e9alables. Dans notre jeu de donn\u00e9es <em>systolic<\/em>, on observe que <em>systolic<\/em> est stock\u00e9e sous un format entier (<em>int<\/em>) tandis que <em>drug<\/em> et <em>disease<\/em> sont enregistr\u00e9es sous un format compact (<em>byte<\/em>), garantissant une gestion optimale de la m\u00e9moire vive.<\/p>\n<p>La commande <em>codebook<\/em> fournit quant \u00e0 elle une analyse structurelle plus fine en d\u00e9nombrant avec exactitude les valeurs manquantes (not\u00e9es sous forme de points d&rsquo;interrogation ou de <em>system missing<\/em> \u00ab\u00a0.\u00a0\u00bb dans Stata), les valeurs uniques prises par chaque champ, ainsi que les quantiles fondamentaux. L&rsquo;analyste doit imp\u00e9rativement s&rsquo;assurer qu&rsquo;aucune valeur n\u00e9gative aberrante n&rsquo;entache la variable d\u00e9pendante (sauf si une baisse de tension est explicitement mod\u00e9lis\u00e9e sous forme de diff\u00e9rentiel n\u00e9gatif) et qu&rsquo;aucune modalit\u00e9 num\u00e9rique impr\u00e9vue (par exemple un code 5 pour la variable <em>drug<\/em>) n&rsquo;a \u00e9t\u00e9 ins\u00e9r\u00e9e par erreur lors de l&rsquo;encodage clinique initial. Pour les 58 observations de l&rsquo;\u00e9chantillon <em>systolic<\/em>, aucune donn\u00e9e manquante n&rsquo;est recens\u00e9e, assurant la compl\u00e9tude de l&rsquo;\u00e9chantillon d&rsquo;analyse.<\/p>\n<h3>4.3 Labellisation et recodage des variables<\/h3>\n<p>La lisibilit\u00e9, la reproductibilit\u00e9 et la conformit\u00e9 aux normes scientifiques des sorties g\u00e9n\u00e9r\u00e9es par Stata reposent directement sur la qualit\u00e9 de l&rsquo;habillage s\u00e9mantique des variables. Des variables cat\u00e9gorielles d\u00e9pourvues de libell\u00e9s textuels explicites contraignent l&rsquo;analyste \u00e0 des allers-retours fastidieux avec le dictionnaire des codes et multiplient le risque d&rsquo;erreur d&rsquo;interpr\u00e9tation lors de la r\u00e9daction des manuscrits acad\u00e9miques.<\/p>\n<p>Le protocole de labellisation sous Stata s&rsquo;articule en deux \u00e9tapes interd\u00e9pendantes : la d\u00e9finition d&rsquo;un masque de libell\u00e9s via la directive <em>label define<\/em>, suivie de l&rsquo;application de ce masque aux modalit\u00e9s de la variable cible au moyen de la directive <em>label values<\/em>. Dans le cadre de notre jeu de donn\u00e9es, les instructions s&rsquo;\u00e9noncent avec une clart\u00e9 exemplaire :<\/p>\n<p><em>label define drug_lbl 1 \u00ab\u00a0Placebo\u00a0\u00bb 2 \u00ab\u00a0Mol\u00e9cule A\u00a0\u00bb 3 \u00ab\u00a0Mol\u00e9cule B\u00a0\u00bb 4 \u00ab\u00a0Mol\u00e9cule C\u00a0\u00bb<\/em><br \/>\n<em>label values drug drug_lbl<\/em><br \/>\n<em>label define disease_lbl 1 \u00ab\u00a0Stade I\u00a0\u00bb 2 \u00ab\u00a0Stade II\u00a0\u00bb 3 \u00ab\u00a0Stade III\u00a0\u00bb<\/em><br \/>\n<em>label values disease disease_lbl<\/em><\/p>\n<p>De surcro\u00eet, il est primordial d&rsquo;adjoindre une description substantielle aux variables elles-m\u00eames afin de documenter leur unit\u00e9 de mesure dans les m\u00e9tadonn\u00e9es du fichier via la commande <em>label variable<\/em> :<\/p>\n<p><em>label variable systolic \u00ab\u00a0Variation de pression art\u00e9rielle systolique (mmHg)\u00a0\u00bb<\/em><br \/>\n<em>label variable drug \u00ab\u00a0Protocole m\u00e9dicamenteux administr\u00e9\u00a0\u00bb<\/em><br \/>\n<em>label variable disease \u00ab\u00a0Stade de s\u00e9v\u00e9rit\u00e9 de la pathologie\u00a0\u00bb<\/em><\/p>\n<p>L&rsquo;ensemble de ces instructions pr\u00e9paratoires doit \u00eatre syst\u00e9matiquement consign\u00e9 dans un script Do-file structur\u00e9, comment\u00e9 et horodat\u00e9. Ce Do-file constitue la cl\u00e9 de vo\u00fbte de la reproductibilit\u00e9 scientifique, garantissant qu&rsquo;un tiers ou un auditeur puisse r\u00e9ex\u00e9cuter l&rsquo;int\u00e9gralit\u00e9 de la cha\u00eene de traitement depuis les micro-donn\u00e9es brutes jusqu&rsquo;aux estimations finales sans divergence de r\u00e9sultat.<\/p>\n<h2>5. Exploration descriptive et visualisation graphique des sous-groupes<\/h2>\n<h3>5.1 Statistiques descriptives par cellule exp\u00e9rimentale<\/h3>\n<p>Avant d&rsquo;engager les calculs d&rsquo;inf\u00e9rence statistique et de tester formellement les rapports de variance, l&rsquo;analyste se doit d&rsquo;explorer minutieusement les propri\u00e9t\u00e9s empiriques de chaque sous-groupe factoriel. Cette immersion descriptive vise \u00e0 appr\u00e9hender la distribution centrale des donn\u00e9es, \u00e0 \u00e9valuer la dispersion interne des mesures et \u00e0 d\u00e9tecter intuitivement les \u00e9ventuelles ruptures de tendance signalant des effets principaux ou des interactions saillantes.<\/p>\n<p>Stata met \u00e0 disposition plusieurs commandes d&rsquo;agr\u00e9gation d&rsquo;une grande efficacit\u00e9 op\u00e9rationnelle. La commande classique <em>tabulate<\/em>, enrichie de l&rsquo;option <em>summarize()<\/em>, permet de g\u00e9n\u00e9rer instantan\u00e9ment un tableau crois\u00e9 matriciel affichant pour chacune des 12 combinaisons factorielles la moyenne arithm\u00e9tique, l&rsquo;\u00e9cart-type ainsi que l&rsquo;effectif valide :<\/p>\n<p><em>tabulate drug disease, summarize(systolic)<\/em><\/p>\n<p>L&rsquo;examen minutieux de la matrice num\u00e9rique issue de cette commande r\u00e9v\u00e8le des \u00e9carts remarquables. Les moyennes de variation de pression systolique oscillent entre des valeurs modestes chez les sujets sous placebo (par exemple une moyenne autour de 6 mmHg au stade I de la maladie) et des r\u00e9ductions substantielles sous principes actifs avanc\u00e9s (atteignant des moyennes exc\u00e9dant 25 \u00e0 30 mmHg pour certains croisements th\u00e9rapeutiques). Cette analyse tabulaire permet \u00e9galement de v\u00e9rifier visuellement le d\u00e9s\u00e9quilibre d&rsquo;\u00e9chantillonnage : on constate que certaines cellules comportent 4 sujets tandis que d&rsquo;autres en comptent 5 ou 6. Pour des investigations descriptives encore plus d\u00e9taill\u00e9es int\u00e9grant les m\u00e9dianes, \u00e9carts interquartiles et erreurs-types, l&rsquo;instruction <em>table (drug) (disease), statistic(mean systolic) statistic(sd systolic) statistic(n systolic)<\/em> g\u00e9n\u00e8re des rapports tabulaires impeccablement calibr\u00e9s pour l&rsquo;exportation vers des traitements de texte scientifiques.<\/p>\n<h3>5.2 Visualisation graphique pr\u00e9liminaire<\/h3>\n<p>L&rsquo;exploration descriptive demeure incompl\u00e8te sans une projection visuelle des distributions empiriques. Le trac\u00e9 de diagrammes en bo\u00eete \u00e0 moustaches (boxplots) repr\u00e9sente la d\u00e9marche diagnostique la plus f\u00e9conde pour appr\u00e9cier simultan\u00e9ment la tendance centrale, la sym\u00e9trie de la distribution, la dispersion interquartile et l&rsquo;existence d&rsquo;observations aberrantes au sein de chaque sous-population d&rsquo;int\u00e9r\u00eat.<\/p>\n<p>Sous Stata, la commande <em>graph box<\/em> permet de segmenter graphiquement la variable d\u00e9pendante en combinant les deux facteurs ind\u00e9pendants via l&rsquo;option <em>over()<\/em> :<\/p>\n<p><em>graph box systolic, over(drug) over(disease) ytitle(\u00ab\u00a0Variation tensionnelle (mmHg)\u00a0\u00bb) title(\u00ab\u00a0Distribution de la pression systolique par groupe\u00a0\u00bb)<\/em><\/p>\n<p>Cette visualisation met en exergue l&rsquo;amplitude comparative des \u00e9carts interquartiles repr\u00e9sent\u00e9s par la hauteur des bo\u00eetes. Une disparit\u00e9 visuelle flagrante entre la hauteur des bo\u00eetes d&rsquo;une cellule \u00e0 l&rsquo;autre fournit un signal d&rsquo;alerte quant \u00e0 une potentielle violation du postulat d&rsquo;homosc\u00e9dasticit\u00e9. De m\u00eame, les observations se positionnant au-del\u00e0 de 1,5 fois l&rsquo;\u00e9cart interquartile (mat\u00e9rialis\u00e9es par des points isol\u00e9s \u00e0 l&rsquo;ext\u00e9rieur des moustaches) identifient imm\u00e9diatement des valeurs extr\u00eames susceptibles d&rsquo;exercer un levier excessif sur l&rsquo;estimation de la moyenne quadratique.<\/p>\n<p>Stata permet une personnalisation esth\u00e9tique pouss\u00e9e de ces figures. L&rsquo;analyste peut modifier les palettes de couleurs, int\u00e9grer des trames monochromes adapt\u00e9es \u00e0 l&rsquo;impression acad\u00e9mique, ajuster l&rsquo;orientation des \u00e9tiquettes d&rsquo;axes pour \u00e9viter les chevauchements textuels et exporter l&rsquo;illustration finale sous des formats vectoriels haute d\u00e9finition (tels que le format EPS ou PDF) via la syntaxe <em>graph export boxplot_initial.pdf, replace<\/em>.<\/p>\n<h2>6. Ex\u00e9cution syntaxique de l&rsquo;ANOVA factorielle avec Stata<\/h2>\n<h3>6.1 Syntaxe fondamentale de la commande anova<\/h3>\n<p>L&rsquo;ajustement du mod\u00e8le factoriel sous Stata repose sur la commande native <em>anova<\/em>, dont la syntaxe modulaire permet de sp\u00e9cifier des architectures d&rsquo;analyse de complexit\u00e9 arbitraire. Dans sa logique structurelle, la commande exige la d\u00e9signation imm\u00e9diate de la variable continue d\u00e9pendante, suivie s\u00e9quentiellement de la liste des facteurs ind\u00e9pendants cat\u00e9goriels et de leurs croisements multiplicatifs.<\/p>\n<p>Dans les versions contemporaines de Stata, la gestion des variables discr\u00e8tes est optimis\u00e9e par le syst\u00e8me d&rsquo;op\u00e9rateurs de variables factorielles (<em>factor variables<\/em>), caract\u00e9ris\u00e9 par l&rsquo;usage du pr\u00e9fixe <em>i.<\/em>. Pour expliciter un plan d&rsquo;exp\u00e9rience crois\u00e9 complet \u00e0 deux facteurs comportant les effets principaux et le terme d&rsquo;interaction, la syntaxe exhaustive s&rsquo;\u00e9nonce traditionnellement :<\/p>\n<p><em>anova systolic drug disease drug#disease<\/em><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/oneWayStata2-1.png\" alt=\"Two-way ANOVA example in Stata\" \/><figcaption>Two-way ANOVA example in Stata<\/figcaption><\/figure>\n<p>Dans cette \u00e9criture, le symbole di\u00e8se unique (<em>#<\/em>) ordonne explicitement au moteur de calcul de construire le terme d&rsquo;interaction entre les deux facteurs sans r\u00e9injecter leurs effets principaux au-del\u00e0 de ce qui a d\u00e9j\u00e0 \u00e9t\u00e9 libell\u00e9. Toutefois, Stata dispose d&rsquo;un raccourci d&rsquo;une grande \u00e9l\u00e9gance syntaxique : l&rsquo;op\u00e9rateur double di\u00e8se (<em>##<\/em>). Cet op\u00e9rateur factoriel instruit automatiquement le logiciel de mod\u00e9liser simultan\u00e9ment l&rsquo;ensemble des effets principaux individuels ainsi que toutes les interactions d&rsquo;ordre inf\u00e9rieur et sup\u00e9rieur qui en d\u00e9coulent. Par cons\u00e9quent, la commande suivante produit un r\u00e9sultat strictement identique \u00e0 la formulation pr\u00e9c\u00e9dente :<\/p>\n<p><em>anova systolic drug##disease<\/em><\/p>\n<p>L&rsquo;usage explicite de la notation factorielle <em>i.drug##i.disease<\/em> bien que facultatif au sein de la commande native <em>anova<\/em> (qui pr\u00e9sume par convention que tous les pr\u00e9dicteurs non continus sont cat\u00e9goriels) devient hautement recommand\u00e9 par mesure de rigueur m\u00e9thodologique, notamment pour garantir une interop\u00e9rabilit\u00e9 sans faille avec l&rsquo;ensemble des suites de commandes post-estimation comme <em>margins<\/em>.<\/p>\n<h3>6.2 Gestion des termes d&rsquo;interaction dans la commande<\/h3>\n<p>La distinction syntaxique entre un mod\u00e8le purement additif et un mod\u00e8le complet avec interaction refl\u00e8te deux postulats th\u00e9oriques fondamentalement divergents quant \u00e0 la nature du ph\u00e9nom\u00e8ne investigu\u00e9. Si l&rsquo;investigateur postule, sur la base de connaissances th\u00e9oriques pr\u00e9alables solides ou d&rsquo;objectifs de mod\u00e9lisation sp\u00e9cifiques, que l&rsquo;efficacit\u00e9 du m\u00e9dicament est strictement ind\u00e9pendante du stade pathologique du patient, il peut contraindre le mod\u00e8le \u00e0 un sch\u00e9ma purement additif en omettant d\u00e9lib\u00e9r\u00e9ment le produit crois\u00e9 :<\/p>\n<p><em>anova systolic drug disease<\/em><\/p>\n<p>Dans cette configuration additive, le mod\u00e8le alloue l&rsquo;int\u00e9gralit\u00e9 de la variance factorielle uniquement aux effets principaux de <em>drug<\/em> et de <em>disease<\/em>, rejetant toute variation conjointe r\u00e9siduelle directement dans le terme d&rsquo;erreur <em>SSE<\/em>. Il s&rsquo;agit l\u00e0 d&rsquo;une contrainte statistique forte : si une interaction r\u00e9elle existe au sein de la population mais se trouve bannie du mod\u00e8le, la variance d&rsquo;erreur r\u00e9siduelle sera artificiellement gonfl\u00e9e, d\u00e9gradant la sensibilit\u00e9 des tests <em>F<\/em> associ\u00e9s aux effets principaux.<\/p>\n<p>Un autre aspect technique crucial r\u00e9side dans le contr\u00f4le de la matrice de plan g\u00e9n\u00e9r\u00e9e en interne par Stata. L&rsquo;analyste doit garder \u00e0 l&rsquo;esprit que, par d\u00e9faut, la commande <em>anova<\/em> traite l&rsquo;ensemble des cellules selon une param\u00e9trisation en sommes des carr\u00e9s partielles de Type III. L&rsquo;ordre d&rsquo;\u00e9criture physique des variables dans la syntaxe (par exemple <em>drug##disease<\/em> versus <em>disease##drug<\/em>) n&rsquo;exerce donc strictement aucune influence sur les valeurs num\u00e9riques des statistiques <em>F<\/em>, des carr\u00e9s moyens ou des seuils de significativit\u00e9 obtenus pour chaque terme. Cette invariance garantit l&rsquo;objectivit\u00e9 des conclusions scientifiques face aux variations arbitraires d&rsquo;\u00e9criture de code.<\/p>\n<h2>7. Lecture et d\u00e9composition de la table de sortie ANOVA<\/h2>\n<h3>7.1 Analyse globale du mod\u00e8le statistique<\/h3>\n<p>\u00c0 la suite de l&rsquo;ex\u00e9cution de la commande factorielle, Stata compile et restitue l&rsquo;ensemble des param\u00e8tres au sein d&rsquo;une table d&rsquo;analyse de variance synth\u00e9tique. La lecture \u00e9clair\u00e9e de cette table requiert une progression ordonn\u00e9e, initi\u00e9e par l&rsquo;\u00e9valuation macroscopique de l&rsquo;ajustement global du mod\u00e8le lin\u00e9aire avant de descendre au niveau des effets analytiques individuels.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/twoWayStata1-1.png\" alt=\"Example of a two-way ANOVA in Stata\" \/><figcaption>Example of a two-way ANOVA in Stata<\/figcaption><\/figure>\n<p>La partie sup\u00e9rieure du tableau d&rsquo;affichage d\u00e9livre les m\u00e9triques r\u00e9sum\u00e9es du syst\u00e8me d&rsquo;\u00e9quations :<\/p>\n<ul>\n<li><strong>Number of obs :<\/strong> Indique la taille effective de l&rsquo;\u00e9chantillon retenu dans le calcul d&rsquo;inf\u00e9rence (dans notre exemple, N = 58).<\/li>\n<li><strong>R-squared (R\u00b2) :<\/strong> Repr\u00e9sente le coefficient de d\u00e9termination multiple global. Il quantifie la proportion de variance de la variable d\u00e9pendante expliqu\u00e9e conjointement par l&rsquo;ensemble des termes incorpor\u00e9s dans le mod\u00e8le (effets principaux et interaction). Une valeur de R\u00b2 \u00e9gale \u00e0 0,54 indique par exemple que 54 % de la variabilit\u00e9 globale de la pression art\u00e9rielle systolique est captur\u00e9e par la prise en compte du traitement, de la maladie et de leur synergie.<\/li>\n<li><strong>Adj R-squared (R\u00b2 ajust\u00e9) :<\/strong> Corrige le coefficient de d\u00e9termination standard pour tenir compte du nombre de param\u00e8tres estim\u00e9s en regard de la taille d&rsquo;\u00e9chantillon disponible. Contrairement au R\u00b2 qui s&rsquo;accro\u00eet m\u00e9caniquement \u00e0 chaque adjonction de param\u00e8tre, le R\u00b2 ajust\u00e9 p\u00e9nalise l&rsquo;introduction de termes non parcimonieux, fournissant une estimation non biais\u00e9e de la validit\u00e9 du mod\u00e8le dans la population parente.<\/li>\n<li><strong>Root MSE :<\/strong> Mat\u00e9rialise l&rsquo;\u00e9cart-type r\u00e9siduel estim\u00e9 (racine carr\u00e9e du carr\u00e9 moyen de l&rsquo;erreur, <em>\u221aMSE<\/em>). Il chiffre l&rsquo;erreur type moyenne de pr\u00e9diction exprim\u00e9e dans l&rsquo;unit\u00e9 native de la variable continue (ici, en mmHg).<\/li>\n<\/ul>\n<p>La premi\u00e8re ligne du corps du tableau, intitul\u00e9e <strong>Model<\/strong>, r\u00e9capitule la variance globale expliqu\u00e9e par l&rsquo;ensemble du dispositif exp\u00e9rimental. On y lit la somme des carr\u00e9s du mod\u00e8le (<em>SS<sub>Model<\/sub> = SSA + SSB + SSAB<\/em>), le total des degr\u00e9s de libert\u00e9 associ\u00e9s (<em>df<sub>Model<\/sub> = (I &#8211; 1) + (J &#8211; 1) + (I &#8211; 1)(J &#8211; 1)<\/em>), et le carr\u00e9 moyen r\u00e9sultant (<em>MS<sub>Model<\/sub> = SS<sub>Model<\/sub> \/ df<sub>Model<\/sub><\/em>). Le test <em>F<\/em> global associ\u00e9 confronte le mod\u00e8le factoriel complet \u00e0 un mod\u00e8le purement al\u00e9atoire n&rsquo;incluant que la constante (la moyenne g\u00e9n\u00e9rale). L&rsquo;obtention d&rsquo;une statistique <em>F<\/em> globale hautement significative avec une p-valeur tr\u00e8s inf\u00e9rieure au seuil critique conventionnel (<em>p &lt; 0,001<\/em>) autorise formellement l&rsquo;analyste \u00e0 rejeter l&rsquo;hypoth\u00e8se nulle g\u00e9n\u00e9rale et \u00e0 affirmer que le dispositif exp\u00e9rimental explique une part substantielle de la variabilit\u00e9 mesur\u00e9e.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/twoWayStata4.png\" alt=\"Two-Way ANOVA results in Stata\" \/><figcaption>Two-Way ANOVA results in Stata<\/figcaption><\/figure>\n<h3>7.2 Interpr\u00e9tation des effets principaux individuels<\/h3>\n<p>Une fois valid\u00e9e la significativit\u00e9 globale du mod\u00e8le, l&rsquo;attention se porte sur les lignes factorielles isolant la contribution sp\u00e9cifique de chaque effet principal : les lignes libell\u00e9es <strong>drug<\/strong> et <strong>disease<\/strong>.<\/p>\n<p>Pour le facteur <em>drug<\/em>, qui dispose de 4 modalit\u00e9s, les degr\u00e9s de libert\u00e9 associ\u00e9s s&rsquo;\u00e9l\u00e8vent logiquement \u00e0 <em>df = 4 &#8211; 1 = 3<\/em>. La table de Stata renseigne la somme des carr\u00e9s partielle de Type III (<em>Partial SS<\/em>), calcul\u00e9e en mesurant la perte d&rsquo;ajustement du mod\u00e8le si le terme <em>drug<\/em> venait \u00e0 \u00eatre retranch\u00e9 alors que <em>disease<\/em> et l&rsquo;interaction demeurent pr\u00e9sents. Le carr\u00e9 moyen (<em>MS<\/em>) est ensuite divis\u00e9 par le carr\u00e9 moyen r\u00e9siduel de la ligne <strong>Residual<\/strong> pour obtenir la statistique de Fisher :<\/p>\n<p><em>F<sub>drug<\/sub> = MS<sub>drug<\/sub> \/ MS<sub>Residual<\/sub><\/em><\/p>\n<p>Dans l&rsquo;\u00e9chantillon d&rsquo;exemple <em>systolic<\/em>, l&rsquo;examen de cette statistique r\u00e9v\u00e8le typiquement un rapport <em>F<\/em> de magnitude tr\u00e8s \u00e9lev\u00e9e (couramment <em>F &gt; 9<\/em> avec <em>df = (3, 46)<\/em>, affichant une p-valeur rigoureusement inf\u00e9rieure \u00e0 0,0001). Ce r\u00e9sultat conduit sans \u00e9quivoque au rejet cat\u00e9gorique de l&rsquo;hypoth\u00e8se nulle <em>H<sub>0(drug)<\/sub><\/em> : l&rsquo;administration des diff\u00e9rentes mol\u00e9cules m\u00e9dicamenteuses induit des variations diff\u00e9rentielles majeures de la pression systolique.<\/p>\n<p>Parall\u00e8lement, la ligne relative au facteur <em>disease<\/em> pr\u00e9sente <em>df = 3 &#8211; 1 = 2<\/em> degr\u00e9s de libert\u00e9. La statistique de Fisher correspondante, calcul\u00e9e selon la m\u00eame d\u00e9marche (<em>F<sub>disease<\/sub> = MS<sub>disease<\/sub> \/ MS<sub>Residual<\/sub><\/em>), produit \u00e9galement une valeur statistiquement significative (par exemple <em>p &lt; 0,05<\/em> ou <em>p &lt; 0,01<\/em>). Cela indique que la gravit\u00e9 pathologique sous-jacente des patients module syst\u00e9matiquement la r\u00e9activit\u00e9 tensionnelle globale, les patients aux stades les plus s\u00e9v\u00e8res pr\u00e9sentant des \u00e9carts de tension structurellement distincts de ceux enregistr\u00e9s aux stades pr\u00e9coces.<\/p>\n<h3>7.3 \u00c9valuation du terme d&rsquo;interaction<\/h3>\n<p>Le troisi\u00e8me palier d&rsquo;investigation statistique converge vers la ligne identifiant le terme combin\u00e9 d&rsquo;interaction : <strong>drug#disease<\/strong>. Ce terme r\u00e9clame <em>df = (4 &#8211; 1) \u00d7 (3 &#8211; 1) = 6<\/em> degr\u00e9s de libert\u00e9. La somme des carr\u00e9s partielle correspondante isole la variation qui ne peut \u00eatre imput\u00e9e ni \u00e0 l&rsquo;action moyenne lin\u00e9aire du traitement, ni \u00e0 la s\u00e9v\u00e9rit\u00e9 moyenne intrins\u00e8que de la pathologie, mais d\u00e9coule de leur synergie non lin\u00e9aire sp\u00e9cifique.<\/p>\n<p>L&rsquo;observation de la statistique de Fisher pour l&rsquo;interaction fournit le verdict m\u00e9thodologique d\u00e9terminant de l&rsquo;\u00e9tude. Dans le cadre des donn\u00e9es <em>systolic<\/em>, on obtient g\u00e9n\u00e9ralement une statistique <em>F<\/em> autour de 1,00 \u00e0 1,50 associ\u00e9e \u00e0 une p-valeur largement sup\u00e9rieure au seuil canonique de 0,05 (typiquement <em>p &gt; 0,30<\/em>). Cette absence de significativit\u00e9 statistique a des implications th\u00e9oriques et pratiques majeures : elle d\u00e9montre que les donn\u00e9es empiriques ne fournissent pas de preuves suffisantes pour rejeter l&rsquo;hypoth\u00e8se nulle d&rsquo;additivit\u00e9 <em>H<sub>0(AB)<\/sub><\/em>.<\/p>\n<p>D&rsquo;un point de vue m\u00e9thodologique, deux sc\u00e9narios antagonistes doivent \u00eatre scrupuleusement distingu\u00e9s selon la p-valeur de l&rsquo;interaction :<\/p>\n<ul>\n<li><strong>Sc\u00e9nario A (Interaction statistiquement non significative, p &ge; 0,05) :<\/strong> Les donn\u00e9es sont compatibles avec un mod\u00e8le strictement additif. L&rsquo;action pharmacologique des mol\u00e9cules m\u00e9dicamenteuses op\u00e8re de fa\u00e7on stable et constante quel que soit le stade clinique d&rsquo;atteinte pathologique du sujet. Dans ce contexte, l&rsquo;investigateur est pleinement autoris\u00e9 \u00e0 interpr\u00e9ter de mani\u00e8re directe, universelle et exhaustive les effets principaux identifi\u00e9s pr\u00e9c\u00e9demment.<\/li>\n<li><strong>Sc\u00e9nario B (Interaction statistiquement significative, p &lt; 0,05) :<\/strong> Le postulat d&rsquo;additivit\u00e9 s&rsquo;effondre. Il devient scientifiquement fallacieux de proclamer qu&rsquo;une mol\u00e9cule donn\u00e9e est \u00ab\u00a0la plus efficace\u00a0\u00bb de fa\u00e7on universelle : son efficacit\u00e9 d\u00e9pend intimement de la maladie du patient. D\u00e8s lors, l&rsquo;interpr\u00e9tation des effets principaux globaux doit \u00eatre suspendue ou subordonn\u00e9e \u00e0 une analyse exhaustive des effets simples (analyse conditionnelle par niveau de cellule).<\/li>\n<\/ul>\n<h2>8. Analyse approfondie de l&rsquo;interaction et moyennes marginales pr\u00e9dictives<\/h2>\n<h3>8.1 Calcul des moyennes marginales avec margins<\/h3>\n<p>L&rsquo;interpr\u00e9tation contemporaine des mod\u00e8les lin\u00e9aires sous Stata ne s&rsquo;arr\u00eate plus \u00e0 l&rsquo;inspection brute des tables de variances ; elle s&rsquo;articule autour du paradigme puissant des moyennes marginales pr\u00e9dictives (estimations lin\u00e9aires ajust\u00e9es), pilot\u00e9 par la suite logicielle post-estimation <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.stata.com\/help.cgi?margins\">margins<\/a>. Cet outil permet de traduire des param\u00e8tres math\u00e9matiques abstraits en grandeurs directement intelligibles exprim\u00e9es dans l&rsquo;\u00e9chelle d&rsquo;origine de la mesure clinique.<\/p>\n<p>Apr\u00e8s l&rsquo;estimation de l&rsquo;ANOVA par la syntaxe factorielle, l&rsquo;analyste ex\u00e9cute l&rsquo;instruction suivante pour calculer les moyennes ajust\u00e9es associ\u00e9es \u00e0 chaque cellule exp\u00e9rimentale :<\/p>\n<p><em>margins drug#disease<\/em><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/twoWayStata2.png\" alt=\"Two-way ANOVA in Stata\" \/><figcaption>Two-way ANOVA in Stata<\/figcaption><\/figure>\n<p>L&rsquo;ex\u00e9cution de cette instruction d\u00e9clenche la production d&rsquo;un tableau hautement structur\u00e9 pr\u00e9sentant, pour chacune des 12 combinaisons factorielles possibles :<\/p>\n<ul>\n<li>La moyenne marginale pr\u00e9dictive (Margin), repr\u00e9sentant la pr\u00e9diction attendue de variation systolique du mod\u00e8le pour ladite condition ;<\/li>\n<li>L&rsquo;erreur-type de pr\u00e9diction associ\u00e9e (Delta-method Standard Error), quantifiant la pr\u00e9cision probabiliste de l&rsquo;estimation lin\u00e9aire ;<\/li>\n<li>La statistique de Wald <em>z<\/em> ou le ratio <em>t<\/em> associ\u00e9 au test de nullit\u00e9 de la moyenne ;<\/li>\n<li>L&rsquo;intervalle de confiance bilat\u00e9ral \u00e0 95 % (95% CI) d\u00e9limitant la plage des valeurs populationnelles plausibles.<\/li>\n<\/ul>\n<p>Il est capital de diff\u00e9rencier conceptuellement les moyennes empiriques brutes calcul\u00e9es lors de la phase descriptive initiale des moyennes marginales pr\u00e9dictives fournies par <em>margins<\/em>. Dans un plan exp\u00e9rimental \u00e9quilibr\u00e9, ces deux ensembles de grandeurs sont rigoureusement identiques. En revanche, au sein d&rsquo;un plan d\u00e9s\u00e9quilibr\u00e9 comme celui du jeu de donn\u00e9es <em>systolic<\/em>, les moyennes marginales calcul\u00e9es par Stata sont des moyennes marginales ajust\u00e9es (souvent qualifi\u00e9es de <em>Least-Squares Means<\/em> ou moyennes des moindres carr\u00e9s). Elles corrigent math\u00e9matiquement l&rsquo;asym\u00e9trie num\u00e9rique des cellules en pond\u00e9rant uniform\u00e9ment chaque modalit\u00e9, restituant une pr\u00e9diction non biais\u00e9e par la surrepr\u00e9sentation ou la sous-repr\u00e9sentation fortuite de certains sous-groupes.<\/p>\n<h3>8.2 Graphique de profil d&rsquo;interaction avec marginsplot<\/h3>\n<p>Si la lecture des estimations num\u00e9riques d\u00e9livr\u00e9es par <em>margins<\/em> fournit l&rsquo;ensemble des m\u00e9triques d&rsquo;inf\u00e9rence, la perception cognitive des patrons d&rsquo;interaction r\u00e9clame imp\u00e9rativement une transcription visuelle. Stata r\u00e9pond \u00e0 cette exigence par l&rsquo;interm\u00e9diaire de la commande post-estimation coupl\u00e9e <em>marginsplot<\/em>, con\u00e7ue pour convertir imm\u00e9diatement les pr\u00e9dictions calcul\u00e9es en figures graphiques professionnelles de haute tenue.<\/p>\n<p>Imm\u00e9diatement apr\u00e8s avoir ex\u00e9cut\u00e9 <em>margins drug#disease<\/em>, l&rsquo;instruction suivante permet de tracer le profil d&rsquo;interaction complet :<\/p>\n<p><em>marginsplot, xdimension(disease) ytitle(\u00ab\u00a0Variation systolique pr\u00e9dite (mmHg)\u00a0\u00bb) xtitle(\u00ab\u00a0Stade de la maladie\u00a0\u00bb) title(\u00ab\u00a0Profil d&rsquo;interaction Traitement x Pathologie\u00a0\u00bb)<\/em><\/p>\n<p>Le graphique ainsi g\u00e9n\u00e9r\u00e9 projette sur l&rsquo;axe des abscisses (axe X) les trois stades pathologiques, positionne sur l&rsquo;axe des ordonn\u00e9es (axe Y) l&rsquo;ampleur de la r\u00e9duction tensionnelle, et repr\u00e9sente les quatre modalit\u00e9s th\u00e9rapeutiques sous forme de profils s\u00e9riels distincts accompagn\u00e9s de leurs intervalles de confiance \u00e0 95 %. L&rsquo;analyste peut intervertir l&rsquo;agencement conceptuel en sp\u00e9cifiant <em>xdimension(drug)<\/em> si l&rsquo;objectif d&rsquo;exposition requiert de visualiser l&rsquo;\u00e9volution continue des mol\u00e9cules selon les diff\u00e9rentes cohortes pathologiques.<\/p>\n<p>L&rsquo;analyse visuelle de cette figure s&rsquo;articule autour d&rsquo;un principe g\u00e9om\u00e9trique fondamental : le crit\u00e8re de parall\u00e9lisme. Si les profils reliant les moyennes de chaque traitement se d\u00e9ploient selon des trajectoires strictement parall\u00e8les, le comportement des mol\u00e9cules est invariant \u00e0 travers les stades de la maladie, ce qui mat\u00e9rialise graphiquement l&rsquo;absence d&rsquo;interaction. Inversement, toute rupture flagrante de parall\u00e9lisme (lignes divergentes, convergentes ou franchement s\u00e9cantes) illustre visuellement la pr\u00e9sence d&rsquo;un effet d&rsquo;interaction. La commande <em>marginsplot<\/em> autorise par ailleurs une personnalisation graphique compl\u00e8te via l&rsquo;ajustement des marqueurs g\u00e9om\u00e9triques, l&rsquo;adjonction de quadrillages de lecture et le param\u00e9trage pr\u00e9cis des l\u00e9gendes explicatives.<\/p>\n<h3>8.3 Analyse des effets simples<\/h3>\n<p>Lorsqu&rsquo;une interaction s&rsquo;av\u00e8re statistiquement significative (ou lorsqu&rsquo;une hypoth\u00e8se clinique hautement sp\u00e9cifique r\u00e9clame d&rsquo;analyser l&rsquo;efficacit\u00e9 diff\u00e9rentielle d&rsquo;un traitement \u00e0 un palier donn\u00e9 de s\u00e9v\u00e9rit\u00e9), l&rsquo;examen global des effets principaux devient inop\u00e9rant. L&rsquo;analyste se doit alors de d\u00e9composer l&rsquo;architecture globale en analysant les effets simples (<em>simple effects<\/em>). L&rsquo;analyse des effets simples consiste \u00e0 examiner l&rsquo;impact d&rsquo;un facteur donn\u00e9 en isolant et en fixant successivement chaque modalit\u00e9 de l&rsquo;autre facteur.<\/p>\n<p>Sous Stata, cette op\u00e9ration analytique fine se programme avec une remarquable fluidit\u00e9 gr\u00e2ce aux options de conditionnement de la commande <em>margins<\/em> et \u00e0 l&rsquo;outil <em>contrast<\/em>. Si l&rsquo;on souhaite examiner l&rsquo;effet du traitement pharmacologique de fa\u00e7on compartiment\u00e9e au sein de chaque stade pathologique sp\u00e9cifique, la syntaxe d\u00e9di\u00e9e s&rsquo;\u00e9nonce :<\/p>\n<p><em>margins disease, over(drug)<\/em><br \/>\n<em>contrast r.drug@disease<\/em><\/p>\n<p>L&rsquo;instruction <em>contrast r.drug@disease<\/em> ordonne \u00e0 Stata d&rsquo;\u00e9valuer le test omnibus de Fisher comparant les mol\u00e9cules th\u00e9rapeutiques isol\u00e9ment au stade I, puis de recalculer ce test de fa\u00e7on ind\u00e9pendante au stade II, et enfin au stade III. Les sorties statistiques fournissent une partition rigoureuse :<\/p>\n<ul>\n<li>Au stade I : le ratio <em>F<\/em> mesure si les traitements divergent chez les patients au stade d\u00e9butant ;<\/li>\n<li>Au stade II : un second ratio <em>F<\/em> quantifie l&rsquo;efficacit\u00e9 diff\u00e9rentielle sous atteinte mod\u00e9r\u00e9e ;<\/li>\n<li>Au stade III : un troisi\u00e8me ratio <em>F<\/em> \u00e9value l&rsquo;impact th\u00e9rapeutique face \u00e0 l&rsquo;affection s\u00e9v\u00e8re.<\/li>\n<\/ul>\n<p>L&rsquo;utilit\u00e9 translationnelle de cette d\u00e9composition est inestimable en pratique m\u00e9dicale. Elle permet de formuler des conclusions th\u00e9rapeutiques diff\u00e9renci\u00e9es : une mol\u00e9cule peut s&rsquo;av\u00e9rer rigoureusement inop\u00e9rante face au stade d\u00e9butant de la pathologie (effet simple non significatif), tout en d\u00e9ployant une sup\u00e9riorit\u00e9 pharmacologique majeure d\u00e8s lors que la maladie bascule dans son stade critique (effet simple hautement significatif). L&rsquo;analyse des effets simples \u00e9vite ainsi d&rsquo;\u00e9carter une mol\u00e9cule prometteuse dont l&rsquo;action cibl\u00e9e aurait \u00e9t\u00e9 dilu\u00e9e dans une moyenne globale d&rsquo;effet principal non significative.<\/p>\n<h2>9. Analyses post-hoc et ajustements pour comparaisons multiples<\/h2>\n<h3>9.1 Comparaisons par paires avec la commande pwcompare<\/h3>\n<p>Lorsqu&rsquo;un facteur cat\u00e9goriel comportant plus de deux modalit\u00e9s d\u00e9ploie un effet principal statistiquement significatif (ou dans le sillage de l&rsquo;identification d&rsquo;effets simples h\u00e9t\u00e9rog\u00e8nes), le test <em>F<\/em> initial informe le chercheur qu&rsquo;au moins deux moyennes divergent entre elles, mais demeure structurellement incapable d&rsquo;identifier pr\u00e9cis\u00e9ment quelles paires de conditions exp\u00e9rimentales sont responsables de cette divergence. Il devient alors indispensable d&rsquo;effectuer des comparaisons multiples par paires (tests post-hoc).<\/p>\n<p>Toutefois, la multiplication non r\u00e9gul\u00e9e de tests d&rsquo;hypoth\u00e8ses s\u00e9quentiels d\u00e9clenche un ph\u00e9nom\u00e8ne pernicieux bien connu des statisticiens : l&rsquo;inflation de l&rsquo;erreur de type I par famille (<em>familywise error rate<\/em>, FWER). Si l&rsquo;on effectue <em>k<\/em> tests statistiques ind\u00e9pendants au seuil nominal usuel de <em>\u03b1 = 0,05<\/em>, la probabilit\u00e9 globale de commettre au moins une fausse d\u00e9couverte (rejet erron\u00e9 d&rsquo;une hypoth\u00e8se nulle authentique) s&rsquo;\u00e9l\u00e8ve \u00e0 <em>1 &#8211; (1 &#8211; \u03b1)<sup>k<\/sup><\/em>. Pour un facteur comprenant 4 modalit\u00e9s (soit 6 comparaisons bilat\u00e9rales possibles), le risque cumul\u00e9 r\u00e9el d&rsquo;erreur de type I franchit le seuil inacceptable des 26 % sans mesure d&rsquo;ajustement appropri\u00e9e.<\/p>\n<p>Stata r\u00e9sout ce d\u00e9fi m\u00e9thodologique par l&rsquo;entremise de la commande post-estimation hautement polyvalente <em>pwcompare<\/em>, qui int\u00e8gre nativement l&rsquo;ensemble des algorithmes contemporains de contr\u00f4le de l&rsquo;erreur globale :<\/p>\n<p><em>pwcompare drug, mcompare(bonferroni) pveffects<\/em><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/twoWayStata3.png\" alt=\"Two-way ANOVA in Stata\" \/><figcaption>Two-way ANOVA in Stata<\/figcaption><\/figure>\n<p>L&rsquo;option <em>mcompare()<\/em> d\u00e9termine l&rsquo;approche d&rsquo;ajustement probabiliste souhait\u00e9e. Parmi les m\u00e9thodes rigoureusement impl\u00e9ment\u00e9es dans Stata, quatre proc\u00e9dures m\u00e9ritent une attention particuli\u00e8re :<\/p>\n<ul>\n<li><strong>Ajustement de Bonferroni :<\/strong> L&rsquo;approche la plus classique et conservatrice. Elle consiste \u00e0 multiplier la p-valeur brute par le nombre total de comparaisons par paires effectu\u00e9es (ou \u00e0 diviser le seuil alpha nominal par ce m\u00eame nombre). Bien que garantissant un contr\u00f4le strict du FWER, cette correction devient excessivement p\u00e9nalisante et s\u00e9v\u00e8re lorsque le nombre de modalit\u00e9s s&rsquo;accro\u00eet, augmentant dramatiquement le risque d&rsquo;erreur de type II (non-d\u00e9tection d&rsquo;effets r\u00e9els).<\/li>\n<li><strong>Ajustement de \u0160id\u00e1k :<\/strong> L\u00e9g\u00e8rement moins conservateur que Bonferroni, ce proc\u00e9d\u00e9 ajuste les probabilit\u00e9s selon la formule exacte <em>\u03b1<sub>ajust\u00e9<\/sub> = 1 &#8211; (1 &#8211; \u03b1)<sup>1\/k<\/sup><\/em> sous postulat d&rsquo;ind\u00e9pendance, offrant un gain marginal mais scientifiquement appr\u00e9ciable de puissance statistique.<\/li>\n<li><strong>Proc\u00e9dure de Tukey-Kramer (Tukey&rsquo;s HSD) :<\/strong> Sollicit\u00e9e via l&rsquo;argument <em>mcompare(tukey)<\/em>, cette m\u00e9thode fond\u00e9e sur la loi de l&rsquo;amplitude studentis\u00e9e (studentized range distribution) s&rsquo;impose comme le standard m\u00e9thodologique absolu pour les comparaisons factorielles de toutes les paires possibles. Elle maintient exactement le FWER au niveau nominal tout en pr\u00e9servant une puissance de d\u00e9tection nettement sup\u00e9rieure \u00e0 celle des m\u00e9thodes de Bonferroni ou \u0160id\u00e1k, et s&rsquo;adapte sans distorsion computationnelle aux plans d\u00e9s\u00e9quilibr\u00e9s.<\/li>\n<li><strong>Proc\u00e9dure de Scheff\u00e9 :<\/strong> Activ\u00e9e par <em>mcompare(scheffe)<\/em>, cette approche est la plus robuste et la plus conservatrice de toutes. Elle prot\u00e8ge le taux d&rsquo;erreur non seulement pour les simples comparaisons par paires, mais pour l&rsquo;infinit\u00e9 de tous les contrastes lin\u00e9aires complexes possibles imaginables \u00e0 travers le plan exp\u00e9rimental.<\/li>\n<\/ul>\n<p>L&rsquo;adjonction de l&rsquo;option <em>pveffects<\/em> ordonne \u00e0 Stata de consigner dans le tableau de sortie la diff\u00e9rence brute entre les deux moyennes compar\u00e9es, son erreur-type estim\u00e9e, la statistique de test associ\u00e9e, la p-valeur corrig\u00e9e finale ainsi que l&rsquo;intervalle de confiance r\u00e9ajust\u00e9 pour le protocole de comparaison multiple.<\/p>\n<h3>9.2 Contrastes planifi\u00e9s a priori<\/h3>\n<p>Si les tests post-hoc s&rsquo;apparentent \u00e0 une d\u00e9marche exploratoire exhaustive conduite <em>a posteriori<\/em> sans hypoth\u00e8se directionnelle initiale, les contrastes planifi\u00e9s rel\u00e8vent d&rsquo;une d\u00e9marche scientifique d\u00e9ductive et confirmatoire sup\u00e9rieure. D\u00e9finis en amont de la collecte des donn\u00e9es sur la base de fondements th\u00e9oriques \u00e9prouv\u00e9s, les contrastes planifi\u00e9s permettent de tester des hypoth\u00e8ses de recherche hautement sp\u00e9cifiques, cibl\u00e9es et focalis\u00e9es, tout en maximisant la puissance statistique disponible.<\/p>\n<p>Dans l&rsquo;environnement Stata, la formulation des contrastes s&rsquo;ex\u00e9cute avec une remarquable souplesse logique via la commande g\u00e9n\u00e9rale <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.stata.com\/help.cgi?contrast\">contrast<\/a>. Supposons qu&rsquo;au sein de notre \u00e9tude clinique sur la pression art\u00e9rielle, l&rsquo;interrogation th\u00e9orique centrale ne vise pas \u00e0 comparer indistinctement toutes les mol\u00e9cules entre elles, mais consiste sp\u00e9cifiquement \u00e0 confronter le groupe t\u00e9moin recevant le placebo (modalit\u00e9 1) \u00e0 l&rsquo;ensemble agr\u00e9g\u00e9 des patients trait\u00e9s par n&rsquo;importe laquelle des trois mol\u00e9cules pharmacologiques actives (modalit\u00e9s 2, 3 et 4 combin\u00e9es). Ce contraste orthogonal d&rsquo;hypoth\u00e8se se formalise par des coefficients de pond\u00e9ration respectant la contrainte arithm\u00e9tique fondamentale d&rsquo;annulation de leur somme :<\/p>\n<p>Poids : <em>[-3, +1, +1, +1]<\/em><\/p>\n<p>Sous Stata, cette hypoth\u00e8se planifi\u00e9e s&rsquo;\u00e9value au moyen de la syntaxe de contraste personnalis\u00e9e suivante :<\/p>\n<p><em>contrast {drug -3 1 1 1}<\/em><\/p>\n<p>Le moteur d&rsquo;estimation g\u00e9n\u00e8re alors un test d&rsquo;hypoth\u00e8se cibl\u00e9 produisant une somme des carr\u00e9s sp\u00e9cifique \u00e0 1 degr\u00e9 de libert\u00e9 et le ratio de Fisher correspondant. Cette approche cibl\u00e9e \u00e9vite de dissiper la puissance statistique dans des corrections p\u00e9nalisantes pour comparaisons multiples inutiles, offrant une sensibilit\u00e9 d&rsquo;inf\u00e9rence tr\u00e8s sup\u00e9rieure face aux questions fondamentales pos\u00e9es par le protocole de recherche.<\/p>\n<h2>10. Quantification de la taille d&rsquo;effet et calcul de puissance statistique<\/h2>\n<h3>10.1 Mesures de la taille d&rsquo;effet : Eta-carr\u00e9 et Eta-carr\u00e9 partiel<\/h3>\n<p>Au sein de la recherche contemporaine, l&rsquo;interpr\u00e9tation scientifique exclusive fond\u00e9e sur la p-valeur fait l&rsquo;objet de vives r\u00e9serves \u00e9pist\u00e9mologiques. Une p-valeur hautement significative atteste uniquement qu&rsquo;une divergence observ\u00e9e est tr\u00e8s improbable sous l&rsquo;hypoth\u00e8se nulle, mais ne renseigne en aucune mani\u00e8re sur l&rsquo;amplitude clinique ou th\u00e9orique du ph\u00e9nom\u00e8ne. L&rsquo;association d&rsquo;un indice standardis\u00e9 de taille d&rsquo;effet (<em>effect size<\/em>) est devenue une obligation d\u00e9ontologique et \u00e9ditoriale dans l&rsquo;ensemble des revues scientifiques de premier rang.<\/p>\n<p>Dans le domaine de l&rsquo;analyse de variance factorielle, l&rsquo;indice d&rsquo;estimation de taille d&rsquo;effet le plus universellement rapport\u00e9 est l&rsquo;Eta-carr\u00e9 partiel (<em>partial \u03b7\u00b2<\/em> ou <em>\u03b7<sub>p<\/sub>\u00b2<\/em>). Math\u00e9matiquement, l&rsquo;Eta-carr\u00e9 partiel repr\u00e9sente la proportion de variance de la variable d\u00e9pendante qui est expliqu\u00e9e par un facteur d\u00e9termin\u00e9, une fois que l&rsquo;on a extrait et contr\u00f4l\u00e9 la variance captur\u00e9e par l&rsquo;ensemble des autres termes et pr\u00e9dicteurs du mod\u00e8le :<\/p>\n<p><em>\u03b7<sub>p<\/sub>\u00b2 = SS<sub>Effet<\/sub> \/ (SS<sub>Effet<\/sub> + SS<sub>Erreur<\/sub>)<\/em><\/p>\n<p>Dans Stata, le calcul automatis\u00e9 de ces m\u00e9triques standardis\u00e9es apr\u00e8s l&rsquo;ex\u00e9cution de la commande <em>anova<\/em> s&rsquo;effectue imm\u00e9diatement \u00e0 l&rsquo;aide de l&rsquo;instruction post-estimation sp\u00e9cialis\u00e9e <em>estat esize<\/em> :<\/p>\n<p><em>estat esize<\/em><\/p>\n<p>Stata g\u00e9n\u00e8re instantan\u00e9ment un tableau exhaustif d\u00e9taillant pour chaque terme factoriel (les effets principaux de chaque variable et leur interaction conjointe) :<\/p>\n<ul>\n<li>La valeur ponctuelle de l&rsquo;Eta-carr\u00e9 partiel (<em>Partial Eta-Squared<\/em>) ;<\/li>\n<li>L&rsquo;intervalle de confiance exact \u00e0 95 % encadrant cet indice, calcul\u00e9 sur la base de la distribution non centrale de Fisher.<\/li>\n<\/ul>\n<p>Il importe toutefois de maintenir une vigilance conceptuelle rigoureuse quant aux propri\u00e9t\u00e9s d&rsquo;\u00e9chantillonnage de cet estimateur. L&rsquo;Eta-carr\u00e9 partiel pr\u00e9sente un biais statistique positif bien document\u00e9 : il tend \u00e0 surestimer l&rsquo;amplitude r\u00e9elle de la taille d&rsquo;effet dans la population parente, particuli\u00e8rement face \u00e0 des \u00e9chantillons de taille modeste. C&rsquo;est pourquoi de nombreux m\u00e9thodologistes pr\u00e9conisent conjointement le calcul de l&rsquo;Om\u00e9ga-carr\u00e9 partiel (<em>partial \u03c9\u00b2<\/em>), obtenu dans Stata en adjoignant l&rsquo;option <em>omega<\/em> \u00e0 la commande : <em>estat esize, omega<\/em>. L&rsquo;Om\u00e9ga-carr\u00e9 int\u00e8gre dans sa formulation le carr\u00e9 moyen de l&rsquo;erreur r\u00e9siduelle afin de corriger le biais d&rsquo;\u00e9chantillonnage, fournissant une estimation populationnelle remarquablement plus conservatrice et exacte.<\/p>\n<p>Pour guider l&rsquo;interpr\u00e9tation substantielle de ces indices, les conventions heuristiques canoniques \u00e9tablies par Jacob Cohen dans ses travaux fondateurs continuent de faire autorit\u00e9 dans la litt\u00e9rature en sciences humaines et biom\u00e9dicales :<\/p>\n<ul>\n<li><strong>Effet faible (small effect) :<\/strong> <em>\u03b7<sub>p<\/sub>\u00b2 &asymp; 0,01<\/em> (le facteur capture approximativement 1 % de la variance r\u00e9siduelle disponible) ;<\/li>\n<li><strong>Effet mod\u00e9r\u00e9 (medium effect) :<\/strong> <em>\u03b7<sub>p<\/sub>\u00b2 &asymp; 0,06<\/em> (le facteur capture environ 6 % de la variance) ;<\/li>\n<li><strong>Effet fort (large effect) :<\/strong> <em>\u03b7<sub>p<\/sub>\u00b2 &ge; 0,14<\/em> (le facteur gouverne 14 % ou plus de la variance r\u00e9siduelle).<\/li>\n<\/ul>\n<h3>10.2 \u00c9valuation de la puissance statistique a posteriori<\/h3>\n<p>La puissance statistique, d\u00e9sign\u00e9e par <em>1 &#8211; \u03b2<\/em>, mat\u00e9rialise la probabilit\u00e9 math\u00e9matique qu&rsquo;un test rejette \u00e0 juste titre une hypoth\u00e8se nulle qui est effectivement fausse au sein de la population. En d&rsquo;autres termes, elle chiffre la capacit\u00e9 op\u00e9rationnelle du protocole exp\u00e9rimental \u00e0 d\u00e9tecter un effet authentique de taille d\u00e9termin\u00e9e. Sous Stata, la suite logicielle d\u00e9di\u00e9e \u00e0 ces calculs s&rsquo;articule autour de la commande g\u00e9n\u00e9rique <em>power<\/em>.<\/p>\n<p>Il est cependant imp\u00e9ratif d&rsquo;\u00e9mettre une mise en garde m\u00e9thodologique fondamentale contre la pratique largement d\u00e9cri\u00e9e du calcul de puissance a posteriori (dite <em>observed power<\/em> ou <em>post-hoc power<\/em>), qui consiste \u00e0 recalculer la puissance d&rsquo;une exp\u00e9rience \u00e0 partir des tailles d&rsquo;effet et des p-valeurs directement observ\u00e9es dans l&rsquo;\u00e9chantillon d&rsquo;analyse. La litt\u00e9rature statistique moderne a formellement d\u00e9montr\u00e9 que la puissance observ\u00e9e post-hoc n&rsquo;est qu&rsquo;une transformation math\u00e9matique circulaire et triviale de la p-valeur : un test affichant une p-valeur non significative produira m\u00e9caniquement une faible puissance observ\u00e9e, sans apporter la moindre information scientifique additionnelle.<\/p>\n<p>En revanche, l&rsquo;utilisation \u00e9clair\u00e9e de la commande <em>power<\/em> trouve tout son sens scientifique dans une d\u00e9marche prospective ou r\u00e9trospectivement projective, visant \u00e0 dimensionner rigoureusement les futurs protocoles de recherche empirique sur la base des donn\u00e9es pr\u00e9liminaires recueillies. En observant par exemple que l&rsquo;effet principal d&rsquo;un nouveau traitement hypotensif s&rsquo;\u00e9tablit \u00e0 un <em>\u03b7<sub>p<\/sub>\u00b2<\/em> de 0,08 dans notre \u00e9tude pilote, l&rsquo;analyste utilise Stata pour d\u00e9terminer la taille d&rsquo;\u00e9chantillon requise afin de reproduire cette d\u00e9couverte avec une puissance statistique robuste conventionnelle fix\u00e9e \u00e0 80 % ou 90 % sous un risque alpha de 5 % :<\/p>\n<p><em>power oneway, k(4) power(0.80) alpha(0.05) varmeans(&#8230;)<\/em><\/p>\n<p>Cette d\u00e9marche de planification prot\u00e8ge le chercheur contre la mise en \u0153uvre de protocoles sous-dimensionn\u00e9s (sous-puissants), dont les r\u00e9sultats n\u00e9gatifs demeureraient syst\u00e9matiquement ind\u00e9cidables, oscillant entre l&rsquo;absence effective d&rsquo;impact pharmacologique et l&rsquo;incapacit\u00e9 statistique du dispositif exp\u00e9rimental \u00e0 capter le signal au milieu du bruit de fond de l&rsquo;erreur r\u00e9siduelle.<\/p>\n<h2>11. Diagnostics post-estimation et mod\u00e9lisation alternative sous Stata<\/h2>\n<h3>11.1 Extraction et examen visuel des r\u00e9sidus<\/h3>\n<p>La validation formelle d&rsquo;un mod\u00e8le factoriel ne saurait s&rsquo;achever sans un contr\u00f4le approfondi du comportement individuel des observations au travers des diagnostics de r\u00e9gression post-estimation. La premi\u00e8re \u00e9tape consiste \u00e0 extraire physiquement les r\u00e9sidus et les valeurs ajust\u00e9es (<em>fitted values<\/em>) de la m\u00e9moire de Stata pour les incorporer sous forme de nouvelles colonnes dans la base active \u00e0 l&rsquo;aide de la commande <em>predict<\/em> :<\/p>\n<p><em>predict yhat, xb<\/em><br \/>\n<em>predict rstud, rstudent<\/em><\/p>\n<p>Dans ces lignes de commande, l&rsquo;option <em>xb<\/em> ordonne le calcul des valeurs pr\u00e9dites moyennes pour chaque sujet selon sa cellule factorielle, tandis que l&rsquo;option <em>rstudent<\/em> extrait les r\u00e9sidus studentis\u00e9s sans suppression (ou r\u00e9sidus jackknife). Les r\u00e9sidus studentis\u00e9s pr\u00e9sentent l&rsquo;avantage d\u00e9cisif d&rsquo;\u00eatre standardis\u00e9s par une estimation de l&rsquo;\u00e9cart-type r\u00e9siduel excluant l&rsquo;observation elle-m\u00eame, ce qui stabilise leur variance unitaire et permet une identification rigoureuse des valeurs aberrantes selon une m\u00e9trique standardis\u00e9e.<\/p>\n<p>L&rsquo;outil d&rsquo;investigation graphique post-estimation central s&rsquo;ex\u00e9cute ensuite via la commande <em>rvfplot<\/em> (Residuals-versus-Fitted plot) :<\/p>\n<p><em>rvfplot, yline(0) mlabel(drug) title(\u00ab\u00a0Graphique des r\u00e9sidus studentis\u00e9s vs valeurs pr\u00e9dites\u00a0\u00bb)<\/em><\/p>\n<p>Ce graphique projette les r\u00e9sidus studentis\u00e9s le long d&rsquo;un axe vertical centr\u00e9 sur z\u00e9ro, en regard des valeurs pr\u00e9dites par le mod\u00e8le factoriel d\u00e9ploy\u00e9es sur l&rsquo;axe horizontal. L&rsquo;analyse s\u00e9miologique de ce nuage de points permet de valider simultan\u00e9ment plusieurs postulats :<\/p>\n<ul>\n<li><strong>Contr\u00f4le de la lin\u00e9arit\u00e9 et centrage :<\/strong> Les r\u00e9sidus doivent se r\u00e9partir de fa\u00e7on \u00e9quilibr\u00e9e et al\u00e9atoire de part et d&rsquo;autre de la ligne horizontale nulle, sans d\u00e9rive ascendante ou descendante syst\u00e9matique ;<\/li>\n<li><strong>Contr\u00f4le de l&rsquo;homosc\u00e9dasticit\u00e9 :<\/strong> L&rsquo;amplitude de dispersion verticale des points ne doit montrer aucune structure en entonnoir (cone-shaped pattern) ou d&rsquo;\u00e9vasement marqu\u00e9 d&rsquo;une valeur ajust\u00e9e \u00e0 l&rsquo;autre ;<\/li>\n<li><strong>D\u00e9tection des valeurs discordantes :<\/strong> Tout point franchissant les seuils critiques horizontaux de <em>&plusmn; 2<\/em> ou <em>&plusmn; 3<\/em> \u00e9carts-types studentis\u00e9s identifie un cas atypique majeur r\u00e9clamant une v\u00e9rification imm\u00e9diate de la transcription clinique d&rsquo;origine.<\/li>\n<\/ul>\n<p>Pour compl\u00e9ter cet audit de diagnostic, l&rsquo;\u00e9valuation des observations influentes et \u00e0 fort effet de levier (leverage) peut \u00eatre quantifi\u00e9e via le calcul des distances de Cook (<em>predict dcook, cooksd<\/em>). Une distance de Cook exc\u00e9dant la valeur seuil usuelle de <em>4 \/ N<\/em> ou s&rsquo;approchant de 1 signale une observation isol\u00e9e dont la pr\u00e9sence modifie substantiellement \u00e0 elle seule l&rsquo;estimation des coefficients du mod\u00e8le factoriel.<\/p>\n<h3>11.2 Alternatives en cas de violation des postulats<\/h3>\n<p>Lorsque les diagnostics r\u00e9v\u00e8lent des violations persistantes et irr\u00e9ductibles des postulats fondamentaux de l&rsquo;ANOVA param\u00e9trique (notamment une h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 s\u00e9v\u00e8re r\u00e9fractaire aux transformations math\u00e9matiques ou une non-normalit\u00e9 r\u00e9siduelle imputable \u00e0 de petits effectifs asym\u00e9triques), la rigueur scientifique impose de d\u00e9laisser la commande classique <em>anova<\/em> au profit d&rsquo;architectures d&rsquo;estimation alternatives plus robustes.<\/p>\n<p>La premi\u00e8re solution de rem\u00e9diation, remarquablement accessible sous Stata, repose sur l&rsquo;\u00e9quivalence formelle unissant l&rsquo;ANOVA factorielle et la r\u00e9gression lin\u00e9aire multiple via le mod\u00e8le lin\u00e9aire g\u00e9n\u00e9ral. Il est parfaitement possible d&rsquo;estimer exactement la m\u00eame structure factorielle en sollicitant la commande <em>regress<\/em> coupl\u00e9e \u00e0 l&rsquo;option d&rsquo;estimation robuste des variances (estimateur en sandwich de Huber-White) :<\/p>\n<p><em>regress systolic i.drug##i.disease, vce(robust)<\/em><br \/>\n<em>testparm i.drug#i.disease<\/em><\/p>\n<p>L&rsquo;utilisation de l&rsquo;option <em>vce(robust)<\/em> lib\u00e8re l&rsquo;inf\u00e9rence statistique de l&rsquo;obligation d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances en recalculant les erreurs-types de chaque contraste de fa\u00e7on consistante face \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9. Les tests d&rsquo;effets principaux et d&rsquo;interaction sont ensuite restitu\u00e9s sous forme de tests de Wald robustes via la commande post-estimation <em>testparm<\/em>, pr\u00e9servant la validit\u00e9 des seuils de significativit\u00e9 nominale m\u00eame en pr\u00e9sence de fortes disparit\u00e9s de dispersion inter-cellules.<\/p>\n<p>La deuxi\u00e8me perspective de recours s&rsquo;inscrit dans le champ des Mod\u00e8les Lin\u00e9aires G\u00e9n\u00e9ralis\u00e9s (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/stats.oarc.ucla.edu\/stata\/dae\/poisson-regression\/\">GLM<\/a>). Si la variable d\u00e9pendante mesur\u00e9e s&rsquo;av\u00e8re structurellement tronqu\u00e9e, strictement positive, ou distribu\u00e9e selon des dynamiques hautement asym\u00e9triques (telles que des dur\u00e9es de survie, des co\u00fbts m\u00e9dicaux ou des taux de rechute), la commande <em>glm<\/em> permet de relier la moyenne factorielle \u00e0 une distribution statistique appropri\u00e9e (famille Gamma, famille Inverse-Gaussienne) via une fonction de lien idoine (lien logarithmique ou lien identit\u00e9), contournant \u00e9l\u00e9gamment l&rsquo;hypoth\u00e8se de normalit\u00e9 :<\/p>\n<p><em>glm systolic i.drug##i.disease, family(gaussian) link(identity) vce(robust)<\/em><\/p>\n<p>Enfin, lorsque la distribution est totalement ind\u00e9termin\u00e9e et que les \u00e9chantillons sont r\u00e9duits, les techniques computationnelles de r\u00e9\u00e9chantillonnage non param\u00e9trique par bootstrap repr\u00e9sentent une alternative puissante. Stata permet d&rsquo;ex\u00e9cuter un bootstrap automatis\u00e9 sur l&rsquo;ANOVA factorielle en r\u00e9\u00e9chantillonnant des milliers de fois l&rsquo;\u00e9chantillon empirique original :<\/p>\n<p><em>bootstrap, reps(2000) seed(12345): anova systolic drug##disease<\/em><\/p>\n<p>Cette mod\u00e9lisation empirique g\u00e9n\u00e8re des intervalles de confiance bootstrapp\u00e9s par centiles affranchis de tout postulat de normalit\u00e9 th\u00e9orique pr\u00e9alable, assurant une solidit\u00e9 sans faille aux conclusions inf\u00e9rentielles de l&rsquo;\u00e9tude.<\/p>\n<h2>12. Restitution acad\u00e9mique des r\u00e9sultats et conformit\u00e9 aux normes APA<\/h2>\n<h3>12.1 Structuration textuelle des r\u00e9sultats statistiques<\/h3>\n<p>La transcription acad\u00e9mique des r\u00e9sultats statistiques doit satisfaire aux exigences formelles \u00e9dict\u00e9es par l&rsquo;American Psychological Association (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\">normes APA 7e \u00e9dition<\/a>), qui constituent le standard de communication scientifique pr\u00e9dominant en sciences du comportement, en neurosciences et en m\u00e9decine translationnelle. La r\u00e9daction doit se caract\u00e9riser par une concision irr\u00e9prochable, une transparence math\u00e9matique totale et une s\u00e9quence logique universelle.<\/p>\n<p>Le format typographique canonique pour restituer un test omnibus issu de l&rsquo;analyse de variance factorielle s&rsquo;\u00e9nonce invariablement selon le formalisme suivant :<\/p>\n<p><em>F(ddl<sub>Effet<\/sub>, ddl<sub>Erreur<\/sub>) = Valeur_F, p = Valeur_p, \u03b7<sub>p<\/sub>\u00b2 = Valeur_Taille_Effet<\/em><\/p>\n<p>Il est imp\u00e9ratif d&rsquo;observer scrupuleusement les conventions de mise en forme impos\u00e9es par l&rsquo;APA :<\/p>\n<ul>\n<li>Les lettres latines symbolisant des grandeurs statistiques (telles que <em>F<\/em>, <em>p<\/em>, <em>t<\/em>, <em>N<\/em>, <em>R\u00b2<\/em>) s&rsquo;\u00e9crivent obligatoirement en caract\u00e8res italiques ;<\/li>\n<li>Les statistiques <em>F<\/em> et les indices de taille d&rsquo;effet comme l&rsquo;Eta-carr\u00e9 partiel sont syst\u00e9matiquement tronqu\u00e9s \u00e0 deux d\u00e9cimales apr\u00e8s la virgule (par exemple <em>F(3, 46) = 9,05<\/em>) ;<\/li>\n<li>La p-valeur s&rsquo;\u00e9nonce avec trois d\u00e9cimales. Si la sortie brute de Stata indique une p-valeur de 0,0000, il est m\u00e9thodologiquement incorrect d&rsquo;\u00e9crire <em>p = .000<\/em> ; il convient imp\u00e9rativement de formaliser <em>p &lt; .001<\/em> ;<\/li>\n<li>Conform\u00e9ment aux normes nord-am\u00e9ricaines anglo-saxonnes (et adapt\u00e9 dans le style typographique francophone), le z\u00e9ro initial pr\u00e9c\u00e9dant la virgule ou le point d\u00e9cimal est omis pour les m\u00e9triques qui ne peuvent math\u00e9matiquement jamais exc\u00e9der la valeur 1 (telles que les probabilit\u00e9s <em>p<\/em> et les coefficients <em>\u03b7<sub>p<\/sub>\u00b2<\/em>, not\u00e9s conventionnellement <em>p = .024<\/em> ou <em>\u03b7<sub>p<\/sub>\u00b2 = .37<\/em>).<\/li>\n<\/ul>\n<p>Sur le plan de l&rsquo;ordonnancement discursif dans le corps du texte acad\u00e9mique, la restitution doit imp\u00e9rativement adopter une d\u00e9marche s\u00e9quentielle hi\u00e9rarchis\u00e9e :<\/p>\n<ol>\n<li>Rapporter en premier lieu le r\u00e9sultat du test d&rsquo;interaction factorielle <em>drug \u00d7 disease<\/em>, car de sa significativit\u00e9 d\u00e9pend l&rsquo;intelligibilit\u00e9 des \u00e9tapes subs\u00e9quentes ;<\/li>\n<li>Si l&rsquo;interaction s&rsquo;av\u00e8re non significative, formaliser le maintien de l&rsquo;hypoth\u00e8se additif, puis expliciter successivement les deux effets principaux globaux accompagn\u00e9s de leurs statistiques respectives ;<\/li>\n<li>Rapporter les r\u00e9sultats des comparaisons post-hoc en mentionnant explicitement la m\u00e9thode de correction retenue (ex. : \u00ab\u00a0selon la proc\u00e9dure de Tukey-Kramer \u00e0 un seuil <em>\u03b1 = .05<\/em>\u00ab\u00a0), en sp\u00e9cifiant les \u00e9carts moyens observ\u00e9s et leurs intervalles de confiance \u00e0 95 % ;<\/li>\n<li>Si l&rsquo;interaction s&rsquo;av\u00e8re significative, suspendre l&rsquo;interpr\u00e9tation des effets principaux globaux pour exposer imm\u00e9diatement la d\u00e9composition ordonn\u00e9e des effets simples au sein de chaque sous-groupe d&rsquo;int\u00e9r\u00eat.<\/li>\n<\/ol>\n<h3>12.2 Conception de tableaux et graphiques publiables<\/h3>\n<p>L&rsquo;int\u00e9gration d&rsquo;un tableau d&rsquo;analyse de variance factorielle conforme aux directives \u00e9ditoriales internationales \u00e9vite de surcharger inutilement le texte narratif de l&rsquo;article tout en offrant une transparence exhaustive sur l&rsquo;ensemble de la matrice factorielle. Selon les standards de l&rsquo;APA 7e \u00e9dition, un tableau scientifique ne doit comporter aucune ligne verticale d&rsquo;encadrement. Seules trois lignes horizontales principales doivent structurer l&rsquo;espace visuel : une ligne sup\u00e9rieure marquant le sommet du tableau, une ligne m\u00e9diane s\u00e9parant les en-t\u00eates de colonnes des donn\u00e9es chiffr\u00e9es, et une ligne inf\u00e9rieure cl\u00f4turant la base du tableau avant les notes explicatives.<\/p>\n<p>Le tableau doit comporter des colonnes d\u00e9di\u00e9es aux sources de variation (Source), \u00e0 la somme des carr\u00e9s (<em>SS<\/em>), aux degr\u00e9s de libert\u00e9 (<em>df<\/em>), au carr\u00e9 moyen (<em>MS<\/em>), \u00e0 la statistique de Fisher (<em>F<\/em>), \u00e0 la p-valeur de d\u00e9cision (<em>p<\/em>) et \u00e0 la taille d&rsquo;effet standardis\u00e9e (<em>\u03b7<sub>p<\/sub>\u00b2<\/em>). En pied de tableau, des notes explicatives documentent la signification des abr\u00e9viations, l&rsquo;usage des sommes des carr\u00e9s partielles de Type III, ainsi que le codage des seuils de significativit\u00e9 par ast\u00e9risques usuels (* <em>p &lt; .05<\/em>, ** <em>p &lt; .01<\/em>, *** <em>p &lt; .001<\/em>).<\/p>\n<p>S&rsquo;agissant des figures graphiques (profils d&rsquo;interaction ou graphiques de barres pr\u00e9dictives), celles-ci doivent imp\u00e9rativement comporter l&rsquo;affichage syst\u00e9matique des barres d&rsquo;erreur, en pr\u00e9cisant sans ambigu\u00eft\u00e9 dans la l\u00e9gende s&rsquo;il s&rsquo;agit d&rsquo;\u00e9carts-types (SD), d&rsquo;erreurs-types de la moyenne (SEM) ou d&rsquo;intervalles de confiance \u00e0 95 % (95% CI). Depuis Stata, l&rsquo;analyste veillera \u00e0 exporter ses figures dans des formats vectoriels encapsul\u00e9s (EPS ou PDF haute r\u00e9solution) ou au format raster non compress\u00e9 (TIFF avec une r\u00e9solution minimale de 300 \u00e0 600 DPI) au moyen de commandes d&rsquo;exportation script\u00e9es :<\/p>\n<p><em>graph export figure_interaction_apa.tif, width(2400) height(1800) replace<\/em><\/p>\n<p>Enfin, l&rsquo;exigence supr\u00eame de la recherche contemporaine r\u00e9side dans l&rsquo;assurance de la r\u00e9plicabilit\u00e9 universelle des analyses. Le Do-file final doit \u00eatre document\u00e9, nettoy\u00e9 de toute commande exploratoire superflue, et archiv\u00e9 aux c\u00f4t\u00e9s de la base de donn\u00e9es brute sur des entrep\u00f4ts de science ouverte (tels que l&rsquo;Open Science Framework, OSF, ou Zenodo). Cette pratique assure la tra\u00e7abilit\u00e9 compl\u00e8te de l&rsquo;investigation, renfor\u00e7ant la cr\u00e9dibilit\u00e9 acad\u00e9mique des r\u00e9sultats publi\u00e9s et participant pleinement au progr\u00e8s des standards m\u00e9thodologiques en sciences comportementales et biom\u00e9dicales.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7th ed.). American Psychological Association. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/li>\n<li>Cohen, J. (1988). <em>Statistical power analysis for the behavioral sciences<\/em> (2nd ed.). Lawrence Erlbaum Associates.<\/li>\n<li>Howell, D. C. (2012). <em>Statistical methods for psychology<\/em> (8th ed.). Cengage Learning.<\/li>\n<li>Kirk, R. E. (2013). <em>Experimental design: Procedures for the behavioral sciences<\/em> (4th ed.). SAGE Publications. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.4135\/9781483384733\">https:\/\/doi.org\/10.4135\/9781483384733<\/a><\/li>\n<li>Maxwell, S. E., Delaney, H. D., &amp; Kelley, K. (2017). <em>Designing experiments and analyzing data: A model comparison perspective<\/em> (3rd ed.). Routledge. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.4324\/9781315642956\">https:\/\/doi.org\/10.4324\/9781315642956<\/a><\/li>\n<li>Mitchell, M. N. (2021). <em>Interpreting and visualizing regression models using Stata<\/em> (2nd ed.). Stata Press.<\/li>\n<li>StataCorp. (2023). <em>Stata base reference manual: Release 18<\/em>. Stata Press. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.stata.com\/manuals\/r.pdf\">https:\/\/www.stata.com\/manuals\/r.pdf<\/a><\/li>\n<li>Tabachnick, B. G., &amp; Fidell, L. S. (2019). <em>Using multivariate statistics<\/em> (7th ed.). Pearson.<\/li>\n<li>Winer, B. J., Brown, D. R., &amp; Michels, K. M. (1991). <em>Statistical principles in experimental design<\/em> (3rd ed.). McGraw-Hill.<\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide m\u00e9thodologique complet pour r\u00e9aliser et interpr\u00e9ter une ANOVA factorielle \u00e0 deux facteurs sous Stata, de la v\u00e9rification des postulats aux normes APA.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2559,1829],"tags":[4814,4815,2304,2216,2655],"class_list":["post-2764","psych_stat","type-psych_stat","status-publish","hentry","category-analyses-de-donnees-quantitatives","category-methodologie-statistique","tag-anova-a-deux-facteurs","tag-effet-d-interaction","tag-normes-apa","tag-stata","tag-tests-post-hoc"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/2764","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=2764"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/2764\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=2764"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=2764"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=2764"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}