{"id":1669,"date":"2026-09-05T10:01:43","date_gmt":"2026-09-05T10:01:43","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/guide-dt-qt-pt-rt-dans-r\/"},"modified":"2026-09-05T10:01:45","modified_gmt":"2026-09-05T10:01:45","slug":"guide-dt-qt-pt-rt-dans-r","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=guide-dt-qt-pt-rt-dans-r","title":{"rendered":"Un guide pour dt, qt, pt et rt dans R"},"content":{"rendered":"<p>L&rsquo;inf\u00e9rence statistique dans les sciences du comportement et les neurosciences cognitives repose sur un \u00e9quilibre d\u00e9licat entre les postulats math\u00e9matiques des mod\u00e8les probabilistes et les contraintes empiriques de la collecte de donn\u00e9es. Lorsque les chercheurs travaillent avec des \u00e9chantillons restreints, l&rsquo;estimation des param\u00e8tres de la population, en particulier la variance, devient particuli\u00e8rement instable. C&rsquo;est pr\u00e9cis\u00e9ment pour surmonter cette instabilit\u00e9 que la loi t de Student a \u00e9t\u00e9 introduite dans le corpus statistique moderne. Dans l&rsquo;\u00e9cosyst\u00e8me de programmation statistique R, la manipulation computationnelle de cette distribution cardinale ne s&rsquo;op\u00e8re pas au moyen d&rsquo;une proc\u00e9dure monolithique et opaque, mais via un quatuor fonctionnel remarquablement \u00e9l\u00e9gant et modulaire constitu\u00e9 par les fonctions <strong>dt<\/strong>, <strong>pt<\/strong>, <strong>qt<\/strong> et <strong>rt<\/strong>.<\/p>\n<p>Ce guide encyclop\u00e9dique propose une immersion exhaustive dans les fondements th\u00e9oriques, les m\u00e9canismes algorithmiques et les applications exp\u00e9rimentales de ces quatre fonctions au sein de l&rsquo;environnement de calcul R. Con\u00e7u sp\u00e9cifiquement pour les chercheurs, statisticiens et \u00e9tudiants avanc\u00e9s en psychologie quantitative, cet ouvrage d\u00e9taille chaque facette du traitement probabiliste de la loi t, depuis le calcul de densit\u00e9s infinit\u00e9simales jusqu&rsquo;\u00e0 la mod\u00e9lisation stochastique par les m\u00e9thodes de Monte-Carlo, en passant par le calcul exact des valeurs p et l&rsquo;estimation d&rsquo;intervalles de confiance non centraux. En abordant \u00e0 la fois la th\u00e9orie sous-jacente et les scripts d&rsquo;impl\u00e9mentation op\u00e9rationnels, ce document \u00e9tablit un pont m\u00e9thodologique indispensable entre la rigueur math\u00e9matique et les imp\u00e9ratifs contemporains de reproductibilit\u00e9 empirique.<\/p>\n<p>\u00c0 travers une exploration syst\u00e9matique des douze sections th\u00e9matiques qui composent cette \u00e9tude, le lecteur d\u00e9couvrira comment ma\u00eetriser la g\u00e9om\u00e9trie des courbes de distribution, comment d\u00e9jouer les pi\u00e8ges classiques d&rsquo;arrondi num\u00e9rique ou d&rsquo;interpr\u00e9tation fallacieuse de la densit\u00e9, et comment exploiter toute la puissance du param\u00e8tre de non-centralit\u00e9 pour les \u00e9tudes de puissance statistique a priori. La ma\u00eetrise approfondie de ce socle d&rsquo;outils analytiques constitue l&rsquo;un des piliers m\u00e9thodologiques fondamentaux pour quiconque souhaite mener des analyses inf\u00e9rentielles robustes, conformes aux exigences les plus strictes de la psychologie scientifique contemporaine.<\/p>\n<h2>1. Introduction th\u00e9orique \u00e0 la loi t de Student en psychologie statistique<\/h2>\n<h3>1.1 Origine historique et fondements math\u00e9matiques de William Gosset<\/h3>\n<p>L&rsquo;\u00e9mergence historique de la loi t de Student est intimement li\u00e9e aux exigences concr\u00e8tes du contr\u00f4le de la qualit\u00e9 industrielle au tournant du vingti\u00e8me si\u00e8cle. C&rsquo;est en 1908 que le chimiste et math\u00e9maticien britannique William Sealy Gosset, alors employ\u00e9 par la brasserie Guinness \u00e0 Dublin, publia sous le pseudonyme d\u00e9sormais universel de \u00ab Student \u00bb son article s\u00e9minal intitul\u00e9 <em>The Probable Error of a Mean<\/em> dans la revue <a rel=\"noopener noreferrer\" href=\"https:\/\/academic.oup.com\/biometrika\">Biometrika<\/a>. La direction de l&rsquo;entreprise interdisait formellement \u00e0 ses scientifiques de divulguer des secrets industriels ou de publier des travaux acad\u00e9miques sous leur v\u00e9ritable identit\u00e9, afin de pr\u00e9server l&rsquo;avantage concurrentiel li\u00e9 aux proc\u00e9d\u00e9s biochimiques de brassage. Gosset \u00e9tait confront\u00e9 \u00e0 un probl\u00e8me m\u00e9thodologique r\u00e9current : comment \u00e9valuer rigoureusement la qualit\u00e9 des mati\u00e8res premi\u00e8res, telles que le houblon ou l&rsquo;orge, \u00e0 partir d&rsquo;\u00e9chantillons drastiquement r\u00e9duits, souvent limit\u00e9s \u00e0 trois ou quatre mesures seulement, l\u00e0 o\u00f9 la th\u00e9orie des grands \u00e9chantillons d\u00e9velopp\u00e9e par Karl Pearson s&rsquo;av\u00e9rait dramatiquement inapplicable ?<\/p>\n<p>Sur le plan math\u00e9matique, Gosset cherchait \u00e0 d\u00e9terminer la distribution de la moyenne d&rsquo;un \u00e9chantillon lorsque la variance de la population g\u00e9n\u00e9rale est inconnue et doit \u00eatre estim\u00e9e \u00e0 partir des donn\u00e9es observ\u00e9es. Dans le cadre de la th\u00e9orie classique d\u00e9riv\u00e9e du th\u00e9or\u00e8me central limite, le quotient de la diff\u00e9rence entre la moyenne de l&rsquo;\u00e9chantillon et la moyenne th\u00e9orique par l&rsquo;erreur type repose sur l&rsquo;hypoth\u00e8se que l&rsquo;\u00e9cart-type de la population est une constante connue. D\u00e8s lors que cette grandeur inconnue est substitu\u00e9e par l&rsquo;\u00e9cart-type de l&rsquo;\u00e9chantillon, le ratio statistique ne suit plus une loi normale centr\u00e9e r\u00e9duite, mais une variable al\u00e9atoire composite. Cette grandeur est d\u00e9finie par le quotient d&rsquo;une variable normale standardis\u00e9e par la racine carr\u00e9e d&rsquo;une variable distribu\u00e9e selon une loi du khi-deux ind\u00e9pendante, elle-m\u00eame divis\u00e9e par ses degr\u00e9s de libert\u00e9 correspondants. Cette formalisation analytique a permis de r\u00e9soudre d\u00e9finitivement le probl\u00e8me de l&rsquo;impr\u00e9cision inh\u00e9rente \u00e0 l&rsquo;estimation conjointe de deux param\u00e8tres gaussiens sur de faibles volumes de donn\u00e9es.<\/p>\n<p>La transition de cette avanc\u00e9e industrielle vers les sciences humaines et comportementales s&rsquo;est op\u00e9r\u00e9e quelques d\u00e9cennies plus tard, principalement sous l&rsquo;impulsion m\u00e9thodologique de Ronald Fisher, qui a g\u00e9n\u00e9ralis\u00e9 et formalis\u00e9 les travaux de Gosset dans ses trait\u00e9s fondamentaux. En psychophysique et en psychom\u00e9trie, o\u00f9 les protocoles exp\u00e9rimentaux imposaient des t\u00e2ches perceptives longues, co\u00fbteuses et \u00e9prouvantes pour les participants, la possibilit\u00e9 de conduire des inf\u00e9rences rigoureuses sans recourir \u00e0 des cohortes de centaines de sujets a profond\u00e9ment boulevers\u00e9 la pratique de recherche. L&rsquo;adoption progressive du ratio t a ainsi permis d&rsquo;asseoir la psychologie exp\u00e9rimentale naissante sur des bases inf\u00e9rentielles solides, en l\u00e9gitimant math\u00e9matiquement l&rsquo;extrapolation de r\u00e9gularit\u00e9s cognitives universelles \u00e0 partir d&rsquo;\u00e9chantillons de taille modeste mais rigoureusement contr\u00f4l\u00e9s.<\/p>\n<h3>1.2 Pertinence de la loi de Student face aux petits \u00e9chantillons en sciences cognitives<\/h3>\n<p>Dans le domaine contemporain des neurosciences cognitives et de la psychologie exp\u00e9rimentale, la contrainte des petits \u00e9chantillons demeure une r\u00e9alit\u00e9 m\u00e9thodologique omnipr\u00e9sente. Que ce soit lors de protocoles d&rsquo;imagerie par r\u00e9sonance magn\u00e9tique fonctionnelle (IRMf) impliquant des populations cliniques sp\u00e9cifiques, d&rsquo;\u00e9tudes neuroendocriniennes \u00e9valuant des dosages hormonaux salivaires rares, ou de recherches sur des l\u00e9sions c\u00e9r\u00e9brales focales, recruter des cohortes massives s&rsquo;av\u00e8re fr\u00e9quemment impossible pour des raisons \u00e9thiques, logistiques ou financi\u00e8res. Dans ces contextes exp\u00e9rimentaux restreints, l&rsquo;application na\u00efve des propri\u00e9t\u00e9s de la loi normale conduit syst\u00e9matiquement \u00e0 une surconfiance d\u00e9sastreuse quant \u00e0 la pr\u00e9cision des r\u00e9sultats obtenus. L&rsquo;\u00e9chantillonnage de petite taille induit en effet une volatilit\u00e9 extr\u00eame de la variance empirique, rendant les estimations ponctuelles de dispersion hautement instables d&rsquo;un \u00e9chantillon \u00e0 l&rsquo;autre.<\/p>\n<p>La loi de Student r\u00e9pond \u00e0 cette vuln\u00e9rabilit\u00e9 \u00e9pist\u00e9mique en int\u00e9grant formellement l&rsquo;incertitude li\u00e9e \u00e0 l&rsquo;estimation de la variance directement dans la g\u00e9om\u00e9trie de sa distribution de probabilit\u00e9. Contrairement \u00e0 la courbe gaussienne, dont la d\u00e9croissance exponentielle des queues de distribution est rapide, la loi t poss\u00e8de des queues structurellement plus lourdes et \u00e9paisses, conf\u00e9rant au mod\u00e8le une propri\u00e9t\u00e9 dite leptokurtique. Cet \u00e9paississement morphologique traduit fid\u00e8lement le fait que l&rsquo;observation de valeurs substantiellement \u00e9loign\u00e9es du centre est nettement plus probable lorsqu&rsquo;on ignore la variance v\u00e9ritable de la population sous-jacente. Ainsi, la loi de Student fonctionne comme un dispositif d&rsquo;ajustement probabiliste qui p\u00e9nalise m\u00e9caniquement le manque d&rsquo;information exp\u00e9rimentale en \u00e9largissant les r\u00e9gions d&rsquo;acceptation du hasard.<\/p>\n<p>Cette particularit\u00e9 math\u00e9matique engendre des cons\u00e9quences directes sur l&rsquo;inf\u00e9rence des variables comportementales sensibles, \u00e0 l&rsquo;image des temps de r\u00e9action ou des scores de discrimination perceptive. Par exemple, lors de la mesure des temps de latence dans une t\u00e2che de type Stroop sur une douzaine de sujets, la survenue occasionnelle d&rsquo;un temps de r\u00e9ponse atypique ne doit pas \u00eatre trait\u00e9e comme un \u00e9cart extr\u00eame improbable selon une m\u00e9trique gaussienne standard. L&rsquo;usage de la distribution de Student permet d&rsquo;absorber ad\u00e9quatement cette dispersion sans rejeter ind\u00fbment l&rsquo;hypoth\u00e8se nulle d&rsquo;absence d&rsquo;effet, prot\u00e9geant ainsi la communaut\u00e9 scientifique contre l&rsquo;inflation pernicieuse des faux positifs qui a largement aliment\u00e9 la crise contemporaine de la r\u00e9plicabilit\u00e9 dans les sciences psychologiques.<\/p>\n<h3>1.3 Les degr\u00e9s de libert\u00e9 comme r\u00e9gulateur morphologique de la densit\u00e9<\/h3>\n<p>L&rsquo;architecture morphologique de la loi de Student est enti\u00e8rement dict\u00e9e par un param\u00e8tre fondamental unique : les degr\u00e9s de libert\u00e9, commun\u00e9ment d\u00e9sign\u00e9s par la lettre grecque nu ou l&rsquo;abr\u00e9viation <em>df<\/em>. D&rsquo;un point de vue conceptuel, ce param\u00e8tre refl\u00e8te le nombre d&rsquo;observations ind\u00e9pendantes disponibles dans le calcul de la statistique de test, diminu\u00e9 du nombre de param\u00e8tres interm\u00e9diaires pr\u00e9alablement estim\u00e9s \u00e0 partir de ces m\u00eames donn\u00e9es. Dans le cas canonique de l&rsquo;estimation d&rsquo;une moyenne unique sur un \u00e9chantillon de taille N, la contrainte impos\u00e9e par le calcul pr\u00e9alable de la moyenne empirique r\u00e9duit d&rsquo;une unit\u00e9 l&rsquo;ind\u00e9pendance r\u00e9siduelle des variations individuelles, fixant ainsi les degr\u00e9s de libert\u00e9 \u00e0 N &#8211; 1. Les degr\u00e9s de libert\u00e9 agissent math\u00e9matiquement comme un curseur modulant avec une pr\u00e9cision chirurgicale l&rsquo;aplatissement de la courbe et l&rsquo;\u00e9paisseur relative de ses queues de distribution.<\/p>\n<p>Sur le plan analytique formel, la fonction de densit\u00e9 de probabilit\u00e9 de la loi de Student s&rsquo;exprime \u00e0 travers l&rsquo;intervention fondamentale de la fonction Gamma d&rsquo;Euler, une g\u00e9n\u00e9ralisation continue de la fonction factorielle aux nombres r\u00e9els et complexes. L&rsquo;\u00e9quation de la densit\u00e9 univari\u00e9e centr\u00e9e standard pour une variable r\u00e9elle x s&rsquo;\u00e9tablit selon la formulation suivante :<\/p>\n<p>f(x; df) = [ Gamma((df + 1) \/ 2) \/ ( sqrt(df * pi) * Gamma(df \/ 2) ) ] * ( 1 + (x^2 \/ df) )^(-(df + 1) \/ 2)<\/p>\n<p>Ce formalisme explicite clairement comment les degr\u00e9s de libert\u00e9 conditionnent l&rsquo;ensemble de la fonction. Lorsque le param\u00e8tre df est tr\u00e8s faible, typiquement inf\u00e9rieur \u00e0 cinq, le terme polynomial au d\u00e9nominateur d\u00e9cro\u00eet beaucoup plus lentement que le terme quadratique exponentiel d&rsquo;une distribution normale, ce qui maintient une densit\u00e9 non n\u00e9gligeable \u00e0 plusieurs \u00e9carts-types du point z\u00e9ro. Cette formulation garantit que le kurtosis d&rsquo;ordre quatre de la distribution est strictement sup\u00e9rieur \u00e0 trois, confirmant l&rsquo;\u00e9tat leptokurtique de la distribution pour tout degr\u00e9 de libert\u00e9 fini sup\u00e9rieur \u00e0 quatre.<\/p>\n<p>\u00c0 mesure que la taille de l&rsquo;\u00e9chantillon cro\u00eet et que les degr\u00e9s de libert\u00e9 tendent vers l&rsquo;infini, la distribution de Student pr\u00e9sente une propri\u00e9t\u00e9 de convergence asymptotique vers la distribution de Gauss. Sur le plan math\u00e9matique, cette convergence s&rsquo;explique par le fait que la variance empirique de l&rsquo;\u00e9chantillon converge en probabilit\u00e9 vers la variance r\u00e9elle de la population, conform\u00e9ment \u00e0 la loi des grands nombres. Dans ces conditions limites, le rapport des fonctions Gamma pond\u00e9r\u00e9 par la racine carr\u00e9e de df se simplifie analytiquement, et la limite du terme polynomial correspond exactement \u00e0 l&rsquo;exponentielle de l&rsquo;oppos\u00e9 du demi-carr\u00e9 de x. En pratique computationnelle dans le logiciel R, lorsque les degr\u00e9s de libert\u00e9 d\u00e9passent une centaine d&rsquo;unit\u00e9s, les divergences entre la loi t et la loi normale deviennent minimes, bien qu&rsquo;il demeure imp\u00e9ratif de conserver la distribution exacte de Student afin d&rsquo;\u00e9viter tout biais infinit\u00e9simal dans le calcul des valeurs p et des bornes de confiance.<\/p>\n<h2>2. Architecture et vue d&rsquo;ensemble des fonctions t dans l&rsquo;environnement R<\/h2>\n<h3>2.1 La nomenclature modulaire d, p, q, r du socle stats<\/h3>\n<p>L&rsquo;une des plus remarquables r\u00e9ussites de l&rsquo;architecture logicielle du langage R r\u00e9side dans l&rsquo;harmonisation fonctionnelle et la coh\u00e9rence syntaxique de son socle statistique standard, contenu dans le paquetage de base <a rel=\"noopener noreferrer\" href=\"https:\/\/stat.ethz.ch\/R-manual\/R-devel\/library\/stats\/html\/00Index.html\">stats<\/a>. Depuis les origines du langage S con\u00e7u aux laboratoires Bell par John Chambers et ses collaborateurs, la manipulation de l&rsquo;ensemble des distributions de probabilit\u00e9 ob\u00e9it \u00e0 un sch\u00e9ma modulaire strict fond\u00e9 sur quatre pr\u00e9fixes alphab\u00e9tiques mn\u00e9motechniques universels : <strong>d<\/strong>, <strong>p<\/strong>, <strong>q<\/strong>, et <strong>r<\/strong>. Ce formalisme matriciel s&rsquo;applique avec une constance rigoureuse \u00e0 toutes les lois impl\u00e9ment\u00e9es dans le syst\u00e8me, qu&rsquo;elles soient continues ou discr\u00e8tes, formant un \u00e9cosyst\u00e8me logique o\u00f9 chaque fonction remplit un r\u00f4le probabiliste univoque et universellement pr\u00e9visible.<\/p>\n<p>Cette standardisation structurelle s&rsquo;articule autour d&rsquo;une racine lexicale d\u00e9signant la famille de lois consid\u00e9r\u00e9e, qui prend la lettre isol\u00e9e <strong>t<\/strong> pour d\u00e9signer la distribution de Student. Ainsi, le syst\u00e8me d\u00e9rive imm\u00e9diatement quatre outils fondamentaux :<\/p>\n<ul>\n<li><strong>dt<\/strong> : le calcul de la densit\u00e9 de probabilit\u00e9 (density) en un point ou le long d&rsquo;une coordonn\u00e9e d&rsquo;abscisses ;<\/li>\n<li><strong>pt<\/strong> : l&rsquo;\u00e9valuation de la fonction de r\u00e9partition empirique ou cumulative (probability distribution function) retournant la surface int\u00e9gr\u00e9e sous la courbe ;<\/li>\n<li><strong>qt<\/strong> : la d\u00e9termination des quantiles (quantile function), qui repr\u00e9sente la bijection r\u00e9ciproque de la fonction de r\u00e9partition ;<\/li>\n<li><strong>rt<\/strong> : la g\u00e9n\u00e9ration de r\u00e9alisations stochastiques pseudo-al\u00e9atoires (random generation) r\u00e9gies par cette loi sp\u00e9cifique.<\/li>\n<\/ul>\n<p>L&rsquo;interop\u00e9rabilit\u00e9 de cette architecture avec le reste du langage est totale et transparente. Les quatre fonctions sont nativement vectoris\u00e9es, ce qui signifie qu&rsquo;elles acceptent indiff\u00e9remment des arguments scalaires isol\u00e9s, des vecteurs de nombres r\u00e9els de dimensions arbitraires, ou des colonnes issues de structures de donn\u00e9es plus complexes comme les data frames ou les tibbles modernes. Cette conception vectorielle \u00e9limine le besoin d&rsquo;\u00e9crire des boucles explicites dans les scripts, favorisant une ex\u00e9cution binaire hautement optimis\u00e9e en langage C sous-jacent et garantissant une fluidit\u00e9 de traitement indispensable lors de l&rsquo;analyse massive de donn\u00e9es en sciences du comportement.<\/p>\n<h3>2.2 Distinguer densit\u00e9, fonction de r\u00e9partition, quantiles et simulation<\/h3>\n<p>La ma\u00eetrise op\u00e9rationnelle de ce quatuor exige une distinction conceptuelle absolue entre les quatre grandeurs math\u00e9matiques qu&rsquo;ils manipulent, une confusion fr\u00e9quente chez les chercheurs non math\u00e9maticiens \u00e9tant d&rsquo;assimiler la densit\u00e9 \u00e0 une probabilit\u00e9 brute. La fonction <strong>dt<\/strong> \u00e9value l&rsquo;ordonn\u00e9e d&rsquo;une courbe continue en un point pr\u00e9cis d&rsquo;abscisse t. Pour une variable continue, la probabilit\u00e9 exacte d&rsquo;observer ponctuellement une valeur sp\u00e9cifique est math\u00e9matiquement nulle. La densit\u00e9 correspond en r\u00e9alit\u00e9 \u00e0 la d\u00e9riv\u00e9e premi\u00e8re de la probabilit\u00e9 cumul\u00e9e au point consid\u00e9r\u00e9 ; elle d\u00e9crit l&rsquo;intensit\u00e9 de concentration relative de la masse probabiliste dans un intervalle infinit\u00e9simal autour de ce point. Sa valeur peut d&rsquo;ailleurs, pour certaines distributions tr\u00e8s resserr\u00e9es, d\u00e9passer l&rsquo;unit\u00e9 sans violer le second axiome de Kolmogorov.<\/p>\n<p>\u00c0 l&rsquo;inverse, la fonction de r\u00e9partition calcul\u00e9e par <strong>pt<\/strong> op\u00e8re une int\u00e9gration continue le long de l&rsquo;axe des r\u00e9els, accumulant la masse de probabilit\u00e9 depuis l&rsquo;infini n\u00e9gatif jusqu&rsquo;\u00e0 une valeur seuil donn\u00e9e. Le r\u00e9sultat renvoy\u00e9 par cette fonction repr\u00e9sente une probabilit\u00e9 stricte, born\u00e9e de fa\u00e7on rigide entre 0 et 1. C&rsquo;est cet outil probabiliste qui intervient dans la prise de d\u00e9cision inf\u00e9rentielle classique, puisqu&rsquo;il permet de d\u00e9duire la probabilit\u00e9 d&rsquo;observer sous l&rsquo;hypoth\u00e8se nulle une statistique de test au moins aussi extr\u00eame que celle relev\u00e9e au cours de l&rsquo;exp\u00e9rience, processus qui conduit \u00e0 la c\u00e9l\u00e8bre valeur p ou p-value.<\/p>\n<p>La fonction <strong>qt<\/strong> r\u00e9sout le probl\u00e8me inverse, formul\u00e9 de mani\u00e8re duale : elle prend en entr\u00e9e une probabilit\u00e9 cumul\u00e9e ou un niveau d&rsquo;erreur alpha pr\u00e9d\u00e9termin\u00e9 et renvoie la valeur critique sur l&rsquo;axe des abscisses correspondant \u00e0 cette superficie int\u00e9gr\u00e9e. Enfin, la fonction <strong>rt<\/strong> op\u00e8re dans le domaine de la simulation stochastique : elle ne calcule pas de propri\u00e9t\u00e9s analytiques th\u00e9oriques, mais \u00e9chantillonne artificiellement des valeurs num\u00e9riques en appliquant des algorithmes de transformation \u00e0 un flux de nombres pseudo-al\u00e9atoires uniformes. Cette articulation quadripartite forme un circuit ferm\u00e9 complet permettant d&rsquo;aborder toutes les facettes de l&rsquo;analyse statistique, de l&rsquo;exploration analytique \u00e0 la mod\u00e9lisation informatique la plus pouss\u00e9e.<\/p>\n<h3>2.3 Param\u00e8tres universels : degr\u00e9s de libert\u00e9 (df) et non-centralit\u00e9 (ncp)<\/h3>\n<p>L&rsquo;ensemble des fonctions du module t partagent une signature de param\u00e8tres canonique qui assure l&rsquo;homog\u00e9n\u00e9it\u00e9 de leur manipulation dans l&rsquo;environnement R. L&rsquo;argument central le plus imp\u00e9ratif est sans conteste <em>df<\/em>, qui sp\u00e9cifie les degr\u00e9s de libert\u00e9 gouvernant la dispersion de la distribution. Dans R, cet argument n&rsquo;est pas restreint aux seuls entiers naturels, bien que les plans d&rsquo;exp\u00e9rience traditionnels reposent sur des d\u00e9nombrements d&rsquo;observations entiers. Le syst\u00e8me accepte des valeurs r\u00e9elles strictement positives pour <em>df<\/em>, ce qui s&rsquo;av\u00e8re particuli\u00e8rement crucial lors de l&rsquo;application de corrections m\u00e9thodologiques avanc\u00e9es, telles que l&rsquo;ajustement de Welch-Satterthwaite pour variances h\u00e9t\u00e9rog\u00e8nes, o\u00f9 les degr\u00e9s de libert\u00e9 r\u00e9sultants sont fr\u00e9quemment fractionnaires.<\/p>\n<p>Le second param\u00e8tre d&rsquo;une importance m\u00e9thodologique capitale, bien que moins fr\u00e9quemment explicit\u00e9 par les praticiens d\u00e9butants, est le param\u00e8tre de non-centralit\u00e9, d\u00e9sign\u00e9 par l&rsquo;acronyme <em>ncp<\/em> (pour non-centrality parameter). Par d\u00e9faut, cet argument est assign\u00e9 \u00e0 la valeur scalaire 0, ce qui positionne la distribution dans sa forme centrale canonique, strictement sym\u00e9trique autour de l&rsquo;origine z\u00e9ro. C&rsquo;est cette forme centrale qui est exploit\u00e9e pour mod\u00e9liser le comportement de la statistique sous l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9galit\u00e9 des moyennes ou d&rsquo;absence d&rsquo;effet exp\u00e9rimental. En revanche, d\u00e8s lors que l&rsquo;on s&rsquo;int\u00e9resse \u00e0 l&rsquo;hypoth\u00e8se alternative d&rsquo;existence d&rsquo;un effet non nul dans la population, le param\u00e8tre <em>ncp<\/em> prend une valeur non nulle proportionnelle \u00e0 la taille d&rsquo;effet th\u00e9orique attendue et \u00e0 la taille de l&rsquo;\u00e9chantillon, brisant la sym\u00e9trie de la courbe et d\u00e9calant son apex.<\/p>\n<p>Il convient de souligner que R g\u00e8re ces arguments avec une rigueur computationnelle intransigeante. Si l&rsquo;utilisateur omet de renseigner le param\u00e8tre <em>df<\/em>, l&rsquo;interpr\u00e9teur interrompt imm\u00e9diatement l&rsquo;ex\u00e9cution en g\u00e9n\u00e9rant une erreur formelle, ce param\u00e8tre ne poss\u00e9dant pas de valeur par d\u00e9faut arbitraire. En outre, la vectorisation de R s&rsquo;\u00e9tend \u00e0 ces param\u00e8tres : il est tout \u00e0 fait envisageable de fournir un vecteur de degr\u00e9s de libert\u00e9 ou un vecteur de param\u00e8tres de non-centralit\u00e9 de m\u00eame longueur que le vecteur de donn\u00e9es \u00e9valu\u00e9. Enfin, toute pr\u00e9sence de donn\u00e9es manquantes (not\u00e9es <em>NA<\/em> ou <em>NaN<\/em>) au sein des arguments en entr\u00e9e se propage de fa\u00e7on d\u00e9terministe dans le r\u00e9sultat retourn\u00e9, garantissant une int\u00e9grit\u00e9 computationnelle totale \u00e9vitant les inf\u00e9rences fallacieuses bas\u00e9es sur des donn\u00e9es incompl\u00e8tes.<\/p>\n<h2>3. La fonction dt() : \u00c9valuation de la fonction de densit\u00e9 de probabilit\u00e9<\/h2>\n<h3>3.1 Formulation formelle et syntaxe canonique de dt()<\/h3>\n<p>La fonction <strong>dt()<\/strong> constitue le noyau computationnel d\u00e9di\u00e9 au calcul de la hauteur de l&rsquo;ordonn\u00e9e de la courbe de distribution pour un point donn\u00e9 de l&rsquo;axe r\u00e9el. La syntaxe canonique de cette fonction, telle que document\u00e9e dans les manuels officiels de R, se d\u00e9ploie selon la signature technique suivante :<\/p>\n<p>dt(x, df, ncp, log = FALSE)<\/p>\n<p>Le premier argument, <em>x<\/em>, repr\u00e9sente la valeur de la statistique t ou le vecteur d&rsquo;abscisses num\u00e9riques auquel la fonction de densit\u00e9 doit \u00eatre \u00e9valu\u00e9e. L&rsquo;argument <em>df<\/em> exige un scalaire ou un vecteur num\u00e9rique positif sp\u00e9cifiant la morphologie des degr\u00e9s de libert\u00e9. L&rsquo;argument optionnel <em>ncp<\/em>, s&rsquo;il est sp\u00e9cifi\u00e9, permet de basculer instantan\u00e9ment de l&rsquo;\u00e9valuation de la loi de Student centrale vers sa variante non-centrale asym\u00e9trique, bien que son omission implicite active automatiquement l&rsquo;algorithme standard centr\u00e9 autour de z\u00e9ro.<\/p>\n<p>L&rsquo;argument logique <em>log<\/em>, configur\u00e9 par d\u00e9faut sur <em>FALSE<\/em>, rev\u00eat une importance capitale pour la stabilit\u00e9 des calculs num\u00e9riques avanc\u00e9s. Lorsqu&rsquo;il est bascul\u00e9 sur <em>TRUE<\/em>, la fonction ne renvoie pas la densit\u00e9 absolue f(x), mais son logarithme n\u00e9p\u00e9rien naturel, ln(f(x)). Dans le cadre d&rsquo;estimations par le maximum de vraisemblance en psychom\u00e9trie computationnelle ou lors de l&rsquo;ajustement de mod\u00e8les hi\u00e9rarchiques complexes, les chercheurs doivent fr\u00e9quemment multiplier des densit\u00e9s conjointes associ\u00e9es \u00e0 des centaines de points d&rsquo;observation. La multiplication directe de densit\u00e9s infinit\u00e9simales conduit in\u00e9vitablement \u00e0 un ph\u00e9nom\u00e8ne mat\u00e9riel d&rsquo;annulation sous le seuil de pr\u00e9cision machine (connu sous le terme anglais d&rsquo;<em>underflow<\/em>), o\u00f9 des valeurs extr\u00eamement proches de z\u00e9ro sont tronqu\u00e9es arbitrairement \u00e0 z\u00e9ro absolu par le processeur. L&rsquo;\u00e9valuation directe dans l&rsquo;espace logarithmique au moyen de <em>log = TRUE<\/em> substitue la somme arithm\u00e9tique aux produits factoriels instables, pr\u00e9servant ainsi l&rsquo;exactitude num\u00e9rique sur plusieurs ordres de grandeur.<\/p>\n<h3>3.2 D\u00e9mystifier la valeur retourn\u00e9e : densit\u00e9 n&rsquo;est pas probabilit\u00e9<\/h3>\n<p>L&rsquo;une des erreurs conceptuelles les plus persistantes chez les \u00e9tudiants et les praticiens de la recherche empirique consiste \u00e0 interpr\u00e9ter le nombre scalaire renvoy\u00e9 par <strong>dt(x, df)<\/strong> comme la probabilit\u00e9 formelle que la variable t prenne pr\u00e9cis\u00e9ment la valeur x. Il est fondamental de r\u00e9affirmer qu&rsquo;au sein de tout mod\u00e8le probabiliste \u00e0 variable continue r\u00e9gi par les axiomes de Lebesgue et de Kolmogorov, la probabilit\u00e9 ponctuelle P(T = x) est strictement \u00e9gale \u00e0 z\u00e9ro pour n&rsquo;importe quelle valeur r\u00e9elle isol\u00e9e. La surface g\u00e9om\u00e9trique correspondant \u00e0 une simple ligne verticale d&rsquo;\u00e9paisseur nulle au-dessus d&rsquo;un point est n\u00e9cessairement vide d&rsquo;aire.<\/p>\n<p>La valeur retourn\u00e9e par la fonction <strong>dt()<\/strong> doit \u00eatre rigoureusement comprise comme un taux de probabilit\u00e9 par unit\u00e9 de mesure, mesurant l&rsquo;intensit\u00e9 de la probabilit\u00e9 marginale locale. Si l&rsquo;on consid\u00e8re un intervalle infinit\u00e9simalement \u00e9troit compris entre le point x et x + dx, la probabilit\u00e9 math\u00e9matique que la statistique de Student se r\u00e9alise \u00e0 l&rsquo;int\u00e9rieur de cet intervalle est donn\u00e9e par l&rsquo;approximation du premier ordre f(x) * dx. Par cons\u00e9quent, la densit\u00e9 f(x) correspond \u00e0 la limite du rapport entre la probabilit\u00e9 cumul\u00e9e sur une tranche arbitrairement mince et la largeur de cette m\u00eame tranche lorsque celle-ci tend vers z\u00e9ro :<\/p>\n<p>f(x) = lim_{dx -&gt; 0} [ P(x &lt;= T &lt;= x + dx) \/ dx ]<\/p>\n<p>Cette distinction conceptuelle n&rsquo;est pas une simple coquetterie de formalisme math\u00e9matique ; elle prot\u00e8ge les chercheurs contre des aberrations d\u00e9ductives majeures lors de l&rsquo;interpr\u00e9tation de mod\u00e8les structurels. Dans certaines distributions continues tr\u00e8s comprim\u00e9es, la densit\u00e9 peut d\u00e9passer largement des valeurs de 1.0 ou 10.0 sans enfreindre aucun principe probabiliste, car seule l&rsquo;int\u00e9grale de cette densit\u00e9 sur l&rsquo;ensemble de l&rsquo;espace des r\u00e9els, allant de moins l&rsquo;infini \u00e0 plus l&rsquo;infini, est tenue de valoir exactement 1.0. Lors de l&rsquo;usage de <strong>dt()<\/strong>, le r\u00e9sultat indique donc uniquement la plausibilit\u00e9 relative d&rsquo;un \u00e9tat statistique par rapport \u00e0 un autre sur l&rsquo;axe continu de dispersion des donn\u00e9es.<\/p>\n<h3>3.3 Impl\u00e9mentation d&rsquo;exemples num\u00e9riques pas \u00e0 pas<\/h3>\n<p>Pour appr\u00e9hender concr\u00e8tement la dynamique calculatoire de <strong>dt()<\/strong>, observons son comportement au sommet g\u00e9om\u00e9trique de la distribution, situ\u00e9 pour une loi centrale au point d&rsquo;abscisse t = 0. En \u00e9valuant l&rsquo;expression <em>dt(x = 0, df = 5)<\/em>, l&rsquo;interpr\u00e9teur R mobilise la formule analytique de la fonction Gamma et d\u00e9livre une ordonn\u00e9e de sommet approximativement \u00e9gale \u00e0 0.3796. Si l&rsquo;on augmente consid\u00e9rablement le param\u00e8tre de dispersion en \u00e9valuant <em>dt(x = 0, df = 30)<\/em>, l&rsquo;apex s&rsquo;\u00e9l\u00e8ve pour atteindre 0.3953. Enfin, lorsque les degr\u00e9s de libert\u00e9 tendent virtuellement vers l&rsquo;infini, par exemple avec <em>dt(x = 0, df = 1000)<\/em>, la hauteur maximale converge asymptotiquement vers la valeur universelle 1 \/ sqrt(2 * pi), qui vaut approximativement 0.3989, soit l&rsquo;ordonn\u00e9e centrale exacte de la courbe normale centr\u00e9e r\u00e9duite.<\/p>\n<p>Examinons \u00e0 pr\u00e9sent ce qui survient lorsque l&rsquo;on interroge la fonction sur des points excentr\u00e9s, repr\u00e9sentatifs des queues de distribution, par exemple \u00e0 une distance de trois \u00e9carts-types du centre (t = 3). L&rsquo;appel de l&rsquo;instruction <em>dt(x = 3, df = 4)<\/em> renvoie une densit\u00e9 r\u00e9siduelle de 0.0207. En revanche, l&rsquo;\u00e9valuation de cette m\u00eame coordonn\u00e9e pour une loi reposant sur cent degr\u00e9s de libert\u00e9 via <em>dt(x = 3, df = 100)<\/em> produit un r\u00e9sultat de 0.0051, soit une densit\u00e9 quatre fois plus faible. Cette divergence num\u00e9rique illustre pr\u00e9cis\u00e9ment la leptokurticit\u00e9 : \u00e0 des distances critiques de l&rsquo;esp\u00e9rance, la loi de Student attribue un poids substantiellement sup\u00e9rieur aux valeurs aberrantes par rapport aux mod\u00e8les fond\u00e9s sur de grands \u00e9chantillons.<\/p>\n<p>Enfin, la puissance de la vectorisation native en R permet d&rsquo;\u00e9valuer simultan\u00e9ment une pluralit\u00e9 de sc\u00e9narios exp\u00e9rimentaux en une unique ligne de code synth\u00e9tique. En injectant un vecteur d&rsquo;abscisses r\u00e9gulier dans la fonction, tel que <em>dt(x = c(-2, -1, 0, 1, 2), df = 10)<\/em>, R retourne instantan\u00e9ment un vecteur num\u00e9rique contenant les cinq hauteurs de densit\u00e9 sym\u00e9triques correspondantes : 0.0611, 0.2360, 0.3891, 0.2360 et 0.0611. Cette propri\u00e9t\u00e9 de sym\u00e9trie parfaite autour de l&rsquo;axe z\u00e9ro se v\u00e9rifie rigoureusement pour toutes les \u00e9valuations \u00e9quidistantes sous hypoth\u00e8se centrale, offrant une base robuste pour la construction de repr\u00e9sentations visuelles des profils statistiques.<\/p>\n<h2>4. Mod\u00e9lisation visuelle des distributions th\u00e9oriques avec dt() et ggplot2<\/h2>\n<h3>4.1 G\u00e9n\u00e9ration de s\u00e9quences coordonn\u00e9es et trac\u00e9s graphiques de base<\/h3>\n<p>L&rsquo;exploration visuelle de la morphologie des fonctions de densit\u00e9 constitue une \u00e9tape p\u00e9dagogique et m\u00e9thodologique fondamentale pour appr\u00e9cier concr\u00e8tement la variation de courbure en fonction des param\u00e8tres d&rsquo;inf\u00e9rence. Le m\u00e9canisme standard pour g\u00e9n\u00e9rer un canevas graphique dans le syst\u00e8me R de base repose sur l&rsquo;articulation de la fonction <strong>seq()<\/strong> et des outils de tra\u00e7age de bas niveau issus de la biblioth\u00e8que graphique native. En premi\u00e8re instance, le chercheur cr\u00e9e un maillage d&rsquo;abscisses continues finement discr\u00e9tis\u00e9, par exemple en \u00e9crivant <em>valeurs_t &lt;- seq(from = -4, to = 4, length.out = 500)<\/em>. Ce vecteur ordonn\u00e9 g\u00e9n\u00e8re une s\u00e9rie dense de points de mesure garantissant une fluidit\u00e9 parfaite du rendu visuel des courbures sans brisure polygonale.<\/p>\n<p>Une fois cette coordonn\u00e9e horizontale \u00e9tablie, il suffit d&rsquo;appliquer la fonction <strong>dt()<\/strong> sur ce vecteur pour en d\u00e9duire les hauteurs verticales correspondantes. La commande <em>densite_t &lt;- dt(x = valeurs_t, df = 5)<\/em> calcule de fa\u00e7on vectoris\u00e9e les 500 ordonn\u00e9es associ\u00e9es. L&rsquo;appel primitif \u00e0 la fonction <em>plot(valeurs_t, densite_t, type = \u00ab\u00a0l\u00a0\u00bb, lwd = 2, col = \u00ab\u00a0blue\u00a0\u00bb)<\/em> trace alors la courbe de Student correspondante. Pour enrichir la visualisation et rendre sensible la notion de convergence morphologique, il est d&rsquo;usage d&rsquo;adjoindre des lignes secondaires au moyen de la primitive <em>lines()<\/em>, en superposant des distributions reposant sur des degr\u00e9s de libert\u00e9 croissants, par exemple df = 2, df = 10 et df = 50, en utilisant des codes couleur distinctifs et un jeu d&rsquo;\u00e9paisseurs de trait hi\u00e9rarchis\u00e9.<\/p>\n<p>Cette approche basique, bien que tr\u00e8s \u00e9conome en ressources de calcul et universellement disponible sans n\u00e9cessiter l&rsquo;installation de biblioth\u00e8ques tierces, pr\u00e9sente n\u00e9anmoins des limitations s\u00e9v\u00e8res en mati\u00e8re de mise en page pour la publication scientifique. L&rsquo;alignement manuel des \u00e9chelles, la configuration minutieuse des l\u00e9gendes et l&rsquo;ajout d&rsquo;ombrages d\u00e9limitant des surfaces critiques exigent un codage verbeux et r\u00e9p\u00e9titif. C&rsquo;est la raison pour laquelle les chercheurs en psychologie privil\u00e9gient massivement l&rsquo;utilisation de l&rsquo;\u00e9cosyst\u00e8me graphique d\u00e9claratif fond\u00e9 sur la grammaire des graphiques.<\/p>\n<h3>4.2 Conception de figures d&rsquo;\u00e9dition avec ggplot2 et stat_function<\/h3>\n<p>L&rsquo;utilisation de la biblioth\u00e8que <a rel=\"noopener noreferrer\" href=\"https:\/\/ggplot2.tidyverse.org\/\">ggplot2<\/a> transcende la production graphique en introduisant une abstraction \u00e9l\u00e9gante o\u00f9 l&rsquo;\u00e9valuation des distributions th\u00e9oriques s&rsquo;affranchit totalement de la cr\u00e9ation manuelle pr\u00e9alable de jeux de donn\u00e9es tabulaires d&rsquo;ordonn\u00e9es. Gr\u00e2ce \u00e0 la couche g\u00e9om\u00e9trique sp\u00e9cialis\u00e9e <strong>stat_function()<\/strong>, il suffit de d\u00e9finir un cadre d&rsquo;observation sur un intervalle d&rsquo;abscisses restreint, puis de d\u00e9l\u00e9guer \u00e0 ggplot2 l&rsquo;\u00e9chantillonnage continu et le trac\u00e9 direct de la fonction analytique de Student.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/tDist1.jpg\" alt=\"Student t Distribution plot in R\" \/><figcaption>Student t Distribution plot in R<\/figcaption><\/figure>\n<p>Pour illustrer la mise en forme rigoureuse d&rsquo;une figure destin\u00e9e \u00e0 l&rsquo;\u00e9dition universitaire selon les strictes directives de l&rsquo;American Psychological Association (APA), le chercheur initialise une toile vide avec <em>ggplot(data = data.frame(x = c(-4, 4)), aes(x = x))<\/em>. Ensuite, l&rsquo;adjonction de couches d\u00e9claratives successives structure le rendu final :<\/p>\n<p>stat_function(fun = dt, args = list(df = 8), color = \u00ab\u00a0steelblue\u00a0\u00bb, size = 1.2)<\/p>\n<p>Cette syntaxe compacte sp\u00e9cifie directement la fonction de densit\u00e9 \u00e0 \u00e9valuer et lui transmet dynamiquement sa liste d&rsquo;arguments structurels. La personnalisation de la figure requiert l&rsquo;\u00e9radication des grilles d&rsquo;arri\u00e8re-plan superflues et des aplats gris\u00e2tres via l&rsquo;application du th\u00e8me formel <em>theme_classic()<\/em>, la configuration soign\u00e9e des polices typographiques institutionnelles, et l&rsquo;\u00e9tiquetage pr\u00e9cis des axes pr\u00e9cisant la m\u00e9trique de la statistique de test sur l&rsquo;axe des x et la densit\u00e9 continue sur l&rsquo;axe des y.<\/p>\n<p>De surcro\u00eet, la mise en exergue des zones de rejet de l&rsquo;hypoth\u00e8se nulle \u00e0 un seuil conventionnel donn\u00e9 (par exemple alpha = 0.05 bilat\u00e9ral) s&rsquo;op\u00e8re par l&rsquo;adjonction de couches g\u00e9om\u00e9triques de remplissage de polygones, au moyen de l&rsquo;argument <em>geom = \u00ab\u00a0ribbon\u00a0\u00bb<\/em> encapsul\u00e9 dans <em>stat_function<\/em>. En bornant les fonctions d&rsquo;ombrage entre les quantiles critiques n\u00e9gatifs et positifs et les extr\u00e9mit\u00e9s de l&rsquo;axe, la figure met graphiquement en relief l&rsquo;aire r\u00e9siduelle de probabilit\u00e9 sous les queues de la distribution. Ce type de repr\u00e9sentation explicative, hautement didactique, garantit une transmission transparente et imm\u00e9diate des r\u00e9sultats inf\u00e9rentiels lors de la r\u00e9daction de manuscrits de recherche ou de pr\u00e9sentations en conf\u00e9rences internationales.<\/p>\n<h3>4.3 Visualisation comparative : Student face \u00e0 la loi normale standard<\/h3>\n<p>L&rsquo;un des exercices visuels les plus instructifs dans l&rsquo;enseignement et la pratique de la statistique comportementale r\u00e9side dans la superposition directe, sur un m\u00eame rep\u00e8re cart\u00e9sien, de la fonction de densit\u00e9 de Student avec celle de la loi normale centr\u00e9e r\u00e9duite. En confrontant <strong>dt(x, df)<\/strong> et <strong>dnorm(x, mean = 0, sd = 1)<\/strong> sur un intervalle d&rsquo;abscisses s&rsquo;\u00e9tendant typiquement de -4 \u00e0 +4, les propri\u00e9t\u00e9s math\u00e9matiques fondamentales \u00e9nonc\u00e9es pr\u00e9c\u00e9demment deviennent intuitivement \u00e9videntes pour l&rsquo;observateur.<\/p>\n<p>Sur une telle composition graphique comparative, la loi normale s&rsquo;\u00e9l\u00e8ve au sommet comme une courbe plus \u00e9troite et plus \u00e9lanc\u00e9e, culminant pr\u00e9cis\u00e9ment \u00e0 0.3989, tandis que la courbe de Student \u00e0 faibles degr\u00e9s de libert\u00e9 (par exemple df = 3) affiche un sommet abaiss\u00e9 et \u00e9mouss\u00e9. Plus significatif encore est l&rsquo;examen attentif des r\u00e9gions lat\u00e9rales situ\u00e9es au-del\u00e0 de deux unit\u00e9s de distance par rapport \u00e0 l&rsquo;origine : \u00e0 cet endroit pr\u00e9cis, les courbes se croisent, et la loi de Student passe r\u00e9solument au-dessus de la courbe gaussienne. Cet exc\u00e9dent de masse de probabilit\u00e9 visible dans les queues de distribution d\u00e9montre pourquoi une statistique t de valeur 2.2 peut se situer largement au-del\u00e0 du seuil critique de significativit\u00e9 sous l&rsquo;assomption normale, mais demeurer parfaitement compatible avec le hasard sous la loi de Student associ\u00e9e \u00e0 un petit \u00e9chantillon.<\/p>\n<p>Pour accentuer la valeur p\u00e9dagogique d&rsquo;une telle illustration comparative, il est particuli\u00e8rement pertinent d&rsquo;appliquer des colorations diff\u00e9rentielles semi-transparentes sous les queues de distribution respectives. La zone de discr\u00e9pance entre l&rsquo;aire normale et l&rsquo;aire de Student mat\u00e9rialise directement ce que l&rsquo;on qualifie d&rsquo;erreur d&rsquo;ajustement \u00e9pist\u00e9mique. Ce contraste visuel saisissant met en garde les investigateurs contre la tentation syst\u00e9matique d&rsquo;invoquer l&rsquo;approximation normale lorsque les effectifs de leurs groupes exp\u00e9rimentaux demeurent modestes, formalisant graphiquement la n\u00e9cessit\u00e9 math\u00e9matique d&rsquo;incorporer l&rsquo;incertitude li\u00e9e \u00e0 l&rsquo;estimation de la variance empirique.<\/p>\n<h2>5. La fonction pt() : Probabilit\u00e9s cumul\u00e9es et calcul de la valeur p<\/h2>\n<h3>5.1 M\u00e9canismes de la fonction de r\u00e9partition empirique<\/h3>\n<p>Alors que la fonction <strong>dt()<\/strong> op\u00e8re une mesure ponctuelle sur la courbe de densit\u00e9, la fonction <strong>pt()<\/strong> impl\u00e9mente le calcul infinit\u00e9simal de l&rsquo;aire sous la courbe, incarnant la fonction de r\u00e9partition cumulative (d\u00e9sign\u00e9e formellement par CDF, pour Cumulative Distribution Function) de la loi de Student. Sur le plan de l&rsquo;analyse math\u00e9matique, la fonction de r\u00e9partition en une coordonn\u00e9e quantile q est d\u00e9finie comme l&rsquo;int\u00e9grale d\u00e9finie de la fonction de densit\u00e9 de probabilit\u00e9 depuis la limite asymptotique n\u00e9gative jusqu&rsquo;au point d&rsquo;\u00e9valuation consid\u00e9r\u00e9 :<\/p>\n<p>F(q; df) = P(T &lt;= q) = int_{-infty}^{q} f(t; df) dt<\/p>\n<p>La signature op\u00e9rationnelle de cette fonction fondamentale dans R s&rsquo;articule comme suit :<\/p>\n<p>pt(q, df, ncp, lower.tail = TRUE, log.p = FALSE)<\/p>\n<p>L&rsquo;argument premier, <em>q<\/em>, d\u00e9signe le vecteur de quantiles ou les valeurs observ\u00e9es de la statistique empirique issue de l&rsquo;analyse des donn\u00e9es comportementales. Le param\u00e8tre <em>df<\/em> fixe de nouveau la courbure via les degr\u00e9s de libert\u00e9 th\u00e9oriques. Par d\u00e9faut, la fonction assigne l&rsquo;argument logique <em>lower.tail = TRUE<\/em>, ce qui signifie que l&rsquo;algorithme proc\u00e8de \u00e0 l&rsquo;int\u00e9gration de la masse de probabilit\u00e9 depuis l&rsquo;extr\u00e9mit\u00e9 gauche (-Inf) jusqu&rsquo;au quantile q sp\u00e9cifi\u00e9. La valeur num\u00e9rique d\u00e9livr\u00e9e en retour repr\u00e9sente fid\u00e8lement la probabilit\u00e9 unilat\u00e9rale gauche P(T &lt;= q), une quantit\u00e9 math\u00e9matique strictement comprise dans l&rsquo;intervalle ouvert ]0, 1[.<\/p>\n<p>L&rsquo;int\u00e9gration analytique de la loi de Student fait intervenir de mani\u00e8re computationnelle des fonctions hyperg\u00e9om\u00e9triques ou la fonction b\u00eata incompl\u00e8te r\u00e9gularis\u00e9e. Sous le capot de l&rsquo;interpr\u00e9teur R, ces calculs ne sont pas ex\u00e9cut\u00e9s au moyen d&rsquo;approximations de Riemann grossi\u00e8res, mais s&rsquo;appuient sur des routines compil\u00e9es de tr\u00e8s haute pr\u00e9cision issues de la c\u00e9l\u00e8bre biblioth\u00e8que Fortran\/C <em>Applied Statistics<\/em>. Cette conception garantit une exactitude arithm\u00e9tique quasi absolue jusqu&rsquo;\u00e0 la limite de pr\u00e9cision des nombres \u00e0 virgule flottante en double pr\u00e9cision standard (norme IEEE 754), assurant aux chercheurs une fiabilit\u00e9 sans faille lors de l&rsquo;extraction des probabilit\u00e9s associ\u00e9es \u00e0 leurs exp\u00e9rimentations.<\/p>\n<h3>5.2 Ma\u00eetrise du param\u00e8tre lower.tail pour les tests unilat\u00e9raux et bilat\u00e9raux<\/h3>\n<p>Une compr\u00e9hension subtile de l&rsquo;argument <em>lower.tail<\/em> est indispensable pour pr\u00e9venir des erreurs d&rsquo;arrondi dramatiques et des contresens m\u00e9thodologiques dans l&rsquo;\u00e9valuation des tests d&rsquo;hypoth\u00e8se. Par d\u00e9finition sym\u00e9trique, si l&rsquo;on cherche \u00e0 \u00e9valuer la probabilit\u00e9 critique dans la queue sup\u00e9rieure droite d&rsquo;une distribution \u2014 situation classique lors d&rsquo;un test directionnel affirmant qu&rsquo;une intervention psychologique accro\u00eet les performances cognitives \u2014, la probabilit\u00e9 d&rsquo;int\u00e9r\u00eat correspond \u00e0 P(T &gt; q). L&rsquo;utilisateur d\u00e9butant pourrait \u00eatre tent\u00e9 d&rsquo;\u00e9crire l&rsquo;expression na\u00efve <em>1 &#8211; pt(q, df)<\/em>. Bien que th\u00e9oriquement exacte sur le plan alg\u00e9brique, cette formulation s&rsquo;av\u00e8re hautement pr\u00e9judiciable sur le plan informatique d\u00e8s lors que le quantile q se situe \u00e0 des distances extr\u00eames dans la queue de distribution.<\/p>\n<p>Dans l&rsquo;architecture informatique moderne, la pr\u00e9cision machine limite le stockage des d\u00e9cimales \u00e0 environ seize chiffres significatifs. Si la valeur de <em>pt(q, df)<\/em> est par exemple \u00e9gale \u00e0 0.9999999999999999, l&rsquo;op\u00e9ration arithm\u00e9tique de soustraction soumise au seuil de pr\u00e9cision machine va tronquer la diff\u00e9rence \u00e0 0 absolu, annihilant toute estimation de la p-valeur. \u00c0 l&rsquo;inverse, l&rsquo;assignation explicite de l&rsquo;argument <em>lower.tail = FALSE<\/em> ordonne \u00e0 l&rsquo;algorithme sous-jacent de calculer directement l&rsquo;int\u00e9grale depuis le point q jusqu&rsquo;\u00e0 l&rsquo;infini positif :<\/p>\n<p>pt(q, df, lower.tail = FALSE)<\/p>\n<p>Cette approche permet d&rsquo;obtenir des grandeurs infinit\u00e9simales pr\u00e9cises, telles que 1.42e-28, sans souffrir du ph\u00e9nom\u00e8ne d&rsquo;annulation soustractive. En outre, si la valeur est encore plus extr\u00eame, l&rsquo;adjonction conjointe de <em>log.p = TRUE<\/em> retourne le logarithme naturel exact de la probabilit\u00e9 de queue sans rupture de pile machine.<\/p>\n<p>Pour le calcul classique de la p-valeur dans le cadre d&rsquo;un test non-directionnel (bilat\u00e9ral), le chercheur postule que la d\u00e9viation empirique peut se manifester indiff\u00e9remment dans l&rsquo;une ou l&rsquo;autre des directions de l&rsquo;espace d&rsquo;\u00e9chantillonnage. En vertu de la sym\u00e9trie absolue de la distribution de Student centrale autour de z\u00e9ro, la p-valeur bilat\u00e9rale s&rsquo;obtient en isolant la probabilit\u00e9 r\u00e9siduelle de la queue la plus extr\u00eame et en la multipliant par deux. La formulation canonique la plus robuste sur le plan algorithmique consiste ainsi \u00e0 extraire la valeur absolue de la statistique observ\u00e9e, \u00e0 calculer son int\u00e9grale sup\u00e9rieure droite et \u00e0 doubler le r\u00e9sultat obtenu :<\/p>\n<p>p_valeur_bilaterale &lt;- 2 * pt(abs(t_observe), df = degres_liberte, lower.tail = FALSE)<\/p>\n<p>Cette formulation \u00e9l\u00e9gante traite de fa\u00e7on rigoureusement identique les statistiques empiriques positives et n\u00e9gatives, \u00e9vitant tout risque d&rsquo;inversion accidentelle de la zone d&rsquo;int\u00e9gration.<\/p>\n<h3>5.3 Conversion manuelle d&rsquo;une statistique t observ\u00e9e en p-valeur exacte<\/h3>\n<p>Pour illustrer la mise en application directe de ces principes, consid\u00e9rons un protocole en psychologie cognitive \u00e9valuant l&rsquo;effet d&rsquo;une privation de sommeil sur une t\u00e2che d&rsquo;attention soutenue. Supposons qu&rsquo;un \u00e9chantillon unique de 15 participants produise une moyenne de scores d&rsquo;erreurs significativement d\u00e9cal\u00e9e par rapport au standard normatif, aboutissant apr\u00e8s calcul alg\u00e9brique \u00e0 une statistique de test empirique <em>t_obs = 2.45<\/em> avec des degr\u00e9s de libert\u00e9 <em>df = 15 &#8211; 1 = 14<\/em>.<\/p>\n<p>Afin de d\u00e9terminer manuellement la p-valeur exacte associ\u00e9e \u00e0 ce r\u00e9sultat sans recourir aveugl\u00e9ment \u00e0 des proc\u00e9dures d&rsquo;analyse automatis\u00e9es, le chercheur ex\u00e9cute dans la console R l&rsquo;instruction unilat\u00e9rale droite <em>pt(2.45, df = 14, lower.tail = FALSE)<\/em>. Le syst\u00e8me calcule instantan\u00e9ment une probabilit\u00e9 de queue unilat\u00e9rale valant approximativement 0.0140. Dans le cadre d&rsquo;une hypoth\u00e8se directionnelle a priori sp\u00e9cifiant que le manque de sommeil d\u00e9t\u00e9riore imp\u00e9rativement la performance, cette valeur de 0.0140 constitue la p-valeur d\u00e9finitive. Si la recherche s&rsquo;inscrit au contraire dans un paradigme bilat\u00e9ral exploratoire, l&rsquo;instruction <em>2 * pt(2.45, df = 14, lower.tail = FALSE)<\/em> d\u00e9livre une p-valeur bilat\u00e9rale exacte de 0.0280.<\/p>\n<p>Pour valider la conformit\u00e9 absolue de ce calcul manuel, il est ais\u00e9 de confronter ce r\u00e9sultat avec la sortie g\u00e9n\u00e9r\u00e9e par la fonction int\u00e9gr\u00e9e <strong>t.test()<\/strong> de R appliqu\u00e9e sur les donn\u00e9es brutes sous-jacentes. L&rsquo;extraction de l&rsquo;\u00e9l\u00e9ment <em>modele$p.value<\/em> retournera rigoureusement cette m\u00eame valeur num\u00e9rique au bit pr\u00e8s. Enfin, consid\u00e9rons le cas d&rsquo;une statistique n\u00e9gative prononc\u00e9e, par exemple <em>t_obs = -3.85<\/em> avec <em>df = 8<\/em>. L&rsquo;ex\u00e9cution de <em>2 * pt(-abs(-3.85), df = 8, lower.tail = TRUE)<\/em> ou de mani\u00e8re \u00e9quivalente <em>2 * pt(abs(-3.85), df = 8, lower.tail = FALSE)<\/em> renvoie une probabilit\u00e9 bilat\u00e9rale hautement significative de 0.0049. Cette d\u00e9composition analytique garantit une tra\u00e7abilit\u00e9 totale du processus de d\u00e9cision statistique, d\u00e9mystifiant enti\u00e8rement les bo\u00eetes noires logicielles.<\/p>\n<h2>6. Applications de pt() aux protocoles exp\u00e9rimentaux en psychologie<\/h2>\n<h3>6.1 Comparaison de deux groupes ind\u00e9pendants avec variances in\u00e9gales (Welch)<\/h3>\n<p>L&rsquo;un des d\u00e9fis m\u00e9thodologiques les plus omnipr\u00e9sents dans la recherche exp\u00e9rimentale concerne la comparaison de deux groupes ind\u00e9pendants lorsque le postulat d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances (homosc\u00e9dasticit\u00e9) est transgress\u00e9. Historiquement, l&rsquo;application du test t de Student standard \u00e0 deux \u00e9chantillons pools la variance des deux groupes, ce qui produit un biais dramatique sur le taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce lorsque les effectifs des groupes sont asym\u00e9triques. Pour pallier cette faiblesse structurelle, la correction de Welch-Satterthwaite modifie le calcul du ratio t et ajuste les degr\u00e9s de libert\u00e9 en une quantit\u00e9 r\u00e9elle souvent non enti\u00e8re.<\/p>\n<p>Supposons une exp\u00e9rience clinique \u00e9valuant l&rsquo;efficacit\u00e9 d&rsquo;une psychoth\u00e9rapie innovante sur la r\u00e9duction de l&rsquo;anxi\u00e9t\u00e9 : un groupe exp\u00e9rimental de n1 = 12 patients est confront\u00e9 \u00e0 un groupe t\u00e9moin sur liste d&rsquo;attente de n2 = 22 participants. Les variances empiriques mesur\u00e9es sont respectivement s1^2 = 18.4 et s2^2 = 5.2. La statistique t ajust\u00e9e de Welch est calcul\u00e9e manuellement \u00e0 t = 2.68. La formule classique de Welch-Satterthwaite pour les degr\u00e9s de libert\u00e9 effectifs se formule selon le quotient :<\/p>\n<p>df_welch = [ (s1^2 \/ n1) + (s2^2 \/ n2) ]^2 \/ [ ( (s1^2 \/ n1)^2 \/ (n1 &#8211; 1) ) + ( (s2^2 \/ n2)^2 \/ (n2 &#8211; 1) ) ]<\/p>\n<p>Apr\u00e8s \u00e9valuation arithm\u00e9tique dans R, cette expression produit une valeur fractionnaire continue, par exemple <em>df_welch = 15.342<\/em>. C&rsquo;est ici que l&rsquo;interop\u00e9rabilit\u00e9 math\u00e9matique de <strong>pt()<\/strong> r\u00e9v\u00e8le toute sa flexibilit\u00e9 conceptuelle : la fonction accepte sans la moindre h\u00e9sitation ce degr\u00e9 de libert\u00e9 non entier :<\/p>\n<p>p_welch &lt;- 2 * pt(abs(2.68), df = 15.342, lower.tail = FALSE)<\/p>\n<p>Le calcul d\u00e9livre une p-valeur exacte de 0.0169. Si l&rsquo;analyste avait arrondi na\u00efvement les degr\u00e9s de libert\u00e9 \u00e0 l&rsquo;entier le plus proche ou utilis\u00e9 les degr\u00e9s de libert\u00e9 nominaux non corrig\u00e9s (12 + 22 &#8211; 2 = 32), la p-valeur aurait \u00e9t\u00e9 sous-estim\u00e9e, illustrant l&rsquo;utilit\u00e9 indispensable de <strong>pt()<\/strong> pour appliquer les corrections contemporaines les plus exigeantes.<\/p>\n<h3>6.2 Analyse de mesures r\u00e9p\u00e9t\u00e9es et devis pr\u00e9-test \/ post-test<\/h3>\n<p>Les devis longitudinaux ou \u00e0 mesures r\u00e9p\u00e9t\u00e9es intrajets repr\u00e9sentent un standard d&rsquo;excellence m\u00e9thodologique en psychologie cognitive et en neuropsychologie, en raison de leur capacit\u00e9 intrins\u00e8que \u00e0 neutraliser la variabilit\u00e9 interindividuelle r\u00e9siduelle. Dans un protocole pr\u00e9-test versus post-test classique mesurant la flexibilit\u00e9 ex\u00e9cutive avant et apr\u00e8s une s\u00e9ance de rem\u00e9diation cognitive sur une cohorte de 20 sujets, les deux s\u00e9ries de scores sont intrins\u00e8quement d\u00e9pendantes et positivement corr\u00e9l\u00e9es.<\/p>\n<p>L&rsquo;analyse ne s&rsquo;op\u00e8re donc pas sur les distributions transversales brutes, mais sur la variable vectorielle des diff\u00e9rences individuelles, d\u00e9finie par D_i = X_{post, i} &#8211; X_{pre, i}. D\u00e8s lors, l&rsquo;inf\u00e9rence se ram\u00e8ne rigoureusement au comportement d&rsquo;un \u00e9chantillon unique de diff\u00e9rences compar\u00e9 \u00e0 une esp\u00e9rance nulle sous l&rsquo;hypoth\u00e8se de non-efficacit\u00e9 (H0 : mu_D = 0). Supposons que la moyenne des gains cognitifs soit de 4.5 points, avec un \u00e9cart-type des diff\u00e9rences s_D = 6.2. L&rsquo;erreur type de la moyenne des diff\u00e9rences s&rsquo;\u00e9tablit \u00e0 s_D \/ sqrt(N) = 6.2 \/ sqrt(20) = 1.386. La statistique observ\u00e9e vaut ainsi t = 4.5 \/ 1.386 = 3.246 avec df = 19.<\/p>\n<p>L&rsquo;utilisation de la fonction <strong>pt()<\/strong> permet de d\u00e9terminer avec pr\u00e9cision le niveau de certitude probabiliste associ\u00e9 \u00e0 ce gain th\u00e9rapeutique :<\/p>\n<p>p_gain &lt;- pt(3.246, df = 19, lower.tail = FALSE)<\/p>\n<p>La commande retourne une probabilit\u00e9 unilat\u00e9rale de 0.0021, confirmant une am\u00e9lioration clinique hautement significative. Dans ce cadre intrajets, la puissance statistique est d\u00e9cupl\u00e9e par la prise en compte de la corr\u00e9lation test-retest r. L&rsquo;absence de corr\u00e9lation aurait gonfl\u00e9 l&rsquo;erreur type et dilu\u00e9 la statistique t, d\u00e9montrant comment l&rsquo;association fonctionnelle entre la m\u00e9trique du devis exp\u00e9rimental et l&rsquo;int\u00e9gration via <strong>pt()<\/strong> permet d&rsquo;asseoir solidement les conclusions sur la dynamique du changement cognitif.<\/p>\n<h3>6.3 V\u00e9rification des zones de rejet et prise de d\u00e9cision bay\u00e9sienne pr\u00e9liminaire<\/h3>\n<p>La pratique contemporaine de l&rsquo;analyse comportementale invite \u00e0 d\u00e9passer la stricte dichotomisation op\u00e9r\u00e9e par le rituel du seuil alpha arbitraire fix\u00e9 \u00e0 0.05. En neurosciences computationnelles et en psychologie quantitative, l&rsquo;interpr\u00e9tation d&rsquo;une p-valeur issue de <strong>pt()<\/strong> est fr\u00e9quemment enrichie par des approches exploratoires confrontant la vraisemblance fr\u00e9quentiste \u00e0 des approximations de facteurs de Bayes pr\u00e9liminaires.<\/p>\n<p>Lorsqu&rsquo;une statistique t empirique est convertie en p-valeur via <strong>pt()<\/strong>, cette probabilit\u00e9 cumulative r\u00e9siduelle refl\u00e8te la raret\u00e9 des donn\u00e9es sous le mod\u00e8le nul exclusif. Cependant, une p-valeur de 0.045, bien que classiquement qualifi\u00e9e de \u00ab statistiquement significative \u00bb, n&rsquo;apporte souvent qu&rsquo;une \u00e9vidence empirique tr\u00e8s mod\u00e9r\u00e9e en faveur de l&rsquo;hypoth\u00e8se alternative r\u00e9elle. Des m\u00e9thodologistes comme Wagenmakers ont popularis\u00e9 des conversions asymptotiques empiriques permettant d&rsquo;estimer une borne sup\u00e9rieure du facteur de Bayes (maximum putative Bayes factor) directement \u00e0 partir de la valeur de la statistique t et de ses degr\u00e9s de libert\u00e9 :<\/p>\n<p>BF_max &lt;- 1 \/ ( -e * p * log(p) )<\/p>\n<p>En employant judicieusement <strong>pt()<\/strong> pour extraire avec une exactitude d\u00e9cimale irr\u00e9prochable la p-valeur associ\u00e9e \u00e0 l&rsquo;exp\u00e9rience, le chercheur peut instantan\u00e9ment quantifier le rapport de cotes bay\u00e9sien maximal compatible avec ses observations. Si cette estimation calibr\u00e9e indique que les donn\u00e9es ne sont que trois fois plus probables sous H1 que sous H0, le chercheur s&rsquo;abstiendra de proclamer une d\u00e9couverte majeure, adoptant une posture \u00e9pist\u00e9mique prudente et rigoureusement document\u00e9e, en parfaite conformit\u00e9 avec les directives m\u00e9thodologiques les plus r\u00e9centes de l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/www.apa.org\/\">American Psychological Association<\/a>.<\/p>\n<h2>7. La fonction qt() : D\u00e9termination des valeurs critiques et des quantiles<\/h2>\n<h3>7.1 L&rsquo;inversion de la fonction de r\u00e9partition cumulative<\/h3>\n<p>Dans l&rsquo;arsenal d&rsquo;analyse statistique, la fonction <strong>qt()<\/strong> constitue le miroir op\u00e9rationnel parfait de <strong>pt()<\/strong>. Sur le plan de la th\u00e9orie math\u00e9matique, elle mat\u00e9rialise la fonction quantile, d\u00e9finie formellement comme l&rsquo;inverse g\u00e9n\u00e9ralis\u00e9 ou la bijection r\u00e9ciproque de la fonction de r\u00e9partition continue :<\/p>\n<p>Q(p; df) = F^{-1}(p; df) = inf { t in R : F(t; df) &gt;= p }<\/p>\n<p>Sa signature syntaxique officielle dans l&rsquo;environnement R se d\u00e9ploie selon le formalisme suivant :<\/p>\n<p>qt(p, df, ncp, lower.tail = TRUE, log.p = FALSE)<\/p>\n<p>Le premier param\u00e8tre, <em>p<\/em>, attend imp\u00e9rativement une valeur ou un vecteur de probabilit\u00e9s num\u00e9riques born\u00e9es strictement dans l&rsquo;intervalle [0, 1]. Si une valeur ext\u00e9rieure \u00e0 cette plage est inject\u00e9e, l&rsquo;interpr\u00e9teur g\u00e9n\u00e8re imm\u00e9diatement un avertissement et assigne la valeur <em>NaN<\/em> (Not a Number). Le r\u00f4le de <strong>qt()<\/strong> est d&rsquo;extraire la coordonn\u00e9e spatiale exacte sur l&rsquo;axe horizontal qui partitionne la distribution th\u00e9orique de sorte que la fraction d&rsquo;aire situ\u00e9e \u00e0 sa gauche (lorsque <em>lower.tail = TRUE<\/em>) corresponde pr\u00e9cis\u00e9ment \u00e0 l&rsquo;argument <em>p<\/em>.<\/p>\n<p>Sur le plan algorithmique, l&rsquo;inversion d&rsquo;une fonction de r\u00e9partition continue qui n&rsquo;admet pas de forme ferm\u00e9e analytique \u00e9l\u00e9mentaire constitue un probl\u00e8me d&rsquo;analyse num\u00e9rique sophistiqu\u00e9. R mobilise des algorithmes d&rsquo;inversion it\u00e9ratifs de tr\u00e8s haut niveau fond\u00e9s sur des approximations rationnelles de type minimax de Chebyshev, compl\u00e9t\u00e9es par des raffinements de Newton-Raphson. Cette m\u00e9canique garantit une localisation instantan\u00e9e et d&rsquo;une pr\u00e9cision microscopique du quantile recherch\u00e9, permettant d&rsquo;\u00e9tablir les fronti\u00e8res critiques de rejet d&rsquo;hypoth\u00e8ses avec une tol\u00e9rance num\u00e9rique quasi nulle.<\/p>\n<h3>7.2 Identification des seuils critiques unilat\u00e9raux et bilat\u00e9raux<\/h3>\n<p>L&rsquo;utilisation la plus courante de <strong>qt()<\/strong> r\u00e9side dans la d\u00e9termination a priori des valeurs critiques n\u00e9cessaires \u00e0 la prise de d\u00e9cision inf\u00e9rentielle de type Neyman-Pearson. Prenons le cas universel d&rsquo;un chercheur planifiant une comparaison d&rsquo;\u00e9chantillons avec un seuil d&rsquo;erreur de premi\u00e8re esp\u00e8ce alpha fix\u00e9 conventionnellement \u00e0 0.05. La morphologie de la r\u00e9gion critique d\u00e9pend crucialement de l&rsquo;orientation directionnelle du mod\u00e8le statistique.<\/p>\n<p>Pour un test directionnel unilat\u00e9ral droit avec df = 20 degr\u00e9s de libert\u00e9, la surface d&rsquo;erreur critique de 0.05 doit \u00eatre int\u00e9gralement projet\u00e9e dans l&rsquo;extr\u00e9mit\u00e9 sup\u00e9rieure de la courbe. Le quantile de coupure s&rsquo;obtient soit par <em>qt(0.95, df = 20, lower.tail = TRUE)<\/em>, soit par son \u00e9criture duelle plus \u00e9l\u00e9gante <em>qt(0.05, df = 20, lower.tail = FALSE)<\/em>. Le syst\u00e8me d\u00e9livre la valeur critique seuil t_crit = 1.7247. Toute statistique t empirique r\u00e9sultant de l&rsquo;exp\u00e9rience qui surpasserait ce seuil de 1.7247 conduira m\u00e9caniquement au rejet de l&rsquo;hypoth\u00e8se nulle.<\/p>\n<p>Dans la situation standard d&rsquo;un test non directionnel (bilat\u00e9ral), le risque d&rsquo;erreur globale alpha = 0.05 doit \u00eatre scind\u00e9 sym\u00e9triquement en deux demi-aires de superficie alpha \/ 2 = 0.025 assign\u00e9es \u00e0 chacune des extr\u00e9mit\u00e9s de la distribution. La borne critique positive est donc identifi\u00e9e par la requ\u00eate :<\/p>\n<p>t_crit_positif &lt;- qt(p = 1 &#8211; (0.05 \/ 2), df = 20, lower.tail = TRUE)<\/p>\n<p>ou plus directement :<\/p>\n<p>t_crit_positif &lt;- qt(p = 0.025, df = 20, lower.tail = FALSE)<\/p>\n<p>Cette instruction renvoie rigoureusement la valeur critique 2.0860. En raison de la sym\u00e9trie absolue de la distribution centrale, la borne inf\u00e9rieure gauche est simplement son oppos\u00e9 arithm\u00e9tique exact, \u00e0 savoir -2.0860, que l&rsquo;on peut extraire directement par <em>qt(0.025, df = 20, lower.tail = TRUE)<\/em>. Cette d\u00e9marche illustre de fa\u00e7on \u00e9clatante comment <strong>qt()<\/strong> formalise les r\u00e8gles de d\u00e9cision g\u00e9om\u00e9triques avant m\u00eame d&rsquo;avoir recueilli la premi\u00e8re mesure empirique.<\/p>\n<h3>7.3 Correction pour comparaisons multiples : Bonferroni et \u0160id\u00e1k via qt()<\/h3>\n<p>En imagerie c\u00e9r\u00e9brale fonctionnelle, en \u00e9lectrophysiologie (potentiels \u00e9voqu\u00e9s ERP) ou lors de protocoles psychom\u00e9triques mesurant des dizaines de facettes de personnalit\u00e9, les chercheurs effectuent couramment des analyses exploratoires massives impliquant une multiplicit\u00e9 de comparaisons statistiques simultan\u00e9es. La r\u00e9p\u00e9tition effr\u00e9n\u00e9e de tests statistiques univari\u00e9s induit une inflation exponentielle du taux d&rsquo;erreur global de premi\u00e8re esp\u00e8ce par famille de tests (FWER, pour Family-Wise Error Rate). Si l&rsquo;on conduit m = 10 comparaisons ind\u00e9pendantes \u00e0 alpha = 0.05 nominal, la probabilit\u00e9 globale d&rsquo;au moins un faux rejet du hasard s&rsquo;\u00e9l\u00e8ve \u00e0 1 &#8211; (1 &#8211; 0.05)^10 = 0.401, soit plus de 40% de risques de fausse d\u00e9couverte.<\/p>\n<p>Pour restaurer l&rsquo;int\u00e9grit\u00e9 de l&rsquo;inf\u00e9rence collective, les corrections classiques ajustent \u00e0 la baisse le seuil alpha allou\u00e9 \u00e0 chaque test univari\u00e9. La proc\u00e9dure conservative de Bonferroni assigne \u00e0 chaque test le seuil alpha_bonf = alpha \/ m. Gr\u00e2ce \u00e0 <strong>qt()<\/strong>, le chercheur automatise la d\u00e9rivation instantan\u00e9e des nouveaux seuils critiques t p\u00e9nalis\u00e9s sans d\u00e9pendre de tables imprim\u00e9es obsol\u00e8tes :<\/p>\n<p>m &lt;- 10; alpha_global &lt;- 0.05; df_exp &lt;- 28<\/p>\n<p>alpha_ajuste &lt;- alpha_global \/ m<\/p>\n<p>t_crit_bonf &lt;- qt(1 &#8211; (alpha_ajuste \/ 2), df = df_exp)<\/p>\n<p>Ce script calcule un seuil t_crit_bonf de 3.047, substantiellement plus s\u00e9v\u00e8re que le seuil standard non p\u00e9nalis\u00e9 de 2.048 associ\u00e9 \u00e0 alpha = 0.05. De mani\u00e8re encore plus rigoureuse, la correction exacte de \u0160id\u00e1k, affranchie de l&rsquo;in\u00e9galit\u00e9 de Boole pour des hypoth\u00e8ses ind\u00e9pendantes, pose alpha_sidak = 1 &#8211; (1 &#8211; alpha_global)^(1 \/ m). L&rsquo;injection de cette probabilit\u00e9 dans <em>qt(1 &#8211; (alpha_sidak \/ 2), df = df_exp)<\/em> d\u00e9livre le seuil critique d&rsquo;exclusion exact \u00e0 appliquer sur les cartographies neurofonctionnelles, prot\u00e9geant ainsi l&rsquo;exp\u00e9rimentateur contre la prolif\u00e9ration pernicieuse d&rsquo;art\u00e9facts d&rsquo;activation cognitive.<\/p>\n<h2>8. Estimation d&rsquo;intervalles de confiance pour la recherche comportementale avec qt()<\/h2>\n<h3>8.1 Calcul analytique de la marge d&rsquo;erreur autour de la moyenne<\/h3>\n<p>L&rsquo;estimation par intervalle de confiance incarne le socle fondamental de la nouvelle statistique (<em>The New Statistics<\/em>) pr\u00e9conis\u00e9e par les m\u00e9thodologistes de la psychologie contemporaine, dans le but d&rsquo;\u00e9manciper la discipline de la tyrannie st\u00e9rile du test de significativit\u00e9 dichotomique. Un intervalle de confiance \u00e0 95% pour une moyenne de population mu repr\u00e9sente l&rsquo;ensemble des valeurs plausibles du param\u00e8tre qui ne seraient pas rejet\u00e9es au seuil alpha = 0.05 par un test bilat\u00e9ral. Sa d\u00e9rivation analytique exacte repose sur la combinaison lin\u00e9aire de la moyenne empirique observ\u00e9e, de l&rsquo;erreur type de la moyenne (SEM, pour Standard Error of the Mean) et du quantile critique fourni par <strong>qt()<\/strong>.<\/p>\n<p>L&rsquo;\u00e9quation classique r\u00e9gissant les bornes d&rsquo;un intervalle de confiance sous loi de Student s&rsquo;\u00e9tablit selon la double in\u00e9galit\u00e9 suivante :<\/p>\n<p>IC_{1 &#8211; alpha} = [ bar{X} &#8211; t_{(1 &#8211; alpha\/2, df)} times frac{s}{sqrt{N}} quad ; quad bar{X} + t_{(1 &#8211; alpha\/2, df)} times frac{s}{sqrt{N}} ]<\/p>\n<p>Consid\u00e9rons un \u00e9chantillon de psychologie du travail mesurant le niveau de charge mentale per\u00e7ue (sur une \u00e9chelle valid\u00e9e) chez 16 op\u00e9rateurs de contr\u00f4le a\u00e9rien : la moyenne mesur\u00e9e est de 68.4 points, avec un \u00e9cart-type empirique corrig\u00e9 de s = 8.2 points. Le protocole comporte df = 16 &#8211; 1 = 15 degr\u00e9s de libert\u00e9. L&rsquo;extraction computationnelle de la marge d&rsquo;erreur dans R s&rsquo;ex\u00e9cute avec une clart\u00e9 limpide :<\/p>\n<p>n &lt;- 16; moy &lt;- 68.4; s &lt;- 8.2; df_charge &lt;- n &#8211; 1<\/p>\n<p>sem &lt;- s \/ sqrt(n)<\/p>\n<p>t_critique &lt;- qt(0.975, df = df_charge)<\/p>\n<p>marge_erreur &lt;- t_critique * sem<\/p>\n<p>borne_inf &lt;- moy &#8211; marge_erreur; borne_sup &lt;- moy + marge_erreur<\/p>\n<p>Ce script calcule un quantile critique de 2.1314, g\u00e9n\u00e9rant une marge d&rsquo;erreur de 4.369 points, et un intervalle de confiance \u00e0 95% born\u00e9 pr\u00e9cis\u00e9ment entre 64.03 et 72.77 points. Si le chercheur avait eu recours \u00e0 l&rsquo;approximation gaussienne z = 1.96, la marge d&rsquo;erreur calcul\u00e9e aurait \u00e9t\u00e9 artificiellement r\u00e9duite \u00e0 4.018 points, produisant un intervalle faussement \u00e9troit qui occulte l&rsquo;incertitude induite par le petit effectif d&rsquo;\u00e9chantillonnage.<\/p>\n<h3>8.2 Intervalles de confiance pour les tailles d&rsquo;effet (d de Cohen)<\/h3>\n<p>Au-del\u00e0 de la simple estimation de la moyenne m\u00e9trique brute, les directives de reporting m\u00e9thodologique imposent d\u00e9sormais la divulgation syst\u00e9matique des tailles d&rsquo;effet standardis\u00e9es accompagn\u00e9es de leurs intervalles d&rsquo;incertitude d&rsquo;\u00e9chantillonnage respectifs. L&rsquo;indice le plus c\u00e9l\u00e8bre en psychologie diff\u00e9rentielle et exp\u00e9rimentale est le d de Cohen, qui exprime l&rsquo;amplitude d&rsquo;un \u00e9cart de moyennes normalis\u00e9 par la d\u00e9viation standard agr\u00e9g\u00e9e du mod\u00e8le (pooled standard deviation).<\/p>\n<p>Toutefois, une erreur fr\u00e9quente consiste \u00e0 appliquer na\u00efvement les formules d&rsquo;intervalles asymptotiques gaussiens sym\u00e9triques autour du d de Cohen calcul\u00e9. L&rsquo;erreur type d&rsquo;un d de Cohen d&rsquo;\u00e9chantillon d\u00e9pend en r\u00e9alit\u00e9 de la magnitude m\u00eame de cet effet dans la population. L&rsquo;approximation asymptotique classique de son erreur type, formul\u00e9e pour un devis \u00e0 deux groupes ind\u00e9pendants de tailles n1 et n2, s&rsquo;\u00e9tablit traditionnellement selon :<\/p>\n<p>SE_d = sqrt( [ (n1 + n2) \/ (n1 * n2) ] + [ d^2 \/ (2 * (n1 + n2)) ] )<\/p>\n<p>En premi\u00e8re approximation lin\u00e9aire, le quantile extrait par <strong>qt(0.975, df = n1 + n2 &#8211; 2)<\/strong> peut \u00eatre conjugu\u00e9 avec cette erreur type pour borner la variation plausible de la taille d&rsquo;effet. N\u00e9anmoins, comme nous l&rsquo;explorerons de fa\u00e7on exhaustive dans la section d\u00e9di\u00e9e \u00e0 la non-centralit\u00e9, la distribution exacte d&rsquo;un indice de taille d&rsquo;effet n&rsquo;est pas sym\u00e9trique, mais r\u00e9gie par une loi t non-centrale. L&rsquo;usage basique de <strong>qt()<\/strong> standard fournit une estimation pr\u00e9liminaire utile mais asymptotique, qui rappelle la n\u00e9cessit\u00e9 absolue d&rsquo;associer la taille de l&rsquo;effet \u00e0 une quantification rigoureuse de son impr\u00e9cision intrins\u00e8que.<\/p>\n<h3>8.3 Visualisation des barres d&rsquo;erreur bas\u00e9es sur la distribution t<\/h3>\n<p>La repr\u00e9sentation graphique des r\u00e9sultats exp\u00e9rimentaux conditionne la validit\u00e9 des inf\u00e9rences intuitives formul\u00e9es par la communaut\u00e9 scientifique. Trop souvent, les graphiques comportementaux publient des barres d&rsquo;erreur correspondant \u00e0 la simple d\u00e9viation standard (qui ne d\u00e9crit que la variabilit\u00e9 de l&rsquo;\u00e9chantillon sans vis\u00e9e inf\u00e9rentielle) ou \u00e0 l&rsquo;erreur type brute (SEM), laquelle correspond visuellement \u00e0 un intervalle de confiance grossier d&rsquo;environ 68%, cr\u00e9ant l&rsquo;illusion trompeuse d&rsquo;une convergence nette.<\/p>\n<p>Pour construire une visualisation d&rsquo;\u00e9dition rigoureuse au sein de <strong>ggplot2<\/strong>, il est vivement recommand\u00e9 d&rsquo;int\u00e9grer des barres d&rsquo;erreur mat\u00e9rialisant exactement l&rsquo;intervalle de confiance \u00e0 95% calcul\u00e9 sous la loi de Student. Pour ce faire, au lieu de d\u00e9pendre d&rsquo;extensions opaques, le chercheur peut programmer sa propre fonction de r\u00e9sum\u00e9 statistique vectoris\u00e9e :<\/p>\n<p>calcul_ic_t &lt;- function(x) {<br \/>\n&nbsp;&nbsp;m &lt;- mean(x, na.rm = TRUE)<br \/>\n&nbsp;&nbsp;n &lt;- sum(!is.na(x))<br \/>\n&nbsp;&nbsp;sem &lt;- sd(x, na.rm = TRUE) \/ sqrt(n)<br \/>\n&nbsp;&nbsp;marge &lt;- qt(0.975, df = n &#8211; 1) * sem<br \/>\n&nbsp;&nbsp;return(c(y = m, ymin = m &#8211; marge, ymax = m + marge))<br \/>\n}<\/p>\n<p>L&rsquo;injection de cette fonction personnalis\u00e9e dans la couche <em>stat_summary(fun.data = calcul_ic_t, geom = \u00ab\u00a0errorbar\u00a0\u00bb, width = 0.2)<\/em> g\u00e9n\u00e8re des segments d&rsquo;erreur parfaitement calibr\u00e9s sur la th\u00e9orie des petits \u00e9chantillons. Il est fondamental de rappeler aux praticiens le pi\u00e8ge cognitif de l&rsquo;interpr\u00e9tation visuelle : la disjonction ou le non-chevauchement des barres d&rsquo;erreur \u00e0 95% entre deux groupes exp\u00e9rimentaux ind\u00e9pendants n&rsquo;est pas une condition n\u00e9cessaire pour \u00e9tablir une diff\u00e9rence significative \u00e0 p &lt; 0.05. Deux moyennes dont les intervalles de confiance \u00e0 95% bas\u00e9s sur la loi t se chevauchent mod\u00e9r\u00e9ment peuvent parfaitement diff\u00e9rer de mani\u00e8re statistiquement significative lors de l&rsquo;ex\u00e9cution formelle du test t bilat\u00e9ral.<\/p>\n<h2>9. La fonction rt() : Simulation stochastique et g\u00e9n\u00e9ration de nombres pseudo-al\u00e9atoires<\/h2>\n<h3>9.1 Algorithmique du g\u00e9n\u00e9rateur pseudo-al\u00e9atoire sous loi t<\/h3>\n<p>La quatri\u00e8me composante du module, la fonction <strong>rt()<\/strong>, op\u00e8re la transition de l&rsquo;analyse d\u00e9terministe vers la mod\u00e9lisation informatique stochastique. Sa mission algorithmique consiste \u00e0 \u00e9chantillonner des r\u00e9alisations num\u00e9riques ind\u00e9pendantes d&rsquo;une variable al\u00e9atoire gouvern\u00e9e par la distribution de Student. La syntaxe canonique de cette fonction s&rsquo;\u00e9nonce selon :<\/p>\n<p>rt(n, df, ncp)<\/p>\n<p>L&rsquo;argument <em>n<\/em> sp\u00e9cifie le volume scalaire de nombres pseudo-al\u00e9atoires \u00e0 synth\u00e9tiser. Si un vecteur est transmis \u00e0 la place d&rsquo;un scalaire pour <em>n<\/em>, la fonction utilise conventionnellement la longueur de ce vecteur comme taille d&rsquo;\u00e9chantillonnage souhait\u00e9e.<\/p>\n<p>Sur le plan algorithmique sous-jacent, le g\u00e9n\u00e9rateur pseudo-al\u00e9atoire ne recourt pas \u00e0 une inversion directe brutale de la fonction de r\u00e9partition par <em>qt(runif(n), df)<\/em>, qui s&rsquo;av\u00e9rerait num\u00e9riquement trop lente et co\u00fbteuse lors d&rsquo;it\u00e9rations massives. R exploite plut\u00f4t les propri\u00e9t\u00e9s math\u00e9matiques structurelles \u00e9nonc\u00e9es dans la premi\u00e8re section de cet article : il g\u00e9n\u00e8re ind\u00e9pendamment un vecteur de variables normales centr\u00e9es r\u00e9duites Z via l&rsquo;algorithme rapide de Box-Muller ou la m\u00e9thode d&rsquo;inversion normale, puis un vecteur de variables V distribu\u00e9es selon une loi du khi-deux \u00e0 <em>df<\/em> degr\u00e9s de libert\u00e9 (g\u00e9n\u00e9r\u00e9es \u00e0 partir de lois Gamma standard). L&rsquo;algorithme assemble ensuite ces grandeurs stochastiques selon le quotient exact :<\/p>\n<p>T = Z \/ sqrt( V \/ df )<\/p>\n<p>Cette d\u00e9rivation physique garantit une vitesse d&rsquo;ex\u00e9cution optimale et des propri\u00e9t\u00e9s stochastiques rigoureusement indiscernables d&rsquo;une distribution th\u00e9orique parfaite. Pour assurer la stricte reproductibilit\u00e9 scientifique des simulations, il est imp\u00e9ratif d&rsquo;initialiser au pr\u00e9alable le g\u00e9n\u00e9rateur de congruence lin\u00e9aire de R \u00e0 l&rsquo;aide de l&rsquo;instruction <strong>set.seed(graine)<\/strong>. La fixation explicite de cette graine garantit que n&rsquo;importe quel chercheur ex\u00e9cutant le m\u00eame script obtiendra une s\u00e9quence pseudo-al\u00e9atoire strictement identique sur sa machine.<\/p>\n<h3>9.2 Simulation de distributions de donn\u00e9es psychologiques \u00e0 queues lourdes<\/h3>\n<p>L&rsquo;utilit\u00e9 primordiale de <strong>rt()<\/strong> en psychologie exp\u00e9rimentale se manifeste dans la mod\u00e9lisation r\u00e9aliste de donn\u00e9es chronom\u00e9triques ou comportementales qui violent les postulats de la normalit\u00e9 gaussienne. L&rsquo;exemple paradigmatique est constitu\u00e9 par les temps de r\u00e9ponse (RT) dans les paradigmes d&rsquo;inhibition ou d&rsquo;amor\u00e7age attentionnel. Les temps de r\u00e9action r\u00e9els ne sont jamais parfaitement gaussiens : ils pr\u00e9sentent une asym\u00e9trie r\u00e9siduelle positive marqu\u00e9e et une proportion non n\u00e9gligeable d&rsquo;essais extr\u00eamement lents caus\u00e9s par des micro-laps attentionnels passagers.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/tDist2.jpg\" alt=\"Histogram of t distributions\" \/><figcaption>Histogram of t distributions<\/figcaption><\/figure>\n<p>La g\u00e9n\u00e9ration synth\u00e9tique de telles dynamiques attentionnelles atypiques peut \u00eatre simul\u00e9e avec une grande \u00e9l\u00e9gance en injectant des r\u00e9alisations issues de la distribution de Student au moyen de <strong>rt()<\/strong>. En configurant des degr\u00e9s de libert\u00e9 faibles, typiquement compris entre 3 et 6, le chercheur injecte artificiellement une propension robuste aux observations extr\u00eames et atypiques (<em>outliers<\/em>) sans devoir inventer de r\u00e8gles ad hoc complexes. Pour caler ces nombres sans dimension sur une \u00e9chelle temporelle empirique concr\u00e8te en millisecondes, il suffit d&rsquo;appliquer une transformation lin\u00e9aire d&rsquo;\u00e9chelle et de position :<\/p>\n<p>set.seed(42)<\/p>\n<p>temps_reaction_simules &lt;- 450 + ( 65 * rt(n = 500, df = 4) )<\/p>\n<p>L&rsquo;observation du vecteur r\u00e9sultant met en \u00e9vidence une moyenne g\u00e9n\u00e9rale centr\u00e9e autour de 450 ms, mais parsem\u00e9e de plusieurs temps de r\u00e9ponse exc\u00e9dant 700 ou 800 ms, reproduisant fid\u00e8lement la signature empirique observ\u00e9e chez des participants humains dans des t\u00e2ches attentionnelles soutenues. La confrontation graphique de cet \u00e9chantillon simul\u00e9 par <strong>rt()<\/strong> avec un \u00e9chantillon de taille \u00e9quivalente g\u00e9n\u00e9r\u00e9 par la fonction gaussienne <em>rnorm()<\/em> r\u00e9v\u00e8le imm\u00e9diatement la capacit\u00e9 de la loi de Student \u00e0 mod\u00e9liser organiquement l&rsquo;incertitude \u00e9pist\u00e9mique et l&rsquo;instabilit\u00e9 comportementale naturelle.<\/p>\n<h3>9.3 Exp\u00e9riences de Monte-Carlo pour tester la robustesse des estimateurs<\/h3>\n<p>L&rsquo;exp\u00e9rimentation de Monte-Carlo repr\u00e9sente le laboratoire virtuel supr\u00eame du statisticien. Elle permet de soumettre les proc\u00e9dures d&rsquo;inf\u00e9rence \u00e0 des conditions exp\u00e9rimentales de stress extr\u00eame afin d&rsquo;\u00e9valuer empiriquement leur r\u00e9sistance aux violations m\u00e9thodologiques. Un questionnement r\u00e9current en psychologie porte sur la robustesse du test t classique \u00e0 deux \u00e9chantillons ind\u00e9pendants face \u00e0 la violation de l&rsquo;hypoth\u00e8se de normalit\u00e9 : l&rsquo;inflation de l&rsquo;erreur de type I (taux de faux rejets) est-elle r\u00e9ellement contenue lorsque les populations r\u00e9elles pr\u00e9sentent des queues lourdes ?<\/p>\n<p>La fonction <strong>rt()<\/strong> permet de concevoir une telle \u00e9tude de Monte-Carlo en quelques lignes de script vectoris\u00e9 exploitant la primitive <strong>replicate()<\/strong>. Le protocole informatique g\u00e9n\u00e8re it\u00e9rativement des milliers d&rsquo;exp\u00e9riences synth\u00e9tiques pr\u00e9levant deux groupes ind\u00e9pendants de petite taille (par exemple n1 = n2 = 10) sous une loi t \u00e0 df = 3 degr\u00e9s de libert\u00e9 sous stricte hypoth\u00e8se nulle (aucune diff\u00e9rence de moyennes r\u00e9elles), puis enregistre la proportion empirique de rejets \u00e0 alpha = 0.05 :<\/p>\n<p>set.seed(123)<\/p>\n<p>B &lt;- 10000<\/p>\n<p>resultats_p &lt;- replicate(B, {<br \/>\n&nbsp;&nbsp;groupe1 &lt;- rt(n = 10, df = 3)<br \/>\n&nbsp;&nbsp;groupe2 &lt;- rt(n = 10, df = 3)<br \/>\n&nbsp;&nbsp;t.test(groupe1, groupe2, var.equal = TRUE)$p.value<br \/>\n})<\/p>\n<p>taux_erreur_type_I &lt;- mean(resultats_p &lt; 0.05)<\/p>\n<p>L&rsquo;ex\u00e9cution de cette simulation massive sur 10 000 r\u00e9plications d\u00e9livre un taux empirique de rejet oscillant g\u00e9n\u00e9ralement autour de 0.047 \u00e0 0.051. Cette observation computationnelle rigoureuse apporte la d\u00e9monstration \u00e9clatante que le test t \u00e0 deux \u00e9chantillons \u00e9quilibr\u00e9s de Student est remarquablement robuste \u00e0 la leptokurticit\u00e9 des queues de distribution, dissipant les angoisses infond\u00e9es relatives aux d\u00e9viations mod\u00e9r\u00e9es de la normalit\u00e9 et illustrant la puissance de <strong>rt()<\/strong> comme outil d&rsquo;investigation m\u00e9ta-scientifique.<\/p>\n<h2>10. Gestion avanc\u00e9e du param\u00e8tre de non-centralit\u00e9 (ncp)<\/h2>\n<h3>10.1 Fondements de la loi t non-centrale<\/h3>\n<p>L&rsquo;immense majorit\u00e9 des manuels d&rsquo;introduction \u00e0 la statistique restreignent leur expos\u00e9 \u00e0 la loi de Student centrale, o\u00f9 la statistique t oscille sym\u00e9triquement autour d&rsquo;une esp\u00e9rance nulle en l&rsquo;absence de tout effet exp\u00e9rimental sous H0. D\u00e8s lors que l&rsquo;on s&rsquo;int\u00e9resse \u00e0 l&rsquo;hypoth\u00e8se alternative (H1), postulant que l&rsquo;effet exp\u00e9rimental dans la population sous-jacente est non nul, la variable al\u00e9atoire ne suit plus une loi centrale sym\u00e9trique. Elle bascule sous le r\u00e9gime plus g\u00e9n\u00e9ral de la loi de Student non-centrale, dont la structure est gouvern\u00e9e par le param\u00e8tre de non-centralit\u00e9 delta, d\u00e9sign\u00e9 dans R par l&rsquo;argument <em>ncp<\/em>.<\/p>\n<p>Sur le plan de sa d\u00e9rivation formelle, la variable t non-centrale est d\u00e9finie par le quotient :<\/p>\n<p>T_{non-centrale} = ( Z + delta ) \/ sqrt{ V \/ df }<\/p>\n<p>o\u00f9 Z repr\u00e9sente une variable normale centr\u00e9e r\u00e9duite, V une variable ind\u00e9pendante distribu\u00e9e selon une loi du khi-deux \u00e0 <em>df<\/em> degr\u00e9s de libert\u00e9, et delta est un nombre r\u00e9el constant d\u00e9calant l&rsquo;esp\u00e9rance du num\u00e9rateur. Dans le contexte classique d&rsquo;une comparaison de deux moyennes ind\u00e9pendantes avec des groupes d&rsquo;effectifs n1 et n2, le param\u00e8tre delta est directement proportionnel \u00e0 la taille d&rsquo;effet standardis\u00e9e de la population (le delta de Cohen sous-jacent) :<\/p>\n<p>delta = d_{population} times sqrt{ frac{n_1 times n_2}{n_1 + n_2} }<\/p>\n<p>Dans le cas d&rsquo;un \u00e9chantillon unique ou de mesures r\u00e9p\u00e9t\u00e9es de taille N, la relation se simplifie en delta = d_{pop} * sqrt(N). L&rsquo;adjonction de ce param\u00e8tre modifie radicalement la g\u00e9om\u00e9trie de la densit\u00e9 de distribution : la courbe perd sa sym\u00e9trie, subit une inclinaison asym\u00e9trique (skewness) marqu\u00e9e dans la direction de delta, et son mode ne co\u00efncide plus exactement avec la valeur de delta en raison de l&rsquo;\u00e9crasement exerc\u00e9 par la division stochastique du d\u00e9nominateur.<\/p>\n<h3>10.2 Calcul de la puissance statistique a priori et a posteriori<\/h3>\n<p>La loi t non-centrale constitue le moteur math\u00e9matique exclusif de l&rsquo;analyse de puissance statistique dans les devis param\u00e9triques univari\u00e9s. La puissance d&rsquo;un test statistique (1 &#8211; b\u00eata) repr\u00e9sente la probabilit\u00e9 de rejeter formellement l&rsquo;hypoth\u00e8se nulle lorsque l&rsquo;hypoth\u00e8se alternative H1 est authentiquement vraie dans la nature. Sur le plan g\u00e9om\u00e9trique, la puissance correspond \u00e0 l&rsquo;int\u00e9grale de la distribution t non-centrale au-del\u00e0 des bornes critiques fix\u00e9es sous la distribution centrale.<\/p>\n<p>L&rsquo;utilisation combin\u00e9e des fonctions <strong>qt()<\/strong> et <strong>pt()<\/strong> permet de programmer un moteur de calcul de puissance sur-mesure d&rsquo;une concision exemplaire, s&rsquo;affranchissant totalement de la d\u00e9pendance \u00e0 des progiciels sp\u00e9cialis\u00e9s comme G*Power ou la biblioth\u00e8que <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=pwr\">pwr<\/a>. Prenons un sc\u00e9nario de recherche o\u00f9 un psychologue planifie une exp\u00e9rience avec deux groupes de N = 25 sujets chacun (df = 25 + 25 &#8211; 2 = 48). Il anticipe une taille d&rsquo;effet moyenne de d = 0.50 au seuil de significativit\u00e9 bilat\u00e9ral alpha = 0.05 :<\/p>\n<p>n1 &lt;- 25; n2 &lt;- 25; df_puiss &lt;- n1 + n2 &#8211; 2; d_effet &lt;- 0.50<\/p>\n<p>seuil_critique &lt;- qt(1 &#8211; (0.05 \/ 2), df = df_puiss)<\/p>\n<p>param_ncp &lt;- d_effet * sqrt((n1 * n2) \/ (n1 + n2))<\/p>\n<p>puissance_calculee &lt;- pt(seuil_critique, df = df_puiss, ncp = param_ncp, lower.tail = FALSE)<\/p>\n<p>Ce script calcule d&rsquo;abord le seuil critique sym\u00e9trique (t_crit = 2.0106) au moyen de <strong>qt()<\/strong> central, puis d\u00e9rive le param\u00e8tre ncp = 0.50 * sqrt(625 \/ 50) = 1.7677. Enfin, <strong>pt()<\/strong> int\u00e8gre la queue sup\u00e9rieure de la distribution non-centrale au-del\u00e0 de 2.0106, retournant une puissance exacte de 0.4038 (soit environ 40.4% de chances de d\u00e9tecter l&rsquo;effet exp\u00e9rimental r\u00e9el). Cette d\u00e9monstration limpide met imm\u00e9diatement en exergue le sous-dimensionnement critique de nombreuses exp\u00e9riences en psychologie cognitive, o\u00f9 des cohortes restreintes de 25 sujets par groupe \u00e9chouent \u00e0 garantir le standard de puissance conventionnel de 80% face \u00e0 des effets mod\u00e9r\u00e9s.<\/p>\n<h3>10.3 Construction d&rsquo;intervalles de confiance exacts non-centraux<\/h3>\n<p>L&rsquo;apog\u00e9e de la rigueur statistique dans le traitement des tailles d&rsquo;effet standardis\u00e9es consiste \u00e0 calculer l&rsquo;intervalle de confiance exact du d de Cohen par l&rsquo;inversion de la distribution t non-centrale, une m\u00e9thodologie formalis\u00e9e par Michael Smithson et Geoff Cumming. Contrairement aux approximations asymptotiques sym\u00e9triques qui violent la morphologie asym\u00e9trique de la distribution d&rsquo;\u00e9chantillonnage de d, cette m\u00e9thode identifie les deux param\u00e8tres de non-centralit\u00e9 uniques qui placent la statistique t empiriquement observ\u00e9e exactement sur leurs quantiles critiques respectifs.<\/p>\n<p>Le principe computationnel proc\u00e8de par recherche num\u00e9rique de z\u00e9ros au moyen de la fonction d&rsquo;optimisation univari\u00e9e <strong>uniroot()<\/strong> disponible nativement dans R. Supposons qu&rsquo;un devis ind\u00e9pendant \u00e0 deux groupes de 20 participants chacun (df = 38) produise une statistique observ\u00e9e t_obs = 2.85 (d d&rsquo;\u00e9chantillon = 2.85 * sqrt(40 \/ 400) = 0.901). Les bornes inf\u00e9rieure et sup\u00e9rieure du param\u00e8tre ncp satisfont les \u00e9quations probabilistes d&rsquo;inversion :<\/p>\n<p>pt(t_obs, df = 38, ncp = ncp_inf) = 0.975<\/p>\n<p>pt(t_obs, df = 38, ncp = ncp_sup) = 0.025<\/p>\n<p>Dans R, la r\u00e9solution num\u00e9rique s&rsquo;op\u00e8re instantan\u00e9ment via le code suivant :<\/p>\n<p>trouver_ncp_inf &lt;- function(ncp) pt(2.85, df = 38, ncp = ncp) &#8211; 0.975<\/p>\n<p>trouver_ncp_sup &lt;- function(ncp) pt(2.85, df = 38, ncp = ncp) &#8211; 0.025<\/p>\n<p>borne_ncp_inf &lt;- uniroot(trouver_ncp_inf, interval = c(-5, 10))$root<\/p>\n<p>borne_ncp_sup &lt;- uniroot(trouver_ncp_sup, interval = c(-5, 10))$root<\/p>\n<p>Une fois les bornes de non-centralit\u00e9 isol\u00e9es (par exemple ncp_inf = 0.778 et ncp_sup = 4.882), il suffit de multiplier ces grandeurs par sqrt((n1 + n2) \/ (n1 * n2)) pour obtenir les bornes exactes et asym\u00e9triques de l&rsquo;intervalle de confiance \u00e0 95% du d de Cohen dans la population : [0.246 ; 1.544]. Cette proc\u00e9dure \u00e9limine toute distorsion math\u00e9matique et garantit une fid\u00e9lit\u00e9 d&rsquo;inf\u00e9rence absolue, incarnant le summum de l&rsquo;exigence m\u00e9thodologique en psychologie exp\u00e9rimentale.<\/p>\n<h2>11. Loi normale vs Loi t : \u00c9tude comparative rigoureuse en R<\/h2>\n<h3>11.1 Analyse comparative formelle : dnorm\/dt, pnorm\/pt, qnorm\/qt, rnorm\/rt<\/h3>\n<p>Pour parfaire la compr\u00e9hension op\u00e9rationnelle du chercheur, il est indispensable de structurer une confrontation syst\u00e9matique entre les quatre fonctions de la famille t et leurs homologues de la famille normale : <strong>dnorm<\/strong>, <strong>pnorm<\/strong>, <strong>qnorm<\/strong> et <strong>rnorm<\/strong>. Bien que ces deux familles partagent des signatures syntaxiques tr\u00e8s proches et une m\u00eame vocation \u00e0 mod\u00e9liser la variabilit\u00e9 unimodale sym\u00e9trique, leurs comportements algorithmiques divergent substantiellement sous les contraintes des petits \u00e9chantillons.<\/p>\n<p>Consid\u00e9rons en premier lieu le couple <strong>dnorm\/dt<\/strong>. La fonction normale d\u00e9pend exclusivement d&rsquo;un param\u00e8tre d&rsquo;\u00e9chelle constant (l&rsquo;\u00e9cart-type sigma), conf\u00e9rant \u00e0 sa courbure une d\u00e9croissance quadratique exponentielle rapide : exp(-x^2 \/ 2). \u00c0 l&rsquo;inverse, <strong>dt<\/strong> int\u00e8gre une d\u00e9croissance polynomiale mod\u00e9r\u00e9e proportionnelle \u00e0 (1 + x^2 \/ df)^(-(df+1)\/2). Cette divergence g\u00e9om\u00e9trique peut \u00eatre formellement quantifi\u00e9e en th\u00e9orie de l&rsquo;information par la divergence de Kullback-Leibler (KL). Lorsque df = 2, la divergence KL entre la loi t et la loi normale est substantielle, traduisant une perte massive d&rsquo;information probabiliste si l&rsquo;on tente de mod\u00e9liser des queues lourdes par une courbe gaussienne. Cette discr\u00e9pance ne devient infinit\u00e9simale qu&rsquo;au-del\u00e0 de df &gt; 50.<\/p>\n<p>Le diptyque <strong>qnorm\/qt<\/strong> met en lumi\u00e8re des divergences pratiques imm\u00e9diates sur les seuils d&rsquo;exclusion. Alors que <em>qnorm(0.975)<\/em> retourne la constante universelle 1.95996 quel que soit l&rsquo;\u00e9chantillon, <em>qt(0.975, df = 3)<\/em> s&rsquo;envole \u00e0 3.1824. Enfin, sur le plan g\u00e9n\u00e9ratif, la confrontation stochastique de <strong>rnorm<\/strong> et <strong>rt<\/strong> montre que <strong>rt()<\/strong> produit une dispersion d&rsquo;\u00e9chantillonnage marqu\u00e9e par une surrepr\u00e9sentation d&rsquo;observations extr\u00eames \u00e0 plus de trois \u00e9carts-types, alors que <strong>rnorm()<\/strong> n&rsquo;en g\u00e9n\u00e8re statistiquement qu&rsquo;environ 27 pour 10 000 tirages. Cette analyse crois\u00e9e formalise de mani\u00e8re incontestable l&rsquo;irrempla\u00e7abilit\u00e9 du quatuor de Student d\u00e8s lors que le param\u00e8tre de variance de population n&rsquo;est pas accessible a priori.<\/p>\n<h3>11.2 Cons\u00e9quences de l&rsquo;assimilation abusive de la loi t \u00e0 la loi normale<\/h3>\n<p>L&rsquo;assimilation aveugle et abusive de la distribution de Student \u00e0 la distribution normale constitue l&rsquo;un des vices m\u00e9thodologiques les plus destructeurs r\u00e9pertori\u00e9s dans l&rsquo;histoire des publications en sciences biom\u00e9dicales et psychologiques. Durant des d\u00e9cennies, des chercheurs ont substitu\u00e9 les quantiles gaussiens aux quantiles de Student, soit par commodit\u00e9 algorithmique \u00e0 l&rsquo;\u00e9poque o\u00f9 les puissances de calcul \u00e9taient rudimentaires, soit par m\u00e9connaissance des fondements math\u00e9matiques de l&rsquo;inf\u00e9rence.<\/p>\n<p>Les cons\u00e9quences analytiques de cette substitution indue sur les cohortes cliniques restreintes sont dramatiques. En utilisant <strong>qnorm()<\/strong> \u00e0 la place de <strong>qt()<\/strong> sur un \u00e9chantillon de taille n = 5 (df = 4), le chercheur applique un seuil critique bilat\u00e9ral de 1.96 au lieu du seuil rigoureux de 2.776. Si la statistique t empirique observ\u00e9e se monte \u00e0 2.30, l&rsquo;analyste abus\u00e9 par l&rsquo;approximation normale d\u00e9clarera le r\u00e9sultat statistiquement significatif au seuil alpha = 0.05 (car 2.30 &gt; 1.96), alors que la p-valeur exacte calcul\u00e9e par <em>2 * pt(2.30, df = 4, lower.tail = FALSE)<\/em> est en r\u00e9alit\u00e9 de 0.083, relevant purement du hasard sous les conventions scientifiques \u00e9tablies. Ce ph\u00e9nom\u00e8ne engendre une inflation d\u00e9vastatrice des erreurs de premi\u00e8re esp\u00e8ce (faux positifs), polluant la litt\u00e9rature scientifique d&rsquo;effets illusoires impossibles \u00e0 r\u00e9pliquer ult\u00e9rieurement.<\/p>\n<p>De surcro\u00eet, cette approximation erron\u00e9e engendre un r\u00e9tr\u00e9cissement trompeur des intervalles de confiance calcul\u00e9s dans les \u00e9tudes de neuro-imagerie fonctionnelle et d&rsquo;\u00e9valuation neuropsychologique. Les marges d&rsquo;incertitude \u00e9tant sous-\u00e9valu\u00e9es, les investigateurs formulent des assertions cliniques d\u00e9mesur\u00e9ment optimistes quant \u00e0 la taille v\u00e9ritable des d\u00e9ficits cognitifs ou \u00e0 la pr\u00e9cision des th\u00e9rapies cibl\u00e9es. La relecture historique des crises de reproductibilit\u00e9 confirme que l&rsquo;oubli syst\u00e9matique de la correction de Student sur les \u00e9chantillons modestes a largement contribu\u00e9 aux errements diagnostiques et inf\u00e9rentiels document\u00e9s au cours des derni\u00e8res d\u00e9cennies.<\/p>\n<h3>11.3 Recommandations m\u00e9thodologiques pour le choix du mod\u00e8le de r\u00e9f\u00e9rence<\/h3>\n<p>Face \u00e0 ces p\u00e9rils d&rsquo;inf\u00e9rence, la communaut\u00e9 des m\u00e9thodologistes contemporains a \u00e9dict\u00e9 des arbres de d\u00e9cision stricts et formalis\u00e9s pour guider le choix du mod\u00e8le de distribution de r\u00e9f\u00e9rence dans le traitement des jeux de donn\u00e9es quantitatifs. La r\u00e8gle d&rsquo;or universelle \u00e9nonce que d\u00e8s lors que l&rsquo;\u00e9cart-type de la population sigma est inconnu \u2014 ce qui repr\u00e9sente la quasi-totalit\u00e9 des situations exp\u00e9rimentales concr\u00e8tes en psychologie \u2014, l&rsquo;usage de la loi de Student via ses fonctions d\u00e9di\u00e9es <strong>dt, pt, qt, rt<\/strong> est formellement obligatoire, quelle que soit la taille de l&rsquo;\u00e9chantillon d&rsquo;analyse.<\/p>\n<p>L&rsquo;argument fallacieux pr\u00e9tendant que \u00ab pour N sup\u00e9rieur \u00e0 30, la loi normale s&rsquo;applique \u00bb doit \u00eatre d\u00e9finitivement banni de la formation des jeunes chercheurs. Si la convergence asymptotique rapproche effectivement les grandeurs num\u00e9riques pour des \u00e9chantillons \u00e9lev\u00e9s, l&rsquo;application continue de la distribution t dans R n&rsquo;impose aucun co\u00fbt de calcul perceptible tout en garantissant une exactitude math\u00e9matique int\u00e9grale \u00e0 tout effectif d&rsquo;\u00e9chantillonnage. En outre, le choix du mod\u00e8le statistique doit incorporer une phase pr\u00e9alable d&rsquo;exploration rigoureuse au moyen de tests de normalit\u00e9 formels, tels que le test de Shapiro-Wilk (<em>shapiro.test()<\/em>), combin\u00e9s \u00e0 l&rsquo;inspection visuelle des diagrammes quantile-quantile (Q-Q plots).<\/p>\n<p>Lorsque les donn\u00e9es empiriques exhibent des asym\u00e9tries prononc\u00e9es ou une pr\u00e9sence irr\u00e9ductible de valeurs extr\u00eames violant le postulat d&rsquo;\u00e9chantillonnage sous-jacent d&rsquo;une population parente normale, le chercheur doit basculer vers des architectures d&rsquo;analyse robustes. Parmi celles-ci figurent les tests bas\u00e9s sur les moyennes tronqu\u00e9es (trimmed means) impl\u00e9ment\u00e9s dans le package <em>WRS2<\/em> de Rand Wilcox, les proc\u00e9dures de r\u00e9\u00e9chantillonnage bootstrap non param\u00e9triques, ou l&rsquo;ajustement de mod\u00e8les lin\u00e9aires mixtes (LMM) via <em>lme4<\/em> lorsque les structures d&rsquo;observation sont embo\u00eet\u00e9es. Dans tous les cas, l&rsquo;ancrage dans la distribution de Student demeure le pivot central autour duquel s&rsquo;organise l&rsquo;appr\u00e9ciation critique de l&rsquo;\u00e9vidence empirique.<\/p>\n<h2>12. Bonnes pratiques, pi\u00e8ges computationnels et reproductibilit\u00e9 scientifique<\/h2>\n<h3>12.1 Les erreurs syntaxiques et conceptuelles classiques dans l&rsquo;usage des fonctions<\/h3>\n<p>L&rsquo;usage op\u00e9rationnel du quatuor des fonctions t dans la pratique quotidienne des laboratoires expose les analystes \u00e0 une s\u00e9rie d&rsquo;erreurs r\u00e9currentes, \u00e0 la fois syntaxiques et conceptuelles, qu&rsquo;il convient de r\u00e9pertorier m\u00e9thodiquement afin d&rsquo;en immuniser l&rsquo;exp\u00e9rimentateur. L&rsquo;une des confusions les plus banales r\u00e9side dans l&rsquo;interversion des arguments cardinaux entre <strong>pt()<\/strong> et <strong>qt()<\/strong> : injecter une probabilit\u00e9 (par exemple 0.05) comme premier argument dans <strong>pt()<\/strong> ou, inversement, transmettre une coordonn\u00e9e spatiale (t = 2.5) \u00e0 <strong>qt()<\/strong> sans g\u00e9n\u00e9rer d&rsquo;erreur formelle si les grandeurs sont comprises dans l&rsquo;intervalle [0, 1], conduisant \u00e0 des inf\u00e9rences totalement aberrantes.<\/p>\n<p>Une autre b\u00e9vue computationnelle majeure concerne l&rsquo;omission ou la mauvaise orientation de l&rsquo;argument <em>lower.tail<\/em>. L&rsquo;\u00e9valuation irr\u00e9fl\u00e9chie d&rsquo;une statistique positive unilat\u00e9rale droite par <em>pt(2.5, df = 10)<\/em> renvoie une probabilit\u00e9 de 0.987, qui correspond \u00e0 la surface cumulative gauche ; l&rsquo;analyste inattentif pourrait d\u00e9clarer le test totalement non significatif en confondant l&rsquo;aire de confiance centrale avec la p-valeur r\u00e9siduelle de queue sup\u00e9rieure. De m\u00eame, dans les tests bilat\u00e9raux, l&rsquo;oubli de doubler la probabilit\u00e9 extraite de la queue unilat\u00e9rale g\u00e9n\u00e8re une p-valeur divis\u00e9e par deux, doublant artificiellement le taux r\u00e9el d&rsquo;erreurs de type I du protocole.<\/p>\n<p>Enfin, la sp\u00e9cification d\u00e9fectueuse des degr\u00e9s de libert\u00e9 repr\u00e9sente un foyer d&rsquo;erreurs chroniques. Dans un devis \u00e0 mesures r\u00e9p\u00e9t\u00e9es intrajets pr\u00e9-test versus post-test portant sur N = 15 sujets (conduisant \u00e0 30 mesures brutes appari\u00e9es), renseigner par m\u00e9garde <em>df = 28<\/em> (comme s&rsquo;il s&rsquo;agissait de deux groupes ind\u00e9pendants de 15 sujets) au lieu de <em>df = 14<\/em> (N &#8211; 1 diff\u00e9rences) gonfle ind\u00fbment les degr\u00e9s de libert\u00e9, faussant la d\u00e9rivation des seuils critiques et la validit\u00e9 de la d\u00e9cision statistique. Une vigilance syntaxique scrupuleuse est donc indispensable \u00e0 chaque \u00e9tape du codage.<\/p>\n<h3>12.2 Normes de rapportage APA 7 pour les statistiques t issues de R<\/h3>\n<p>La transcription des analyses quantitatives issues de l&rsquo;environnement R dans les publications scientifiques exige une conformit\u00e9 absolue avec les canons r\u00e9dactionnels de la 7e \u00e9dition du manuel de style de l&rsquo;APA. Selon ces directives internationales, la communication d&rsquo;un r\u00e9sultat d&rsquo;inf\u00e9rence bas\u00e9 sur la loi de Student ne saurait se r\u00e9duire \u00e0 une mention lapidaire de significativit\u00e9 binaire ; elle doit pr\u00e9senter l&rsquo;ensemble des param\u00e8tres articul\u00e9s de mani\u00e8re formalis\u00e9e et transparente.<\/p>\n<p>Le standard d&rsquo;\u00e9criture typographique impose d&rsquo;indiquer la lettre t en italique minuscule, suivie imm\u00e9diatement entre parenth\u00e8ses des degr\u00e9s de libert\u00e9 effectifs non arrondis (ou arrondis \u00e0 deux d\u00e9cimales dans le cas d&rsquo;ajustements continus de Welch), de la valeur empirique exacte de la statistique de test, du niveau exact de la p-valeur, et de la m\u00e9trique de taille d&rsquo;effet standardis\u00e9e escort\u00e9e de son intervalle de confiance \u00e0 95% calcul\u00e9 sous distribution exacte. La structure narrative canonique s&rsquo;organise selon le sch\u00e9ma d&rsquo;\u00e9dition suivant :<\/p>\n<p><em>t<\/em>(28) = 2.45, <em>p<\/em> = .021, <em>d<\/em> = 0.89, IC 95% [0.13, 1.64]<\/p>\n<p>Il convient de souligner deux subtilit\u00e9s formelles prescrites par l&rsquo;APA : la valeur p ne doit jamais \u00eatre pr\u00e9c\u00e9d\u00e9e d&rsquo;un z\u00e9ro initial lorsqu&rsquo;elle est born\u00e9e entre 0 et 1 (\u00e9crire <em>p<\/em> = .021 et non <em>p<\/em> = 0.021), et d\u00e8s lors que l&rsquo;algorithme de R retourne une probabilit\u00e9 infinit\u00e9simale inf\u00e9rieure au seuil de pr\u00e9cision standard, elle doit \u00eatre rapport\u00e9e sous la forme <em>p<\/em> &lt; .001 plut\u00f4t que <em>p<\/em> = .000. Enfin, l&rsquo;auteur a l&rsquo;obligation \u00e9thique d&rsquo;indiquer explicitement si le test statistique ex\u00e9cut\u00e9 \u00e9tait directionnel (unilat\u00e9ral) ou non directionnel (bilat\u00e9ral), garantissant une reproductibilit\u00e9 cognitive parfaite de la lecture du manuscrit.<\/p>\n<h3>12.3 Structuration d&rsquo;un script d&rsquo;analyse reproductible et auto-document\u00e9<\/h3>\n<p>L&rsquo;id\u00e9al moderne de la science ouverte (<em>Open Science<\/em>) postule que toute publication scientifique doit \u00eatre accompagn\u00e9e de son mat\u00e9riel m\u00e9thodologique complet et de son code d&rsquo;analyse source v\u00e9rifiable. Pour garantir une reproductibilit\u00e9 computationnelle sans faille de l&rsquo;\u00e9valuation des fonctions <strong>dt, pt, qt, rt<\/strong>, le script R doit \u00eatre architectur\u00e9 de mani\u00e8re modulaire, auto-document\u00e9e et ind\u00e9pendante de l&rsquo;\u00e9tat volatile de la session locale de l&rsquo;analyste.<\/p>\n<p>Une bonne pratique fondamentale consiste \u00e0 scapsuler les analyses au sein d&rsquo;environnements programmatiques reproductibles fond\u00e9s sur des documents informatiques dynamiques tels que <a rel=\"noopener noreferrer\" href=\"https:\/\/rmarkdown.rstudio.com\/\">R Markdown<\/a> ou sa version contemporaine <a rel=\"noopener noreferrer\" href=\"https:\/\/quarto.org\/\">Quarto<\/a>. Ces canevas informatiques unifient le texte acad\u00e9mique r\u00e9dig\u00e9, le code ex\u00e9cutable et les sorties graphiques ou num\u00e9riques en un artefact unique infalsifiable. Au sommet du document, un bloc d&rsquo;initialisation rigoureux configure les graines pseudo-al\u00e9atoires n\u00e9cessaires \u00e0 la stabilit\u00e9 des fonctions g\u00e9n\u00e9ratives (<em>set.seed<\/em>), verrouille la version des d\u00e9pendances logicielles via des outils d&rsquo;instanciation comme <em>renv<\/em>, et charge les biblioth\u00e8ques d&rsquo;\u00e9dition requises.<\/p>\n<p>Chaque manipulation des fonctions t doit \u00eatre assortie de commentaires descriptifs pr\u00e9cisant l&rsquo;intention m\u00e9thodologique sous-jacente : formulation de l&rsquo;hypoth\u00e8se nulle, justification des degr\u00e9s de libert\u00e9 retenus, et tra\u00e7abilit\u00e9 des transformations alg\u00e9briques conduisant aux indices de taille d&rsquo;effet. Enfin, l&rsquo;archivage p\u00e9renne du script sur des d\u00e9p\u00f4ts institutionnels publics certifi\u00e9s, tels que l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/osf.io\/\">Open Science Framework<\/a> (OSF) ou Zenodo, permet \u00e0 la communaut\u00e9 des pairs de r\u00e9ex\u00e9cuter l&rsquo;int\u00e9gralit\u00e9 du pipeline computationnel au bit pr\u00e8s, consolidant l&rsquo;\u00e9difice scientifique collectif et incarnant l&rsquo;\u00e9tat de l&rsquo;art de l&rsquo;investigation quantitative en psychologie 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>Chambers, J. M. (2008). <em>Software for data analysis: Programming with R<\/em>. Springer. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-0-387-75936-4\">https:\/\/doi.org\/10.1007\/978-0-387-75936-4<\/a><\/li>\n<li>Cohen, J. (1988). <em>Statistical power analysis for the behavioral sciences<\/em> (2nd ed.). Lawrence Erlbaum Associates.<\/li>\n<li>Cumming, G. (2012). <em>Understanding the new statistics: Effect sizes, confidence intervals, and meta-analysis<\/em>. Routledge. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.4324\/9780203807002\">https:\/\/doi.org\/10.4324\/9780203807002<\/a><\/li>\n<li>Fisher, R. A. (1925). <em>Statistical methods for research workers<\/em>. Oliver and Boyd.<\/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>Smithson, M. (2003). <em>Confidence intervals<\/em>. SAGE Publications. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.4135\/9781412983761\">https:\/\/doi.org\/10.4135\/9781412983761<\/a><\/li>\n<li>Student. (1908). The probable error of a mean. <em>Biometrika<\/em>, 6(1), 1\u201325. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/6.1.1\">https:\/\/doi.org\/10.1093\/biomet\/6.1.1<\/a><\/li>\n<li>Wagenmakers, E.-J. (2007). A practical solution to the pervasive problems of p values. <em>Psychonomic Bulletin &amp; Review<\/em>, 14(5), 779\u2013804. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.3758\/BF03194105\">https:\/\/doi.org\/10.3758\/BF03194105<\/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>Wilcox, R. R. (2017). <em>Introduction to robust estimation and hypothesis testing<\/em> (4th ed.). Academic Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/C2015-0-02263-2\">https:\/\/doi.org\/10.1016\/C2015-0-02263-2<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour ma\u00eetriser les fonctions dt, qt, pt et rt dans R, appliqu\u00e9es \u00e0 l&rsquo;analyse statistique et aux donn\u00e9es en psychologie exp\u00e9rimentale.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2193],"tags":[3503,1940,2956,2886,3502],"class_list":["post-1669","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-programmation-r","tag-fonctions-statistiques-r","tag-loi-de-student","tag-programmation-r","tag-statistiques-psychologie","tag-tests-d-hypotheses"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1669","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=1669"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1669\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1669"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1669"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1669"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}