{"id":4887,"date":"2026-10-14T18:01:19","date_gmt":"2026-10-14T18:01:19","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-appliquer-la-regle-empirique-dans-r\/"},"modified":"2026-10-14T18:01:19","modified_gmt":"2026-10-14T18:01:19","slug":"comment-appliquer-la-regle-empirique-dans-r","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-appliquer-la-regle-empirique-dans-r","title":{"rendered":"Comment appliquer la r\u00e8gle empirique dans R"},"content":{"rendered":"<p>Dans le domaine de l&rsquo;\u00e9valuation psychologique et des sciences quantitatives du comportement, la mod\u00e9lisation statistique constitue le socle indispensable sur lequel reposent l&rsquo;interpr\u00e9tation diagnostique et la validation des instruments de mesure. La compr\u00e9hension de la distribution des traits latents, qu&rsquo;il s&rsquo;agisse de dimensions cognitives complexes ou d&rsquo;indicateurs de psychopathologie, d\u00e9pend fondamentalement de mod\u00e8les probabilistes robustes. Parmi ces abstractions math\u00e9matiques, la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Loi_normale\">distribution normale de Gauss-Laplace<\/a> occupe une position pr\u00e9\u00e9minente, servant \u00e0 la fois de r\u00e9f\u00e9rence normative pour le calibrage des scores standardis\u00e9s et de pilier th\u00e9orique pour les m\u00e9thodologies inf\u00e9rentielles modernes.<\/p>\n<p>Au c\u0153ur de cette architecture se trouve un th\u00e9or\u00e8me fondamental de dispersion souvent vulgaris\u00e9 sous la d\u00e9signation de \u00ab r\u00e8gle empirique \u00bb, ou r\u00e8gle des 68-95-99,7. Cette formulation condens\u00e9e exprime les propri\u00e9t\u00e9s g\u00e9om\u00e9triques et int\u00e9grales immuables d&rsquo;une courbe en cloche parfaite, fixant avec une r\u00e9gularit\u00e9 remarquable la proportion d&rsquo;observations attendues \u00e0 l&rsquo;int\u00e9rieur d&rsquo;un, deux et trois \u00e9carts-types autour de l&rsquo;esp\u00e9rance math\u00e9matique. Pour le psychom\u00e8tre, le neuropsychologue ou le chercheur en neurosciences computationnelles, cette r\u00e8gle ne repr\u00e9sente pas une simple curiosit\u00e9 arithm\u00e9tique, mais une grille de lecture op\u00e9rationnelle servant \u00e0 discerner la typicit\u00e9 de la d\u00e9viance clinique et \u00e0 quantifier la raret\u00e9 statistique d&rsquo;un profil cognitif donn\u00e9.<\/p>\n<p>Toutefois, la transposition de ce mod\u00e8le th\u00e9orique id\u00e9alis\u00e9 vers les donn\u00e9es empiriques recueillies sur le terrain n\u00e9cessite une ma\u00eetrise technique rigoureuse. L&rsquo;environnement de calcul statistique <a rel=\"noopener noreferrer\" href=\"https:\/\/www.r-project.org\/\">R<\/a> offre une infrastructure logicielle d&rsquo;une flexibilit\u00e9 in\u00e9gal\u00e9e pour impl\u00e9menter, tester et visualiser les postulats de la r\u00e8gle empirique. Ce guide propose une exploration approfondie, th\u00e9orique, algorithmique et pratique de l&rsquo;application de cette r\u00e8gle fondamentale sous R, en articulant les d\u00e9rivations math\u00e9matiques sous-jacentes avec les imp\u00e9ratifs cliniques de la mesure psychologique moderne.<\/p>\n<h2>1. Fondements th\u00e9oriques de la r\u00e8gle empirique en psychom\u00e9trie<\/h2>\n<h3>1.1 D\u00e9finition formelle de la r\u00e8gle 68-95-99,7<\/h3>\n<p>La r\u00e8gle empirique, fr\u00e9quemment identifi\u00e9e dans la litt\u00e9rature anglo-saxonne sous l&rsquo;appellation de <em>three-sigma rule<\/em>, stipule que pour toute variable al\u00e9atoire continue suivant rigoureusement une loi normale de moyenne th\u00e9orique &mu; et d&rsquo;\u00e9cart-type &sigma;, la quasi-totalit\u00e9 des r\u00e9alisations se r\u00e9partit de mani\u00e8re d\u00e9terministe au sein de trois intervalles centr\u00e9s. De fa\u00e7on formelle, la proportion d&rsquo;observations comprises entre &mu; &minus; &sigma; et &mu; + &sigma; est approximativement \u00e9gale \u00e0 68,27 % ; la proportion situ\u00e9e entre &mu; &minus; 2&sigma; et &mu; + 2&sigma; atteint 95,45 % ; et l&rsquo;intervalle d\u00e9limit\u00e9 par &mu; &minus; 3&sigma; et &mu; + 3&sigma; englobe 99,73 % de la masse de probabilit\u00e9 totale de la distribution.<\/p>\n<p>D&rsquo;un point de vue conceptuel, ce d\u00e9coupage r\u00e9gulier fournit une m\u00e9trique standardis\u00e9e universelle pour segmenter la variabilit\u00e9 d&rsquo;un ph\u00e9nom\u00e8ne biologique ou psychologique. L&rsquo;\u00e9cart-type ne se contente pas de mesurer l&rsquo;\u00e9talement brut des valeurs ; il devient une unit\u00e9 de mesure fonctionnelle \u00e0 part enti\u00e8re. Le passage de l&rsquo;intervalle \u00e0 un \u00e9cart-type \u00e0 celui \u00e0 deux \u00e9carts-types ne traduit pas seulement une augmentation lin\u00e9aire de la largeur de la fen\u00eatre observ\u00e9e, mais une rar\u00e9faction abrupte de la densit\u00e9 marginale, illustrant la convergence rapide des queues de la distribution vers l&rsquo;asymptote horizontale z\u00e9ro.<\/p>\n<p>Sur le plan de l&rsquo;histoire des sciences, cette formulation r\u00e9sulte d&rsquo;une s\u00e9dimentation conceptuelle majeure reliant les travaux d&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Abraham_de_Moivre\">Abraham de Moivre<\/a> au d\u00e9but du XVIIIe si\u00e8cle, qui d\u00e9couvrit l&rsquo;approximation de la loi binomiale par la loi normale, aux d\u00e9veloppements analytiques ult\u00e9rieurs de Carl Friedrich Gauss et de Pierre-Simon de Laplace concernant la th\u00e9orie des erreurs d&rsquo;observation. \u00c0 la fin du XIXe si\u00e8cle, le biom\u00e9tricien Karl Pearson a formalis\u00e9 l&rsquo;int\u00e9gration syst\u00e9matique de l&rsquo;\u00e9cart-type comme mesure de r\u00e9f\u00e9rence, conf\u00e9rant \u00e0 la r\u00e8gle son statut de standard pour les sciences empiriques naissantes.<\/p>\n<p>Il est fondamental de maintenir une distinction \u00e9pist\u00e9mologique stricte entre la r\u00e8gle empirique en tant qu&rsquo;expression math\u00e9matique exacte appliqu\u00e9e \u00e0 un continuum th\u00e9orique infini et son comportement lorsqu&rsquo;elle est confront\u00e9e \u00e0 un \u00e9chantillonnage fini. Dans le monde th\u00e9orique, les valeurs 68,27 %, 95,45 % et 99,73 % sont des constantes immuables d\u00e9riv\u00e9es du calcul diff\u00e9rentiel et int\u00e9gral. Dans les jeux de donn\u00e9es r\u00e9els, les fluctuations d&rsquo;\u00e9chantillonnage, les erreurs de mesure inh\u00e9rentes aux questionnaires et l&rsquo;impr\u00e9cision des instruments font que les pourcentages observ\u00e9s ne constituent que des approximations fluctuant autour de ces bornes canoniques.<\/p>\n<h3>1.2 Pertinence de la loi normale dans les mesures psychologiques<\/h3>\n<p>L&rsquo;omnipr\u00e9sence du mod\u00e8le gaussien dans les sciences de l&rsquo;esprit ne d\u00e9coule pas d&rsquo;un choix arbitraire, mais s&rsquo;enracine dans les cons\u00e9quences directes du <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Th%C3%A9or%C3%A8me_central_limite\">th\u00e9or\u00e8me central limite<\/a>. Ce th\u00e9or\u00e8me math\u00e9matique d\u00e9montre que la somme ou la r\u00e9sultante lin\u00e9aire d&rsquo;un grand nombre de variables al\u00e9atoires ind\u00e9pendantes, ou faiblement corr\u00e9l\u00e9es, tend asymptotiquement vers une loi normale, quelles que soient les distributions sous-jacentes de ces composantes \u00e9l\u00e9mentaires. En psychologie cognitive et diff\u00e9rentielle, un score global \u00e0 un test n&rsquo;est jamais l&rsquo;expression d&rsquo;un facteur isol\u00e9, mais l&rsquo;aboutissement macroscopique d&rsquo;une myriade de processus physiologiques, g\u00e9n\u00e9tiques, environnementaux et attentionnels interagissant de concert.<\/p>\n<p>Cette r\u00e9alit\u00e9 biologique et psychique sous-tend la construction des batteries psychom\u00e9triques contemporaines les plus \u00e9prouv\u00e9es. Les \u00e9chelles d&rsquo;intelligence de David Wechsler, telles que la WAIS-IV ou le WISC-V, sont express\u00e9ment \u00e9talonn\u00e9es pour que le Quotient Intellectuel (QI) se distribue normalement avec une moyenne de population fix\u00e9e \u00e0 100 et un \u00e9cart-type de 15. De m\u00eame, les inventaires de personnalit\u00e9 multidimensionnels, \u00e0 l&rsquo;image du <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/NEO-PI-R\">NEO-PI-3<\/a> fond\u00e9 sur le mod\u00e8le en cinq facteurs, utilisent des m\u00e9triques d\u00e9riv\u00e9es o\u00f9 les traits latents adoptent cette g\u00e9om\u00e9trie sym\u00e9trique. L&rsquo;utilisation de ces rep\u00e8res permet d&rsquo;\u00e9tablir une correspondance directe et univoque entre la position d&rsquo;un individu dans la cohorte et sa raret\u00e9 probabiliste.<\/p>\n<p>Les propri\u00e9t\u00e9s intrins\u00e8ques de sym\u00e9trie parfaite et d&rsquo;unimodalit\u00e9 de la fonction gaussienne conditionnent l&rsquo;ensemble de l&rsquo;inf\u00e9rence clinique. L&rsquo;absence de distorsion directionnelle implique que la moyenne, la m\u00e9diane et le mode co\u00efncident rigoureusement au centre de la distribution. En recherche clinique, cette convergence asymptotique garantit que les d\u00e9viations positives et n\u00e9gatives d&rsquo;amplitude identique poss\u00e8dent rigoureusement la m\u00eame probabilit\u00e9 d&rsquo;occurrence, simplifiant ainsi le calcul des intervalles de confiance et la mod\u00e9lisation des erreurs r\u00e9siduelles dans les \u00e9quations structurelles.<\/p>\n<p>Le recours aux scores standardis\u00e9s Z d\u00e9coule directement de cette structure. En retranchant la moyenne empirique d&rsquo;un score brut et en divisant le r\u00e9sidu obtenu par l&rsquo;\u00e9cart-type de l&rsquo;\u00e9chantillon de r\u00e9f\u00e9rence, le chercheur convertit une m\u00e9trique brute h\u00e9t\u00e9rog\u00e8ne en un \u00e9cart norm\u00e9 sans unit\u00e9. Un score Z de +2,00 indique imm\u00e9diatement, en vertu de la r\u00e8gle empirique, que le sujet se situe au-del\u00e0 d&rsquo;environ 97,7 % des individus de sa population d&rsquo;\u00e9talonnage. Cette standardisation unitaire est la condition <em>sine qua non<\/em> de toute comparaison inter-instruments valide dans les bilans neuropsychologiques.<\/p>\n<h3>1.3 Crit\u00e8res d&rsquo;applicabilit\u00e9 et hypoth\u00e8ses pr\u00e9alables<\/h3>\n<p>L&rsquo;utilisation non critique de la r\u00e8gle 68-95-99,7 expose le statisticien \u00e0 des b\u00e9vues interpr\u00e9tatives majeures d\u00e8s lors que les donn\u00e9es sous-jacentes transgressent les postulats fondateurs du mod\u00e8le continu. Le premier crit\u00e8re fondamental concerne l&rsquo;\u00e9chelle de mesure : la variable \u00e9tudi\u00e9e doit imp\u00e9rativement relever d&rsquo;une \u00e9chelle d&rsquo;intervalles constants ou d&rsquo;une \u00e9chelle de rapports, au sens de la classification m\u00e9trologique de Stanley Smith Stevens. Appliquer machinalement la moyenne et l&rsquo;\u00e9cart-type sur des scores purement ordinaux, tels que des r\u00e9ponses isol\u00e9es sur une \u00e9chelle de Likert \u00e0 quatre modalit\u00e9s, constitue une violation m\u00e9thodologique qui vicie l&rsquo;estimation des proportions de dispersion.<\/p>\n<p>Par ailleurs, la sensibilit\u00e9 de la r\u00e8gle empirique aux anomalies distributionnelles est aigu\u00eb. Deux descripteurs morphologiques de forme doivent faire l&rsquo;objet d&rsquo;un examen vigilant : le coefficient d&rsquo;asym\u00e9trie (ou <em>skewness<\/em>) et le coefficient d&rsquo;aplatissement (ou <em>kurtosis<\/em>). La moindre asym\u00e9trie d\u00e9truit l&rsquo;\u00e9quiprobabilit\u00e9 des queues de distribution ; si un \u00e9chantillon pr\u00e9sente une asym\u00e9trie positive marqu\u00e9e, l&rsquo;intervalle moyen plus ou moins un \u00e9cart-type n&rsquo;englobera plus 68 % des scores, mais concentrera artificiellement une proportion d\u00e9s\u00e9quilibr\u00e9e d&rsquo;observations dans sa moiti\u00e9 inf\u00e9rieure tout en \u00e9tirant de mani\u00e8re disproportionn\u00e9e sa queue sup\u00e9rieure.<\/p>\n<p>De surcro\u00eet, l&rsquo;architecture d&rsquo;aplatissement joue un r\u00f4le tout aussi critique. Une distribution leptocurtique, caract\u00e9ris\u00e9e par un pic central ac\u00e9r\u00e9 et des queues \u00e9paisses (exc\u00e8s de kurtosis positif), pr\u00e9sentera une concentration de valeurs autour de la moyenne bien sup\u00e9rieure aux 68,27 % th\u00e9oriques dans le premier \u00e9cart-type, tout en h\u00e9bergeant un contingent d&rsquo;observations extr\u00eames anormalement \u00e9lev\u00e9 au-del\u00e0 de trois \u00e9carts-types. \u00c0 l&rsquo;inverse, une distribution platycurtique, dont la morphologie est aplatie et les queues fines, sous-estimera l&rsquo;effectif contenu dans les paliers interm\u00e9diaires tout en r\u00e9duisant presque \u00e0 n\u00e9ant la survenue de scores d\u00e9viants.<\/p>\n<p>Enfin, l&rsquo;existence de troncatures d&rsquo;\u00e9chantillonnage repr\u00e9sente un \u00e9cueil classique en psychologie appliqu\u00e9e. Dans le cas de protocoles \u00e9valuant des pathologies s\u00e9v\u00e8res au sein de populations non cliniques, on assiste couramment \u00e0 des effets de plancher massifs : une majorit\u00e9 \u00e9crasante de participants obtient le score minimal de z\u00e9ro. Appliquer la r\u00e8gle des trois sigmas sur de telles donn\u00e9es tronqu\u00e9es rel\u00e8ve du non-sens statistique absolu, la variance calcul\u00e9e ne refl\u00e9tant plus la dispersion d&rsquo;un continuum latent mais la borne physique impos\u00e9e par le protocole exp\u00e9rimental.<\/p>\n<h2>2. Architecture math\u00e9matique et d\u00e9rivation de la distribution normale<\/h2>\n<h3>2.1 La fonction de densit\u00e9 de probabilit\u00e9 normale<\/h3>\n<p>La transcription formelle de la courbe de Gauss est r\u00e9gie par sa fonction de densit\u00e9 de probabilit\u00e9 continue, not\u00e9e <em>f(x)<\/em>. Cette \u00e9quation diff\u00e9rentie la probabilit\u00e9 ponctuelle &mdash; qui est nulle pour toute valeur exacte d&rsquo;un continuum &mdash; de la densit\u00e9 relative de probabilit\u00e9 en un point donn\u00e9 de la droite r\u00e9elle. Pour une variable al\u00e9atoire <em>X<\/em> caract\u00e9ris\u00e9e par une moyenne &mu; et une variance &sigma;<sup>2<\/sup>, la fonction analytique s&rsquo;\u00e9nonce selon l&rsquo;expression classique :<\/p>\n<p><em>f(x) = (1 \/ (&sigma; &radic;(2&pi;))) &times; exp(&minus;(x &minus; &mu;)<sup>2<\/sup> \/ (2&sigma;<sup>2<\/sup>))<\/em><\/p>\n<p>Au sein de cette \u00e9quation, les constantes &pi; et <em>e<\/em> jouent un r\u00f4le g\u00e9om\u00e9trique central. Le terme multiplicatif <em>1 \/ (&sigma; &radic;(2&pi;))<\/em> assure la fonction essentielle de normalisation : il garantit que l&rsquo;int\u00e9grale d\u00e9finie de la fonction sur l&rsquo;intervalle ouvert de moins l&rsquo;infini \u00e0 plus l&rsquo;infini soit rigoureusement \u00e9gale \u00e0 1, satisfaisant ainsi le premier axiome de Kolmogorov. L&rsquo;argument exponentiel contr\u00f4le la d\u00e9croissance quadratique de la fonction \u00e0 mesure que l&rsquo;on s&rsquo;\u00e9loigne du barycentre &mu;.<\/p>\n<p>Une propri\u00e9t\u00e9 g\u00e9om\u00e9trique fondamentale de cette courbe r\u00e9side dans la position exacte de ses points d&rsquo;inflexion, c&rsquo;est-\u00e0-dire les coordonn\u00e9es o\u00f9 la d\u00e9riv\u00e9e seconde de la fonction s&rsquo;annule, marquant le passage d&rsquo;une concavit\u00e9 tourn\u00e9e vers le bas \u00e0 une concavit\u00e9 tourn\u00e9e vers le haut. En calculant la d\u00e9riv\u00e9e seconde <em>f\u00a0\u00bb(x)<\/em> et en r\u00e9solvant l&rsquo;\u00e9quation <em>f\u00a0\u00bb(x) = 0<\/em>, on d\u00e9montre que ces points d&rsquo;inflexion se situent rigoureusement \u00e0 <em>x = &mu; &minus; &sigma;<\/em> et <em>x = &mu; + &sigma;<\/em>. Ainsi, le premier intervalle de la r\u00e8gle empirique correspond pr\u00e9cis\u00e9ment \u00e0 la zone morphologique o\u00f9 la courbe ralentit son effondrement vertical pour amortir sa courbure vers les queues asymptotiques.<\/p>\n<p>Pour obtenir la probabilit\u00e9 que la variable <em>X<\/em> prenne une valeur comprise dans un intervalle ferm\u00e9 <em>[a, b]<\/em>, il est math\u00e9matiquement imp\u00e9ratif de calculer l&rsquo;int\u00e9grale d\u00e9finie de cette fonction de densit\u00e9 entre ces deux bornes :<\/p>\n<p><em>P(a &le; X &le; b) = &int;<sub>a<\/sub><sup>b<\/sup> f(x) dx<\/em><\/p>\n<p>Toutefois, la fonction <em>exp(&minus;t<sup>2<\/sup>)<\/em> ne poss\u00e8de pas de primitive analytique pouvant s&rsquo;exprimer au moyen des fonctions \u00e9l\u00e9mentaires standard (polyn\u00f4mes, trigonom\u00e9trie, logarithmes). Cette impossibilit\u00e9 fondamentale, d\u00e9montr\u00e9e par Joseph Liouville au XIXe si\u00e8cle, impose le recours imp\u00e9ratif \u00e0 des techniques d&rsquo;analyse num\u00e9rique avanc\u00e9es ou \u00e0 des d\u00e9veloppements en s\u00e9ries de Taylor pour approximer la valeur de l&rsquo;aire sous la courbe avec une pr\u00e9cision arbitrairement \u00e9lev\u00e9e.<\/p>\n<h3>2.2 La fonction de r\u00e9partition cumulative (CDF)<\/h3>\n<p>Face \u00e0 l&rsquo;inexistence de primitive \u00e9l\u00e9mentaire, la formalisation des probabilit\u00e9s normales repose sur la fonction de r\u00e9partition cumulative, universellement d\u00e9sign\u00e9e par la lettre grecque majuscule &Phi;(z) pour la variable centr\u00e9e r\u00e9duite. Par d\u00e9finition, cette fonction exprime la probabilit\u00e9 non conditionnelle qu&rsquo;une variable normale standardis\u00e9e <em>Z<\/em> (ayant une moyenne &mu; = 0 et un \u00e9cart-type &sigma; = 1) prenne une valeur inf\u00e9rieure ou \u00e9gale \u00e0 un seuil donn\u00e9 <em>z<\/em> :<\/p>\n<p><em>&Phi;(z) = P(Z &le; z) = (1 \/ &radic;(2&pi;)) &times; &int;<sub>&minus;&infin;<\/sub><sup>z<\/sup> exp(&minus;t<sup>2<\/sup> \/ 2) dt<\/em><\/p>\n<p>Cette fonction cumulative est monotone strictement croissante et se d\u00e9ploie asymptotiquement de la limite 0 lorsque <em>z<\/em> tend vers moins l&rsquo;infini, jusqu&rsquo;\u00e0 la limite 1 lorsque <em>z<\/em> tend vers plus l&rsquo;infini. En vertu de la sym\u00e9trie sp\u00e9culaire parfaite de la distribution normale autour de z\u00e9ro, la fonction de r\u00e9partition poss\u00e8de la propri\u00e9t\u00e9 fondamentale de compl\u00e9mentarit\u00e9 suivante :<\/p>\n<p><em>&Phi;(&minus;z) = 1 &minus; &Phi;(z)<\/em><\/p>\n<p>Cette relation d&rsquo;\u00e9quivalence permet de d\u00e9river avec une grande \u00e9l\u00e9gance la probabilit\u00e9 associ\u00e9e \u00e0 n&rsquo;importe quel intervalle sym\u00e9trique born\u00e9 par <em>&minus;z<\/em> et <em>+z<\/em>. En appliquant les propri\u00e9t\u00e9s de base des mesures d&rsquo;intervalles probabilistes, nous pouvons \u00e9tablir l&rsquo;\u00e9galit\u00e9 structurelle :<\/p>\n<p><em>P(&minus;z &le; Z &le; z) = &Phi;(z) &minus; &Phi;(&minus;z) = &Phi;(z) &minus; (1 &minus; &Phi;(z)) = 2&Phi;(z) &minus; 1<\/em><\/p>\n<p>Historiquement, cette op\u00e9ration n\u00e9cessitait la consultation fastidieuse de volumineuses tables num\u00e9riques imprim\u00e9es en fin d&rsquo;ouvrages statistiques, dans lesquelles les valeurs de &Phi;(z) \u00e9taient calcul\u00e9es avec quatre \u00e0 six d\u00e9cimales au moyen d&rsquo;int\u00e9grations num\u00e9riques manuelles bas\u00e9es sur la m\u00e9thode de Simpson ou les polyn\u00f4mes orthogonaux de Gauss-Legendre. Dans les environnements computationnels contemporains, ces approximations sont ex\u00e9cut\u00e9es instantan\u00e9ment par des algorithmes optimis\u00e9s de haute pr\u00e9cision, tels que l&rsquo;algorithme d&rsquo;approximation polynomiale de Cody ou des approximations rationnelles de Chebyshev, garantissant une pr\u00e9cision jusqu&rsquo;au seizi\u00e8me chiffre apr\u00e8s la virgule.<\/p>\n<h3>2.3 Quantification exacte des probabilit\u00e9s sous la courbe<\/h3>\n<p>L&rsquo;appellation \u00ab r\u00e8gle 68-95-99,7 \u00bb rel\u00e8ve d&rsquo;une convention mn\u00e9motechnique commode, mais dissimule les v\u00e9ritables valeurs transcendantes issues du calcul infinit\u00e9simal. L&rsquo;\u00e9valuation exacte de ces seuils r\u00e9v\u00e8le la pr\u00e9cision sous-jacente du mod\u00e8le de Gauss-Laplace.<\/p>\n<p>Pour le premier \u00e9cart-type, l&rsquo;\u00e9valuation de <em>2&Phi;(1) &minus; 1<\/em> avec des outils d&rsquo;analyse \u00e0 haute r\u00e9solution d\u00e9cimale fournit une valeur de :<\/p>\n<p><em>P(&minus;1 &le; Z &le; 1) &asymp; 0,682689492137&#8230;<\/em><\/p>\n<p>Soit tr\u00e8s exactement 68,26895 % de la masse de probabilit\u00e9 globale. La part r\u00e9siduelle totale exclue de cet intervalle est de 31,73105 %, ce qui implique, par sym\u00e9trie bilat\u00e9rale, que pr\u00e9cis\u00e9ment 15,8655 % des observations se situent au-dessous de &mu; &minus; &sigma;, et un pourcentage rigoureusement \u00e9quivalent au-dessus de &mu; + &sigma;.<\/p>\n<p>Pour le deuxi\u00e8me palier d&rsquo;\u00e9cart-type, le calcul de l&rsquo;expression <em>2&Phi;(2) &minus; 1<\/em> livre :<\/p>\n<p><em>P(&minus;2 &le; Z &le; 2) &asymp; 0,954499736103&#8230;<\/em><\/p>\n<p>Ce r\u00e9sultat d\u00e9montre que la valeur arrondie de 95 % sous-estime l\u00e9g\u00e8rement la r\u00e9alit\u00e9 math\u00e9matique, le v\u00e9ritable pourcentage \u00e9tant d&rsquo;environ 95,45 %. Il en d\u00e9coule que le seuil standard de significativit\u00e9 clinique et statistique fix\u00e9 \u00e0 &alpha; = 0,05 bilat\u00e9ral (soit exactement 95,000 % de couverture centrale) ne correspond pas pr\u00e9cis\u00e9ment \u00e0 deux \u00e9carts-types entiers, mais au quantile bien connu de 1,959964 \u00e9cart-type. Le reliquat hors de l&rsquo;intervalle &plusmn;2&sigma; s&rsquo;\u00e9tablit \u00e0 4,5500 %, r\u00e9parti en deux queues sym\u00e9triques de 2,2750 % chacune.<\/p>\n<p>Enfin, pour le troisi\u00e8me palier de dispersion, l&rsquo;\u00e9valuation num\u00e9rique de <em>2&Phi;(3) &minus; 1<\/em> produit :<\/p>\n<p><em>P(&minus;3 &le; Z &le; 3) &asymp; 0,997300203937&#8230;<\/em><\/p>\n<p>Soit environ 99,7300 %. Seule une fraction minuscule de 0,26998 % de la population totale \u00e9chappe \u00e0 cet intervalle tridimensionnel, r\u00e9partie en deux zones d&rsquo;extr\u00e9mit\u00e9 \u00e9gales de 0,1350 % (approximativement un cas sur 740 dans chaque queue). Ces grandeurs d\u00e9montrent la radicalit\u00e9 de l&rsquo;effondrement probabiliste : franchir la fronti\u00e8re du troisi\u00e8me \u00e9cart-type n&rsquo;est pas un \u00e9v\u00e9nement ordinaire, mais l&rsquo;indicateur statistique d&rsquo;un profil exceptionnel ou d&rsquo;une anomalie m\u00e9trologique majeure.<\/p>\n<h2>3. Environnement de travail et outillage computationnel sous R<\/h2>\n<h3>3.1 Configuration de la session R et d\u00e9pendances n\u00e9cessaires<\/h3>\n<p>L&rsquo;utilisation de R pour mod\u00e9liser et v\u00e9rifier la r\u00e8gle empirique repose sur une double architecture : la puissance statistique brute fournie par les composants internes du langage et l&rsquo;\u00e9cosyst\u00e8me enrichi apport\u00e9 par les biblioth\u00e8ques tierces. Bien que la distribution standard de R contienne toutes les routines n\u00e9cessaires pour effectuer des calculs de densit\u00e9 et de probabilit\u00e9 cumulative, l&rsquo;adoption d&rsquo;un environnement de programmation moderne facilite l&rsquo;\u00e9laboration de pipelines d&rsquo;analyse reproductibles et esth\u00e9tiquement rigoureux.<\/p>\n<p>Pour mener \u00e0 bien ce travail, l&rsquo;installation et le chargement de la suite <a rel=\"noopener noreferrer\" href=\"https:\/\/www.tidyverse.org\/\">tidyverse<\/a> constituent une \u00e9tape pr\u00e9alable recommand\u00e9e. Cette collection coh\u00e9rente d&rsquo;outils, qui englobe notamment les packages <em>readr<\/em>, <em>dplyr<\/em>, <em>tidyr<\/em> et <em>ggplot2<\/em>, transforme la manipulation de donn\u00e9es en une s\u00e9rie d&rsquo;op\u00e9rations fluides et lisibles. Parall\u00e8lement, des packages sp\u00e9cialis\u00e9s en analyse psychologique et m\u00e9trologique s&rsquo;av\u00e8rent pr\u00e9cieux. Le package <em>psych<\/em> d\u00e9velopp\u00e9 par William Revelle fournit des descripteurs de synth\u00e8se enrichis, adapt\u00e9s aux sp\u00e9cificit\u00e9s des \u00e9chelles composites. Le package <em>moments<\/em> offre un acc\u00e8s imm\u00e9diat aux estimateurs formels de sym\u00e9trie et de voussure, tandis que le package <em>car<\/em> apporte des outils de diagnostic avanc\u00e9s pour tester la lin\u00e9arit\u00e9 et l&rsquo;homog\u00e9n\u00e9it\u00e9 des r\u00e9sidus.<\/p>\n<p>Dans la conduite de ces exp\u00e9rimentations num\u00e9riques, la gestion des param\u00e8tres globaux de session garantit l&rsquo;int\u00e9grit\u00e9 des r\u00e9sultats. La configuration du formatage num\u00e9rique via la commande globale d&rsquo;options permet d&rsquo;\u00e9tendre l&rsquo;affichage jusqu&rsquo;\u00e0 sept ou huit d\u00e9cimales, \u00e9vitant ainsi les arrondis pr\u00e9matur\u00e9s g\u00e9n\u00e9r\u00e9s par l&rsquo;interface par d\u00e9faut de la console. De m\u00eame, l&rsquo;initialisation d\u00e9terministe du g\u00e9n\u00e9rateur de nombres pseudo-al\u00e9atoires au moyen de l&rsquo;instruction <strong>set.seed()<\/strong> constitue une exigence de reproductibilit\u00e9 scientifique : elle permet de s&rsquo;assurer que les simulations de cohortes psychom\u00e9triques produisent des jeux de donn\u00e9es strictement identiques d&rsquo;une machine \u00e0 une autre.<\/p>\n<h3>3.2 Pr\u00e9sentation des primitives statistiques du noyau R<\/h3>\n<p>Le noyau de base du langage R incorpore une famille coh\u00e9rente et hautement vectoris\u00e9e de quatre fonctions fondamentales d\u00e9di\u00e9es \u00e0 l&rsquo;exploration de la distribution normale. Cette nomenclature quadripartite s&rsquo;organise selon un pr\u00e9fixe d&rsquo;une seule lettre qui d\u00e9termine la nature math\u00e9matique de l&rsquo;op\u00e9ration, suivi de la cha\u00eene racine invariable <em>norm<\/em>. La compr\u00e9hension lucide de leurs fonctions respectives est essentielle :<\/p>\n<ul>\n<li><strong>dnorm()<\/strong> : Cette primitive calcule la densit\u00e9 de probabilit\u00e9 au point <em>x<\/em>. Elle correspond rigoureusement \u00e0 l&rsquo;application de l&rsquo;\u00e9quation analytique de la courbe de Gauss d\u00e9crite pr\u00e9c\u00e9demment. Elle ne fournit pas une probabilit\u00e9 cumul\u00e9e, mais la hauteur g\u00e9om\u00e9trique de l&rsquo;ordonn\u00e9e sur la courbe pour une valeur donn\u00e9e en abscisse. Son utilit\u00e9 computationnelle principale r\u00e9side dans le tra\u00e7age continu de profils th\u00e9oriques au-dessus de distributions de fr\u00e9quences empiriques.<\/li>\n<li><strong>pnorm()<\/strong> : Il s&rsquo;agit de la fonction de r\u00e9partition cumulative &Phi;(x). Elle calcule la surface int\u00e9gr\u00e9e sous la courbe depuis la limite th\u00e9orique inf\u00e9rieure jusqu&rsquo;\u00e0 une coordonn\u00e9e quantile sp\u00e9cifique. C&rsquo;est l&rsquo;instrument de pr\u00e9dilection pour v\u00e9rifier la r\u00e8gle empirique, puisqu&rsquo;elle permet d&rsquo;\u00e9valuer directement les masses de probabilit\u00e9 circonscrites entre n&rsquo;importe quelle paire de bornes.<\/li>\n<li><strong>qnorm()<\/strong> : Fonction r\u00e9ciproque de la pr\u00e9c\u00e9dente, ou fonction quantile. \u00c9tant donn\u00e9 une probabilit\u00e9 cumul\u00e9e <em>p<\/em> comprise imp\u00e9rativement dans l&rsquo;intervalle ]0, 1[, <strong>qnorm()<\/strong> d\u00e9termine la coordonn\u00e9e <em>x<\/em> exacte telle que <em>pnorm(x) = p<\/em>. Elle est indispensable pour identifier les seuils critiques cliniques, par exemple le niveau de score isolant les 2,5 % les plus faibles d&rsquo;une population de r\u00e9f\u00e9rence.<\/li>\n<li><strong>rnorm()<\/strong> : Ce g\u00e9n\u00e9rateur pseudo-al\u00e9atoire produit des tirages ind\u00e9pendants et identiquement distribu\u00e9s selon une sp\u00e9cification gaussienne d\u00e9finie par sa moyenne et son \u00e9cart-type. S&rsquo;appuyant sur des algorithmes ultra-rapides comme la transformation de Box-Muller ou la m\u00e9thode du zigourrat, cette primitive permet de concevoir des simulations de Monte-Carlo pour \u00e9prouver la robustesse des mod\u00e8les psychom\u00e9triques face aux contraintes d&rsquo;\u00e9chantillonnage.<\/li>\n<\/ul>\n<h2>4. D\u00e9monstration formelle de la r\u00e8gle empirique avec la fonction pnorm()<\/h2>\n<h3>4.1 Syntaxe et param\u00e8tres op\u00e9rationnels de pnorm()<\/h3>\n<p>La fonction <strong>pnorm()<\/strong> encapsule une infrastructure computationnelle d&rsquo;une remarquable sophistication derri\u00e8re une interface syntaxique intuitive. Sa signature standard comprend cinq arguments fonctionnels dont la ma\u00eetrise conditionne la validit\u00e9 des calculs :<\/p>\n<p><em>pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)<\/em><\/p>\n<p>L&rsquo;argument central <strong>q<\/strong> repr\u00e9sente le vecteur de quantiles dont on d\u00e9sire quantifier l&rsquo;aire cumul\u00e9e. Les param\u00e8tres <strong>mean<\/strong> et <strong>sd<\/strong> caract\u00e9risent l&rsquo;esp\u00e9rance math\u00e9matique et l&rsquo;\u00e9cart-type de la distribution cible ; lorsqu&rsquo;ils sont omis, R applique par d\u00e9faut la configuration standardis\u00e9e (&mu; = 0, &sigma; = 1). Le quatri\u00e8me param\u00e8tre, <strong>lower.tail<\/strong>, accepte une valeur logique : s&rsquo;il est configur\u00e9 sur <em>TRUE<\/em> (valeur par d\u00e9faut), la fonction int\u00e8gre la surface depuis l&rsquo;extr\u00eame gauche, renvoyant <em>P(X &le; q)<\/em>. \u00c0 l&rsquo;inverse, si le statisticien le commute sur <em>FALSE<\/em>, l&rsquo;algorithme int\u00e8gre la queue sup\u00e9rieure droite, \u00e9valuant directement <em>P(X &gt; q)<\/em> sans n\u00e9cessiter une soustraction compl\u00e9mentaire manuelle \u00e0 l&rsquo;unit\u00e9, minimisant ainsi les impr\u00e9cisions d&rsquo;arrondi num\u00e9rique sur les valeurs extr\u00eames.<\/p>\n<p>Le param\u00e8tre <strong>log.p<\/strong> permet, en cas de commutation sur <em>TRUE<\/em>, de retourner les probabilit\u00e9s sous forme de logarithmes n\u00e9p\u00e9riens naturels. Cette option s&rsquo;av\u00e8re vitale lors du traitement d&rsquo;\u00e9v\u00e9nements excessivement improbables survenant au-del\u00e0 de cinq ou six \u00e9carts-types, o\u00f9 les probabilit\u00e9s directes approcheraient la limite de pr\u00e9cision flottante de la machine (ph\u00e9nom\u00e8ne de <em>floating-point underflow<\/em>). Enfin, la nature intrins\u00e8quement vectoris\u00e9e de cette primitive implique que si l&rsquo;argument <strong>q<\/strong> re\u00e7oit un vecteur de dix coordonn\u00e9es distinctes, <strong>pnorm()<\/strong> retournera instantan\u00e9ment un vecteur de dix probabilit\u00e9s calcul\u00e9es de mani\u00e8re concurrente.<\/p>\n<h3>4.2 Calcul pas \u00e0 pas des probabilit\u00e9s d&rsquo;intervalles standards<\/h3>\n<p>L&rsquo;\u00e9valuation empirique des fondements math\u00e9matiques sous R s&rsquo;op\u00e8re par la mise en \u0153uvre de soustractions entre probabilit\u00e9s cumulatives bilat\u00e9rales. L&rsquo;ex\u00e9cution pas \u00e0 pas de ces commandes dans la console d\u00e9montre la convergence exacte vers les constantes de la r\u00e8gle 68-95-99,7.<\/p>\n<p>Pour d\u00e9terminer l&rsquo;aire circonscrite dans le premier \u00e9cart-type de la loi normale centr\u00e9e r\u00e9duite, l&rsquo;op\u00e9ration consiste \u00e0 retrancher la probabilit\u00e9 cumul\u00e9e au seuil <em>z = &minus;1<\/em> de celle observ\u00e9e au seuil <em>z = 1<\/em> :<\/p>\n<p><em>pnorm(1) &#8211; pnorm(-1)<\/em><\/p>\n<p>L&rsquo;interrogation de l&rsquo;interpr\u00e9teur R renvoie imm\u00e9diatement la valeur 0,6826895. Si l&rsquo;on souhaite observer la d\u00e9composition interne de ce calcul, <strong>pnorm(1)<\/strong> produit 0,8413447, repr\u00e9sentant l&rsquo;ensemble de la masse comprise entre moins l&rsquo;infini et +1 \u00e9cart-type, tandis que <strong>pnorm(-1)<\/strong> isole les 0,1586553 r\u00e9siduels situ\u00e9s dans la queue inf\u00e9rieure. La diff\u00e9rence nette consolide le fait que plus des deux tiers d&rsquo;une cohorte normale se stabilisent obligatoirement \u00e0 l&rsquo;int\u00e9rieur de cet intervalle pivot.<\/p>\n<p>L&rsquo;extension de cette logique computationnelle au deuxi\u00e8me palier d&rsquo;\u00e9cart-type requiert l&rsquo;\u00e9valuation de l&rsquo;\u00e9cart sym\u00e9trique :<\/p>\n<p><em>pnorm(2) &#8211; pnorm(-2)<\/em><\/p>\n<p>Le r\u00e9sultat g\u00e9n\u00e9r\u00e9 par R atteint 0,9544997. L&rsquo;inspection des segments isol\u00e9s montre que la probabilit\u00e9 cumul\u00e9e jusqu&rsquo;\u00e0 +2 sigmas est de 0,9772499, alors que la queue basse n&rsquo;abrite que 0,02275013 de la distribution totale. La soustraction confirme la couverture d&rsquo;environ 95,45 % de la population.<\/p>\n<p>Enfin, l&rsquo;exploration du troisi\u00e8me palier de dispersion met en \u0153uvre la soustraction canonique suivante :<\/p>\n<p><em>pnorm(3) &#8211; pnorm(-3)<\/em><\/p>\n<p>La console affiche le coefficient flottant 0,9973002. En compl\u00e9ment, l&rsquo;interrogation directe de la queue r\u00e9siduelle d&rsquo;un seul c\u00f4t\u00e9 via l&rsquo;instruction <em>pnorm(3, lower.tail = FALSE)<\/em> renvoie pr\u00e9cis\u00e9ment 0,001349898. Ce calcul d\u00e9montre avec clart\u00e9 pourquoi la r\u00e8gle empirique constitue le thermom\u00e8tre statistique par excellence : tout score situ\u00e9 au-del\u00e0 de trois \u00e9carts-types repr\u00e9sente un fait clinique dont la fr\u00e9quence d&rsquo;apparition fortuite est inf\u00e9rieure \u00e0 1,35 cas pour mille observations sous l&rsquo;hypoth\u00e8se de normalit\u00e9.<\/p>\n<h3>4.3 G\u00e9n\u00e9ralisation \u00e0 des distributions non standardis\u00e9es<\/h3>\n<p>L&rsquo;une des propri\u00e9t\u00e9s structurales les plus remarquables de la distribution normale est son invariance par transformation affine. La r\u00e8gle empirique ne s&rsquo;applique pas uniquement aux variables normalis\u00e9es de moyenne nulle et de dispersion unitaire, mais se conserve sur l&rsquo;ensemble des distributions r\u00e9elles, quels que soient les param\u00e8tres d&rsquo;\u00e9chelle &mu; et &sigma;.<\/p>\n<p>Consid\u00e9rons, \u00e0 titre d&rsquo;illustration psychom\u00e9trique concr\u00e8te, l&rsquo;\u00e9chelle d&rsquo;efficience intellectuelle globale de la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Wechsler_Adult_Intelligence_Scale\">WAIS-IV<\/a>, caract\u00e9ris\u00e9e par une moyenne &mu; = 100 et un \u00e9cart-type &sigma; = 15. Pour tester si la r\u00e8gle empirique conserve son int\u00e9grit\u00e9 probabiliste sur les scores bruts composites de cet instrument, il suffit de configurer dynamiquement les arguments de la primitive R sans standardiser pr\u00e9alablement les valeurs :<\/p>\n<p>Pour le premier \u00e9cart-type, les bornes physiques s&rsquo;\u00e9tablissent \u00e0 <em>100 &minus; 15 = 85<\/em> et <em>100 + 15 = 115<\/em>. Le calcul ex\u00e9cut\u00e9 sous R s&rsquo;\u00e9crit :<\/p>\n<p><em>pnorm(115, mean = 100, sd = 15) &#8211; pnorm(85, mean = 100, sd = 15)<\/em><\/p>\n<p>La valeur retourn\u00e9e par le compilateur est strictement identique : 0,6826895. De la m\u00eame mani\u00e8re, le calcul de l&rsquo;intervalle \u00e0 deux \u00e9carts-types d\u00e9limit\u00e9 par les scores de QI 70 (seuil classique d&rsquo;alerte pour le d\u00e9ficit intellectuel) et 130 (seuil traditionnel pour le haut potentiel intellectuel) s&rsquo;exprime comme suit :<\/p>\n<p><em>pnorm(130, mean = 100, sd = 15) &#8211; pnorm(70, mean = 100, sd = 15)<\/em><\/p>\n<p>Le syst\u00e8me restitue imperturbablement 0,9544997. Ce r\u00e9sultat formalise l&rsquo;\u00e9quivalence num\u00e9rique absolue entre le travail op\u00e9r\u00e9 sur les scores bruts param\u00e9tr\u00e9s et celui conduit sur les scores Z r\u00e9duits. L&rsquo;automatisation de ce principe en R peut \u00eatre encapsul\u00e9e dans une formule o\u00f9 l&rsquo;on sp\u00e9cifie de mani\u00e8re dynamique le quantile sup\u00e9rieur sous la forme <em>mean + k * sd<\/em> et le quantile inf\u00e9rieur sous la forme <em>mean &#8211; k * sd<\/em>, d\u00e9montrant que la constante g\u00e9om\u00e9trique transcende les unit\u00e9s de mesure physiques du ph\u00e9nom\u00e8ne explor\u00e9.<\/p>\n<h2>5. Application concr\u00e8te sur des jeux de donn\u00e9es psychom\u00e9triques<\/h2>\n<h3>5.1 Pr\u00e9paration et inspection d&rsquo;un \u00e9chantillon psychologique simul\u00e9<\/h3>\n<p>Pour appr\u00e9hender le comportement empirique de ces lois de dispersion face aux contingences du recueil de donn\u00e9es, la simulation d&rsquo;une cohorte clinique standardis\u00e9e offre un cadre d&rsquo;exp\u00e9rimentation p\u00e9dagogique id\u00e9al. Imaginons un protocole de recherche \u00e9valuant l&rsquo;anxi\u00e9t\u00e9-trait au moyen d&rsquo;un inventaire clinique calibr\u00e9 sur un score moyen attendu de 50 points avec un \u00e9cart-type de population de 10 points. Au moyen du noyau computationnel de R, nous pouvons g\u00e9n\u00e9rer un \u00e9chantillon repr\u00e9sentatif de 1000 participants en invoquant l&rsquo;instruction :<\/p>\n<p><em>scores_anxiete &lt;- rnorm(n = 1000, mean = 50, sd = 10)<\/em><\/p>\n<p>Avant d&rsquo;engager toute v\u00e9rification distributionnelle, le contr\u00f4le d&rsquo;int\u00e9grit\u00e9 du vecteur num\u00e9rique constitue le premier geste m\u00e9thodologique du chercheur. Ce processus implique l&rsquo;inspection des valeurs manquantes potentielles \u00e0 l&rsquo;aide de la commande <em>anyNA(scores_anxiete)<\/em> et la v\u00e9rification de la classe structurelle via <em>is.numeric(scores_anxiete)<\/em>. L&rsquo;\u00e9valuation descriptive \u00e9l\u00e9mentaire est ensuite produite par l&rsquo;agr\u00e9gation de la moyenne arithm\u00e9tique d&rsquo;\u00e9chantillon <em>mean(scores_anxiete)<\/em> et de l&rsquo;\u00e9cart-type d&rsquo;\u00e9chantillon corrig\u00e9 des degr\u00e9s de libert\u00e9 <em>sd(scores_anxiete)<\/em>.<\/p>\n<p>Sous l&rsquo;effet des fluctuations stochastiques d&rsquo;\u00e9chantillonnage, la moyenne calcul\u00e9e ne sera pas rigoureusement \u00e9gale \u00e0 50,000, mais oscillera par exemple autour de 49,87 ou 50,12, tandis que l&rsquo;\u00e9cart-type empirique se stabilisera \u00e0 proximit\u00e9 de 9,94 ou 10,08. L&rsquo;ex\u00e9cution de la commande <em>summary(scores_anxiete)<\/em> fournit un diagnostic pr\u00e9liminaire de sym\u00e9trie en comparant la moyenne et la m\u00e9diane : un \u00e9cart minime entre ces deux param\u00e8tres de tendance centrale apporte un premier indice en faveur de l&rsquo;absence de biais de distribution directionnel dans notre cohorte clinique simul\u00e9e.<\/p>\n<h3>5.2 D\u00e9termination empirique des bornes d&rsquo;intervalles<\/h3>\n<p>D\u00e8s lors que les descripteurs de localisation et d&rsquo;\u00e9talement sont extraits de l&rsquo;\u00e9chantillon, la v\u00e9rification pragmatique de la r\u00e8gle empirique impose de d\u00e9couper le continuum des scores r\u00e9els en tranches de variabilit\u00e9 fond\u00e9es sur l&rsquo;\u00e9cart-type empirique calcul\u00e9. Ce processus commence par la vectorisation des seuils critiques :<\/p>\n<p>La borne inf\u00e9rieure \u00e0 un \u00e9cart-type est assign\u00e9e \u00e0 <em>b1_inf &lt;- mean(scores_anxiete) &#8211; sd(scores_anxiete)<\/em>, et son homologue sup\u00e9rieure \u00e0 <em>b1_sup &lt;- mean(scores_anxiete) + sd(scores_anxiete)<\/em>. Les couples de bornes pour deux et trois \u00e9carts-types sont d\u00e9termin\u00e9s de fa\u00e7on analogue en multipliant le terme de dispersion par deux et par trois. La puissance fonctionnelle de R permet d&rsquo;\u00e9valuer le nombre de sujets tombant \u00e0 l&rsquo;int\u00e9rieur de ces intervalles au moyen d&rsquo;op\u00e9rateurs bool\u00e9ens relationnels combin\u00e9s :<\/p>\n<p><em>dans_1_sigma &lt;- scores_anxiete &gt;= b1_inf &amp; scores_anxiete &le; b1_sup<\/em><\/p>\n<p>En appliquant la fonction de sommation <em>sum(dans_1_sigma)<\/em>, R effectue une conversion implicite des valeurs logiques <em>TRUE<\/em> en entiers unitaires, retournant le d\u00e9nombrement pr\u00e9cis des observations capt\u00e9es. La division de ce d\u00e9compte par l&rsquo;effectif global de l&rsquo;\u00e9chantillon, formalis\u00e9e par <em>mean(dans_1_sigma) * 100<\/em>, extrait instantan\u00e9ment le pourcentage observ\u00e9 r\u00e9el. Sur un \u00e9chantillon bien calibr\u00e9 de 1000 sujets, on observera usuellement des proportions empiriques de l&rsquo;ordre de 68,1 % ou 68,9 % pour le premier niveau, 95,2 % \u00e0 95,8 % pour le deuxi\u00e8me niveau, et 99,6 % \u00e0 99,9 % pour le troisi\u00e8me palier, corroborant le comportement attendu du mod\u00e8le de Gauss appliqu\u00e9 \u00e0 une cohorte finie.<\/p>\n<h3>5.3 Analyse de la convergence empirique selon la taille d&rsquo;\u00e9chantillon<\/h3>\n<p>L&rsquo;\u00e9cart syst\u00e9matique constat\u00e9 entre les proportions observ\u00e9es au sein d&rsquo;une cohorte et les constantes th\u00e9oriques parfaites (68,27 %, 95,45 %, 99,73 %) est directement modul\u00e9 par la taille de l&rsquo;\u00e9chantillon <em>N<\/em>. Cette dynamique illustre l&rsquo;action de la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Loi_des_grands_nombres\">loi des grands nombres<\/a>, selon laquelle les estimateurs d&rsquo;\u00e9chantillon convergent en probabilit\u00e9 vers les param\u00e8tres de la population \u00e0 mesure que le volume d&rsquo;observations tend vers l&rsquo;infini.<\/p>\n<p>Sur une cohorte clinique restreinte de taille <em>N = 30<\/em> &mdash; configuration fr\u00e9quente dans les essais pilotes neuropsychologiques &mdash; l&rsquo;erreur d&rsquo;\u00e9chantillonnage s&rsquo;av\u00e8re consid\u00e9rable. L&rsquo;exclusion accidentelle de deux ou trois individus atypiques due \u00e0 la faible taille de l&rsquo;\u00e9chantillon peut faire fluctuer la proportion observ\u00e9e dans le premier \u00e9cart-type entre 60 % et 76 %. Le calcul de l&rsquo;erreur relative d\u00e9montre que l&rsquo;impr\u00e9cision statistique peut atteindre des seuils incompatibles avec une standardisation rigoureuse.<\/p>\n<p>\u00c0 l&rsquo;inverse, si l&rsquo;on \u00e9tend l&rsquo;exp\u00e9rience \u00e0 une m\u00e9ga-cohorte \u00e9pid\u00e9miologique de taille <em>N = 10 000<\/em> participants au moyen d&rsquo;une it\u00e9ration de <em>rnorm()<\/em>, la variance d&rsquo;\u00e9chantillonnage se comprime de mani\u00e8re spectaculaire. L&rsquo;erreur standard des proportions, gouvern\u00e9e par l&rsquo;expression math\u00e9matique <em>&radic;(p(1 &minus; p) \/ N)<\/em>, devient infinit\u00e9simale. L&rsquo;\u00e9cart absolu entre la proportion simul\u00e9e et les valeurs canoniques s&rsquo;amenuise pour ne repr\u00e9senter que quelques centi\u00e8mes de pourcent. Cette stabilit\u00e9 asymptotique d\u00e9montre empiriquement pourquoi les \u00e9talonnages des tests psychom\u00e9triques nationaux exigent des cohortes de normalisation comprenant au minimum plusieurs centaines, voire plusieurs milliers de sujets, pour que les seuils d&rsquo;\u00e9cart-type acqui\u00e8rent une validit\u00e9 diagnostique p\u00e9renne.<\/p>\n<h2>6. Visualisation avanc\u00e9e de la r\u00e8gle empirique avec ggplot2<\/h2>\n<h3>6.1 Construction de la courbe de Gauss th\u00e9orique<\/h3>\n<p>La transcription graphique de la r\u00e8gle empirique exige un degr\u00e9 de pr\u00e9cision visuelle conforme aux standards acad\u00e9miques internationaux. La biblioth\u00e8que <a rel=\"noopener noreferrer\" href=\"https:\/\/ggplot2.tidyverse.org\/\">ggplot2<\/a>, fond\u00e9e sur la grammaire des graphiques de Leland Wilkinson, constitue le standard de r\u00e9f\u00e9rence pour produire des visualisations scientifiques \u00e9l\u00e9gantes sous R. Contrairement aux syst\u00e8mes graphiques conventionnels, ggplot2 s\u00e9pare la mod\u00e9lisation math\u00e9matique du rendu esth\u00e9tique des couches.<\/p>\n<p>Pour tracer une distribution gaussienne th\u00e9orique pure sans d\u00e9pendre d&rsquo;un \u00e9chantillon discret particulier, l&rsquo;instruction <strong>stat_function()<\/strong> offre une m\u00e9thode directe. En instanciant un canevas vierge d\u00e9fini sur un intervalle d&rsquo;abscisses \u00e9tendu, typiquement de &minus;4 \u00e0 +4 sur l&rsquo;\u00e9chelle standardis\u00e9e, la commande permet d&rsquo;invoquer analytiquement la primitive <strong>dnorm<\/strong> :<\/p>\n<p><em>ggplot(data.frame(x = c(-4, 4)), aes(x = x)) + stat_function(fun = dnorm, geom = \u00ab\u00a0line\u00a0\u00bb, size = 1.2, color = \u00ab\u00a0#2C3E50\u00a0\u00bb)<\/em><\/p>\n<p>Ce proc\u00e9d\u00e9 garantit que la courbe trac\u00e9e correspond \u00e0 la r\u00e9solution infinit\u00e9simale de la fonction de densit\u00e9, \u00e9vitant les artefacts de segmentation polygonale propres aux histogrammes de faible r\u00e9solution. Pour satisfaire aux exigences de publication de l&rsquo;American Psychological Association (APA 7e \u00e9dition), l&rsquo;environnement graphique doit \u00eatre \u00e9pur\u00e9 : suppression du quadrillage lourd d&rsquo;arri\u00e8re-plan au moyen du th\u00e8me pr\u00e9d\u00e9fini <strong>theme_classic()<\/strong>, typographie proportionn\u00e9e sans empattement, et labellisation explicite des axes cart\u00e9siens indiquant la \u00ab Densit\u00e9 de probabilit\u00e9 \u00bb sur l&rsquo;ordonn\u00e9e et le \u00ab Score standardis\u00e9 (Z) \u00bb sur l&rsquo;abscisse.<\/p>\n<h3>6.2 Ombrage s\u00e9lectif des aires sous la courbe<\/h3>\n<p>La communication visuelle de la r\u00e8gle empirique prend tout son sens lorsque les intervalles &plusmn;1&sigma;, &plusmn;2&sigma; et &plusmn;3&sigma; sont mis en valeur par des jeux d&rsquo;ombres chromatiques \u00e9tag\u00e9s. Dans l&rsquo;\u00e9cosyst\u00e8me de ggplot2, cette op\u00e9ration s&rsquo;effectue avec \u00e9l\u00e9gance via l&rsquo;adjonction de couches successives exploitant l&rsquo;argument <strong>geom = \u00ab\u00a0ribbon\u00a0\u00bb<\/strong> ou l&rsquo;appel r\u00e9p\u00e9t\u00e9 \u00e0 <strong>stat_function()<\/strong> coupl\u00e9 \u00e0 l&rsquo;argument spatial <strong>xlim<\/strong>.<\/p>\n<p>En d\u00e9finissant une palette chromatique ordonn\u00e9e &mdash; par exemple trois nuances d\u00e9grad\u00e9es de bleu &mdash; le chercheur configure trois rubans int\u00e9gr\u00e9s : le premier circonscrit entre &minus;3 et +3, le deuxi\u00e8me restreint de &minus;2 \u00e0 +2, et le troisi\u00e8me focalis\u00e9 au c\u0153ur du profil entre &minus;1 et +1. Le contr\u00f4le rigoureux du canal alpha de transparence (g\u00e9n\u00e9ralement configur\u00e9 entre 0,25 et 0,40) est crucial : il assure que les superpositions de pigments ne saturent pas l&rsquo;espace visuel et laissent transpara\u00eetre la trajectoire de la courbe directrice principale trac\u00e9e au premier plan.<\/p>\n<p>Pour parachever la lisibilit\u00e9 g\u00e9om\u00e9trique, l&rsquo;insertion de lignes de d\u00e9marcation verticales en pointill\u00e9s, au moyen de la g\u00e9om\u00e9trie <strong>geom_vline()<\/strong>, mat\u00e9rialise avec pr\u00e9cision les fronti\u00e8res discr\u00e8tes aux abscisses enti\u00e8res &minus;3, &minus;2, &minus;1, 0, 1, 2 et 3. Ces marqueurs verticaux relient visuellement les points d&rsquo;inflexion morphologiques aux zones de dispersion probabiliste, offrant une lisibilit\u00e9 imm\u00e9diate des segments de risque clinique pour les praticiens non statisticiens.<\/p>\n<h3>6.3 Superposition d&rsquo;un histogramme empirique et d&rsquo;annotations textuelles<\/h3>\n<p>L&rsquo;un des d\u00e9fis p\u00e9dagogiques et diagnostiques majeurs consiste \u00e0 confronter graphiquement un recueil de donn\u00e9es cliniques r\u00e9elles aux exigences id\u00e9alis\u00e9es de la courbe normale th\u00e9orique. Cette confrontation impose une manipulation rigoureuse des \u00e9chelles d&rsquo;ordonn\u00e9es. Si un histogramme classique affiche les effectifs bruts (fr\u00e9quences absolues), il ne peut pas \u00eatre superpos\u00e9 directement \u00e0 la fonction <strong>dnorm<\/strong>, dont l&rsquo;int\u00e9grale unitaire produit des ordonn\u00e9es d\u00e9cimales tr\u00e8s faibles.<\/p>\n<p>La solution technique consiste \u00e0 forcer l&rsquo;histogramme \u00e0 normaliser son \u00e9chelle en densit\u00e9 relative d&rsquo;\u00e9chantillon en renseignant la correspondance interne <strong>aes(y = after_stat(density))<\/strong> au sein de l&rsquo;op\u00e9rateur <strong>geom_histogram()<\/strong>. Une fois cette homog\u00e9n\u00e9it\u00e9 m\u00e9trique \u00e9tablie, le statisticien peut projeter concomitamment :<\/p>\n<ul>\n<li>L&rsquo;histogramme d&rsquo;\u00e9chantillon, segment\u00e9 en classes de largeur optimis\u00e9e selon la r\u00e8gle de Freedman-Diaconis pour \u00e9viter tout biais de lissage arbitraire.<\/li>\n<li>L&rsquo;estimateur de densit\u00e9 par noyau non param\u00e9trique via <strong>geom_density()<\/strong>, permettant de capter les d\u00e9viations morphologiques r\u00e9elles de la cohorte.<\/li>\n<li>La trajectoire id\u00e9ale de la densit\u00e9 gaussienne th\u00e9orique via <strong>stat_function()<\/strong>, agissant comme mod\u00e8le de r\u00e9f\u00e9rence normatif.<\/li>\n<\/ul>\n<p>Enfin, l&rsquo;incorporation d&rsquo;annotations textuelles contextuelles au moyen de l&rsquo;op\u00e9rateur <strong>annotate()<\/strong> finalise la planche graphique. En fixant des cha\u00eenes de caract\u00e8res pr\u00e9cises au-dessus de chaque compartiment d&rsquo;\u00e9cart-type &mdash; par exemple en inscrivant \u00ab ~68,3 % \u00bb, \u00ab ~95,5 % \u00bb et \u00ab ~99,7 % \u00bb associ\u00e9es \u00e0 des fl\u00e8ches directionnelles fines &mdash; le graphique acquiert une valeur d\u00e9monstrative imm\u00e9diate, indispensable \u00e0 l&rsquo;\u00e9laboration de bilans neuropsychologiques ou de publications scientifiques \u00e9valu\u00e9es par les pairs.<\/p>\n<h2>7. V\u00e9rification rigoureuse des postulats de normalit\u00e9<\/h2>\n<h3>7.1 M\u00e9thodes d&rsquo;inspection graphique de la distribution<\/h3>\n<p>L&rsquo;application pr\u00e9cipit\u00e9e de la r\u00e8gle des trois sigmas sans v\u00e9rification empirique pr\u00e9alable de la normalit\u00e9 de la distribution constitue une faute m\u00e9thodologique av\u00e9r\u00e9e. Avant d&rsquo;engager des calculs de d\u00e9tection ou de classification fond\u00e9s sur l&rsquo;\u00e9cart-type, le statisticien se doit de conduire un examen visuel approfondi de la morphologie de son \u00e9chantillon.<\/p>\n<p>L&rsquo;outil le plus sensible et informatif pour cette mission est incontestablement le diagramme quantile-quantile standardis\u00e9, ou <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Diagramme_quantile-quantile\">Q-Q plot<\/a>. Ce dispositif cartographie les quantiles empiriques observ\u00e9s dans le jeu de donn\u00e9es face aux quantiles th\u00e9oriques qui seraient produits par une distribution gaussienne parfaite de param\u00e8tres identiques. Sous R et l&rsquo;architecture ggplot2, la construction de cet art\u00e9fact s&rsquo;op\u00e8re au moyen de la combinaison conjointe de <strong>stat_qq()<\/strong> pour le nuage de points et de <strong>stat_qq_line()<\/strong> pour la droite de r\u00e9f\u00e9rence asymptotique.<\/p>\n<p>L&rsquo;interpr\u00e9tation clinique d&rsquo;un Q-Q plot repose sur l&rsquo;examen rigoureux de la trajectoire des observations par rapport \u00e0 cette diagonale id\u00e9ale \u00e0 45 degr\u00e9s :<\/p>\n<ul>\n<li>Si les donn\u00e9es proviennent d&rsquo;une population gaussienne, les points s&rsquo;alignent de fa\u00e7on rectiligne et continue sur l&rsquo;ensemble de la m\u00e9diane et des queues.<\/li>\n<li>Une incurvation globale en forme de vo\u00fbte signale une asym\u00e9trie marqu\u00e9e du ph\u00e9nom\u00e8ne sous-jacent.<\/li>\n<li>Une configuration sigmo\u00efdale en forme de \u00ab S \u00bb invers\u00e9 d\u00e9note des queues de distribution excessivement lourdes (leptocurticit\u00e9), avertissant imm\u00e9diatement le chercheur que la proportion r\u00e9elle d&rsquo;observations aberrantes exc\u00e9dera substantiellement le seuil des 0,27 % th\u00e9oriques dict\u00e9 par la r\u00e8gle empirique.<\/li>\n<\/ul>\n<p>En compl\u00e9ment, l&rsquo;exploitation conjointe de bo\u00eetes \u00e0 moustaches (<em>boxplots<\/em>) et de diagrammes en violon permet de mat\u00e9rialiser instantan\u00e9ment l&rsquo;\u00e9cart interquartile et d&rsquo;identifier la pr\u00e9sence de points isol\u00e9s au-del\u00e0 des lignes de Tukey, affinant ainsi le diagnostic morphologique pr\u00e9liminaire.<\/p>\n<h3>7.2 Tests statistiques formels d&rsquo;ad\u00e9quation<\/h3>\n<p>Si l&rsquo;exploration graphique renseigne efficacement sur la nature de la distribution, elle doit \u00eatre \u00e9tay\u00e9e par des tests d&rsquo;hypoth\u00e8se formels destin\u00e9s \u00e0 rejeter ou \u00e0 conforter l&rsquo;hypoth\u00e8se nulle d&rsquo;ad\u00e9quation gaussienne (<em>H0 : la distribution sous-jacente suit une loi normale<\/em>).<\/p>\n<p>Le test statistique le plus puissant et le plus couramment mobilis\u00e9 dans la pratique psychom\u00e9trique est le test de Shapiro-Wilk, accessible nativement dans R via la fonction <strong>shapiro.test()<\/strong>. Fond\u00e9 sur le rapport entre l&rsquo;estimateur des moindres carr\u00e9s g\u00e9n\u00e9ralis\u00e9s de la variance et l&rsquo;estimateur non biais\u00e9 standard, ce test d\u00e9montre une sensibilit\u00e9 remarquable aux distorsions intervenant au c\u0153ur de la distribution comme dans ses zones p\u00e9riph\u00e9riques. Cependant, le praticien doit rester conscient d&rsquo;une limitation algorithmique propre \u00e0 R : la fonction rejette syst\u00e9matiquement les vecteurs dont la taille est inf\u00e9rieure \u00e0 3 observations ou exc\u00e8de 5000 individus.<\/p>\n<p>Pour contourner cette borne sur les vastes cohortes cliniques, le package sp\u00e9cialis\u00e9 <em>nortest<\/em> met \u00e0 disposition des alternatives inf\u00e9rentielles de premier plan, notamment le test d&rsquo;Anderson-Darling accessible par l&rsquo;instruction <strong>ad.test()<\/strong>. Ce dernier applique un sch\u00e9ma de pond\u00e9ration quadratique sp\u00e9cifique conf\u00e9rant un poids accru aux queues de distribution, le rendant sup\u00e9rieur au test de Kolmogorov-Smirnov corrig\u00e9 par Lilliefors (<strong>lillie.test()<\/strong>) pour \u00e9valuer la fiabilit\u00e9 des bornes extr\u00eames associ\u00e9es aux deuxi\u00e8me et troisi\u00e8me \u00e9carts-types.<\/p>\n<p>Cependant, l&rsquo;interpr\u00e9tation des seuils de significativit\u00e9 <em>p<\/em> doit \u00eatre men\u00e9e avec discernement \u00e9pist\u00e9mologique. Les tests formels d&rsquo;ad\u00e9quation sont soumis au paradoxe de la sur-puissance statistique sur les grands \u00e9chantillons : d\u00e8s lors que l&rsquo;effectif <em>N<\/em> franchit le millier d&rsquo;individus, la moindre d\u00e9viation microscopique et sans port\u00e9e clinique r\u00e9elle vis-\u00e0-vis d&rsquo;une sym\u00e9trie parfaite induit m\u00e9caniquement une valeur de <em>p<\/em> inf\u00e9rieure au seuil fatidique de 0,05. Ainsi, le rejet formel de la normalit\u00e9 par un test statistique ne disqualifie pas n\u00e9cessairement l&rsquo;utilisation op\u00e9rationnelle de la r\u00e8gle empirique si les d\u00e9viations restent d&rsquo;amplitude n\u00e9gligeable.<\/p>\n<h3>7.3 \u00c9valuation quantitative de l&rsquo;asym\u00e9trie et de l&rsquo;aplatissement<\/h3>\n<p>Afin de temp\u00e9rer les arbitrages des tests formels sur grands \u00e9chantillons, la quantification de l&rsquo;asym\u00e9trie (<em>skewness<\/em>) et de la voussure (<em>kurtosis<\/em>) fournit des indicateurs morphologiques objectifs. L&rsquo;infrastructure computationnelle du package <em>moments<\/em> permet d&rsquo;extraire directement ces grandeurs via les instructions <strong>skewness()<\/strong> et <strong>kurtosis()<\/strong>.<\/p>\n<p>Pour une loi normale th\u00e9orique parfaite, le coefficient d&rsquo;asym\u00e9trie standardis\u00e9 de Fisher s&rsquo;\u00e9tablit rigoureusement \u00e0 z\u00e9ro. Dans la pratique de l&rsquo;\u00e9valuation psychologique et selon les conventions m\u00e9thodologiques admises (notamment celles document\u00e9es par Kline ou West), une valeur de <em>skewness<\/em> comprise dans l&rsquo;intervalle ferm\u00e9 [&minus;1,0, +1,0] atteste d&rsquo;une quasi-sym\u00e9trie compatible avec l&rsquo;application de la r\u00e8gle empirique. D\u00e8s lors que cette grandeur s&rsquo;\u00e9chappe au-del\u00e0 de &plusmn;2,0, l&rsquo;asym\u00e9trie est qualifi\u00e9e de s\u00e9v\u00e8re, invalidant formellement l&rsquo;utilisation des \u00e9carts-types sym\u00e9triques.<\/p>\n<p>En ce qui concerne l&rsquo;aplatissement, la commande standard de R renvoie l&rsquo;indice non centr\u00e9 de Pearson, dont la valeur de r\u00e9f\u00e9rence normale est \u00e9gale \u00e0 3,0 (contrairement \u00e0 l&rsquo;exc\u00e8s de kurtosis qui soustrait 3 pour atteindre une r\u00e9f\u00e9rence de z\u00e9ro). Une m\u00e9trique de Pearson substantiellement sup\u00e9rieure \u00e0 3 traduit une distribution leptocurtique aux queues \u00e9paisses, tandis qu&rsquo;une valeur inf\u00e9rieure caract\u00e9rise un patron platycurtique. En r\u00e8gle g\u00e9n\u00e9rale, une fourchette de kurtosis situ\u00e9e entre 2,0 et 4,0 pr\u00e9serve une approximation satisfaisante des seuils 68 % et 95 %.<\/p>\n<p>En pr\u00e9sence d&rsquo;anomalies morphologiques av\u00e9r\u00e9es, le chercheur dispose d&rsquo;alternatives de pr\u00e9traitement par transformation univari\u00e9e continue. L&rsquo;application d&rsquo;une transformation logarithmique <em>log(x + c)<\/em> ou par racine carr\u00e9e <em>sqrt(x)<\/em> permet fr\u00e9quemment de r\u00e9duire une asym\u00e9trie positive sur des temps de r\u00e9action cognitive. De mani\u00e8re plus g\u00e9n\u00e9rale, l&rsquo;algorithme d&rsquo;optimisation param\u00e9trique de Box-Cox, accessible dans la biblioth\u00e8que <em>MASS<\/em>, identifie automatiquement l&rsquo;exposant &lambda; permettant de stabiliser la variance et de rapprocher les donn\u00e9es brutes d&rsquo;une morphologie normale exploitable.<\/p>\n<h2>8. D\u00e9tection et traitement des valeurs aberrantes via la r\u00e8gle empirique<\/h2>\n<h3>8.1 Cadre conceptuel de l&rsquo;identification des outliers statistiques<\/h3>\n<p>L&rsquo;une des applications op\u00e9rationnelles les plus imm\u00e9diates de la r\u00e8gle empirique r\u00e9side dans la d\u00e9tection des valeurs aberrantes univari\u00e9es (les <em>outliers<\/em>). Sur le plan statistique, un score qui s&rsquo;\u00e9loigne excessivement du centre de gravit\u00e9 d&rsquo;une distribution normale poss\u00e8de une probabilit\u00e9 d&rsquo;occurrence spontan\u00e9e extr\u00eamement faible. En fixant arbitrairement une fronti\u00e8re de d\u00e9marcation au-del\u00e0 du seuil des trois \u00e9carts-types, le chercheur postule que toute observation v\u00e9rifiant <em>|X &minus; &mu;| &gt; 3&sigma;<\/em> a moins de 2,7 chances sur 1000 d&rsquo;\u00e9merger du processus stochastique standard.<\/p>\n<p>Il est fondamental d&rsquo;introduire une distinction clinique et m\u00e9thodologique nette entre les diverses natures d&rsquo;observations extr\u00eames :<\/p>\n<ul>\n<li>Les anomalies de saisie ou d&rsquo;instrumentation : erreurs typographiques lors du codage d&rsquo;un questionnaire, art\u00e9facts \u00e9lectriques sur un signal \u00e9lectroenc\u00e9phalographique ou d\u00e9faut de calibrage d&rsquo;un capteur de temps de r\u00e9ponse. Ces donn\u00e9es repr\u00e9sentent du bruit parasite et n&rsquo;ont aucune validit\u00e9 scientifique.<\/li>\n<li>Les cas cliniques authentiques : individus poss\u00e9dant un profil neuropsychologique singulier, une symptomatologie d&rsquo;une s\u00e9v\u00e9rit\u00e9 exceptionnelle ou des aptitudes cognitives hors normes. Ces valeurs appartiennent l\u00e9gitimement \u00e0 la population \u00e9tudi\u00e9e et leur rejet aveugle biaiserait l&rsquo;analyse \u00e9pid\u00e9miologique.<\/li>\n<\/ul>\n<p>La vuln\u00e9rabilit\u00e9 majeure de cette approche de d\u00e9tection tient \u00e0 la nature des estimateurs mobilis\u00e9s. La moyenne arithm\u00e9tique et l&rsquo;\u00e9cart-type sont des mesures dot\u00e9es d&rsquo;un point de rupture (<em>breakdown point<\/em>) nul : la pr\u00e9sence d&rsquo;une unique observation hautement aberrante suffit \u00e0 d\u00e9caler artificiellement la moyenne vers elle et \u00e0 gonfler d\u00e9mesur\u00e9ment l&rsquo;\u00e9cart-type. Ce ph\u00e9nom\u00e8ne d&rsquo;inflation des variances peut engendrer un effet d&rsquo;occultation (<em>masking<\/em>), o\u00f9 l&rsquo;amplitude gonfl\u00e9e de &sigma; emp\u00eache d&rsquo;autres valeurs d\u00e9viantes d&rsquo;atteindre le seuil critique d&rsquo;exclusion.<\/p>\n<h3>8.2 Impl\u00e9mentation d&rsquo;un algorithme de filtrage dans R<\/h3>\n<p>Malgr\u00e9 ces r\u00e9serves m\u00e9thodologiques, le filtrage fond\u00e9 sur la r\u00e8gle des trois sigmas demeure une proc\u00e9dure standard en psychom\u00e9trie exploratoire. L&rsquo;environnement R, conjugu\u00e9 aux verbes de manipulation du package <em>dplyr<\/em>, simplifie l&rsquo;extraction et l&rsquo;\u00e9tiquetage de ces points critiques au sein des tables de donn\u00e9es.<\/p>\n<p>Pour isoler les observations d\u00e9viantes au sein d&rsquo;une cohorte, la syntaxe op\u00e9rationnelle s&rsquo;appuie sur la g\u00e9n\u00e9ration de variables indicatrices calcul\u00e9es \u00e0 la vol\u00e9e. En exploitant la fonction <strong>mutate()<\/strong> combin\u00e9e \u00e0 l&rsquo;instruction conditionnelle <strong>case_when()<\/strong>, le statisticien peut formaliser la r\u00e8gle de d\u00e9tection vectoris\u00e9e suivante :<\/p>\n<p><em>donnees &lt;- donnees %&gt;% mutate(score_z = (score &#8211; mean(score, na.rm = TRUE)) \/ sd(score, na.rm = TRUE), statut_outlier = case_when(abs(score_z) &gt; 3 ~ \u00ab\u00a0Aberrant\u00a0\u00bb, TRUE ~ \u00ab\u00a0Normal\u00a0\u00bb))<\/em><\/p>\n<p>Ce script enrichit instantan\u00e9ment le tableau d&rsquo;analyse d&rsquo;un vecteur cat\u00e9goriel permettant de d\u00e9nombrer les cas probl\u00e9matiques via <em>count(statut_outlier)<\/em> ou de filtrer les identifiants concern\u00e9s au moyen du verbe <em>filter(statut_outlier == \u00ab\u00a0Aberrant\u00a0\u00bb)<\/em>.<\/p>\n<p>Sur le plan comparatif, il convient de confronter ce calibrage gaussien \u00e0 la c\u00e9l\u00e8bre m\u00e9thode non param\u00e9trique de John Tukey fond\u00e9e sur l&rsquo;\u00e9cart interquartile (<em>Interquartile Range<\/em> ou <em>IQR<\/em>). La r\u00e8gle de Tukey qualifie d&rsquo;aberrant tout point situ\u00e9 au-del\u00e0 de 1,5 fois l&rsquo;IQR en dehors du premier ou du troisi\u00e8me quartile. Sur une distribution rigoureusement normale, les barri\u00e8res de Tukey \u00e0 1,5 IQR correspondent approximativement \u00e0 &plusmn;2,698 \u00e9carts-types, isolant th\u00e9oriquement 0,70 % des observations, ce qui en fait un crible plus conservateur que le seuil des trois \u00e9carts-types. La cartographie crois\u00e9e des deux m\u00e9thodes sous ggplot2, en surlignant les outliers par une couleur contrast\u00e9e sur un nuage de points ou un histogramme, permet de documenter visuellement la coh\u00e9rence des deux approches.<\/p>\n<h3>8.3 Proc\u00e9d\u00e9s de traitement : exclusion, winsorisation ou imputation<\/h3>\n<p>D\u00e8s lors que des observations sont formellement isol\u00e9es au-del\u00e0 du seuil critique des trois \u00e9carts-types, trois options m\u00e9thodologiques distinctes s&rsquo;offrent au chercheur, chacune engageant sa responsabilit\u00e9 quant \u00e0 la validit\u00e9 des conclusions statistiques :<\/p>\n<p>La premi\u00e8re voie &mdash; l&rsquo;exclusion unilat\u00e9rale des donn\u00e9es (ou <em>trimming<\/em>) &mdash; consiste \u00e0 \u00e9liminer purement et simplement les lignes correspondantes du jeu d&rsquo;analyse. Si cette d\u00e9cision se justifie pleinement face \u00e0 une erreur instrumentale manifeste (tel qu&rsquo;un temps de r\u00e9action physiologique enregistr\u00e9 \u00e0 12 millisecondes, biologiquement impossible), elle doit \u00eatre proscrite s&rsquo;il s&rsquo;agit d&rsquo;une observation clinique authentique. Supprimer arbitrairement les sujets atypiques tronque artificiellement la variabilit\u00e9 humaine r\u00e9elle et sous-estime l&rsquo;erreur standard des estimations ult\u00e9rieures.<\/p>\n<p>La deuxi\u00e8me alternative, \u00e9l\u00e9gante et largement reconnue dans les sciences sociales, est la technique de <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Winsorisation\">winsorisation<\/a>. Conceptualis\u00e9e par le biostatisticien Charles Winsor, cette proc\u00e9dure consiste non pas \u00e0 supprimer les scores extr\u00eames, mais \u00e0 ramener leur valeur \u00e0 la borne maximale tol\u00e9r\u00e9e dans l&rsquo;intervalle normatif d\u00e9fini. Dans le cadre de la r\u00e8gle empirique, tout score sup\u00e9rieur \u00e0 &mu; + 3&sigma; se voit r\u00e9assign\u00e9 \u00e0 la valeur exacte de &mu; + 3&sigma;, et sym\u00e9triquement pour la borne inf\u00e9rieure. L&rsquo;impl\u00e9mentation sous R s&rsquo;op\u00e8re directement par l&rsquo;appel de la fonction <strong>Winsorize()<\/strong> issue du package <em>DescTools<\/em>, neutralisant l&rsquo;impact d\u00e9stabilisant des valeurs extr\u00eames sur la variance tout en pr\u00e9servant l&rsquo;effectif total <em>N<\/em> de l&rsquo;\u00e9chantillon.<\/p>\n<p>La troisi\u00e8me approche repose sur l&rsquo;imputation statistique. Les observations aberrantes av\u00e9r\u00e9es sont converties en valeurs manquantes (<em>NA<\/em>), puis r\u00e9estim\u00e9es au moyen d&rsquo;algorithmes robustes : imputation par la m\u00e9diane pour les d\u00e9marches univari\u00e9es simples, ou imputation multiple par \u00e9quations cha\u00een\u00e9es (package <em>mice<\/em>) exploitant les corr\u00e9lations inter-items pour restituer le score le plus probable du patient. Quelle que soit la strat\u00e9gie retenue, l&rsquo;\u00e9thique scientifique impose de conduire syst\u00e9matiquement une analyse de sensibilit\u00e9 comparant les r\u00e9sultats statistiques obtenus avant et apr\u00e8s traitement des outliers afin de garantir la transparence des inf\u00e9rences cliniques.<\/p>\n<h2>9. Automatisation et programmation fonctionnelle en R<\/h2>\n<h3>9.1 Conception d&rsquo;une fonction modulaire pour la r\u00e8gle empirique<\/h3>\n<p>L&rsquo;automatisation des flux de travail analytiques constitue l&rsquo;un des atouts majeurs de l&rsquo;\u00e9cosyst\u00e8me R. Plut\u00f4t que de r\u00e9p\u00e9ter manuellement des s\u00e9ries de commandes pour chaque variable d&rsquo;une batterie de tests, l&rsquo;\u00e9criture d&rsquo;une fonction personnalis\u00e9e robuste et r\u00e9utilisable optimise la productivit\u00e9 tout en \u00e9liminant les risques d&rsquo;erreur de copier-coller.<\/p>\n<p>Une fonction d\u00e9di\u00e9e \u00e0 la v\u00e9rification de la r\u00e8gle empirique doit accepter en entr\u00e9e un vecteur num\u00e9rique et comporter des m\u00e9canismes de contr\u00f4le d\u00e9fensifs (gestion des valeurs manquantes, validation de la nature num\u00e9rique des entr\u00e9es). Elle doit calculer dynamiquement la moyenne, l&rsquo;\u00e9cart-type, les seuils th\u00e9oriques et les pourcentages r\u00e9els capt\u00e9s, pour enfin retourner une structure tabulaire claire. La logique algorithmique de cette fonction peut \u00eatre structur\u00e9e comme suit :<\/p>\n<p>La fonction initialise d&rsquo;abord une validation d&rsquo;int\u00e9grit\u00e9 via l&rsquo;instruction <em>stopifnot(is.numeric(x))<\/em>. Elle extrait ensuite les statistiques descriptives essentielles en configurant syst\u00e9matiquement l&rsquo;argument <em>na.rm = TRUE<\/em>. Dans un second temps, elle calcule de mani\u00e8re vectoris\u00e9e les proportions d&rsquo;observations incluses dans les fen\u00eatres &plusmn;1&sigma;, &plusmn;2&sigma; et &plusmn;3&sigma;. Enfin, elle assemble ces donn\u00e9es dans un objet tibble ordonn\u00e9 r\u00e9capitulant, pour chaque palier de dispersion :<\/p>\n<ul>\n<li>Le niveau de dispersion analys\u00e9 (&plusmn;1&sigma;, &plusmn;2&sigma;, &plusmn;3&sigma;).<\/li>\n<li>La borne num\u00e9rique inf\u00e9rieure calcul\u00e9e.<\/li>\n<li>La borne num\u00e9rique sup\u00e9rieure calcul\u00e9e.<\/li>\n<li>Le pourcentage observ\u00e9 r\u00e9el au sein du vecteur fourni.<\/li>\n<li>Le pourcentage th\u00e9orique gaussien exact de r\u00e9f\u00e9rence.<\/li>\n<li>L&rsquo;\u00e9cart r\u00e9siduel absolu constatant le degr\u00e9 de conformit\u00e9 empirique.<\/li>\n<\/ul>\n<p>Pour parfaire l&rsquo;utilit\u00e9 clinique du script, l&rsquo;adjonction d&rsquo;un message d&rsquo;alerte contextuel d\u00e9clench\u00e9 par une condition d&rsquo;asym\u00e9trie (par exemple si la d\u00e9viation absolue d&rsquo;un palier exc\u00e8de 4 %) avertit instantan\u00e9ment l&rsquo;utilisateur d&rsquo;une violation flagrante des hypoth\u00e8ses de distribution, pr\u00e9venant ainsi toute erreur d&rsquo;interpr\u00e9tation diagnostique.<\/p>\n<h3>9.2 Vectorisation et ex\u00e9cution sur plusieurs variables simultan\u00e9es<\/h3>\n<p>En contexte psychom\u00e9trique r\u00e9el, le praticien ou le chercheur n&rsquo;analyse que rarement une variable isol\u00e9e. Une batterie neuropsychologique compl\u00e8te &mdash; comme la NEPSY-II ou une \u00e9chelle multidimensionnelle de personnalit\u00e9 &mdash; g\u00e9n\u00e8re fr\u00e9quemment des dizaines de sous-scores corr\u00e9l\u00e9s qu&rsquo;il convient de cartographier simultan\u00e9ment.<\/p>\n<p>L&rsquo;association de la programmation fonctionnelle issue du package <em>purrr<\/em> et des s\u00e9lecteurs de colonnes du package <em>dplyr<\/em> permet de d\u00e9ployer l&rsquo;\u00e9valuation de la r\u00e8gle empirique sur l&rsquo;ensemble d&rsquo;un classeur clinique en une unique ligne de commande expressive. Gr\u00e2ce au verbe <strong>summarise()<\/strong> combin\u00e9 \u00e0 la directive <strong>across()<\/strong>, l&rsquo;analyste peut sp\u00e9cifier un ensemble de pr\u00e9dicats :<\/p>\n<p><em>resultats_batterie &lt;- donnees_cliniques %&gt;% summarise(across(where(is.numeric), ~ list(calculer_regle_empirique(.x))))<\/em><\/p>\n<p>L&rsquo;utilisation de la fonction <strong>map_dfr()<\/strong> de <em>purrr<\/em> permet ensuite de \u00ab d\u00e9baller \u00bb ces listes imbriqu\u00e9es pour former un tableau de synth\u00e8se pr\u00eat pour l&rsquo;analyse. Cette architecture tabulaire peut \u00eatre imm\u00e9diatement format\u00e9e pour les rapports acad\u00e9miques gr\u00e2ce au package <em>gt<\/em> ou \u00e0 la commande <strong>kable()<\/strong> de <em>knitr<\/em>, produisant des tableaux conformes aux normes APA incluant l&rsquo;alignement d\u00e9cimal automatique et la mise en gras des d\u00e9viations pathologiques.<\/p>\n<p>Dans une perspective de d\u00e9ploiement en laboratoire de recherche, la robustesse de cette infrastructure peut \u00eatre s\u00e9curis\u00e9e au moyen de tests unitaires formels programm\u00e9s avec le package <em>testthat<\/em>. Ces routines valident par exemple que la somme des pourcentages observ\u00e9s ne d\u00e9passe jamais 100 %, ou que la fonction r\u00e9agit conform\u00e9ment aux attentes lorsqu&rsquo;un vecteur constant (variance nulle) lui est inject\u00e9 en entr\u00e9e, garantissant ainsi l&rsquo;int\u00e9grit\u00e9 logicielle de l&rsquo;outil au fil des mises \u00e0 jour.<\/p>\n<h3>9.3 Cr\u00e9ation d&rsquo;un pipeline graphique reproductible<\/h3>\n<p>L&rsquo;aboutissement de cette d\u00e9marche d&rsquo;automatisation r\u00e9side dans l&rsquo;encapsulation de la cha\u00eene de production visuelle au sein d&rsquo;une fonction graphique param\u00e9trable. En fournissant simplement le nom d&rsquo;un jeu de donn\u00e9es et l&rsquo;\u00e9tiquette de la variable cible, ce pipeline g\u00e9n\u00e8re un document de synth\u00e8se visuelle complet pr\u00eat pour la publication.<\/p>\n<p>Cette fonction graphique standardis\u00e9e orchestre de fa\u00e7on s\u00e9quentielle :<\/p>\n<ul>\n<li>L&rsquo;extraction des descripteurs de localisation et d&rsquo;\u00e9chelle de la variable s\u00e9lectionn\u00e9e.<\/li>\n<li>Le calcul automatis\u00e9 des coordonn\u00e9es des rubans d&rsquo;ombrage correspondant aux trois seuils de la r\u00e8gle empirique.<\/li>\n<li>La g\u00e9n\u00e9ration de l&rsquo;histogramme empirique normalis\u00e9 en densit\u00e9, adoss\u00e9 \u00e0 la courbe de r\u00e9f\u00e9rence th\u00e9orique.<\/li>\n<li>L&rsquo;int\u00e9gration dynamique d&rsquo;\u00e9tiquettes textuelles renseignant la moyenne empirique, l&rsquo;\u00e9cart-type et les pourcentages r\u00e9els constat\u00e9s sur l&rsquo;\u00e9chantillon.<\/li>\n<\/ul>\n<p>L&rsquo;avantage technique d&rsquo;une telle encapsulation r\u00e9side dans la standardisation de l&rsquo;esth\u00e9tique sur l&rsquo;ensemble des graphiques d&rsquo;un projet clinique, assurant une coh\u00e9rence visuelle parfaite entre les diff\u00e9rents chapitres d&rsquo;un rapport neuropsychologique. Coupl\u00e9 \u00e0 l&rsquo;instruction d&rsquo;exportation <strong>ggsave()<\/strong>, ce pipeline enregistre les visualisations en formats vectoriels haute r\u00e9solution (PDF ou SVG) sans aucune pixellisation, garantissant une int\u00e9gration fluide dans les environnements de r\u00e9daction dynamique tels que Quarto ou R Markdown.<\/p>\n<h2>10. \u00c9tude de cas psychologique : Analyse comparative d&rsquo;un inventaire d\u00e9pressif<\/h2>\n<h3>10.1 Description du protocole et de la cohorte d&rsquo;\u00e9tude<\/h3>\n<p>Afin d&rsquo;ancrer ces d\u00e9veloppements dans la r\u00e9alit\u00e9 clinique de l&rsquo;\u00e9valuation psychologique, examinons une \u00e9tude de cas portant sur l&rsquo;inventaire de d\u00e9pression de Beck de deuxi\u00e8me \u00e9dition (<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Inventaire_de_d%C3%A9pression_de_Beck\">BDI-II<\/a>). Cet auto-questionnaire r\u00e9put\u00e9 comprend 21 items \u00e9valuant l&rsquo;intensit\u00e9 de la symptomatologie d\u00e9pressive, g\u00e9n\u00e9rant un score composite continu born\u00e9 th\u00e9oriquement entre 0 et 63 points. Selon les conventions cliniques d&rsquo;Aaron Beck, les scores de 0 \u00e0 13 d\u00e9notent une d\u00e9pression minimale, de 14 \u00e0 19 une d\u00e9pression l\u00e9g\u00e8re, de 20 \u00e0 28 une d\u00e9pression mod\u00e9r\u00e9e, et de 29 \u00e0 63 une symptomatologie s\u00e9v\u00e8re.<\/p>\n<p>Notre protocole d&rsquo;investigation porte sur une cohorte de validation constitu\u00e9e de <em>N = 850<\/em> patients admis cons\u00e9cutivement au sein d&rsquo;un centre de consultations psychologiques sp\u00e9cialis\u00e9. L&rsquo;objectif m\u00e9trologique de la recherche consiste \u00e0 d\u00e9terminer si les scores globaux au BDI-II au sein de cette population symptomatique suivent une distribution suffisamment gaussienne pour autoriser l&rsquo;application de la r\u00e8gle empirique dans le diagnostic de s\u00e9v\u00e9rit\u00e9 individuelle.<\/p>\n<p>La formulation des hypoth\u00e8ses de travail repose sur une tension th\u00e9orique classique : bien que le th\u00e9or\u00e8me central limite s&rsquo;applique \u00e0 la somme des 21 items individuels de l&rsquo;\u00e9chelle, la nature m\u00eame d&rsquo;un \u00e9chantillon clinique admis en consultation biaise la distribution vers les scores \u00e9lev\u00e9s, cr\u00e9ant un risque d&rsquo;asym\u00e9trie n\u00e9gative ou d&rsquo;effets de concentration s\u00e9lective qui menaceraient la r\u00e9gularit\u00e9 des paliers 68-95-99,7.<\/p>\n<h3>10.2 Ex\u00e9cution int\u00e9grale du flux d&rsquo;analyse sous R<\/h3>\n<p>L&rsquo;ex\u00e9cution de la s\u00e9quence analytique d\u00e9bute par le chargement des biblioth\u00e8ques logicielles et l&rsquo;importation de la table de donn\u00e9es cliniques au moyen du package <em>readr<\/em>. Apr\u00e8s v\u00e9rification du nettoyage \u00e9l\u00e9mentaire des donn\u00e9es via <em>tidyr::drop_na()<\/em>, le profil num\u00e9rique des scores globaux au BDI-II est quantifi\u00e9 :<\/p>\n<p>L&rsquo;estimation des moments fournit une moyenne d&rsquo;\u00e9chantillon <em>m = 26,45<\/em> points et un \u00e9cart-type <em>s = 7,82<\/em> points. Les bornes d&rsquo;intervalles associ\u00e9es \u00e0 la r\u00e8gle empirique se d\u00e9ploient par cons\u00e9quent selon les segments suivants :<\/p>\n<ul>\n<li>Intervalle &plusmn;1&sigma; : de 18,63 \u00e0 34,27 points.<\/li>\n<li>Intervalle &plusmn;2&sigma; : de 10,81 \u00e0 42,09 points.<\/li>\n<li>Intervalle &plusmn;3&sigma; : de 2,99 \u00e0 49,91 points.<\/li>\n<\/ul>\n<p>La confrontation de ces seuils avec les donn\u00e9es r\u00e9elles au moyen du filtrage logique vectoris\u00e9 r\u00e9v\u00e8le des r\u00e9sultats concrets : 70,82 % des patients de la cohorte se concentrent dans le premier intervalle, 94,35 % dans le deuxi\u00e8me, et 99,53 % dans le troisi\u00e8me intervalle de dispersion.<\/p>\n<p>L&rsquo;interpr\u00e9tation clinique de ces d\u00e9viations relatives s&rsquo;av\u00e8re particuli\u00e8rement instructive. La concentration de pr\u00e8s de 71 % des patients au sein du premier \u00e9cart-type (soit 2,55 % de plus que la norme th\u00e9orique de Gauss) combin\u00e9e \u00e0 une l\u00e9g\u00e8re sous-repr\u00e9sentation dans le deuxi\u00e8me palier trahit un profil mod\u00e9r\u00e9ment leptocurtique. Les patients orient\u00e9s vers ce centre sp\u00e9cialis\u00e9 pr\u00e9sentent une homog\u00e9n\u00e9it\u00e9 diagnostique autour des niveaux mod\u00e9r\u00e9s \u00e0 s\u00e9v\u00e8res (scores de 20 \u00e0 32), comprimant la variabilit\u00e9 centrale de l&rsquo;\u00e9chantillon par rapport \u00e0 une population tout-venant.<\/p>\n<h3>10.3 Synth\u00e8se et communication des r\u00e9sultats statistiques<\/h3>\n<p>La formalisation du compte-rendu d&rsquo;\u00e9valuation dans une publication ou un dossier m\u00e9dical informatis\u00e9 doit se conformer aux standards r\u00e9dactionnels de l&rsquo;APA. La transcription narrative de cette \u00e9tude de cas s&rsquo;articule comme suit :<\/p>\n<p>\u00ab L&rsquo;analyse distributionnelle des scores au BDI-II au sein de notre cohorte clinique (<em>N = 850<\/em>, <em>M = 26,45<\/em>, <em>SD = 7,82<\/em>) met en \u00e9vidence une ad\u00e9quation globale satisfaisante avec le mod\u00e8le de la r\u00e8gle empirique, malgr\u00e9 une l\u00e9g\u00e8re leptocurticit\u00e9 (kurtosis empirique = 3,42 ; skewness = &minus;0,18). Les taux de couverture r\u00e9els pour un, deux et trois \u00e9carts-types s&rsquo;\u00e9tablissent respectivement \u00e0 70,82 %, 94,35 % et 99,53 %, affichant des \u00e9carts absolus mod\u00e9r\u00e9s (&le; 2,55 %) vis-\u00e0-vis des constantes gaussiennes th\u00e9oriques de 68,27 %, 95,45 % et 99,73 %. \u00bb<\/p>\n<p>Sur le plan diagnostique individuel, ce mod\u00e8le offre un rep\u00e8re normatif pour positionner chaque patient. Par exemple, un patient obtenant un score brut de 43 points se situe au-del\u00e0 de la borne &mu; + 2&sigma; (fix\u00e9e \u00e0 42,09). Le clinicien peut affirmer avec un haut degr\u00e9 de certitude statistique que ce profil symptomatique ne concerne que moins de 2,5 % des consultants du service, justifiant une intensification imm\u00e9diate de la prise en charge th\u00e9rapeutique.<\/p>\n<p>Cette approche assist\u00e9e par R d\u00e9montre que si les postulats de distribution ne sont jamais parfaitement respect\u00e9s sur les \u00e9chantillons humains, la r\u00e8gle empirique fournit un cadre de calibrage op\u00e9ratoire d&rsquo;une grande valeur pragmatique, \u00e0 la condition expresse que le praticien conserve un regard critique sur les biais d&rsquo;orientation propres \u00e0 son contexte de recrutement.<\/p>\n<h2>11. Limites m\u00e9thodologiques et controverses en psychologie quantitative<\/h2>\n<h3>11.1 L&rsquo;illusion de la normalit\u00e9 dans les variables comportementales<\/h3>\n<p>L&rsquo;adh\u00e9sion inconditionnelle au mod\u00e8le gaussien en psychologie fait l&rsquo;objet de contestations m\u00e9thodologiques r\u00e9currentes depuis plusieurs d\u00e9cennies. Dans un article devenu classique intitul\u00e9 <em>The Micceri Fallacy<\/em>, le statisticien Theodore Micceri a proc\u00e9d\u00e9 \u00e0 l&rsquo;analyse empirique syst\u00e9matique de 440 distributions issues d&rsquo;\u00e9preuves d&rsquo;\u00e9valuation \u00e9ducative et psychologique de grande envergure. Ses conclusions furent sans appel : aucune des distributions \u00e9tudi\u00e9es ne pr\u00e9sentait une normalit\u00e9 rigoureuse, et la majorit\u00e9 exhibait des asym\u00e9tries prononc\u00e9es, des bimodalit\u00e9s sous-jacentes ou des queues \u00e0 d\u00e9croissance non exponentielle.<\/p>\n<p>L&rsquo;une des causes fondamentales de cette distorsion r\u00e9side dans l&rsquo;existence d&rsquo;effets de plancher et de plafond inh\u00e9rents \u00e0 la structure des \u00e9preuves. Dans les tests mesurant des comp\u00e9tences cognitives exceptionnelles ou des phobies sp\u00e9cifiques au sein de la population g\u00e9n\u00e9rale, une part pr\u00e9pond\u00e9rante des participants se regroupe de mani\u00e8re contrainte sur le score plancher z\u00e9ro ou sur le score plafond maximal, d\u00e9truisant l&rsquo;asymptotique gaussienne. Dans ces conditions m\u00e9trologiques, postuler qu&rsquo;environ 68 % des individus gravitent autour de la moyenne arithm\u00e9tique constitue une erreur d&rsquo;interpr\u00e9tation majeure, la moyenne arithm\u00e9tique perdant toute repr\u00e9sentativit\u00e9 face \u00e0 l&rsquo;asym\u00e9trie des observations.<\/p>\n<p>De surcro\u00eet, le recrutement traditionnel des \u00e9chantillons en psychologie universitaire &mdash; reposant historiquement sur des cohortes d&rsquo;\u00e9tudiants de premier cycle issues de soci\u00e9t\u00e9s dites <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/WEIRD\">WEIRD<\/a> (<em>Western, Educated, Industrialized, Rich, Democratic<\/em>) &mdash; induit une homog\u00e9n\u00e9isation artificielle des donn\u00e9es. Cette r\u00e9duction artificielle de la variance naturelle engendre une illusion de distribution normale au sein de l&rsquo;\u00e9chantillon d&rsquo;\u00e9tude, qui s&rsquo;effondre lorsque l&rsquo;outil est d\u00e9ploy\u00e9 sur des populations \u00e9pid\u00e9miologiques plus h\u00e9t\u00e9rog\u00e8nes.<\/p>\n<h3>11.2 Inad\u00e9quation face aux distributions multimodales et mixtes<\/h3>\n<p>La r\u00e8gle empirique postule formellement l&rsquo;unimodalit\u00e9 stricte de la variable continue. Or, de nombreuses variables comportementales et psychiatriques dissimulent des m\u00e9langes h\u00e9t\u00e9rog\u00e8nes de sous-populations sous une \u00e9tiquette diagnostique globale unique.<\/p>\n<p>En neuropsychiatrie, la mesure des temps de r\u00e9action psychomoteurs complexes illustre parfaitement cette impasse. La distribution de ces latences temporelles ne suit jamais une loi normale, mais adopte des profils de type ex-Gaussienne (convolution d&rsquo;une loi normale et d&rsquo;une loi exponentielle) ou de loi log-normale, caract\u00e9ris\u00e9es par une asym\u00e9trie droite massive due aux d\u00e9faillances attentionnelles ponctuelles. Appliquer m\u00e9caniquement le crit\u00e8re des trois sigmas sur de telles latences aboutit \u00e0 fixer une borne inf\u00e9rieure th\u00e9orique &mu; &minus; 3&sigma; souvent n\u00e9gative &mdash; ce qui constitue une impossibilit\u00e9 physique &mdash; tout en sous-estimant dramatiquement le volume r\u00e9el de temps de r\u00e9action anormalement lents dans la queue droite.<\/p>\n<p>De mani\u00e8re similaire, lorsqu&rsquo;un \u00e9chantillon clinique rassemble sans distinction des sujets sains et des patients atteints de troubles neurod\u00e9g\u00e9n\u00e9ratifs d\u00e9butants, la distribution sous-jacente forme un m\u00e9lange de densit\u00e9s (<em>mixture distribution<\/em>) potentiellement bimodal. L&rsquo;estimation d&rsquo;une moyenne et d&rsquo;un \u00e9cart-type combin\u00e9s globaux produit une moyenne qui tombe dans un \u00ab creux \u00bb statistique entre les deux populations, ne repr\u00e9sentant aucun patient r\u00e9el. Dans ce cas de figure, l&rsquo;intervalle &plusmn;1&sigma; englobe un pourcentage d&rsquo;individus nettement inf\u00e9rieur aux 68,27 % th\u00e9oriques, invalidant le triage clinique.<\/p>\n<h3>11.3 Divergence entre pertinence statistique et significativit\u00e9 clinique<\/h3>\n<p>L&rsquo;\u00e9cueil le plus insidieux dans l&rsquo;usage clinique des rep\u00e8res gaussiens r\u00e9side dans la confusion r\u00e9currente entre d\u00e9viance statistique normative et pertinence diagnostique individuelle. Le fait qu&rsquo;un sujet se situe \u00e0 plus de deux \u00e9carts-types de la moyenne d&rsquo;une population de r\u00e9f\u00e9rence n&rsquo;atteste en rien, de mani\u00e8re univoque, de la pr\u00e9sence d&rsquo;une psychopathologie ou d&rsquo;un d\u00e9ficit fonctionnel.<\/p>\n<p>Un score isol\u00e9 \u00e0 &minus;2,1&sigma; sur une \u00e9preuve de fluidit\u00e9 verbale peut simplement refl\u00e9ter une variabilit\u00e9 interindividuelle normale, un temp\u00e9rament r\u00e9serv\u00e9 ou un niveau de motivation temporairement affaibli lors de la passation, plut\u00f4t qu&rsquo;une atteinte l\u00e9sionnelle focale du lobe frontal. C&rsquo;est pourquoi la psychom\u00e9trie contemporaine subordonne la r\u00e8gle des sigmas \u00e0 des indices m\u00e9trologiques plus int\u00e9grateurs, tels que le changement clinique significatif (<em>Reliable Change Index<\/em> ou <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Reliable_change_index\">RCI<\/a> formalis\u00e9 par Jacobson et Truax) et l&rsquo;estimation de la taille d&rsquo;effet mesur\u00e9e par le d de Cohen.<\/p>\n<p>L&rsquo;indice RCI int\u00e8gre l&rsquo;erreur standard de mesure (<em>Standard Error of Measurement<\/em> ou <em>SEM<\/em>), calcul\u00e9e \u00e0 partir du coefficient de fid\u00e9lit\u00e9 test-retest de l&rsquo;instrument. Il rappelle utilement que l&rsquo;\u00e9cart-type d&rsquo;\u00e9chantillon amalgame la v\u00e9ritable variabilit\u00e9 biologique du trait et l&rsquo;impr\u00e9cision m\u00e9trologique du test. D\u00e8s lors, le jugement clinique \u00e9clair\u00e9 du praticien, adoss\u00e9 \u00e0 l&rsquo;anamn\u00e8se et \u00e0 l&rsquo;observation comportementale, doit pr\u00e9valoir sur les d\u00e9coupages automatis\u00e9s fournis par les logiciels d&rsquo;analyse statistique.<\/p>\n<h2>12. Alternatives non param\u00e9triques et approches computationnelles robustes<\/h2>\n<h3>12.1 L&rsquo;in\u00e9galit\u00e9 de Bienaym\u00e9-Tchebychev pour distributions quelconques<\/h3>\n<p>Lorsque les donn\u00e9es recueillies enfreignent de fa\u00e7on irr\u00e9ductible les postulats de normalit\u00e9, le chercheur n&rsquo;est pas d\u00e9muni pour autant. Il peut se tourner vers un th\u00e9or\u00e8me universel fondamental de la th\u00e9orie des probabilit\u00e9s : l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/In%C3%A9galit%C3%A9_de_Bienaym%C3%A9-Tchebychev\">in\u00e9galit\u00e9 de Bienaym\u00e9-Tchebychev<\/a>.<\/p>\n<p>Formul\u00e9 par Ir\u00e9n\u00e9e-Jules Bienaym\u00e9 et formalis\u00e9 par Pafnouti Tchebychev, ce th\u00e9or\u00e8me garantit que pour toute variable al\u00e9atoire poss\u00e9dant une moyenne &mu; et une variance finie &sigma;<sup>2<\/sup>, <em>quelle que soit la forme g\u00e9om\u00e9trique de sa distribution<\/em>, la proportion d&rsquo;observations s&rsquo;\u00e9cartant de la moyenne d&rsquo;au moins <em>k<\/em> \u00e9carts-types v\u00e9rifie l&rsquo;in\u00e9galit\u00e9 conservatrice universelle :<\/p>\n<p><em>P(|X &minus; &mu;| &ge; k&sigma;) &le; 1 \/ k<sup>2<\/sup><\/em><\/p>\n<p>Par compl\u00e9mentarit\u00e9, la masse minimale de donn\u00e9es obligatoirement contenue au sein de l&rsquo;intervalle centr\u00e9 &mu; &plusmn; k&sigma; s&rsquo;exprime par la borne :<\/p>\n<p><em>P(|X &minus; &mu;| &lt; k&sigma;) &ge; 1 &minus; (1 \/ k<sup>2<\/sup>)<\/em><\/p>\n<p>L&rsquo;application de ce th\u00e9or\u00e8me d\u00e9livre des garanties math\u00e9matiques strictes, ind\u00e9pendantes de toute hypoth\u00e8se de normalit\u00e9 :<\/p>\n<ul>\n<li>Pour <em>k = 1<\/em> \u00e9cart-type : la garantie math\u00e9matique est nulle (<em>1 &minus; 1 = 0 %<\/em>), ce qui traduit l&rsquo;existence possible de distributions o\u00f9 l&rsquo;ensemble des points se situent pr\u00e9cis\u00e9ment sur les bornes d&rsquo;un sigma.<\/li>\n<li>Pour <em>k = 2<\/em> \u00e9carts-types : la borne garantit qu&rsquo;au minimum <em>1 &minus; (1 \/ 4) = 75,0 %<\/em> des observations se trouvent imp\u00e9rativement dans l&rsquo;intervalle, loin des 95,45 % de la loi normale.<\/li>\n<li>Pour <em>k = 3<\/em> \u00e9carts-types : au minimum <em>1 &minus; (1 \/ 9) &asymp; 88,89 %<\/em> des donn\u00e9es sont n\u00e9cessairement circonscrites, face aux 99,73 % attendus sous hypoth\u00e8se gaussienne.<\/li>\n<\/ul>\n<p>Dans R, la comparaison directe entre les pourcentages th\u00e9oriques de Gauss, les observations empiriques r\u00e9elles et les bornes conservatrices de Tchebychev permet de cadrer l&rsquo;incertitude lors de l&rsquo;exploration de variables psychologiques fortement asym\u00e9triques ou plurimodales.<\/p>\n<h3>12.2 Statistiques non param\u00e9triques bas\u00e9es sur les quantiles<\/h3>\n<p>Face aux limites des estimateurs de dispersion param\u00e9triques, la statistique computationnelle non param\u00e9trique privil\u00e9gie les m\u00e9triques bas\u00e9es sur les quantiles et l&rsquo;ordonnancement des observations. Ces descripteurs poss\u00e8dent un niveau \u00e9lev\u00e9 de robustesse \u00e0 la rupture et s&rsquo;affranchissent totalement des pr\u00e9suppos\u00e9s morphologiques gaussiens.<\/p>\n<p>La primitive centrale sous R pour impl\u00e9menter cette d\u00e9marche est la fonction <strong>quantile()<\/strong>. Elle accepte en argument le vecteur d&rsquo;analyse et un vecteur de probabilit\u00e9s cibles <em>probs<\/em>, offrant neuf algorithmes distincts de calcul des quantiles continus (le type 7 \u00e9tant le standard par d\u00e9faut du logiciel). L&rsquo;\u00e9cart interquartile (<em>IQR<\/em>), extrait nativement par la commande <strong>IQR()<\/strong>, \u00e9value l&rsquo;\u00e9tendue couverte par les 50 % centraux de l&rsquo;\u00e9chantillon (diff\u00e9rence entre le 75e et le 25e percentile), fournissant un \u00e9quivalent direct et r\u00e9sistant \u00e0 l&rsquo;\u00e9cart-type classique.<\/p>\n<p>Parall\u00e8lement, la d\u00e9viation absolue m\u00e9diane (<em>Median Absolute Deviation<\/em> ou <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/%C3%89cart_m%C3%A9dian_absolu\">MAD<\/a>) constitue l&rsquo;estimateur de dispersion robuste par excellence. D\u00e9finie math\u00e9matiquement comme la m\u00e9diane des \u00e9carts absolus par rapport \u00e0 la m\u00e9diane de l&rsquo;\u00e9chantillon :<\/p>\n<p><em>MAD = median(|X<sub>i<\/sub> &minus; median(X)|)<\/em><\/p>\n<p>Dans l&rsquo;environnement de base de R, la fonction <strong>mad()<\/strong> int\u00e8gre par d\u00e9faut un facteur d&rsquo;ajustement asymptotique multiplicatif fix\u00e9 \u00e0 1,4826. Cette constante normalisatrice sp\u00e9cifique est choisie de telle sorte que, si la variable sous-jacente suit rigoureusement une loi normale, la valeur de la MAD estim\u00e9e co\u00efncide exactement avec l&rsquo;\u00e9cart-type classique &sigma;. Ainsi, la commande <em>mad()<\/em> permet d&rsquo;estimer l&rsquo;\u00e9chelle d&rsquo;une distribution sans subir la contamination des valeurs aberrantes extr\u00eames, offrant une alternative robuste pour d\u00e9finir des intervalles normatifs de dispersion clinique.<\/p>\n<h3>12.3 M\u00e9thodes de r\u00e9\u00e9chantillonnage par bootstrap<\/h3>\n<p>La troisi\u00e8me alternative moderne &mdash; rendue accessible par la puissance de calcul des processeurs contemporains &mdash; repose sur les m\u00e9thodes de r\u00e9\u00e9chantillonnage non param\u00e9trique par <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Bootstrap_(statistiques)\">bootstrap<\/a>, conceptualis\u00e9es par Bradley Efron \u00e0 la fin des ann\u00e9es 1970. Le principe directeur du bootstrap consiste \u00e0 traiter l&rsquo;\u00e9chantillon empirique initial comme une pseudo-population source, au sein de laquelle on effectue un grand nombre de tirages it\u00e9ratifs avec remise (typiquement <em>B = 2000<\/em> \u00e0 <em>10 000<\/em> r\u00e9p\u00e9titions) de taille <em>N<\/em> identique \u00e0 l&rsquo;\u00e9chantillon initial.<\/p>\n<p>L&rsquo;impl\u00e9mentation de cette strat\u00e9gie computationnelle sous R s&rsquo;appuie principalement sur le package standard <em>boot<\/em>. L&rsquo;utilisateur d\u00e9finit en premier lieu une fonction d&rsquo;extraction statistique renvoyant la m\u00e9trique d&rsquo;int\u00e9r\u00eat (par exemple les bornes d&rsquo;un intervalle de dispersion empirique donn\u00e9), puis d\u00e9l\u00e8gue au moteur <strong>boot()<\/strong> la g\u00e9n\u00e9ration de la distribution d&rsquo;\u00e9chantillonnage empirique. Cette m\u00e9thode permet de cartographier la pr\u00e9cision d&rsquo;un intervalle de confiance sans formuler la moindre hypoth\u00e8se a priori sur l&rsquo;\u00e9quation math\u00e9matique g\u00e9n\u00e9ratrice des donn\u00e9es.<\/p>\n<p>L&rsquo;atout majeur de cette biblioth\u00e8que r\u00e9side dans sa capacit\u00e9 \u00e0 d\u00e9river des intervalles de confiance de type <em>BCa<\/em> (<em>bias-corrected and accelerated<\/em>) via la fonction <strong>boot.ci()<\/strong>. Ces intervalles corrigent de mani\u00e8re non param\u00e9trique les asym\u00e9tries et les biais d&rsquo;acc\u00e9l\u00e9ration distributionnels observ\u00e9s dans les donn\u00e9es d&rsquo;origine. La comparaison finale entre l&rsquo;intervalle param\u00e9trique d\u00e9riv\u00e9 de la r\u00e8gle empirique de Gauss, la borne conservatively large garantie par Bienaym\u00e9-Tchebychev et les intervalles empiriques ajust\u00e9s issus du bootstrap dote le psychom\u00e8tre d&rsquo;un panorama analytique complet, alliant rigueur th\u00e9orique formelle et pragmatisme computationnel moderne.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7th ed.). American Psychological Association. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/li>\n<li>Beck, A. T., Steer, R. A., &amp; Brown, G. K. (1996). <em>Manual for the Beck Depression Inventory-II<\/em>. Psychological Corporation. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.pearsonassessments.com\">https:\/\/www.pearsonassessments.com<\/a><\/li>\n<li>Box, G. E. P., &amp; Cox, D. R. (1964). An analysis of transformations. <em>Journal of the Royal Statistical Society: Series B (Methodological)<\/em>, 26(2), 211-252. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x\">https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x<\/a><\/li>\n<li>Cantelli, F. P. (1910). Intorno ad un teorema fondamentale per la teoria degli errori. <em>Rendiconti del Circolo Matematico di Palermo<\/em>, 30(1), 321-326. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/BF03014884\">https:\/\/doi.org\/10.1007\/BF03014884<\/a><\/li>\n<li>Davison, A. C., &amp; Hinkley, D. V. (1997). <em>Bootstrap methods and their application<\/em>. Cambridge University Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1017\/CBO9780511802843\">https:\/\/doi.org\/10.1017\/CBO9780511802843<\/a><\/li>\n<li>Efron, B. (1979). Bootstrap methods: Another look at the jackknife. <em>The Annals of Statistics<\/em>, 7(1), 1-26. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aos\/1176344552\">https:\/\/doi.org\/10.1214\/aos\/1176344552<\/a><\/li>\n<li>Gauss, C. F. (1809). <em>Theoria motus corporum coelestium in sectionibus conicis solem ambientium<\/em>. Perthes et Besser. <a rel=\"noopener noreferrer\" href=\"https:\/\/archive.org\/details\/theoriamotuscorp00gausuoft\">https:\/\/archive.org\/details\/theoriamotuscorp00gausuoft<\/a><\/li>\n<li>Jacobson, N. S., &amp; Truax, P. (1991). Clinical significance: A statistical approach to defining meaningful change in psychotherapy research. <em>Journal of Consulting and Clinical Psychology<\/em>, 59(1), 12-19. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0022-006X.59.1.12\">https:\/\/doi.org\/10.1037\/0022-006X.59.1.12<\/a><\/li>\n<li>Micceri, T. (1989). The unicorn, the normal curve, and other improbable creatures. <em>Psychological Bulletin<\/em>, 105(1), 156-166. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0033-2909.105.1.156\">https:\/\/doi.org\/10.1037\/0033-2909.105.1.156<\/a><\/li>\n<li>Pearson, K. (1895). Contributions to the mathematical theory of evolution: II. Skew variation in homogeneous material. <em>Philosophical Transactions of the Royal Society of London. A<\/em>, 186, 343-414. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1098\/rsta.1895.0010\">https:\/\/doi.org\/10.1098\/rsta.1895.0010<\/a><\/li>\n<li>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><\/li>\n<li>Revelle, W. (2023). <em>psych: Procedures for psychological, psychometric, and personality research<\/em> (Version 2.3.9). Northwestern University. <a rel=\"noopener noreferrer\" href=\"https:\/\/CRAN.R-project.org\/package=psych\">https:\/\/CRAN.R-project.org\/package=psych<\/a><\/li>\n<li>Shapiro, S. S., &amp; Wilk, M. B. (1965). An analysis of variance test for normality (complete samples). <em>Biometrika<\/em>, 52(3\/4), 591-611. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2333709\">https:\/\/doi.org\/10.2307\/2333709<\/a><\/li>\n<li>Tchebychev, P. (1867). Des valeurs moyennes. <em>Journal de math\u00e9matiques pures et appliqu\u00e9es<\/em>, 2e s\u00e9rie, 12, 177-184. <a rel=\"noopener noreferrer\" href=\"https:\/\/eudml.org\/doc\/234479\">https:\/\/eudml.org\/doc\/234479<\/a><\/li>\n<li>Tukey, J. W. (1977). <em>Exploratory data analysis<\/em>. Addison-Wesley. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.pearson.com\">https:\/\/www.pearson.com<\/a><\/li>\n<li>Wechsler, D. (2008). <em>Wechsler Adult Intelligence Scale \u2013 Fourth Edition (WAIS-IV)<\/em>. Psychological Corporation. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.pearsonassessments.com\">https:\/\/www.pearsonassessments.com<\/a><\/li>\n<li>Wickham, H. (2016). <em>ggplot2: Elegant graphics for data analysis<\/em>. Springer-Verlag. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-3-319-24277-4\">https:\/\/doi.org\/10.1007\/978-3-319-24277-4<\/a><\/li>\n<li>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><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour comprendre, calculer et visualiser la r\u00e8gle empirique (68-95-99,7) dans R \u00e0 l&rsquo;aide de pnorm() et ggplot2 en psychom\u00e9trie.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2193,2],"tags":[1953,2688,16,2658,1951],"class_list":["post-4887","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-programmation-r","category-psychometrie","tag-distribution-normale","tag-ggplot2","tag-psychometrie","tag-r","tag-regle-empirique"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/4887","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=4887"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/4887\/revisions"}],"predecessor-version":[{"id":16703,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/4887\/revisions\/16703"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=4887"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=4887"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=4887"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}