{"id":3547,"date":"2026-09-10T22:20:14","date_gmt":"2026-09-10T22:20:14","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-effectuer-test-grubbs-python\/"},"modified":"2026-09-10T22:20:14","modified_gmt":"2026-09-10T22:20:14","slug":"comment-effectuer-test-grubbs-python","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-effectuer-test-grubbs-python","title":{"rendered":"Comment effectuer le test de Grubbs en Python"},"content":{"rendered":"<p>L&rsquo;int\u00e9grit\u00e9 des conclusions scientifiques repose fondamentalement sur la rigueur du traitement appliqu\u00e9 aux observations empiriques. Dans le champ des sciences computationnelles, de la psychom\u00e9trie et de l&rsquo;analyse statistique contemporaine, la pr\u00e9sence d&rsquo;observations aberrantes \u2014 commun\u00e9ment d\u00e9sign\u00e9es sous le vocable d&rsquo;<em>outliers<\/em> \u2014 constitue une pr\u00e9occupation m\u00e9thodologique majeure. Ces points de donn\u00e9es atypiques, qu&rsquo;ils r\u00e9sultent d&rsquo;erreurs de mesure instrumentales, d&rsquo;artefacts d&rsquo;\u00e9chantillonnage ou de manifestations extr\u00eames mais authentiques du ph\u00e9nom\u00e8ne \u00e9tudi\u00e9, exercent une influence disproportionn\u00e9e sur les estimateurs classiques. Lorsque des m\u00e9triques param\u00e9triques fondamentales telles que la moyenne arithm\u00e9tique ou la variance d&rsquo;\u00e9chantillon sont mobilis\u00e9es, une unique valeur d\u00e9viante peut alt\u00e9rer la significativit\u00e9 statistique, inverser des corr\u00e9lations ou invalider les postulats des mod\u00e8les lin\u00e9aires g\u00e9n\u00e9raux.<\/p>\n<p>Face \u00e0 cet enjeu, le recours \u00e0 des crit\u00e8res d\u00e9cisionnels objectifs et math\u00e9matiquement formalis\u00e9s s&rsquo;av\u00e8re indispensable pour substituer une inf\u00e9rence probabiliste reproductible \u00e0 l&rsquo;intuition visuelle subjective des chercheurs. D\u00e9velopp\u00e9 au milieu du vingti\u00e8me si\u00e8cle, le test de Grubbs s&rsquo;est impos\u00e9 comme l&rsquo;une des proc\u00e9dures param\u00e9triques les plus \u00e9l\u00e9gantes et les plus robustes pour isoler une anomalie univari\u00e9e au sein d&rsquo;une distribution pr\u00e9sum\u00e9e gaussienne. En quantifiant l&rsquo;\u00e9cart studentis\u00e9 maximal entre une observation suspecte et le barycentre des donn\u00e9es, cette m\u00e9thode permet de statuer formellement sur l&rsquo;hypoth\u00e8se d&rsquo;une contamination de l&rsquo;\u00e9chantillon, tout en contr\u00f4lant rigoureusement le risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce.<\/p>\n<p>L&rsquo;\u00e9cosyst\u00e8me Python contemporain offre une infrastructure logicielle de premier plan pour ex\u00e9cuter, automatiser et personnaliser de telles proc\u00e9dures inf\u00e9rentielles. \u00c0 travers l&rsquo;exploitation conjointe de modules sp\u00e9cialis\u00e9s comme <em>outlier_utils<\/em> et des biblioth\u00e8ques scientifiques fondamentales que sont <em>NumPy<\/em>, <em>SciPy<\/em> et <em>Pandas<\/em>, le chercheur dispose d&rsquo;un outillage algorithmique complet pour scruter ses s\u00e9ries de donn\u00e9es. Cet article propose un guide exhaustif et approfondi d\u00e9di\u00e9 \u00e0 l&rsquo;ex\u00e9cution du test de Grubbs en Python, articulant les fondements th\u00e9oriques, les exigences math\u00e9matiques, la validation empirique des postulats pr\u00e9alables, l&rsquo;impl\u00e9mentation logicielle et les implications \u00e9pist\u00e9mologiques inh\u00e9rentes au traitement des donn\u00e9es d\u00e9viantes.<\/p>\n<h2>1. Introduction th\u00e9orique au test de Grubbs et \u00e0 la d\u00e9tection des valeurs aberrantes<\/h2>\n<h3>1.1 D\u00e9finition et origines du test de Grubbs<\/h3>\n<p>Le test de Grubbs, formalis\u00e9 en 1950 par le statisticien am\u00e9ricain Frank Ephraim Grubbs dans son article s\u00e9minal publi\u00e9 au sein des <a rel=\"noopener noreferrer\" href=\"https:\/\/www.jstor.org\/journal\/annmathstat\">Annals of Mathematical Statistics<\/a>, r\u00e9pond \u00e0 la n\u00e9cessit\u00e9 d&rsquo;\u00e9tablir une r\u00e8gle de rejet math\u00e9matiquement justifiable pour les observations hautement divergentes. Historiquement d\u00e9sign\u00e9 sous le terme de test de l&rsquo;\u00e9cart studentis\u00e9 maximal (\u00ab maximum normed residual test \u00bb), ce protocole a \u00e9t\u00e9 con\u00e7u pour d\u00e9terminer si la valeur la plus \u00e9loign\u00e9e de la moyenne dans un \u00e9chantillon de taille finie peut \u00eatre raisonnablement consid\u00e9r\u00e9e comme issue de la m\u00eame population gaussienne que les autres observations, ou si elle t\u00e9moigne d&rsquo;une perturbation externe.<\/p>\n<p>Sur le plan taxinomique, le test de Grubbs s&rsquo;inscrit au sein de la famille des proc\u00e9dures d&rsquo;\u00e9cart studentis\u00e9 maximal. Contrairement aux approches fond\u00e9es sur des seuils arbitraires non standardis\u00e9s, il op\u00e8re en divisant la d\u00e9viation maximale observ\u00e9e par l&rsquo;\u00e9cart-type de l&rsquo;\u00e9chantillon, int\u00e9grant ainsi la variabilit\u00e9 inh\u00e9rente aux donn\u00e9es dans la statistique de test. Cette normalisation interne conf\u00e8re au test une propri\u00e9t\u00e9 fondamentale : sa distribution d&rsquo;\u00e9chantillonnage ne d\u00e9pend pas des param\u00e8tres de position et d&rsquo;\u00e9chelle de la population parente, pourvu que l&rsquo;hypoth\u00e8se de normalit\u00e9 sous-jacente soit satisfaite.<\/p>\n<p>Il importe d&rsquo;\u00e9tablir une diff\u00e9renciation conceptuelle rigoureuse entre trois manifestations de variance atypique : le bruit de mesure al\u00e9atoire, la variabilit\u00e9 extr\u00eame biologique ou psychologique, et l&rsquo;erreur exp\u00e9rimentale syst\u00e9matique. Le bruit de mesure d\u00e9coule de l&rsquo;impr\u00e9cision in\u00e9vitable des capteurs ou des \u00e9chelles d&rsquo;\u00e9valuation, demeurant centr\u00e9 autour de la valeur vraie. La variabilit\u00e9 extr\u00eame traduit quant \u00e0 elle l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 naturelle d&rsquo;un syst\u00e8me complexe, documentant des \u00e9tats rares mais parfaitement valides. Enfin, l&rsquo;erreur exp\u00e9rimentale ou art\u00e9factuelle proc\u00e8de d&rsquo;un dysfonctionnement m\u00e9thodologique, d&rsquo;un d\u00e9faut de calibration mat\u00e9rielle ou d&rsquo;une erreur de saisie. Le test de Grubbs fournit une assise probabiliste permettant d&rsquo;isoler une observation dont la divergence d\u00e9passe ce que le hasard \u00e9chantillonnal peut expliquer, guidant l&rsquo;analyste dans l&rsquo;attribution de cette anomalie \u00e0 une population sous-jacente distincte.<\/p>\n<h3>1.2 Pertinence du test de Grubbs en recherche psychologique et comportementale<\/h3>\n<p>Dans les sciences du comportement, la psychom\u00e9trie et la psycholinguistique exp\u00e9rimentale, les donn\u00e9es quantitatives pr\u00e9sentent une vuln\u00e9rabilit\u00e9 chronique aux art\u00e9facts. L&rsquo;analyse des latences cognitives ou des temps de r\u00e9action chronom\u00e9triques illustre parfaitement cette vuln\u00e9rabilit\u00e9. Lors d&rsquo;une t\u00e2che d&rsquo;attention soutenue ou de d\u00e9cision lexicale, un sujet humain peut manifester un \u00e9pisode d&rsquo;inattention transitoire, d\u00e9clencher une r\u00e9ponse anticip\u00e9e involontaire en pressant une touche avant la perception du stimulus, ou subir une micro-distraction environnementale. De tels incidents g\u00e9n\u00e8rent des valeurs aberrantes qui ne refl\u00e8tent en rien les processus neurocognitifs cibl\u00e9s par le paradigme exp\u00e9rimental.<\/p>\n<p>L&rsquo;inclusion non r\u00e9gul\u00e9e de ces latences extr\u00eames produit une distorsion imm\u00e9diate de la moyenne arithm\u00e9tique et une inflation dramatique de la variance au sein des groupes exp\u00e9rimentaux. Dans le cadre de l&rsquo;application de mod\u00e8les lin\u00e9aires, tels que le test t de Student ou l&rsquo;analyse de variance (<a rel=\"noopener noreferrer\" href=\"https:\/\/www.sciencedirect.com\/topics\/mathematics\/analysis-of-variance\">ANOVA<\/a>), cette inflation r\u00e9siduelle affaiblit la puissance statistique globale du test en augmentant le d\u00e9nominateur des ratios de variance de Fisher. Par cons\u00e9quent, des effets exp\u00e9rimentaux authentiques risquent d&rsquo;\u00eatre masqu\u00e9s, conduisant \u00e0 une erreur de type II. Inversement, si une valeur extr\u00eame s&rsquo;agr\u00e8ge asym\u00e9triquement dans une condition exp\u00e9rimentale sp\u00e9cifique, elle peut induire artificiellement une diff\u00e9rence intergroupe statistiquement significative, provoquant une redoutable erreur de type I.<\/p>\n<p>Les r\u00e9percussions sur la validit\u00e9 scientifique globale sont consid\u00e9rables. La validit\u00e9 interne se trouve compromise d\u00e8s lors que la variance observ\u00e9e est imputable \u00e0 des art\u00e9facts plut\u00f4t qu&rsquo;\u00e0 la manipulation exp\u00e9rimentale de la variable ind\u00e9pendante. Parall\u00e8lement, la validit\u00e9 externe et la reproductibilit\u00e9 des conclusions souffrent d&rsquo;une telle d\u00e9gradation : une th\u00e9orie psychologique construite sur des estimations param\u00e9triques biais\u00e9es par des valeurs r\u00e9siduelles anormales \u00e9chouera syst\u00e9matiquement lors des tentatives de r\u00e9plication ind\u00e9pendantes men\u00e9es au sein d&rsquo;autres laboratoires universitaires.<\/p>\n<h3>1.3 Typologie des tests de d\u00e9tection : Grubbs vs autres proc\u00e9dures classiques<\/h3>\n<p>Le test de Grubbs coexiste avec une multiplicit\u00e9 de techniques d\u00e9di\u00e9es au traitement des valeurs atypiques, chacune reposant sur des postulats distincts. L&rsquo;une des heuristiques les plus r\u00e9pandues est la r\u00e8gle de l&rsquo;\u00e9cart interquartile conceptualis\u00e9e par John Tukey, qui d\u00e9finit des seuils d\u00e9cisionnels \u00e0 une distance de une fois et demie ou trois fois l&rsquo;intervalle interquartile au-del\u00e0 des premier et troisi\u00e8me quartiles. Si l&rsquo;approche de Tukey pr\u00e9sente l&rsquo;avantage d&rsquo;\u00eatre non param\u00e9trique et hautement r\u00e9sistante aux asym\u00e9tries distributionnelles, elle demeure une heuristique d\u00e9nu\u00e9e de niveau de signification probabiliste formel. Elle ne permet pas d&rsquo;\u00e9valuer le risque alpha exact associ\u00e9 \u00e0 la qualification d&rsquo;un point comme aberrant, contrairement au test de Grubbs qui calcule une probabilit\u00e9 formelle d&rsquo;erreur.<\/p>\n<p>Une autre alternative classique est constitu\u00e9e par les tests de ratio de Dixon, fr\u00e9quemment d\u00e9sign\u00e9s sous l&rsquo;appellation de tests Q de Dixon. Sp\u00e9cifiquement optimis\u00e9s pour les micro-\u00e9chantillons compris entre trois et trente observations, les tests de Dixon \u00e9valuent le rapport entre l&rsquo;\u00e9cart s\u00e9parant la valeur suspecte de sa voisine imm\u00e9diate et l&rsquo;\u00e9tendue totale de la s\u00e9rie. N\u00e9anmoins, lorsque la taille de l&rsquo;\u00e9chantillon cro\u00eet, l&rsquo;efficacit\u00e9 relative du test de Dixon d\u00e9cro\u00eet rapidement au profit du test de Grubbs, qui mobilise l&rsquo;ensemble de l&rsquo;information contenue dans la distribution par l&rsquo;entremise de la variance globale de l&rsquo;\u00e9chantillon.<\/p>\n<p>Enfin, un crit\u00e8re comparatif majeur r\u00e9side dans la vuln\u00e9rabilit\u00e9 au ph\u00e9nom\u00e8ne de masquage statistique (\u00ab masking effect \u00bb). Le masquage survient lorsque la pr\u00e9sence simultan\u00e9e de deux ou plusieurs valeurs aberrantes de magnitudes similaires dans le m\u00eame sens a pour effet de gonfler artificiellement la variance de l&rsquo;\u00e9chantillon. Cette inflation de l&rsquo;\u00e9cart-type r\u00e9duit m\u00e9caniquement la valeur du score standardis\u00e9 calcul\u00e9 par le test de Grubbs, rendant l&rsquo;anomalie ind\u00e9tectable. Alors que les m\u00e9thodes d&rsquo;estimation hautement robustes comme l&rsquo;\u00e9cart absolu \u00e0 la m\u00e9diane r\u00e9sistent intrins\u00e8quement au masquage, le test de Grubbs univari\u00e9 standard est originellement calibr\u00e9 pour la d\u00e9tection d&rsquo;une unique aberration, ce qui impose une vigilance analytique extr\u00eame lors de son ex\u00e9cution pratique.<\/p>\n<h2>2. Les hypoth\u00e8ses statistiques fondamentales et conditions d&rsquo;application<\/h2>\n<h3>2.1 Le postulat imp\u00e9ratif de normalit\u00e9 univari\u00e9e<\/h3>\n<p>L&rsquo;application l\u00e9gitime du test de Grubbs est subordonn\u00e9e au respect rigoureux de l&rsquo;hypoth\u00e8se de normalit\u00e9 univari\u00e9e de la distribution parente. Sur le plan math\u00e9matique, la distribution th\u00e9orique de la statistique de Grubbs d\u00e9rive directement des propri\u00e9t\u00e9s distributionnelles d&rsquo;un \u00e9chantillon al\u00e9atoire extrait d&rsquo;une loi normale d&rsquo;esp\u00e9rance math\u00e9matique mu et de variance sigma carr\u00e9. La d\u00e9rivation analytique de ses valeurs critiques repose sur l&rsquo;ind\u00e9pendance stochastique existant, sous l&rsquo;hypoth\u00e8se gaussienne, entre la moyenne empirique et la variance d&rsquo;\u00e9chantillon, conform\u00e9ment au c\u00e9l\u00e8bre th\u00e9or\u00e8me de Basu.<\/p>\n<p>D\u00e8s lors que la population d&rsquo;origine s&rsquo;\u00e9carte substantiellement de la sym\u00e9trie gaussienne, le cadre inf\u00e9rentiel du test s&rsquo;effondre. En pr\u00e9sence d&rsquo;une distribution asym\u00e9trique positive, ph\u00e9nom\u00e8ne extr\u00eamement fr\u00e9quent lors de la mesure de variables chronom\u00e9triques, les observations situ\u00e9es dans la longue queue de distribution g\u00e9n\u00e9reront de mani\u00e8re syst\u00e9matique des statistiques de Grubbs artificiellement gonfl\u00e9es. Le test conclura alors \u00e0 tort \u00e0 la pr\u00e9sence d&rsquo;une valeur aberrante avec un taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce largement sup\u00e9rieur au seuil nominal alpha fix\u00e9 par le chercheur. De m\u00eame, face \u00e0 des distributions leptokurtiques caract\u00e9ris\u00e9es par des queues lourdes, le test manifestera une instabilit\u00e9 critique.<\/p>\n<p>Il est donc m\u00e9thodologiquement proscrit d&rsquo;appliquer le test de Grubbs de mani\u00e8re m\u00e9canique sans avoir au pr\u00e9alable valid\u00e9 formellement l&rsquo;ad\u00e9quation de la s\u00e9rie de donn\u00e9es aux propri\u00e9t\u00e9s d&rsquo;une loi normale. Si la normalit\u00e9 est rejet\u00e9e, l&rsquo;analyste doit imp\u00e9rativement renoncer \u00e0 cette proc\u00e9dure param\u00e9trique au profit de transformations math\u00e9matiques stabilisatrices \u2014 telles que la transformation logarithmique ou la transformation de Box-Cox \u2014 ou s&rsquo;orienter vers des proc\u00e9dures non param\u00e9triques fond\u00e9es sur les rangs ou les estimateurs robustes.<\/p>\n<h3>2.2 Exigence relative \u00e0 la taille de l&rsquo;\u00e9chantillon<\/h3>\n<p>Outre la normalit\u00e9 distributionnelle, le test de Grubbs requiert un volume d&rsquo;observations minimal pour d\u00e9ployer sa validit\u00e9 math\u00e9matique. Les tables statistiques historiques et les travaux m\u00e9thodologiques contemporains du <a rel=\"noopener noreferrer\" href=\"https:\/\/www.nist.gov\/itl\/sed\/handbook\/eda\/section3\/eda35h.htm\">NIST\/SEMATECH e-Handbook of Statistical Methods<\/a> s&rsquo;accordent sur le crit\u00e8re formel imposant un \u00e9chantillon comportant au minimum sept observations distinctes (N sup\u00e9rieur ou \u00e9gal \u00e0 7). En de\u00e7\u00e0 de ce seuil critique, les degr\u00e9s de libert\u00e9 associ\u00e9s \u00e0 l&rsquo;estimation conjointe de la moyenne et de la variance d&rsquo;\u00e9chantillon deviennent trop t\u00e9nus pour distinguer math\u00e9matiquement une valeur v\u00e9ritablement d\u00e9viante d&rsquo;une simple fluctuation d&rsquo;\u00e9chantillonnage.<\/p>\n<p>Sur le plan asymptotique, lorsque la taille de l&rsquo;\u00e9chantillon s&rsquo;accro\u00eet consid\u00e9rablement pour atteindre plusieurs centaines ou milliers d&rsquo;individus, le comportement de la statistique de Grubbs converge vers les propri\u00e9t\u00e9s des valeurs extr\u00eames de Gumbel. Toutefois, dans les \u00e9chantillons massifs caract\u00e9ristiques du forage de donn\u00e9es contemporain, la puissance statistique du test devient tellement \u00e9crasante qu&rsquo;elle tend \u00e0 rejeter l&rsquo;hypoth\u00e8se nulle pour des d\u00e9viations marginales sans pertinence pratique ou substantielle. \u00c0 l&rsquo;inverse, au sein des micro-\u00e9chantillons fr\u00e9quemment rencontr\u00e9s en neurobiologie ou en psychologie exp\u00e9rimentale sp\u00e9cialis\u00e9e (N compris entre 3 et 6), le test souffre d&rsquo;un d\u00e9ficit d&rsquo;efficacit\u00e9 drastique, la borne th\u00e9orique sup\u00e9rieure de la statistique ne pouvant pas math\u00e9matiquement d\u00e9passer le seuil critique requis au niveau alpha usuel.<\/p>\n<p>La gestion des risques d&rsquo;erreurs inf\u00e9rentielles est intrins\u00e8quement li\u00e9e \u00e0 cette taille d&rsquo;\u00e9chantillon. Un \u00e9chantillon de dimension inad\u00e9quate expose l&rsquo;exp\u00e9rimentateur soit \u00e0 une inflation dramatique de l&rsquo;erreur de type I \u2014 rejet d&rsquo;observations parfaitement l\u00e9gitimes provoquant une r\u00e9duction artificielle de la variabilit\u00e9 r\u00e9elle du ph\u00e9nom\u00e8ne \u2014 soit \u00e0 une accumulation d&rsquo;erreurs de type II, maintenant dans le corpus analytique des donn\u00e9es hautement contamin\u00e9es qui d\u00e9t\u00e9riorent la pr\u00e9cision des estimations subs\u00e9quentes.<\/p>\n<h3>2.3 Ind\u00e9pendance et distribution identique des observations<\/h3>\n<p>L&rsquo;\u00e9difice probabiliste du test de Grubbs postule que l&rsquo;ensemble des observations vectoris\u00e9es constitue une suite de variables al\u00e9atoires ind\u00e9pendantes et identiquement distribu\u00e9es (i.i.d.). Ce postulat d&rsquo;ind\u00e9pendance signifie que l&rsquo;occurrence ou la magnitude d&rsquo;une mesure donn\u00e9e n&rsquo;exerce aucune influence causale, corr\u00e9lationnelle ou temporelle sur les autres observations du vecteur. Chaque point de donn\u00e9e doit r\u00e9sulter d&rsquo;une unit\u00e9 exp\u00e9rimentale autonome soumise aux m\u00eames conditions environnementales et instrumentales.<\/p>\n<p>Cette condition interdit l&rsquo;application directe du test de Grubbs standard aux structures de donn\u00e9es temporelles, aux signaux \u00e9lectrophysiologiques s\u00e9quentiels (tels que les enregistrements \u00e9lectroenc\u00e9phalographiques) ou aux devis de mesures r\u00e9p\u00e9t\u00e9es non contr\u00f4l\u00e9s. Dans une s\u00e9rie chronologique, la pr\u00e9sence d&rsquo;une autocorr\u00e9lation s\u00e9rielle viole le principe fondamental de la variance r\u00e9siduelle ind\u00e9pendante. Une valeur d\u00e9viante observ\u00e9e \u00e0 l&rsquo;instant t engendre une tra\u00een\u00e9e de d\u00e9pendance sur les observations t+1 et t+2, ce qui distord le calcul de l&rsquo;\u00e9cart-type global et biaise la m\u00e9trique de l&rsquo;\u00e9cart studentis\u00e9 maximal.<\/p>\n<p>Lorsque des mesures r\u00e9p\u00e9t\u00e9es proviennent du m\u00eame participant au sein d&rsquo;un protocole cognitif, les corr\u00e9lations intra-sujets violent le crit\u00e8re d&rsquo;\u00e9chantillonnage al\u00e9atoire simple. Dans ces configurations complexes, l&rsquo;analyste doit d&rsquo;abord extraire les r\u00e9sidus d&rsquo;un mod\u00e8le lin\u00e9aire mixte ou d&rsquo;une d\u00e9composition autor\u00e9gressive avant de pouvoir envisager de tester les d\u00e9viations standardis\u00e9es sur ces r\u00e9sidus d\u00e9corr\u00e9l\u00e9s, sous peine de formuler des conclusions inf\u00e9rentielles fallacieuses.<\/p>\n<h2>3. Formulation math\u00e9matique et cadre d\u00e9cisionnel du test de Grubbs<\/h2>\n<h3>3.1 \u00c9quations de la statistique G<\/h3>\n<p>La formalisation math\u00e9matique de la statistique de Grubbs repose sur l&rsquo;\u00e9valuation de la d\u00e9viation absolue maximale d&rsquo;une valeur par rapport \u00e0 la moyenne empirique de l&rsquo;\u00e9chantillon, mise \u00e0 l&rsquo;\u00e9chelle par l&rsquo;\u00e9cart-type non biais\u00e9. Soit un \u00e9chantillon univari\u00e9 constitu\u00e9 de n observations ind\u00e9pendantes ordonn\u00e9es ou non ordonn\u00e9es, not\u00e9 y indice 1, y indice 2, jusqu&rsquo;\u00e0 y indice n. La moyenne arithm\u00e9tique de cet \u00e9chantillon, not\u00e9e y barre, et son \u00e9cart-type d&rsquo;\u00e9chantillon corrig\u00e9, not\u00e9 s, sont formul\u00e9s selon les estimateurs non biais\u00e9s usuels.<\/p>\n<p>Pour un test bilat\u00e9ral ou non directionnel, visant \u00e0 identifier si l&rsquo;observation la plus distante du barycentre \u2014 qu&rsquo;elle soit exceptionnellement \u00e9lev\u00e9e ou exceptionnellement basse \u2014 constitue une valeur aberrante, la statistique G est d\u00e9finie par l&rsquo;\u00e9quation suivante :<\/p>\n<p><strong>G = max(|y_i &#8211; y_barre|) \/ s<\/strong>, pour i variant de 1 \u00e0 n.<\/p>\n<p>Cette expression capture le ratio de d\u00e9viation normalis\u00e9e le plus \u00e9lev\u00e9 pr\u00e9sent au sein du jeu de donn\u00e9es. Toutefois, lorsque l&rsquo;exp\u00e9rimentateur formule des hypoth\u00e8ses scientifiques directionnelles a priori, la statistique de Grubbs peut \u00eatre d\u00e9clin\u00e9e de mani\u00e8re unilat\u00e9rale. Pour tester si la valeur minimale observ\u00e9e est anormalement basse, la statistique unilat\u00e9rale minimale s&rsquo;exprime par le ratio :<\/p>\n<p><strong>G_min = (y_barre &#8211; y_min) \/ s<\/strong><\/p>\n<p>Inversement, pour \u00e9valuer si la valeur maximale d&rsquo;un vecteur est anormalement \u00e9lev\u00e9e par rapport \u00e0 la dispersion collective des donn\u00e9es, la statistique unilat\u00e9rale maximale est formul\u00e9e par :<\/p>\n<p><strong>G_max = (y_max &#8211; y_barre) \/ s<\/strong><\/p>\n<p>Dans les deux configurations unilat\u00e9rales, le num\u00e9rateur pr\u00e9serve le signe alg\u00e9brique de la d\u00e9viation d&rsquo;int\u00e9r\u00eat afin de sonder sp\u00e9cifiquement l&rsquo;extr\u00e9mit\u00e9 sup\u00e9rieure ou inf\u00e9rieure de la distribution empirique.<\/p>\n<h3>3.2 Hypoth\u00e8ses nulle et alternative<\/h3>\n<p>Le cadre inf\u00e9rentiel n\u00e9o-fissh\u00e9rien ou de Neyman-Pearson exige la sp\u00e9cification formelle d&rsquo;un couple d&rsquo;hypoth\u00e8ses mutuellement exclusives avant le calcul de la statistique de test. L&rsquo;hypoth\u00e8se nulle, d\u00e9sign\u00e9e conventionnellement par H0, postule l&rsquo;homog\u00e9n\u00e9it\u00e9 distributionnelle int\u00e9grale de l&rsquo;\u00e9chantillon. Formellement, sous H0, l&rsquo;ensemble des n observations observ\u00e9es sont issues d&rsquo;une population parente unique distribu\u00e9e normalement, avec une moyenne mu et une variance sigma carr\u00e9 invariantes :<\/p>\n<p><strong>H0 : Toutes les observations y_i proviennent de la m\u00eame distribution normale N(mu, sigma^2). Il n&rsquo;existe aucune valeur aberrante au sein de l&rsquo;\u00e9chantillon.<\/strong><\/p>\n<p>L&rsquo;hypoth\u00e8se alternative, d\u00e9sign\u00e9e par H1, postule une contamination de la distribution par un m\u00e9canisme g\u00e9n\u00e9rateur perturbateur. La nature formelle de H1 diverge selon la directionnalit\u00e9 s\u00e9lectionn\u00e9e :<\/p>\n<ul>\n<li>Dans le cas d&rsquo;un test bilat\u00e9ral non directionnel, H1 \u00e9nonce qu&rsquo;au moins une observation (celle qui maximise l&rsquo;\u00e9cart studentis\u00e9 absolu) ne provient pas de la distribution parente normale r\u00e9gissant le reste de l&rsquo;\u00e9chantillon, pr\u00e9sentant un glissement de moyenne vers une valeur diff\u00e9rente.<\/li>\n<li>Dans le cas d&rsquo;un test unilat\u00e9ral maximal, H1 sp\u00e9cifie formellement que l&rsquo;observation maximale y_max provient d&rsquo;une distribution distincte caract\u00e9ris\u00e9e par une esp\u00e9rance math\u00e9matique strictement sup\u00e9rieure \u00e0 mu.<\/li>\n<li>Dans le cas d&rsquo;un test unilat\u00e9ral minimal, H1 sp\u00e9cifie que l&rsquo;observation minimale y_min est issue d&rsquo;une distribution distincte affichant une esp\u00e9rance math\u00e9matique strictement inf\u00e9rieure \u00e0 mu.<\/li>\n<\/ul>\n<h3>3.3 Calcul de la valeur critique bas\u00e9e sur la distribution de Student<\/h3>\n<p>La distribution d&rsquo;\u00e9chantillonnage exacte de la statistique de Grubbs sous l&rsquo;hypoth\u00e8se nulle n&rsquo;est pas gaussienne ; elle est intimement li\u00e9e \u00e0 la loi de Student par l&rsquo;interm\u00e9diaire d&rsquo;une transformation non lin\u00e9aire qui tient compte de la corr\u00e9lation inh\u00e9rente entre les \u00e9carts individuels et l&rsquo;\u00e9cart-type global de l&rsquo;\u00e9chantillon. Frank Grubbs a d\u00e9montr\u00e9 que la valeur critique th\u00e9orique, d\u00e9sign\u00e9e par G_crit, pour un seuil de signification alpha donn\u00e9 et une taille d&rsquo;\u00e9chantillon n, est calcul\u00e9e exactement par l&rsquo;\u00e9quation math\u00e9matique suivante :<\/p>\n<p><strong>G_crit = ((n &#8211; 1) \/ sqrt(n)) * sqrt(t^2 \/ (n &#8211; 2 + t^2))<\/strong><\/p>\n<p>Dans cette formulation analytique, le terme t repr\u00e9sente la valeur critique sup\u00e9rieure issue d&rsquo;une loi de Student \u00e0 (n &#8211; 2) degr\u00e9s de libert\u00e9. L&rsquo;\u00e9l\u00e9ment d\u00e9terminant r\u00e9side dans l&rsquo;ajustement du seuil de probabilit\u00e9 utilis\u00e9 pour d\u00e9terminer le quantile de Student. Pour un test de Grubbs bilat\u00e9ral, le niveau de signification unitaire appliqu\u00e9 au calcul de la valeur critique de Student est ajust\u00e9 par une correction de type <a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/BonferroniCorrection.html\">Bonferroni<\/a> pour neutraliser le risque d&rsquo;inflation alpha r\u00e9sultant de la recherche du maximum parmi n observations :<\/p>\n<p><strong>t = t_critique(alpha \/ (2n), n &#8211; 2)<\/strong><\/p>\n<p>Pour un test de Grubbs unilat\u00e9ral, l&rsquo;ajustement du quantile se formule en retenant un niveau d&rsquo;erreur \u00e9gal \u00e0 alpha divis\u00e9 par n :<\/p>\n<p><strong>t = t_critique(alpha \/ n, n &#8211; 2)<\/strong><\/p>\n<p>La r\u00e8gle de d\u00e9cision statistique s&rsquo;\u00e9nonce de mani\u00e8re univoque : si la statistique G calcul\u00e9e \u00e0 partir des donn\u00e9es empiriques est strictement sup\u00e9rieure \u00e0 la valeur critique th\u00e9orique G_crit calcul\u00e9e pour le niveau alpha choisi, l&rsquo;hypoth\u00e8se nulle H0 est rejet\u00e9e avec un niveau de confiance de 1 moins alpha. L&rsquo;observation associ\u00e9e est alors d\u00e9clar\u00e9e formellement aberrante au sens de Grubbs.<\/p>\n<h2>4. Configuration de l&rsquo;environnement Python et installation des biblioth\u00e8ques<\/h2>\n<h3>4.1 Architecture logicielle et gestion des d\u00e9pendances<\/h3>\n<p>L&rsquo;impl\u00e9mentation rigoureuse d&rsquo;analyses statistiques en Python requiert une isolation m\u00e9ticuleuse de l&rsquo;environnement de travail afin de garantir la p\u00e9rennit\u00e9 et la reproductibilit\u00e9 des calculs. L&rsquo;utilisation d&rsquo;environnements virtuels standardis\u00e9s constitue une exigence fondamentale de l&rsquo;ing\u00e9nierie logicielle scientifique. Il est vivement recommand\u00e9 d&rsquo;initialiser un environnement d\u00e9di\u00e9 par l&rsquo;entremise de l&rsquo;utilitaire venv natif de Python ou via le gestionnaire d&rsquo;environnements Conda. Cette pratique pr\u00e9vient tout conflit de versions entre les d\u00e9pendances sous-jacentes des biblioth\u00e8ques num\u00e9riques compil\u00e9es en langage C ou Fortran.<\/p>\n<p>La conduite du test de Grubbs repose principalement sur le paquet sp\u00e9cialis\u00e9 <em>outlier_utils<\/em>, accessible publiquement via l&rsquo;index de paquets officiel PyPI. Son installation s&rsquo;effectue ais\u00e9ment \u00e0 l&rsquo;aide de l&rsquo;outil pip au sein du terminal d&rsquo;instructions. Parall\u00e8lement, l&rsquo;environnement doit imp\u00e9rativement h\u00e9berger les trois piliers du calcul scientifique en Python : <em>NumPy<\/em> pour l&rsquo;alg\u00e8bre vectorielle et la manipulation d&rsquo;\u00e9chantillons multidimensionnels, <em>SciPy<\/em> pour l&rsquo;\u00e9valuation des lois de probabilit\u00e9s et des quantiles statistiques th\u00e9oriques, et <em>Pandas<\/em> pour la structuration tubulaire et le nettoyage des donn\u00e9es exp\u00e9rimentales.<\/p>\n<p>L&rsquo;int\u00e9grit\u00e9 de la cha\u00eene analytique doit \u00eatre valid\u00e9e d\u00e8s l&rsquo;initialisation de l&rsquo;interpr\u00e9teur Python. Il convient de s&rsquo;assurer que les biblioth\u00e8ques import\u00e9es affichent des versions compatibles avec les sp\u00e9cifications de calcul matriciel modernes. La synchronisation pr\u00e9cise entre SciPy et NumPy est particuli\u00e8rement critique, car le calcul des degr\u00e9s de libert\u00e9 et l&rsquo;\u00e9valuation de la distribution cumulative de Student reposent sur les sous-routines de physique statistique de la biblioth\u00e8que Cephes int\u00e9gr\u00e9e au c\u0153ur de SciPy.<\/p>\n<h3>4.2 Structure du paquet outlier_utils et du module smirnov_grubbs<\/h3>\n<p>Le paquet <em>outlier_utils<\/em> structure ses fonctionnalit\u00e9s de d\u00e9tection au sein d&rsquo;un sous-module embl\u00e9matique nomm\u00e9 <em>smirnov_grubbs<\/em>. Ce nommage rend hommage \u00e0 la contribution math\u00e9matique parall\u00e8le de Nikola\u00ef Smirnov, qui a \u00e9galement document\u00e9 les propri\u00e9t\u00e9s de ces \u00e9carts standardis\u00e9s extr\u00eames dans la litt\u00e9rature sovi\u00e9tique des ann\u00e9es 1940. Ce module encapsule une s\u00e9rie de fonctions fonctionnelles optimis\u00e9es, permettant d&rsquo;\u00e9viter une programmation manuelle fastidieuse tout en minimisant les risques de bugs num\u00e9riques.<\/p>\n<p>La signature standard de la fonction ma\u00eetresse s&rsquo;articule autour de la m\u00e9thode test du module smirnov_grubbs. Cette fonction accepte deux arguments directeurs : le vecteur de donn\u00e9es empiriques sous forme d&rsquo;it\u00e9rable num\u00e9rique, et le param\u00e8tre scalaire alpha d\u00e9finissant le seuil de signification de l&rsquo;erreur de premi\u00e8re esp\u00e8ce, calibr\u00e9 par d\u00e9faut \u00e0 la valeur canonique de 0.05. En interne, l&rsquo;algorithme proc\u00e8de \u00e0 l&rsquo;\u00e9valuation de la distribution de Grubbs, compare la valeur extr\u00eame observ\u00e9e au seuil critique th\u00e9orique, puis retourne une nouvelle structure de donn\u00e9es purg\u00e9e des \u00e9l\u00e9ments d\u00e9clar\u00e9s aberrants.<\/p>\n<p>Un atout architectural majeur de cette impl\u00e9mentation r\u00e9side dans son agilit\u00e9 face aux formats d&rsquo;entr\u00e9e h\u00e9t\u00e9rog\u00e8nes. La m\u00e9thode accepte de mani\u00e8re transparente des listes standards Python, des tableaux vectoriels unidimensionnels cr\u00e9\u00e9s via NumPy ou des objets de type Series issus de l&rsquo;\u00e9cosyst\u00e8me Pandas. En coulisses, l&rsquo;algorithme convertit l&rsquo;entr\u00e9e en un tableau NumPy contigu en m\u00e9moire vive, optimisant ainsi les calculs de somme des carr\u00e9s et d&rsquo;\u00e9carts centrom\u00e9triques sans surcharge de calcul pour l&rsquo;analyste.<\/p>\n<h2>5. Pr\u00e9traitement et structuration des donn\u00e9es exp\u00e9rimentales<\/h2>\n<h3>5.1 Importation de jeux de donn\u00e9es psychom\u00e9triques<\/h3>\n<p>Le point de d\u00e9part de toute investigation r\u00e9side dans l&rsquo;acquisition et le chargement m\u00e9thodique des matrices de donn\u00e9es collect\u00e9es lors d&rsquo;exp\u00e9rimentations comportementales. \u00c0 l&rsquo;aide de la biblioth\u00e8que Pandas, les fichiers de formats usuels, tels que les valeurs s\u00e9par\u00e9es par des virgules (CSV) ou les feuilles de calcul Microsoft Excel (XLSX), sont instanci\u00e9s sous la forme d&rsquo;un objet DataFrame. Cet objet structure les variables sous forme de colonnes typ\u00e9es, facilitant grandement les op\u00e9rations subs\u00e9quentes de d\u00e9coupage et de filtrage statistique.<\/p>\n<p>L&rsquo;extraction de la variable continue d&rsquo;int\u00e9r\u00eat, comme les latences cognitives exprim\u00e9es en millisecondes lors de t\u00e2ches d&rsquo;\u00e9vocation s\u00e9mantique, n\u00e9cessite un ciblage vectoriel strict. \u00c0 ce stade pr\u00e9liminaire, une vigilance toute particuli\u00e8re doit \u00eatre accord\u00e9e au traitement des donn\u00e9es manquantes, encod\u00e9es conventionnellement sous la valeur NaN (Not a Number) ou None. L&rsquo;existence de cellules vides au sein d&rsquo;un vecteur de calcul math\u00e9matique perturbe irr\u00e9m\u00e9diablement l&rsquo;\u00e9valuation de la moyenne arithm\u00e9tique et de l&rsquo;\u00e9cart-type, aboutissant \u00e0 la propagation de r\u00e9sultats ind\u00e9termin\u00e9s.<\/p>\n<p>L&rsquo;analyste doit mettre en \u0153uvre une politique de gestion des valeurs manquantes explicite et justifi\u00e9e m\u00e9thodologiquement. Dans le cadre univari\u00e9 propre au test de Grubbs, la suppression par liste des observations non renseign\u00e9es constitue l&rsquo;approche conventionnelle la plus robuste, sous r\u00e9serve que ces omissions r\u00e9pondent au m\u00e9canisme de donn\u00e9es manquantes de mani\u00e8re totalement al\u00e9atoire (Missing Completely at Random, MCAR). Les commandes de filtrage de Pandas permettent d&rsquo;\u00e9purer instantan\u00e9ment les s\u00e9ries des valeurs nulles avant toute conversion vers un vecteur num\u00e9rique pur destin\u00e9 au calcul inf\u00e9rentiel.<\/p>\n<h3>5.2 Exploration descriptive univari\u00e9e<\/h3>\n<p>Avant d&rsquo;initier toute proc\u00e9dure inf\u00e9rentielle formalis\u00e9e, une analyse exploratoire des donn\u00e9es univari\u00e9es est indispensable. Cette phase empirique fondamentale permet au chercheur de s&rsquo;impr\u00e9gner de l&rsquo;architecture distributionnelle de ses mesures. Elle d\u00e9bute par l&rsquo;extraction des estimateurs de tendance centrale classiques (moyenne arithm\u00e9tique, m\u00e9diane d&rsquo;\u00e9chantillon) et des m\u00e9triques de dispersion canoniques (variance d&rsquo;\u00e9chantillon non biais\u00e9e, \u00e9cart-type, \u00e9tendue interquartile). La comparaison imm\u00e9diate entre la moyenne et la m\u00e9diane offre un premier indicateur pr\u00e9cieux de l&rsquo;asym\u00e9trie potentielle de la s\u00e9rie.<\/p>\n<p>Cette description num\u00e9rique doit imp\u00e9rativement \u00eatre consolid\u00e9e par des repr\u00e9sentations visuelles graphiques adapt\u00e9es. La g\u00e9n\u00e9ration d&rsquo;histogrammes combin\u00e9s \u00e0 une estimation de densit\u00e9 par noyau (Kernel Density Estimation, KDE) permet d&rsquo;observer la topologie de la distribution empirique et d&rsquo;identifier la pr\u00e9sence \u00e9ventuelle de bimodalit\u00e9 ou d&rsquo;\u00e9talements anormaux des queues de distribution. L&rsquo;estimation par noyau offre une approximation continue de la fonction de densit\u00e9 sous-jacente, lissant les irr\u00e9gularit\u00e9s artificielles cr\u00e9\u00e9es par le choix arbitraire de la largeur des classes de l&rsquo;histogramme.<\/p>\n<p>Parall\u00e8lement, la conception de diagrammes en bo\u00eete, ou boxplots, constitue un instrument privil\u00e9gi\u00e9 pour rep\u00e9rer intuitivement les singularit\u00e9s num\u00e9riques. En visualisant graphiquement les quartiles, la m\u00e9diane et les moustaches d\u00e9limit\u00e9es par les crit\u00e8res de Tukey, le chercheur localise imm\u00e9diatement les candidats potentiels \u00e0 l&rsquo;aberrance. Cette exploration visuelle ne saurait en aucun cas se substituer au test formel de Grubbs, mais elle remplit une mission d&rsquo;alerte diagnostique cruciale, permettant de d\u00e9tecter d&rsquo;\u00e9ventuelles anomalies globales avant de soumettre les donn\u00e9es aux \u00e9quations d\u00e9cisionnelles.<\/p>\n<h2>6. Validation empirique des postulats pr\u00e9alables en Python<\/h2>\n<h3>6.1 \u00c9valuation de la normalit\u00e9 avec SciPy<\/h3>\n<p>La validit\u00e9 des conclusions issues du test de Grubbs \u00e9tant subordonn\u00e9e \u00e0 la normalit\u00e9 des donn\u00e9es, une validation empirique pr\u00e9alable s&rsquo;impose de mani\u00e8re cat\u00e9gorique. Le sous-module <em>scipy.stats<\/em> int\u00e8gre un \u00e9ventail complet d&rsquo;outils probabilistes permettant de tester cette ad\u00e9quation. Parmi eux, le test de Shapiro-Wilk, accessible via la fonction <em>scipy.stats.shapiro<\/em>, fait figure de r\u00e9f\u00e9rence m\u00e9thodologique pour les \u00e9chantillons de taille petite \u00e0 mod\u00e9r\u00e9e (typiquement pour N compris entre 10 et 50 observations). Ce test compare la corr\u00e9lation empirique des quantiles de l&rsquo;\u00e9chantillon avec ceux d&rsquo;une loi normale th\u00e9orique ; une p-valeur inf\u00e9rieure au seuil critique alpha (g\u00e9n\u00e9ralement 0.05) commande le rejet cat\u00e9gorique de l&rsquo;hypoth\u00e8se de normalit\u00e9.<\/p>\n<p>Pour des jeux de donn\u00e9es plus cons\u00e9quents, le test omnibus de D&rsquo;Agostino-Pearson, mobilisable par la fonction <em>scipy.stats.normaltest<\/em>, pr\u00e9sente une robustesse th\u00e9orique remarquable. Cette proc\u00e9dure statistique op\u00e8re en \u00e9valuant simultan\u00e9ment le coefficient d&rsquo;asym\u00e9trie (skewness) et le coefficient d&rsquo;aplatissement (kurtosis) de la distribution empirique. En combinant les transformations normalisantes de ces deux estimateurs dans un test du chi-deux \u00e0 deux degr\u00e9s de libert\u00e9, le test de D&rsquo;Agostino-Pearson d\u00e9tecte avec une sensibilit\u00e9 \u00e9lev\u00e9e les d\u00e9viations structurelles par rapport \u00e0 la courbe en cloche gaussienne.<\/p>\n<p>Cette v\u00e9rification num\u00e9rique doit \u00eatre compl\u00e9t\u00e9e par l&rsquo;inspection visuelle d&rsquo;un diagramme Quantile-Quantile (Q-Q plot), g\u00e9n\u00e9r\u00e9 conventionnellement via la fonction <em>scipy.stats.probplot<\/em>. Ce graphique projette les quantiles observ\u00e9s de l&rsquo;\u00e9chantillon empirique contre les quantiles th\u00e9oriques attendus d&rsquo;une loi normale standardis\u00e9e id\u00e9ale. Si les donn\u00e9es d\u00e9rivent d&rsquo;une distribution gaussienne, les points s&rsquo;alignent rigoureusement le long de la bissectrice th\u00e9orique. Des d\u00e9crochages marqu\u00e9s aux extr\u00e9mit\u00e9s de la ligne traduisent sans \u00e9quivoque des queues de distribution \u00e9paisses ou une asym\u00e9trie r\u00e9dhibitoire pour l&rsquo;application l\u00e9gitime du test de Grubbs.<\/p>\n<h3>6.2 V\u00e9rification algorithmique de la taille d&rsquo;\u00e9chantillon<\/h3>\n<p>Dans une perspective de programmation d\u00e9fensive et de rigueur m\u00e9thodologique, il est imp\u00e9ratif d&rsquo;int\u00e9grer au sein des scripts Python une barri\u00e8re de validation algorithmique contr\u00f4lant formellement la dimension du vecteur de donn\u00e9es avant toute sollicitation des fonctions de Grubbs. Comme \u00e9tabli lors de l&rsquo;examen des fondements th\u00e9oriques, un effectif d&rsquo;au moins sept observations (N sup\u00e9rieur ou \u00e9gal \u00e0 7) est requis pour assurer la consistance des calculs d&rsquo;\u00e9carts studentis\u00e9s.<\/p>\n<p>Cette validation s&rsquo;impl\u00e9mente par le biais d&rsquo;instructions conditionnelles ou d&rsquo;assertions structurelles levant des exceptions explicites lorsque le vecteur analys\u00e9 ne satisfait pas au crit\u00e8re dimensionnel requis. Si la fonction len() appliqu\u00e9e au tableau de donn\u00e9es nettoy\u00e9 renvoie une valeur strictement inf\u00e9rieure \u00e0 sept, le programme doit imm\u00e9diatement suspendre l&rsquo;ex\u00e9cution du test param\u00e9trique et \u00e9mettre un message d&rsquo;erreur ou d&rsquo;avertissement document\u00e9, pr\u00e9venant ainsi la production de r\u00e9sultats statistiques d\u00e9nu\u00e9s de toute signification probabiliste.<\/p>\n<p>Face \u00e0 un \u00e9chantillon ne satisfaisant pas aux crit\u00e8res de taille minimale ou en situation d&rsquo;\u00e9chec av\u00e9r\u00e9 des tests de normalit\u00e9 de Shapiro-Wilk, le script d&rsquo;analyse doit basculer automatiquement vers des strat\u00e9gies alternatives non param\u00e9triques. L&rsquo;analyste se tournera avantageusement vers le test de Dixon pour les micro-\u00e9chantillons normaux, ou adoptera des approches fond\u00e9es sur l&rsquo;\u00e9cart absolu \u00e0 la m\u00e9diane (MAD) ou la r\u00e8gle de Tukey, \u00e9vitant ainsi d&rsquo;imposer un mod\u00e8le math\u00e9matique inadapt\u00e9 \u00e0 une r\u00e9alit\u00e9 empirique non conforme.<\/p>\n<h2>7. Ex\u00e9cution pratique du test de Grubbs bilat\u00e9ral avec outlier_utils<\/h2>\n<h3>7.1 Impl\u00e9mentation pas \u00e0 pas de smirnov_grubbs.test()<\/h3>\n<p>L&rsquo;ex\u00e9cution standard d&rsquo;un test de Grubbs bilat\u00e9ral via la biblioth\u00e8que <em>outlier_utils<\/em> s&rsquo;orchestre avec une remarquable concision syntaxique, tout en mobilisant en sous-main des calculs matriciels \u00e9labor\u00e9s. La premi\u00e8re \u00e9tape op\u00e9ratoire consiste \u00e0 convertir la s\u00e9rie de mesures exp\u00e9rimentales \u2014 pr\u00e9alablement valid\u00e9e quant \u00e0 sa normalit\u00e9 et son volume \u2014 en un vecteur de calcul standardis\u00e9, g\u00e9n\u00e9ralement sous la forme d&rsquo;un tableau unidimensionnel NumPy.<\/p>\n<p>La d\u00e9tection de l&rsquo;aberration est alors d\u00e9clench\u00e9e en sollicitant directement la fonction test h\u00e9berg\u00e9e au sein du sous-module <em>smirnov_grubbs<\/em>. L&rsquo;analyste transmet \u00e0 cette fonction deux param\u00e8tres fondamentaux : la collection num\u00e9rique cibl\u00e9e et le param\u00e8tre alpha gouvernant le niveau de rigueur inf\u00e9rentielle, usuellement fix\u00e9 \u00e0 0.05 pour tol\u00e9rer une probabilit\u00e9 de cinq pour cent de fausse alarme statistique. Lors de cet appel synchrone, le module ex\u00e9cute les calculs analytiques de la moyenne, de l&rsquo;\u00e9cart-type non biais\u00e9, localise le point maximisant la d\u00e9viation absolue, \u00e9value la statistique G contre la valeur critique issue de la distribution de Student, et formule le verdict inf\u00e9rentiel.<\/p>\n<p>L&rsquo;objet retourn\u00e9 en sortie par l&rsquo;ex\u00e9cution de la m\u00e9thode <em>smirnov_grubbs.test()<\/em> est une nouvelle s\u00e9quence de donn\u00e9es, typ\u00e9e \u00e0 l&rsquo;identique de la structure inject\u00e9e en entr\u00e9e. Si aucune observation n&rsquo;exc\u00e8de la valeur critique th\u00e9orique au seuil alpha sp\u00e9cifi\u00e9, la collection retourn\u00e9e est rigoureusement identique au jeu de donn\u00e9es initial. En revanche, si une observation pr\u00e9sente un \u00e9cart studentis\u00e9 statistiquement significatif, elle est automatiquement extirp\u00e9e du vecteur, offrant ainsi au chercheur le sous-ensemble imm\u00e9diatement purg\u00e9 pour ses analyses ult\u00e9rieures.<\/p>\n<h3>7.2 Extraction programmatique des valeurs et indices aberrants<\/h3>\n<p>Si l&rsquo;obtention directe du vecteur assaini s&rsquo;av\u00e8re commode pour les cha\u00eenes de traitement automatis\u00e9es, la d\u00e9marche scientifique exige une tra\u00e7abilit\u00e9 documentaire int\u00e9grale des \u00e9l\u00e9ments rejet\u00e9s. Le chercheur doit pouvoir documenter pr\u00e9cis\u00e9ment quelle observation a \u00e9t\u00e9 qualifi\u00e9e d&rsquo;aberrante, quelle \u00e9tait sa magnitude m\u00e9trique exacte, et \u00e0 quel sujet ou essai exp\u00e9rimental elle correspondait au sein du registre d&rsquo;origine. La m\u00e9thode <em>test()<\/em> \u00e9l\u00e9mentaire s&rsquo;av\u00e8re insuffisante pour cette t\u00e2che puisqu&rsquo;elle occulte l&rsquo;identit\u00e9 contextuelle du point \u00e9lagu\u00e9.<\/p>\n<p>Pour surmonter cette limitation, le module fournit une m\u00e9thode jumelle indispensable : <em>two_sided_test_indices()<\/em>. Au lieu de renvoyer le sous-ensemble \u00e9pur\u00e9 des valeurs num\u00e9riques, cette fonction interroge le vecteur d&rsquo;origine et renvoie sous forme de liste Python l&rsquo;ensemble des coordonn\u00e9es scalaires, c&rsquo;est-\u00e0-dire les index positionnels, correspondant aux observations ayant provoqu\u00e9 le rejet de l&rsquo;hypoth\u00e8se nulle. Cette disposition algorithmique est d&rsquo;une grande \u00e9l\u00e9gance, car elle permet de relier instantan\u00e9ment l&rsquo;anomalie aux m\u00e9tadonn\u00e9es associ\u00e9es dans un DataFrame Pandas parent.<\/p>\n<p>Gr\u00e2ce \u00e0 ces coordonn\u00e9es positionnelles extraites par <em>two_sided_test_indices()<\/em>, il devient trivial d&rsquo;interroger la table de donn\u00e9es initiale par indexation s\u00e9lective pour isoler la ligne d&rsquo;enregistrement d\u00e9faillante. Le statisticien peut d\u00e8s lors archiver la valeur exacte de l&rsquo;anomalie, analyser l&rsquo;identifiant du participant concern\u00e9, l&rsquo;heure de passation de l&rsquo;\u00e9preuve et les param\u00e8tres instrumentaux annexes, assurant une conformit\u00e9 exemplaire avec les exigences de tra\u00e7abilit\u00e9 m\u00e9thodologique.<\/p>\n<h3>7.3 Analyse du comportement face \u00e0 diff\u00e9rents niveaux de confiance<\/h3>\n<p>Le choix du niveau de signification alpha constitue un arbitrage statistique d&rsquo;une port\u00e9e d\u00e9terminante sur l&rsquo;agressivit\u00e9 du test de Grubbs. La modulation intentionnelle de ce seuil permet de mesurer la sensibilit\u00e9 du mod\u00e8le face au conservatisme m\u00e9thodologique. Traditionnellement calibr\u00e9 \u00e0 0.05 dans les sciences humaines et comportementales, ce param\u00e8tre peut \u00eatre restreint \u00e0 0.01 pour exiger un niveau de preuve sup\u00e9rieur, ou assoupli \u00e0 0.10 dans le cadre d&rsquo;\u00e9tudes exploratoires pilotes cherchant \u00e0 maximiser le pouvoir de d\u00e9tection des signaux faibles.<\/p>\n<p>Lorsque le param\u00e8tre alpha est positionn\u00e9 au seuil strict de 0.01, le quantile th\u00e9orique de la distribution de Student mobilis\u00e9 dans l&rsquo;\u00e9quation de Grubbs augmente substantiellement. La valeur critique G_crit qui en r\u00e9sulte est repouss\u00e9e vers des grandeurs plus extr\u00eames. Ce resserrement diminue drastiquement le risque d&rsquo;erreur de type I : seules les observations pr\u00e9sentant une d\u00e9viation colossale par rapport \u00e0 la moyenne seront d\u00e9clar\u00e9es aberrantes. Ce profil conservateur est particuli\u00e8rement pr\u00e9conis\u00e9 dans les contextes o\u00f9 l&rsquo;exclusion d&rsquo;une donn\u00e9e engendre des cons\u00e9quences th\u00e9oriques ou \u00e9thiques majeures, ou lorsque la variabilit\u00e9 naturelle de la population est r\u00e9put\u00e9e importante.<\/p>\n<p>\u00c0 l&rsquo;inverse, une configuration \u00e0 alpha \u00e9gal \u00e0 0.10 abaisse le seuil critique d&rsquo;acceptation, augmentant la propension du test \u00e0 rejeter l&rsquo;hypoth\u00e8se nulle. Si cette sensibilit\u00e9 accrue minimise les risques de conserver des mesures contamin\u00e9es (r\u00e9duction de l&rsquo;erreur de type II), elle accro\u00eet corr\u00e9lativement la probabilit\u00e9 d&rsquo;\u00e9laguer de mani\u00e8re injustifi\u00e9e des donn\u00e9es valides situ\u00e9es dans les marges physiologiques de la population. L&rsquo;analyste doit imp\u00e9rativement expliciter et justifier a priori la s\u00e9lection de son seuil alpha, \u00e9vitant tout ajustement opportuniste post-hoc destin\u00e9 \u00e0 conformer les donn\u00e9es aux attentes th\u00e9oriques de l&rsquo;exp\u00e9rimentateur.<\/p>\n<h2>8. Ex\u00e9cution du test de Grubbs unilat\u00e9ral pour extr\u00eames minimaux et maximaux<\/h2>\n<h3>8.1 D\u00e9tection cibl\u00e9e de la valeur maximale<\/h3>\n<p>Dans de nombreuses configurations exp\u00e9rimentales, l&rsquo;asym\u00e9trie fonctionnelle des processus psychologiques sous-jacents conf\u00e8re un statut privil\u00e9gi\u00e9 \u00e0 l&rsquo;une des extr\u00e9mit\u00e9s du spectre d&rsquo;observation. L&rsquo;\u00e9tude des latences d&rsquo;ex\u00e9cution dans les paradigmes cognitifs en fournit une illustration frappante. Alors qu&rsquo;un temps de r\u00e9action ne peut pas descendre en de\u00e7\u00e0 d&rsquo;un seuil physiologique incompressible d\u00e9termin\u00e9 par la conduction nerveuse aff\u00e9rente et eff\u00e9rente, il peut en revanche s&rsquo;\u00e9tirer de mani\u00e8re consid\u00e9rable sous l&rsquo;effet d&rsquo;un rel\u00e2chement attentionnel brutal ou d&rsquo;une divagation mentale du sujet exp\u00e9rimental.<\/p>\n<p>Dans ce contexte th\u00e9orique pr\u00e9cis, tester l&rsquo;\u00e9mergence d&rsquo;une valeur anormalement basse est d\u00e9nu\u00e9 de pertinence substantielle : seule la queue sup\u00e9rieure de la distribution abrite les art\u00e9facts d&rsquo;inattention suspect\u00e9s. Le test unilat\u00e9ral ciblant la valeur maximale r\u00e9pond pr\u00e9cis\u00e9ment \u00e0 cet imp\u00e9ratif d&rsquo;ajustement conceptuel. La biblioth\u00e8que <em>outlier_utils<\/em> encapsule cette proc\u00e9dure d\u00e9di\u00e9e par le biais de la fonction sp\u00e9cialis\u00e9e <em>smirnov_grubbs.max_test()<\/em>, compl\u00e9t\u00e9e par son homologue positionnelle <em>smirnov_grubbs.max_test_indices()<\/em>.<\/p>\n<p>En concentrant l&rsquo;int\u00e9gralit\u00e9 de la r\u00e9gion de rejet probabiliste sur la seule d\u00e9viation positive extr\u00eame, la fonction <em>max_test()<\/em> compare la statistique G_max exclusivement au quantile de Student unilat\u00e9ral. Ce dispositif focalis\u00e9 dispense l&rsquo;analyse de fractionner le risque alpha sur les deux extr\u00e9mit\u00e9s de la distribution, optimisant ainsi substantiellement la puissance du test pour isoler une latence anormale sans subir la p\u00e9nalit\u00e9 statistique d&rsquo;une incertitude directionnelle non fond\u00e9e.<\/p>\n<h3>8.2 D\u00e9tection cibl\u00e9e de la valeur minimale<\/h3>\n<p>\u00c0 l&rsquo;autre extr\u00e9mit\u00e9 de la dynamique comportementale se situe le sc\u00e9nario des r\u00e9ponses motrices pr\u00e9matur\u00e9es, des anticipations involontaires ou des dysfonctionnements instrumentaux provoquant un d\u00e9clenchement quasi instantan\u00e9 des capteurs temporels. Lors d&rsquo;une \u00e9preuve psychophysique imposant un choix forc\u00e9 complexe, l&rsquo;observation d&rsquo;un temps de r\u00e9ponse infinit\u00e9simal, par exemple cinquante millisecondes, est physiologiquement incompatible avec la dur\u00e9e incompressible du traitement sensoriel visuel et de la programmation motrice intracorticale. Ces points repr\u00e9sentent de toute \u00e9vidence des anticipations accidentelles survenues avant m\u00eame le traitement cognitif de l&rsquo;information.<\/p>\n<p>Pour traquer et \u00e9liminer s\u00e9lectivement ces anomalies par le bas sans sonder la variabilit\u00e9 de la queue sup\u00e9rieure, le chercheur mobilise le test de Grubbs unilat\u00e9ral minimal. Au sein du paquet <em>outlier_utils<\/em>, cette op\u00e9ration est orchestr\u00e9e par les fonctions jumelles <em>smirnov_grubbs.min_test()<\/em> et <em>smirnov_grubbs.min_test_indices()<\/em>. L&rsquo;algorithme \u00e9value alors sp\u00e9cifiquement la m\u00e9trique G_min, calculant l&rsquo;\u00e9cart entre la moyenne d&rsquo;\u00e9chantillon et la borne minimale observ\u00e9e.<\/p>\n<p>L&rsquo;int\u00e9r\u00eat de cette approche unilat\u00e9rale minimale r\u00e9side dans sa neutralit\u00e9 totale vis-\u00e0-vis des \u00e9ventuelles valeurs \u00e9lev\u00e9es de l&rsquo;\u00e9chantillon. M\u00eame si une observation haute pr\u00e9sente un \u00e9tirement marqu\u00e9, celle-ci n&rsquo;interf\u00e9rera aucunement avec le verdict port\u00e9 sur la valeur minimale suspecte, hormis par sa contribution m\u00e9canique \u00e0 la variance d&rsquo;\u00e9chantillon. Cette focalisation chirurgicale garantit que l&rsquo;analyste teste une hypoth\u00e8se de recherche rigoureusement circonscrite aux m\u00e9canismes d&rsquo;invalidation motrice pr\u00e9coce.<\/p>\n<h3>8.3 Comparaison fonctionnelle entre d\u00e9marches unilat\u00e9rales et bilat\u00e9rales<\/h3>\n<p>La d\u00e9cision m\u00e9thodologique d&rsquo;opter pour un test de Grubbs bilat\u00e9ral ou unilat\u00e9ral ne rel\u00e8ve pas d&rsquo;une simple commodit\u00e9 computationnelle, mais engage des consid\u00e9rations \u00e9pist\u00e9mologiques profondes relatives \u00e0 la formulation des hypoth\u00e8ses scientifiques. Sur le plan purement statistique, le test bilat\u00e9ral est une proc\u00e9dure conservatrice omnibus : en r\u00e9partissant la probabilit\u00e9 de rejet alpha \u00e0 parts \u00e9gales sur les deux queues de la distribution (alpha divis\u00e9 par deux \u00e0 chaque extr\u00e9mit\u00e9), il impose une barri\u00e8re critique th\u00e9orique plus \u00e9lev\u00e9e pour d\u00e9clarer la significativit\u00e9 d&rsquo;un \u00e9cart.<\/p>\n<p>Par cons\u00e9quent, lorsque l&rsquo;exp\u00e9rimentateur poss\u00e8de des justifications th\u00e9oriques inattaquables prescrivant qu&rsquo;une valeur aberrante ne peut survenir que dans une direction donn\u00e9e, le recours au test bilat\u00e9ral engendre une perte de puissance statistique \u00e9vidente. La statistique observ\u00e9e pourrait \u00e9chouer de peu \u00e0 franchir le seuil critique bilat\u00e9ral plus exigeant, g\u00e9n\u00e9rant une erreur de type II alors m\u00eame que le test unilat\u00e9ral concordant aurait l\u00e9gitimement conclu au rejet de l&rsquo;hypoth\u00e8se nulle d&rsquo;homog\u00e9n\u00e9it\u00e9.<\/p>\n<p>N\u00e9anmoins, l&rsquo;adoption d&rsquo;un test unilat\u00e9ral exige une probit\u00e9 scientifique absolue. Le choix de la directionnalit\u00e9 doit imp\u00e9rativement \u00eatre arr\u00eat\u00e9 a priori, dans le cadre d&rsquo;un protocole d&rsquo;analyse pr\u00e9enregistr\u00e9 ou fond\u00e9 sur des contraintes physiques et physiologiques universelles. Basculer vers un test unilat\u00e9ral a posteriori, apr\u00e8s avoir constat\u00e9 empiriquement qu&rsquo;un point d\u00e9viant extr\u00eame ne franchit pas le seuil bilat\u00e9ral canonique, constitue une infraction m\u00e9thodologique majeure assimilable au dragage de donn\u00e9es (<a rel=\"noopener noreferrer\" href=\"https:\/\/www.nature.com\/articles\/d41586-019-00857-9\">p-hacking<\/a>). Cette pratique inflationniste invalide le contr\u00f4le nominal de l&rsquo;erreur de premi\u00e8re esp\u00e8ce et vicie l&rsquo;ensemble des d\u00e9ductions inf\u00e9rentielles de la recherche.<\/p>\n<h2>9. Impl\u00e9mentation personnalis\u00e9e du test de Grubbs avec NumPy et SciPy<\/h2>\n<h3>9.1 D\u00e9veloppement de la fonction de calcul de la statistique G<\/h3>\n<p>Bien que le recours \u00e0 des modules sp\u00e9cialis\u00e9s tiers simplifie les flux de travail standards, la programmation personnalis\u00e9e du test de Grubbs \u00e0 partir des biblioth\u00e8ques fondamentales <em>NumPy<\/em> et <em>SciPy<\/em> pr\u00e9sente des vertus heuristiques et fonctionnelles majeures. Elle permet \u00e0 l&rsquo;analyste de ma\u00eetriser l&rsquo;int\u00e9gralit\u00e9 de la machinerie math\u00e9matique sous-jacente, de s&rsquo;affranchir de d\u00e9pendances externes parfois mal maintenues, et d&rsquo;adapter sur mesure le traitement des m\u00e9triques annexes.<\/p>\n<p>La premi\u00e8re composante de cette impl\u00e9mentation algorithmique consiste \u00e0 structurer le calcul vectoris\u00e9 de la statistique G. \u00c0 partir d&rsquo;un tableau NumPy d&rsquo;observations num\u00e9riques continues, il convient d&rsquo;\u00e9valuer la moyenne arithm\u00e9tique par la m\u00e9thode <em>mean()<\/em>, puis d&rsquo;estimer l&rsquo;\u00e9cart-type d&rsquo;\u00e9chantillon par la m\u00e9thode <em>std()<\/em> en veillant scrupuleusement \u00e0 fixer le param\u00e8tre des degr\u00e9s de libert\u00e9 diff\u00e9rentiels, not\u00e9 <em>ddof<\/em>, \u00e0 la valeur enti\u00e8re 1. Ce r\u00e9glage est capital : il instancie le diviseur n moins un garantissant l&rsquo;absence de biais de l&rsquo;estimateur de dispersion, condition sine qua non de l&rsquo;exactitude de la th\u00e9orie de Grubbs.<\/p>\n<p>Une fois ces scalaires obtenus, le code \u00e9value le vecteur des \u00e9carts absolus en retranchant la moyenne arithm\u00e9tique de chaque observation, puis en appliquant la fonction de valeur absolue <em>abs()<\/em> de NumPy. La recherche de l&rsquo;\u00e9cart maximal s&rsquo;op\u00e8re par l&rsquo;instruction <em>max()<\/em>, tandis que la fonction <em>argmax()<\/em> identifie l&rsquo;index positionnel pr\u00e9cis de l&rsquo;observation responsable de cette d\u00e9viation centrom\u00e9trique. Le ratio de cette d\u00e9viation maximale sur l&rsquo;\u00e9cart-type non biais\u00e9 d\u00e9livre alors la statistique scalaire G empirique avec une pr\u00e9cision en virgule flottante maximale.<\/p>\n<h3>9.2 Programmation de la valeur critique \u00e0 partir de la loi de Student<\/h3>\n<p>La seconde phase de la construction algorithmique autonome concerne la d\u00e9rivation exacte de la valeur critique th\u00e9orique. Cette \u00e9tape requiert l&rsquo;\u00e9valuation de la fonction quantile \u2014 \u00e9galement d\u00e9sign\u00e9e sous le terme de fonction de point de pourcentage (Percent Point Function, PPF) \u2014 de la loi de Student univari\u00e9e. Le module <em>scipy.stats<\/em> h\u00e9berge cette m\u00e9thode sous la fonction <em>t.ppf<\/em>, qui calcule la valeur d&rsquo;abscisse correspondant \u00e0 une probabilit\u00e9 cumul\u00e9e sp\u00e9cifi\u00e9e pour un nombre d\u00e9fini de degr\u00e9s de libert\u00e9.<\/p>\n<p>Pour transposer fid\u00e8lement la formule de Frank Grubbs en code Python, l&rsquo;analyste programme l&rsquo;ajustement probabiliste du seuil alpha. Dans la configuration d&rsquo;un test bilat\u00e9ral \u00e0 un niveau de risque unitaire alpha, l&rsquo;argument de probabilit\u00e9 transmis \u00e0 la fonction <em>t.ppf<\/em> s&rsquo;\u00e9l\u00e8ve \u00e0 1 moins alpha divis\u00e9 par le double de la taille d&rsquo;\u00e9chantillon n. Les degr\u00e9s de libert\u00e9 sont quant \u00e0 eux rigoureusement param\u00e9tr\u00e9s \u00e0 n moins deux. Le scalaire r\u00e9sultant correspond au quantile t exact de l&rsquo;\u00e9quation.<\/p>\n<p>Le calcul de la valeur critique G_crit s&rsquo;ach\u00e8ve par l&rsquo;application de la transformation alg\u00e9brique compl\u00e8te : multiplication du terme (n moins un) divis\u00e9 par la racine carr\u00e9e de n, par la racine carr\u00e9e du ratio articulant le carr\u00e9 du quantile t et la somme des degr\u00e9s de libert\u00e9 n moins deux et du carr\u00e9 de ce m\u00eame quantile t. Parall\u00e8lement, l&rsquo;impl\u00e9mentation de la fonction de survie de Student (<em>t.sf<\/em>) permet de calculer par inversion analytique la p-valeur exacte associ\u00e9e \u00e0 la statistique G observ\u00e9e, fournissant une m\u00e9trique probabiliste continue d&rsquo;une pr\u00e9cision totale.<\/p>\n<h3>9.3 Encapsulation dans une classe Python orient\u00e9e objet<\/h3>\n<p>Afin de conf\u00e9rer \u00e0 cette architecture une modularit\u00e9 logicielle irr\u00e9prochable conforme aux standards modernes de la programmation orient\u00e9e objet, l&rsquo;ensemble des m\u00e9thodes fonctionnelles gagne \u00e0 \u00eatre encapsul\u00e9 au sein d&rsquo;une classe d\u00e9di\u00e9e, que l&rsquo;on peut d\u00e9signer sous le nom de <em>GrubbsDetector<\/em>. Cette abstraction structurelle r\u00e9unit l&rsquo;\u00e9tat interne des donn\u00e9es, les param\u00e8tres de calibration statistique et les r\u00e9sultats d&rsquo;inf\u00e9rence au sein d&rsquo;une m\u00eame entit\u00e9 coh\u00e9rente.<\/p>\n<p>Le constructeur de cette classe prend en charge l&rsquo;initialisation du seuil de signification alpha ainsi que la sp\u00e9cification de l&rsquo;orientation directionnelle du test (bilat\u00e9rale, unilat\u00e9rale maximale ou unilat\u00e9rale minimale). Une m\u00e9thode ma\u00eetresse, traditionnellement d\u00e9nomm\u00e9e <em>fit()<\/em> ou <em>detect()<\/em>, accepte en entr\u00e9e la collection de donn\u00e9es empiriques, d\u00e9clenche la validation d\u00e9fensive pr\u00e9alable des conditions d&rsquo;\u00e9chantillonnage, orchestre le calcul de la statistique G, d\u00e9rive la valeur critique th\u00e9orique correspondante et calcule la p-valeur de l&rsquo;observation suspecte.<\/p>\n<p>Cette structure de classe int\u00e8gre des m\u00e9thodes d&rsquo;acc\u00e8s permettant d&rsquo;extraire un bilan statistique d\u00e9taill\u00e9 sous la forme d&rsquo;un dictionnaire ou d&rsquo;un rapport textuel format\u00e9, articulant clairement la taille de l&rsquo;\u00e9chantillon, la moyenne, la variance, la statistique G calcul\u00e9e, le seuil critique th\u00e9orique, la p-valeur exacte et la conclusion logique d&rsquo;acceptation ou de rejet de l&rsquo;hypoth\u00e8se nulle. La confrontation syst\u00e9matique des r\u00e9sultats produits par cette classe maison avec les sorties standard du paquet <em>outlier_utils<\/em> offre une proc\u00e9dure de validation crois\u00e9e garantissant l&rsquo;exactitude absolue du moteur computationnel.<\/p>\n<h2>10. Proc\u00e9dures it\u00e9ratives pour la d\u00e9tection de multiples valeurs aberrantes<\/h2>\n<h3>10.1 M\u00e9canisme de r\u00e9application s\u00e9quentielle du test de Grubbs<\/h3>\n<p>Dans les contextes d&rsquo;investigation empirique r\u00e9elle, il est fr\u00e9quent que les donn\u00e9es soient contamin\u00e9es non pas par une anomalie isol\u00e9e, mais par plusieurs perturbations concurrentes. Face \u00e0 cette situation, une pratique analytique largement document\u00e9e dans la litt\u00e9rature computationnelle consiste \u00e0 recourir \u00e0 une proc\u00e9dure it\u00e9rative ou r\u00e9cursive du test de Grubbs. Ce protocole it\u00e9ratif vise \u00e0 purger s\u00e9quentiellement l&rsquo;\u00e9chantillon de ses donn\u00e9es polluantes jusqu&rsquo;\u00e0 ce que la distribution r\u00e9siduelle satisfasse pleinement \u00e0 l&rsquo;hypoth\u00e8se d&rsquo;homog\u00e9n\u00e9it\u00e9.<\/p>\n<p>Le d\u00e9roulement algorithmique de cette proc\u00e9dure s&rsquo;ordonne selon une boucle s\u00e9quentielle conditionnelle contr\u00f4l\u00e9e par le verdict statistique. \u00c0 chaque cycle de calcul, le test de Grubbs standard est appliqu\u00e9 \u00e0 l&rsquo;\u00e9chantillon courant de dimension n. Si la statistique G observ\u00e9e d\u00e9passe le seuil critique th\u00e9orique G_crit calcul\u00e9 pour cet effectif sp\u00e9cifique, l&rsquo;observation suspecte est d\u00e9clar\u00e9e aberrante et formellement extraite du vecteur de donn\u00e9es. Le sous-ensemble restant, dont l&rsquo;effectif s&rsquo;\u00e9tablit d\u00e9sormais \u00e0 n moins un, fait alors l&rsquo;objet d&rsquo;un recalcul int\u00e9gral de sa moyenne arithm\u00e9tique et de son \u00e9cart-type d&rsquo;\u00e9chantillon.<\/p>\n<p>Cette r\u00e9\u00e9valuation compl\u00e8te des estimateurs est absolument indispensable : l&rsquo;extraction d&rsquo;un point extr\u00eame modifie substantiellement le barycentre et compresse imm\u00e9diatement la variance globale du groupe, ce qui modifie la sensibilit\u00e9 du test pour les points restants. Une nouvelle valeur critique th\u00e9orique est alors calcul\u00e9e, r\u00e9ajust\u00e9e sur les nouveaux degr\u00e9s de libert\u00e9 n moins trois. La boucle d&rsquo;it\u00e9ration poursuit son ex\u00e9cution s\u00e9quentielle jusqu&rsquo;\u00e0 ce qu&rsquo;aucune observation r\u00e9siduelle ne franchisse le seuil critique th\u00e9orique recalcul\u00e9, ou jusqu&rsquo;\u00e0 ce que la taille de l&rsquo;\u00e9chantillon r\u00e9siduel atteigne la borne minimale de recevabilit\u00e9 statistique.<\/p>\n<h3>10.2 Le pi\u00e8ge statistique du masquage (Masking) et de l&rsquo;inondation (Swamping)<\/h3>\n<p>Bien que l&rsquo;automatisation d&rsquo;un test de Grubbs r\u00e9cursif soit s\u00e9duisante de prime abord, cette pratique expose l&rsquo;analyste \u00e0 deux pathologies statistiques redoutables intrins\u00e8ques aux proc\u00e9dures d&rsquo;extr\u00eames studentis\u00e9s : le ph\u00e9nom\u00e8ne de masquage (masking) et le ph\u00e9nom\u00e8ne d&rsquo;inondation (swamping). Ces dysfonctionnements conceptuels imposent une circonspection absolue quant \u00e0 la g\u00e9n\u00e9ralisation aveugle de la boucle de Grubbs it\u00e9rative.<\/p>\n<p>Le masquage statistique survient typiquement en pr\u00e9sence d&rsquo;un cluster d&rsquo;observations aberrantes de magnitudes comparables. Par exemple, si une exp\u00e9rience chronom\u00e9trique enregistre deux latences art\u00e9factuelles massives de 4500 et 4600 millisecondes au sein d&rsquo;une distribution gravitant habituellement autour de 500 millisecondes, la pr\u00e9sence conjointe de ces deux extr\u00eames gonfle dramatiquement l&rsquo;\u00e9cart-type global de l&rsquo;\u00e9chantillon. En raison de ce d\u00e9nominateur d\u00e9mesur\u00e9ment large, le score G calcul\u00e9 pour la valeur la plus extr\u00eame (4600 ms) se trouve comprim\u00e9 et \u00e9choue \u00e0 franchir le seuil critique G_crit. La seconde aberration masque litt\u00e9ralement la premi\u00e8re, conduisant l&rsquo;analyste \u00e0 accepter \u00e0 tort l&rsquo;hypoth\u00e8se nulle et \u00e0 conserver les deux art\u00e9facts.<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;inondation ou surd\u00e9tection (swamping) d\u00e9signe la situation o\u00f9 la pr\u00e9sence d&rsquo;une ou plusieurs valeurs aberrantes majeures d\u00e9place le barycentre et distord tellement la structure d&rsquo;\u00e9chelle qu&rsquo;une observation parfaitement saine et l\u00e9gitime se retrouve, apr\u00e8s l&rsquo;extraction s\u00e9quentielle d&rsquo;un premier point, qualifi\u00e9e d&rsquo;aberrante lors du cycle suivant. L&rsquo;effet de cascade expulse alors injustement des donn\u00e9es authentiques, r\u00e9duisant artificiellement la dispersion intrins\u00e8que du ph\u00e9nom\u00e8ne biologique ou cognitif mod\u00e9lis\u00e9.<\/p>\n<h3>10.3 Alternative recommand\u00e9e : Le test ESD g\u00e9n\u00e9ralis\u00e9 de Rosner<\/h3>\n<p>Pour contourner de mani\u00e8re math\u00e9matiquement rigoureuse les impasses inh\u00e9rentes au masquage et \u00e0 l&rsquo;inondation, Bernard Rosner a d\u00e9velopp\u00e9 en 1983 une extension fondamentale de la th\u00e9orie de Grubbs : le test de d\u00e9viation extr\u00eame studentis\u00e9e g\u00e9n\u00e9ralis\u00e9e, universellement d\u00e9sign\u00e9 par l&rsquo;acronyme GESD (Generalized Extreme Studentized Deviate). Cette proc\u00e9dure est express\u00e9ment con\u00e7ue pour d\u00e9tecter jusqu&rsquo;\u00e0 un nombre pr\u00e9d\u00e9termin\u00e9 k de valeurs aberrantes au sein d&rsquo;une distribution normale, tout en r\u00e9sistant remarquablement au masquage par paires.<\/p>\n<p>La sup\u00e9riorit\u00e9 algorithmique du test de Rosner repose sur sa structure d\u00e9cisionnelle invers\u00e9e. Au lieu de tester et d&rsquo;arr\u00eater la proc\u00e9dure d\u00e8s qu&rsquo;un point n&rsquo;est pas significatif, l&rsquo;algorithme de Rosner extrait s\u00e9quentiellement les k observations les plus distantes de la moyenne, r\u00e9estimant les param\u00e8tres de dispersion \u00e0 chaque \u00e9tape, pour ensuite \u00e9valuer la significativit\u00e9 statistique en partant de la k-i\u00e8me anomalie candidate en remontant vers la premi\u00e8re. Si la statistique calcul\u00e9e pour la transition d&rsquo;ordre r est sup\u00e9rieure \u00e0 sa valeur critique associ\u00e9e, alors exactement r valeurs sont formellement d\u00e9clar\u00e9es aberrantes d&rsquo;un seul bloc, neutralisant ainsi totalement l&rsquo;effet de masquage.<\/p>\n<p>En Python, l&rsquo;impl\u00e9mentation du test GESD de Rosner peut \u00eatre ais\u00e9ment programm\u00e9e \u00e0 l&rsquo;aide de boucles matricielles mobilisant <em>SciPy<\/em>, ou ex\u00e9cut\u00e9e par l&rsquo;entremise de fonctions d\u00e9di\u00e9es issues de biblioth\u00e8ques sp\u00e9cialis\u00e9es comme <em>scikit-posthocs<\/em>. D\u00e8s lors qu&rsquo;un jeu de donn\u00e9es exp\u00e9rimentales complexes laisse soup\u00e7onner la pr\u00e9sence potentielle de contaminations multiples ou asym\u00e9triques, l&rsquo;adoption du protocole GESD de Rosner doit imp\u00e9rativement \u00eatre pr\u00e9f\u00e9r\u00e9e \u00e0 la r\u00e9application na\u00efve du test de Grubbs simple.<\/p>\n<h2>11. Interpr\u00e9tation psychom\u00e9trique et gouvernance des donn\u00e9es aberrantes<\/h2>\n<h3>11.1 Dilemme analytique : suppression, imputation ou r\u00e9tention<\/h3>\n<p>L&rsquo;identification formelle d&rsquo;une valeur aberrante par le test de Grubbs marque le commencement d&rsquo;une d\u00e9lib\u00e9ration m\u00e9thodologique complexe, et non son aboutissement. La qualification statistique d&rsquo;une anomalie n&#8217;emporte pas automatiquement l&rsquo;injonction de son effacement pur et simple du fichier de donn\u00e9es. L&rsquo;exp\u00e9rimentateur fait face \u00e0 un dilemme analytique fondamental articul\u00e9 autour de trois options strat\u00e9giques : l&rsquo;\u00e9lagage (suppression d\u00e9finitive), l&rsquo;imputation probabiliste et la r\u00e9tention assortie d&rsquo;un traitement de r\u00e9duction d&rsquo;influence.<\/p>\n<p>L&rsquo;\u00e9limination d\u00e9finitive d&rsquo;une observation ne doit \u00eatre envisag\u00e9e que si l&rsquo;analyste dispose d&rsquo;\u00e9l\u00e9ments contextuels tangibles \u00e9tablissant que la valeur r\u00e9sulte d&rsquo;une panne instrumentale av\u00e9r\u00e9e, d&rsquo;une erreur mat\u00e9rielle de saisie ou d&rsquo;une violation indiscutable du protocole exp\u00e9rimental par le participant. Supprimer une observation simplement parce qu&rsquo;elle d\u00e9vie statistiquement de la trajectoire attendue constitue un raisonnement circulaire destructeur pour la d\u00e9couvrabilit\u00e9 scientifique, amputant l&rsquo;\u00e9chantillon de sa variance naturelle authentique.<\/p>\n<p>Lorsque la pr\u00e9servation de la taille d&rsquo;\u00e9chantillon initiale est imp\u00e9rative, par exemple dans le cadre de devis exp\u00e9rimentaux r\u00e9p\u00e9t\u00e9s balanc\u00e9s, les techniques d&rsquo;imputation adapt\u00e9es aux donn\u00e9es psychom\u00e9triques offrent une alternative mesur\u00e9e. Le remplacement par la m\u00e9diane de la condition ou l&rsquo;imputation multiple par \u00e9quations cha\u00een\u00e9es permettent de combler l&rsquo;absence tout en minimisant la distorsion de la variance. Parall\u00e8lement, les m\u00e9thodes de <a rel=\"noopener noreferrer\" href=\"https:\/\/online.stat.psu.edu\/stat500\/lesson\/1\/1.5\/1.5.1\">winsorisation<\/a> \u2014 consistant \u00e0 ramener la valeur extr\u00eame \u00e0 la borne du quantile acceptable le plus proche (par exemple le 95e ou 99e centile) \u2014 permettent de neutraliser l&rsquo;effet de levier disproportionn\u00e9 de l&rsquo;anomalie sur la moyenne et l&rsquo;\u00e9cart-type, sans sacrifier l&rsquo;existence informationnelle du point d&rsquo;\u00e9valuation.<\/p>\n<h3>11.2 Typologie des causes d&rsquo;aberrance en sciences du comportement<\/h3>\n<p>Pour \u00e9clairer cette d\u00e9cision de gouvernance de donn\u00e9es, le chercheur en sciences comportementales doit proc\u00e9der \u00e0 une analyse \u00e9tiologique fine de l&rsquo;aberrance statistique. Trois grandes cat\u00e9gories explicatives pr\u00e9valent dans les laboratoires empiriques : les d\u00e9faillances technologiques mat\u00e9rielles, la non-conformit\u00e9 comportementale des agents humains, et l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 biologique ou psychologique intrins\u00e8que.<\/p>\n<p>Les d\u00e9faillances d&rsquo;enregistrement mat\u00e9riel englobent les dysfonctionnements \u00e9lectroniques des interfaces de r\u00e9ponse (rebonds de contacts m\u00e9caniques sur les claviers de r\u00e9ponse chronom\u00e9triques), les d\u00e9synchronisations logicielles entre le rafra\u00eechissement de l&rsquo;affichage visuel et l&rsquo;horloge syst\u00e8me du microprocesseur, ou les micro-coupures de capteurs physiologiques. De tels artefacts doivent imp\u00e9rativement \u00eatre exclus du corpus analytique puisqu&rsquo;ils n&rsquo;entretiennent aucun lien de causalit\u00e9 avec l&rsquo;architecture psychologique du sujet.<\/p>\n<p>La non-conformit\u00e9 aux consignes de la t\u00e2che repr\u00e9sente la seconde source de perturbation. Elle documente les \u00e9pisodes o\u00f9 le participant cesse d&rsquo;allouer ses ressources attentionnelles au protocole : somnolence passag\u00e8re, interruption verbale avec l&rsquo;exp\u00e9rimentateur, ou s\u00e9lection syst\u00e9matique et aveugle de la m\u00eame r\u00e9ponse motrice pour achever l&rsquo;\u00e9preuve au plus vite. Si ces comportements invalident la repr\u00e9sentativit\u00e9 du proc\u00e8s cognitif, leur rejet doit \u00eatre adoss\u00e9 \u00e0 des crit\u00e8res de contr\u00f4le de validit\u00e9 standardis\u00e9s.<\/p>\n<p>Enfin, la troisi\u00e8me cause \u2014 et la plus d\u00e9licate \u00e0 arbitrer \u2014 r\u00e9side dans la variabilit\u00e9 humaine authentique. Un participant peut manifester des capacit\u00e9s cognitives prodigieuses ou, \u00e0 l&rsquo;inverse, appartenir \u00e0 une sous-population clinique non diagnostiqu\u00e9e caract\u00e9ris\u00e9e par des d\u00e9ficits sp\u00e9cifiques. Dans cette configuration, la donn\u00e9e extr\u00eame constitue pr\u00e9cis\u00e9ment le c\u0153ur de la d\u00e9couverte scientifique : l&rsquo;exclure reviendrait \u00e0 op\u00e9rer une censure m\u00e9thodologique masquant l&rsquo;existence d&rsquo;endoph\u00e9notypes cognitifs originaux.<\/p>\n<h3>11.3 Exigences de transparence et de reproductibilit\u00e9 selon les normes APA<\/h3>\n<p>L&rsquo;\u00e9thique scientifique contemporaine, renforc\u00e9e par la crise de la reproductibilit\u00e9 en psychologie et les recommandations formelles de l&rsquo;American Psychological Association (<a rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\/\">normes APA<\/a>, 7e \u00e9dition), proscrit cat\u00e9goriquement les pratiques d&rsquo;exclusion opaque de donn\u00e9es. Tout protocole d&rsquo;assainissement statistique mobilisant le test de Grubbs doit faire l&rsquo;objet d&rsquo;une divulgation exhaustive au sein de la section M\u00e9thode des manuscrits soumis \u00e0 \u00e9valuation.<\/p>\n<p>Les directives de notification imposent au chercheur de sp\u00e9cifier sans ambigu\u00eft\u00e9 :<\/p>\n<ul>\n<li>Le mod\u00e8le formel de test mobilis\u00e9 (test de Grubbs univari\u00e9 simple, bilat\u00e9ral ou unilat\u00e9ral, ou test de Rosner) ;<\/li>\n<li>Le seuil de d\u00e9cision probabiliste alpha fix\u00e9 a priori, assorti de ses justifications m\u00e9thodologiques ;<\/li>\n<li>Le nombre exact d&rsquo;observations identifi\u00e9es comme d\u00e9viantes, leur magnitude m\u00e9trique individuelle et le pourcentage de donn\u00e9es initiales que repr\u00e9sente cette exclusion par rapport \u00e0 l&rsquo;\u00e9chantillon brut ;<\/li>\n<li>Le sort r\u00e9serv\u00e9 aux anomalies d\u00e9tect\u00e9es (suppression d\u00e9finitive, imputation, winsorisation ou mod\u00e9lisation robuste).<\/li>\n<\/ul>\n<p>De surcro\u00eet, la d\u00e9monstration de la robustesse des conclusions exp\u00e9rimentales exige la conduite syst\u00e9matique d&rsquo;une analyse de sensibilit\u00e9. L&rsquo;auteur doit pr\u00e9senter dans ses synth\u00e8ses inf\u00e9rentielles une comparaison synoptique des r\u00e9sultats statistiques fondamentaux calcul\u00e9s avec l&rsquo;ensemble des donn\u00e9es brutes versus calcul\u00e9s apr\u00e8s exclusion des valeurs aberrantes de Grubbs. Si la significativit\u00e9 statistique d&rsquo;un effet majeur s&rsquo;\u00e9vanouit lors de l&rsquo;inclusion d&rsquo;une seule observation d\u00e9viante, ce constat doit \u00eatre honn\u00eatement discut\u00e9. Enfin, le partage public du script de traitement Python au sein d&rsquo;un entrep\u00f4t de science ouverte certifi\u00e9, tel que l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/osf.io\/\">Open Science Framework (OSF)<\/a>, parach\u00e8ve la cha\u00eene de transparence et de reproductibilit\u00e9 computationnelle exig\u00e9e par la communaut\u00e9 savante.<\/p>\n<h2>12. Comparaison critique du test de Grubbs avec d&rsquo;autres m\u00e9thodes en Python<\/h2>\n<h3>12.1 Test de Grubbs versus Z-score standardis\u00e9 et Z-score robuste (MAD)<\/h3>\n<p>Pour appr\u00e9cier la singularit\u00e9 du test de Grubbs dans la panoplie des outils de pr\u00e9traitement, il convient de le confronter analytiquement aux m\u00e9thodes d&rsquo;\u00e9cart standardis\u00e9 courantes. L&rsquo;heuristique universelle du Z-score standard, qui qualifie d&rsquo;aberrante toute observation dont l&rsquo;\u00e9cart \u00e0 la moyenne d\u00e9passe arbitrairement deux ou trois \u00e9carts-types, souffre d&rsquo;une aporie circulaire majeure. La moyenne et l&rsquo;\u00e9cart-type mobilis\u00e9s pour calculer le Z-score int\u00e8grent pr\u00e9cis\u00e9ment en leur sein la valeur aberrante que l&rsquo;on cherche \u00e0 traquer. Cette contamination r\u00e9ciproque induit une r\u00e9sistance fallacieuse du Z-score, qui sous-estime syst\u00e9matiquement la v\u00e9ritable anormalit\u00e9 du point dans les petits \u00e9chantillons.<\/p>\n<p>Le test de Grubbs r\u00e9sout pr\u00e9cis\u00e9ment cette limite du Z-score conventionnel en mod\u00e9lisant la loi de distribution exacte de l&rsquo;\u00e9cart standardis\u00e9 maximal sous forme d&rsquo;une distribution de Student ajust\u00e9e pour la taille d&rsquo;\u00e9chantillon. Toutefois, face \u00e0 une contamination volumineuse, l&rsquo;estimation de la moyenne et de la variance de Grubbs s&rsquo;effondre \u00e9galement en raison d&rsquo;un point de rupture (\u00ab breakdown point \u00bb) th\u00e9orique extr\u00eamement faible : une unique valeur d\u00e9viante suffit \u00e0 influencer les deux estimateurs de base.<\/p>\n<p>C&rsquo;est ici qu&rsquo;intervient la sup\u00e9riorit\u00e9 des m\u00e9thodes robustes fond\u00e9es sur l&rsquo;\u00e9cart absolu \u00e0 la m\u00e9diane (Median Absolute Deviation, MAD), accessible dans l&rsquo;\u00e9cosyst\u00e8me Python via la fonction <em>scipy.stats.median_abs_deviation<\/em>. Le score Z modifi\u00e9 fond\u00e9 sur le MAD remplace la moyenne arithm\u00e9tique par la m\u00e9diane et l&rsquo;\u00e9cart-type par la m\u00e9diane des \u00e9carts absolus \u00e0 la m\u00e9diane. Dot\u00e9e d&rsquo;un point de rupture optimal atteignant 50 %, cette m\u00e9thodologie r\u00e9siste imperturbablement \u00e0 la pr\u00e9sence simultan\u00e9e de multiples valeurs extr\u00eames, offrant une r\u00e9silience math\u00e9matique que le cadre param\u00e9trique strict de Grubbs ne peut \u00e9galer.<\/p>\n<h3>12.2 Test de Grubbs versus approches d&rsquo;apprentissage automatique non supervis\u00e9<\/h3>\n<p>L&rsquo;essor de la science des donn\u00e9es et de l&rsquo;apprentissage machine (machine learning) a d\u00e9mocratis\u00e9 des algorithmes sophistiqu\u00e9s de d\u00e9tection d&rsquo;anomalies en environnement non supervis\u00e9, parmi lesquels <em>Isolation Forest<\/em> et l&rsquo;algorithme des facteurs d&rsquo;aberration locale (Local Outlier Factor, LOF), largement impl\u00e9ment\u00e9s au sein du module <em>sklearn.neighbors<\/em> et <em>sklearn.ensemble<\/em> de la biblioth\u00e8que Scikit-Learn.<\/p>\n<p>L&rsquo;algorithme <em>Isolation Forest<\/em> op\u00e8re par partitionnement r\u00e9cursif de l&rsquo;espace g\u00e9om\u00e9trique des donn\u00e9es \u00e0 l&rsquo;aide d&rsquo;arbres de d\u00e9cision al\u00e9atoires. Les anomalies structurelles, \u00e9tant quantitativement rares et distantes des concentrations denses d&rsquo;observations, sont isol\u00e9es \u00e0 des profondeurs d&rsquo;arbres nettement moindres que les donn\u00e9es nominales. De son c\u00f4t\u00e9, l&rsquo;algorithme LOF mesure l&rsquo;isolement local d&rsquo;un point en \u00e9valuant la densit\u00e9 volumique de ses k plus proches voisins dans l&rsquo;espace vectoriel d&rsquo;\u00e9chantillonnage.<\/p>\n<p>Bien que ces algorithmes de machine learning excellent dans le traitement de donn\u00e9es massives hautement multidimensionnelles o\u00f9 les anomalies r\u00e9sultent d&rsquo;interactions non lin\u00e9aires complexes entre de multiples variables, leur mobilisation face \u00e0 un vecteur univari\u00e9 issu d&rsquo;un devis psychologique exp\u00e9rimental s&rsquo;av\u00e8re souvent inappropri\u00e9e et surdimensionn\u00e9e. Ces mod\u00e8les fonctionnent comme des \u00ab bo\u00eetes noires \u00bb d\u00e9nu\u00e9es de formulation de p-valeur probabiliste exacte, interdisant le contr\u00f4le formel du risque alpha exig\u00e9 par la statistique inf\u00e9rentielle classique. Pour une s\u00e9rie univari\u00e9e de petite \u00e0 moyenne dimension r\u00e9pondant au postulat gaussien, le test de Grubbs demeure infiniment sup\u00e9rieur en termes de parcimonie, de tractabilit\u00e9 math\u00e9matique et de clart\u00e9 \u00e9pist\u00e9mologique.<\/p>\n<h3>12.3 Synth\u00e8se d\u00e9cisionnelle pour le chercheur en psychologie<\/h3>\n<p>Pour guider le praticien et le chercheur \u00e0 travers la complexit\u00e9 des choix algorithmiques de d\u00e9tection d&rsquo;anomalies en Python, l&rsquo;adoption d&rsquo;un arbre de d\u00e9cision m\u00e9thodologique standardis\u00e9 s&rsquo;av\u00e8re hautement b\u00e9n\u00e9fique. Cette synth\u00e8se op\u00e9rationnelle structure les \u00e9tapes incontournables garantissant une gouvernance de donn\u00e9es rigoureuse et reproductible :<\/p>\n<ol>\n<li><strong>\u00c9valuation de la dimensionnalit\u00e9 :<\/strong> Si l&rsquo;investigation porte sur une structure multivari\u00e9e non r\u00e9ductible, orienter la cha\u00eene de traitement vers des algorithmes non param\u00e9triques multidimensionnels (ex. Isolation Forest). Si l&rsquo;investigation est univari\u00e9e, poursuivre vers l&rsquo;\u00e9tape de diagnostic distributionnel.<\/li>\n<li><strong>Diagnostic de la normalit\u00e9 et de la taille d&rsquo;\u00e9chantillon :<\/strong> Soumettre la s\u00e9rie univari\u00e9e aux tests de Shapiro-Wilk et \u00e0 l&rsquo;examen du diagramme Q-Q plot via SciPy.\n<ul>\n<li>Si l&rsquo;\u00e9chantillon compte moins de sept observations (N &lt; 7), exclure formellement le test de Grubbs au profit du test de Dixon ou de m\u00e9thodes heuristiques prudentes.<\/li>\n<li>Si l&rsquo;\u00e9chantillon pr\u00e9sente une asym\u00e9trie marqu\u00e9e ou une leptokurticit\u00e9 violant le postulat gaussien, renoncer aux tests param\u00e9triques et d\u00e9ployer la d\u00e9tection robuste par le Z-score modifi\u00e9 fond\u00e9 sur le MAD ou la r\u00e8gle interquartile de Tukey.<\/li>\n<\/ul>\n<\/li>\n<li><strong>Nombre d&rsquo;anomalies pr\u00e9sum\u00e9es et directionnalit\u00e9 :<\/strong>\n<ul>\n<li>Si la s\u00e9rie est confirm\u00e9e gaussienne et que l&rsquo;hypoth\u00e8se m\u00e9thodologique ne cible qu&rsquo;une seule aberration suspecte, ex\u00e9cuter le test de Grubbs via <em>outlier_utils.smirnov_grubbs<\/em> (en choisissant rigoureusement le mode bilat\u00e9ral ou unilat\u00e9ral selon la nature des postulats a priori).<\/li>\n<li>Si des anomalies multiples sont soup\u00e7onn\u00e9es ou si un risque de masquage statistique est identifi\u00e9, \u00e9carter la boucle de Grubbs it\u00e9rative na\u00efve et mobiliser sans r\u00e9serve le test Generalized ESD de Rosner.<\/li>\n<\/ul>\n<\/li>\n<li><strong>Gouvernance \u00e9thique et archivage reproductible :<\/strong> Documenter l&rsquo;int\u00e9gralit\u00e9 du processus de d\u00e9tection, conduire obligatoirement une analyse de sensibilit\u00e9 comparant les mod\u00e8les inf\u00e9rentiels avec et sans les donn\u00e9es atypiques, expliciter les arbitrages (suppression, imputation, winsorisation) et d\u00e9poser les scripts Python de pr\u00e9traitement sur un d\u00e9p\u00f4t public en libre acc\u00e8s conform\u00e9ment aux normes scientifiques APA.<\/li>\n<\/ol>\n<p>En inscrivant l&rsquo;usage du test de Grubbs dans ce cadre m\u00e9thodologique \u00e9quilibr\u00e9, le chercheur s&rsquo;assure de conjuguer la puissance d&rsquo;analyse du langage Python avec l&rsquo;int\u00e9grit\u00e9 d\u00e9ontologique requise pour la production d&rsquo;un savoir scientifique robuste et p\u00e9renne.<\/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 rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/li>\n<li>D&rsquo;Agostino, R. B., &amp; Pearson, E. S. (1973). Tests for departure from normality. Empirical results for the distributions of b2 and sqrt(b1). <em>Biometrika<\/em>, 60(3), 613\u2013622. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/60.3.613\">https:\/\/doi.org\/10.1093\/biomet\/60.3.613<\/a><\/li>\n<li>Dixon, W. J. (1950). Analysis of extreme values. <em>The Annals of Mathematical Statistics<\/em>, 21(4), 488\u2013506. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aoms\/1177729747\">https:\/\/doi.org\/10.1214\/aoms\/1177729747<\/a><\/li>\n<li>Grubbs, F. E. (1950). Procedures for detecting outlying observations in samples. <em>The Annals of Mathematical Statistics<\/em>, 21(1), 27\u201358. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aoms\/1177729885\">https:\/\/doi.org\/10.1214\/aoms\/1177729885<\/a><\/li>\n<li>Grubbs, F. E. (1969). Procedures for detecting outlying observations in samples. <em>Technometrics<\/em>, 11(1), 1\u201321. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00401706.1969.10490657\">https:\/\/doi.org\/10.1080\/00401706.1969.10490657<\/a><\/li>\n<li>Leys, C., Ley, C., Klein, O., Bernard, P., &amp; Licata, L. (2013). Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median. <em>Journal of Experimental Social Psychology<\/em>, 49(4), 764\u2013766. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/j.jesp.2013.03.013\">https:\/\/doi.org\/10.1016\/j.jesp.2013.03.013<\/a><\/li>\n<li>National Institute of Standards and Technology. (2012). <em>NIST\/SEMATECH e-Handbook of Statistical Methods<\/em>. U.S. Department of Commerce. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.18434\/M32189\">https:\/\/doi.org\/10.18434\/M32189<\/a><\/li>\n<li>Ratcliff, R. (1993). Methods for dealing with reaction time outliers. <em>Psychological Bulletin<\/em>, 114(3), 510\u2013532. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0033-2909.114.3.510\">https:\/\/doi.org\/10.1037\/0033-2909.114.3.510<\/a><\/li>\n<li>Rosner, B. (1983). Percentage points for a generalized ESD many-outlier procedure. <em>Technometrics<\/em>, 25(2), 165\u2013172. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00401706.1983.10487848\">https:\/\/doi.org\/10.1080\/00401706.1983.10487848<\/a><\/li>\n<li>Shapiro, S. S., &amp; Wilk, M. B. (1965). An analysis of variance test for normality (complete samples). <em>Biometrika<\/em>, 52(3\/4), 591\u2013611. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/52.3-4.591\">https:\/\/doi.org\/10.1093\/biomet\/52.3-4.591<\/a><\/li>\n<li>Tukey, J. W. (1977). <em>Exploratory data analysis<\/em>. Addison-Wesley.<\/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., \u2026 SciPy 1.0 Contributors. (2020). SciPy 1.0: Fundamental algorithms for scientific computing in Python. <em>Nature Methods<\/em>, 17(3), 261\u2013272. <a 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<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour d\u00e9tecter les valeurs aberrantes en psychologie : th\u00e9orie, formulation et ex\u00e9cution pratique du test de Grubbs avec Python.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2465,2169],"tags":[16,2468,2231,4227,1963],"class_list":["post-3547","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-programmation-python","category-psychologie-quantitative","tag-psychometrie","tag-python","tag-statistiques-inferentielles","tag-test-de-grubbs","tag-valeurs-aberrantes"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/3547","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=3547"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/3547\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=3547"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=3547"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=3547"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}