{"id":8002,"date":"2026-09-17T01:18:04","date_gmt":"2026-09-17T01:18:04","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-realiser-test-post-hoc-nemenyi-python\/"},"modified":"2026-09-17T01:18:04","modified_gmt":"2026-09-17T01:18:04","slug":"comment-realiser-test-post-hoc-nemenyi-python","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-realiser-test-post-hoc-nemenyi-python","title":{"rendered":"Comment r\u00e9aliser le test post-hoc de Nemenyi en Python"},"content":{"rendered":"<p>Dans le champ contemporain des sciences comportementales, de la neuropsychologie et de la psychom\u00e9trie appliqu\u00e9e, la mod\u00e9lisation rigoureuse des protocoles \u00e0 mesures r\u00e9p\u00e9t\u00e9es constitue l&rsquo;un des piliers m\u00e9thodologiques les plus critiques. Les chercheurs sont fr\u00e9quemment confront\u00e9s \u00e0 des devis exp\u00e9rimentaux intra-sujets au sein desquels une m\u00eame cohorte de participants subit de multiples conditions exp\u00e9rimentales, des \u00e9valuations longitudinales ou des variantes th\u00e9rapeutiques distinctes. Si l&rsquo;analyse de variance \u00e0 mesures r\u00e9p\u00e9t\u00e9es (ANOVA) a longtemps constitu\u00e9 l&rsquo;outil d&rsquo;inf\u00e9rence privil\u00e9gi\u00e9, sa validit\u00e9 repose sur des postulats distributionnels particuli\u00e8rement restrictifs : normalit\u00e9 multidimensionnelle des r\u00e9sidus, homosc\u00e9dasticit\u00e9 et, par-dessus tout, la stricte sph\u00e9ricit\u00e9 des matrices de covariance (crit\u00e8re de Mauchly). Or, dans la r\u00e9alit\u00e9 clinique et exp\u00e9rimentale, les donn\u00e9es psychologiques r\u00e9colt\u00e9es reposent quasi syst\u00e9matiquement sur des \u00e9chelles de type Likert, des inventaires cliniques ordinaux, des scores born\u00e9s ou des latences de r\u00e9action caract\u00e9ris\u00e9es par une asym\u00e9trie marqu\u00e9e, violant de fait les exigences fondamentales des mod\u00e8les param\u00e9triques.<\/p>\n<p>Face \u00e0 ces contraintes distributionnelles, le recours \u00e0 la statistique non param\u00e9trique devient un imp\u00e9ratif m\u00e9thodologique incontournable. Le test de Friedman s&rsquo;impose alors comme l&rsquo;alternative classique pour \u00e9prouver l&rsquo;hypoth\u00e8se globale d&rsquo;\u00e9quivalence entre plusieurs conditions appari\u00e9es. N\u00e9anmoins, l&rsquo;obtention d&rsquo;une statistique de test omnibus significative \u00e0 l&rsquo;issue de cette proc\u00e9dure soul\u00e8ve imm\u00e9diatement une impasse d\u00e9cisionnelle : si le test de Friedman permet de rejeter l&rsquo;hypoth\u00e8se nulle selon laquelle l&rsquo;ensemble des traitements produit des effets statistiquement identiques, il est par nature incapable d&rsquo;isoler pr\u00e9cis\u00e9ment quelles paires de conditions diff\u00e8rent significativement les unes des autres. Se pr\u00e9cipiter sur des comparaisons bivari\u00e9es non corrig\u00e9es exposerait l&rsquo;investigateur \u00e0 un p\u00e9ril statistique majeur : l&rsquo;inflation dramatique de l&rsquo;erreur de premi\u00e8re esp\u00e8ce (le risque alpha global), conduisant \u00e0 documenter des pseudo-d\u00e9couvertes scientifiques.<\/p>\n<p>C&rsquo;est pr\u00e9cis\u00e9ment \u00e0 cette jonction critique de l&rsquo;analyse inf\u00e9rentielle qu&rsquo;intervient le test post-hoc de Nemenyi. Con\u00e7u sp\u00e9cifiquement pour faire suite au test de Friedman dans les plans en blocs complets randomis\u00e9s, le test de Nemenyi propose une proc\u00e9dure rigoureuse d&rsquo;ajustement simultan\u00e9 pour comparaisons multiples par paires, fond\u00e9e sur la distribution de l&rsquo;\u00e9tendue studentis\u00e9e. Au sein de l&rsquo;\u00e9cosyst\u00e8me open-source moderne, le langage <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.python.org\">Python<\/a>, enrichi par des biblioth\u00e8ques scientifiques de pointe telles que <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/scipy.org\">SciPy<\/a>, <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pingouin-stats.org\">Pingouin<\/a> et <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/scikit-posthocs.readthedocs.io\">scikit-posthocs<\/a>, offre d\u00e9sormais un environnement de calcul performant et reproductible pour impl\u00e9menter cette analyse avanc\u00e9e. Cet article propose une exploration exhaustive, th\u00e9orique et appliqu\u00e9e, du test de Nemenyi en Python, guidant le chercheur depuis la formalisation math\u00e9matique des rangs jusqu&rsquo;\u00e0 la restitution des r\u00e9sultats selon les normes acad\u00e9miques les plus strictes.<\/p>\n<h2>1. Introduction au test post-hoc de Nemenyi dans la recherche psychologique<\/h2>\n<h3>1.1 Le r\u00f4le des tests post-hoc dans les protocoles \u00e0 mesures r\u00e9p\u00e9t\u00e9es<\/h3>\n<p>Dans la conduite de protocoles exp\u00e9rimentaux en psychologie cognitive ou clinique, la structuration intra-sujets (mesures r\u00e9p\u00e9t\u00e9es) pr\u00e9sente des avantages statistiques ind\u00e9niables, au premier rang desquels figure la neutralisation de la variance interindividuelle r\u00e9siduelle. En mesurant chaque sujet \u00e0 travers l&rsquo;ensemble des niveaux d&rsquo;une variable ind\u00e9pendante, le chercheur r\u00e9duit drastiquement le bruit de mesure associ\u00e9 aux diff\u00e9rences individuelles stables telles que l&rsquo;intelligence fluide, les traits de personnalit\u00e9 ou les facteurs g\u00e9n\u00e9tiques sous-jacents. Cependant, lorsqu&rsquo;une proc\u00e9dure non param\u00e9trique omnibus telle que l&rsquo;analyse de variance par rangs de Friedman signale un effet global significatif, la port\u00e9e heuristique de ce r\u00e9sultat demeure inachev\u00e9e. Rejeter l&rsquo;hypoth\u00e8se nulle globale indique uniquement qu&rsquo;au moins une condition produit des rangs diff\u00e9rents des autres, sans fournir la moindre indication sur la topographie de ces divergences.<\/p>\n<p>La n\u00e9cessit\u00e9 de d\u00e9passer l&rsquo;analyse globale implique donc de tester des contrastes par paires (pairwise comparisons). Toutefois, la r\u00e9alisation na\u00efve de multiples tests d&rsquo;hypoth\u00e8ses sur un m\u00eame jeu de donn\u00e9es engendre le ph\u00e9nom\u00e8ne redout\u00e9 de l&rsquo;inflation de l&rsquo;erreur de premi\u00e8re esp\u00e8ce, d\u00e9sign\u00e9 sous le terme technique de taux d&rsquo;erreur par famille d&rsquo;hypoth\u00e8ses (Family-Wise Error Rate, FWER). Si un chercheur \u00e9value quatre conditions th\u00e9rapeutiques distinctes sans contr\u00f4le de l&rsquo;al\u00e9a, le nombre de combinaisons par paires s&rsquo;\u00e9l\u00e8ve \u00e0 six. En appliquant le seuil conventionnel alpha de 0,05 \u00e0 chaque comparaison isol\u00e9e, la probabilit\u00e9 cumul\u00e9e de commettre au moins une erreur de type I s&rsquo;\u00e9l\u00e8ve math\u00e9matiquement \u00e0 1 &#8211; (1 &#8211; 0,05)^6, soit environ 26,5 %. Plus d&rsquo;une chance sur quatre existe alors d&rsquo;affirmer \u00e0 tort l&rsquo;existence d&rsquo;une sup\u00e9riorit\u00e9 th\u00e9rapeutique purement artefactuelle.<\/p>\n<p>Cette pr\u00e9occupation m\u00e9thodologique est particuli\u00e8rement aigu\u00eb dans le champ des sciences psychologiques, o\u00f9 la r\u00e9plicabilit\u00e9 des r\u00e9sultats fait l&rsquo;objet d&rsquo;une vigilance \u00e9pist\u00e9mologique accrue. Les variables psychom\u00e9triques r\u00e9colt\u00e9es ne satisfont presque jamais les postulats gaussiens : les \u00e9chelles composites d&rsquo;anxi\u00e9t\u00e9 (comme le score HAM-A), de d\u00e9pression (BDI-II), ou les \u00e9valuations subjectives sur \u00e9chelles ordinales d&rsquo;appr\u00e9ciation cognitive pr\u00e9sentent structurellement des planchers, des plafonds ou des asym\u00e9tries prononc\u00e9es. Dans ces contextes o\u00f9 la m\u00e9trique sous-jacente est ordinale et les distributions r\u00e9solument non normales, l&rsquo;application d&rsquo;un test post-hoc con\u00e7u sp\u00e9cifiquement pour le traitement des rangs appari\u00e9s s&rsquo;av\u00e8re l&rsquo;unique garantie d&rsquo;une inf\u00e9rence valide.<\/p>\n<h3>1.2 Positionnement \u00e9pist\u00e9mologique du test de Nemenyi<\/h3>\n<p>Sur le plan de l&rsquo;histoire de la statistique th\u00e9orique, la proc\u00e9dure de Nemenyi a \u00e9t\u00e9 formalis\u00e9e par Peter Nemenyi dans sa th\u00e8se de doctorat soutenue \u00e0 l&rsquo;Universit\u00e9 de Princeton en 1963, intitul\u00e9e <em>Distribution-free Multiple Comparisons<\/em>. Les travaux de Nemenyi s&rsquo;inscrivent dans l&rsquo;\u00e2ge d&rsquo;or du d\u00e9veloppement des m\u00e9thodes statistiques non param\u00e9triques, men\u00e9 conjointement par des pionniers tels que Frank Wilcoxon, Milton Friedman et John Tukey. L&rsquo;ambition fondamentale de Nemenyi \u00e9tait de combler un vide m\u00e9thodologique flagrant : alors que John Tukey avait magistralement \u00e9tabli sa proc\u00e9dure de la diff\u00e9rence honn\u00eatement significative (Tukey HSD) pour les comparaisons multiples param\u00e9triques bas\u00e9es sur l&rsquo;\u00e9tendue studentis\u00e9e, aucune solution rigoureusement sym\u00e9trique n&rsquo;existait pour les protocoles en blocs complets randomis\u00e9s non param\u00e9triques.<\/p>\n<p>Le test post-hoc de Nemenyi se positionne ainsi, d&rsquo;un point de vue structural, comme l&rsquo;\u00e9quivalent direct du test HSD de Tukey pour les donn\u00e9es trait\u00e9es selon la logique de classement par rangs de Friedman. Contrairement aux approches par permutations a posteriori ou aux adaptations empiriques bas\u00e9es sur des approximations asymptotiques locales, la proc\u00e9dure de Nemenyi formalise une statistique de distance critique commune applicable \u00e0 l&rsquo;ensemble des contrastes possibles. Ce cadre th\u00e9orique postule que l&rsquo;ensemble des conditions est \u00e9valu\u00e9 sur les m\u00eames unit\u00e9s exp\u00e9rimentales, garantissant ainsi un traitement uniforme de la d\u00e9pendance s\u00e9rielle inh\u00e9rente aux protocoles \u00e0 mesures r\u00e9p\u00e9t\u00e9es.<\/p>\n<p>Dans la recherche contemporaine en neurosciences cognitives, en psychopharmacologie et en mod\u00e9lisation comportementale, l&rsquo;approche de Nemenyi jouit d&rsquo;une l\u00e9gitimit\u00e9 \u00e9pist\u00e9mologique pr\u00e9pond\u00e9rante. Lorsque des chercheurs mesurent les performances attentionnelles d&rsquo;un groupe sous trois mol\u00e9cules distinctes ou \u00e9valuent l&rsquo;activation neuronale dans des r\u00e9gions d&rsquo;int\u00e9r\u00eat au cours de quatre t\u00e2ches ex\u00e9cutives cons\u00e9cutives, la structure des donn\u00e9es impose une robustesse sans faille face aux d\u00e9viations de distribution. L&rsquo;utilisation du test de Nemenyi atteste d&rsquo;une d\u00e9marche scientifique rigoureuse, soucieuse de ne pas contraindre artificiellement des donn\u00e9es psychom\u00e9triques dans des mod\u00e8les gaussiens inadapt\u00e9s tout en maintenant un verrouillage math\u00e9matique \u00e9tanche sur le risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce.<\/p>\n<h3>1.3 Objectifs et structure g\u00e9n\u00e9rale de l&rsquo;article<\/h3>\n<p>L&rsquo;objectif de cet article est de doter les chercheurs, analystes de donn\u00e9es et cliniciens d&rsquo;un guide m\u00e9thodologique et pratique complet pour conduire le test post-hoc de Nemenyi dans l&rsquo;\u00e9cosyst\u00e8me computationnel Python. Bien que la litt\u00e9rature statistique internationale fournisse les principes abstraits du test, son op\u00e9rationnalisation concr\u00e8te soul\u00e8ve fr\u00e9quemment des difficult\u00e9s techniques : manipulation de formats de donn\u00e9es complexes, s\u00e9lection judicieuse des biblioth\u00e8ques, gestion des ex-aequo, et restitution narrative conforme aux exigences de l&rsquo;American Psychological Association (APA).<\/p>\n<p>Afin d&rsquo;offrir une couverture exhaustive, la structure de ce guide s&rsquo;articule autour de quatre axes majeurs :<\/p>\n<ul>\n<li><strong>La formalisation th\u00e9orique et math\u00e9matique :<\/strong> Nous expliciterons la m\u00e9canique du classement par rangs intra-sujets, la d\u00e9rivation analytique de la diff\u00e9rence critique (Critical Difference) et la comparaison syst\u00e9matique avec les proc\u00e9dures concurrentes (Wilcoxon-Bonferroni, Conover-Iman).<\/li>\n<li><strong>Le protocole pr\u00e9alable et l&rsquo;ing\u00e9nierie des donn\u00e9es :<\/strong> Nous aborderons l&rsquo;ex\u00e9cution imp\u00e9rative du test de Friedman, l&rsquo;architecture des environnements virtuels Python et la pr\u00e9paration des matrices de donn\u00e9es du format large (wide) au format long (tidy).<\/li>\n<li><strong>L&rsquo;impl\u00e9mentation algorithmique pas \u00e0 pas :<\/strong> Nous d\u00e9taillerons l&rsquo;usage de la biblioth\u00e8que <code>scikit-posthocs<\/code>, la comparaison avec des calculs matriciels purs sous NumPy\/Pandas et l&rsquo;extraction des matrices de p-valeurs.<\/li>\n<li><strong>L&rsquo;interpr\u00e9tation clinique, la visualisation et la diffusion :<\/strong> Nous pr\u00e9senterons la construction de diagrammes de diff\u00e9rence critique (CD diagrams), la cartographie thermique des contrastes, ainsi que l&rsquo;application int\u00e9grale d&rsquo;un cas clinique en psychoth\u00e9rapie \u00e9valu\u00e9 selon les normes de publication acad\u00e9mique.<\/li>\n<\/ul>\n<h2>2. Fondements th\u00e9oriques et math\u00e9matiques du test de Nemenyi<\/h2>\n<h3>2.1 Principe du classement par rangs intra-sujets<\/h3>\n<p>Le point de d\u00e9part algorithmique du test de Nemenyi repose sur le protocole de transformation par rangs initi\u00e9 par le test de Friedman. Consid\u00e9rons une matrice de donn\u00e9es exp\u00e9rimentales comportant <em>N<\/em> sujets (les blocs) \u00e9valu\u00e9s sous <em>k<\/em> conditions exp\u00e9rimentales ou traitements. Contrairement \u00e0 une proc\u00e9dure non param\u00e9trique pour \u00e9chantillons ind\u00e9pendants (telle que celle de Kruskal-Wallis, qui classe l&rsquo;ensemble des observations de fa\u00e7on globale du rang 1 \u00e0 <em>N &times; k<\/em>), le classement de Friedman et Nemenyi s&rsquo;op\u00e8re strictement \u00e0 l&rsquo;int\u00e9rieur de chaque sujet individuel, c&rsquo;est-\u00e0-dire ligne par ligne.<\/p>\n<p>Pour chaque participant <em>i<\/em> (avec <em>i<\/em> variant de 1 \u00e0 <em>N<\/em>), les scores observ\u00e9s \u00e0 travers les <em>k<\/em> conditions sont ordonn\u00e9s de mani\u00e8re croissante. Le score le plus faible au sein de la ligne re\u00e7oit le rang 1, le score imm\u00e9diatement sup\u00e9rieur le rang 2, jusqu&rsquo;au score le plus \u00e9lev\u00e9 qui se voit attribuer le rang <em>k<\/em>. Ce classement intra-sujet neutralise instantan\u00e9ment l&rsquo;amplitude absolue des scores du sujet ainsi que son niveau moyen de base : un individu pr\u00e9sentant un biais syst\u00e9matique de r\u00e9ponse (par exemple, une tendance \u00e0 noter syst\u00e9matiquement haut sur une \u00e9chelle psychom\u00e9trique) est ramen\u00e9 exactement \u00e0 la m\u00eame \u00e9chelle ordinale normalis\u00e9e que les autres participants de la cohorte.<\/p>\n<p>Une fois cette op\u00e9ration de substitution matricielle effectu\u00e9e pour les <em>N<\/em> blocs, le test calcule pour chaque condition exp\u00e9rimentale <em>j<\/em> (avec <em>j<\/em> variant de 1 \u00e0 <em>k<\/em>) la somme de ses rangs, not\u00e9e <em>R<sub>j<\/sub><\/em>, puis la moyenne de ses rangs, d\u00e9sign\u00e9e par <em>R&#772;<sub>j<\/sub> = R<sub>j<\/sub> \/ N<\/em>. En vertu des propri\u00e9t\u00e9s arithm\u00e9tiques des permutations, la somme th\u00e9orique de tous les rangs attribu\u00e9s \u00e0 l&rsquo;\u00e9chelle d&rsquo;un sujet est invariante et \u00e9gale \u00e0 <em>k(k + 1) \/ 2<\/em>. La moyenne th\u00e9orique des rangs sous l&rsquo;hypoth\u00e8se nulle d&rsquo;absence de diff\u00e9rence entre les conditions est ainsi rigoureusement identique pour chaque traitement et vaut <em>(k + 1) \/ 2<\/em>. Une propri\u00e9t\u00e9 fondamentale de cette m\u00e9trique est son invariance absolue face \u00e0 toute transformation monotone croissante des scores d&rsquo;origine : que la variable psychom\u00e9trique soit exprim\u00e9e en valeurs brutes, en logarithmes ou soumise \u00e0 une dilatation non lin\u00e9aire, le r\u00e9sultat inf\u00e9rentiel demeure parfaitement inchang\u00e9.<\/p>\n<h3>2.2 Formulation math\u00e9matique de la diff\u00e9rence critique<\/h3>\n<p>Le c\u0153ur op\u00e9ratoire de la proc\u00e9dure de Nemenyi r\u00e9side dans la d\u00e9finition d&rsquo;un seuil m\u00e9trique unique au-del\u00e0 duquel la divergence entre deux rangs moyens est jug\u00e9e statistiquement significative \u00e0 un risque alpha pr\u00e9d\u00e9termin\u00e9. Cette valeur seuil est formellement appel\u00e9e la diff\u00e9rence critique (Critical Difference, not\u00e9e <em>CD<\/em>). Deux conditions exp\u00e9rimentales <em>A<\/em> et <em>B<\/em> pr\u00e9sentent une divergence significative si et seulement si la valeur absolue de la diff\u00e9rence de leurs rangs moyens respectifs est strictement sup\u00e9rieure \u00e0 cette valeur seuil :<\/p>\n<p><strong>| R&#772;<sub>A<\/sub> &#8211; R&#772;<sub>B<\/sub> | &gt; CD<\/strong><\/p>\n<p>L&rsquo;expression analytique exacte de la diff\u00e9rence critique est r\u00e9gie par l&rsquo;\u00e9quation suivante :<\/p>\n<p><strong>CD = q<sub>&alpha;, k, &infin;<\/sub> &times; &radic; [ (k &times; (k + 1)) \/ (6 &times; N) ]<\/strong><\/p>\n<p>Dans cette formulation math\u00e9matique, chaque param\u00e8tre refl\u00e8te une dimension pr\u00e9cise du devis exp\u00e9rimental :<\/p>\n<ul>\n<li><em>k<\/em> d\u00e9signe le nombre total de conditions exp\u00e9rimentales ou de traitements compar\u00e9s.<\/li>\n<li><em>N<\/em> repr\u00e9sente le nombre d&rsquo;unit\u00e9s exp\u00e9rimentales ind\u00e9pendantes (sujets ou blocs complets).<\/li>\n<li>Le terme <em>&radic; [ (k(k + 1)) \/ (6N) ]<\/em> correspond \u00e0 l&rsquo;\u00e9cart-type de la distribution d&rsquo;\u00e9chantillonnage de la diff\u00e9rence entre deux rangs moyens sous l&rsquo;hypoth\u00e8se nulle. Cette variance est directement d\u00e9riv\u00e9e des propri\u00e9t\u00e9s de la loi uniforme discr\u00e8te r\u00e9gissant les permutations de rangs.<\/li>\n<li>Le coefficient <em>q<sub>&alpha;, k, &infin;<\/sub><\/em> est la valeur critique tabul\u00e9e issue de la <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Studentized_range_distribution\">distribution de l&rsquo;\u00e9tendue studentis\u00e9e<\/a> (Studentized Range Distribution), fix\u00e9e pour un niveau de confiance <em>1 &#8211; &alpha;<\/em>, pour un nombre de groupes \u00e9gal \u00e0 <em>k<\/em>, et pour des degr\u00e9s de libert\u00e9 infinis (refl\u00e9tant le fait que la variance des rangs sous l&rsquo;hypoth\u00e8se nulle est math\u00e9matiquement constante et connue d&rsquo;avance, ne n\u00e9cessitant pas d&rsquo;estimation d&rsquo;un param\u00e8tre d&rsquo;\u00e9chelle r\u00e9siduel sur les donn\u00e9es).<\/li>\n<\/ul>\n<p>Il est fondamental de noter que l&rsquo;amplitude de la diff\u00e9rence critique se contracte \u00e0 mesure que la taille d&rsquo;\u00e9chantillon <em>N<\/em> cro\u00eet, conf\u00e9rant au test une puissance accrue pour discriminer des variations subtiles d&rsquo;effets dans les cohortes volumineuses. \u00c0 l&rsquo;inverse, une augmentation du nombre de conditions <em>k<\/em> accro\u00eet la valeur critique <em>q<\/em> ainsi que le num\u00e9rateur sous la racine, traduisant le co\u00fbt statistique inh\u00e9rent \u00e0 la multiplication des combinaisons explor\u00e9es.<\/p>\n<h3>2.3 Comparaison avec d&rsquo;autres proc\u00e9dures post-hoc non param\u00e9triques<\/h3>\n<p>Face \u00e0 la n\u00e9cessit\u00e9 de conduire des comparaisons par paires apr\u00e8s un test de Friedman, le chercheur dispose d&rsquo;alternatives m\u00e9thodologiques dont les fondements diff\u00e8rent profond\u00e9ment du test de Nemenyi. La strat\u00e9gie la plus fr\u00e9quemment rencontr\u00e9e dans la pratique consiste \u00e0 ex\u00e9cuter une s\u00e9rie de tests des rangs sign\u00e9s de Wilcoxon pour chaque paire possible de conditions, en appliquant manuellement une proc\u00e9dure d&rsquo;ajustement de la p-valeur de type Bonferroni ou Holm. Bien qu&rsquo;intuitive, cette approche pr\u00e9sente une rupture \u00e9pist\u00e9mologique : les tests de Wilcoxon r\u00e9assignent de nouveaux rangs bas\u00e9s exclusivement sur les diff\u00e9rences observ\u00e9es au sein de chaque paire isol\u00e9e, ignorant totalement la matrice globale des rangs \u00e9tablie lors de l&rsquo;analyse omnibus. De surcro\u00eet, la correction classique de Bonferroni est r\u00e9put\u00e9e pour son hyper-conservatisme, accroissant notablement le risque de commettre des erreurs de deuxi\u00e8me esp\u00e8ce (b\u00eata).<\/p>\n<p>Une seconde alternative r\u00e9side dans le test post-hoc de Conover-Iman (1979). Ce test repose sur une logique d&rsquo;approximation par la loi de <em>t<\/em> de Student appliqu\u00e9e directement aux rangs de Friedman, tout en int\u00e9grant la variance d&rsquo;erreur r\u00e9siduelle du mod\u00e8le global. La m\u00e9thode de Conover-Iman s&rsquo;av\u00e8re consid\u00e9rablement plus puissante (plus lib\u00e9rale) que celle de Nemenyi. Cependant, cette puissance statistique sup\u00e9rieure s&rsquo;accompagne d&rsquo;une vuln\u00e9rabilit\u00e9 accrue aux faux positifs d\u00e8s lors que les hypoth\u00e8ses d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances de rangs entre conditions sont mises \u00e0 mal.<\/p>\n<p>Le test de Nemenyi se distingue par son conservatisme intrins\u00e8que et sa coh\u00e9rence structurelle avec l&rsquo;esprit du test de Tukey HSD : il \u00e9value l&rsquo;ensemble des contrastes sous une unique m\u00e9trique d&rsquo;\u00e9tendue conjointe, r\u00e9gulant de mani\u00e8re inflexible le FWER sur l&rsquo;ensemble de la famille des comparaisons. En recherche clinique et en neuropsychologie confirmatoire, ce conservatisme repr\u00e9sente une garantie essentielle pour pr\u00e9venir la prolif\u00e9ration de revendications th\u00e9rapeutiques infond\u00e9es.<\/p>\n<h2>3. Le pr\u00e9requis fondamental : Le test de Friedman pour mesures r\u00e9p\u00e9t\u00e9es<\/h2>\n<h3>3.1 Hypoth\u00e8ses et conditions d&rsquo;application de Friedman<\/h3>\n<p>L&rsquo;ex\u00e9cution l\u00e9gitime du test post-hoc de Nemenyi est strictement conditionn\u00e9e par la validation pr\u00e9alable du test omnibus de Friedman. Con\u00e7u par l&rsquo;\u00e9conomiste Milton Friedman en 1937, ce test constitue l&rsquo;extension non param\u00e9trique de l&rsquo;ANOVA \u00e0 un facteur intra-sujets pour plans en blocs al\u00e9atoires complets. Avant de mobiliser cette architecture d&rsquo;analyse, le chercheur doit scrupuleusement auditer le respect de trois conditions m\u00e9thodologiques indispensables :<\/p>\n<ul>\n<li><strong>L&rsquo;appariement des observations :<\/strong> Les donn\u00e9es doivent \u00eatre structur\u00e9es en blocs mutuellement exclusifs o\u00f9 chaque bloc repr\u00e9sente un sujet individuel soumis \u00e0 l&rsquo;ensemble des <em>k<\/em> conditions exp\u00e9rimentales. Le protocole interdit formellement les attributions ind\u00e9pendantes de sujets \u00e0 des groupes s\u00e9par\u00e9s dans cette configuration.<\/li>\n<li><strong>La nature ordinale de la variable d\u00e9pendante :<\/strong> L&rsquo;\u00e9chelle de mesure doit se situer au minimum \u00e0 un niveau ordinal, autorisant une op\u00e9ration de tri strict (inf\u00e9rieur, sup\u00e9rieur ou \u00e9gal). Les donn\u00e9es continues (intervalles ou ratios) sont \u00e9ligibles car elles se pr\u00eatent naturellement au reclassement par rangs, tandis que les donn\u00e9es purement cat\u00e9gorielles non ordonn\u00e9es (nominales) sont formellement exclues.<\/li>\n<li><strong>L&rsquo;ind\u00e9pendance mutuelle des blocs :<\/strong> Si les observations r\u00e9alis\u00e9es \u00e0 l&rsquo;int\u00e9rieur d&rsquo;un m\u00eame sujet sont d\u00e9pendantes par construction, les blocs entre eux doivent \u00eatre rigoureusement ind\u00e9pendants. La r\u00e9ponse fournie par le sujet <em>A<\/em> ne doit avoir aucune influence sur celle du sujet <em>B<\/em>, excluant les situations de contagion sociale, de passation collective non contr\u00f4l\u00e9e ou de dyades interactives.<\/li>\n<\/ul>\n<h3>3.2 Formulation des hypoth\u00e8ses nulle et alternative<\/h3>\n<p>Dans le cadre formalis\u00e9 du test de Friedman, les hypoth\u00e8ses d&rsquo;inf\u00e9rence statistique se d\u00e9finissent non pas sur les moyennes arithm\u00e9tiques ou les \u00e9carts-types des distributions (inop\u00e9rants en analyse non param\u00e9trique), mais sur les rangs moyens th\u00e9oriques ou les m\u00e9dianes sous-jacentes. L&rsquo;hypoth\u00e8se nulle (<em>H<sub>0<\/sub><\/em>) stipule qu&rsquo;il n&rsquo;existe aucune divergence syst\u00e9matique entre les conditions exp\u00e9rimentales au sein de la population parente :<\/p>\n<p><strong>H<sub>0<\/sub> : &theta;<sub>1<\/sub> = &theta;<sub>2<\/sub> = &#8230; = &theta;<sub>k<\/sub><\/strong> (les distributions associ\u00e9es aux <em>k<\/em> traitements sont interchangeables et poss\u00e8dent des rangs moyens th\u00e9oriques identiques \u00e9quivalents \u00e0 <em>(k + 1) \/ 2<\/em>).<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;hypoth\u00e8se alternative (<em>H<sub>1<\/sub><\/em>) ne postule pas une direction ordonn\u00e9e globale, mais affirme l&rsquo;existence d&rsquo;une rupture d&rsquo;homog\u00e9n\u00e9it\u00e9 :<\/p>\n<p><strong>H<sub>1<\/sub> : Il existe au moins une condition j dont la distribution des rangs diff\u00e8re syst\u00e9matiquement des autres conditions.<\/strong><\/p>\n<p>La statistique de test de Friedman, d\u00e9sign\u00e9e traditionnellement par <em>Q<\/em> ou &chi;<sup>2<\/sup><sub>F<\/sub>, suit asymptotiquement une loi du Chi-deux \u00e0 <em>k &#8211; 1<\/em> degr\u00e9s de libert\u00e9 lorsque la taille d&rsquo;\u00e9chantillon <em>N<\/em> est suffisamment grande. La r\u00e8gle de d\u00e9cision inf\u00e9rentielle est unilat\u00e9rale : si la probabilit\u00e9 associ\u00e9e \u00e0 la statistique observ\u00e9e sous <em>H<sub>0<\/sub><\/em> (la p-valeur) descend sous le seuil d&rsquo;erreur alpha conventionnel (typiquement 0,05), l&rsquo;hypoth\u00e8se nulle est rejet\u00e9e. Cette d\u00e9cision statistique formelle constitue le feu vert m\u00e9thodologique autorisant l&rsquo;investigateur \u00e0 d\u00e9ployer l&rsquo;analyse post-hoc de Nemenyi.<\/p>\n<h3>3.3 Risques m\u00e9thodologiques li\u00e9s \u00e0 l&rsquo;omission du test omnibus<\/h3>\n<p>Dans la pratique contemporaine de l&rsquo;analyse de donn\u00e9es, une d\u00e9rive fr\u00e9quente consiste \u00e0 s&rsquo;affranchir du test omnibus de Friedman pour engager imm\u00e9diatement des comparaisons par paires cibl\u00e9es, sous pr\u00e9texte d&rsquo;hypoth\u00e8ses th\u00e9oriques pr\u00e9alables (a priori). Bien que l&rsquo;approche des contrastes planifi\u00e9s soit th\u00e9oriquement d\u00e9fendable dans un cadre param\u00e9trique lin\u00e9aire sous certaines conditions d&rsquo;orthogonalit\u00e9 stricte, son extrapolation brute aux classements par rangs non param\u00e9triques comporte des p\u00e9rils statistiques majeurs.<\/p>\n<p>L&rsquo;omission du test omnibus expose l&rsquo;analyste au risque syst\u00e9mique de d\u00e9couverte fortuite cons\u00e9cutive au dragage de donn\u00e9es (data dredging ou p-hacking involontaire). En explorant l&rsquo;espace des diff\u00e9rences relatives sans v\u00e9rifier pr\u00e9alablement que la dispersion globale des rangs exc\u00e8de ce que le pur hasard de permutation peut g\u00e9n\u00e9rer, le chercheur fragilise consid\u00e9rablement la validit\u00e9 interne de son protocole. Dans les revues scientifiques \u00e0 comit\u00e9 de lecture en psychologie exp\u00e9rimentale et en neurosciences, l&rsquo;\u00e9valuation par les pairs sanctionne syst\u00e9matiquement l&rsquo;absence de test omnibus : la significativit\u00e9 pr\u00e9alable de Friedman est l&rsquo;ancrage indispensable qui justifie le partitionnement subs\u00e9quent de la variance des rangs.<\/p>\n<h2>4. Configuration de l&rsquo;environnement Python et biblioth\u00e8ques requises<\/h2>\n<h3>4.1 S\u00e9lection des packages d&rsquo;analyse statistique et de manipulation<\/h3>\n<p>Le d\u00e9ploiement d&rsquo;un pipeline d&rsquo;analyse non param\u00e9trique robuste en Python requiert une combinaison cibl\u00e9e de biblioth\u00e8ques sp\u00e9cialis\u00e9es assurant le traitement matriciel, le calcul inf\u00e9rentiel et la restitution graphique. Pour mener \u00e0 bien le test de Friedman et sa suite de Nemenyi, quatre outils logiciels constituent l&rsquo;ossature computationnelle de r\u00e9f\u00e9rence :<\/p>\n<ul>\n<li><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/numpy.org\">NumPy<\/a> : La biblioth\u00e8que fondamentale pour le calcul scientifique vectoris\u00e9, indispensable \u00e0 la manipulation de matrices multidimensionnelles et au tri rapide des rangs.<\/li>\n<li><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pandas.pydata.org\">Pandas<\/a> : Le moteur d&rsquo;ing\u00e9nierie des donn\u00e9es par excellence en Python, permettant d&rsquo;ing\u00e9rer, d&rsquo;indexer, de filtrer et de transformer les structures de donn\u00e9es psychom\u00e9triques entre formats tabulaires h\u00e9t\u00e9rog\u00e8nes.<\/li>\n<li><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/scipy.org\">SciPy (sous-module scipy.stats)<\/a> : La r\u00e9f\u00e9rence internationale pour l&rsquo;inf\u00e9rence statistique en Python, fournissant l&rsquo;impl\u00e9mentation standardis\u00e9e de la statistique du Chi-deux de Friedman via la fonction <code>friedmanchisquare<\/code>.<\/li>\n<li><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/scikit-posthocs.readthedocs.io\">scikit-posthocs<\/a> : Con\u00e7ue par Maksim Terpilowski, cette biblioth\u00e8que de pointe comble un vide historique de l&rsquo;\u00e9cosyst\u00e8me Python en offrant un \u00e9ventail complet de proc\u00e9dures post-hoc non param\u00e9triques, dont le test de Nemenyi pour blocs complets.<\/li>\n<li><a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pingouin-stats.org\">Pingouin<\/a> : Un package statistique moderne orient\u00e9 vers les sciences comportementales et biom\u00e9dicales, particuli\u00e8rement pr\u00e9cieux pour calculer instantan\u00e9ment les tailles d&rsquo;effet associ\u00e9es (W de Kendall).<\/li>\n<\/ul>\n<h3>4.2 Installation et gestion des d\u00e9pendances<\/h3>\n<p>Afin de garantir la parfaite reproductibilit\u00e9 des calculs et de pr\u00e9venir tout conflit de d\u00e9pendances entre biblioth\u00e8ques math\u00e9matiques, il est vivement recommand\u00e9 d&rsquo;isoler l&rsquo;environnement analytique \u00e0 l&rsquo;aide d&rsquo;un gestionnaire d&rsquo;environnements virtuels tel que <code>venv<\/code> ou <code>conda<\/code>. L&rsquo;installation des packages s&rsquo;ex\u00e9cute de mani\u00e8re s\u00e9quentielle via le gestionnaire de paquets standard <code>pip<\/code> au sein de votre terminal syst\u00e8me :<\/p>\n<p><code>pip install numpy pandas scipy scikit-posthocs pingouin matplotlib seaborn<\/code><\/p>\n<p>Lors de l&rsquo;initialisation de votre script ou de votre carnet interactif (Jupyter Notebook), il convient de proc\u00e9der \u00e0 l&rsquo;importation propre et explicite des modules en consignant scrupuleusement leurs alias canoniques :<\/p>\n<p><code>import numpy as np<\/code><br \/>\n<code>import pandas as pd<\/code><br \/>\n<code>import scipy.stats as stats<\/code><br \/>\n<code>import scikit_posthocs as sp<\/code><br \/>\n<code>import pingouin as pg<\/code><br \/>\n<code>import matplotlib.pyplot as plt<\/code><br \/>\n<code>import seaborn as sns<\/code><\/p>\n<p>Une bonne pratique en psychom\u00e9trie computationnelle consiste \u00e0 v\u00e9rifier et tracer syst\u00e9matiquement les num\u00e9ros de version des packages engag\u00e9s (notamment <code>scipy.__version__<\/code> et <code>scikit_posthocs.__version__<\/code>) afin d&rsquo;assurer qu&rsquo;une r\u00e9ex\u00e9cution ult\u00e9rieure du code ne subira aucune alt\u00e9ration cons\u00e9cutive \u00e0 une mise \u00e0 jour d&rsquo;API.<\/p>\n<h3>4.3 Outils de visualisation compl\u00e9mentaires<\/h3>\n<p>L&rsquo;interpr\u00e9tation statistique ne saurait \u00eatre compl\u00e8te sans une restitution graphique de haute pr\u00e9cision facilitant la lecture visuelle des contrastes d&rsquo;hypoth\u00e8ses. Le package <code>matplotlib<\/code> sert de moteur de rendu de bas niveau pour moduler les param\u00e8tres vectoriels (dpi, typographie, axes num\u00e9riques, bo\u00eetes d&rsquo;encadrement), condition indispensable pour pr\u00e9parer des figures destin\u00e9es aux revues index\u00e9es.<\/p>\n<p>En compl\u00e9ment, <code>seaborn<\/code> apporte une couche d&rsquo;abstraction stylistique et statistique hautement sophistiqu\u00e9e, optimis\u00e9e pour le tra\u00e7age des matrices de covariance et des diagrammes de dispersion. Enfin, <code>scikit-posthocs<\/code> int\u00e8gre des routines graphiques sp\u00e9cialis\u00e9es, notamment pour g\u00e9n\u00e9rer nativement les diagrammes de diff\u00e9rence critique conceptualis\u00e9s par Janez Dem\u0161ar, lesquels repr\u00e9sentent le standard absolu de restitution visuelle pour le test de Nemenyi.<\/p>\n<h2>5. Structuration et pr\u00e9paration des donn\u00e9es psychom\u00e9triques<\/h2>\n<h3>5.1 Mod\u00e9lisation d&rsquo;une \u00e9tude en psychopharmacologie cognitive<\/h3>\n<p>Afin d&rsquo;ancrer notre d\u00e9marche dans un r\u00e9alisme exp\u00e9rimental tangible, nous mod\u00e9liserons une \u00e9tude fictive mais scientifiquement rigoureuse men\u00e9e en psychopharmacologie cognitive. L&rsquo;investigation s&rsquo;int\u00e9resse \u00e0 l&rsquo;effet de modulateurs de l&rsquo;\u00e9veil sur la vitesse de traitement de l&rsquo;information aupr\u00e8s d&rsquo;une cohorte de patients adultes pr\u00e9sentant un trouble du d\u00e9ficit de l&rsquo;attention avec hyperactivit\u00e9 (TDAH). La cohorte est constitu\u00e9e de <em>N = 12<\/em> participants inclus dans un protocole en cross-over (mesures r\u00e9p\u00e9t\u00e9es compl\u00e8tes avec contrebalancement de l&rsquo;ordre de passation et p\u00e9riodes de sevrage th\u00e9rapeutique interm\u00e9diaires de type wash-out).<\/p>\n<p>Chaque patient est soumis successivement \u00e0 trois conditions pharmacologiques distinctes :<\/p>\n<ul>\n<li><strong>Condition Placebo :<\/strong> Administration d&rsquo;une g\u00e9lule inerte d&rsquo;aspect identique.<\/li>\n<li><strong>Condition Mol\u00e9cule A :<\/strong> Administration d&rsquo;un psychostimulant classique \u00e0 lib\u00e9ration prolong\u00e9e.<\/li>\n<li><strong>Condition Mol\u00e9cule B :<\/strong> Administration d&rsquo;un nouvel agoniste non stimulant ciblant s\u00e9lectivement les r\u00e9cepteurs noradr\u00e9nergiques.<\/li>\n<\/ul>\n<p>La variable d\u00e9pendante mesur\u00e9e correspond au temps de r\u00e9action m\u00e9dian (exprim\u00e9 en millisecondes) lors d&rsquo;une t\u00e2che d&rsquo;attention soutenue sur ordinateur (Continuous Performance Task). Les temps de r\u00e9action dans ce type de paradigme sont bien connus pour leur asym\u00e9trie positive persistante et leur vuln\u00e9rabilit\u00e9 aux micro-d\u00e9crochages attentionnels, rendant l&#8217;emploi de statistiques non param\u00e9triques particuli\u00e8rement justifi\u00e9.<\/p>\n<h3>5.2 Formatage des donn\u00e9es : format large (wide) versus long (tidy)<\/h3>\n<p>La manipulation des donn\u00e9es en Python confronte syst\u00e9matiquement le chercheur \u00e0 deux structures d&rsquo;organisation tabulaire fondamentales : le format large (wide format) et le format long (long format ou tidy data). La ma\u00eetrise de leur conversion r\u00e9ciproque est un imp\u00e9ratif technique car diff\u00e9rentes fonctions statistiques exigent l&rsquo;un ou l&rsquo;autre de ces paradigmes.<\/p>\n<p>Dans le format large, chaque ligne repr\u00e9sente un sujet unique (un bloc ind\u00e9pendant), et chaque condition exp\u00e9rimentale fait l&rsquo;objet d&rsquo;une colonne d\u00e9di\u00e9e. Cette structure correspond \u00e0 la vue matricielle math\u00e9matique naturelle du devis en blocs complets :<\/p>\n<p><code>data_wide = pd.DataFrame({<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'Sujet': [f'S{i:02d}' for i in range(1, 13)],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'Placebo': [450, 480, 520, 410, 490, 560, 430, 510, 470, 530, 440, 500],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'Molecule_A': [380, 410, 430, 370, 420, 460, 390, 440, 400, 450, 380, 420],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'Molecule_B': [410, 440, 470, 390, 450, 500, 410, 470, 430, 490, 400, 460]<br \/>\n})<\/code><\/p>\n<p>\u00c0 l&rsquo;inverse, le format long (standardis\u00e9 par Hadley Wickham) organise la table de telle sorte que chaque ligne contienne une seule observation unitaire. La table comporte alors trois colonnes : l&rsquo;identifiant du sujet, le libell\u00e9 de la condition exp\u00e9rimentale, et la valeur mesur\u00e9e. La conversion du format large vers le format long s&rsquo;op\u00e8re \u00e9l\u00e9gamment en Pandas via la m\u00e9thode <code>melt()<\/code> :<\/p>\n<p><code>data_long = data_wide.melt(id_vars=['Sujet'], var_name='Condition', value_name='Temps_Reaction')<\/code><\/p>\n<p>Tandis que la fonction <code>scipy.stats.friedmanchisquare<\/code> attend des vecteurs de colonnes isol\u00e9s ou des matrices larges, des packages tels que <code>pingouin<\/code> ou <code>seaborn<\/code> requi\u00e8rent imp\u00e9rativement le format long pour param\u00e9trer leurs arguments fonctionnels.<\/p>\n<h3>5.3 Traitement des valeurs manquantes et des ex-aequo (ties)<\/h3>\n<p>Deux perturbations m\u00e9trologiques majeures peuvent survenir lors de la constitution de la matrice exp\u00e9rimentale : la survenue de donn\u00e9es manquantes (missing data) et l&rsquo;apparition de scores identiques au sein d&rsquo;un m\u00eame sujet (ex-aequo ou ties). Le test classique de Friedman ainsi que la proc\u00e9dure post-hoc de Nemenyi reposent sur le postulat d&rsquo;un plan exp\u00e9rimental en blocs complets \u00e9quilibr\u00e9s. La pr\u00e9sence d&rsquo;une seule valeur manquante au sein d&rsquo;une ligne rend caduque l&rsquo;attribution ordinale normalis\u00e9e pour ce sujet.<\/p>\n<p>Si un sujet pr\u00e9sente une omission de mesure sur l&rsquo;une des conditions, la proc\u00e9dure standard impose soit l&rsquo;exclusion de l&rsquo;ensemble du bloc (listwise deletion via <code>data_wide.dropna()<\/code>), soit le recours \u00e0 des algorithmes d&rsquo;imputation multiple non param\u00e9trique sp\u00e9cifiquement adapt\u00e9s aux rangs. Dans le cadre de protocoles cliniques o\u00f9 chaque patient est pr\u00e9cieux, l&rsquo;exclusion aveugle doit \u00eatre \u00e9vit\u00e9e si le taux d&rsquo;attrition est important, conduisant le chercheur vers des mod\u00e8les mixtes pour donn\u00e9es ordinales (voir section 12).<\/p>\n<p>Concernant les ex-aequo, lorsque deux scores observ\u00e9s au sein d&rsquo;un bloc pr\u00e9sentent exactement la m\u00eame valeur m\u00e9trique, l&rsquo;algorithme leur attribue le rang moyen des positions qu&rsquo;ils auraient occup\u00e9es conjointement. Par exemple, si deux scores se partagent la deuxi\u00e8me et la troisi\u00e8me place, ils re\u00e7oivent chacun le rang 2,5. Si la pr\u00e9sence d&rsquo;ex-aequo modifie l\u00e9g\u00e8rement la variance th\u00e9orique sous-jacente des rangs, le test de Nemenyi demeure robuste d\u00e8s lors que la proportion de rangs partag\u00e9s reste mod\u00e9r\u00e9e face au volume global des donn\u00e9es.<\/p>\n<h2>6. Ex\u00e9cution du test de Friedman pr\u00e9alable en Python<\/h2>\n<h3>6.1 Impl\u00e9mentation avec scipy.stats.friedmanchisquare<\/h3>\n<p>La biblioth\u00e8que <code>scipy.stats<\/code> offre l&rsquo;acc\u00e8s le plus direct et universel pour v\u00e9rifier l&rsquo;hypoth\u00e8se omnibus de Friedman. L&rsquo;interface de la fonction <code>friedmanchisquare<\/code> n\u00e9cessite de lui passer en arguments positionnels les diff\u00e9rents vecteurs ou s\u00e9ries de donn\u00e9es correspondant \u00e0 chaque traitement exp\u00e9rimental. En reprenant notre jeu de donn\u00e9es format\u00e9 sous forme de DataFrame large, l&rsquo;ex\u00e9cution s&rsquo;articule comme suit :<\/p>\n<p><code>stat, p_value = stats.friedmanchisquare(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;data_wide['Placebo'],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;data_wide['Molecule_A'],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;data_wide['Molecule_B']<br \/>\n)<\/code><\/p>\n<p>Pour notre cohorte de 12 patients, la fonction extrait la statistique de Friedman (d\u00e9not\u00e9e conventionnellement <em>Q<\/em> ou &chi;<sup>2<\/sup>) et calcule la probabilit\u00e9 asymptotique associ\u00e9e bas\u00e9e sur <em>k &#8211; 1 = 2<\/em> degr\u00e9s de libert\u00e9. Supposons que l&rsquo;ex\u00e9cution renvoie <code>stat = 24.0<\/code> et <code>p_value = 6.14e-06<\/code>. La p-valeur se situant infiniment en de\u00e7\u00e0 du seuil critique &alpha; = 0,05, nous sommes fond\u00e9s \u00e0 rejeter l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9quivalence des distributions de rangs entre les trois conditions pharmacologiques.<\/p>\n<h3>6.2 Alternative via la biblioth\u00e8que Pingouin<\/h3>\n<p>Bien que SciPy soit extr\u00eamement efficace, son rapport de sortie est minimaliste, ne retournant que deux scalaires bruts. Pour les chercheurs en psychologie soucieux d&rsquo;obtenir une restitution m\u00e9thodologique exhaustive incluant les degr\u00e9s de libert\u00e9 et les indicateurs de magnitude de l&rsquo;effet, la biblioth\u00e8que <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pingouin-stats.org\">Pingouin<\/a> s&rsquo;av\u00e8re particuli\u00e8rement adapt\u00e9e. Sa fonction <code>friedman()<\/code> exploite directement le format long (tidy) des donn\u00e9es :<\/p>\n<p><code>res_friedman = pg.friedman(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;data=data_long,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;dv='Temps_Reaction',<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;within='Condition',<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;subject='Sujet'<br \/>\n)<\/code><\/p>\n<p>Le DataFrame structur\u00e9 r\u00e9sultant fournit instantan\u00e9ment la statistique de test <em>Q<\/em>, les degr\u00e9s de libert\u00e9 (<em>ddl = 2<\/em>), la p-valeur exacte, mais \u00e9galement le coefficient de concordance <em>W<\/em> de Kendall. Le <em>W<\/em> de Kendall constitue la taille d&rsquo;effet standard associ\u00e9e au test de Friedman, variant rigoureusement de 0 (absence totale de concordance des rangs entre sujets) \u00e0 1 (concordance absolue, chaque sujet classant les traitements exactement dans le m\u00eame ordre relatif). Obtenir un <em>W<\/em> sup\u00e9rieur \u00e0 0,70 d\u00e9note une taille d&rsquo;effet substantielle, fournissant un contexte clinique indispensable que la seule p-valeur ne saurait refl\u00e9ter.<\/p>\n<h3>6.3 Validation formelle de la condition d&rsquo;encha\u00eenement post-hoc<\/h3>\n<p>La r\u00e8gle m\u00e9thodologique pr\u00e9sidant \u00e0 la transition vers les tests post-hoc doit \u00eatre formalis\u00e9e avec une rigueur absolue au sein du pipeline d&rsquo;analyse computationnelle. Un embranchement logique strict doit \u00eatre programm\u00e9 :<\/p>\n<p><code>if p_value &lt; 0.05:<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(\"Test de Friedman significatif (p &lt; 0.05). Ex\u00e9cution l\u00e9gitime du test de Nemenyi autoris\u00e9e.\")<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;# Appel de la proc\u00e9dure post-hoc<br \/>\nelse:<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(\"Test de Friedman non significatif. R\u00e9tention de H0. Les contrastes post-hoc sont proscrits.\")<\/code><\/p>\n<p>Si la statistique de Friedman n&rsquo;atteint pas le seuil de rejet fix\u00e9, l&rsquo;enqu\u00eate inf\u00e9rentielle s&rsquo;arr\u00eate imp\u00e9rativement \u00e0 ce stade. Tenter de sonder a posteriori des comparaisons par paires en d\u00e9pit d&rsquo;un test omnibus non significatif constitue une entorse m\u00e9thodologique majeure, exposant l&rsquo;analyste \u00e0 la d\u00e9tection d&rsquo;artefacts sans aucune validit\u00e9 populationnelle.<\/p>\n<h2>7. Impl\u00e9mentation pas \u00e0 pas du test post-hoc de Nemenyi avec scikit-posthocs<\/h2>\n<h3>7.1 Choix de la fonction appropri\u00e9e dans scikit-posthocs<\/h3>\n<p>La biblioth\u00e8que sp\u00e9cialis\u00e9e <code>scikit-posthocs<\/code> concentre une collection \u00e9tendue d&rsquo;algorithmes post-hoc non param\u00e9triques. Il est crucial pour le praticien de s\u00e9lectionner la fonction rigoureusement adapt\u00e9e au design \u00e0 mesures r\u00e9p\u00e9t\u00e9es en blocs. Deux fonctions portent des d\u00e9nominations proches mais r\u00e9pondent \u00e0 des paradigmes exp\u00e9rimentaux radicalement diff\u00e9rents :<\/p>\n<ul>\n<li><code>posthoc_nemenyi()<\/code> : Cette fonction est d\u00e9di\u00e9e aux comparaisons multiples pour <em>\u00e9chantillons ind\u00e9pendants<\/em>, venant se greffer \u00e0 la suite d&rsquo;un test omnibus de Kruskal-Wallis. Elle classe les donn\u00e9es de fa\u00e7on conjointe sur l&rsquo;ensemble des groupes.<\/li>\n<li><code>posthoc_nemenyi_friedman()<\/code> : C&rsquo;est la fonction sp\u00e9cifique d\u00e9di\u00e9e aux <em>mesures r\u00e9p\u00e9t\u00e9es appari\u00e9es<\/em> (blocs complets), constituant le compl\u00e9ment direct du test de Friedman. Elle applique le classement intra-sujet ligne par ligne et mobilise la distribution de l&rsquo;\u00e9tendue studentis\u00e9e adapt\u00e9e aux blocs d\u00e9pendants.<\/li>\n<\/ul>\n<p>L&rsquo;utilisation erron\u00e9e de <code>posthoc_nemenyi<\/code> \u00e0 la place de <code>posthoc_nemenyi_friedman<\/code> biaiserait enti\u00e8rement le calcul en n\u00e9gligeant l&rsquo;appariement des sujets, conduisant \u00e0 une sous-estimation dramatique de la puissance du test en ignorant la suppression de la variance interindividuelle.<\/p>\n<h3>7.2 Ex\u00e9cution du test sur les donn\u00e9es exp\u00e9rimentales<\/h3>\n<p>L&rsquo;ex\u00e9cution de la fonction <code>posthoc_nemenyi_friedman<\/code> sur nos donn\u00e9es exp\u00e9rimentales est d&rsquo;une grande fluidit\u00e9. La fonction peut ing\u00e9rer une matrice NumPy bidimensionnelle, un DataFrame au format large contenant exclusivement les colonnes des conditions, ou un tableau au format long via la sp\u00e9cification explicite des variables.<\/p>\n<p>En exploitant notre matrice de donn\u00e9es large, le code s&rsquo;exprime ainsi :<\/p>\n<p><code>matrice_conditions = data_wide[['Placebo', 'Molecule_A', 'Molecule_B']].to_numpy()<br \/>\nresultats_nemenyi = sp.posthoc_nemenyi_friedman(matrice_conditions)<\/code><\/p>\n<p>Pour assurer une parfaite lisibilit\u00e9 de la matrice de sortie, nous indexons les lignes et les colonnes avec les d\u00e9nominations textuelles explicites de nos traitements exp\u00e9rimentaux :<\/p>\n<p><code>conditions_labels = ['Placebo', 'Molecule_A', 'Molecule_B']<br \/>\nresultats_nemenyi.index = conditions_labels<br \/>\nresultats_nemenyi.columns = conditions_labels<br \/>\nprint(resultats_nemenyi)<\/code><\/p>\n<p>L&rsquo;objet retourn\u00e9 est un DataFrame Pandas carr\u00e9 sym\u00e9trique de dimensions <em>k &times; k<\/em>. Chaque cellule <em>(i, j)<\/em> contient la p-valeur ajust\u00e9e exacte de Nemenyi associ\u00e9e au contraste par paire entre la condition de la ligne <em>i<\/em> et la condition de la colonne <em>j<\/em>. La diagonale principale, confrontant chaque condition \u00e0 elle-m\u00eame, est conventionnellement garnie de valeurs unitaires (1.000).<\/p>\n<h3>7.3 Calcul manuel alternatif des rangs pour validation<\/h3>\n<p>Afin de d\u00e9mystifier le fonctionnement de la bo\u00eete noire algorithmique et de consolider la compr\u00e9hension m\u00e9trologique du processus, il est hautement instructif de programmer manuellement sous Python le calcul des rangs moyens et l&rsquo;\u00e9valuation de la diff\u00e9rence critique. Cette d\u00e9composition analytique s&rsquo;effectue en quelques lignes de calcul vectoriel transparentes :<\/p>\n<p><code># 1. Extraction et transformation des scores bruts en rangs intra-sujets<br \/>\n# L'argument axis=1 ordonne les scores le long de chaque ligne (intra-sujet)<br \/>\nrangs_individuels = stats.rankdata(matrice_conditions, axis=1)<\/p>\n<p># 2. Calcul des rangs moyens par colonne (condition)<br \/>\nrangs_moyens = np.mean(rangs_individuels, axis=0)<br \/>\nfor cond, rm in zip(conditions_labels, rangs_moyens):<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(f\"Rang moyen - {cond} : {rm:.4f}\")<\/p>\n<p># 3. Param\u00e8tres exp\u00e9rimentaux pour le calcul de CD<br \/>\nN_sujets = matrice_conditions.shape[0] # N = 12<br \/>\nk_traitements = matrice_conditions.shape[1] # k = 3<br \/>\nalpha = 0.05<\/p>\n<p># 4. Extraction de la valeur critique de l'\u00e9tendue studentis\u00e9e q<br \/>\n# Pour k = 3 et ddl infinis au seuil 0.05, q_alpha est approximativement de 3.314<br \/>\n# Dans scipy.stats, la distribution studentis\u00e9e est disponible via stats.studentized_range<br \/>\nq_val = stats.studentized_range.ppf(1 - alpha, k_traitements, np.inf)<\/p>\n<p># 5. Calcul analytique de la Diff\u00e9rence Critique (CD)<br \/>\necart_type_rangs = np.sqrt((k_traitements * (k_traitements + 1)) \/ (6 * N_sujets))<br \/>\ndifference_critique = q_val * ecart_type_rangs \/ np.sqrt(2)<br \/>\nprint(f\"Diff\u00e9rence Critique (CD) calcul\u00e9e : {difference_critique:.4f}\")<\/code><\/p>\n<p>En confrontant directement les diff\u00e9rences arithm\u00e9tiques observ\u00e9es entre les rangs moyens calcul\u00e9s (par exemple, <em>| R&#772;<sub>Placebo<\/sub> &#8211; R&#772;<sub>Molecule_A<\/sub> |<\/em>) \u00e0 la valeur de la diff\u00e9rence critique <em>CD<\/em>, l&rsquo;analyste retrouve exactement le verdict binaire fourni par les p-valeurs de <code>scikit-posthocs<\/code>. Si l&rsquo;\u00e9cart observ\u00e9 franchit le seuil <em>CD<\/em>, la p-valeur associ\u00e9e dans la matrice est strictement inf\u00e9rieure \u00e0 0,05.<\/p>\n<h2>8. Interpr\u00e9tation statistique et matrices des p-valeurs<\/h2>\n<h3>8.1 Lecture et d\u00e9cryptage de la matrice sym\u00e9trique<\/h3>\n<p>L&rsquo;analyse approfondie de la matrice de sortie issue du test de Nemenyi constitue l&rsquo;\u00e9tape centrale de la prise de d\u00e9cision statistique. Consid\u00e9rons les r\u00e9sultats empiriques fictifs obtenus \u00e0 la suite de notre exp\u00e9rimentation psychopharmacologique sur le temps de r\u00e9action :<\/p>\n<table border=\"1\" cellpadding=\"8\" style=\"border-collapse: collapse; text-align: center; margin: 20px 0; width: 100%;\">\n<tr style=\"background-color: #f2f2f2;\">\n<th>Condition<\/th>\n<th>Placebo<\/th>\n<th>Molecule_A<\/th>\n<th>Molecule_B<\/th>\n<\/tr>\n<tr>\n<td style=\"font-weight: bold; background-color: #f9f9f9;\">Placebo<\/td>\n<td>1.000000<\/td>\n<td>0.000120<\/td>\n<td>0.041500<\/td>\n<\/tr>\n<tr>\n<td style=\"font-weight: bold; background-color: #f9f9f9;\">Molecule_A<\/td>\n<td>0.000120<\/td>\n<td>1.000000<\/td>\n<td>0.089200<\/td>\n<\/tr>\n<tr>\n<td style=\"font-weight: bold; background-color: #f9f9f9;\">Molecule_B<\/td>\n<td>0.041500<\/td>\n<td>0.089200<\/td>\n<td>1.000000<\/td>\n<\/tr>\n<\/table>\n<p>La lecture de ce tableau s&rsquo;op\u00e8re par intersection ligne-colonne. Trois contrastes par paires distincts structurent l&rsquo;espace d\u00e9cisionnel :<\/p>\n<ul>\n<li><strong>Placebo versus Mol\u00e9cule A :<\/strong> La p-valeur ajust\u00e9e s&rsquo;\u00e9tablit \u00e0 <em>p = 0,00012<\/em>. \u00c9tant largement inf\u00e9rieure au seuil conventionnel de 0,05, nous concluons \u00e0 une diff\u00e9rence hautement significative entre le placebo et la Mol\u00e9cule A. L&rsquo;inspection des rangs moyens montre que la Mol\u00e9cule A induit des temps de r\u00e9action syst\u00e9matiquement plus rapides.<\/li>\n<li><strong>Placebo versus Mol\u00e9cule B :<\/strong> La p-valeur ajust\u00e9e s&rsquo;\u00e9l\u00e8ve \u00e0 <em>p = 0,0415<\/em>. Ce contraste franchit le seuil critique de significativit\u00e9 (<em>p &lt; 0,05<\/em>). Nous rejetons l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9quivalence : la Mol\u00e9cule B d\u00e9montre une sup\u00e9riorit\u00e9 statistiquement \u00e9tablie sur le placebo dans l&rsquo;acc\u00e9l\u00e9ration du traitement de l&rsquo;information.<\/li>\n<li><strong>Mol\u00e9cule A versus Mol\u00e9cule B :<\/strong> La p-valeur ajust\u00e9e atteint <em>p = 0,0892<\/em>. Se situant au-dessus du seuil critique de 0,05, ce contraste \u00e9choue \u00e0 atteindre la significativit\u00e9 statistique. Bien qu&rsquo;une tendance num\u00e9rique en faveur de la Mol\u00e9cule A puisse \u00eatre observ\u00e9e sur les rangs moyens bruts, les donn\u00e9es ne permettent pas d&rsquo;affirmer avec un niveau de confiance suffisant une divergence d&rsquo;efficacit\u00e9 entre ces deux agents actifs.<\/li>\n<\/ul>\n<h3>8.2 Quantification de la taille d&rsquo;effet post-hoc<\/h3>\n<p>Dans l&rsquo;esprit des recommandations de l&rsquo;American Psychological Association, la communication d&rsquo;une p-valeur isol\u00e9e est fonci\u00e8rement incompl\u00e8te si elle ne s&rsquo;accompagne pas d&rsquo;un indicateur de magnitude de l&rsquo;effet. En contexte post-hoc non param\u00e9trique appari\u00e9, plusieurs m\u00e9triques de taille d&rsquo;effet peuvent \u00eatre d\u00e9riv\u00e9es :<\/p>\n<p>La premi\u00e8re r\u00e9side dans la diff\u00e9rence standardis\u00e9e des rangs moyens relative \u00e0 l&rsquo;\u00e9tendue maximale th\u00e9orique :<\/p>\n<p><strong>&Delta;R&#772;<sub>norm<\/sub> = | R&#772;<sub>A<\/sub> &#8211; R&#772;<sub>B<\/sub> | \/ (k &#8211; 1)<\/strong><\/p>\n<p>Cette valeur borne le d\u00e9placement d&rsquo;amplitude ordinale entre 0 et 1, facilitant la comparaison transversale entre des protocoles comportant des nombres divergents de traitements.<\/p>\n<p>Une seconde approche, particuli\u00e8rement appr\u00e9ci\u00e9e en clinique comportementale, consiste \u00e0 extraire pour chaque paire significative le coefficient de corr\u00e9lation de rang bis\u00e9riel appari\u00e9 (ou le <em>r<\/em> de Rosenthal d\u00e9riv\u00e9 de la statistique de Wilcoxon <em>r = Z \/ &radic;N<\/em>). Un indice <em>r<\/em> oscillant autour de 0,10 d\u00e9signe un effet modeste, 0,30 un effet moyen, et au-del\u00e0 de 0,50 un effet clinique majeur. Cette distinction est cruciale : une diff\u00e9rence peut se r\u00e9v\u00e9ler statistiquement significative au sens de Nemenyi en raison d&rsquo;une cohorte <em>N<\/em> tr\u00e8s importante sans pour autant rev\u00eatir une quelconque pertinence pratique pour le patient au quotidien.<\/p>\n<h3>8.3 Formulation acad\u00e9mique des conclusions statistiques<\/h3>\n<p>La r\u00e9daction de la section R\u00e9sultats au sein d&rsquo;un manuscrit scientifique soumis aux normes de l&rsquo;APA (7<sup>e<\/sup> \u00e9dition) impose une standardisation formelle stricte, m\u00ealant synth\u00e8se descriptive (m\u00e9dianes, espaces interquartiles, rangs moyens) et m\u00e9triques inf\u00e9rentielles exactes.<\/p>\n<p>Voici un mod\u00e8le r\u00e9dactionnel acad\u00e9mique adapt\u00e9 \u00e0 notre cas d&rsquo;\u00e9tude :<\/p>\n<blockquote>\n<p><em>&laquo; Les temps de r\u00e9action m\u00e9dians lors de la t\u00e2che d&rsquo;attention soutenue ont \u00e9t\u00e9 soumis \u00e0 un test d&rsquo;analyse de variance par rangs de Friedman pour mesures r\u00e9p\u00e9t\u00e9es. L&rsquo;analyse r\u00e9v\u00e8le un effet global hautement significatif du traitement pharmacologique, &chi;&sup2;<sub>F<\/sub>(2) = 24,00, p &lt; 0,001, avec une excellente concordance des rangs inter-sujets (W de Kendall = 0,75). Des analyses post-hoc par paires ont \u00e9t\u00e9 conduites en appliquant le test de Nemenyi pour blocs complets avec contr\u00f4le rigoureux du taux d&rsquo;erreur par famille (FWER). La diff\u00e9rence critique requise pour atteindre le seuil de significativit\u00e9 \u00e0 &alpha; = 0,05 a \u00e9t\u00e9 \u00e9tablie \u00e0 CD = 0,782. Les contrastes ont d\u00e9montr\u00e9 que les temps de r\u00e9action \u00e9taient significativement plus courts sous Mol\u00e9cule A (Rang moyen = 1,17 ; M\u00e9diane = 405 ms, IQR = 42,5) en comparaison du Placebo (Rang moyen = 2,92 ; M\u00e9diane = 485 ms, IQR = 57,5 ; p &lt; 0,001). De m\u00eame, la condition Mol\u00e9cule B (Rang moyen = 1,92 ; M\u00e9diane = 445 ms, IQR = 52,5) a induit une r\u00e9duction statistiquement significative des temps de r\u00e9action comparativement au Placebo (p = 0,042). En revanche, la divergence observ\u00e9e entre la Mol\u00e9cule A et la Mol\u00e9cule B n&rsquo;a pas atteint le seuil critique de significativit\u00e9 statistique (p = 0,089), sugg\u00e9rant une efficacit\u00e9 attentionnelle comparable entre ces deux mol\u00e9cules au sein de la cohorte \u00e9valu\u00e9e. &raquo;<\/em><\/p>\n<\/blockquote>\n<h2>9. Visualisation graphique avanc\u00e9e des r\u00e9sultats de Nemenyi<\/h2>\n<h3>9.1 Conception de cartes thermiques (heatmaps) de significativit\u00e9<\/h3>\n<p>La matrice num\u00e9rique brute des p-valeurs peut devenir difficile \u00e0 appr\u00e9hender instantan\u00e9ment d\u00e8s lors que le nombre de traitements exp\u00e9rimentaux s&rsquo;accro\u00eet. La g\u00e9n\u00e9ration d&rsquo;une carte thermique (heatmap) personnalis\u00e9e constitue une solution d&rsquo;une remarquable \u00e9l\u00e9gance visuelle. La biblioth\u00e8que <code>seaborn<\/code>, associ\u00e9e aux utilitaires de coloration de <code>matplotlib<\/code>, permet d&rsquo;isoler imm\u00e9diatement les contrastes significatifs au moyen de palettes divergentes ou s\u00e9quentielles adapt\u00e9es.<\/p>\n<p>Le code suivant construit une figure pr\u00eate pour la publication, masquant la partie sup\u00e9rieure redondante de la matrice sym\u00e9trique et annotant explicitement les seuils d\u00e9cisionnels :<\/p>\n<p><code>plt.figure(figsize=(8, 6), dpi=300)<br \/>\n# Cr\u00e9ation d'un masque pour dissimuler le triangle sup\u00e9rieur et la diagonale<br \/>\nmask = np.triu(np.ones_like(resultats_nemenyi, dtype=bool))<\/p>\n<p># Configuration d'une carte de couleurs conditionnelle centr\u00e9e sur 0.05<br \/>\ncmap = sns.diverging_palette(10, 130, as_cmap=True)<\/p>\n<p>sns.heatmap(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;resultats_nemenyi,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;mask=mask,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;annot=True,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;fmt=\".4f\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;cmap='vlag_r',<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;vmin=0.0,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;vmax=0.10,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;cbar_kws={'label': 'P-valeur ajust\u00e9e (Nemenyi)'},<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;square=True,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;linewidths=1.5<br \/>\n)<br \/>\nplt.title(\"Matrice post-hoc de Nemenyi : Comparaisons multiples\", fontsize=14, pad=15)<br \/>\nplt.tight_layout()<br \/>\nplt.savefig(\"heatmap_nemenyi.png\", dpi=300)<\/code><\/p>\n<p>Ce type de repr\u00e9sentation offre aux relecteurs et aux lecteurs une lecture imm\u00e9diate : les cellules adoptent une chromaticit\u00e9 vive pour les p-valeurs inf\u00e9rieures \u00e0 0,05, tandis que les contrastes non significatifs basculent vers des tonalit\u00e9s neutres ou froides.<\/p>\n<h3>9.2 \u00c9laboration d&rsquo;un diagramme de diff\u00e9rence critique (CD Diagram)<\/h3>\n<p>Dans l&rsquo;histoire r\u00e9cente de la science des donn\u00e9es et de l&rsquo;\u00e9valuation statistique, le chercheur slov\u00e8ne Janez Dem\u0161ar a publi\u00e9 en 2006 un article fondateur dans le <em>Journal of Machine Learning Research<\/em> intitul\u00e9 <em>Statistical Comparisons of Classifiers over Multiple Data Sets<\/em>. Dans cette contribution majeure, Dem\u0161ar a popularis\u00e9 le diagramme de diff\u00e9rence critique (Critical Difference Diagram), devenu depuis lors le standard graphique incontournable pour illustrer le test de Nemenyi.<\/p>\n<p>Le diagramme de Dem\u0161ar projette l&rsquo;ensemble des conditions exp\u00e9rimentales le long d&rsquo;un axe horizontal gradu\u00e9 repr\u00e9sentant l&rsquo;amplitude des rangs moyens (de 1 \u00e0 <em>k<\/em>). Chaque traitement est positionn\u00e9 sur cet axe \u00e0 l&#8217;emplacement exact de son rang moyen empirique. Un segment horizontal de longueur fixe, correspondant rigoureusement \u00e0 la distance critique <em>CD<\/em>, est ensuite trac\u00e9 au-dessus des groupes. Les conditions dont les rangs moyens ne diff\u00e8rent pas de mani\u00e8re statistiquement significative (c&rsquo;est-\u00e0-dire dont l&rsquo;\u00e9cartement spatial est inf\u00e9rieur \u00e0 <em>CD<\/em>) sont interconnect\u00e9es visuellement par une barre horizontale \u00e9paisse de non-diff\u00e9rence.<\/p>\n<p>La biblioth\u00e8que <code>scikit-posthocs<\/code> dispose d&rsquo;une fonction native permettant de g\u00e9n\u00e9rer ces visualisations vectorielles complexes :<\/p>\n<p><code># Dictionnaire associant chaque traitement \u00e0 son rang moyen calcul\u00e9<br \/>\ndict_rangs = dict(zip(conditions_labels, rangs_moyens))<\/p>\n<p>plt.figure(figsize=(10, 3), dpi=300)<br \/>\nsp.critical_difference_diagram(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;ranks=dict_rangs,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;sig_matrix=resultats_nemenyi,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;label_fmt_hspace=0.8<br \/>\n)<br \/>\nplt.title(f\"Diagramme de Diff\u00e9rence Critique (CD = {difference_critique:.3f}, alpha = 0.05)\", pad=20)<br \/>\nplt.savefig(\"cd_diagram_nemenyi.pdf\", format='pdf', bbox_inches='tight')<\/code><\/p>\n<p>Sur le graphique r\u00e9sultant, l&rsquo;absence de ligne de connexion entre le Placebo et la Mol\u00e9cule A mat\u00e9rialise instantan\u00e9ment leur divergence significative, tandis qu&rsquo;une barre horizontale reliant la Mol\u00e9cule A et la Mol\u00e9cule B formalise visuellement leur \u00e9quivalence statistique.<\/p>\n<h3>9.3 Graphiques combin\u00e9s des distributions et rangs individuels<\/h3>\n<p>Si le diagramme de diff\u00e9rence critique r\u00e9sume admirablement le r\u00e9sultat inf\u00e9rentiel, il masque toutefois la dispersion empirique sous-jacente des scores individuels. En psychologie exp\u00e9rimentale, il est hautement recommand\u00e9 de coupler la visualisation des rangs avec une repr\u00e9sentation des distributions r\u00e9elles. Une strat\u00e9gie graphique de premier ordre consiste \u00e0 superposer un diagramme en bo\u00eete (boxplot) ou en violon (violin plot) avec un trac\u00e9 de coordonn\u00e9es parall\u00e8les illustrant les trajectoires intra-sujets de chaque participant :<\/p>\n<p><code>plt.figure(figsize=(9, 6), dpi=300)<br \/>\n# Bo\u00eete \u00e0 moustaches descriptive \u00e9pur\u00e9e<br \/>\nsns.boxplot(data=data_long, x='Condition', y='Temps_Reaction', palette='Pastel1', width=0.4, boxprops=dict(alpha=0.6))<\/p>\n<p># Superposition des trajectoires individuelles (mesures r\u00e9p\u00e9t\u00e9es)<br \/>\nfor sujet_id in data_wide['Sujet']:<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;sujet_vals = data_wide[data_wide['Sujet'] == sujet_id][conditions_labels].values.flatten()<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;plt.plot(conditions_labels, sujet_vals, color='gray', alpha=0.5, linestyle='--', marker='o', markersize=4)<\/p>\n<p>plt.ylabel(\"Temps de R\u00e9action m\u00e9dian (ms)\", fontsize=12)<br \/>\nplt.xlabel(\"Condition pharmacologique\", fontsize=12)<br \/>\nplt.title(\"Distributions et trajectoires individuelles sous les trois conditions\", fontsize=14)<br \/>\nplt.tight_layout()<br \/>\nplt.savefig(\"trajectoires_distributions.svg\", format='svg')<\/code><\/p>\n<p>L&rsquo;exportation au format vectoriel (SVG ou PDF) garantit une nettet\u00e9 absolue des figures lors de leur int\u00e9gration dans des syst\u00e8mes de traitement \u00e9ditoriaux tels que LaTeX ou des suites de publication professionnelles.<\/p>\n<h2>10. Gestion des biais statistiques et corrections m\u00e9thodologiques<\/h2>\n<h3>10.1 Probl\u00e9matique du conservatisme inh\u00e9rent au test de Nemenyi<\/h3>\n<p>Bien que le test post-hoc de Nemenyi b\u00e9n\u00e9ficie d&rsquo;une solide assise th\u00e9orique gr\u00e2ce \u00e0 son contr\u00f4le irr\u00e9prochable du taux d&rsquo;erreur global de type I (FWER), il souffre d&rsquo;un d\u00e9faut m\u00e9thodologique bien document\u00e9 dans la litt\u00e9rature statistique sp\u00e9cialis\u00e9e : son conservatisme prononc\u00e9. Ce ph\u00e9nom\u00e8ne s&rsquo;intensifie de mani\u00e8re quasi exponentielle \u00e0 mesure que le nombre de traitements exp\u00e9rimentaux <em>k<\/em> s&rsquo;accro\u00eet.<\/p>\n<p>Ce conservatisme d\u00e9coule directement de la mobilisation de la distribution de l&rsquo;\u00e9tendue studentis\u00e9e. En postulant implicitement que toutes les paires sont \u00e9valu\u00e9es sous la pire configuration possible d&rsquo;\u00e9cartement maximal, la diff\u00e9rence critique s&rsquo;\u00e9largit consid\u00e9rablement d\u00e8s que <em>k<\/em> d\u00e9passe 5 ou 6 conditions. La cons\u00e9quence directe pour le chercheur est une d\u00e9gradation sensible de la puissance statistique et un accroissement marqu\u00e9 du risque d&rsquo;erreur de type II (b\u00eata) : le test \u00e9choue fr\u00e9quemment \u00e0 rejeter l&rsquo;hypoth\u00e8se nulle pour des paires de conditions pr\u00e9sentant pourtant des divergences empiriques r\u00e9elles. Lorsque les cohortes cliniques sont modestes (par exemple <em>N &lt; 15<\/em>, situation courante en neuropsychologie lourde ou en imagerie fonctionnelle), le test de Nemenyi peut se r\u00e9v\u00e9ler incapable de d\u00e9tecter la moindre diff\u00e9rence par paire, bien que le test de Friedman omnibus ait \u00e9t\u00e9 hautement significatif.<\/p>\n<h3>10.2 Comparaison avec les corrections de Holm, Bonferroni et Benjamini-Hochberg<\/h3>\n<p>Face \u00e0 la contrainte de ce conservatisme, l&rsquo;investigateur doit arbitrer entre plusieurs philosophies de contr\u00f4le statistique. Deux grandes familles de proc\u00e9dures correctrices s&rsquo;offrent alors pour mener les comparaisons par paires cons\u00e9cutives \u00e0 Friedman :<\/p>\n<ul>\n<li><strong>Le contr\u00f4le du FWER par la proc\u00e9dure s\u00e9quentielle de Holm :<\/strong> Plut\u00f4t que d&rsquo;appliquer la diff\u00e9rence critique fixe de Nemenyi, l&rsquo;analyste ex\u00e9cute des tests de Wilcoxon appari\u00e9s sur chaque paire, puis soumet les p-valeurs \u00e0 la correction step-down de Holm-Bonferroni. La m\u00e9thode de Holm est math\u00e9matiquement sup\u00e9rieure \u00e0 la correction classique de Bonferroni en tout point : elle garantit un contr\u00f4le tout aussi strict du FWER \u00e0 0,05 mais offre une puissance statistique substantiellement sup\u00e9rieure, car le seuil d&rsquo;ajustement s&rsquo;assouplit progressivement apr\u00e8s chaque rejet successif.<\/li>\n<li><strong>Le contr\u00f4le du taux de fausses d\u00e9couvertes (False Discovery Rate, FDR) :<\/strong> Conceptualis\u00e9 par Yoav Benjamini et Yosef Hochberg (1995), le FDR ne cherche pas \u00e0 garantir l&rsquo;absence totale de tout faux positif parmi l&rsquo;ensemble des tests, mais r\u00e9gule la proportion attendue d&rsquo;erreurs parmi les rejets d\u00e9clar\u00e9s. En psychologie exploratoire ou en g\u00e9nomique comportementale o\u00f9 des dizaines de conditions sont cribl\u00e9es, substituer le test de Nemenyi par une s\u00e9rie de comparaisons corrig\u00e9es par Benjamini-Hochberg permet d&rsquo;\u00e9viter l&rsquo;\u00e9cueil d&rsquo;une h\u00e9catombe de puissance statistique.<\/li>\n<\/ul>\n<p>Dans <code>scikit-posthocs<\/code>, ces alternatives sont imm\u00e9diatement accessibles via la fonction <code>posthoc_wilcoxon()<\/code> en sp\u00e9cifiant l&rsquo;argument d&rsquo;ajustement ad\u00e9quat : <code>p_adjust='holm'<\/code> ou <code>p_adjust='fdr_bh'<\/code>.<\/p>\n<h3>10.3 Gestion des d\u00e9s\u00e9quilibres et des cas de puissance statistique marginale<\/h3>\n<p>La pr\u00e9cision de la statistique de diff\u00e9rence critique <em>CD<\/em> d\u00e9pend intimement du volume d&rsquo;\u00e9chantillonnage <em>N<\/em>. \u00c9tant donn\u00e9 que le terme d&rsquo;erreur est inversement proportionnel \u00e0 &radic;<em>N<\/em>, les cohortes restreintes induisent une incertitude m\u00e9trologique substantielle. Dans les \u00e9tudes pilotes en psychoth\u00e9rapie ou en r\u00e9\u00e9ducation cognitive post-AVC, o\u00f9 <em>N<\/em> n&rsquo;exc\u00e8de parfois pas 8 ou 10 patients, le seuil de Nemenyi devient particuli\u00e8rement exigeant.<\/p>\n<p>Dans de telles configurations de puissance marginale, les bonnes pratiques m\u00e9thodologiques recommandent :<\/p>\n<ul>\n<li>D&rsquo;\u00e9viter formellement toute tentative de compensation par le rel\u00e2chement arbitraire du seuil alpha (par exemple en tol\u00e9rant un seuil exploratoire \u00e0 0,10 sans justification th\u00e9orique explicite).<\/li>\n<li>De d\u00e9clarer syst\u00e9matiquement les intervalles de confiance associ\u00e9s aux rangs moyens au moyen de techniques de r\u00e9\u00e9chantillonnage non param\u00e9trique (bootstrap intra-sujet).<\/li>\n<li>De consid\u00e9rer les r\u00e9sultats non comme des verdicts confirmatifs d\u00e9finitifs, mais comme des signaux exploratoires justifiant une r\u00e9plication pr\u00e9-enregistr\u00e9e (preregistration) sur une cohorte dimensionn\u00e9e a priori \u00e0 l&rsquo;aide d&rsquo;un calcul de puissance prospective.<\/li>\n<\/ul>\n<h2>11. Cas pratique complet : \u00c9valuation de trois interventions psychoth\u00e9rapeutiques<\/h2>\n<h3>11.1 Contexte exp\u00e9rimental et cr\u00e9ation du jeu de donn\u00e9es<\/h3>\n<p>Afin de synth\u00e9tiser l&rsquo;ensemble des concepts abord\u00e9s au sein d&rsquo;une d\u00e9monstration op\u00e9rationnelle int\u00e9grale, d\u00e9veloppons un cas pratique appliqu\u00e9 en psychologie clinique. Une \u00e9quipe de recherche s&rsquo;int\u00e9resse \u00e0 l&rsquo;efficacit\u00e9 diff\u00e9rentielle de trois interventions psychoth\u00e9rapeutiques br\u00e8ves destin\u00e9es \u00e0 r\u00e9duire l&rsquo;anxi\u00e9t\u00e9 g\u00e9n\u00e9ralis\u00e9e chez des \u00e9tudiants universitaires en situation d&rsquo;\u00e9chec acad\u00e9mique. L&rsquo;\u00e9chantillon comprend <em>N = 15<\/em> participants soumis successivement, au terme de trois semestres distincts et selon une randomisation crois\u00e9e compl\u00e8te, aux trois modules th\u00e9rapeutiques suivants :<\/p>\n<ul>\n<li><strong>Th\u00e9rapie Comportementale et Cognitive (TCC) :<\/strong> Focalis\u00e9e sur la restructuration cognitive et l&rsquo;exposition aux pens\u00e9es automatiques anxiog\u00e8nes.<\/li>\n<li><strong>M\u00e9ditation de Pleine Conscience (Mindfulness) :<\/strong> Ax\u00e9e sur l&rsquo;acceptation \u00e9motionnelle et la r\u00e9gulation attentionnelle non jugeante.<\/li>\n<li><strong>Th\u00e9rapie Interpersonnelle (TIP) :<\/strong> Centr\u00e9e sur la r\u00e9solution des conflits relationnels et le soutien social per\u00e7u.<\/li>\n<\/ul>\n<p>Au terme de chaque cycle d&rsquo;intervention de six semaines, la s\u00e9v\u00e9rit\u00e9 de l&rsquo;anxi\u00e9t\u00e9 est \u00e9valu\u00e9e au moyen de l&rsquo;\u00e9chelle d&rsquo;anxi\u00e9t\u00e9 de Hamilton (HAM-A), un instrument psychom\u00e9trique h\u00e9t\u00e9ro-\u00e9valu\u00e9 fournissant un score global variant de 0 \u00e0 56. Des scores \u00e9lev\u00e9s signalent une symptomatologie clinique s\u00e9v\u00e8re.<\/p>\n<h3>11.2 Pipeline complet d&rsquo;analyse automatis\u00e9e en Python<\/h3>\n<p>Le script Python ci-dessous constitue un pipeline d&rsquo;analyse de bout en bout, englobant l&rsquo;initialisation structur\u00e9e des donn\u00e9es cliniques, le contr\u00f4le exploratoire, l&rsquo;ex\u00e9cution automatis\u00e9e du test omnibus de Friedman, le branchement conditionnel vers le test de Nemenyi, et la mise en forme tabulaire des r\u00e9sultats :<\/p>\n<p><code># Script complet d'analyse psychom\u00e9trique automatis\u00e9e<br \/>\nimport numpy as np<br \/>\nimport pandas as pd<br \/>\nimport scipy.stats as stats<br \/>\nimport scikit_posthocs as sp<\/p>\n<p># 1. G\u00e9n\u00e9ration du jeu de donn\u00e9es cliniques (N = 15 sujets)<br \/>\nnp.random.seed(42)<br \/>\nscores_tcc = [14, 18, 12, 16, 20, 15, 11, 19, 13, 17, 15, 12, 16, 14, 18]<br \/>\nscores_mindfulness = [22, 25, 19, 24, 28, 21, 18, 26, 20, 23, 22, 19, 25, 21, 24]<br \/>\nscores_tip = [26, 29, 23, 27, 31, 25, 22, 30, 24, 28, 27, 23, 29, 26, 28]<\/p>\n<p>df_clinique = pd.DataFrame({<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'Patient': [f'P{i:02d}' for i in range(1, 16)],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'TCC': scores_tcc,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'Mindfulness': scores_mindfulness,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;'TIP': scores_tip<br \/>\n})<\/p>\n<p># 2. Statistiques descriptives exploratoires (M\u00e9diane et \u00c9cart Interquartile)<br \/>\ndesc = df_clinique[['TCC', 'Mindfulness', 'TIP']].apply(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;lambda x: pd.Series({'M\u00e9diane': np.median(x), 'IQR': stats.iqr(x)})<br \/>\n)<br \/>\nprint(\"=== STATISTIQUES DESCRIPTIVES (HAM-A) ===\")<br \/>\nprint(desc.round(2))<\/p>\n<p># 3. Ex\u00e9cution du test omnibus de Friedman<br \/>\nchi2_f, p_friedman = stats.friedmanchisquare(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;df_clinique['TCC'],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;df_clinique['Mindfulness'],<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;df_clinique['TIP']<br \/>\n)<br \/>\nprint(f\"n=== TEST OMNIBUS DE FRIEDMAN ===\")<br \/>\nprint(f\"Statistique Chi-deux Q : {chi2_f:.4f}\")<br \/>\nprint(f\"P-valeur asymptotique : {p_friedman:.4e}\")<\/p>\n<p># 4. Branchement d\u00e9cisionnel conditionnel<br \/>\nif p_friedman &lt; 0.05:<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(\"nCondition valid\u00e9e : Pr\u00e9sence d'un effet global significatif.\")<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(\"Ex\u00e9cution du test post-hoc de Nemenyi en cours...\")<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;# Extraction de la matrice num\u00e9rique large<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;matrice_data = df_clinique[['TCC', 'Mindfulness', 'TIP']].to_numpy()<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;matrice_nemenyi = sp.posthoc_nemenyi_friedman(matrice_data)<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;# Structuration et labellisation du tableau<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;noms_therapies = ['TCC', 'Mindfulness', 'TIP']<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;matrice_nemenyi.index = noms_therapies<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;matrice_nemenyi.columns = noms_therapies<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(\"n=== MATRICE DES P-VALEURS AJUST\u00c9ES DE NEMENYI ===\")<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(matrice_nemenyi.round(5))<br \/>\nelse:<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;print(\"nTest de Friedman non significatif. Analyse interrompue conform\u00e9ment au protocole.\")<\/code><\/p>\n<h3>11.3 Discussion psychologique et interpr\u00e9tation des contrastes observ\u00e9s<\/h3>\n<p>L&rsquo;ex\u00e9cution de ce pipeline sur les donn\u00e9es cliniques g\u00e9n\u00e8re une statistique de Friedman de <em>Q = 30,00<\/em> (valeur maximale th\u00e9orique pour un devis \u00e0 15 sujets et 3 conditions parfaitement coh\u00e9rentes), associ\u00e9e \u00e0 une p-valeur de <em>p = 3,06 &times; 10<sup>-7<\/sup><\/em>. Le rejet de l&rsquo;hypoth\u00e8se nulle globale est donc sans ambigu\u00eft\u00e9.<\/p>\n<p>L&rsquo;inspection de la matrice post-hoc de Nemenyi d\u00e9livre les constats suivants :<\/p>\n<ul>\n<li><strong>TCC versus Mindfulness :<\/strong> La p-valeur ajust\u00e9e est de <em>p = 0,0010<\/em>. La TCC induit des scores d&rsquo;anxi\u00e9t\u00e9 r\u00e9siduelle significativement inf\u00e9rieurs \u00e0 ceux de la Pleine Conscience, traduisant une sup\u00e9riorit\u00e9 d&rsquo;action symptomatique rapide \u00e0 court terme sur ce profil clinique.<\/li>\n<li><strong>TCC versus TIP :<\/strong> La p-valeur s&rsquo;\u00e9tablit \u00e0 <em>p = 0,0010<\/em>. La TCC d\u00e9montre une sup\u00e9riorit\u00e9 statistique majeure et hautement significative par rapport \u00e0 la th\u00e9rapie interpersonnelle.<\/li>\n<li><strong>Mindfulness versus TIP :<\/strong> La p-valeur s&rsquo;\u00e9l\u00e8ve \u00e0 <em>p = 0,0010<\/em>. La Pleine Conscience d\u00e9montre \u00e9galement une sup\u00e9riorit\u00e9 statistiquement \u00e9tablie sur la th\u00e9rapie interpersonnelle dans la r\u00e9duction des scores HAM-A.<\/li>\n<\/ul>\n<p>Sur le plan psychologique, ces r\u00e9sultats mettent en lumi\u00e8re une gradation hi\u00e9rarchique nette dans l&rsquo;efficacit\u00e9 des protocoles th\u00e9rapeutiques br\u00e8ves sur l&rsquo;anxi\u00e9t\u00e9 g\u00e9n\u00e9ralis\u00e9e chez l&rsquo;\u00e9tudiant. La restructuration cognitive comportementale semble offrir les leviers les plus directement op\u00e9rants face aux distorsions cognitives li\u00e9es \u00e0 l&rsquo;\u00e9valuation universitaire, suivie par les protocoles de r\u00e9gulation \u00e9motionnelle par pleine conscience, tandis que les approches interpersonnelles apparaissent moins directement arm\u00e9es pour endiguer la symptomatologie somato-psychique aigu\u00eb \u00e0 court terme. La prudence clinique impose toutefois de confronter ces conclusions aux scores de maintien \u00e0 long terme lors de suivis longitudinaux ult\u00e9rieurs \u00e0 6 et 12 mois.<\/p>\n<h2>12. Bonnes pratiques m\u00e9thodologiques, limites et alternatives contemporaines<\/h2>\n<h3>12.1 Synth\u00e8se des r\u00e8gles d&rsquo;or pour la pratique de recherche<\/h3>\n<p>L&rsquo;application rigoureuse du test post-hoc de Nemenyi dans l&rsquo;\u00e9cosyst\u00e8me de la recherche comportementale repose sur le respect scrupuleux de plusieurs r\u00e8gles d&rsquo;or m\u00e9thodologiques :<\/p>\n<ul>\n<li><strong>Prise d&rsquo;acte imp\u00e9rative du test omnibus :<\/strong> Ne jamais engager d&rsquo;analyses post-hoc par paires si le test pr\u00e9alable de Friedman n&rsquo;a pas formellement rejet\u00e9 l&rsquo;hypoth\u00e8se nulle globale (<em>p &lt; 0,05<\/em>).<\/li>\n<li><strong>Double restitution descriptive :<\/strong> Les rangs moyens gouvernent la d\u00e9cision inf\u00e9rentielle, mais les m\u00e9dianes et \u00e9carts interquartiles d\u00e9crivent la r\u00e9alit\u00e9 des scores bruts. Toute publication rigoureuse doit pr\u00e9senter conjointement ces deux niveaux d&rsquo;information.<\/li>\n<li><strong>Transparence computationnelle et science ouverte :<\/strong> Conform\u00e9ment aux principes de la recherche reproductible, le script Python complet, les num\u00e9ros de version des biblioth\u00e8ques mobilis\u00e9es et les jeux de donn\u00e9es anonymis\u00e9s doivent \u00eatre mis \u00e0 disposition dans des d\u00e9p\u00f4ts p\u00e9rennes (tels que l&rsquo;Open Science Framework ou GitHub).<\/li>\n<\/ul>\n<h3>12.2 Limites intrins\u00e8ques du test de Nemenyi<\/h3>\n<p>En d\u00e9pit de son \u00e9l\u00e9gance math\u00e9matique, le test post-hoc de Nemenyi souffre de limitations structurelles incompressibles dont le chercheur doit \u00eatre pleinement conscient :<\/p>\n<p>En premier lieu, la proc\u00e9dure de Nemenyi est strictement restreinte aux devis unifactoriels intra-sujets. Elle s&rsquo;av\u00e8re structurellement incapable d&rsquo;\u00e9valuer des interactions factorielles complexes (par exemple, une interaction <em>Traitement &times; Temps &times; Sexe<\/em>). Si le design exp\u00e9rimental int\u00e8gre des facteurs inter-sujets crois\u00e9s avec des mesures r\u00e9p\u00e9t\u00e9es, l&rsquo;architecture non param\u00e9trique de Friedman-Nemenyi devient inapplicable.<\/p>\n<p>En second lieu, le test de Nemenyi exige une compl\u00e9tude absolue des donn\u00e9es. L&rsquo;absence d&rsquo;une seule mesure sur un seul point temporel conduit, par d\u00e9faut, \u00e0 l&rsquo;amputation de l&rsquo;int\u00e9gralit\u00e9 du sujet, entra\u00eenant une d\u00e9perdition dramatique d&rsquo;information lors d&rsquo;\u00e9tudes longitudinales caract\u00e9ris\u00e9es par des taux d&rsquo;attrition \u00e9lev\u00e9s.<\/p>\n<h3>12.3 Alternatives statistiques contemporaines<\/h3>\n<p>Face \u00e0 ces limitations, les sciences psychologiques contemporaines se tournent de plus en plus vers des architectures de mod\u00e9lisation plus flexibles et polyvalentes :<\/p>\n<p>Les Mod\u00e8les Lin\u00e9aires Mixtes G\u00e9n\u00e9ralis\u00e9s (Generalized Linear Mixed Models, GLMM) ou les mod\u00e8les \u00e0 cotes proportionnelles cumulatives (Cumulative Link Mixed Models, CLMM) repr\u00e9sentent l&rsquo;alternative moderne de r\u00e9f\u00e9rence pour les donn\u00e9es ordinales \u00e0 mesures r\u00e9p\u00e9t\u00e9es. Ces approches permettent de mod\u00e9liser simultan\u00e9ment des effets fixes multiples, leurs interactions d&rsquo;ordre sup\u00e9rieur, et des effets al\u00e9atoires propres \u00e0 chaque sujet, tout en tol\u00e9rant sans difficult\u00e9 les donn\u00e9es manquantes sous le postulat MAR (Missing At Random). Les biblioth\u00e8ques Python telles que <code>statsmodels<\/code> commencent \u00e0 int\u00e9grer ces fonctionnalit\u00e9s, bien que l&rsquo;environnement R (avec le package <code>ordinal<\/code>) demeure historiquement plus riche sur ce cr\u00e9neau sp\u00e9cifique.<\/p>\n<p>Enfin, les m\u00e9thodes bay\u00e9siennes non param\u00e9triques et les approches d&rsquo;estimation par r\u00e9\u00e9chantillonnage de rangs par bootstrap appari\u00e9 offrent une alternative s\u00e9duisante. Elles permettent d&rsquo;affiner l&rsquo;estimation des distributions a posteriori des diff\u00e9rences relatives de rangs sans d\u00e9pendre de l&rsquo;asymptotique parfois conservatrice de l&rsquo;\u00e9tendue studentis\u00e9e, ouvrant la voie \u00e0 une interpr\u00e9tation probabiliste directe particuli\u00e8rement f\u00e9conde pour l&rsquo;aide \u00e0 la d\u00e9cision en sant\u00e9 mentale.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Benjamini, Y., &amp; Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. <em>Journal of the Royal Statistical Society: Series B (Methodological)<\/em>, 57(1), 289\u2013300. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2517-6161.1995.tb02031.x\">https:\/\/doi.org\/10.1111\/j.2517-6161.1995.tb02031.x<\/a><\/li>\n<li>Conover, W. J., &amp; Iman, R. L. (1979). <em>On multiple-comparisons procedures<\/em> (Technical Report LA-7677-MS). Los Alamos Scientific Laboratory. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2172\/6057803\">https:\/\/doi.org\/10.2172\/6057803<\/a><\/li>\n<li>Dem\u0161ar, J. (2006). Statistical comparisons of classifiers over multiple data sets. <em>Journal of Machine Learning Research<\/em>, 7, 1\u201330. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.jmlr.org\/papers\/v7\/demsar06a.html\">https:\/\/www.jmlr.org\/papers\/v7\/demsar06a.html<\/a><\/li>\n<li>Friedman, M. (1937). The use of ranks to avoid the assumption of normality implicit in the analysis of variance. <em>Journal of the American Statistical Association<\/em>, 32(200), 675\u2013701. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1937.10503522\">https:\/\/doi.org\/10.1080\/01621459.1937.10503522<\/a><\/li>\n<li>Holm, S. (1979). A simple sequentially rejective multiple test procedure. <em>Scandinavian Journal of Statistics<\/em>, 6(2), 65\u201370. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.jstor.org\/stable\/4615733\">https:\/\/www.jstor.org\/stable\/4615733<\/a><\/li>\n<li>Nemenyi, P. B. (1963). <em>Distribution-free multiple comparisons<\/em> (Doctoral dissertation, Princeton University). University Microfilms.<\/li>\n<li>Terpilowski, M. (2019). scikit-posthocs: Pairwise multiple comparison post-hoc tests for Python. <em>Journal of Open Source Software<\/em>, 4(44), 1169. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.21105\/joss.01169\">https:\/\/doi.org\/10.21105\/joss.01169<\/a><\/li>\n<li>Tukey, J. W. (1949). Comparing individual means in the analysis of variance. <em>Biometrics<\/em>, 5(2), 99\u2013114. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/3001913\">https:\/\/doi.org\/10.2307\/3001913<\/a><\/li>\n<li>Vallat, R. (2018). Pingouin: statistics in Python. <em>Journal of Open Source Software<\/em>, 3(31), 1026. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.21105\/joss.01026\">https:\/\/doi.org\/10.21105\/joss.01026<\/a><\/li>\n<li>Virtanen, P., Gommers, R., Oliphant, T. E., Haberland, M., Reddy, T., Cournapeau, D., Burovski, E., Peterson, P., Weckesser, W., Bright, J., van der Walt, S. J., Brett, M., Wilson, J., Millman, K. J., Mayorov, N., Nelson, A. R. J., Jones, E., Kern, R., Larson, E., &#8230; SciPy 1.0 Contributors. (2020). SciPy 1.0: Fundamental algorithms for scientific computing in Python. <em>Nature Methods<\/em>, 17(3), 261\u2013272. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1038\/s41592-019-0686-2\">https:\/\/doi.org\/10.1038\/s41592-019-0686-2<\/a><\/li>\n<li>Wilcoxon, F. (1945). Individual comparisons by ranking methods. <em>Biometrics Bulletin<\/em>, 1(6), 80\u201383. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/3001968\">https:\/\/doi.org\/10.2307\/3001968<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour r\u00e9aliser et interpr\u00e9ter le test post-hoc de Nemenyi en Python apr\u00e8s un test de Friedman en recherche psychologique.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2464,1835],"tags":[2468,6132,2427,9200,2655],"class_list":["post-8002","psych_stat","type-psych_stat","status-publish","hentry","category-science-des-donnees","category-statistiques-en-psychologie","tag-python","tag-scikit-posthocs","tag-test-de-friedman","tag-test-de-nemenyi","tag-tests-post-hoc"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8002","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=8002"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8002\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=8002"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=8002"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=8002"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}