{"id":9985,"date":"2026-10-14T18:02:35","date_gmt":"2026-10-14T18:02:35","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-identifier-valeurs-aberrantes-sas-exemple\/"},"modified":"2026-10-14T18:02:35","modified_gmt":"2026-10-14T18:02:35","slug":"comment-identifier-valeurs-aberrantes-sas-exemple","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-identifier-valeurs-aberrantes-sas-exemple","title":{"rendered":"Comment identifier les valeurs aberrantes dans SAS (avec exemple)"},"content":{"rendered":"<p>L&rsquo;int\u00e9grit\u00e9 analytique d&rsquo;une recherche quantitative repose fondamentalement sur la qualit\u00e9 m\u00e9trologique des observations recueillies et sur la rigueur du traitement statistique appliqu\u00e9 aux distributions empiriques. Dans les sciences empiriques contemporaines, qu&rsquo;il s&rsquo;agisse de la biom\u00e9trie, de l&rsquo;\u00e9conom\u00e9trie, de l&rsquo;\u00e9pid\u00e9miologie ou des sciences comportementales, l&rsquo;\u00e9mergence d&rsquo;observations atypiques &mdash; commun\u00e9ment d\u00e9sign\u00e9es sous le vocable de valeurs aberrantes ou <em>outliers<\/em> &mdash; constitue un ph\u00e9nom\u00e8ne omnipr\u00e9sent dont l&rsquo;impact sur l&rsquo;inf\u00e9rence statistique peut s&rsquo;av\u00e9rer d\u00e9vastateur. Une valeur aberrante ne se r\u00e9sume pas \u00e0 une simple curiosit\u00e9 descriptive ; elle incarne une d\u00e9viation substantielle par rapport \u00e0 la structure probabiliste sous-jacente du reste de l&rsquo;\u00e9chantillon, mena\u00e7ant la validit\u00e9 interne des estimations, distordant les erreurs-types et conduisant potentiellement \u00e0 des conclusions scientifiques fallacieuses.<\/p>\n<p>Le progiciel <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.sas.com\">SAS (Statistical Analysis System)<\/a> s&rsquo;est impos\u00e9 depuis plusieurs d\u00e9cennies comme le standard de r\u00e9f\u00e9rence au sein des industries r\u00e9glement\u00e9es, des agences de sant\u00e9 publique et des laboratoires acad\u00e9miques pour l&rsquo;analyse de donn\u00e9es complexes et massives. Gr\u00e2ce \u00e0 son architecture logicielle \u00e0 double composante &mdash; combinant la flexibilit\u00e9 proc\u00e9durale de l&rsquo;\u00e9tape <code>DATA<\/code> et la puissance calculatoire de proc\u00e9dures optimis\u00e9es telles que <code>PROC UNIVARIATE<\/code>, <code>PROC STDIZE<\/code> ou <code>PROC ROBUSTREG<\/code> &mdash; l&rsquo;environnement SAS offre un arsenal m\u00e9thodologique sans \u00e9quivalent pour isoler, diagnostiquer et traiter les observations discordantes. La ma\u00eetrise de ces outils computationnels exige n\u00e9anmoins une compr\u00e9hension approfondie des fondements math\u00e9matiques qui sous-tendent chaque algorithme d&rsquo;identification, afin d&rsquo;\u00e9viter l&rsquo;\u00e9cueil d&rsquo;une \u00e9puration aveugle des bases de donn\u00e9es.<\/p>\n<p>Le pr\u00e9sent trait\u00e9 m\u00e9thodologique a pour vocation de fournir un guide exhaustif, th\u00e9orique et op\u00e9rationnel d\u00e9di\u00e9 \u00e0 l&rsquo;identification rigoureuse des valeurs aberrantes univari\u00e9es et multivari\u00e9es sous SAS. En articulant les d\u00e9monstrations formelles des cl\u00f4tures de John Tukey avec des impl\u00e9mentations de code reproductibles pas \u00e0 pas sur un jeu de donn\u00e9es empirique \u00e9talonn\u00e9, ce document explore l&rsquo;ensemble du cycle de vie du diagnostic statistique : du partitionnement quantile \u00e0 l&rsquo;automatisation par macro-programmation, en passant par la mod\u00e9lisation robuste et les arbitrages d\u00e9ontologiques post-d\u00e9tection conformes aux normes internationales de publication scientifique.<\/p>\n<h2>1. Cadre th\u00e9orique et d\u00e9finition des valeurs aberrantes en analyse statistique<\/h2>\n<h3>1.1 Typologie et nature des valeurs atypiques dans les donn\u00e9es empiriques<\/h3>\n<p>Dans la litt\u00e9rature statistique fondamentale, telle qu&rsquo;\u00e9tablie par des auteurs pr\u00e9curseurs comme Vic Barnett et Toby Lewis, une valeur aberrante se d\u00e9finit formellement comme une observation qui semble diverger de mani\u00e8re marqu\u00e9e du sch\u00e9ma structurel caract\u00e9risant les autres \u00e9l\u00e9ments de la distribution au sein de laquelle elle s&rsquo;inscrit. Cependant, l&rsquo;\u00e9tiquette d&rsquo;\u00ab observation aberrante \u00bb recouvre des r\u00e9alit\u00e9s ontologiques et g\u00e9n\u00e9ratives profond\u00e9ment divergentes qu&rsquo;il convient de discriminer avec la plus haute pr\u00e9cision conceptuelle avant d&rsquo;envisager la moindre intervention algorithmique.<\/p>\n<p>On distingue traditionnellement trois cat\u00e9gories m\u00e8res d&rsquo;observations discordantes. La premi\u00e8re rel\u00e8ve des <strong>erreurs de mesure et d&rsquo;encodage proc\u00e9durales<\/strong>. Celles-ci regroupent les fautes de frappe humaine (comme l&rsquo;adjonction involontaire d&rsquo;un z\u00e9ro terminal transformant un score de 22 en 220), les d\u00e9faillances mat\u00e9rielles de capteurs \u00e9lectroniques, les ruptures de calibration d&rsquo;appareillages de laboratoire ou les erreurs de fusion de tables relationnelles lors des \u00e9tapes d&rsquo;extraction de donn\u00e9es. Ces valeurs ne poss\u00e8dent aucune validit\u00e9 \u00e9cologique et doivent \u00eatre imp\u00e9rativement purg\u00e9es ou corrig\u00e9es.<\/p>\n<p>La deuxi\u00e8me cat\u00e9gorie englobe les manifestations extr\u00eames de la <strong>variabilit\u00e9 biologique ou ph\u00e9nom\u00e9nologique intrins\u00e8que<\/strong>. Dans ce cas de figure, l&rsquo;observation r\u00e9sulte fid\u00e8lement d&rsquo;un processus g\u00e9n\u00e9rateur valide, mais refl\u00e8te la queue lourde d&rsquo;une distribution naturelle ou l&rsquo;occurrence d&rsquo;un \u00e9v\u00e9nement rare mais parfaitement plausible. Dans les \u00e9tudes \u00e9pid\u00e9miologiques ou psychom\u00e9triques, un sujet peut pr\u00e9senter un temps de r\u00e9action exceptionnellement long ou un profil immunologique singulier sans que la mesure ne soit intrins\u00e8quement vici\u00e9e. Traiter ces observations comme de simples anomalies techniques reviendrait \u00e0 appauvrir artificiellement la variabilit\u00e9 du ph\u00e9nom\u00e8ne mod\u00e9lis\u00e9.<\/p>\n<p>La troisi\u00e8me cat\u00e9gorie correspond aux <strong>artefacts de contamination d&rsquo;\u00e9chantillonnage<\/strong>, o\u00f9 l&rsquo;individu mesur\u00e9 n&rsquo;appartient pas, en r\u00e9alit\u00e9, \u00e0 la population cible postul\u00e9e par le chercheur. Par exemple, l&rsquo;inclusion accidentelle d&rsquo;un patient souffrant d&rsquo;une pathologie comorbide s\u00e9v\u00e8re non diagnostiqu\u00e9e dans un essai clinique ax\u00e9 sur des sujets sains g\u00e9n\u00e8re une observation atypique issue d&rsquo;un mod\u00e8le de m\u00e9lange sous-jacent.<\/p>\n<p>D&rsquo;un point de vue structural, les observations discordantes se classifient selon leur dimensionalit\u00e9 g\u00e9om\u00e9trique :<\/p>\n<ul>\n<li><strong>Anomalies univari\u00e9es :<\/strong> Elles se manifestent par un \u00e9cart substantiel le long de l&rsquo;axe r\u00e9el d&rsquo;une variable quantitative isol\u00e9e. Leur rep\u00e9rage s&rsquo;appuie sur la dispersion marginale de cette seule variable, ind\u00e9pendamment de toute corr\u00e9lation externe.<\/li>\n<li><strong>Anomalies bivari\u00e9es et multivari\u00e9es :<\/strong> Elles correspondent \u00e0 des configurations vectorielles atypiques o\u00f9 chaque variable prise isol\u00e9ment pr\u00e9sente une valeur parfaitement plausible au sein de sa distribution marginale, mais dont la combinaison conjointe viole la structure de covariance ou de corr\u00e9lation de l&rsquo;\u00e9chantillon multidimensionnel.<\/li>\n<\/ul>\n<p>L&rsquo;asym\u00e9trie distributionnelle (mesur\u00e9e par le coefficient de <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/Skewness.html\"><em>skewness<\/em><\/a> de Fisher-Pearson) exerce une influence trompeuse majeure sur l&rsquo;\u00e9mergence apparente d&rsquo;anomalies. Dans une distribution fortement asym\u00e9trique \u00e0 droite &mdash; telle que celle des revenus financiers, des temps d&rsquo;acc\u00e8s serveur ou de la concentration plasmatique de biomarqueurs &mdash;, les valeurs situ\u00e9es dans la queue sup\u00e9rieure peuvent \u00eatre \u00e9tiquet\u00e9es \u00e0 tort comme des aberrations par les m\u00e9triques classiques alors qu&rsquo;elles sont strictement conformes \u00e0 la loi th\u00e9orique sous-jacente (par exemple, une loi log-normale ou de Pareto). Les cons\u00e9quences directes sur les estimateurs classiques sont d\u00e9vastatrices : la moyenne arithm\u00e9tique et la variance d&rsquo;\u00e9chantillon, \u00e9tant d\u00e9pourvues de toute r\u00e9sistance, sont litt\u00e9ralement tract\u00e9es vers les extr\u00eames, gonflant artificiellement les intervalles de confiance et alt\u00e9rant l&rsquo;ensemble des diagnostics param\u00e9triques d\u00e9riv\u00e9s.<\/p>\n<h3>1.2 Enjeux \u00e9pist\u00e9mologiques et psychom\u00e9triques de la d\u00e9tection d&rsquo;observations extr\u00eames<\/h3>\n<p>La d\u00e9tection et le traitement des observations extr\u00eames ne rel\u00e8vent pas uniquement d&rsquo;une m\u00e9canique algorithmique ; ils soul\u00e8vent des dilemmes \u00e9pist\u00e9mologiques critiques quant \u00e0 la nature m\u00eame de l&rsquo;induction scientifique. Dans le champ des sciences comportementales, de la psychom\u00e9trie et de la recherche biom\u00e9dicale, l&rsquo;exclusion arbitraire ou non document\u00e9e de donn\u00e9es atypiques constitue l&rsquo;un des vecteurs majeurs de la crise de la reproductibilit\u00e9 scientifique. Le chercheur se trouve confront\u00e9 \u00e0 une tension dialectique permanente entre la pr\u00e9servation de la <em>validit\u00e9 \u00e9cologique<\/em> &mdash; le respect de la complexit\u00e9 intrins\u00e8que et de la variabilit\u00e9 du r\u00e9el &mdash; et la sauvegarde de l&rsquo;<em>int\u00e9grit\u00e9 statistique<\/em> des mod\u00e9lisations lin\u00e9aires g\u00e9n\u00e9rales.<\/p>\n<p>Au sein des mod\u00e8les lin\u00e9aires g\u00e9n\u00e9raux (r\u00e9gression multiple, analyse de variance factorielle, mod\u00e9lisation par \u00e9quations structurelles), la pr\u00e9sence de valeurs aberrantes non contr\u00f4l\u00e9es viole les postulats cardinaux de normalit\u00e9 des r\u00e9sidus et d&rsquo;homosc\u00e9dasticit\u00e9. Ces violations induisent des biais s\u00e9v\u00e8res dans l&rsquo;estimation des coefficients de r\u00e9gression par les moindres carr\u00e9s ordinaires (MCO), dont le point de rupture (<em>breakdown point<\/em>) th\u00e9orique est \u00e9gal \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"1\/n\">$1\/n$<\/span>, ce qui signifie qu&rsquo;une seule observation infinit\u00e9simale ou gigantesque suffit \u00e0 faire diverger arbitrairement la droite d&rsquo;ajustement. De surcro\u00eet, les erreurs-types des param\u00e8tres estim\u00e9s se trouvent artificiellement gonfl\u00e9es ou comprim\u00e9es selon la position g\u00e9om\u00e9trique de l&rsquo;anomalie, provoquant une distorsion impr\u00e9visible des taux d&rsquo;erreur de type I (rejet \u00e0 tort de l&rsquo;hypoth\u00e8se nulle d&rsquo;absence d&rsquo;effet) et de type II (incapacit\u00e9 \u00e0 d\u00e9tecter un effet r\u00e9el existant).<\/p>\n<p>Face \u00e0 ces risques, les standards internationaux de publication scientifique &mdash; notamment codifi\u00e9s par l&rsquo;<em>American Psychological Association<\/em> (normes APA 7e \u00e9dition) et le r\u00e9seau EQUATOR via les lignes directrices CONSORT et STROBE &mdash; imposent une transparence m\u00e9thodologique absolue. Il est d\u00e9sormais formellement proscrit d&rsquo;\u00e9liminer des observations en coulisses sans rapporter explicitement les crit\u00e8res math\u00e9matiques pr\u00e9alablement fix\u00e9s ayant guid\u00e9 l&rsquo;exclusion, le nombre exact d&rsquo;observations \u00e9cart\u00e9es, leur profil descriptif et une analyse de sensibilit\u00e9 d\u00e9montrant la stabilit\u00e9 des conclusions post-exclusion.<\/p>\n<h3>1.3 Panorama comparatif des approches de d\u00e9tection univari\u00e9es<\/h3>\n<p>L&rsquo;arsenal statistique d\u00e9di\u00e9 au diagnostic univari\u00e9 des valeurs atypiques se scinde en deux paradigmes fondamentaux : les approches param\u00e9triques fond\u00e9es sur l&rsquo;hypoth\u00e8se de normalit\u00e9 et les m\u00e9thodes non param\u00e9triques ax\u00e9es sur la position relative et les quantiles d&rsquo;\u00e9chantillonnage.<\/p>\n<p>L&rsquo;approche param\u00e9trique la plus intuitive repose sur la transformation en <strong>scores standardis\u00e9s (scores Z)<\/strong>, o\u00f9 chaque observation <span class=\"arabpsy-math-inline\" data-latex=\"x_i\">$x_i$<\/span> est recentr\u00e9e par rapport \u00e0 la moyenne arithm\u00e9tique <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{x}\">$\\bar{x}$<\/span> et r\u00e9duite par l&rsquo;\u00e9cart-type d&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Z_i = \\frac{x_i - \\bar{x}}{s}\">$$Z_i = \\frac{x_i &#8211; \\bar{x}}{s}$$<\/span><\/p>\n<p>La r\u00e8gle empirique conventionnelle stipule qu&rsquo;une observation dont la valeur absolue du score Z exc\u00e8de le seuil critique de 2,5 ou 3,0 doit \u00eatre consid\u00e9r\u00e9e comme potentiellement anormale, en vertu des propri\u00e9t\u00e9s de la loi normale standard o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"99,73%\">$99,73%$<\/span> de la masse de probabilit\u00e9 est confin\u00e9e dans l&rsquo;intervalle <span class=\"arabpsy-math-inline\" data-latex=\"[-3\\sigma, +3\\sigma]\">$[-3\\sigma, +3\\sigma]$<\/span>. N\u00e9anmoins, cette m\u00e9thode souffre d&rsquo;une aporie fondamentale connue sous le nom d&rsquo;<strong>effet de masque<\/strong> (<em>masking effect<\/em>). La moyenne arithm\u00e9tique et l&rsquo;\u00e9cart-type \u00e9tant eux-m\u00eames extr\u00eamement sensibles aux valeurs aberrantes, la pr\u00e9sence de plusieurs anomalies de grande magnitude gonfle d\u00e9mesur\u00e9ment l&rsquo;\u00e9cart-type <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span> et attire <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{x}\">$\\bar{x}$<\/span> vers l&rsquo;ext\u00e9rieur. Par cons\u00e9quent, les d\u00e9nominateurs des scores Z deviennent anormalement \u00e9lev\u00e9s, ce qui compresse artificiellement la valeur calcul\u00e9e de <span class=\"arabpsy-math-inline\" data-latex=\"Z_i\">$Z_i$<\/span> pour des anomalies secondaires qui \u00e9chappent alors totalement \u00e0 la d\u00e9tection.<\/p>\n<p>\u00c0 l&rsquo;inverse, les m\u00e9thodologies non param\u00e9triques fond\u00e9es sur les quantiles &mdash; au premier rang desquelles figure la m\u00e9thode de l&rsquo;\u00e9cart interquartile conceptualis\u00e9e par John Tukey &mdash; s&rsquo;affranchissent int\u00e9gralement de toute hypoth\u00e8se distributionnelle gaussienne. Les quantiles d&rsquo;\u00e9chantillonnage poss\u00e8dent une robustesse intrins\u00e8que remarquable, leur point de rupture atteignant 25 % pour les quartiles inf\u00e9rieur et sup\u00e9rieur. Le tableau analytique suivant synth\u00e9tise les propri\u00e9t\u00e9s op\u00e9rationnelles de ces grandes familles de d\u00e9tection univari\u00e9e :<\/p>\n<ul>\n<li><strong>R\u00e8gle empirique des scores Z :<\/strong> N\u00e9cessite une stricte normalit\u00e9 sous-jacente ; point de rupture th\u00e9orique nul (0 %) ; hautement vuln\u00e9rable au masquage et \u00e0 l&rsquo;asym\u00e9trie ; applicable pr\u00e9f\u00e9rentiellement aux tr\u00e8s grands \u00e9chantillons strictement gaussiens.<\/li>\n<li><strong>\u00c9cart Absolu M\u00e9dian (MAD) :<\/strong> Mesure de dispersion hautement robuste centr\u00e9e sur la m\u00e9diane ; point de rupture optimal de 50 % ; immunis\u00e9e contre le masquage ; particuli\u00e8rement adapt\u00e9e aux distributions unimodales contamin\u00e9es.<\/li>\n<li><strong>\u00c9cart Interquartile (IQR \/ Tukey) :<\/strong> Absence d&rsquo;hypoth\u00e8se de normalit\u00e9 ; point de rupture de 25 % ; excellente r\u00e9sistance aux queues lourdes ; id\u00e9ale pour les jeux de donn\u00e9es de taille petite \u00e0 mod\u00e9r\u00e9e dans les cadres de recherche exploratoire et appliqu\u00e9e.<\/li>\n<\/ul>\n<h2>2. Fondements math\u00e9matiques de la m\u00e9thode de l&rsquo;\u00e9cart interquartile (IQR)<\/h2>\n<h3>2.1 D\u00e9finition analytique des quartiles et de la dispersion m\u00e9diane<\/h3>\n<p>L&rsquo;analyse exploratoire de donn\u00e9es (<em>Exploratory Data Analysis<\/em>, EDA), th\u00e9oris\u00e9e par <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/John_Tukey\">John W. Tukey<\/a> en 1977, privil\u00e9gie syst\u00e9matiquement les statistiques d&rsquo;ordre (<em>order statistics<\/em>) par rapport aux moments d&rsquo;\u00e9chantillonnage pour caract\u00e9riser la tendance centrale et l&rsquo;\u00e9talement des distributions empiriques. Les quartiles constituent les trois points de coupure cardinaux divisant un ensemble d&rsquo;observations tri\u00e9es par ordre croissant en quatre sous-groupes d&rsquo;effectifs rigoureusement \u00e9quiprobables.<\/p>\n<p>Soit un \u00e9chantillon ordonn\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"X = (x_{(1)}, x_{(2)}, dots, x_{(n)})\">$X = (x_{(1)}, x_{(2)}, dots, x_{(n)})$<\/span> issu d&rsquo;une variable quantitative continue, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"x_{(1)} le x_{(2)} le dots le x_{(n)}\">$x_{(1)} le x_{(2)} le dots le x_{(n)}$<\/span>. Le premier quartile, not\u00e9 formalement <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"P_{25}\">$P_{25}$<\/span>, correspond au 25e centile empirique : il repr\u00e9sente la valeur seuil en de\u00e7\u00e0 de laquelle se situe exactement un quart (25 %) des donn\u00e9es observ\u00e9es. Le troisi\u00e8me quartile, not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"P_{75}\">$P_{75}$<\/span>, d\u00e9signe le 75e centile empirique, marquant le seuil en dessous duquel sont contenues 75 % des observations, ou de fa\u00e7on sym\u00e9trique, au-del\u00e0 duquel ne subsiste que le quart sup\u00e9rieur de la distribution. Le deuxi\u00e8me quartile (<span class=\"arabpsy-math-inline\" data-latex=\"Q_2\">$Q_2$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"P_{50}\">$P_{50}$<\/span>) correspond rigoureusement \u00e0 la m\u00e9diane empirique de l&rsquo;\u00e9chantillon.<\/p>\n<p>L&rsquo;<strong>\u00e9cart interquartile<\/strong>, universellement d\u00e9sign\u00e9 par l&rsquo;acronyme <strong>IQR<\/strong> (pour <em>Interquartile Range<\/em>), constitue la mesure non param\u00e9trique canonique de la dispersion statistique. Il se d\u00e9finit analytiquement par la diff\u00e9rence lin\u00e9aire directe entre le troisi\u00e8me et le premier quartile :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{IQR} = Q_3 - Q_1\">$$\\text{IQR} = Q_3 &#8211; Q_1$$<\/span><\/p>\n<p>La grandeur <span class=\"arabpsy-math-inline\" data-latex=\"\\text{IQR}\">$\\text{IQR}$<\/span> quantifie avec une pr\u00e9cision robuste l&rsquo;\u00e9tendue spatiale couverte par les 50 % centraux de la population observ\u00e9e. Contrairement \u00e0 l&rsquo;\u00e9tendue globale de l&rsquo;\u00e9chantillon (<span class=\"arabpsy-math-inline\" data-latex=\"R = x_{(n)} - x_{(1)}\">$R = x_{(n)} &#8211; x_{(1)}$<\/span>) qui est enti\u00e8rement asservie aux deux observations les plus extr\u00eames, l&rsquo;\u00e9cart interquartile poss\u00e8de une propri\u00e9t\u00e9 de r\u00e9sistance structurelle absolue face aux perturbations marginales. Toute alt\u00e9ration, quelle qu&rsquo;en soit la d\u00e9mesure math\u00e9matique, apport\u00e9e aux valeurs contenues dans le premier quart inf\u00e9rieur ($x <Q> Q_3$) laisse l&rsquo;\u00e9cart interquartile strictement invariant, pourvu que les rangs relatifs des observations ne basculent pas de part et d&rsquo;autre de ces fronti\u00e8res fractiles.<\/p>\n<h3>2.2 Formalisation des cl\u00f4tures de Tukey (Tukey&rsquo;s Fences)<\/h3>\n<p>Pour operationaliser la d\u00e9marcation objective entre la variabilit\u00e9 attendue d&rsquo;une distribution et les observations manifestant une d\u00e9viance anormale, John Tukey a d\u00e9velopp\u00e9 un dispositif algorithmique g\u00e9om\u00e9trique reposant sur la construction de seuils d&rsquo;exclusion ext\u00e9rieurs, conventionnellement d\u00e9sign\u00e9s sous le terme de <strong>cl\u00f4tures de Tukey<\/strong> (<em>Tukey&rsquo;s Fences<\/em>). Ces barri\u00e8res statistiques sont g\u00e9n\u00e9r\u00e9es en projetant une distance proportionnelle \u00e0 l&rsquo;\u00e9cart interquartile de part et d&rsquo;autre du c\u0153ur de la distribution mat\u00e9rialis\u00e9 par <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span>.<\/p>\n<p>Analytiquement, les cl\u00f4tures int\u00e9rieure inf\u00e9rieure (<span class=\"arabpsy-math-inline\" data-latex=\"B_{\\text{\\inf}}\">$B_{\\text{\\inf}}$<\/span>) et int\u00e9rieure sup\u00e9rieure (<span class=\"arabpsy-math-inline\" data-latex=\"B_{\\text{\\sup}}\">$B_{\\text{\\sup}}$<\/span>) sont d\u00e9finies par les \u00e9quations fondamentales suivantes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"B_{\\text{\\inf}} = Q_1 - (k \\times \\text{IQR})\">$$B_{\\text{\\inf}} = Q_1 &#8211; (k \\times \\text{IQR})$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"B_{\\text{\\sup}} = Q_3 + (k \\times \\text{IQR})\">$$B_{\\text{\\sup}} = Q_3 + (k \\times \\text{IQR})$$<\/span><\/p>\n<p>Dans la formalisation canonique \u00e9tablie par Tukey pour la d\u00e9tection standard des valeurs aberrantes mod\u00e9r\u00e9es, le coefficient multiplicateur <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> est fix\u00e9 rigoureusement \u00e0 la constante scalaire $1,5$ :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Cl\u00f4ture Inf\u00e9rieure} = Q_1 - 1,5 \\times \\text{IQR}\">$$\\text{Cl\u00f4ture Inf\u00e9rieure} = Q_1 &#8211; 1,5 \\times \\text{IQR}$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Cl\u00f4ture Sup\u00e9rieure} = Q_3 + 1,5 \\times \\text{IQR}\">$$\\text{Cl\u00f4ture Sup\u00e9rieure} = Q_3 + 1,5 \\times \\text{IQR}$$<\/span><\/p>\n<p>La justification probabiliste de la s\u00e9lection du coefficient multiplicateur <span class=\"arabpsy-math-inline\" data-latex=\"k = 1,5\">$k = 1,5$<\/span> s&rsquo;enracine dans le comportement th\u00e9orique d&rsquo;une distribution gaussienne id\u00e9ale <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{N}(\\mu, \\sigma^2)\">$\\mathcal{N}(\\mu, \\sigma^2)$<\/span>. Sous une loi normale standardis\u00e9e, les quartiles th\u00e9oriques se d\u00e9duisent analytiquement de la fonction de r\u00e9partition inverse <span class=\"arabpsy-math-inline\" data-latex=\"\\Phi^{-1}\">$\\Phi^{-1}$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Q_1 = \\mu - 0,6745\\sigma \\quad \\text{et} \\quad Q_3 = \\mu + 0,6745\\sigma\">$$Q_1 = \\mu &#8211; 0,6745\\sigma \\quad \\text{et} \\quad Q_3 = \\mu + 0,6745\\sigma$$<\/span><\/p>\n<p>Il en d\u00e9coule que l&rsquo;\u00e9cart interquartile th\u00e9orique d&rsquo;une loi normale vaut approximativement :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{IQR} = 1,349\\sigma\">$$\\text{IQR} = 1,349\\sigma$$<\/span><\/p>\n<p>En substituant ces valeurs dans l&rsquo;expression des cl\u00f4tures de Tukey avec <span class=\"arabpsy-math-inline\" data-latex=\"k = 1,5\">$k = 1,5$<\/span>, les bornes critiques se d\u00e9ploient spatialement \u00e0 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Bornes} = \\mu \\pm (0,6745\\sigma + 1,5 \\times 1,349\\sigma) \\approx \\mu \\pm 2,698\\sigma\">$$\\text{Bornes} = \\mu \\pm (0,6745\\sigma + 1,5 \\times 1,349\\sigma) \\approx \\mu \\pm 2,698\\sigma$$<\/span><\/p>\n<p>Sur une variable al\u00e9atoire parfaitement normale, la probabilit\u00e9 th\u00e9orique pour qu&rsquo;une observation valide franchisse l&rsquo;une ou l&rsquo;autre de ces cl\u00f4tures est d&rsquo;exactement <span class=\"arabpsy-math-inline\" data-latex=\"2 \\times \\Phi(-2,698) \\approx 0,00698\">$2 \\times \\Phi(-2,698) \\approx 0,00698$<\/span>, soit un taux d&rsquo;\u00e9tiquetage d&rsquo;environ <strong>0,7 %<\/strong> (environ 7 observations sur 1000). Ce calibrage offre un \u00e9quilibre remarquable entre sensibilit\u00e9 diagnostique et protection contre le faux \u00e9tiquetage.<\/p>\n<p>Pour discriminer les observations d\u00e9viantes ordinaires des anomalies catastrophiques, Tukey a \u00e9galement formalis\u00e9 les <strong>cl\u00f4tures ext\u00e9rieures<\/strong> (<em>outer fences<\/em>) en \u00e9levant le multiplicateur \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"k = 3,0\">$k = 3,0$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Cl\u00f4ture Ext\u00e9rieure Inf\u00e9rieure} = Q_1 - 3,0 \\times \\text{IQR}\">$$\\text{Cl\u00f4ture Ext\u00e9rieure Inf\u00e9rieure} = Q_1 &#8211; 3,0 \\times \\text{IQR}$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Cl\u00f4ture Ext\u00e9rieure Sup\u00e9rieure} = Q_3 + 3,0 \\times \\text{IQR}\">$$\\text{Cl\u00f4ture Ext\u00e9rieure Sup\u00e9rieure} = Q_3 + 3,0 \\times \\text{IQR}$$<\/span><\/p>\n<p>Une observation situ\u00e9e entre la cl\u00f4ture int\u00e9rieure (1,5) et la cl\u00f4ture ext\u00e9rieure (3,0) est qualifi\u00e9e de <em>valeur aberrante mod\u00e9r\u00e9e<\/em> (ou valeur suspecte), tandis qu&rsquo;une observation exc\u00e9dant la cl\u00f4ture ext\u00e9rieure de 3,0 IQR est classifi\u00e9e comme <em>valeur aberrante extr\u00eame<\/em>. Dans le cadre normal th\u00e9orique, une valeur extr\u00eame (<span class=\"arabpsy-math-inline\" data-latex=\"k = 3,0\">$k = 3,0$<\/span>) correspond \u00e0 un \u00e9cartement de <span class=\"arabpsy-math-inline\" data-latex=\"\\pm 4,72\\sigma\">$\\pm 4,72\\sigma$<\/span>, ce qui correspond \u00e0 une probabilit\u00e9 th\u00e9orique infime de l&rsquo;ordre de <span class=\"arabpsy-math-inline\" data-latex=\"2,3 \\times 10^{-6}\">$2,3 \\times 10^{-6}$<\/span> (deux observations par million).<\/p>\n<h3>2.3 Robustesse et propri\u00e9t\u00e9s asymptotiques sous distributions asym\u00e9triques<\/h3>\n<p>Si les cl\u00f4tures de Tukey avec un coefficient <span class=\"arabpsy-math-inline\" data-latex=\"k = 1,5\">$k = 1,5$<\/span> pr\u00e9sentent des vertus heuristiques ind\u00e9niables, leur comportement asymptotique face \u00e0 des distributions empiriques hautement asym\u00e9triques &mdash; telles que les lois exponentielles, log-normales, Gamma ou de Weibull &mdash; expose l&rsquo;analyste \u00e0 un risque s\u00e9v\u00e8re de <strong>sur-d\u00e9tection unilat\u00e9rale<\/strong>.<\/p>\n<p>Lorsque la distribution pr\u00e9sente une asym\u00e9trie positive marqu\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"Skewness &gt; 0\">$Skewness &gt; 0$<\/span>), la distance physique s\u00e9parant la m\u00e9diane de <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span> s&rsquo;av\u00e8re nettement sup\u00e9rieure \u00e0 celle s\u00e9parant <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span> de la m\u00e9diane. Bien que l&rsquo;IQR capte une mesure d&rsquo;\u00e9talement global, l&rsquo;application d&rsquo;un intervalle strictement sym\u00e9trique (<span class=\"arabpsy-math-inline\" data-latex=\"1,5 \\times \\text{IQR}\">$1,5 \\times \\text{IQR}$<\/span>) au-del\u00e0 de <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span> ne compense pas l&rsquo;allongement naturel de la queue droite de la loi. En cons\u00e9quence math\u00e9matique directe, le taux empirique d&rsquo;observations l\u00e9gitimes exc\u00e9dant la cl\u00f4ture sup\u00e9rieure peut s&rsquo;\u00e9lever jusqu&rsquo;\u00e0 5 %, 10 % voire 15 % dans des \u00e9chantillons log-normaux \u00e0 forte variance r\u00e9siduelle, conduisant l&rsquo;analyste \u00e0 d\u00e9truire des fractions massives de donn\u00e9es parfaitement valides.<\/p>\n<p>Pour pallier cette distorsion structurelle sans renoncer \u00e0 la philosophie des quantiles, la statistique computationnelle moderne propose l&rsquo;utilisation de l&rsquo;<strong>\u00e9cart interquartile ajust\u00e9 par la m\u00e9dcouple<\/strong> (not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"MC\">$MC$<\/span>). D\u00e9velopp\u00e9e par Mia Hubert et Ellen Vandervieren en 2008, cette m\u00e9thodologie int\u00e8gre un coefficient de sym\u00e9trie robuste born\u00e9 dans l&rsquo;intervalle <span class=\"arabpsy-math-inline\" data-latex=\"[-1, 1]\">$[-1, 1]$<\/span>. Les cl\u00f4tures asym\u00e9triques sont alors red\u00e9finies dynamiquement par des fonctions exponentielles de <span class=\"arabpsy-math-inline\" data-latex=\"MC\">$MC$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"B_{\\text{\\sup}} = Q_3 + 1,5 \\times e^{3 \\times MC} \\times \\text{IQR} \\quad \\text{si } MC ge 0\">$$B_{\\text{\\sup}} = Q_3 + 1,5 \\times e^{3 \\times MC} \\times \\text{IQR} \\quad \\text{si } MC ge 0$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"B_{\\text{\\inf}} = Q_1 - 1,5 \\times e^{-4 \\times MC} \\times \\text{IQR} \\quad \\text{si } MC ge 0\">$$B_{\\text{\\inf}} = Q_1 &#8211; 1,5 \\times e^{-4 \\times MC} \\times \\text{IQR} \\quad \\text{si } MC ge 0$$<\/span><\/p>\n<p>Lorsque <span class=\"arabpsy-math-inline\" data-latex=\"MC = 0\">$MC = 0$<\/span>, la formule se r\u00e9duit rigoureusement aux cl\u00f4tures classiques de Tukey. N\u00e9anmoins, pour des variables continues issues de protocoles psychologiques standardis\u00e9s o\u00f9 l&rsquo;asym\u00e9trie demeure mod\u00e9r\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"|Skewness| &lt; 1,0\">$|Skewness| &lt; 1,0$<\/span>), la r\u00e8gle empirique stricte de Tukey (<span class=\"arabpsy-math-inline\" data-latex=\"1,5 \\times \\text{IQR}\">$1,5 \\times \\text{IQR}$<\/span>) demeure un \u00e9talon analytique fiable, universellement accept\u00e9 et directement programmable sous les architectures SAS.<\/p>\n<h2>3. Configuration de l&rsquo;environnement SAS et g\u00e9n\u00e9ration du jeu de donn\u00e9es<\/h2>\n<h3>3.1 Initialisation de l&rsquo;\u00e9tape DATA et saisie des observations<\/h3>\n<p>La mise en \u0153uvre pratique de l&rsquo;identification des observations atypiques sous SAS requiert la ma\u00eetrise fondamentale de l&rsquo;\u00e9tape <code>DATA<\/code>. Cette composante programmatique primaire permet d&rsquo;initialiser l&rsquo;espace m\u00e9moire vectoriel, de structurer les descripteurs de m\u00e9tadonn\u00e9es et d&rsquo;injecter s\u00e9quentiellement les flux de donn\u00e9es brutes via un tampon d&rsquo;entr\u00e9e (<em>Input Buffer<\/em>) vers le vecteur d&rsquo;observation du programme (<em>Program Data Vector<\/em>, PDV).<\/p>\n<p>Afin d&rsquo;illustrer la m\u00e9canique algorithmique de d\u00e9tection sur un cas didactique reproductible et scientifiquement explicite, nous construisons une table SAS nomm\u00e9e <code>DONNEES_ORIGINALES<\/code>. Ce jeu de donn\u00e9es synth\u00e9tise les performances d&rsquo;\u00e9quipes de recherche ou de collectifs sportifs soumises \u00e0 une m\u00e9trique d&rsquo;\u00e9valuation quantitative continue &mdash; la variable <code>Points<\/code> &mdash; index\u00e9e par une variable nominale identifiante nomm\u00e9e <code>Equipe<\/code>.<\/p>\n<p>L&rsquo;impl\u00e9mentation de l&rsquo;\u00e9tape <code>DATA<\/code> s&rsquo;ex\u00e9cute selon le code SAS certifi\u00e9 suivant :<\/p>\n<p><!-- Syntax-clean SAS Step Data Block --><\/p>\n<p><code>DATA DONNEES_ORIGINALES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;LENGTH Equipe $ 1;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;INPUT Equipe $ Points;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DATALINES;<br \/>\nA 18<br \/>\nB 24<br \/>\nC 26<br \/>\nD 34<br \/>\nE 38<br \/>\nF 45<br \/>\nG 48<br \/>\nH 54<br \/>\nI 60<br \/>\nJ 73<br \/>\nK 79<br \/>\nL 85<br \/>\nM 94<br \/>\nN 98<br \/>\nO 221<br \/>\nP 223<br \/>\n;<br \/>\nRUN;<\/code><\/p>\n<p>La ligne <code>LENGTH Equipe $ 1;<\/code> dimensionne explicitement la variable textuelle \u00e0 un seul octet alphanum\u00e9rique, optimisant ainsi l&rsquo;allocation des pages m\u00e9moire du moteur de stockage SAS (V9 engine). L&rsquo;instruction <code>INPUT Equipe $ Points;<\/code> informe le compilateur que la premi\u00e8re colonne doit \u00eatre interpr\u00e9t\u00e9e comme un scalaire textuel (indiqu\u00e9 par le modificateur mon\u00e9taire <code>$<\/code>), tandis que la deuxi\u00e8me colonne <code>Points<\/code> est compil\u00e9e sous le format num\u00e9rique standard \u00e0 virgule flottante double pr\u00e9cision (8 octets selon la norme IEEE 754). Le d\u00e9limiteur <code>DATALINES;<\/code> ouvre le segment d&rsquo;ingestion s\u00e9quentielle des 16 observations empiriques, cl\u00f4tur\u00e9 formellement par le point-virgule isol\u00e9.<\/p>\n<h3>3.2 Inspection pr\u00e9liminaire du jeu de donn\u00e9es via PROC PRINT<\/h3>\n<p>D\u00e8s l&rsquo;ex\u00e9cution de l&rsquo;\u00e9tape <code>DATA<\/code>, l&rsquo;analyste rigoureux doit proc\u00e9der \u00e0 la validation de la table physique g\u00e9n\u00e9r\u00e9e. Le journal d&rsquo;ex\u00e9cution de SAS (<em>SAS Log<\/em>) doit faire l&rsquo;objet d&rsquo;un examen minutieux : il convient de v\u00e9rifier l&rsquo;absence absolue de messages de type <code>NOTE: Invalid data for Points in line...<\/code> ou d&rsquo;avertissements (<code>WARNING<\/code>) signalant des conversions implicites de variables ou des troncatures de cha\u00eenes.<\/p>\n<p>L&rsquo;instruction <code>PROC PRINT<\/code> permet d&rsquo;\u00e9diter une vue tabulaire exhaustive du contenu du vecteur m\u00e9moire pour s&rsquo;assurer de la parfaite concordance s\u00e9quentielle des observations et valider l&rsquo;absence de valeurs manquantes (<em>missing values<\/em> symbolis\u00e9es par un point en SAS num\u00e9rique) :<\/p>\n<p><code>PROC PRINT DATA=DONNEES_ORIGINALES NOOBS;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"Examen pr\u00e9liminaire du jeu de donn\u00e9es exp\u00e9rimentales\";<br \/>\nRUN;<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/outlier1.jpg\" alt=\"\" \/><\/figure>\n<p>La consultation conjointe de la proc\u00e9dure <code>PROC CONTENTS<\/code> assure la conformit\u00e9 du dictionnaire de donn\u00e9es et des attributs de table :<\/p>\n<p><code>PROC CONTENTS DATA=DONNEES_ORIGINALES;<br \/>\nRUN;<\/code><\/p>\n<p>Cette v\u00e9rification confirme la pr\u00e9sence d&rsquo;exactement 16 observations r\u00e9parties sur deux variables statistiques distinctes, sans compression de blocs, scellant ainsi l&rsquo;int\u00e9grit\u00e9 de la table avant toute manipulation calculatoire.<\/p>\n<h3>3.3 Structuration formelle de l&rsquo;exemple didactique<\/h3>\n<p>La structure num\u00e9rique d\u00e9lib\u00e9r\u00e9ment instill\u00e9e dans cet \u00e9chantillon de 16 modalit\u00e9s (index\u00e9es de <code>A<\/code> \u00e0 <code>P<\/code>) permet d&rsquo;observer l&rsquo;ensemble des dynamiques associ\u00e9es aux queues de distribution. Pour les quatorze premi\u00e8res observations (de l&rsquo;\u00e9quipe <code>A<\/code> avec 18 points \u00e0 l&rsquo;\u00e9quipe <code>N<\/code> avec 98 points), la variable quantitative progresse de fa\u00e7on mod\u00e9r\u00e9e et r\u00e9guli\u00e8re, simulant une distribution empirique standard dont l&rsquo;\u00e9talement s&rsquo;av\u00e8re continu.<\/p>\n<p>En revanche, les deux derni\u00e8res observations &mdash; l&rsquo;\u00e9quipe <code>O<\/code> (221 points) et l&rsquo;\u00e9quipe <code>P<\/code> (223 points) &mdash; introduisent une discontinuit\u00e9 quantitative disproportionn\u00e9e. Ce saut d&rsquo;amplitude num\u00e9rique, mat\u00e9rialis\u00e9 par plus d&rsquo;un doublement de valeur par rapport au score pr\u00e9c\u00e9dent (passant de 98 \u00e0 221), configure un sc\u00e9nario critique typique de donn\u00e9es de terrain pollu\u00e9es par des aberrations extr\u00eames. L&rsquo;objectif m\u00e9thodologique consiste \u00e0 d\u00e9terminer si les cl\u00f4tures de Tukey calcul\u00e9es de fa\u00e7on automatis\u00e9e sous SAS permettront d&rsquo;isoler objectivement les \u00e9quipes <code>O<\/code> et <code>P<\/code> sans rejeter abusivement les scores l\u00e9gitimement \u00e9lev\u00e9s de la s\u00e9rie continue (tels que les \u00e9quipes <code>M<\/code> et <code>N<\/code>).<\/p>\n<h2>4. Extraction computationnelle des quantiles avec PROC UNIVARIATE<\/h2>\n<h3>4.1 Param\u00e9trage avanc\u00e9 de la proc\u00e9dure PROC UNIVARIATE<\/h3>\n<p>La proc\u00e9dure <code>PROC UNIVARIATE<\/code> constitue le moteur analytique par excellence sous SAS pour l&rsquo;\u00e9valuation descriptive approfondie, l&rsquo;estimation des fractiles et l&rsquo;analyse de conformation distributionnelle. Contrairement \u00e0 la proc\u00e9dure <code>PROC MEANS<\/code> ou <code>PROC SUMMARY<\/code> &mdash; orient\u00e9es vers l&rsquo;agr\u00e9gation de masse rapide &mdash;, <code>PROC UNIVARIATE<\/code> d\u00e9ploie des routines algorithmiques avanc\u00e9es pour le tri et l&rsquo;interpolation pr\u00e9cise des quantiles.<\/p>\n<p>Une probl\u00e9matique technique fondamentale r\u00e9side dans l&rsquo;algorithme sous-jacent employ\u00e9 pour localiser et interpoler les quantiles lorsque la taille de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> ne permet pas une division enti\u00e8re stricte par 4. SAS r\u00e9sout cette contingence computationnelle en mettant \u00e0 disposition de l&rsquo;utilisateur l&rsquo;option syst\u00e8me <code>PCTLDEF=<\/code> (<em>Percentile Definition<\/em>), qui supporte cinq d\u00e9finitions math\u00e9matiques distinctes d&rsquo;interpolation :<\/p>\n<ul>\n<li><strong>PCTLDEF=1 :<\/strong> Moyenne pond\u00e9r\u00e9e au point empirique (m\u00e9thode standard des probabilit\u00e9s d&rsquo;ordre).<\/li>\n<li><strong>PCTLDEF=2 :<\/strong> Observation empirique la plus proche de la valeur <span class=\"arabpsy-math-inline\" data-latex=\"n \\times p\">$n \\times p$<\/span>.<\/li>\n<li><strong>PCTLDEF=3 :<\/strong> R\u00e8gle empirique de troncature sans interpolation (statistique d&rsquo;ordre brute).<\/li>\n<li><strong>PCTLDEF=4 :<\/strong> Interpolation lin\u00e9aire standardis\u00e9e entre les observations adjacentes (option par d\u00e9faut sous de nombreux tableurs).<\/li>\n<li><strong>PCTLDEF=5 :<\/strong> Interpolation lin\u00e9aire empirique bas\u00e9e sur la fonction de r\u00e9partition continue (d\u00e9finition canonique retenue par d\u00e9faut par la proc\u00e9dure <code>PROC UNIVARIATE<\/code> sous SAS).<\/li>\n<\/ul>\n<p>Sous la d\u00e9finition par d\u00e9faut <code>PCTLDEF=5<\/code>, la position du quantile d&rsquo;ordre <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> (avec <span class=\"arabpsy-math-inline\" data-latex=\"p = 0,25\">$p = 0,25$<\/span> pour <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"p = 0,75\">$p = 0,75$<\/span> pour <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span>) est calcul\u00e9e sur la base de la valeur indicatrice <span class=\"arabpsy-math-inline\" data-latex=\"j = n \\times p\">$j = n \\times p$<\/span>. Si <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span> est un entier, le quantile est la moyenne de <span class=\"arabpsy-math-inline\" data-latex=\"x_{(j)}\">$x_{(j)}$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"x_{(j+1)}\">$x_{(j+1)}$<\/span>. L&rsquo;appel optimis\u00e9 de la proc\u00e9dure s&rsquo;effectue comme suit :<\/p>\n<p><code>PROC UNIVARIATE DATA=DONNEES_ORIGINALES NOPRINT;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR Points;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;OUTPUT OUT=METRIQUES_QUANTILES<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Q1=Premier_Quartile<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Q3=Troisieme_Quartile<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;QRANGE=Ecart_Interquartile;<br \/>\nRUN;<\/code><\/p>\n<p>L&rsquo;argument syst\u00e8me <code>NOPRINT<\/code> joue un r\u00f4le d\u00e9terminant dans les environnements de production industrielle : il neutralise la g\u00e9n\u00e9ration volumineuse des tables d&rsquo;\u00e9dition dans le visualiseur ODS, d\u00e9chargeant le sous-syst\u00e8me graphique et acc\u00e9l\u00e9rant drastiquement le temps d&rsquo;ex\u00e9cution CPU sur les m\u00e9ga-bases de donn\u00e9es.<\/p>\n<h3>4.2 Exportation des m\u00e9triques vers une table SAS via l&rsquo;instruction OUTPUT<\/h3>\n<p>L&rsquo;instruction <code>OUTPUT OUT=<\/code> permet de d\u00e9tourner le flux des r\u00e9sultats scalaires calcul\u00e9s par le moteur d&rsquo;inf\u00e9rence de <code>PROC UNIVARIATE<\/code> pour les persister structurellement au sein d&rsquo;une nouvelle table SAS exploitable par des \u00e9tapes avales.<\/p>\n<p>Dans la syntaxe expos\u00e9e ci-dessus, les mots-cl\u00e9s r\u00e9serv\u00e9s du langage SAS sont mobilis\u00e9s avec une rigueur absolue :<\/p>\n<ul>\n<li><code>Q1=Premier_Quartile<\/code> commande la capture de la valeur num\u00e9rique calcul\u00e9e pour le 25e centile empirique.<\/li>\n<li><code>Q3=Troisieme_Quartile<\/code> intercepte la valeur correspondant au 75e centile empirique.<\/li>\n<li><code>QRANGE=Ecart_Interquartile<\/code> ordonne au moteur computationnel de soustraire nativement <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span> de <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span>, \u00e9liminant le risque d&rsquo;erreurs d&rsquo;arrondi interm\u00e9diaires.<\/li>\n<\/ul>\n<p>La table SAS r\u00e9sultante, baptis\u00e9e <code>METRIQUES_QUANTILES<\/code>, est constitu\u00e9e d&rsquo;un enregistrement unique (une seule ligne vectorielle) h\u00e9bergeant trois variables scalaires en double pr\u00e9cision. L&rsquo;ex\u00e9cution d&rsquo;un <code>PROC PRINT<\/code> de contr\u00f4le permet de valider empiriquement les valeurs g\u00e9n\u00e9r\u00e9es pour notre \u00e9chantillon de 16 observations :<\/p>\n<p><code>PROC PRINT DATA=METRIQUES_QUANTILES NOOBS;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"M\u00e9triques de dispersion robuste issues de PROC UNIVARIATE\";<br \/>\nRUN;<\/code><\/p>\n<p>Pour notre s\u00e9rie de donn\u00e9es didactique, le calcul manuel direct sous la d\u00e9finition standard confirme un premier quartile <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span> s&rsquo;\u00e9tablissant \u00e0 <strong>30,0<\/strong>, un troisi\u00e8me quartile <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span> positionn\u00e9 \u00e0 <strong>91,75<\/strong> (ou 89,5 selon la r\u00e8gle d&rsquo;interpolation exacte appliqu\u00e9e sur <span class=\"arabpsy-math-inline\" data-latex=\"N=16\">$N=16$<\/span>), fixant un \u00e9cart interquartile robuste <span class=\"arabpsy-math-inline\" data-latex=\"\\text{IQR}\">$\\text{IQR}$<\/span> pr\u00e9cis\u00e9ment mesur\u00e9 autour de <strong>61,75<\/strong>.<\/p>\n<h3>4.3 Alternative d&rsquo;extraction via le syst\u00e8me ODS OUTPUT<\/h3>\n<p>Une alternative de programmation moderne et hautement modulaire au sein de l&rsquo;\u00e9cosyst\u00e8me SAS r\u00e9side dans l&rsquo;exploitation de l&rsquo;<strong>Output Delivery System (ODS)<\/strong>. L&rsquo;instruction <code>ODS OUTPUT<\/code> intercepte directement les tables d&rsquo;objets graphiques et statistiques g\u00e9n\u00e9r\u00e9es nativement en arri\u00e8re-plan par SAS avant leur formatage visuel final.<\/p>\n<p>L&rsquo;impl\u00e9mentation via ODS s&rsquo;articule comme suit :<\/p>\n<p><code>ODS EXCLUDE ALL;<br \/>\nODS OUTPUT Quantiles=TABLE_QUANTILES_ODS;<br \/>\nPROC UNIVARIATE DATA=DONNEES_ORIGINALES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR Points;<br \/>\nRUN;<br \/>\nODS EXCLUDE NONE;<\/code><\/p>\n<p>L&rsquo;instruction <code>ODS EXCLUDE ALL;<\/code> d\u00e9sactive la restitution visuelle de toutes les tables, \u00e9quivalant fonctionnellement \u00e0 l&rsquo;option <code>NOPRINT<\/code> tout en conservant la capacit\u00e9 du syst\u00e8me d&rsquo;interception d&rsquo;extraire la table d\u00e9nomm\u00e9e <code>Quantiles<\/code>. Cette table ODS pr\u00e9sente une topologie verticale o\u00f9 chaque centile (Min, 1%, 5%, 10%, Q1, M\u00e9diane, Q3, 90%, 95%, 99%, Max) occupe un enregistrement distinct index\u00e9 par une variable descriptive. Bien que cette approche soit d&rsquo;une souplesse incomparable pour construire des m\u00e9tadonn\u00e9es dynamiques multi-percentiles, elle n\u00e9cessite une \u00e9tape de transposition ult\u00e9rieure (via <code>PROC TRANSPOSE<\/code>) pour aligner les quartiles horizontalement sur une seule ligne. L&rsquo;instruction native <code>OUTPUT OUT=<\/code> demeure ainsi la m\u00e9thode la plus directe et la plus efficiente en m\u00e9moire vive pour le calcul cons\u00e9cutif des bornes d&rsquo;exclusion.<\/p>\n<h2>5. Calcul algorithmique des bornes de Tukey dans une \u00e9tape DATA<\/h2>\n<h3>5.1 Fusion cart\u00e9sienne et transmission des seuils critiques<\/h3>\n<p>L&rsquo;obtention de la table scalaire <code>METRIQUES_QUANTILES<\/code> pose un d\u00e9fi d&rsquo;architecture logicielle classique sous SAS : comment propager des statistiques synth\u00e9tiques r\u00e9sum\u00e9es sur une ligne unique vers l&rsquo;int\u00e9gralit\u00e9 des observations de la table originale (qui compte <span class=\"arabpsy-math-inline\" data-latex=\"n = 16\">$n = 16$<\/span> lignes) sans g\u00e9n\u00e9rer une surcharge computationnelle disproportionn\u00e9e ?<\/p>\n<p>Deux strat\u00e9gies algorithmiques pr\u00e9valent dans l&rsquo;ing\u00e9nierie SAS. La premi\u00e8re consiste \u00e0 ex\u00e9cuter une jointure cart\u00e9sienne via le langage SQL (<code>PROC SQL<\/code>). N\u00e9anmoins, l&rsquo;approche la plus \u00e9conome en ressources syst\u00e8me et la plus rapide en temps machine repose sur l&rsquo;utilisation imbriqu\u00e9e de deux instructions <code>SET<\/code> s\u00e9quentielles au sein d&rsquo;une seule \u00e9tape <code>DATA<\/code>, exploitant la persistance des registres dans le <em>Program Data Vector<\/em>.<\/p>\n<p>Consid\u00e9rons l&rsquo;\u00e9tape d&rsquo;agr\u00e9gation suivante :<\/p>\n<p><code>DATA DONNEES_ET_BORNES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;IF _N_ = 1 THEN SET METRIQUES_QUANTILES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SET DONNEES_ORIGINALES;<br \/>\nRUN;<\/code><\/p>\n<p>La logique interne de cette syntaxe repose sur la variable automatique du syst\u00e8me SAS <code>_N_<\/code>, qui trace le num\u00e9ro d&rsquo;it\u00e9ration de la boucle d&rsquo;ex\u00e9cution de l&rsquo;\u00e9tape <code>DATA<\/code>. Lors de la toute premi\u00e8re it\u00e9ration (<code>_N_ = 1<\/code>), l&rsquo;instruction conditionnelle <code>SET METRIQUES_QUANTILES;<\/code> est ex\u00e9cut\u00e9e. SAS lit l&rsquo;unique enregistrement de la table des quantiles et verse les variables <code>Premier_Quartile<\/code>, <code>Troisieme_Quartile<\/code> et <code>Ecart_Interquartile<\/code> dans le PDV. Fondamentalement, parce que cette instruction <code>SET<\/code> n&rsquo;est plus jamais r\u00e9ex\u00e9cut\u00e9e pour les it\u00e9rations suivantes (<span class=\"arabpsy-math-inline\" data-latex=\"_N_ ge 2\">$_N_ ge 2$<\/span>), SAS <strong>ne r\u00e9initialise pas<\/strong> ces scalaires \u00e0 des valeurs manquantes : ils demeurent automatiquement pr\u00e9serv\u00e9s en m\u00e9moire tout au long du balayage de la table principale command\u00e9 par la seconde instruction <code>SET DONNEES_ORIGINALES;<\/code>. Cette m\u00e9thode \u00e9limine le besoin de cl\u00e9s d&rsquo;appariement complexes et optimise le temps de traitement sur des millions d&rsquo;enregistrements.<\/p>\n<h3>5.2 Impl\u00e9mentation math\u00e9matique des seuils inf\u00e9rieur et sup\u00e9rieur<\/h3>\n<p>D\u00e8s lors que les scalaires de dispersion robuste sont inject\u00e9s de fa\u00e7on continue dans le flux vectoriel de chaque observation, le calcul algorithmique des cl\u00f4tures int\u00e9rieure et sup\u00e9rieure de Tukey s&rsquo;effectue par des d\u00e9clarations d&rsquo;assignation lin\u00e9aire directes.<\/p>\n<p>Nous formalisons le calcul dans le prolongement direct de l&rsquo;\u00e9tape <code>DATA<\/code> :<\/p>\n<p><code>DATA DONNEES_CALCULEES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;IF _N_ = 1 THEN SET METRIQUES_QUANTILES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SET DONNEES_ORIGINALES;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* Calcul des cl\u00f4tures formelles de Tukey avec k = 1.5 *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;Borne_Inf = Premier_Quartile - (1.5 * Ecart_Interquartile);<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;Borne_Sup = Troisieme_Quartile + (1.5 * Ecart_Interquartile);<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* M\u00e9triques analytiques d'\u00e9cartement relatif *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;Distance_Inf = Points - Borne_Inf;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;Distance_Sup = Points - Borne_Sup;<br \/>\nRUN;<\/code><\/p>\n<p>Dans ce segment de code, les variables <code>Borne_Inf<\/code> et <code>Borne_Sup<\/code> fixent rigoureusement les seuils critiques de conformit\u00e9 statistique. L&rsquo;introduction des variables d&rsquo;\u00e9cartement <code>Distance_Inf<\/code> et <code>Distance_Sup<\/code> fournit une mesure imm\u00e9diate de l&rsquo;amplitude du d\u00e9bordement \u00e9ventuel : un score <span class=\"arabpsy-math-inline\" data-latex=\"Points\">$Points$<\/span> affichant une valeur positive pour <code>Distance_Sup<\/code> signale arithm\u00e9tiquement une violation de la cl\u00f4ture sup\u00e9rieure de Tukey.<\/p>\n<p>Il importe de noter que les processeurs modernes effectuent les calculs en arithm\u00e9tique flottante 64 bits selon le standard IEEE. Dans certains cas limites rarissimes o\u00f9 une observation est rigoureusement \u00e9gale \u00e0 la cl\u00f4ture calcul\u00e9e aux d\u00e9cimales pr\u00e8s, des impr\u00e9cisions de l&rsquo;ordre de <span class=\"arabpsy-math-inline\" data-latex=\"10^{-16}\">$10^{-16}$<\/span> peuvent induire des erreurs de classification. L&rsquo;analyste pr\u00e9cautionneux peut int\u00e9grer la fonction d&rsquo;arrondi native de SAS, par exemple <code>ROUND(Borne_Sup, 1e-9)<\/code>, pour pr\u00e9munir son pipeline de tout artefact machine.<\/p>\n<h3>5.3 Cr\u00e9ation d&rsquo;indicateurs binaires d&rsquo;exclusion<\/h3>\n<p>La derni\u00e8re \u00e9tape du partitionnement logique consiste \u00e0 encoder l&rsquo;\u00e9tiquetage des observations via un indicateur dichotomique formel (<em>flag<\/em> binaire). Cette variable indicatrice permettra de piloter de mani\u00e8re dynamique l&rsquo;ensemble des requ\u00eates de filtrage, des visualisations conditionnelles et des mod\u00e8les d&rsquo;inf\u00e9rence en aval.<\/p>\n<p>Le code SAS d&rsquo;affectation conditionnelle s&rsquo;articule comme suit :<\/p>\n<p><code>DATA DONNEES_FLAGGEES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SET DONNEES_CALCULEES;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* D\u00e9claration conditionnelle d'exclusion *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;IF Points &lt; Borne_Inf OR Points &gt; Borne_Sup THEN DO;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Valeur_Aberrante = 1;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;IF Points &lt; Borne_Inf THEN Type_Anomalie = \"Inf\u00e9rieure\";<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ELSE IF Points &gt; Borne_Sup THEN Type_Anomalie = \"Sup\u00e9rieure\";<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;END;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;ELSE DO;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Valeur_Aberrante = 0;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Type_Anomalie = \"Conforme  \";<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;END;<br \/>\nRUN;<\/code><\/p>\n<p>Cette impl\u00e9mentation permet non seulement de classer instantan\u00e9ment les individus statistiques en deux groupes herm\u00e9tiques (<span class=\"arabpsy-math-inline\" data-latex=\"0 = \\text{Distributionnelle}\">$0 = \\text{Distributionnelle}$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"1 = \\text{Aberrante}\">$1 = \\text{Aberrante}$<\/span>), mais \u00e9galement de qualifier s\u00e9miologiquement la directionnalit\u00e9 de l&rsquo;anomalie au moyen de la variable textuelle <code>Type_Anomalie<\/code>. Le partitionnement logique est ainsi rendu totalement transparent et auditable pour toute expertise statistique ult\u00e9rieure.<\/p>\n<h2>6. Extraction, filtrage et rapport analytique des donn\u00e9es aberrantes<\/h2>\n<h3>6.1 Filtrage syst\u00e9matique des observations cibles avec l&rsquo;instruction WHERE<\/h3>\n<p>Une fois les indicateurs binaires d&rsquo;exclusion grav\u00e9s dans le mod\u00e8le relationnel de la table, la phase d&rsquo;extraction analytique peut \u00eatre engag\u00e9e. L&rsquo;instruction <code>WHERE<\/code> repr\u00e9sente la commande la plus efficiente du langage SAS pour op\u00e9rer un partitionnement strict du flux de donn\u00e9es, celle-ci agissant directement au niveau du sous-syst\u00e8me d&rsquo;E\/S (Entr\u00e9es\/Sorties) avant le chargement effectif des enregistrements dans le tampon m\u00e9moire.<\/p>\n<p>Nous proc\u00e9dons \u00e0 la scission architectur\u00e9e du jeu de donn\u00e9es en deux tables distinctes : une table d\u00e9di\u00e9e \u00e0 l&rsquo;audit clinique des aberrations isol\u00e9es (<code>TABLE_ANOMALIES<\/code>) et une table assainie r\u00e9serv\u00e9e aux inf\u00e9rences statistiques ult\u00e9rieures (<code>TABLE_PURGEE<\/code>) :<\/p>\n<p><code>DATA TABLE_ANOMALIES TABLE_PURGEE;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SET DONNEES_FLAGGEES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;IF Valeur_Aberrante = 1 THEN OUTPUT TABLE_ANOMALIES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;ELSE OUTPUT TABLE_PURGEE;<br \/>\nRUN;<\/code><\/p>\n<p>Cette syntaxe \u00e0 double sortie s\u00e9pare herm\u00e9tiquement les sous-populations sans n\u00e9cessiter deux lectures s\u00e9quentielles distinctes du disque dur. La table <code>TABLE_ANOMALIES<\/code> contient exclusivement les enregistrements ayant enfreint les cl\u00f4tures de Tukey, tandis que <code>TABLE_PURGEE<\/code> pr\u00e9serve l&rsquo;ensemble des profils r\u00e9guliers tout en conservant la tra\u00e7abilit\u00e9 int\u00e9grale de l&rsquo;op\u00e9ration de s\u00e9lection.<\/p>\n<h3>6.2 \u00c9dition de rapports acad\u00e9miques avec PROC PRINT et PROC REPORT<\/h3>\n<p>L&rsquo;\u00e9dition des r\u00e9sultats de la d\u00e9tection statistique doit ob\u00e9ir \u00e0 des standards visuels et documentaires irr\u00e9prochables pour figurer au sein d&rsquo;un rapport d&rsquo;\u00e9valuation clinique, d&rsquo;une soumission r\u00e9glementaire (FDA\/EMA) ou d&rsquo;un manuscrit scientifique destin\u00e9 \u00e0 une revue \u00e0 comit\u00e9 de lecture.<\/p>\n<p>La proc\u00e9dure <code>PROC REPORT<\/code> permet de g\u00e9n\u00e9rer un tableau formel hautement typographi\u00e9 des valeurs discordantes d\u00e9tect\u00e9es :<\/p>\n<p><code>PROC REPORT DATA=TABLE_ANOMALIES NOWD HEADLINE HEADSKIP;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;COLUMN Equipe Points Premier_Quartile Troisieme_Quartile Borne_Sup Distance_Sup Type_Anomalie;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DEFINE Equipe \/ \"Identifiant \u00c9quipe\" WIDTH=18 CENTER;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DEFINE Points \/ \"Score Empirique\" FORMAT=8.2 WIDTH=15 CENTER;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DEFINE Premier_Quartile \/ \"Q1\" FORMAT=8.2 WIDTH=10 CENTER;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DEFINE Troisieme_Quartile \/ \"Q3\" FORMAT=8.2 WIDTH=10 CENTER;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DEFINE Borne_Sup \/ \"Seuil Tukey (Sup)\" FORMAT=8.2 WIDTH=18 CENTER;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DEFINE Distance_Sup \/ \"D\u00e9passement\" FORMAT=8.2 WIDTH=15 CENTER;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DEFINE Type_Anomalie \/ \"Classification\" WIDTH=15 CENTER;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"Rapport d'audit statistique des observations aberrantes identifi\u00e9es\";<br \/>\nRUN;<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/outlier2.jpg\" alt=\"\" \/><\/figure>\n<p>L&rsquo;ex\u00e9cution de ce rapport r\u00e9v\u00e8le formellement que sur l&rsquo;\u00e9chantillon initial de 16 \u00e9quipes, exactement <strong>deux observations<\/strong> violent la cl\u00f4ture sup\u00e9rieure de Tukey :<\/p>\n<ul>\n<li><strong>\u00c9quipe O :<\/strong> Score empirique de <strong>221<\/strong>, affichant un d\u00e9passement massif sup\u00e9rieur \u00e0 35 points au-del\u00e0 du seuil critique sup\u00e9rieur fix\u00e9 par l&rsquo;algorithme.<\/li>\n<li><strong>\u00c9quipe P :<\/strong> Score empirique de <strong>223<\/strong>, surpassant la borne de cl\u00f4ture avec un niveau d&rsquo;amplitude analogue.<\/li>\n<\/ul>\n<p>Toutes les autres observations (de l&rsquo;\u00e9quipe <code>A<\/code> jusqu&rsquo;\u00e0 l&rsquo;\u00e9quipe <code>N<\/code> dont le score culmine \u00e0 98) demeurent strictement confin\u00e9es \u00e0 l&rsquo;int\u00e9rieur de la zone de validit\u00e9 probabiliste <span class=\"arabpsy-math-inline\" data-latex=\"[B_{\\text{\\inf}}, B_{\\text{\\sup}}]\">$[B_{\\text{\\inf}}, B_{\\text{\\sup}}]$<\/span>. L&rsquo;algorithme a ainsi op\u00e9r\u00e9 une discrimination chirurgicale parfaite, n&rsquo;\u00e9tiquetant que les perturbations extr\u00eames sans amputer la variabilit\u00e9 continue de la s\u00e9rie initiale.<\/p>\n<h3>6.3 Quantification de l&rsquo;impact de l&rsquo;exclusion sur les param\u00e8tres de distribution<\/h3>\n<p>L&rsquo;\u00e9valuation \u00e9pist\u00e9mologique de l&rsquo;extraction des valeurs aberrantes exige de mesurer avec pr\u00e9cision les alt\u00e9rations param\u00e9triques induites par l&rsquo;\u00e9puration de l&rsquo;\u00e9chantillon. L&rsquo;ex\u00e9cution comparative de <code>PROC UNIVARIATE<\/code> sur la table globale non filtr\u00e9e et sur la table purg\u00e9e permet de mettre en lumi\u00e8re l&rsquo;ampleur ph\u00e9nom\u00e9nale des distorsions que faisaient peser les \u00e9quipes <code>O<\/code> et <code>P<\/code> sur l&rsquo;\u00e9difice statistique global.<\/p>\n<p>Le code d&rsquo;analyse comparative s&rsquo;ex\u00e9cute ainsi :<\/p>\n<p><code>PROC UNIVARIATE DATA=DONNEES_ORIGINALES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR Points;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"Param\u00e8tres distributionnels - \u00c9chantillon Brut Int\u00e9gral (N=16)\";<br \/>\nRUN;<\/p>\n<p>PROC UNIVARIATE DATA=TABLE_PURGEE;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR Points;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"Param\u00e8tres distributionnels - \u00c9chantillon Purge (N=14)\";<br \/>\nRUN;<\/code><\/p>\n<p>L&rsquo;examen des estimations fournit un enseignement saisissant quant \u00e0 la vuln\u00e9rabilit\u00e9 des statistiques descriptives traditionnelles :<\/p>\n<ul>\n<li><strong>Moyenne arithm\u00e9tique :<\/strong> Dans l&rsquo;\u00e9chantillon global brut (<span class=\"arabpsy-math-inline\" data-latex=\"N=16\">$N=16$<\/span>), la moyenne arithm\u00e9tique culmine \u00e0 <strong>73,75 points<\/strong>. D\u00e8s lors que les deux observations discordantes sont retir\u00e9es (<span class=\"arabpsy-math-inline\" data-latex=\"N=14\">$N=14$<\/span>), la moyenne s&rsquo;effondre \u00e0 <strong>52,43 points<\/strong>, r\u00e9v\u00e9lant un d\u00e9placement artificiel de plus de 21 points (soit une sur\u00e9valuation initiale de plus de 40 %) imputable \u00e0 seulement deux individus !<\/li>\n<li><strong>Variance et \u00c9cart-type :<\/strong> L&rsquo;\u00e9cart-type d&rsquo;\u00e9chantillon passe de <strong>66,88 points<\/strong> dans la table compl\u00e8te \u00e0 <strong>26,55 points<\/strong> dans la table purg\u00e9e. La dispersion a ainsi \u00e9t\u00e9 artificiellement gonfl\u00e9e d&rsquo;un facteur 2,5 par la pr\u00e9sence des anomalies, ce qui aurait rendu tout test param\u00e9trique subs\u00e9quent (test <span class=\"arabpsy-math-inline\" data-latex=\"t\">$t$<\/span> de Student ou ANOVA) totalement incapable de d\u00e9tecter un effet exp\u00e9rimental significatif du fait de l&rsquo;explosion du d\u00e9nominateur d&rsquo;erreur r\u00e9siduelle.<\/li>\n<li><strong>Coefficient d&rsquo;asym\u00e9trie (Skewness) :<\/strong> Le coefficient de Fisher-Pearson plonge d&rsquo;un niveau critique de <strong>+1,78<\/strong> (d\u00e9non\u00e7ant une asym\u00e9trie droite extr\u00eame incompatible avec les tests gaussiens) \u00e0 une valeur parfaitement admissible de <strong>+0,21<\/strong>, restaurant une configuration remarquablement sym\u00e9trique conforme aux exigences de la mod\u00e9lisation param\u00e9trique standard.<\/li>\n<\/ul>\n<h2>7. Validation et diagnostic visuel avec PROC SGPLOT<\/h2>\n<h3>7.1 Construction et param\u00e9trage du diagramme en bo\u00eete (Boxplot)<\/h3>\n<p>Le diagnostic visuel constitue le prolongement heuristique indispensable de la d\u00e9tection num\u00e9rique. En vertu du principe d&rsquo;Anscombe selon lequel des r\u00e9sum\u00e9s statistiques identiques peuvent masquer des architectures empiriques totalement dissemblables, la visualisation graphique permet \u00e0 l&rsquo;analyste de valider la pertinence ph\u00e9nom\u00e9nologique des seuils calcul\u00e9s. Sous SAS, la proc\u00e9dure <code>PROC SGPLOT<\/code> (<em>Statistical Graphics Plot<\/em>) incarne l&rsquo;\u00e9tat de l&rsquo;art du moteur graphique ODS Graphics.<\/p>\n<p>Le <strong>diagramme en bo\u00eete<\/strong> &mdash; ou bo\u00eete \u00e0 moustaches de Tukey (<em>boxplot<\/em>) &mdash; est l&rsquo;expression g\u00e9om\u00e9trique directe de l&rsquo;\u00e9cart interquartile. Au sein de ce graphique, le rectangle central d\u00e9limite le premier quartile <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span>, la m\u00e9diane <span class=\"arabpsy-math-inline\" data-latex=\"Q_2\">$Q_2$<\/span> (marqu\u00e9e par un segment transversal) et le troisi\u00e8me quartile <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span>. Les segments lin\u00e9aires ext\u00e9rieurs &mdash; les \u00ab moustaches \u00bb &mdash; ne s&rsquo;\u00e9tendent pas arbitrairement jusqu&rsquo;au minimum et au maximum de l&rsquo;\u00e9chantillon ; sous la convention stricte de Tukey respect\u00e9e par SAS, elles se d\u00e9ploient jusqu&rsquo;\u00e0 l&rsquo;observation empirique la plus extr\u00eame qui demeure <em>\u00e0 l&rsquo;int\u00e9rieur<\/em> des cl\u00f4tures d&rsquo;exclusion (<span class=\"arabpsy-math-inline\" data-latex=\"[Q_1 - 1,5\\times \\text{IQR} ,;, Q_3 + 1,5\\times \\text{IQR}]\">$[Q_1 &#8211; 1,5\\times \\text{IQR} ,;, Q_3 + 1,5\\times \\text{IQR}]$<\/span>). Toute donn\u00e9e localis\u00e9e au-del\u00e0 des moustaches est dessin\u00e9e sous la forme de points discrets isol\u00e9s.<\/p>\n<p>Le code SAS pour concevoir un boxplot conforme aux normes graphiques acad\u00e9miques les plus strictes s&rsquo;\u00e9tablit comme suit :<\/p>\n<p><code>ODS GRAPHICS ON \/ RESET IMAGENAME=\"Boxplot_Tukey\" IMAGEFMT=PNG WIDTH=12CM HEIGHT=10CM;<br \/>\nPROC SGPLOT DATA=DONNEES_ORIGINALES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VBOX Points \/ BOXWIDTH=0.4<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;DATALABEL=Equipe<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;MEANATTRS=(SYMBOL=Diamond COLOR=DarkRed SIZE=9)<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;OUTLIERATTRS=(SYMBOL=CircleFilled COLOR=Red SIZE=9);<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;YAXIS LABEL=\"Distribution des Scores (Points)\" GRID;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;XAXIS DISPLAY=(NOLABEL);<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"Repr\u00e9sentation sch\u00e9matique par bo\u00eete \u00e0 moustaches de Tukey\";<br \/>\nRUN;<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/outlier3.jpg\" alt=\"\" \/><\/figure>\n<p>L&rsquo;attribut <code>DATALABEL=Equipe<\/code> ordonne l&rsquo;\u00e9tiquetage direct des points extr\u00eames par leur identifiant cat\u00e9goriel, mat\u00e9rialisant instantan\u00e9ment les \u00e9quipes <code>O<\/code> et <code>P<\/code> comme deux marqueurs flottant \u00e0 une distance spectaculaire au-dessus de la moustache terminale sup\u00e9rieure.<\/p>\n<h3>7.2 Visualisation conjointe par nuage de points et lignes de r\u00e9f\u00e9rence<\/h3>\n<p>Afin de procurer une lisibilit\u00e9 cognitive maximale aux auditeurs de l&rsquo;analyse, il est hautement avantageux de superposer la distribution univari\u00e9e brute des observations individuelles avec les cl\u00f4tures formelles calcul\u00e9es pr\u00e9c\u00e9demment. Cette mise en regard conjointe s&rsquo;op\u00e8re par la combinaison des instructions <code>SCATTER<\/code> et <code>REFLINE<\/code> au sein de <code>PROC SGPLOT<\/code> :<\/p>\n<p><code>PROC SGPLOT DATA=DONNEES_FLAGGEES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;\/* Nuage de dispersion univari\u00e9e avec diff\u00e9rentiation chromatique *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SCATTER X=Equipe Y=Points \/ GROUP=Valeur_Aberrante<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;MARKERATTRS=(SYMBOL=CircleFilled SIZE=11)<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;DATALABEL=Points;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* Tra\u00e7age des cl\u00f4tures de Tukey comme seuils d'exclusion *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;REFLINE 184.375 \/ AXIS=Y LINEATTRS=(PATTERN=Dash COLOR=DarkRed THICKNESS=2)<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;LABEL=\"Cl\u00f4ture Sup\u00e9rieure de Tukey (Q3 + 1.5*IQR)\"<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;LABELPOS=Top;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;REFLINE -62.625 \/ AXIS=Y LINEATTRS=(PATTERN=Dash COLOR=Blue THICKNESS=2)<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;LABEL=\"Cl\u00f4ture Inf\u00e9rieure de Tukey (Q1 - 1.5*IQR)\"<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;LABELPOS=Bottom;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;YAXIS LABEL=\"Valeurs observ\u00e9es et cl\u00f4tures th\u00e9oriques\" GRID;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;XAXIS LABEL=\"Modalit\u00e9s cat\u00e9gorielles (\u00c9quipes)\";<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"Cartographie diagnostique de l'espace de d\u00e9cision univari\u00e9\";<br \/>\nRUN;<\/code><\/p>\n<p>Ce graphique produit une d\u00e9monstration p\u00e9remptoire : la ligne rouge tiret\u00e9e mat\u00e9rialisant la cl\u00f4ture sup\u00e9rieure s\u00e9pare visuellement le continuum r\u00e9gulier des points exp\u00e9rimentaux des deux observations fautives, justifiant sans ambigu\u00eft\u00e9 les d\u00e9cisions de traitement en aval.<\/p>\n<h3>7.3 Diagnostic graphique de la normalit\u00e9 avec PROC UNIVARIATE et QQ-Plot<\/h3>\n<p>L&rsquo;\u00e9valuation de la conformit\u00e9 distributionnelle s&rsquo;ach\u00e8ve rigoureusement par l&rsquo;examen du <strong>graphique quantile-quantile (QQ-Plot)<\/strong>. Cet outil confronte directement les quantiles de la distribution empirique observ\u00e9e aux quantiles th\u00e9oriques attendus d&rsquo;une loi normale id\u00e9ale de m\u00eame esp\u00e9rance et de m\u00eame variance (droite de r\u00e9f\u00e9rence de Henry).<\/p>\n<p>L&rsquo;instruction <code>QQPLOT<\/code> sous <code>PROC UNIVARIATE<\/code> automatise cette confrontation analytique :<\/p>\n<p><code>PROC UNIVARIATE DATA=DONNEES_ORIGINALES NOPRINT;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;QQPLOT Points \/ NORMAL(MU=EST SIGMA=EST COLOR=DarkRed THICKNESS=2)<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;SQUARE;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;TITLE \"Diagramme Quantile-Quantile (QQ-Plot) de conformit\u00e9 gaussienne\";<br \/>\nRUN;<\/code><\/p>\n<p>L&rsquo;interpr\u00e9tation diagnostique d&rsquo;un QQ-plot ob\u00e9it \u00e0 des r\u00e9gularit\u00e9s g\u00e9om\u00e9triques bien document\u00e9es en psychom\u00e9trie :<\/p>\n<ul>\n<li>Si les donn\u00e9es proc\u00e8dent d&rsquo;un processus g\u00e9n\u00e9rateur rigoureusement gaussien, les points empiriques s&rsquo;alignent le long de la diagonale th\u00e9orique d&rsquo;ajustement.<\/li>\n<li>La pr\u00e9sence de valeurs aberrantes sup\u00e9rieures se traduit par un d\u00e9collement vertical brutal des points situ\u00e9s dans l&rsquo;extr\u00e9mit\u00e9 droite sup\u00e9rieure du rep\u00e8re, les points s&rsquo;\u00e9cartant vers le haut par rapport \u00e0 la droite de normalit\u00e9. C&rsquo;est exactement la structure cin\u00e9matique observ\u00e9e avec les points 221 et 223, attestant d&rsquo;une queue de distribution anormalement lourde (leptokurticit\u00e9 pathologique).<\/li>\n<\/ul>\n<h2>8. M\u00e9thodes univari\u00e9es alternatives sous SAS : Scores Z et \u00c9cart Absolu M\u00e9dian (MAD)<\/h2>\n<h3>8.1 D\u00e9tection standardis\u00e9e via la proc\u00e9dure PROC STDIZE<\/h3>\n<p>Si la m\u00e9thode de l&rsquo;\u00e9cart interquartile s&rsquo;av\u00e8re particuli\u00e8rement adapt\u00e9e \u00e0 l&rsquo;exploration robuste de donn\u00e9es asym\u00e9triques, l&rsquo;usage des scores standardis\u00e9s demeure un standard industriel pr\u00e9pond\u00e9rant dans de nombreux protocoles d&rsquo;assurance qualit\u00e9. La proc\u00e9dure <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/support.sas.com\/documentation\/onlinedoc\/stat\/142\/stdize.pdf\"><code>PROC STDIZE<\/code><\/a> permet d&rsquo;effectuer des transformations d&rsquo;\u00e9chelle et des calculs de d\u00e9viance \u00e0 une vitesse de traitement optimale.<\/p>\n<p>L&rsquo;impl\u00e9mentation de la r\u00e8gle empirique des scores Z traditionnels s&rsquo;ex\u00e9cute selon la syntaxe suivante :<\/p>\n<p><code>PROC STDIZE DATA=DONNEES_ORIGINALES OUT=SCORES_Z METHOD=STD;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR Points;<br \/>\nRUN;<\/p>\n<p>DATA SCORES_Z_ANALYSE;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SET SCORES_Z (RENAME=(Points=Score_Z));<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;IF ABS(Score_Z) &gt; 2.5 THEN Aberrant_Z = 1;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;ELSE Aberrant_Z = 0;<br \/>\nRUN;<\/code><\/p>\n<p>L&rsquo;option <code>METHOD=STD<\/code> soustrait la moyenne arithm\u00e9tique d&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{x}\">$\\bar{x}$<\/span> et divise par l&rsquo;\u00e9cart-type standard <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>. N\u00e9anmoins, l&rsquo;exp\u00e9rimentation met ici imm\u00e9diatement en lumi\u00e8re l&rsquo;<strong>effet de masque<\/strong> th\u00e9oris\u00e9 dans la section 1.3 :<\/p>\n<p>Dans notre \u00e9chantillon didactique, du fait de l&rsquo;explosion de l&rsquo;\u00e9cart-type d&rsquo;\u00e9chantillon \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"s = 66,88\">$s = 66,88$<\/span> sous l&rsquo;effet conjugu\u00e9 des deux anomalies, le score Z de l&rsquo;\u00e9quipe <code>O<\/code> (221 points) s&rsquo;\u00e9tablit \u00e0 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Z_O = \\frac{221 - 73,75}{66,88} \\approx +2,20\">$$Z_O = \\frac{221 &#8211; 73,75}{66,88} \\approx +2,20$$<\/span><\/p>\n<p>De m\u00eame, le score Z de l&rsquo;\u00e9quipe <code>P<\/code> (223 points) atteint seulement :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Z_P = \\frac{223 - 73,75}{66,88} \\approx +2,23\">$$Z_P = \\frac{223 &#8211; 73,75}{66,88} \\approx +2,23$$<\/span><\/p>\n<p><strong>Constat critique :<\/strong> Si le chercheur applique aveugl\u00e9ment le seuil d\u00e9cisionnel conventionnel de <span class=\"arabpsy-math-inline\" data-latex=\"|Z| &gt; 2,5\">$|Z| &gt; 2,5$<\/span> (ou <em>a fortiori<\/em> <span class=\"arabpsy-math-inline\" data-latex=\"|Z| &gt; 3,0\">$|Z| &gt; 3,0$<\/span>), <strong>aucune<\/strong> des deux observations aberrantes n&rsquo;est d\u00e9tect\u00e9e ! Les valeurs 221 et 223 ont tellement gonfl\u00e9 le d\u00e9nominateur de standardisation qu&rsquo;elles se sont auto-dissimul\u00e9es derri\u00e8re la variance qu&rsquo;elles ont engendr\u00e9e. Cet \u00e9chec d\u00e9monstratif souligne la dangerosit\u00e9 m\u00e9thodologique de l&#8217;emploi des scores Z classiques en pr\u00e9sence d&rsquo;anomalies s\u00e9v\u00e8res.<\/p>\n<h3>8.2 Approche robuste non param\u00e9trique avec l&rsquo;\u00c9cart Absolu M\u00e9dian (MAD)<\/h3>\n<p>Pour restaurer l&rsquo;efficience de la standardisation sans subir le pi\u00e8ge du masquage, les statisticiens Boris Iglewicz et David Hoaglin (1993) pr\u00e9conisent le recours au <strong>score Z modifi\u00e9<\/strong>, fond\u00e9 sur la m\u00e9diane et l&rsquo;<strong>\u00c9cart Absolu M\u00e9dian (MAD &mdash; <em>Median Absolute Deviation<\/em>)<\/strong>. La MAD est math\u00e9matiquement formul\u00e9e comme la m\u00e9diane des \u00e9carts absolus par rapport \u00e0 la m\u00e9diane de l&rsquo;\u00e9chantillon :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{MAD} = \\text{m\u00e9diane}(|x_i - \\tilde{x}|)\">$$\\text{MAD} = \\text{m\u00e9diane}(|x_i &#8211; \\tilde{x}|)$$<\/span><\/p>\n<p>Pour assurer la convergence asymptotique de la MAD vers l&rsquo;\u00e9cart-type <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma\">$\\sigma$<\/span> sous l&rsquo;hypoth\u00e8se d&rsquo;une loi normale, la statistique est multipli\u00e9e par une constante d&rsquo;\u00e9chelle universelle <span class=\"arabpsy-math-inline\" data-latex=\"b = 1,4826\">$b = 1,4826$<\/span> (<span class=\"arabpsy-math-inline\" data-latex=\"1 \/ \\Phi^{-1}(0,75) \\approx 1,4826\">$1 \/ \\Phi^{-1}(0,75) \\approx 1,4826$<\/span>). Le score Z modifi\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"M_i\">$M_i$<\/span>) s&rsquo;\u00e9nonce alors :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"M_i = \\frac{0,6745 \\times (x_i - \\tilde{x})}{\\text{MAD}}\">$$M_i = \\frac{0,6745 \\times (x_i &#8211; \\tilde{x})}{\\text{MAD}}$$<\/span><\/p>\n<p>Sous SAS, <code>PROC STDIZE<\/code> supporte nativement la normalisation par la MAD gr\u00e2ce \u00e0 l&rsquo;option de m\u00e9thode robuste <code>METHOD=MAD<\/code> :<\/p>\n<p><code>PROC STDIZE DATA=DONNEES_ORIGINALES OUT=SCORES_MAD METHOD=MAD;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR Points;<br \/>\nRUN;<\/code><\/p>\n<p>Sous cette transformation, la m\u00e9diane de l&rsquo;\u00e9chantillon (<span class=\"arabpsy-math-inline\" data-latex=\"\\tilde{x} = 51,0\">$\\tilde{x} = 51,0$<\/span>) et la MAD r\u00e9sistent totalement aux valeurs 221 et 223. Le d\u00e9nominateur demeure stable et repr\u00e9sentatif de la variabilit\u00e9 centrale des 14 observations r\u00e9guli\u00e8res. En appliquant la r\u00e8gle d&rsquo;Iglewicz et Hoaglin recommandant un seuil critique d&rsquo;exclusion de <span class=\"arabpsy-math-inline\" data-latex=\"|M_i| &gt; 3,5\">$|M_i| &gt; 3,5$<\/span>, les \u00e9quipes <code>O<\/code> et <code>P<\/code> affichent des scores modifi\u00e9s exc\u00e9dant largement 5,0, assurant une d\u00e9tection sans faille, exempte de tout effet de masque.<\/p>\n<h3>8.3 Synth\u00e8se d\u00e9cisionnelle pour le choix de la m\u00e9thode univari\u00e9e<\/h3>\n<p>Le choix de l&rsquo;instrument computationnel pour l&rsquo;analyse univari\u00e9e ne doit pas s&rsquo;effectuer au hasard mais d\u00e9couler d&rsquo;un arbitrage formel int\u00e9grant la taille d&rsquo;\u00e9chantillon, la sym\u00e9trie distributionnelle et le point de rupture des m\u00e9triques. Le tableau d\u00e9cisionnel ci-apr\u00e8s synth\u00e9tise les crit\u00e8res d&rsquo;orientation pour l&rsquo;analyste SAS :<\/p>\n<ul>\n<li><strong>R\u00e8gle des scores Z classiques (PROC STDIZE METHOD=STD) :<\/strong> Applicable uniquement si <span class=\"arabpsy-math-inline\" data-latex=\"n &gt; 1000\">$n &gt; 1000$<\/span> avec une preuve formelle de normalit\u00e9 sans contamination pr\u00e9alable ; Point de rupture nul (<span class=\"arabpsy-math-inline\" data-latex=\"0%\">$0%$<\/span>) ; Risque maximal d&rsquo;effet de masque.<\/li>\n<li><strong>Cl\u00f4tures de Tukey \/ IQR (PROC UNIVARIATE Q1 Q3) :<\/strong> Applicable pour toutes tailles d&rsquo;\u00e9chantillon (<span class=\"arabpsy-math-inline\" data-latex=\"n ge 10\">$n ge 10$<\/span>) ; Hypoth\u00e8se de normalit\u00e9 non requise ; Point de rupture robuste de <span class=\"arabpsy-math-inline\" data-latex=\"25%\">$25%$<\/span> ; M\u00e9thode de r\u00e9f\u00e9rence pour l&rsquo;analyse exploratoire et les flux de production standard.<\/li>\n<li><strong>Score Z modifi\u00e9 \/ MAD (PROC STDIZE METHOD=MAD) :<\/strong> Id\u00e9al pour les \u00e9chantillons mod\u00e9r\u00e9s (<span class=\"arabpsy-math-inline\" data-latex=\"20 le n le 500\">$20 le n le 500$<\/span>) avec suspicion de contaminations multiples s\u00e9v\u00e8res ; Point de rupture optimal (<span class=\"arabpsy-math-inline\" data-latex=\"50%\">$50%$<\/span>) ; R\u00e9sistance absolue au masquage.<\/li>\n<\/ul>\n<h2>9. Extension aux anomalies multivari\u00e9es : Distance de Mahalanobis et PROC ROBUSTREG<\/h2>\n<h3>9.1 Limites intrins\u00e8ques de l&rsquo;\u00e9valuation univari\u00e9e en pr\u00e9sence de covariations<\/h3>\n<p>L&rsquo;\u00e9valuation marginale variable par variable &mdash; aussi perfectionn\u00e9e soit-elle via les cl\u00f4tures de Tukey ou la MAD &mdash; se r\u00e9v\u00e8le structurellement incapable d&rsquo;intercepter les <strong>valeurs aberrantes multivari\u00e9es<\/strong>. Une observation multivari\u00e9e aberrante est un vecteur multidimensionnel qui, lorsqu&rsquo;il est projet\u00e9 individuellement sur chacun des axes orthogonaux pris isol\u00e9ment, se fond parfaitement au c\u0153ur de chaque distribution univari\u00e9e, mais dont la position conjointe viole radicalement la matrice de corr\u00e9lation ou de covariance du syst\u00e8me.<\/p>\n<p>Consid\u00e9rons l&rsquo;exemple classique d&rsquo;une \u00e9tude biom\u00e9trique combinant l&rsquo;\u00e2ge et la tension art\u00e9rielle systolique. Un \u00e2ge de 12 ans constitue une observation courante ; une tension art\u00e9rielle systolique de 160 mmHg est une mesure fr\u00e9quemment observ\u00e9e chez des patients hypertendus \u00e2g\u00e9s. Les deux grandeurs se situent parfaitement \u00e0 l&rsquo;int\u00e9rieur des cl\u00f4tures de Tukey univari\u00e9es de leurs distributions respectives. En revanche, l&rsquo;occurrence simultan\u00e9e d&rsquo;un individu \u00e2g\u00e9 de 12 ans pr\u00e9sentant une tension systolique de 160 mmHg repr\u00e9sente une anomalie m\u00e9dicale et relationnelle majeure. Les algorithmes univari\u00e9s d\u00e9crits jusqu&rsquo;ici sont aveugles \u00e0 ce ph\u00e9nom\u00e8ne.<\/p>\n<h3>9.2 Calcul matriciel de la distance de Mahalanobis avec PROC CANDISC ou PROC PRINCOMP<\/h3>\n<p>Pour mesurer l&rsquo;\u00e9loignement multidimensionnel d&rsquo;un vecteur d&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_i = (x_{i1}, dots, x_{ip})&#039;\">$\\mathbf{x}_i = (x_{i1}, dots, x_{ip})&rsquo;$<\/span> par rapport au centre de gravit\u00e9 vectoriel de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{\\mathbf{x}}\">$\\bar{\\mathbf{x}}$<\/span>, la distance euclidienne standard est inappropri\u00e9e car elle ignore les corr\u00e9lations inter-variables et traite chaque \u00e9chelle de fa\u00e7on isotrope. Le statisticien Prasanta Chandra Mahalanobis a r\u00e9solu ce probl\u00e8me en forgeant la <strong>distance de Mahalanobis<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"D_i^2\">$D_i^2$<\/span>), qui pond\u00e8re l&rsquo;\u00e9cart multidimensionnel par l&rsquo;inverse de la matrice de variance-covariance d&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}\">$\\mathbf{S}$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"D_i^2 = (\\mathbf{x}_i - \\bar{\\mathbf{x}})&#039; \\mathbf{S}^{-1} (\\mathbf{x}_i - \\bar{\\mathbf{x}})\">$$D_i^2 = (\\mathbf{x}_i &#8211; \\bar{\\mathbf{x}})&rsquo; \\mathbf{S}^{-1} (\\mathbf{x}_i &#8211; \\bar{\\mathbf{x}})$$<\/span><\/p>\n<p>Sous l&rsquo;hypoth\u00e8se de multinormalit\u00e9 des donn\u00e9es, le carr\u00e9 de la distance de Mahalanobis converge asymptotiquement vers une loi du Khi-deux dont les degr\u00e9s de libert\u00e9 correspondent au nombre <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> de dimensions vectorielles (<span class=\"arabpsy-math-inline\" data-latex=\"D^2 \\sim \\chi^2_p\">$D^2 \\sim \\chi^2_p$<\/span>). Une observation pr\u00e9sentant une valeur de <span class=\"arabpsy-math-inline\" data-latex=\"D_i^2\">$D_i^2$<\/span> sup\u00e9rieure au centile critique <span class=\"arabpsy-math-inline\" data-latex=\"\\chi^2_{p, 1 - \\alpha}\">$\\chi^2_{p, 1 &#8211; \\alpha}$<\/span> (avec typiquement <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,001\">$\\alpha = 0,001$<\/span>) est diagnostiqu\u00e9e comme une anomalie multivari\u00e9e.<\/p>\n<p>Sous SAS, la distance de Mahalanobis peut \u00eatre calcul\u00e9e de fa\u00e7on \u00e9l\u00e9gante en exploitant la proc\u00e9dure <code>PROC PRINCOMP<\/code> (analyse en composantes principales) ou <code>PROC CANDISC<\/code> :<\/p>\n<p><code>\/* Exemple d'extraction multivari\u00e9e sur une table bivari\u00e9e *\/<br \/>\nPROC PRINCOMP DATA=DONNEES_MULTIVAR STD OUT=SCORES_PCA NOPRINT;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR Variable_X Variable_Y;<br \/>\nRUN;<\/p>\n<p>DATA DISTANCES_MAHALANOBIS;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SET SCORES_PCA;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;\/* La somme des carr\u00e9s des scores factoriels standardis\u00e9s \u00e9quivaut \u00e0 D^2 *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;Distance_Mahalanobis_Sq = Prin1**2 + Prin2**2;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;P_Valeur = 1 - PROBCHI(Distance_Mahalanobis_Sq, 2);<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;IF P_Valeur &lt; 0.001 THEN Anomalie_Multivariee = 1;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;ELSE Anomalie_Multivariee = 0;<br \/>\nRUN;<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/outlier4.jpg\" alt=\"\" \/><\/figure>\n<p>Cette m\u00e9thodologie permet de sonder en profondeur l&rsquo;espace factoriel et d&rsquo;extraire les observations multivari\u00e9es discordantes qui \u00e9chappent aux filtres univari\u00e9s s\u00e9quentiels.<\/p>\n<h3>9.3 R\u00e9gression robuste et M-estimateurs avec PROC ROBUSTREG<\/h3>\n<p>Lorsque le chercheur cherche \u00e0 mod\u00e9liser une variable r\u00e9ponse en pr\u00e9sence simultan\u00e9e de points aberrants dans les variables explicatives (<em>points de levier<\/em> &mdash; <em>leverage points<\/em>) et dans la dimension r\u00e9siduelle (<em>outliers de r\u00e9gression<\/em>), les MCO \u00e9chouent catastrophiquement. SAS d\u00e9ploie une solution de niveau mondial au travers de la proc\u00e9dure <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/support.sas.com\/documentation\/onlinedoc\/stat\/142\/robustreg.pdf\"><code>PROC ROBUSTREG<\/code><\/a>.<\/p>\n<p><code>PROC ROBUSTREG<\/code> fournit quatre classes d&rsquo;estimateurs robustes de pointe :<\/p>\n<ul>\n<li><strong>M-estimateurs (METHOD=M) :<\/strong> Introduits par Peter Huber (1973), ils modifient la fonction de co\u00fbt quadratique des MCO en appliquant une fonction de perte lin\u00e9aire au-del\u00e0 d&rsquo;un certain seuil d&rsquo;\u00e9cart r\u00e9siduel, limitant l&rsquo;influence des grandes anomalies r\u00e9siduelles.<\/li>\n<li><strong>LTS-estimateurs (METHOD=LTS) :<\/strong> La m\u00e9thode des moindres carr\u00e9s tronqu\u00e9s (<em>Least Trimmed Squares<\/em>) de Peter Rousseeuw minimise la somme des carr\u00e9s ordonn\u00e9s des seuls <span class=\"arabpsy-math-inline\" data-latex=\"h\">$h$<\/span> plus petits r\u00e9sidus, offrant un point de rupture de 50 %.<\/li>\n<li><strong>S-estimateurs (METHOD=S) :<\/strong> Estimateurs fond\u00e9s sur la minimisation d&rsquo;une \u00e9chelle robuste des r\u00e9sidus.<\/li>\n<li><strong>MM-estimateurs (METHOD=MM) :<\/strong> M\u00e9thodologie combinant la haute r\u00e9sistance des estimateurs S avec la haute efficacit\u00e9 asymptotique des M-estimateurs sous le mod\u00e8le normal.<\/li>\n<\/ul>\n<p>L&rsquo;impl\u00e9mentation op\u00e9rationnelle de la r\u00e9gression robuste sous SAS s&rsquo;effectue via le code suivant :<\/p>\n<p><code>PROC ROBUSTREG DATA=DONNEES_COMPLEXES METHOD=MM PLOTS=(RDPLOT LEVERAGE);<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;MODEL Reponse = Predicteur1 Predicteur2 \/ DIAGNOSTICS;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;OUTPUT OUT=REG_ROBUSTE_DIAG RESIDUAL=Residu_Robuste WEIGHT=Poids_Final;<br \/>\nRUN;<\/code><\/p>\n<p>L&rsquo;instruction <code>OUTPUT<\/code> g\u00e9n\u00e8re une variable fondamentale : <code>Poids_Final<\/code>. Les algorithmes it\u00e9ratifs de <code>PROC ROBUSTREG<\/code> assignent dynamiquement un poids compris entre 0 et 1 \u00e0 chaque individu statistique : une observation totalement conforme re\u00e7oit un poids de 1, tandis qu&rsquo;une anomalie extr\u00eame voit son poids \u00e9cras\u00e9 vers 0, la neutralisant int\u00e9gralement lors de l&rsquo;estimation des coefficients sans exiger de suppression manuelle subjective de la part du chercheur.<\/p>\n<h2>10. Traitement m\u00e9thodologique et \u00e9thique des donn\u00e9es aberrantes identifi\u00e9es<\/h2>\n<h3>10.1 Typologie des arbitrages post-d\u00e9tection : Suppression, Imputation, Winsorisation<\/h3>\n<p>L&rsquo;identification algorithmique rigoureuse des observations aberrantes ne marque pas le terme de la d\u00e9marche d&rsquo;ing\u00e9nierie statistique ; elle inaugure la phase critique de l&rsquo;arbitrage m\u00e9thodologique. L&rsquo;analyste se trouve confront\u00e9 \u00e0 quatre alternatives cardinales de traitement, dont chacune emporte des cons\u00e9quences th\u00e9oriques et d\u00e9ontologiques majeures :<\/p>\n<p>La premi\u00e8re voie r\u00e9side dans la <strong>d\u00e9l\u00e9tion pure et simple (suppression cas par cas)<\/strong>. Cette option n&rsquo;est scientifiquement recevable que si l&rsquo;anomalie r\u00e9sulte sans \u00e9quivoque d&rsquo;une erreur d&rsquo;encodage irr\u00e9cup\u00e9rable, d&rsquo;une rupture physique d&rsquo;instrumentation ou d&rsquo;une violation patente des crit\u00e8res d&rsquo;admissibilit\u00e9 de l&rsquo;\u00e9chantillonnage. En dehors de ces justifications m\u00e9caniques document\u00e9es, l&rsquo;exclusion arbitraire est assimilable \u00e0 une falsification opportuniste des distributions (<em>cherry-picking<\/em>).<\/p>\n<p>La deuxi\u00e8me alternative est la <strong>Winsorisation (ou rabattement de Tukey)<\/strong>, nomm\u00e9e en l&rsquo;honneur du bio-statisticien Charles Winsor. Cette approche conserve l&rsquo;int\u00e9gralit\u00e9 de l&rsquo;effectif d&rsquo;\u00e9chantillonnage <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span>, mais remplace la valeur d&rsquo;une observation exc\u00e9dant un centile donn\u00e9 (ou une cl\u00f4ture de Tukey) par la valeur exacte de ce seuil critique. Par exemple, un score situ\u00e9 au-del\u00e0 de la cl\u00f4ture sup\u00e9rieure <span class=\"arabpsy-math-inline\" data-latex=\"B_{\\text{\\sup}}\">$B_{\\text{\\sup}}$<\/span> est arbitrairement tronqu\u00e9 \u00e0 la valeur de <span class=\"arabpsy-math-inline\" data-latex=\"B_{\\text{\\sup}}\">$B_{\\text{\\sup}}$<\/span>. La winsorisation pr\u00e9serve la cardinalit\u00e9 des donn\u00e9es tout en jugulant l&rsquo;influence disproportionn\u00e9e des queues lourdes sur la variance.<\/p>\n<p>La troisi\u00e8me approche rel\u00e8ve de l&rsquo;<strong>imputation statistique robuste<\/strong>, o\u00f9 la valeur anormale av\u00e9r\u00e9e est convertie en donn\u00e9e manquante, puis r\u00e9estim\u00e9e par des algorithmes d&rsquo;imputation multiple par \u00e9quations cha\u00een\u00e9es (MICE via <code>PROC MI<\/code> et <code>PROC MIANALYZE<\/code> sous SAS).<\/p>\n<p>La quatri\u00e8me voie, hautement pr\u00e9conis\u00e9e par l&rsquo;\u00e9pist\u00e9mologie moderne, consiste en la <strong>conservation int\u00e9grale<\/strong> des donn\u00e9es brutes, conjugu\u00e9e \u00e0 l&rsquo;abandon pur et simple des estimateurs gaussiens au profit de mod\u00e9lisations non param\u00e9triques (tests de Wilcoxon-Mann-Whitney via <code>PROC NPAR1WAY<\/code>) ou de mod\u00e9lisations robustes (<code>PROC ROBUSTREG<\/code>).<\/p>\n<h3>10.2 Impl\u00e9mentation programmatique de la Winsorisation sous SAS<\/h3>\n<p>L&rsquo;impl\u00e9mentation algorithmique de la Winsorisation s&rsquo;effectue avec une remarquable concision au sein d&rsquo;une \u00e9tape <code>DATA<\/code> aliment\u00e9e par les bornes extraites de <code>PROC UNIVARIATE<\/code>. En reprenant notre table de calcul didactique, le script de plafonnement s&rsquo;\u00e9tablit comme suit :<\/p>\n<p><code>DATA DONNEES_WINSORISEES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;SET DONNEES_CALCULEES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;LENGTH Statut_Winsorisation $ 15;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;Points_Winsorises = Points;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;IF Points &gt; Borne_Sup THEN DO;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Points_Winsorises = Borne_Sup;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Statut_Winsorisation = \"Plafonn\u00e9 (Sup)\";<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;END;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;ELSE IF Points &lt; Borne_Inf THEN DO;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Points_Winsorises = Borne_Inf;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Statut_Winsorisation = \"Planch\u00e9 (Inf)\";<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;END;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;ELSE DO;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Statut_Winsorisation = \"Inchang\u00e9\";<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;END;<br \/>\nRUN;<\/code><\/p>\n<p>Dans ce mod\u00e8le r\u00e9assign\u00e9, les \u00e9quipes <code>O<\/code> (221) et <code>P<\/code> (223) voient leur score num\u00e9rique plafonn\u00e9 exactement \u00e0 la valeur num\u00e9rique de la cl\u00f4ture sup\u00e9rieure de Tukey, neutralisant leur effet de traction sur la variance d&rsquo;\u00e9chantillon tout en maintenant un effectif de <span class=\"arabpsy-math-inline\" data-latex=\"N = 16\">$N = 16$<\/span>. Cette technique offre un compromis math\u00e9matique de haut niveau pour stabiliser les mod\u00e8les de pr\u00e9vision.<\/p>\n<h3>10.3 Analyses de sensibilit\u00e9 et transparence dans le reporting scientifique<\/h3>\n<p>Le standard d&rsquo;or (<em>gold standard<\/em>) de l&rsquo;int\u00e9grit\u00e9 m\u00e9thodologique contemporaine r\u00e9side dans la r\u00e9alisation syst\u00e9matique d&rsquo;une <strong>analyse de sensibilit\u00e9<\/strong> (<em>sensitivity analysis<\/em>). Cette approche impose de soumettre l&rsquo;ensemble du pipeline inf\u00e9rentiel \u00e0 une double ex\u00e9cution en parall\u00e8le : une premi\u00e8re fois sur le jeu de donn\u00e9es brut int\u00e9gral comprenant la totalit\u00e9 des observations, et une seconde fois sur le jeu de donn\u00e9es \u00e9pur\u00e9 des observations discordantes identifi\u00e9es par les cl\u00f4tures de Tukey.<\/p>\n<p>La communication scientifique conforme aux exigences APA 7e \u00e9dition implique de rapporter textuellement les deux \u00e9tats du mod\u00e8le :<\/p>\n<blockquote>\n<p><em>\u00ab Le diagnostic pr\u00e9liminaire des valeurs univari\u00e9es extr\u00eames au moyen de la m\u00e9thode de l&rsquo;\u00e9cart interquartile de Tukey (k = 1,5) a mis en \u00e9vidence deux observations discordantes majeures au sein des scores de groupe (modalit\u00e9 O = 221 et modalit\u00e9 P = 223, exc\u00e9dant la cl\u00f4ture sup\u00e9rieure fix\u00e9e \u00e0 184,38). L&rsquo;ajustement du mod\u00e8le de r\u00e9gression lin\u00e9aire initial incluant la totalit\u00e9 de l&rsquo;\u00e9chantillon (N = 16) r\u00e9v\u00e8le un param\u00e8tre de pente non significatif (beta = 0,24, p = 0,18, R2 = 0,12). La r\u00e9estimation du mod\u00e8le apr\u00e8s exclusion syst\u00e9matique des deux valeurs aberrantes (N = 14) atteste d&rsquo;une inversion structurelle de la significativit\u00e9 statistique (beta = 0,68, p &lt; 0,001, R2 = 0,58). Conform\u00e9ment aux directives d&#039;int\u00e9grit\u00e9 m\u00e9thodologique, les deux mod\u00e8les sont document\u00e9s de mani\u00e8re concurrente. \u00bb<\/em><\/p>\n<\/blockquote>\n<p>L&rsquo;annexion publique du script SAS reproductible &mdash; contenant l&rsquo;horodatage exact et la table de tra\u00e7abilit\u00e9 des exclusions &mdash; constitue la condition sine qua non de la confiance \u00e9ditoriale et de la reproductibilit\u00e9 computationnelle des travaux de recherche.<\/p>\n<h2>11. Automatisation du pipeline de d\u00e9tection via le langage Macro SAS<\/h2>\n<h3>11.1 Conception modulaire d&rsquo;une macro param\u00e9trable %DETECT_IQR<\/h3>\n<p>Dans les environnements industriels et les centres de recherche clinique traitant des dizaines de tables comprenant des centaines de variables quantitatives, le codage s\u00e9quentiel manuel des \u00e9tapes <code>DATA<\/code> et <code>PROC UNIVARIATE<\/code> devient rapidement impraticable, source d&rsquo;erreurs d&rsquo;encodage et contraire au principe de non-redondance (DRY &mdash; <em>Don&rsquo;t Repeat Yourself<\/em>). Le processeur de macro-programmation de SAS fournit un outil d&rsquo;abstraction puissant permettant d&rsquo;encapsuler la logique de Tukey au sein d&rsquo;une macro universelle param\u00e9trable.<\/p>\n<p>La conception logicielle de la macro, baptis\u00e9e <code>%DETECT_IQR<\/code>, int\u00e8gre quatre param\u00e8tres dynamiques indispensables :<\/p>\n<ul>\n<li><code>DATA_IN :<\/code> Nom de la table SAS source \u00e0 inspecter.<\/li>\n<li><code>VAR_ANALYSE :<\/code> Nom de la variable quantitative soumise au filtrage.<\/li>\n<li><code>COEF_IQR :<\/code> Coefficient multiplicateur des cl\u00f4tures de Tukey (avec 1,5 par d\u00e9faut, extensible \u00e0 3,0 pour les anomalies extr\u00eames).<\/li>\n<li><code>DATA_OUT_FLAG :<\/code> Nom de la table finale g\u00e9n\u00e9r\u00e9e, h\u00e9bergeant les indicateurs binaires de classification.<\/li>\n<\/ul>\n<h3>11.2 D\u00e9veloppement du code macro et gestion des variables macro\u00e9conomiques<\/h3>\n<p>La compilation et l&rsquo;ex\u00e9cution du pipeline modulaire s&rsquo;impl\u00e9mentent au travers de la macro-commande suivante, exploitant les interfaces de m\u00e9moire globale et la routine <code>CALL SYMPUTX<\/code> :<\/p>\n<p><code>%MACRO DETECT_IQR(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DATA_IN=<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;,VAR_ANALYSE=<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;,COEF_IQR=1.5<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;,DATA_OUT_FLAG=TABLE_AUDITEE<br \/>\n);<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;%PUT NOTE: ====================================================;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;%PUT NOTE: INITIALISATION DE LA MACRO DETECT_IQR SUR &amp;VAR_ANALYSE ;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;%PUT NOTE: ====================================================;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* 1. Extraction computationnelle des quantiles sans affichage *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;PROC UNIVARIATE DATA=&amp;DATA_IN NOPRINT;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;VAR &amp;VAR_ANALYSE;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;OUTPUT OUT=_WORK_QUANTILES<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Q1=_M_Q1<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Q3=_M_Q3<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;QRANGE=_M_IQR;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;RUN;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* 2. \u00c9criture des m\u00e9triques dans des macro-variables globales *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DATA _NULL_;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;SET _WORK_QUANTILES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;CALL SYMPUTX('M_Q1', _M_Q1, 'G');<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;CALL SYMPUTX('M_Q3', _M_Q3, 'G');<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;CALL SYMPUTX('M_IQR', _M_IQR, 'G');<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;CALL SYMPUTX('M_B_INF', _M_Q1 - (&amp;COEF_IQR * _M_IQR), 'G');<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;CALL SYMPUTX('M_B_SUP', _M_Q3 + (&amp;COEF_IQR * _M_IQR), 'G');<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;RUN;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;%PUT NOTE: VALEURS CALCUL\u00c9ES : Q1=&amp;M_Q1 | Q3=&amp;M_Q3 | IQR=&amp;M_IQR ;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;%PUT NOTE: CL\u00d4TURE INF\u00c9RIEURE : &amp;M_B_INF ;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;%PUT NOTE: CL\u00d4TURE SUP\u00c9RIEURE : &amp;M_B_SUP ;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* 3. Injection et balisage conditionnel *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DATA &amp;DATA_OUT_FLAG;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;SET &amp;DATA_IN;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Borne_Inf_Seuil = &amp;M_B_INF;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Borne_Sup_Seuil = &amp;M_B_SUP;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;IF NOT MISSING(&amp;VAR_ANALYSE) THEN DO;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;IF &amp;VAR_ANALYSE &lt; &amp;M_B_INF OR &amp;VAR_ANALYSE &gt; &amp;M_B_SUP THEN Aberrant_Flag = 1;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ELSE Aberrant_Flag = 0;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;END;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ELSE Aberrant_Flag = .;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;RUN;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;\/* 4. Purge des tables de travail temporaires *\/<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;PROC DATASETS LIB=WORK NOLIST;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;DELETE _WORK_QUANTILES;<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;QUIT;<\/p>\n<p>&nbsp;&nbsp;&nbsp;&nbsp;%PUT NOTE: TRAITEMENT MACRO %DETECT_IQR ACHEV\u00c9 AVEC SUCC\u00c8S.;<br \/>\n%MEND DETECT_IQR;<\/code><\/p>\n<p>L&rsquo;utilisation de la routine <code>CALL SYMPUTX<\/code> avec le troisi\u00e8me argument positionn\u00e9 \u00e0 <code>'G'<\/code> garantit l&rsquo;affectation formelle dans la table des symboles globaux, \u00e9vitant toute collision de port\u00e9e lexicale. L&rsquo;instruction d\u00e9fensive <code>IF NOT MISSING(&amp;VAR_ANALYSE)<\/code> prot\u00e8ge le script contre la classification erron\u00e9e des valeurs absentes.<\/p>\n<h3>11.3 Application concr\u00e8te de la macro sur des volum\u00e9tries massives<\/h3>\n<p>L&rsquo;ex\u00e9cution de la macro sur notre \u00e9chantillon didactique s&rsquo;op\u00e8re par une ligne de commande d&rsquo;une extr\u00eame simplicit\u00e9 :<\/p>\n<p><code>%DETECT_IQR(<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DATA_IN=DONNEES_ORIGINALES,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;VAR_ANALYSE=Points,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;COEF_IQR=1.5,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;DATA_OUT_FLAG=RESULTAT_AUTOMATISE<br \/>\n);<\/code><\/p>\n<p>Pour industrialiser le diagnostic sur l&rsquo;int\u00e9gralit\u00e9 des variables num\u00e9riques d&rsquo;une base de donn\u00e9es complexe comprenant des centaines d&rsquo;items, l&rsquo;analyste peut coupler cette macro avec le dictionnaire de m\u00e9tadonn\u00e9es interne de SAS (<code>SASHELP.VCOLUMN<\/code>). Au moyen d&rsquo;une boucle g\u00e9n\u00e9r\u00e9e en SQL via l&rsquo;instruction <code>SELECT INTO :MACRO_LIST SEPARATED BY ' '<\/code>, la macro <code>%DETECT_IQR<\/code> peut \u00eatre it\u00e9r\u00e9e sur des milliers de colonnes de fa\u00e7on totalement autonome, stockant chaque observation discordante dans une table d&rsquo;audit consolid\u00e9e pr\u00eate pour les analyses de conformit\u00e9 r\u00e9glementaire.<\/p>\n<h2>12. Synth\u00e8se pratique, pi\u00e8ges courants et recommandations op\u00e9rationnelles<\/h2>\n<h3>12.1 \u00c9cueils programmatiques fr\u00e9quents sous SAS et solutions<\/h3>\n<p>Au fil du d\u00e9ploiement de pipelines industriels sous SAS, un ensemble d&rsquo;erreurs computationnelles et de pi\u00e8ges conceptuels r\u00e9currents peuvent compromettre la validit\u00e9 des d\u00e9tections. L&rsquo;analyste doit imp\u00e9rativement ma\u00eetriser les solutions pr\u00e9ventives associ\u00e9es :<\/p>\n<ul>\n<li><strong>Le pi\u00e8ge de l&rsquo;\u00e9valuation logique des valeurs manquantes :<\/strong> Dans le langage SAS, une valeur manquante num\u00e9rique (repr\u00e9sent\u00e9e par un point <code>.<\/code>) est encod\u00e9e en m\u00e9moire physique sous la forme du plus petit nombre \u00e0 virgule flottante possible (\u00e9quivalant \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"-\\infty\">$-\\infty$<\/span>). Par cons\u00e9quent, si un analyste r\u00e9dige na\u00efvement la clause conditionnelle :<br \/>\n    <br \/><code>IF Points &lt; Borne_Inf THEN Valeur_Aberrante = 1;<\/code><br \/>\n    toutes les valeurs manquantes de la base de donn\u00e9es voient la condition arithm\u00e9tique valid\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"-\\infty &lt; B_{\\text{\\inf}}\">$-\\infty &lt; B_{\\text{\\inf}}$<\/span>) et sont <strong>\u00e9tiquet\u00e9es \u00e0 tort comme des anomalies inf\u00e9rieures extr\u00eames<\/strong> ! Il est imp\u00e9ratif d&rsquo;inclure syst\u00e9matiquement la garde conditionnelle <code>IF NOT MISSING(Points) AND Points &lt; Borne_Inf THEN...<\/code>.<\/li>\n<li><strong>La confusion des algorithmes d&rsquo;interpolation (PCTLDEF) :<\/strong> Si une \u00e9quipe compare des r\u00e9sultats obtenus sous SAS avec ceux calcul\u00e9s sous R ou Python, des divergences d\u00e9cimales peuvent \u00e9merger si l&rsquo;option <code>PCTLDEF=<\/code> n&rsquo;est pas harmonis\u00e9e. SAS utilise par d\u00e9faut <code>PCTLDEF=5<\/code>, tandis que la fonction <code>quantile()<\/code> par d\u00e9faut sous R applique la d\u00e9finition de type 7 (correspondant \u00e0 <code>PCTLDEF=4<\/code> sous SAS). Dans le cadre d&rsquo;\u00e9tudes cliniques crois\u00e9es, sp\u00e9cifier formellement <code>PCTLDEF=4<\/code> ou <code>PCTLDEF=5<\/code> garantit la reproductibilit\u00e9 inter-logicielle.<\/li>\n<li><strong>L&rsquo;invariance d&rsquo;arrondi sur les bornes calcul\u00e9es :<\/strong> Lors des jointures ou des transferts via des macro-variables \u00e9crites sous forme de texte brut, des troncatures de d\u00e9cimales peuvent faire basculer une valeur rigoureusement \u00e9gale \u00e0 la cl\u00f4ture du c\u00f4t\u00e9 des observations exclues. L&rsquo;usage de la fonction d&rsquo;arrondi <code>ROUND()<\/code> ou la conservation stricte des variables sous forme num\u00e9rique continue dans les \u00e9tapes <code>DATA<\/code> \u00e9vite ce risque.<\/li>\n<\/ul>\n<h3>12.2 Guide d&rsquo;audit m\u00e9thodique pour l&rsquo;analyste de donn\u00e9es<\/h3>\n<p>Pour institutionnaliser un processus de d\u00e9tection robuste et auditable au sein d&rsquo;un laboratoire de recherche ou d&rsquo;un service de biostatistique, il est vivement recommand\u00e9 d&rsquo;appliquer une check-list s\u00e9quentielle invariable avant toute exclusion d\u00e9finitive d&rsquo;observations :<\/p>\n<p><strong>\u00c9tape 1 : Audit d&rsquo;encodage et de conformit\u00e9 m\u00e9trologique :<\/strong> V\u00e9rifier syst\u00e9matiquement si l&rsquo;observation atypique r\u00e9sulte d&rsquo;une erreur mat\u00e9rielle d&rsquo;acquisition. Contr\u00f4ler les bordereaux sources (CRF), interroger le syst\u00e8me de capture et corriger l&rsquo;erreur \u00e0 la racine si une preuve mat\u00e9rielle d&rsquo;invalidation est \u00e9tablie.<\/p>\n<p><strong>\u00c9tape 2 : Exploration univari\u00e9e visuelle et distributionnelle :<\/strong> \u00c9valuer le coefficient d&rsquo;asym\u00e9trie de Fisher-Pearson (<span class=\"arabpsy-math-inline\" data-latex=\"Skewness\">$Skewness$<\/span>). Examiner la bo\u00eete \u00e0 moustaches (<code>PROC SGPLOT VBOX<\/code>) et le diagramme quantile-quantile (<code>PROC UNIVARIATE QQPLOT<\/code>). Si l&rsquo;asym\u00e9trie est s\u00e9v\u00e8re (<span class=\"arabpsy-math-inline\" data-latex=\"Skewness &gt; 1,5\">$Skewness &gt; 1,5$<\/span>), privil\u00e9gier l&rsquo;\u00e9cart interquartile ajust\u00e9 par la m\u00e9dcouple ou op\u00e9rer une transformation logarithmique pr\u00e9alable avant d&rsquo;appliquer les cl\u00f4tures de Tukey.<\/p>\n<p><strong>\u00c9tape 3 : Calcul robuste des cl\u00f4tures de Tukey :<\/strong> Ex\u00e9cuter l&rsquo;\u00e9tape <code>PROC UNIVARIATE<\/code> coupl\u00e9e \u00e0 l&rsquo;\u00e9tape <code>DATA<\/code> pour borner la distribution (<span class=\"arabpsy-math-inline\" data-latex=\"k = 1,5\">$k = 1,5$<\/span>). V\u00e9rifier l&rsquo;absence de masquage en comparant les r\u00e9sultats aux scores Z modifi\u00e9s issus de <code>PROC STDIZE METHOD=MAD<\/code>.<\/p>\n<p><strong>\u00c9tape 4 : Diagnostic d&rsquo;influence multivari\u00e9e :<\/strong> V\u00e9rifier si les observations univari\u00e9es conformes ne cachent pas des anomalies dimensionnelles conjointes en estimant les distances de Mahalanobis ou en ex\u00e9cutant une r\u00e9gression robuste sous <code>PROC ROBUSTREG<\/code>.<\/p>\n<p><strong>\u00c9tape 5 : Analyse de sensibilit\u00e9 syst\u00e9matique :<\/strong> Estimer les mod\u00e8les statistiques cibles (r\u00e9gressions, ANOVA, tests d&rsquo;hypoth\u00e8se) sur la table brute int\u00e9grale, puis sur la table \u00e9pur\u00e9e ou winsoris\u00e9e. Si les conclusions scientifiques divergent, documenter exhaustivement les deux \u00e9tats au sein du rapport final.<\/p>\n<h3>12.3 Perspectives computationnelles futures et int\u00e9gration en science des donn\u00e9es<\/h3>\n<p>L&rsquo;identification des valeurs aberrantes conna\u00eet une transformation profonde \u00e0 l&rsquo;\u00e8re du <em>Big Data<\/em> et de l&rsquo;intelligence artificielle. Si les cl\u00f4tures de Tukey et l&rsquo;\u00e9cart interquartile demeurent des fondations analytiques in\u00e9branlables pour l&rsquo;explicabilit\u00e9 m\u00e9thodologique et les environnements r\u00e9glement\u00e9s, les proc\u00e9dures traditionnelles de SAS sont d\u00e9sormais enrichies par l&rsquo;\u00e9cosyst\u00e8me <strong>SAS Viya<\/strong> et son moteur in-memory <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/support.sas.com\/en\/software\/cloud-analytic-services-cas.html\">Cloud Analytic Services (CAS)<\/a>.<\/p>\n<p>Dans cet espace massivement parall\u00e9lis\u00e9, des proc\u00e9dures d&rsquo;apprentissage statistique non supervis\u00e9 &mdash; telles que <code>PROC FASTCLUS<\/code> pour le partitionnement spatial, <code>PROC SVMACHINE<\/code> pour les s\u00e9parateurs \u00e0 vaste marge \u00e0 une classe (<em>One-Class Support Vector Machines<\/em>) et les algorithmes de <strong>for\u00eats d&rsquo;isolement (Isolation Forests)<\/strong> &mdash; permettent d&rsquo;isoler des anomalies g\u00e9om\u00e9triques hautement non lin\u00e9aires au sein d&rsquo;espaces \u00e0 plusieurs milliers de dimensions. N\u00e9anmoins, ces approches algorithmiques avanc\u00e9es ne dispensent aucunement l&rsquo;analyste de donn\u00e9es de la rigueur math\u00e9matique et conceptuelle fondamentale : le traitement des observations discordantes demeurera toujours une discipline \u00e0 la crois\u00e9e de l&rsquo;alg\u00e8bre probabiliste, de la sagacit\u00e9 informatique et de l&rsquo;\u00e9thique de la recherche.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Barnett, V., &amp; Lewis, T. (1994). <em>Outliers in Statistical Data<\/em> (3rd ed.). John Wiley &amp; Sons.<\/li>\n<li>Hoaglin, D. C., Mosteller, F., &amp; Tukey, J. W. (Eds.). (1983). <em>Understanding Robust and Exploratory Data Analysis<\/em>. John Wiley &amp; Sons.<\/li>\n<li>Hubert, M., &amp; Vandervieren, E. (2008). An adjusted boxplot for skewed distributions. <em>Computational Statistics &amp; Data Analysis<\/em>, 52(12), 5186\u20135201. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/j.csda.2007.11.008\">https:\/\/doi.org\/10.1016\/j.csda.2007.11.008<\/a><\/li>\n<li>Iglewicz, B., &amp; Hoaglin, D. C. (1993). <em>How to Detect and Handle Outliers<\/em> (Vol. 16, The ASQC Basic References in Quality Control). American Society for Quality Control.<\/li>\n<li>Mahalanobis, P. C. (1936). On the generalized distance in statistics. <em>Proceedings of the National Institute of Sciences of India<\/em>, 2(1), 49\u201355.<\/li>\n<li>Rousseeuw, P. J., &amp; Leroy, A. M. (1987). <em>Robust Regression and Outlier Detection<\/em>. John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471725382\">https:\/\/doi.org\/10.1002\/0471725382<\/a><\/li>\n<li>SAS Institute Inc. (2020). <em>SAS\/STAT 15.2 User&rsquo;s Guide: The UNIVARIATE Procedure<\/em>. SAS Institute Inc. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/support.sas.com\/documentation\/onlinedoc\/stat\/152\/univariate.pdf\">https:\/\/support.sas.com\/documentation\/onlinedoc\/stat\/152\/univariate.pdf<\/a><\/li>\n<li>SAS Institute Inc. (2020). <em>SAS\/STAT 15.2 User&rsquo;s Guide: The ROBUSTREG Procedure<\/em>. SAS Institute Inc. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/support.sas.com\/documentation\/onlinedoc\/stat\/152\/robustreg.pdf\">https:\/\/support.sas.com\/documentation\/onlinedoc\/stat\/152\/robustreg.pdf<\/a><\/li>\n<li>Tukey, J. W. (1977). <em>Exploratory Data Analysis<\/em>. Addison-Wesley.<\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour identifier les valeurs aberrantes dans SAS avec la m\u00e9thode de l&rsquo;\u00e9cart interquartile (IQR), codes et exemples pratiques.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[9499,359],"tags":[11408,11407,8812,9575,11406],"class_list":["post-9985","psych_stat","type-psych_stat","status-publish","hentry","category-programmation-sas","category-statistiques-appliquees","tag-clotures-de-tukey","tag-ecart-interquartile-iqr","tag-nettoyage-des-donnees","tag-proc-univariate","tag-valeurs-aberrantes-sas"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/9985","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=9985"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/9985\/revisions"}],"predecessor-version":[{"id":9990,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/9985\/revisions\/9990"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=9985"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=9985"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=9985"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}