{"id":2226,"date":"2026-10-14T18:00:50","date_gmt":"2026-10-14T18:00:50","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-ecart-interquartile-dans-r-exemples\/"},"modified":"2026-10-14T18:00:50","modified_gmt":"2026-10-14T18:00:50","slug":"comment-calculer-ecart-interquartile-dans-r-exemples","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-ecart-interquartile-dans-r-exemples","title":{"rendered":"Comment calculer l&rsquo;\u00e9cart interquartile dans R (avec des exemples)"},"content":{"rendered":"<p>Dans le paysage contemporain des sciences comportementales, de la neuropsychologie et de la psychom\u00e9trie quantitative, l&rsquo;int\u00e9grit\u00e9 de l&rsquo;analyse exploratoire conditionne la validit\u00e9 des inf\u00e9rences statistiques. Les chercheurs manipulent fr\u00e9quemment des distributions de donn\u00e9es empiriques qui s&rsquo;\u00e9cartent des postulats th\u00e9oriques d&rsquo;une courbe en cloche id\u00e9ale. Qu&rsquo;il s&rsquo;agisse de latences de r\u00e9ponses cognitives, de scores d&rsquo;inventaires de d\u00e9tresse psychologique caract\u00e9ris\u00e9s par une asym\u00e9trie marqu\u00e9e, ou de biomarqueurs comportementaux sujets \u00e0 des contaminations sporadiques, l&rsquo;usage non critique de la moyenne arithm\u00e9tique et de l&rsquo;\u00e9cart-type expose l&rsquo;analyste \u00e0 des distorsions m\u00e9trologiques substantielles. L&rsquo;\u00e9valuation rigoureuse de la dispersion centrale d&rsquo;une s\u00e9rie d&rsquo;observations n\u00e9cessite alors des outils non param\u00e9triques r\u00e9sistants, capables de r\u00e9sumer la variabilit\u00e9 sans \u00eatre corrompus par la morphologie aberrante des extr\u00e9mit\u00e9s de distribution.<\/p>\n<p>L&rsquo;\u00e9cart interquartile (commun\u00e9ment d\u00e9sign\u00e9 sous l&rsquo;acronyme anglophone <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Interquartile_range\">IQR pour <em>Interquartile Range<\/em><\/a>) s&rsquo;impose comme la m\u00e9trique de r\u00e9f\u00e9rence pour quantifier l&rsquo;\u00e9talement des cinquante pour cent m\u00e9dians d&rsquo;une distribution statistique. Th\u00e9oris\u00e9 dans les travaux fondateurs de l&rsquo;analyse exploratoire des donn\u00e9es initi\u00e9s par John Tukey, cet indicateur d&rsquo;ordre isole le noyau informationnel d&rsquo;un \u00e9chantillon en s&rsquo;affranchissant du comportement stochastique des queues de distribution. Alors que l&rsquo;\u00e9cart-type pr\u00e9sente un point de rupture nul \u2014 une seule valeur infiniment \u00e9loign\u00e9e suffisant \u00e0 d\u00e9stabiliser irr\u00e9m\u00e9diablement son estimation \u2014, l&rsquo;\u00e9cart interquartile offre une r\u00e9sistance intrins\u00e8que remarquable face aux observations extr\u00eames, garantissant une stabilit\u00e9 m\u00e9thodologique essentielle lors de l&rsquo;investigation de populations cliniques h\u00e9t\u00e9rog\u00e8nes ou d&rsquo;exp\u00e9rimentations psychophysiques exigeantes.<\/p>\n<p>L&rsquo;environnement statistique R constitue aujourd&rsquo;hui la plateforme pr\u00e9dominante pour la mod\u00e9lisation et l&rsquo;\u00e9valuation psychom\u00e9trique en libre acc\u00e8s. Bien que la mise en \u0153uvre de la commande native <code>IQR()<\/code> paraisse triviale au premier abord, sa manipulation experte requiert une compr\u00e9hension aigu\u00eb de ses fondements formels, de sa gestion rigide des donn\u00e9es manquantes et, surtout, de ses algorithmes d&rsquo;interpolation sous-jacents. Le pr\u00e9sent trait\u00e9 propose une d\u00e9construction exhaustive du calcul de l&rsquo;\u00e9cart interquartile au sein de l&rsquo;\u00e9cosyst\u00e8me R, alliant bases th\u00e9oriques, cas pratiques issus de protocoles de recherche et recommandations \u00e9ditoriales strictes, conform\u00e9ment aux standards de l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\/\">American Psychological Association<\/a>.<\/p>\n<h2>1. Introduction conceptuelle \u00e0 l&rsquo;\u00e9cart interquartile (IQR) en psychologie quantitative<\/h2>\n<h3>1.1 D\u00e9finition th\u00e9orique et pertinence de l&rsquo;intervalle interquartile<\/h3>\n<p>L&rsquo;\u00e9cart interquartile se d\u00e9finit formellement comme la distance s\u00e9parant le troisi\u00e8me quartile empirique (<span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span>, \u00e9quivalant au 75e centile) du premier quartile empirique (<span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span>, \u00e9quivalant au 25e centile) au sein d&rsquo;un ensemble ordonn\u00e9 d&rsquo;observations r\u00e9elles. Sur le plan arithm\u00e9tique, la formulation s&rsquo;exprime sous la forme :<\/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>Cette diff\u00e9rence absolue d\u00e9limite l&rsquo;\u00e9tendue du domaine num\u00e9rique h\u00e9bergeant exactement la moiti\u00e9 centrale des observations collect\u00e9es. En psychom\u00e9trie, o\u00f9 les instruments de mesure quantifient des \u00e9tats latents tels que l&rsquo;anxi\u00e9t\u00e9 situationnelle, le sentiment d&rsquo;auto-efficacit\u00e9 ou les capacit\u00e9s de m\u00e9moire de travail, la distribution des scores observ\u00e9s ne suit que tr\u00e8s rarement une fonction de densit\u00e9 normale pure. L&rsquo;apparition d&rsquo;effets de plancher (caract\u00e9ris\u00e9s par un regroupement massif de participants aux valeurs minimales de l&rsquo;\u00e9chelle) ou d&rsquo;effets de plafond rend le recours \u00e0 la variance et \u00e0 l&rsquo;\u00e9cart-type th\u00e9oriquement trompeur. Dans ces conditions non gaussiennes, l&rsquo;\u00e9cart-type surestime l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des sujets typiques en \u00e9tant math\u00e9matiquement gonfl\u00e9 par l&rsquo;\u00e9loignement quadratique des scores p\u00e9riph\u00e9riques.<\/p>\n<p>L&rsquo;application typique de l&rsquo;IQR se manifeste avec une acuit\u00e9 particuli\u00e8re dans la mesure des temps de r\u00e9action informatis\u00e9s (exprim\u00e9s en millisecondes) au sein de paradigmes attentionnels informatis\u00e9s tels que la t\u00e2che de Stroop ou la t\u00e2che de d\u00e9tection de signaux de Posner. Dans ces protocoles, les latences brutes d\u00e9montrent invariablement une distribution fortement \u00e9tir\u00e9e vers la droite, d\u00e9sign\u00e9e techniquement sous le nom de distribution ex-gaussienne ou log-normale. L&rsquo;\u00e9cart interquartile refl\u00e8te de mani\u00e8re fid\u00e8le la dispersion des vitesses d&rsquo;ex\u00e9cution chez les participants sans \u00eatre fauss\u00e9 par les inattentions momentan\u00e9es \u2014 qui g\u00e9n\u00e8rent des temps artificiellement longs \u2014 ou par les r\u00e9ponses anticipatoires accidentelles. De m\u00eame, lors de l&rsquo;administration d&rsquo;inventaires cliniques, tels que le <a rel=\"noopener noreferrer\" href=\"https:\/\/www.apa.org\/pi\/about\/publications\/caregivers\/practice-settings\/assessment\/tools\/beck-depression\">Beck Depression Inventory<\/a> au sein d&rsquo;un \u00e9chantillon issu de la population g\u00e9n\u00e9rale, la majeure partie des r\u00e9pondants obtient des scores nuls ou marginaux, g\u00e9n\u00e9rant une asym\u00e9trie positive abrupte o\u00f9 seul l&rsquo;\u00e9cart interquartile associ\u00e9 \u00e0 la m\u00e9diane peut fournir une synth\u00e8se fid\u00e8le des tendances du groupe.<\/p>\n<h3>1.2 Propri\u00e9t\u00e9s m\u00e9trologiques et robustesse statistique<\/h3>\n<p>La valeur fondamentale de l&rsquo;\u00e9cart interquartile repose sur ses remarquables propri\u00e9t\u00e9s de robustesse statistique, plus pr\u00e9cis\u00e9ment articul\u00e9es autour de son point de rupture empirique. En th\u00e9orie de l&rsquo;estimation robuste, le point de rupture d\u00e9signe la proportion minimale d&rsquo;observations arbitrairement erron\u00e9es qu&rsquo;un \u00e9chantillon peut contenir avant que l&rsquo;estimateur consid\u00e9r\u00e9 ne produise une valeur infinie ou d\u00e9nu\u00e9e de sens. L&rsquo;\u00e9cart-type classique poss\u00e8de un point de rupture asymptotique de <span class=\"arabpsy-math-inline\" data-latex=\"0%\">$0%$<\/span>, ce qui signifie qu&rsquo;une seule observation contamin\u00e9e \u2014 issue par exemple d&rsquo;une erreur d&rsquo;encodage manuel transformant un score de d\u00e9pression de 12 en 120 \u2014 suffit \u00e0 vicier sans limite l&rsquo;indice de dispersion. \u00c0 l&rsquo;inverse, l&rsquo;IQR tol\u00e8re jusqu&rsquo;\u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"25%\">$25%$<\/span> de donn\u00e9es extr\u00eames \u00e0 chaque extr\u00e9mit\u00e9 de la distribution sans que son estimation des cinquante pour cent centraux ne soit corrompue, lui conf\u00e9rant un point de rupture effectif de <span class=\"arabpsy-math-inline\" data-latex=\"25%\">$25%$<\/span>.<\/p>\n<p>Cette insensibilit\u00e9 relative permet \u00e0 l&rsquo;IQR de conserver toute sa validit\u00e9 descriptive en pr\u00e9sence d&rsquo;une asym\u00e9trie de distribution (ou <em>skewness<\/em>) prononc\u00e9e. Lorsque la masse des probabilit\u00e9s se d\u00e9place asym\u00e9triquement vers une borne de l&rsquo;\u00e9chelle, les quartiles <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> se d\u00e9calent en coh\u00e9rence avec la densit\u00e9 empirique, maintenant une mesure d&rsquo;\u00e9talement qui n&rsquo;est pas distendue de mani\u00e8re non lin\u00e9aire par les r\u00e9sidus distaux. Ce comportement justifie son int\u00e9gration syst\u00e9matique dans les protocoles d&rsquo;Analyse Exploratoire des Donn\u00e9es (<em>Exploratory Data Analysis<\/em> ou EDA) promus par le statisticien John Tukey d\u00e8s les ann\u00e9es 1970. L&rsquo;EDA pr\u00e9conise d&rsquo;\u00e9couter la structure brute des donn\u00e9es avant d&rsquo;imposer des mod\u00e8les probabilistes rigides aux ph\u00e9nom\u00e8nes observ\u00e9s.<\/p>\n<p>N\u00e9anmoins, la robustesse m\u00e9trologique s&rsquo;accompagne d&rsquo;une contrepartie \u00e9pist\u00e9mologique inh\u00e9rente : l&rsquo;occultation compl\u00e8te du comportement des queues de distribution. Deux jeux de donn\u00e9es cliniques distincts peuvent pr\u00e9senter des valeurs strictement identiques pour <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span>, la m\u00e9diane et <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span> \u2014 et par voie de cons\u00e9quence un IQR rigoureusement identique \u2014, tout en affichant des structures de dispersion radicalement disparates au-del\u00e0 du 75e centile ou en de\u00e7\u00e0 du 25e centile. L&rsquo;analyste doit donc garder \u00e0 l&rsquo;esprit que l&rsquo;\u00e9cart interquartile sacrifie d\u00e9lib\u00e9r\u00e9ment l&rsquo;information relative \u00e0 l&rsquo;ampleur absolue des extr\u00eames afin d&rsquo;offrir une stabilit\u00e9 optimale sur le segment interm\u00e9diaire.<\/p>\n<h3>1.3 Pertinence m\u00e9thodologique dans la recherche exp\u00e9rimentale et clinique<\/h3>\n<p>Dans le champ de la validation psychom\u00e9trique des instruments de mesure, l&rsquo;\u00e9valuation de la dispersion via l&rsquo;IQR intervient pour contr\u00f4ler les biais de r\u00e9ponse extr\u00eames, tels que le biais de complaisance ou l&rsquo;acquiescement syst\u00e9matique. Lorsque des participants remplissent de longues batteries d&rsquo;\u00e9chelles de Likert sans engagement cognitif r\u00e9el, ils peuvent saturer les ancres extr\u00eames de l&rsquo;\u00e9chelle. L&rsquo;\u00e9valuation de l&rsquo;\u00e9cart interquartile des items permet alors de distinguer la variabilit\u00e9 r\u00e9elle des traits de personnalit\u00e9 sous-jacents de la dispersion artificielle g\u00e9n\u00e9r\u00e9e par des art\u00e9facts de passation.<\/p>\n<p>Au niveau des cohortes cliniques, l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 syndromique repr\u00e9sente un d\u00e9fi m\u00e9thodologique constant. Consid\u00e9rons par exemple une recherche \u00e9valuant l&rsquo;efficacit\u00e9 d&rsquo;une psychoth\u00e9rapie cognitive et comportementale sur une cohorte de patients pr\u00e9sentant un Trouble du D\u00e9ficit de l&rsquo;Attention avec ou sans Hyperactivit\u00e9 (TDAH). Les performances attentionnelles de ces sujets se caract\u00e9risent intrins\u00e8quement par une forte labilit\u00e9 intra- et interindividuelle. L&rsquo;utilisation de l&rsquo;\u00e9cart interquartile pour calibrer les scores de dispersion pr\u00e9- et post-intervention prot\u00e8ge le chercheur contre les biais inf\u00e9rentiels induits par les profils atypiques, garantissant que l&rsquo;estimation de l&rsquo;efficacit\u00e9 th\u00e9rapeutique repose sur le d\u00e9placement du groupe central plut\u00f4t que sur la normalisation fortuite d&rsquo;un patient singulier ayant r\u00e9pondu de mani\u00e8re atypique lors du pr\u00e9-test.<\/p>\n<p>De surcro\u00eet, les agences de r\u00e9gulation et les comit\u00e9s de lecture des revues scientifiques de premier rang exigent d\u00e9sormais une transparence m\u00e9thodologique accrue dans le traitement de la non-normalit\u00e9. L&rsquo;alignement sur les standards m\u00e9thodologiques contemporains implique d&rsquo;abandonner le recours automatique aux statistiques gaussiennes d\u00e8s lors que les postulats sous-jacents sont viol\u00e9s. L&rsquo;adoption de l&rsquo;IQR s&rsquo;inscrit directement dans cette d\u00e9marche de scientificit\u00e9 reproductible en fournissant une m\u00e9trique standardis\u00e9e, exempte de biais de manipulation, pour quantifier l&rsquo;incertitude et la variabilit\u00e9 comportementale.<\/p>\n<h2>2. Fondements math\u00e9matiques et m\u00e9thodes d&rsquo;interpolation des quartiles<\/h2>\n<h3>2.1 Calcul manuel et ordonnancement des observations<\/h3>\n<p>L&rsquo;estimation des quartiles commence imp\u00e9rativement par le tri s\u00e9quentiel croissant des observations collect\u00e9es. Soit un \u00e9chantillon de taille <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span>, not\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>. L&rsquo;ordonnancement math\u00e9matique produit la suite des statistiques d&rsquo;ordre not\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"x_{(1)} le x_{(2)} le dots le x_{(n)}\">$$x_{(1)} le x_{(2)} le dots le x_{(n)}$$<\/span><\/p>\n<p>Dans ce cadre ordonn\u00e9, les quartiles <span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"Q_2\">$Q_2$<\/span> (la m\u00e9diane) et <span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span> repr\u00e9sentent les points de coupure qui divisent la distribution en quatre segments contenant chacun un quart (<span class=\"arabpsy-math-inline\" data-latex=\"25%\">$25%$<\/span>) de l&rsquo;effectif total. Bien que conceptuellement limpide, la d\u00e9termination de l&rsquo;indice ou du rang empirique associ\u00e9 \u00e0 chaque quartile pr\u00e9sente une ambigu\u00eft\u00e9 arithm\u00e9tique d\u00e8s lors que le produit de la proportion vis\u00e9e <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>) par la taille de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> ne correspond pas \u00e0 un entier naturel.<\/p>\n<p>Il importe d&rsquo;\u00e9tablir ici une distinction conceptuelle nette entre quartiles, d\u00e9ciles et centiles au sein de l&rsquo;\u00e9talonnage psychom\u00e9trique. Si les centiles d\u00e9coupent la distribution en cent fractions \u00e9gales et les d\u00e9ciles en dix, les quartiles synth\u00e9tisent les transitions macro-structurelles. Consid\u00e9rons un \u00e9chantillon empirique restreint de <span class=\"arabpsy-math-inline\" data-latex=\"n = 9\">$n = 9$<\/span> scores observ\u00e9s sur une \u00e9chelle de phobie sociale :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S = {14, 17, 18, 20, 22, 25, 29, 31, 35}\">$$S = {14, 17, 18, 20, 22, 25, 29, 31, 35}$$<\/span><\/p>\n<p>Ces donn\u00e9es \u00e9tant ordonn\u00e9es de <span class=\"arabpsy-math-inline\" data-latex=\"x_{(1)} = 14\">$x_{(1)} = 14$<\/span> \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"x_{(9)} = 35\">$x_{(9)} = 35$<\/span>, la m\u00e9diane (<span class=\"arabpsy-math-inline\" data-latex=\"Q_2\">$Q_2$<\/span>) correspond directement \u00e0 la cinqui\u00e8me valeur, soit <span class=\"arabpsy-math-inline\" data-latex=\"x_{(5)} = 22\">$x_{(5)} = 22$<\/span>. Pour d\u00e9terminer <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>, plusieurs approches manuelles traditionnelles s&rsquo;affrontent : l&rsquo;approche inclusive (qui int\u00e8gre la m\u00e9diane dans les sous-ensembles inf\u00e9rieurs et sup\u00e9rieurs) et l&rsquo;approche exclusive (qui l&rsquo;exclut). Selon l&rsquo;approche exclusive de Tukey, le premier quartile est la m\u00e9diane de la moiti\u00e9 inf\u00e9rieure <span class=\"arabpsy-math-inline\" data-latex=\"{14, 17, 18, 20}\">${14, 17, 18, 20}$<\/span>, soit la moyenne arithm\u00e9tique entre $17$ et $18$ :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Q_1 = \\frac{17 + 18}{2} = 17{,}5\">$$Q_1 = \\frac{17 + 18}{2} = 17{,}5$$<\/span><\/p>\n<p>De m\u00eame, le troisi\u00e8me quartile correspond \u00e0 la m\u00e9diane de la moiti\u00e9 sup\u00e9rieure <span class=\"arabpsy-math-inline\" data-latex=\"{25, 29, 31, 35}\">${25, 29, 31, 35}$<\/span>, soit :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Q_3 = \\frac{29 + 31}{2} = 30{,}0\">$$Q_3 = \\frac{29 + 31}{2} = 30{,}0$$<\/span><\/p>\n<p>L&rsquo;\u00e9cart interquartile manuel selon cette m\u00e9thode \u00e9l\u00e9mentaire s&rsquo;\u00e9tablit donc \u00e0 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{IQR} = 30{,}0 - 17{,}5 = 12{,}5\">$$\\text{IQR} = 30{,}0 &#8211; 17{,}5 = 12{,}5$$<\/span><\/p>\n<h3>2.2 Probl\u00e9matique de la discontinuit\u00e9 et n\u00e9cessit\u00e9 de l&rsquo;interpolation<\/h3>\n<p>La simplicit\u00e9 apparente de l&rsquo;exemple pr\u00e9c\u00e9dent masque un probl\u00e8me th\u00e9orique fondamental de la statistique non param\u00e9trique : la discontinuit\u00e9 intrins\u00e8que des donn\u00e9es empiriques discr\u00e8tes compar\u00e9e \u00e0 la continuit\u00e9 suppos\u00e9e du continuum latent sous-jacent. D\u00e8s lors que <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> varie, ou que l&rsquo;on manipule des distributions pr\u00e9sentant des \u00e9galit\u00e9s de rangs (d\u00e9sign\u00e9es sous le terme technique d&rsquo;<em>ex \u00e6quo<\/em> ou <em>ties<\/em>), la d\u00e9finition des rangs fractionnaires devient tributaire d&rsquo;hypoth\u00e8ses g\u00e9om\u00e9triques sur la mani\u00e8re d&rsquo;interpoler l&rsquo;espace s\u00e9parant deux statistiques d&rsquo;ordre cons\u00e9cutives.<\/p>\n<p>Si l&rsquo;on cherche le centile <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> d&rsquo;une distribution, l&rsquo;indice th\u00e9orique du rang peut \u00eatre calcul\u00e9 selon diff\u00e9rentes formulations, telles que <span class=\"arabpsy-math-inline\" data-latex=\"k = p \\times n\">$k = p \\times n$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"k = p \\times (n + 1)\">$k = p \\times (n + 1)$<\/span>, ou encore <span class=\"arabpsy-math-inline\" data-latex=\"k = p \\times (n - 1) + 1\">$k = p \\times (n &#8211; 1) + 1$<\/span>. Lorsque cet indice <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> est d\u00e9compos\u00e9 en sa partie enti\u00e8re <span class=\"arabpsy-math-inline\" data-latex=\"lfloor k rfloor\">$lfloor k rfloor$<\/span> et sa partie fractionnaire <span class=\"arabpsy-math-inline\" data-latex=\"\\gamma = k - \\lfloor k \\rfloor\">$\\gamma = k &#8211; \\lfloor k \\rfloor$<\/span>, l&rsquo;estimation du quantile <span class=\"arabpsy-math-inline\" data-latex=\"Q(p)\">$Q(p)$<\/span> proc\u00e8de g\u00e9n\u00e9ralement par interpolation lin\u00e9aire entre deux valeurs adjacentes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Q(p) = (1 - \\gamma) x_{(\\lfloor k \\rfloor)} + \\gamma x_{(\\lfloor k \\rfloor + 1)}\">$$Q(p) = (1 &#8211; \\gamma) x_{(\\lfloor k \\rfloor)} + \\gamma x_{(\\lfloor k \\rfloor + 1)}$$<\/span><\/p>\n<p>Le choix de la pond\u00e9ration <span class=\"arabpsy-math-inline\" data-latex=\"\\gamma\">$\\gamma$<\/span> et de la fonction d&rsquo;indexation n&rsquo;est pas anodin. Dans l&rsquo;exercice clinique, notamment lors de bilans neuropsychologiques utilisant les \u00e9chelles d&rsquo;intelligence de Wechsler (comme la WAIS-IV) ou des tests de fluence verbale, une variation infinit\u00e9simale de la m\u00e9thode d&rsquo;interpolation peut d\u00e9placer un seuil diagnostique. Si le seuil d&rsquo;un d\u00e9ficit cognitif est d\u00e9fini par une performance inf\u00e9rieure au premier quartile (<span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span>) d&rsquo;une cohorte de r\u00e9f\u00e9rence, le basculement d&rsquo;une m\u00e9thode d&rsquo;interpolation continue \u00e0 une m\u00e9thode par paliers discrets est susceptible de requalifier un profil de la normalit\u00e9 vers la zone pathologique, illustrant la port\u00e9e clinique de ces subtilit\u00e9s math\u00e9matiques.<\/p>\n<h2>3. La fonction de base IQR() dans R : Syntaxe et param\u00e8tres essentiels<\/h2>\n<h3>3.1 Structure syntaxique g\u00e9n\u00e9rale et arguments de la fonction<\/h3>\n<p>L&rsquo;environnement statistique R impl\u00e9mente nativement le calcul de cette m\u00e9trique par l&rsquo;interm\u00e9diaire de la fonction <code>IQR()<\/code>, h\u00e9berg\u00e9e au sein du paquet fondamental <code>stats<\/code>, automatiquement charg\u00e9 \u00e0 l&rsquo;initialisation de chaque session de travail. La signature officielle de cette fonction se structure de la mani\u00e8re suivante :<\/p>\n<p><code>IQR(x, na.rm = FALSE, type = 7)<\/code><\/p>\n<p>L&rsquo;argument <code>x<\/code> constitue le vecteur num\u00e9rique primaire contenant les donn\u00e9es quantitatives dont l&rsquo;analyste souhaite extraire la dispersion. Ce vecteur doit imp\u00e9rativement poss\u00e9der un mode de stockage conforme aux types de donn\u00e9es num\u00e9riques de R, \u00e0 savoir <code>numeric<\/code>, <code>double<\/code> ou <code>integer<\/code>. L&rsquo;introduction accidentelle d&rsquo;une variable non quantitative d\u00e9clenche des comportements d&rsquo;\u00e9chec d&rsquo;ex\u00e9cution stricts.<\/p>\n<p>Une erreur fr\u00e9quente chez les chercheurs manipulant des donn\u00e9es import\u00e9es depuis des logiciels tiers r\u00e9side dans le passage involontaire de colonnes typ\u00e9es comme facteurs (<code>factor<\/code>) ou cha\u00eenes de caract\u00e8res (<code>character<\/code>). Cette situation survient r\u00e9guli\u00e8rement lorsque des codes de modalit\u00e9 textuelle (par exemple <code>\"Absence\"<\/code>) sont ins\u00e9r\u00e9s au sein d&rsquo;une variable num\u00e9rique de comptage de sympt\u00f4mes. L&rsquo;\u00e9valuation de l&rsquo;instruction <code>IQR(x)<\/code> sur un tel vecteur g\u00e9n\u00e8re invariablement l&rsquo;erreur classique :<\/p>\n<p><code>Error in x[!is.na(x)] : target of assignment expands to non-language object<\/code> ou <code>Error in quantile.default(as.numeric(x)) : 'x' must be numeric<\/code>.<\/p>\n<p>Le chercheur doit donc s&rsquo;assurer du typage pr\u00e9alable de ses variables \u00e0 l&rsquo;aide des fonctions diagnostiques <code>is.numeric(x)<\/code> et, le cas \u00e9ch\u00e9ant, proc\u00e9der \u00e0 des coercitions explicites via <code>as.numeric(as.character(x))<\/code> lorsque des facteurs num\u00e9riques doivent \u00eatre convertis en valeurs continues r\u00e9elles.<\/p>\n<h3>3.2 Le param\u00e8tre type et l&rsquo;alignement algorithmique<\/h3>\n<p>Le troisi\u00e8me argument de la fonction, <code>type = 7<\/code>, repr\u00e9sente l&rsquo;un des param\u00e8tres les plus cruciaux et pourtant les plus fr\u00e9quemment ignor\u00e9s par les praticiens de la recherche quantitative. R int\u00e8gre au sein de son moteur d&rsquo;interpolation pas moins de neuf algorithmes distincts de calcul des quantiles, document\u00e9s de fa\u00e7on magistrale dans l&rsquo;article de r\u00e9f\u00e9rence de Rob J. Hyndman et Yanan Fan publi\u00e9 en 1996. Le comportement par d\u00e9faut de R consiste \u00e0 assigner la valeur <code>type = 7<\/code>.<\/p>\n<p>Cette d\u00e9cision architecturale distingue imm\u00e9diatement R d&rsquo;autres logiciels propri\u00e9taires omnipr\u00e9sents dans les laboratoires de psychologie. Notamment, le logiciel commercial SPSS (IBM Corp.) ainsi que le progiciel SAS emploient par d\u00e9faut un algorithme d&rsquo;interpolation lin\u00e9aire correspondant au <code>type = 6<\/code> de R. En cons\u00e9quence, si un chercheur calcule l&rsquo;\u00e9cart interquartile sur un m\u00eame jeu de donn\u00e9es cliniques sans modifier les param\u00e8tres par d\u00e9faut, les r\u00e9sultats num\u00e9riques fournis par la console R et par la fen\u00eatre de sortie de SPSS afficheront une divergence syst\u00e9matique sur les petits et moyens \u00e9chantillons.<\/p>\n<p>Cette disparit\u00e9 algorithmique met en p\u00e9ril direct la r\u00e9plicabilit\u00e9 exacte des protocoles exp\u00e9rimentaux. Pour assurer une coh\u00e9rence stricte lors d&rsquo;analyses crois\u00e9es ou lors de la soumission d&rsquo;\u00e9tudes multicentriques collaboratives, l&rsquo;analyste doit expliciter syst\u00e9matiquement le param\u00e8tre <code>type<\/code> employ\u00e9 au sein de son script R. L&rsquo;alignement rigoureux sur les standards de calcul retenus par l&rsquo;ensemble des partenaires d&rsquo;un consortium scientifique conditionne l&rsquo;exactitude des comparaisons inter-laboratoires.<\/p>\n<h2>4. Exemple 1 : Calcul de l&rsquo;\u00e9cart interquartile sur un vecteur simple de donn\u00e9es<\/h2>\n<h3>4.1 Cr\u00e9ation du vecteur et ex\u00e9cution de l&rsquo;instruction de base<\/h3>\n<p>Afin d&rsquo;illustrer la m\u00e9canique fondamentale de la fonction, simulons un vecteur repr\u00e9sentatif des scores d&rsquo;anxi\u00e9t\u00e9-\u00e9tat collect\u00e9s aupr\u00e8s de douze \u00e9tudiants universitaires avant le passage d&rsquo;une \u00e9preuve d&rsquo;\u00e9valuation psychologique standardis\u00e9e. Les scores sont mesur\u00e9s sur une \u00e9chelle continue gradu\u00e9e de 10 \u00e0 50 points :<\/p>\n<p><code>scores_anxiete &lt;- c(18, 21, 22, 24, 25, 26, 29, 32, 35, 41, 44, 48)<\/code><\/p>\n<p>L&rsquo;\u00e9valuation directe de la commande dans l&rsquo;interpr\u00e9teur R s&rsquo;effectue sans param\u00e8tre additionnel :<\/p>\n<p><code>iqr_anxiete &lt;- IQR(scores_anxiete)<\/code><br \/>\n<code>print(iqr_anxiete)<\/code><\/p>\n<p>La console renvoie imm\u00e9diatement la valeur <code>13.75<\/code>. Sur le plan de l&rsquo;interpr\u00e9tation psychologique, ce r\u00e9sultat indique formellement que l&rsquo;amplitude de l&rsquo;\u00e9chelle d&rsquo;anxi\u00e9t\u00e9 occup\u00e9e par les 50 % d&rsquo;\u00e9tudiants se situant au centre de la distribution s&rsquo;\u00e9l\u00e8ve \u00e0 13,75 points. Les participants centraux ne s&rsquo;\u00e9talent pas sur l&rsquo;ensemble de la gamme th\u00e9orique de l&rsquo;instrument, mais pr\u00e9sentent une variabilit\u00e9 mod\u00e9r\u00e9e autour du niveau m\u00e9dian.<\/p>\n<p>Pour valider l&rsquo;int\u00e9grit\u00e9 de cette mesure, l&rsquo;analyste recourt traditionnellement \u00e0 la fonction <code>summary(scores_anxiete)<\/code>, laquelle renvoie un tableau synth\u00e9tique incluant le minimum, le premier quartile, la m\u00e9diane, la moyenne, le troisi\u00e8me quartile et le maximum. L&rsquo;ex\u00e9cution de cette commande fournit les rep\u00e8res suivants :<\/p>\n<p><code>Min. : 18.00<\/code><br \/>\n<code>1st Qu. : 23.50<\/code><br \/>\n<code>Median : 27.50<\/code><br \/>\n<code>Mean : 30.42<\/code><br \/>\n<code>3rd Qu. : 37.25<\/code><br \/>\n<code>Max. : 48.00<\/code><\/p>\n<p>Nous observons sans ambigu\u00eft\u00e9 que la diff\u00e9rence directe entre le troisi\u00e8me quartile (<span class=\"arabpsy-math-inline\" data-latex=\"37{,}25\">$37{,}25$<\/span>) et le premier quartile (<span class=\"arabpsy-math-inline\" data-latex=\"23{,}50\">$23{,}50$<\/span>) correspond rigoureusement \u00e0 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"37{,}25 - 23{,}50 = 13{,}75\">$$37{,}25 &#8211; 23{,}50 = 13{,}75$$<\/span><\/p>\n<h3>4.2 V\u00e9rification manuelle via la fonction quantile()<\/h3>\n<p>Pour ma\u00eetriser pleinement la machinerie analytique sous-jacente, il est fortement recommand\u00e9 aux chercheurs de d\u00e9composer l&rsquo;appel fonctionnel d&rsquo;<code>IQR()<\/code> en exploitant la primitive matricielle sous-jacente : la fonction <code>quantile()<\/code>. En effet, la fonction <code>IQR()<\/code> n&rsquo;est qu&rsquo;une surcouche syntaxique (un <em>wrapper<\/em>) optimis\u00e9e encapsulant la commande <code>diff(quantile(x, probs = c(0.25, 0.75)))<\/code>.<\/p>\n<p>Nous pouvons expliciter cette extraction quartilaire cibl\u00e9e en d\u00e9finissant les probabilit\u00e9s cumul\u00e9es vectorielles dans l&rsquo;argument <code>probs<\/code> :<\/p>\n<p><code>quartiles_extraits &lt;- quantile(scores_anxiete, probs = c(0.25, 0.75), type = 7)<\/code><br \/>\n<code>print(quartiles_extraits)<\/code><\/p>\n<p>La console renvoie un vecteur nomm\u00e9 contenant les pourcentages correspondants :<\/p>\n<p><code>  25%   75% <\/code><br \/>\n<code>23.50 37.25<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de la diff\u00e9rence arithm\u00e9tique explicite confirme l&rsquo;ad\u00e9quation logique :<\/p>\n<p><code>difference_manuelle &lt;- quartiles_extraits[2] - quartiles_extraits[1]<\/code><br \/>\n<code>names(difference_manuelle) &lt;- NULL<\/code><br \/>\n<code>print(difference_manuelle)<\/code><\/p>\n<p>Ce calcul d\u00e9compos\u00e9 retourne exactement <code>13.75<\/code>. L&rsquo;insertion de ces \u00e9tapes explicites au sein de scripts document\u00e9s favorise une auditabilit\u00e9 compl\u00e8te des donn\u00e9es de laboratoire, garantissant que chaque membre d&rsquo;une \u00e9quipe de recherche puisse inspecter visuellement les points d&rsquo;ancrage de la dispersion sans d\u00e9pendre d&rsquo;abstractions implicites.<\/p>\n<h2>5. Exemple 2 : Gestion rigoureuse des donn\u00e9es manquantes (NA)<\/h2>\n<h3>5.1 Comportement par d\u00e9faut face aux valeurs NA dans les questionnaires psychom\u00e9triques<\/h3>\n<p>Dans la pratique r\u00e9elle de la recherche quantitative en sciences humaines, les jeux de donn\u00e9es parfaits et int\u00e8gres constituent une exception rarissime. Lors de la passation de questionnaires auto-rapport\u00e9s ou d&rsquo;inventaires d&rsquo;\u00e9valuation comportementale, les chercheurs sont constamment confront\u00e9s \u00e0 des donn\u00e9es manquantes, mod\u00e9lis\u00e9es dans le langage R sous la constante logique <code>NA<\/code> (signifiant <em>Not Available<\/em>). Les causes de ces absences d&rsquo;information sont multiples : omissions involontaires d&rsquo;un participant tournant deux pages simultan\u00e9ment, refus explicite de r\u00e9pondre \u00e0 une question relative \u00e0 un trauma personnel, ou interruption technique lors d&rsquo;une transmission de donn\u00e9es via une plateforme en ligne.<\/p>\n<p>Consid\u00e9rons l&rsquo;insertion d\u00e9lib\u00e9r\u00e9e d&rsquo;une donn\u00e9e manquante au sein de notre vecteur d&rsquo;anxi\u00e9t\u00e9 initial :<\/p>\n<p><code>scores_lacunaires &lt;- c(18, 21, NA, 24, 25, 26, 29, 32, 35, 41, 44, 48)<\/code><\/p>\n<p>Si l&rsquo;analyste ex\u00e9cute naivement la fonction par d\u00e9faut :<\/p>\n<p><code>IQR(scores_lacunaires)<\/code><\/p>\n<p>Le syst\u00e8me renvoie syst\u00e9matiquement et inexorablement la valeur <code>NA<\/code>. Ce comportement, bien que source de frustration fr\u00e9quente pour les d\u00e9butants, d\u00e9coule d&rsquo;une rigueur philosophique propre \u00e0 la logique computationnelle de R. Par d\u00e9faut, si un \u00e9l\u00e9ment demeure inconnu dans une s\u00e9rie math\u00e9matique, la valeur r\u00e9elle de l&rsquo;\u00e9tendue des 50 % centraux de cette m\u00eame s\u00e9rie devient formellement ind\u00e9termin\u00e9e. R refuse d&rsquo;\u00e9mettre une hypoth\u00e8se implicite sur la valeur de la donn\u00e9e absente.<\/p>\n<p>D&rsquo;un point de vue m\u00e9thodologique, cette neutralit\u00e9 algorithmique contraint l&rsquo;exp\u00e9rimentateur \u00e0 s&rsquo;interroger sur la nature \u00e9pist\u00e9mologique de l&rsquo;absence selon la typologie de Donald Rubin : les donn\u00e9es sont-elles Manquantes Compl\u00e8tement au Hasard (<a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Missing_data#Missing_completely_at_random\">Missing Completely at Random ou MCAR<\/a>), Manquantes au Hasard (MAR), ou Manquantes de Fa\u00e7on Non-Al\u00e9atoire (Missing Not at Random ou MNAR) ? L&rsquo;ignorance aveugle de cette classification peut induire des biais d&rsquo;\u00e9chantillonnage massifs lors de la quantification de la dispersion.<\/p>\n<h3>5.2 Impl\u00e9mentation du param\u00e8tre na.rm = TRUE<\/h3>\n<p>Pour forcer le calcul de l&rsquo;\u00e9cart interquartile sur les observations effectivement disponibles en \u00e9cartant de mani\u00e8re cibl\u00e9e les valeurs manquantes, il convient d&rsquo;activer l&rsquo;argument logique <code>na.rm = TRUE<\/code> (d\u00e9signant <em>NA removal<\/em>) :<\/p>\n<p><code>iqr_nettoye &lt;- IQR(scores_lacunaires, na.rm = TRUE)<\/code><br \/>\n<code>print(iqr_nettoye)<\/code><\/p>\n<p>Sous cette instruction, R commence par appliquer un filtre logique vectoriel en interne, \u00e9liminant tous les \u00e9l\u00e9ments pour lesquels <code>is.na(x)<\/code> retourne la valeur <code>TRUE<\/code>. Le calcul s&rsquo;ex\u00e9cute ensuite sur le sous-vecteur r\u00e9duit compos\u00e9 de <span class=\"arabpsy-math-inline\" data-latex=\"n = 11\">$n = 11$<\/span> participants effectifs, renvoyant un nouvel \u00e9cart interquartile recalcul\u00e9 sur ces donn\u00e9es r\u00e9siduelles : <code>14.5<\/code>.<\/p>\n<p>Cette op\u00e9ration m\u00e9canique d&rsquo;omission n&rsquo;est pas sans cons\u00e9quences m\u00e9thodologiques. L&rsquo;exclusion unilat\u00e9rale d&rsquo;observations modifie la taille r\u00e9elle de l&rsquo;\u00e9chantillon disponible et alt\u00e8re la puissance statistique des analyses futures. Avant toute application g\u00e9n\u00e9ralis\u00e9e de l&rsquo;argument <code>na.rm = TRUE<\/code>, les bonnes pratiques de la recherche exp\u00e9rimentale imposent de quantifier et de consigner explicitement le taux d&rsquo;attrition au sein des rapports de laboratoire :<\/p>\n<p><code>taux_manquants &lt;- mean(is.na(scores_lacunaires)) * 100<\/code><br \/>\n<code>cat(\"Pourcentage de donn\u00e9es manquantes :\", round(taux_manquants, 2), \"%n\")<\/code><\/p>\n<p>Si le taux d&rsquo;attrition d\u00e9passe un seuil critique (commun\u00e9ment fix\u00e9 \u00e0 5 % ou 10 % selon les paradigmes), le chercheur doit obligatoirement justifier l&rsquo;usage de la suppression par liste ou envisager des approches d&rsquo;imputation multiple adapt\u00e9es aux donn\u00e9es non gaussiennes (telles que l&rsquo;imputation pr\u00e9dictive par les plus proches voisins ou via des for\u00eats al\u00e9atoires avec le paquet <code>mice<\/code>) plut\u00f4t que d&rsquo;appliquer silencieusement une troncature univari\u00e9e.<\/p>\n<h2>6. Exemple 3 : Calcul de l&rsquo;IQR au sein de data frames structur\u00e9s<\/h2>\n<h3>6.1 Acc\u00e8s aux colonnes et utilisation de l&rsquo;indexation de base R<\/h3>\n<p>Dans un contexte exp\u00e9rimental r\u00e9aliste, les variables psychom\u00e9triques ne sont que rarement manipul\u00e9es sous la forme de vecteurs isol\u00e9s flottant dans l&rsquo;environnement de travail. Elles sont structur\u00e9es sous forme de matrices tabulaires ou de cadres de donn\u00e9es (les fameux <code>data.frame<\/code> ou <code>tibble<\/code>) regroupant plusieurs dizaines de variables op\u00e9rationnelles pour des centaines de participants. G\u00e9n\u00e9rons une structure de donn\u00e9es typique issue d&rsquo;un protocole d&rsquo;\u00e9valuation clinique :<\/p>\n<p><code>donnees_cliniques &lt;- data.frame(<br \/>\n  identifiant = paste0(\"SUJET_\", 101:108),<br \/>\n  score_depressif = c(14, 28, 19, 42, 22, 16, 31, 25),<br \/>\n  latence_attention = c(420, 510, 480, 890, 460, 430, 620, 540),<br \/>\n  statut_clinique = factor(c(\"Temoin\", \"Patient\", \"Temoin\", \"Patient\", \"Temoin\", \"Temoin\", \"Patient\", \"Patient\"))<br \/>\n)<\/code><\/p>\n<p>Pour calculer l&rsquo;\u00e9cart interquartile sur une colonne sp\u00e9cifique en utilisant la syntaxe historique de base de R, l&rsquo;op\u00e9rateur dollar <code>$<\/code> repr\u00e9sente le moyen d&rsquo;indexation le plus direct :<\/p>\n<p><code>iqr_depression &lt;- IQR(donnees_cliniques$score_depressif, na.rm = TRUE)<\/code><br \/>\n<code>print(iqr_depression)<\/code><\/p>\n<p>Lorsque le protocole de recherche requiert le calcul simultan\u00e9 de l&rsquo;\u00e9cart interquartile \u00e0 travers une multitude de colonnes continues, l&rsquo;\u00e9valuation manuelle instruction par instruction devient inop\u00e9rante et source d&rsquo;erreurs. Il est alors judicieux d&rsquo;exploiter la famille des fonctions de vectorisation de base, plus particuli\u00e8rement <code>sapply()<\/code> ou <code>lapply()<\/code> :<\/p>\n<p><code>colonnes_cibles &lt;- donnees_cliniques[, c(\"score_depressif\", \"latence_attention\")]<\/code><br \/>\n<code>iqr_multiples &lt;- sapply(colonnes_cibles, IQR, na.rm = TRUE)<\/code><br \/>\n<code>print(iqr_multiples)<\/code><\/p>\n<p>Cette commande it\u00e8re sur chaque colonne s\u00e9lectionn\u00e9e et renvoie un vecteur nomm\u00e9 clair indiquant imm\u00e9diatement la dispersion de l&rsquo;\u00e9chelle de d\u00e9pression (<span class=\"arabpsy-math-inline\" data-latex=\"10{,}25\">$10{,}25$<\/span>) et celle de la latence attentionnelle (<span class=\"arabpsy-math-inline\" data-latex=\"127{,}5\">$127{,}5$<\/span> millisecondes).<\/p>\n<h3>6.2 Approche moderne et fluide avec le Tidyverse (dplyr)<\/h3>\n<p>Bien que les commandes de base de R demeurent fonctionnelles, l&rsquo;adoption pr\u00e9dominante du m\u00e9tapaquet <a rel=\"noopener noreferrer\" href=\"https:\/\/www.tidyverse.org\/\">Tidyverse<\/a> et plus singuli\u00e8rement de l&rsquo;extension <code>dplyr<\/code> offre une lisibilit\u00e9 accrue et une int\u00e9gration directe dans des pipelines de traitement analytique complexes. L&rsquo;approche repose sur l&rsquo;encha\u00eenement logique d&rsquo;instructions via l&rsquo;op\u00e9rateur de tuyau (le <em>pipe<\/em> natif de R <code>|&gt;<\/code> ou le pipe historique de magrittr <code>%&gt;%<\/code>).<\/p>\n<p>Pour extraire une synth\u00e8se robuste au moyen du verbe <code>summarise()<\/code>, le pipeline s&rsquo;agence comme suit :<\/p>\n<p><code>library(dplyr)<\/code><\/p>\n<p><code>synthese_robuste &lt;- donnees_cliniques |&gt;<br \/>\n  summarise(<br \/>\n    iqr_dep = IQR(score_depressif, na.rm = TRUE),<br \/>\n    mediane_dep = median(score_depressif, na.rm = TRUE),<br \/>\n    iqr_latence = IQR(latence_attention, na.rm = TRUE),<br \/>\n    mediane_latence = median(latence_attention, na.rm = TRUE),<br \/>\n    n_effectif = n()<br \/>\n  )<\/code><br \/>\n<code>print(synthese_robuste)<\/code><\/p>\n<p>Cette approche prend tout son sens lors du traitement de collectes massives de donn\u00e9es en ligne h\u00e9berg\u00e9es sur des plateformes comme Qualtrics ou Prolific, o\u00f9 des centaines de variables brutes doivent \u00eatre agr\u00e9g\u00e9es. Gr\u00e2ce \u00e0 la fonction auxiliaire <code>across()<\/code>, il est possible d&rsquo;automatiser le calcul de l&rsquo;IQR sur toutes les variables num\u00e9riques du tableau en une unique formulation \u00e9l\u00e9gante :<\/p>\n<p><code>synthese_globale &lt;- donnees_cliniques |&gt;<br \/>\n  summarise(across(where(is.numeric), list(iqr = ~IQR(.x, na.rm = TRUE), med = ~median(.x, na.rm = TRUE))))<\/code><br \/>\n<code>print(synthese_globale)<\/code><\/p>\n<p>Cette flexibilit\u00e9 programmatique limite la redondance du code (principe <em>DRY: Don&rsquo;t Repeat Yourself<\/em>) et r\u00e9duit drastiquement les risques d&rsquo;omission lors de l&rsquo;exploration de r\u00e9pertoires volumineux de donn\u00e9es comportementales.<\/p>\n<h2>7. Approfondissement des 9 algorithmes de calcul des quantiles (argument type)<\/h2>\n<h3>7.1 Classification des m\u00e9thodes d&rsquo;estimation disponibles dans R<\/h3>\n<p>La taxonomie \u00e9tablie par Hyndman et Fan en 1996 structure les neuf algorithmes disponibles au sein de R en deux grandes familles math\u00e9matiques : les estimateurs non continus ou discrets (types 1 \u00e0 3) et les estimateurs continus fond\u00e9s sur une interpolation lin\u00e9aire par morceaux (types 4 \u00e0 9). Comprendre cette typologie s&rsquo;av\u00e8re indispensable pour naviguer entre les diff\u00e9rents logiciels de traitement statistique employ\u00e9s en recherche exp\u00e9rimentale.<\/p>\n<p>Les m\u00e9thodes non continues s&rsquo;appliquent originellement aux distributions empiriques pures :<\/p>\n<ul>\n<li><strong>Type 1 :<\/strong> Fonction empirique inverse par paliers discrets sans interpolation (standard historique des quantiles empiriques).<\/li>\n<li><strong>Type 2 :<\/strong> Fonction empirique invers\u00e9e avec interpolation par palier aux points de discontinuit\u00e9 (moyenne arithm\u00e9tique aux points d&rsquo;\u00e9galit\u00e9).<\/li>\n<li><strong>Type 3 :<\/strong> Algorithme des plus proches entiers, historiquement impl\u00e9ment\u00e9 dans le langage SAS pour l&rsquo;estimation de quantiles discrets.<\/li>\n<\/ul>\n<p>Les m\u00e9thodes continues interpolent lin\u00e9airement les rangs th\u00e9oriques pour d\u00e9river une estimation continue :<\/p>\n<ul>\n<li><strong>Type 4 :<\/strong> Interpolation lin\u00e9aire directe bas\u00e9e sur les probabilit\u00e9s cumul\u00e9es <span class=\"arabpsy-math-inline\" data-latex=\"p(k) = k \/ n\">$p(k) = k \/ n$<\/span>.<\/li>\n<li><strong>Type 5 :<\/strong> Formule populaire en hydrologie et m\u00e9t\u00e9orologie avec <span class=\"arabpsy-math-inline\" data-latex=\"p(k) = (k - 0{,}5) \/ n\">$p(k) = (k &#8211; 0{,}5) \/ n$<\/span>, fournissant un estimateur sans biais sous certaines distributions sym\u00e9triques.<\/li>\n<li><strong>Type 6 :<\/strong> Utilise la formulation de positionnement <span class=\"arabpsy-math-inline\" data-latex=\"p(k) = k \/ (n + 1)\">$p(k) = k \/ (n + 1)$<\/span>. C&rsquo;est le standard immuable de logiciels d&rsquo;ing\u00e9nierie et de statistique r\u00e9put\u00e9s tels que <em>Minitab<\/em>, <em>SPSS<\/em> et les modules de base de <em>SAS<\/em>.<\/li>\n<li><strong>Type 7 :<\/strong> Utilise la relation de positionnement <span class=\"arabpsy-math-inline\" data-latex=\"p(k) = (k - 1) \/ (n - 1)\">$p(k) = (k &#8211; 1) \/ (n &#8211; 1)$<\/span>. Il constitue le param\u00e8tre <strong>par d\u00e9faut sous R<\/strong>, l&rsquo;ancien langage <em>S<\/em>, ainsi que le standard des biblioth\u00e8ques scientifiques <em>Python (NumPy, SciPy)<\/em>.<\/li>\n<li><strong>Type 8 :<\/strong> Repose sur les recommandations de rang m\u00e9dian non biais\u00e9 issues de la statistique d\u00e9cisionnelle : <span class=\"arabpsy-math-inline\" data-latex=\"p(k) = (k - 1\/3) \/ (n + 1\/3)\">$p(k) = (k &#8211; 1\/3) \/ (n + 1\/3)$<\/span>.<\/li>\n<li><strong>Type 9 :<\/strong> Fond\u00e9 sur l&rsquo;approximation de Blom pour les scores normaux : <span class=\"arabpsy-math-inline\" data-latex=\"p(k) = (k - 3\/8) \/ (n + 1\/4)\">$p(k) = (k &#8211; 3\/8) \/ (n + 1\/4)$<\/span>.<\/li>\n<\/ul>\n<p>L&rsquo;existence de cette multiplicit\u00e9 d\u00e9montre qu&rsquo;il n&rsquo;existe pas d&rsquo;estimateur de quantile uniform\u00e9ment sup\u00e9rieur dans toutes les conditions. Toutefois, la divergence structurelle entre le <code>type = 7<\/code> (R) et le <code>type = 6<\/code> (SPSS) est la cause premi\u00e8re des incoh\u00e9rences de calcul constat\u00e9es par les chercheurs passant d&rsquo;un logiciel \u00e0 l&rsquo;autre.<\/p>\n<h3>7.2 Sensibilit\u00e9 des estimations selon la taille de l&rsquo;\u00e9chantillon psychologique<\/h3>\n<p>L&rsquo;impact de la s\u00e9lection de l&rsquo;algorithme d&rsquo;interpolation est inversement proportionnel \u00e0 la taille de l&rsquo;\u00e9chantillon analys\u00e9. Dans des cohortes \u00e9pid\u00e9miologiques massives ou des \u00e9tudes en g\u00e9n\u00e9tique comportementale comprenant des milliers d&rsquo;enregistrements (<span class=\"arabpsy-math-inline\" data-latex=\"n &gt; 1000\">$n > 1000$<\/span>), les diff\u00e9rences observ\u00e9es entre les neuf algorithmes deviennent totalement n\u00e9gligeables en vertu du th\u00e9or\u00e8me de Glivenko-Cantelli et de la convergence asymptotique de la fonction de distribution empirique vers la fonction de r\u00e9partition th\u00e9orique.<\/p>\n<p>\u00c0 l&rsquo;oppos\u00e9, dans les \u00e9tudes cliniques neuropsychologiques ou les protocoles de neuro-imagerie fonctionnelle o\u00f9 la taille d&rsquo;\u00e9chantillon est restreinte en raison des co\u00fbts logistiques ($n < 30$), le choix du type algorithmique produit des variations sensibles sur la valeur absolue de l'IQR. Illustrons cette sensibilit\u00e9 par une simulation concr\u00e8te sous R :<\/p>\n<p><code>echantillon_restreint &lt;- c(10, 12, 15, 18, 22, 27, 35, 42)<\/code><br \/>\n<code>comparaison_types &lt;- sapply(1:9, function(m) IQR(echantillon_restreint, type = m))<\/code><br \/>\n<code>names(comparaison_types) &lt;- paste0(\"Type_\", 1:9)<\/code><br \/>\n<code>print(comparaison_types)<\/code><\/p>\n<p>L&rsquo;inspection des r\u00e9sultats montre des \u00e9carts flagrants : le <code>Type_1<\/code> renvoie un IQR de <code>17<\/code>, le standard R <code>Type_7<\/code> fournit <code>17.25<\/code>, tandis que le standard SPSS <code>Type_6<\/code> produit <code>19.25<\/code>. Une diff\u00e9rence de deux points entiers sur une \u00e9chelle restreinte repr\u00e9sente une variation relative non n\u00e9gligeable. En cons\u00e9quence, les chercheurs manipulant de petites cohortes doivent documenter explicitement le param\u00e8tre <code>type<\/code> afin de pr\u00e9server l&rsquo;interop\u00e9rabilit\u00e9 et la fid\u00e9lit\u00e9 de leurs conclusions.<\/p>\n<h2>8. Calcul comparatif de l&rsquo;IQR par groupes exp\u00e9rimentaux<\/h2>\n<h3>8.1 Segmentation par conditions avec tapply() et aggregate()<\/h3>\n<p>Dans la recherche exp\u00e9rimentale en psychologie, l&rsquo;objectif fondamental consiste rarement \u00e0 \u00e9valuer une variable \u00e0 l&rsquo;\u00e9chelle globale ; il s&rsquo;agit plut\u00f4t de contraster des conditions exp\u00e9rimentales distinctes (par exemple, Groupe Exp\u00e9rimental recevant un protocole de rem\u00e9diation cognitive versus Groupe T\u00e9moin plac\u00e9 sous condition d&rsquo;attente). Segmenter la dispersion pour chaque strate est donc une n\u00e9cessit\u00e9 routini\u00e8re.<\/p>\n<p>R de base met \u00e0 disposition la fonction historique <code>tapply()<\/code>, con\u00e7ue pour appliquer une fonction scalaire sur un vecteur num\u00e9rique ventil\u00e9 par les facteurs d&rsquo;un tableau indic\u00e9 :<\/p>\n<p><code>donnees_experience &lt;- data.frame(<br \/>\n  score_vigilance = c(85, 88, 72, 60, 45, 90, 84, 82, 55, 40, 92, 87, 80, 50, 42),<br \/>\n  condition = factor(rep(c(\"Controle\", \"Privation_Sommeil_Moderee\", \"Privation_Sommeil_Severe\"), each = 5))<br \/>\n)<\/code><\/p>\n<p><code>iqr_par_condition &lt;- tapply(donnees_experience<span class=\"arabpsy-math-inline\" data-latex=\"score_vigilance, donnees_experience\">$score_vigilance, donnees_experience$<\/span>condition, IQR, na.rm = TRUE)<\/code><br \/>\n<code>print(iqr_par_condition)<\/code><\/p>\n<p>Cette commande synth\u00e9tise instantan\u00e9ment l&rsquo;\u00e9talement de la vigilance au sein de chaque niveau du facteur exp\u00e9rimental. Si l&rsquo;analyste d\u00e9sire extraire cette information dans une structure de donn\u00e9es tabulaire formelle pr\u00eate pour une inclusion dans un pipeline de transformation, la fonction <code>aggregate()<\/code> constitue une alternative particuli\u00e8rement lisible :<\/p>\n<p><code>tableau_recapitulatif &lt;- aggregate(score_vigilance ~ condition, data = donnees_experience, FUN = function(x) c(IQR = IQR(x), Mediane = median(x)))<\/code><br \/>\n<code>print(tableau_recapitulatif)<\/code><\/p>\n<p>Cette syntaxe \u00e0 base de formules (<code>y ~ x<\/code>) conf\u00e8re au code une structure quasi-conceptuelle analogue aux mod\u00e9lisations inf\u00e9rentielles (comme les mod\u00e8les lin\u00e9aires g\u00e9n\u00e9alogiques), ce qui facilite sa m\u00e9morisation et son d\u00e9ploiement rapide lors des phases pr\u00e9liminaires d&rsquo;audit de donn\u00e9es.<\/p>\n<h3>8.2 Analyse stratifi\u00e9e multifactorielle avec group_by() de dplyr<\/h3>\n<p>D\u00e8s lors que le protocole exp\u00e9rimental s&rsquo;enrichit et adopte un plan factoriel crois\u00e9 (par exemple un plan factoriel <span class=\"arabpsy-math-inline\" data-latex=\"2 \\times 2\">$2 \\times 2$<\/span> \u00e9valuant simultan\u00e9ment le Type de Traitement et le Sexe Biologique des participants), les outils de base peuvent rapidement conduire \u00e0 une syntaxe verbeuse. L&rsquo;utilisation conjointe des fonctions <code>group_by()<\/code> et <code>summarise()<\/code> issues de l&rsquo;extension <code>dplyr<\/code> offre une fluidit\u00e9 d&rsquo;analyse optimale :<\/p>\n<p><code>donnees_croisees &lt;- data.frame(<br \/>\n  anxiete = c(12, 14, 15, 20, 19, 24, 28, 30, 18, 22, 25, 27, 32, 35, 40, 42),<br \/>\n  traitement = factor(rep(c(\"Placebo\", \"Molecule_Active\"), each = 8)),<br \/>\n  sexe = factor(rep(c(\"Femme\", \"Homme\"), times = 8))<br \/>\n)<\/code><\/p>\n<p><code>synthese_multifactorielle &lt;- donnees_croisees |&gt;<br \/>\n  group_by(traitement, sexe) |&gt;<br \/>\n  summarise(<br \/>\n    iqr_anxiete = IQR(anxiete, na.rm = TRUE),<br \/>\n    mediane_anxiete = median(anxiete, na.rm = TRUE),<br \/>\n    taille_groupe = n(),<br \/>\n    .groups = \"drop\"<br \/>\n  )<\/code><br \/>\n<code>print(synthese_multifactorielle)<\/code><\/p>\n<p>L&rsquo;argument <code>.groups = \"drop\"<\/code> garantit la d\u00e9sactivation explicite du regroupement dans la table r\u00e9sultante, \u00e9vitant ainsi des erreurs d&rsquo;agr\u00e9gation non d\u00e9sir\u00e9es lors des \u00e9tapes ult\u00e9rieures de transformation de donn\u00e9es. La table produite peut ensuite \u00eatre directement export\u00e9e vers des formats de publication acad\u00e9mique (tels que des tableaux LaTeX ou des documents Word formatt\u00e9s via les biblioth\u00e8ques <code>knitr::kable()<\/code>, <code>gt<\/code> ou <code>flextable<\/code>), assurant un flux de travail continu de la donn\u00e9e brute au manuscrit final.<\/p>\n<h2>9. Identification et traitement des valeurs aberrantes via la m\u00e9thode de Tukey<\/h2>\n<h3>9.1 R\u00e8gle op\u00e9rationnelle des cl\u00f4tures de Tukey (Tukey&rsquo;s Fences)<\/h3>\n<p>Au-del\u00e0 de sa fonction de descripteur univari\u00e9 de la dispersion, l&rsquo;\u00e9cart interquartile constitue la pierre angulaire de l&rsquo;une des m\u00e9thodes les plus c\u00e9l\u00e8bres pour l&rsquo;identification objective des valeurs aberrantes (ou <em>outliers<\/em>) : la m\u00e9thode des cl\u00f4tures de Tukey. Face \u00e0 l&rsquo;arbitraire consistant \u00e0 \u00e9carter des observations sur de simples jugements subjectifs, John Tukey a propos\u00e9 une r\u00e8gle g\u00e9om\u00e9trique formelle reposant exclusivement sur la position des quartiles et l&rsquo;amplitude de l&rsquo;\u00e9cart interquartile.<\/p>\n<p>Les fronti\u00e8res d\u00e9limitant les observations acceptables, d\u00e9sign\u00e9es sous le terme de cl\u00f4tures (<em>fences<\/em>), sont d\u00e9finies math\u00e9matiquement par les expressions suivantes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Cl\u00f4ture Inf\u00e9rieure} = Q_1 - (k \\times \\text{IQR})\">$$\\text{Cl\u00f4ture Inf\u00e9rieure} = Q_1 &#8211; (k \\times \\text{IQR})$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Cl\u00f4ture Sup\u00e9rieure} = Q_3 + (k \\times \\text{IQR})\">$$\\text{Cl\u00f4ture Sup\u00e9rieure} = Q_3 + (k \\times \\text{IQR})$$<\/span><\/p>\n<p>Le coefficient multiplicateur <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> d\u00e9termine la s\u00e9v\u00e9rit\u00e9 du crit\u00e8re de d\u00e9tection :<\/p>\n<ul>\n<li>Lorsque <span class=\"arabpsy-math-inline\" data-latex=\"k = 1{,}5\">$k = 1{,}5$<\/span>, la proc\u00e9dure isole les <strong>valeurs aberrantes mod\u00e9r\u00e9es<\/strong>. Dans une distribution gaussienne th\u00e9orique, cette zone englobe environ <span class=\"arabpsy-math-inline\" data-latex=\"99{,}3%\">$99{,}3%$<\/span> des donn\u00e9es, rel\u00e9guant approximativement <span class=\"arabpsy-math-inline\" data-latex=\"0{,}7%\">$0{,}7%$<\/span> des scores au statut d&rsquo;anomalies potentielles.<\/li>\n<li>Lorsque <span class=\"arabpsy-math-inline\" data-latex=\"k = 3{,}0\">$k = 3{,}0$<\/span>, la proc\u00e9dure identifie les <strong>valeurs aberrantes extr\u00eames<\/strong>. La probabilit\u00e9 d&rsquo;observer une donn\u00e9e au-del\u00e0 de ces bornes sous l&rsquo;hypoth\u00e8se de normalit\u00e9 devient infinit\u00e9simale (environ <span class=\"arabpsy-math-inline\" data-latex=\"0{,}0002%\">$0{,}0002%$<\/span>).<\/li>\n<\/ul>\n<p>Sur le plan d\u00e9ontologique et \u00e9pist\u00e9mologique, l&rsquo;\u00e9radication aveugle des valeurs aberrantes repr\u00e9sente un travers m\u00e9thodologique majeur. En psychologie cognitive et en psychiatrie exp\u00e9rimentale, une valeur extr\u00eame ne constitue pas syst\u00e9matiquement une erreur de mesure ou un bruit exp\u00e9rimental ind\u00e9sirable ; elle peut mat\u00e9rialiser l&rsquo;expression clinique authentique d&rsquo;un individu neuroatypique ou un ph\u00e9nom\u00e8ne de rupture attentionnelle d&rsquo;un int\u00e9r\u00eat scientifique sup\u00e9rieur. L&rsquo;analyste a l&rsquo;obligation \u00e9thique d&rsquo;examiner chaque valeur \u00e9tiquet\u00e9e comme atypique avant d&rsquo;envisager son exclusion ou sa neutralisation.<\/p>\n<h3>9.2 Automatisation de la d\u00e9tection et du nettoyage sous R<\/h3>\n<p>L&rsquo;impl\u00e9mentation de la r\u00e8gle de Tukey peut \u00eatre automatis\u00e9e de fa\u00e7on robuste par la cr\u00e9ation d&rsquo;une fonction R personnalis\u00e9e, capable d&rsquo;\u00e9tiqueter chaque observation d&rsquo;un vecteur selon son statut m\u00e9trologique :<\/p>\n<p><code>detecter_aberrations &lt;- function(vecteur, facteur_k = 1.5) {<br \/>\n  q1 &lt;- quantile(vecteur, 0.25, na.rm = TRUE)<br \/>\n  q3 &lt;- quantile(vecteur, 0.75, na.rm = TRUE)<br \/>\n  iqr_valeur &lt;- q3 - q1<\/p>\n<p>  borne_inf &lt;- q1 - (facteur_k * iqr_valeur)<br \/>\n  borne_sup &lt;- q3 + (facteur_k * iqr_valeur)<\/p>\n<p>  est_aberrant &lt;- vecteur &lt; borne_inf | vecteur &gt; borne_sup<br \/>\n  return(est_aberrant)<br \/>\n}<\/code><\/p>\n<p>Appliquons cette fonction sur un jeu de donn\u00e9es simul\u00e9 repr\u00e9sentant des temps de r\u00e9ponse (en millisecondes) au sein d&rsquo;une t\u00e2che de m\u00e9moire de travail, pollu\u00e9 par deux scores de distraction extr\u00eame :<\/p>\n<p><code>temps_reponse &lt;- c(450, 480, 510, 525, 530, 545, 560, 580, 610, 1450, 1820)<\/code><br \/>\n<code>aberrations &lt;- detecter_aberrations(temps_reponse, facteur_k = 1.5)<\/code><br \/>\n<code>donnees_chronometriques &lt;- data.frame(TR = temps_reponse, Aberrant = aberrations)<\/code><br \/>\n<code>print(donnees_chronometriques)<\/code><\/p>\n<p>La console signale imm\u00e9diatement les deux derni\u00e8res observations ($1450$ et $1820$) avec l&rsquo;\u00e9tiquette logique <code>TRUE<\/code>. Pour proc\u00e9der \u00e0 un filtrage conditionnel propre dans le cadre d&rsquo;une analyse de sensibilit\u00e9 au sein du Tidyverse, nous pouvons combiner <code>filter()<\/code> et la n\u00e9gation logique de la condition :<\/p>\n<p><code>donnees_filtrees &lt;- donnees_chronometriques |&gt;<br \/>\n  filter(!Aberrant)<\/code><\/p>\n<p>Il est imp\u00e9ratif, dans les rapports de recherche scientifique, de mener une analyse de sensibilit\u00e9 syst\u00e9matique : l&rsquo;analyste doit ex\u00e9cuter ses mod\u00e9lisations inf\u00e9rentielles avec <em>et<\/em> sans les valeurs aberrantes d\u00e9tect\u00e9es via l&rsquo;IQR. Si les conclusions et la significativit\u00e9 statistique (<span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>-valeurs, tailles d&rsquo;effet) divergent substantiellement entre les deux approches, cette instabilit\u00e9 doit \u00eatre int\u00e9gralement divulgu\u00e9e et discut\u00e9e dans le manuscrit soumis pour publication.<\/p>\n<h2>10. Visualisation graphique de l&rsquo;IQR avec Base R et ggplot2<\/h2>\n<h3>10.1 La bo\u00eete \u00e0 moustaches (Boxplot) en R de base<\/h3>\n<p>Le mode de repr\u00e9sentation graphique privil\u00e9gi\u00e9 de l&rsquo;\u00e9cart interquartile est la bo\u00eete \u00e0 moustaches, con\u00e7ue par John Tukey sous l&rsquo;appellation anglophone de <em>box-and-whisker plot<\/em>. L&rsquo;environnement de base de R dispose d&rsquo;un moteur graphique ultra-performant permettant de g\u00e9n\u00e9rer ces diagrammes sans n\u00e9cessiter le chargement de packages tiers :<\/p>\n<p><code>scores_test &lt;- c(12, 15, 18, 19, 21, 22, 23, 24, 25, 27, 28, 30, 45)<\/code><br \/>\n<code>boxplot(scores_test, main = \"Distribution des scores au test\",<br \/>\n        ylab = \"Points obtenus\", col = \"lightgray\", border = \"black\")<\/code><\/p>\n<p>L&rsquo;anatomie graphique de cette bo\u00eete constitue une incarnation visuelle directe de l&rsquo;\u00e9cart interquartile :<\/p>\n<ul>\n<li>La ligne horizontale \u00e9paisse mat\u00e9rialis\u00e9e au c\u0153ur de la bo\u00eete repr\u00e9sente la <strong>m\u00e9diane<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"Q_2\">$Q_2$<\/span>).<\/li>\n<li>Le bord inf\u00e9rieur de la bo\u00eete correspond pr\u00e9cis\u00e9ment au <strong>premier quartile<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"Q_1\">$Q_1$<\/span>).<\/li>\n<li>Le bord sup\u00e9rieur de la bo\u00eete correspond au <strong>troisi\u00e8me quartile<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"Q_3\">$Q_3$<\/span>).<\/li>\n<li>Par voie de cons\u00e9quence directe, la <strong>hauteur totale du corps de la bo\u00eete<\/strong> mat\u00e9rialise visuellement la magnitude exacte de l&rsquo;<strong>\u00e9cart interquartile (IQR)<\/strong>. Plus la bo\u00eete est \u00e9tir\u00e9e verticalement, plus la dispersion des 50 % centraux de l&rsquo;\u00e9chantillon est prononc\u00e9e.<\/li>\n<\/ul>\n<p>Les moustaches s&rsquo;\u00e9tendent g\u00e9n\u00e9ralement jusqu&rsquo;\u00e0 la valeur observ\u00e9e la plus \u00e9loign\u00e9e qui demeure \u00e0 l&rsquo;int\u00e9rieur des cl\u00f4tures de Tukey (<span class=\"arabpsy-math-inline\" data-latex=\"1{,}5 \\times \\text{IQR}\">$1{,}5 \\times \\text{IQR}$<\/span>). Les points isol\u00e9s figurant au-del\u00e0 des extr\u00e9mit\u00e9s des moustaches visualisent explicitement les observations atypiques individuelles (comme la valeur 45 dans l&rsquo;exemple ci-dessus), permettant un diagnostic perceptuel instantan\u00e9 de l&rsquo;\u00e9tat de la distribution.<\/p>\n<h3>10.2 Repr\u00e9sentations avanc\u00e9es et esth\u00e9tiques de l&rsquo;IQR avec ggplot2<\/h3>\n<p>Si la fonction graphique de base excelle pour l&rsquo;exploration rapide, la conception de visualisations de haute qualit\u00e9 typographique destin\u00e9es aux manuscrits finaux requiert l&rsquo;\u00e9cosyst\u00e8me <a rel=\"noopener noreferrer\" href=\"https:\/\/ggplot2.tidyverse.org\/\">ggplot2<\/a>. La puissance de la grammaire des graphiques permet de superposer l&rsquo;IQR \u00e0 d&rsquo;autres couches morphologiques pour surmonter le principal \u00e9cueil du boxplot standard : son incapacit\u00e9 \u00e0 r\u00e9v\u00e9ler la bimodalit\u00e9 sous-jacente.<\/p>\n<p>Le script suivant assemble une visualisation combinant un diagramme en violon (r\u00e9v\u00e9lant la fonction d&rsquo;estimation de densit\u00e9 par noyau), une bo\u00eete \u00e0 moustaches int\u00e9gr\u00e9e (mettant en exergue l&rsquo;IQR et la m\u00e9diane), et un nuage de points individualis\u00e9s transparents soumis \u00e0 une gigue al\u00e9atoire (<em>jitter<\/em>) :<\/p>\n<p><code>library(ggplot2)<\/code><\/p>\n<p><code>set.seed(42)<\/code><br \/>\n<code>donnees_visu &lt;- data.frame(<br \/>\n  groupe = factor(rep(c(\"Controle\", \"Intervention\"), each = 60)),<br \/>\n  score = c(rnorm(60, mean = 50, sd = 10), rnorm(60, mean = 62, sd = 18))<br \/>\n)<\/code><\/p>\n<p><code>graphique_apa &lt;- ggplot(donnees_visu, aes(x = groupe, y = score)) +<br \/>\n  geom_violin(trim = FALSE, fill = \"gray95\", color = \"gray60\") +<br \/>\n  geom_boxplot(width = 0.2, fill = \"white\", color = \"black\", outlier.shape = NA, alpha = 0.8) +<br \/>\n  geom_jitter(width = 0.08, alpha = 0.35, color = \"black\", size = 1.5) +<br \/>\n  theme_classic() +<br \/>\n  labs(<br \/>\n    title = \"Comparaison de la dispersion centrale inter-groupes\",<br \/>\n    subtitle = \"Superposition de l'estimation de densit\u00e9, de l'IQR et des points individuels\",<br \/>\n    x = \"Condition exp\u00e9rimentale\",<br \/>\n    y = \"Score psychom\u00e9trique standardis\u00e9\"<br \/>\n  ) +<br \/>\n  theme(<br \/>\n    plot.title = element_text(face = \"bold\", size = 14),<br \/>\n    axis.title = element_text(face = \"bold\", size = 12),<br \/>\n    axis.text = element_text(color = \"black\", size = 10)<br \/>\n  )<\/code><\/p>\n<p><code>print(graphique_apa)<\/code><\/p>\n<p>Ce type de figure s&rsquo;aligne fid\u00e8lement sur les chartes graphiques de l&rsquo;APA : \u00e9limination des grilles d&rsquo;arri\u00e8re-plan redondantes, utilisation d&rsquo;un contraste monochrome \u00e9l\u00e9gant et lisible en noir et blanc, et conservation de l&rsquo;int\u00e9grit\u00e9 brute des donn\u00e9es gr\u00e2ce aux points sous-jacents. Cette d\u00e9marche offre une transparence m\u00e9thodologique totale, salu\u00e9e par les \u00e9valuateurs universitaires.<\/p>\n<h2>11. Arbitrage statistique : \u00c9cart interquartile versus \u00c9cart-type<\/h2>\n<h3>11.1 Crit\u00e8res de s\u00e9lection fond\u00e9s sur la forme de la distribution<\/h3>\n<p>Le dilemme m\u00e9trologique opposant l&rsquo;\u00e9cart interquartile \u00e0 l&rsquo;\u00e9cart-type n&rsquo;est pas une simple affaire de pr\u00e9f\u00e9rence stylistique ; il proc\u00e8de d&rsquo;une \u00e9valuation math\u00e9matique rigoureuse de la structure de l&rsquo;\u00e9chantillon. L&rsquo;analyste doit engager un diagnostic de normalit\u00e9 formalis\u00e9 avant de statuer sur le couple d&rsquo;indicateurs de position et de dispersion \u00e0 privil\u00e9gier :<\/p>\n<ul>\n<li><strong>Tests d&rsquo;hypoth\u00e8se formels :<\/strong> Le test de normalit\u00e9 de Shapiro-Wilk (impl\u00e9ment\u00e9 sous R via <code>shapiro.test(x)<\/code>) \u00e9value l&rsquo;ad\u00e9quation des donn\u00e9es \u00e0 une loi normale. Une valeur $p < 0{,}05$ indique une d\u00e9viation significative rejetant l'hypoth\u00e8se de normalit\u00e9 univari\u00e9e.<\/li>\n<li><strong>Coefficients de forme :<\/strong> L&rsquo;analyse conjointe de l&rsquo;asym\u00e9trie (<em>skewness<\/em>) et de l&rsquo;aplatissement (<em>kurtosis<\/em>). Des valeurs absolues du coefficient d&rsquo;asym\u00e9trie d\u00e9passant le seuil de <span class=\"arabpsy-math-inline\" data-latex=\"1{,}0\">$1{,}0$<\/span> (ou <span class=\"arabpsy-math-inline\" data-latex=\"2{,}0\">$2{,}0$<\/span> selon les standards conservateurs) signalent une dissym\u00e9trie incompatible avec l&rsquo;usage fid\u00e8le de l&rsquo;\u00e9cart-type.<\/li>\n<li><strong>Inspections visuelles :<\/strong> L&rsquo;examen des diagrammes quantile-quantile (<em>Q-Q plots<\/em> via <code>qqnorm(x) ; qqline(x)<\/code>) permet de d\u00e9tecter les d\u00e9crochages structurels des observations par rapport \u00e0 la droite de normalit\u00e9 th\u00e9orique.<\/li>\n<\/ul>\n<p>Dans l&rsquo;\u00e9ventualit\u00e9 o\u00f9 la distribution satisfait aux conditions de stricte normalit\u00e9, l&rsquo;\u00e9cart-type constitue sans conteste l&rsquo;estimateur param\u00e9trique le plus efficace au sens de Fisher, exploitant l&rsquo;int\u00e9gralit\u00e9 des \u00e9carts quadratiques \u00e0 la moyenne et maximisant la puissance statistique (atteinte de la borne de Cram\u00e9r-Rao). Toutefois, d\u00e8s lors qu&rsquo;apparaissent des contaminations de donn\u00e9es, une asym\u00e9trie marqu\u00e9e ou des distributions multimodales, la sup\u00e9riorit\u00e9 op\u00e9rationnelle du tandem M\u00e9diane \/ \u00c9cart interquartile devient indiscutable pour pr\u00e9server la v\u00e9rit\u00e9 descriptive des donn\u00e9es psychologiques observ\u00e9es.<\/p>\n<h3>11.2 Implications pour l&rsquo;analyse inf\u00e9rentielle ult\u00e9rieure<\/h3>\n<p>Le diagnostic de dispersion via l&rsquo;IQR poss\u00e8de un impact d\u00e9terminant sur le choix de la mod\u00e9lisation inf\u00e9rentielle ult\u00e9rieure. Lorsqu&rsquo;une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 massive des \u00e9carts interquartiles est mise en \u00e9vidence entre plusieurs groupes exp\u00e9rimentaux, cela indique fr\u00e9quemment une violation s\u00e9v\u00e8re du postulat d&rsquo;homosc\u00e9dasticit\u00e9 (homog\u00e9n\u00e9it\u00e9 des variances). Face \u00e0 un tel constat, les mod\u00e8les param\u00e9triques lin\u00e9aires classiques (comme le test <span class=\"arabpsy-math-inline\" data-latex=\"t\">$t$<\/span> de Student ou l&rsquo;Analyse de Variance univari\u00e9e \/ ANOVA) perdent leur calibration du risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce (<span class=\"arabpsy-math-inline\" data-latex=\"\\alpha\">$\\alpha$<\/span>). Le chercheur est alors orient\u00e9 vers des m\u00e9thodes non param\u00e9triques (test de Wilcoxon-Mann-Whitney, test de Kruskal-Wallis) ou vers des mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s (GLM) autorisant des structures de variance flexibles.<\/p>\n<p>De plus, l&rsquo;\u00e9cart interquartile intervient de mani\u00e8re pr\u00e9pond\u00e9rante dans le pr\u00e9traitement des donn\u00e9es massives via les techniques de standardisation robuste (<em>Robust Scaling<\/em>). Contrairement \u00e0 la standardisation <span class=\"arabpsy-math-inline\" data-latex=\"Z\">$Z$<\/span> classique \u2014 qui soustrait la moyenne empirique et divise par l&rsquo;\u00e9cart-type, propageant ainsi les distorsions induites par les valeurs aberrantes \u2014, la mise \u00e0 l&rsquo;\u00e9chelle robuste transforme les scores selon la fonction suivante :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"x_{\\text{robuste}} = \\frac{x - \\text{M\u00e9diane}(X)}{\\text{IQR}(X)}\">$$x_{\\text{robuste}} = \\frac{x &#8211; \\text{M\u00e9diane}(X)}{\\text{IQR}(X)}$$<\/span><\/p>\n<p>Cette transformation, facilement scriptable dans R via <code>(x - median(x)) \/ IQR(x)<\/code>, stabilise les distributions avant leur injection dans des algorithmes d&rsquo;apprentissage automatique (comme les machines \u00e0 vecteurs de support ou les r\u00e9seaux de neurones appliqu\u00e9s au neuro-diagnostic) en \u00e9vitant que les participants atypiques n&rsquo;\u00e9crasent le poids computationnel des covariables centrales.<\/p>\n<h2>12. Directives de r\u00e9daction acad\u00e9mique et pi\u00e8ges fr\u00e9quents dans R<\/h2>\n<h3>12.1 Normes de publication APA pour la communication de l&rsquo;IQR<\/h3>\n<p>La transmission fid\u00e8le des donn\u00e9es quantitatives dans les revues internationales exige une adh\u00e9sion scrupuleuse aux prescriptions stylistiques \u00e9dict\u00e9es par la 7e \u00e9dition du manuel de publication de l&rsquo;American Psychological Association (APA). Lorsqu&rsquo;une distribution est r\u00e9sum\u00e9e \u00e0 l&rsquo;aide d&rsquo;indicateurs non param\u00e9triques, l&rsquo;APA prescrit d&rsquo;associer invariablement la m\u00e9diane \u00e0 l&rsquo;\u00e9cart interquartile, de la m\u00eame mani\u00e8re que la moyenne arithm\u00e9tique s&rsquo;accompagne syst\u00e9matiquement de l&rsquo;\u00e9cart-type.<\/p>\n<p>Dans le corps du texte, les acronymes statistiques doivent \u00eatre typographi\u00e9s avec exactitude. L&rsquo;abr\u00e9viation de la m\u00e9diane s&rsquo;\u00e9crit avec une majuscule suivie de lettres minuscules en italique (<em>Mdn<\/em>), tandis que l&rsquo;\u00e9cart interquartile s&rsquo;abr\u00e8ge en lettres capitales romaines (sans italique) sous la forme <em>IQR<\/em>. Voici des exemples d&rsquo;insertion textuelle conformes aux normes :<\/p>\n<p><em>\u00ab L&rsquo;\u00e9valuation de l&rsquo;anxi\u00e9t\u00e9 pr\u00e9-op\u00e9ratoire met en \u00e9vidence une dispersion mod\u00e9r\u00e9e au sein du groupe trait\u00e9 (Mdn = 24,50 ; IQR = 8,00), contrastant avec l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 prononc\u00e9e document\u00e9e dans le groupe t\u00e9moin (Mdn = 38,00 ; IQR = 17,25). \u00bb<\/em><\/p>\n<p><em>\u00ab En raison d&rsquo;une asym\u00e9trie positive prononc\u00e9e des latences d&rsquo;inhibition attentionnelle (skewness = 2,41), les donn\u00e9es sont r\u00e9sum\u00e9es par leur m\u00e9diane et leur intervalle interquartile : les participants sous condition de privation de sommeil ont affich\u00e9 une dispersion accrue (Mdn = 612 ms, IQR = 148 ms) comparativement aux sujets repos\u00e9s (Mdn = 445 ms, IQR = 52 ms). \u00bb<\/em><\/p>\n<p>Lors de la pr\u00e9sentation de donn\u00e9es synth\u00e9tiques au sein d&rsquo;un tableau acad\u00e9mique, une colonne sp\u00e9cifique intitul\u00e9e <em>IQR<\/em> doit \u00eatre ins\u00e9r\u00e9e, et les notes de bas de tableau doivent expliciter les modalit\u00e9s de calcul : <em>\u00ab Note. Mdn = m\u00e9diane ; IQR = \u00e9cart interquartile calcul\u00e9 selon l&rsquo;algorithme Type 7 de R (Hyndman &#038; Fan, 1996). \u00bb<\/em> Cette pr\u00e9cision d\u00e9montre une rigueur m\u00e9trologique sans faille et garantit la r\u00e9plicabilit\u00e9 exacte des r\u00e9sultats pr\u00e9sent\u00e9s.<\/p>\n<h3>12.2 Erreurs courantes commises par les chercheurs et solutions de d\u00e9bogage<\/h3>\n<p>L&rsquo;exploitation de l&rsquo;IQR dans l&rsquo;environnement R expose le chercheur \u00e0 divers pi\u00e8ges m\u00e9thodologiques et informatiques document\u00e9s dans les pratiques de laboratoire :<\/p>\n<ul>\n<li><strong>Confusion conceptuelle entre IQR et SIQR :<\/strong> Une erreur r\u00e9currente consiste \u00e0 confondre l&rsquo;\u00e9cart interquartile total (<span class=\"arabpsy-math-inline\" data-latex=\"\\text{IQR} = Q_3 - Q_1\">$\\text{IQR} = Q_3 &#8211; Q_1$<\/span>) avec l&rsquo;\u00e9cart <em>semi-interquartile<\/em> (<span class=\"arabpsy-math-inline\" data-latex=\"\\text{SIQR} = (Q_3 - Q_1) \/ 2\">$\\text{SIQR} = (Q_3 &#8211; Q_1) \/ 2$<\/span>). Le semi-interquartile est parfois utilis\u00e9 comme un analogue visuel d&rsquo;un demi-intervalle sym\u00e9trique autour de la m\u00e9diane, mais son assimilation trompeuse \u00e0 l&rsquo;IQR induit une sous-estimation syst\u00e9matique de 50 % de la dispersion r\u00e9elle rapport\u00e9e.<\/li>\n<li><strong>Erreurs de typage silencieuses :<\/strong> L&rsquo;ex\u00e9cution d&rsquo;op\u00e9rations d&rsquo;agr\u00e9gation group\u00e9e sur des variables import\u00e9es au format date, facteur ordinal ou cha\u00eene de caract\u00e8res peut provoquer des coercitions invisibles o\u00f9 R utilise les rangs entiers internes du facteur au lieu des valeurs num\u00e9riques sous-jacentes. L&rsquo;application pr\u00e9alable de <code>str(donnees)<\/code> permet de neutraliser ce risque.<\/li>\n<li><strong>Omission de na.rm = TRUE dans les fonctions automatis\u00e9es :<\/strong> L&rsquo;\u00e9criture de fonctions d&rsquo;analyse par lots sans inclusion explicite de l&rsquo;exclusion des valeurs manquantes entra\u00eene l&rsquo;\u00e9chec silencieux de scripts entiers, renvoyant des tables d&rsquo;analyse pollu\u00e9es de valeurs <code>NA<\/code>.<\/li>\n<li><strong>N\u00e9gligence de l&rsquo;environnement computationnel :<\/strong> Ne pas consigner la version pr\u00e9cise de R et des packages mobilis\u00e9s met en p\u00e9ril la reproductibilit\u00e9. L&rsquo;encapsulation de l&rsquo;analyse dans un flux contr\u00f4l\u00e9 via le package <code>renv<\/code> et l&rsquo;insertion syst\u00e9matique de la commande <code>sessionInfo()<\/code> au terme du script garantissent la tra\u00e7abilit\u00e9 des environnements num\u00e9riques.<\/li>\n<\/ul>\n<p>En observant ces recommandations formelles et en ma\u00eetrisant les m\u00e9canismes algorithmiques de la fonction <code>IQR()<\/code>, les chercheurs en psychologie et sciences du comportement s&rsquo;assurent d&rsquo;une analyse quantitative rigoureuse, transparente et conforme aux plus hauts standards m\u00e9thodologiques contemporains.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<p>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7th ed.). American Psychological Association. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/p>\n<p>Beck, A. T., Steer, R. A., &#038; Brown, G. K. (1996). <em>Manual for the Beck Depression Inventory-II<\/em>. Psychological Corporation.<\/p>\n<p>Hyndman, R. J., &#038; Fan, Y. (1996). Sample quantiles in statistical packages. <em>The American Statistician<\/em>, 50(4), 361-365. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00031305.1996.10473566\">https:\/\/doi.org\/10.1080\/00031305.1996.10473566<\/a><\/p>\n<p>Little, R. J. A., &#038; Rubin, D. B. (2019). <em>Statistical analysis with missing data<\/em> (3rd ed.). John Wiley &#038; Sons. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/9781119482260\">https:\/\/doi.org\/10.1002\/9781119482260<\/a><\/p>\n<p>R Core Team. (2023). <em>R: A language and environment for statistical computing<\/em>. R Foundation for Statistical Computing. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.R-project.org\/\">https:\/\/www.R-project.org\/<\/a><\/p>\n<p>Tukey, J. W. (1977). <em>Exploratory data analysis<\/em>. Addison-Wesley.<\/p>\n<p>Wickham, H., Averick, M., Bryan, J., Chang, W., McGowan, L. D., Fran\u00e7ois, R., Grolemund, G., Hayes, A., Henry, L., Hester, J., Kuhn, M., Pedersen, T. L., Miller, E., Bache, S. M., M\u00fcller, K., Ooms, J., Robinson, D., Seidel, D. P., Spinu, V., \u2026 Yutani, H. (2019). Welcome to the Tidyverse. <em>Journal of Open Source Software<\/em>, 4(43), Article 1686. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.21105\/joss.01686\">https:\/\/doi.org\/10.21105\/joss.01686<\/a><\/p>\n<p>Wilcox, R. R. (2017). <em>Introduction to robust estimation and hypothesis testing<\/em> (4th ed.). Academic Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/C2015-0-02384-5\">https:\/\/doi.org\/10.1016\/C2015-0-02384-5<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour calculer et interpr\u00e9ter l&rsquo;\u00e9cart interquartile (IQR) dans R, avec applications concr\u00e8tes en recherche psychologique.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2193,2169],"tags":[1982,2211,16,2658,2212],"class_list":["post-2226","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-programmation-r","category-psychologie-quantitative","tag-ecart-interquartile","tag-iqr","tag-psychometrie","tag-r","tag-statistiques-descriptives"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/2226","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=2226"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/2226\/revisions"}],"predecessor-version":[{"id":16343,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/2226\/revisions\/16343"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=2226"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=2226"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=2226"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}