{"id":1627,"date":"2026-09-05T09:28:00","date_gmt":"2026-09-05T09:28:00","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-la-p-valeur-d-une-statistique-f-dans-r\/"},"modified":"2026-09-05T09:28:00","modified_gmt":"2026-09-05T09:28:00","slug":"comment-calculer-la-p-valeur-d-une-statistique-f-dans-r","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-la-p-valeur-d-une-statistique-f-dans-r","title":{"rendered":"Comment calculer la p-valeur d&rsquo;une statistique F dans R"},"content":{"rendered":"<p>Dans l&rsquo;arsenal m\u00e9thodologique du chercheur contemporain, de l&rsquo;\u00e9conom\u00e8tre et du biostatisticien, l&rsquo;analyse de la variance et la mod\u00e9lisation lin\u00e9aire occupent une place pr\u00e9pond\u00e9rante. Au c\u0153ur de ces d\u00e9marches inf\u00e9rentielles r\u00e9side une quantit\u00e9 pivotale : la statistique <em>F<\/em>, con\u00e7ue initialement pour quantifier le rapport de deux variances estim\u00e9es et d\u00e9terminer si la dispersion imputable \u00e0 un mod\u00e8le explicatif d\u00e9passe de mani\u00e8re significative le bruit r\u00e9siduel al\u00e9atoire. Cependant, l&rsquo;obtention de cette valeur empirique ne constitue qu&rsquo;une \u00e9tape interm\u00e9diaire du processus d&rsquo;\u00e9valuation scientifique. Pour transformer cette mesure brute en une d\u00e9cision statistique rigoureuse, il est indispensable de d\u00e9terminer la probabilit\u00e9 d&rsquo;observer une valeur au moins aussi extr\u00eame sous l&rsquo;hypoth\u00e8se de l&rsquo;absence d&rsquo;effet, couramment appel\u00e9e la p-valeur (ou <em>p-value<\/em>).<\/p>\n<p>Le langage de programmation <a rel=\"noopener noreferrer\" href=\"https:\/\/www.r-project.org\/\">R<\/a> s&rsquo;est impos\u00e9 au fil des d\u00e9cennies comme la r\u00e9f\u00e9rence internationale pour l&rsquo;ex\u00e9cution d&rsquo;analyses statistiques avanc\u00e9es, le traitement de donn\u00e9es complexes et la recherche reproductible. Gr\u00e2ce \u00e0 ses fonctions internes hautement optimis\u00e9es bas\u00e9es sur des algorithmes num\u00e9riques de pointe, R permet non seulement d&rsquo;ajuster des mod\u00e8les sophistiqu\u00e9s, mais aussi de calculer avec une pr\u00e9cision analytique absolue les probabilit\u00e9s associ\u00e9es aux distributions th\u00e9oriques. La ma\u00eetrise du calcul de la p-valeur d&rsquo;une statistique <em>F<\/em> au sein de cet environnement ne se limite pas \u00e0 la simple saisie d&rsquo;une commande informatique \u00e9l\u00e9mentaire ; elle implique une compr\u00e9hension intime de la distribution sous-jacente de Fisher-Snedecor, de l&rsquo;alg\u00e8bre des degr\u00e9s de libert\u00e9, ainsi que des particularit\u00e9s algorithmiques propres au traitement des nombres en virgule flottante.<\/p>\n<p>Le pr\u00e9sent trait\u00e9 propose une exploration exhaustive, m\u00e9thodique et didactique des fondements th\u00e9oriques, des impl\u00e9mentations algorithmiques et des applications pratiques du test <em>F<\/em> dans l&rsquo;environnement R. Con\u00e7u \u00e0 l&rsquo;attention des chercheurs, analystes de donn\u00e9es et \u00e9tudiants avanc\u00e9s en m\u00e9thodes quantitatives, ce guide d\u00e9taille chaque \u00e9tape du processus, de l&rsquo;\u00e9quation math\u00e9matique originelle \u00e0 l&rsquo;interpr\u00e9tation des sorties de mod\u00e8les de r\u00e9gression multiple et d&rsquo;analyse de variance (ANOVA), en passant par la visualisation graphique vectorielle et la r\u00e9daction acad\u00e9mique normalis\u00e9e selon les standards de l&rsquo;American Psychological Association (APA).<\/p>\n<h2>1. Introduction aux fondements du test F et de la statistique F<\/h2>\n<h3>1.1 Origine et utilit\u00e9 de la statistique F en analyse quantitative<\/h3>\n<p>L&rsquo;histoire de la statistique moderne est indissociable des contributions pionni\u00e8res de Sir Ronald Aylmer Fisher. Au cours des ann\u00e9es 1920, alors qu&rsquo;il travaillait \u00e0 la station exp\u00e9rimentale agricole de Rothamsted, Fisher s&rsquo;est heurt\u00e9 \u00e0 la n\u00e9cessit\u00e9 d&rsquo;\u00e9valuer simultan\u00e9ment l&rsquo;effet de multiples traitements agronomiques appliqu\u00e9s \u00e0 des parcelles exp\u00e9rimentales, tout en isolant la variabilit\u00e9 naturelle inh\u00e9rente \u00e0 la composition des sols. Les tests usuels de comparaison de deux moyennes, popularis\u00e9s par William Sealy Gosset sous le pseudonyme de Student via la distribution <em>t<\/em>, s&rsquo;av\u00e9raient inadapt\u00e9s pour comparer plus de deux groupes sans gonfler dramatiquement le taux global d&rsquo;erreur de premi\u00e8re esp\u00e8ce (l&rsquo;erreur de type I). Pour surmonter cette limitation, Fisher d\u00e9veloppa l&rsquo;analyse de la variance (ANOVA) et formalisa le rapport de variance, nomm\u00e9 ult\u00e9rieurement \u00ab statistique <em>F<\/em> \u00bb en son honneur par George Snedecor.<\/p>\n<p>Le principe sous-jacent de la statistique <em>F<\/em> repose sur une intuition math\u00e9matique d&rsquo;une rare \u00e9l\u00e9gance : comparer deux estimations ind\u00e9pendantes de la variance d&rsquo;une population sous l&rsquo;hypoth\u00e8se de leur \u00e9galit\u00e9. Plus pr\u00e9cis\u00e9ment, dans un contexte exp\u00e9rimental ou observationnel, la statistique <em>F<\/em> \u00e9value le ratio entre la variance inter-groupe (la variabilit\u00e9 mesur\u00e9e entre les diff\u00e9rents groupes ou conditions, attribuable aux facteurs \u00e9tudi\u00e9s) et la variance intra-groupe (la variabilit\u00e9 mesur\u00e9e au sein de chaque groupe, refl\u00e9tant l&rsquo;erreur de mesure ou le bruit r\u00e9siduel al\u00e9atoire). Lorsque l&rsquo;action des variables explicatives est substantielle, la variance inter-groupe surpasse consid\u00e9rablement la variance intra-groupe, propulsant le ratio <em>F<\/em> vers des valeurs largement sup\u00e9rieures \u00e0 l&rsquo;unit\u00e9.<\/p>\n<p>Dans le champ des sciences comportementales, de la psychologie exp\u00e9rimentale, de la sociologie quantitative et de l&rsquo;\u00e9conom\u00e9trie, cette d\u00e9marche constitue la pierre angulaire de l&rsquo;inf\u00e9rence. Elle autorise le chercheur \u00e0 tester la structure globale de relations complexes entre des covariables continues ou cat\u00e9gorielles et une variable r\u00e9ponse quantitative. Qu&rsquo;il s&rsquo;agisse de mesurer l&rsquo;impact de plusieurs protocoles psychoth\u00e9rapeutiques sur des scores d&rsquo;anxi\u00e9t\u00e9 ou de juger de la pertinence conjointe de multiples variables macro\u00e9conomiques sur l&rsquo;inflation, la statistique <em>F<\/em> fournit le crit\u00e8re quantitatif unifi\u00e9 permettant d&rsquo;\u00e9valuer la pr\u00e9sence d&rsquo;une structure pr\u00e9dictive non fortuite au sein des donn\u00e9es empiriques.<\/p>\n<h3>1.2 Formulation des hypoth\u00e8ses statistiques nulle et alternative<\/h3>\n<p>La rigueur de l&rsquo;inf\u00e9rence statistique repose sur une dichotomie formelle entre deux propositions exclusives : l&rsquo;hypoth\u00e8se nulle, traditionnellement d\u00e9sign\u00e9e par la notation standardis\u00e9e <em>H<\/em><sub>0<\/sub>, et l&rsquo;hypoth\u00e8se alternative, d\u00e9sign\u00e9e par <em>H<\/em><sub>1<\/sub> ou <em>H<\/em><sub>a<\/sub>. Dans le cadre g\u00e9n\u00e9ral de l&rsquo;analyse de variance univari\u00e9e \u00e0 un facteur comprenant <em>k<\/em> groupes distincts, l&rsquo;hypoth\u00e8se nulle postule l&rsquo;\u00e9galit\u00e9 stricte de l&rsquo;ensemble des esp\u00e9rances math\u00e9matiques sous-jacentes. Elle s&rsquo;\u00e9nonce math\u00e9matiquement sous la forme :<\/p>\n<p><em>H<\/em><sub>0<\/sub> : &mu;<sub>1<\/sub> = &mu;<sub>2<\/sub> = &#8230; = &mu;<sub><em>k<\/em><\/sub><\/p>\n<p>Ce postulat implique que l&rsquo;appartenance \u00e0 un groupe particulier n&rsquo;exerce aucune influence syst\u00e9matique sur la variable d\u00e9pendante mesur\u00e9e ; toute divergence observ\u00e9e entre les moyennes \u00e9chantillonnales n&rsquo;est alors imputable qu&rsquo;aux fluctuations d&rsquo;\u00e9chantillonnage al\u00e9atoires. Transpos\u00e9e au contexte de la mod\u00e9lisation lin\u00e9aire multiple comportant <em>p<\/em> variables explicatives, l&rsquo;hypoth\u00e8se nulle stipule que l&rsquo;ensemble des coefficients de r\u00e9gression associ\u00e9s aux pr\u00e9dicteurs sont conjointement nuls :<\/p>\n<p><em>H<\/em><sub>0<\/sub> : &beta;<sub>1<\/sub> = &beta;<sub>2<\/sub> = &#8230; = &beta;<sub><em>p<\/em><\/sub> = 0<\/p>\n<p>Cette assertion \u00e9quivaut \u00e0 affirmer que le mod\u00e8le restreint, ne comportant que la constante (l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine), d\u00e9crit les donn\u00e9es aussi ad\u00e9quatement que le mod\u00e8le complet int\u00e9grant l&rsquo;ensemble des covariables.<\/p>\n<p>\u00c0 l&rsquo;oppos\u00e9, l&rsquo;hypoth\u00e8se alternative <em>H<\/em><sub>1<\/sub> rev\u00eat un caract\u00e8re qualifi\u00e9 d&rsquo;omnibus. Elle ne requiert nullement que toutes les moyennes diff\u00e8rent mutuellement les unes des autres, ni que l&rsquo;int\u00e9gralit\u00e9 des coefficients de r\u00e9gression soient significativement non nuls. Elle pose simplement l&rsquo;existence d&rsquo;au moins une rupture d&rsquo;\u00e9galit\u00e9 ou d&rsquo;au moins une contribution pr\u00e9dictive non nulle :<\/p>\n<p><em>H<\/em><sub>1<\/sub> : Il existe au moins une paire (<em>i<\/em>, <em>j<\/em>) telle que &mu;<sub><em>i<\/em><\/sub> &ne; &mu;<sub><em>j<\/em><\/sub> (pour l&rsquo;ANOVA), ou il existe au moins un indice <em>j<\/em> &isin; {1, &#8230;, <em>p<\/em>} tel que &beta;<sub><em>j<\/em><\/sub> &ne; 0 (pour la r\u00e9gression multiple).<\/p>\n<p>Le rejet statistique de l&rsquo;hypoth\u00e8se nulle entra\u00eene des implications conceptuelles d\u00e9cisives pour le chercheur. Il autorise l&rsquo;abandon formel du mod\u00e8le de pure dispersion al\u00e9atoire au profit d&rsquo;un mod\u00e8le explicatif structur\u00e9. N\u00e9anmoins, en raison de sa nature omnibus, ce rejet n&rsquo;indique en aucun cas la localisation pr\u00e9cise des divergences au sein du dispositif exp\u00e9rimental. Il sert de verrou inf\u00e9rentiel pr\u00e9liminaire : ce n&rsquo;est qu&rsquo;une fois <em>H<\/em><sub>0<\/sub> rejet\u00e9e au seuil global que le praticien est th\u00e9oriquement et m\u00e9thodologiquement fond\u00e9 \u00e0 d\u00e9ployer des proc\u00e9dures de comparaisons multiples planifi\u00e9es ou <em>post-hoc<\/em> (telles que les tests de Tukey, de Bonferroni ou de Scheff\u00e9) ou \u00e0 examiner individuellement la significativit\u00e9 partielle de chaque coefficient &beta; via des tests <em>t<\/em> de Student.<\/p>\n<h3>1.3 Importance de la p-valeur dans l&rsquo;inf\u00e9rence statistique<\/h3>\n<p>La p-valeur occupe une position centrale dans le paradigme fr\u00e9quentiste issu de la fusion historique des approches de Fisher d&rsquo;une part, et de Jerzy Neyman et Egon Pearson d&rsquo;autre part. Math\u00e9matiquement, la p-valeur se d\u00e9finit de mani\u00e8re rigoureuse comme la probabilit\u00e9 conditionnelle d&rsquo;obtenir, sous l&rsquo;hypoth\u00e8se nulle <em>H<\/em><sub>0<\/sub>, une statistique de test au moins aussi extr\u00eame que celle effectivement observ\u00e9e dans l&rsquo;\u00e9chantillon d&rsquo;\u00e9tude. Exprim\u00e9e formellement pour le test <em>F<\/em>, elle correspond \u00e0 l&rsquo;int\u00e9grale sup\u00e9rieure de la densit\u00e9 de probabilit\u00e9 :<\/p>\n<p><em>p<\/em> = <em>P<\/em>(<em>F<\/em> &ge; <em>F<\/em><sub>obs<\/sub> | <em>H<\/em><sub>0<\/sub>)<\/p>\n<p>o\u00f9 <em>F<\/em><sub>obs<\/sub> repr\u00e9sente la valeur num\u00e9rique calcul\u00e9e \u00e0 partir des donn\u00e9es empiriques et <em>F<\/em> la variable al\u00e9atoire th\u00e9orique r\u00e9gie par la distribution de Fisher-Snedecor.<\/p>\n<p>Il importe de clarifier une confusion r\u00e9currente observ\u00e9e dans la litt\u00e9rature scientifique : la p-valeur ne mesure absolument pas la probabilit\u00e9 que l&rsquo;hypoth\u00e8se nulle soit vraie compte tenu des donn\u00e9es recueillies, soit <em>P<\/em>(<em>H<\/em><sub>0<\/sub> | Donn\u00e9es), ce qui constituerait une inversion fallacieuse des probabilit\u00e9s conditionnelles relevant de l&rsquo;inf\u00e9rence bay\u00e9sienne. De m\u00eame, la p-valeur ne quantifie en rien la magnitude ou la port\u00e9e pratique d&rsquo;un ph\u00e9nom\u00e8ne exp\u00e9rimental. Un \u00e9chantillon de tr\u00e8s grande taille dot\u00e9 d&rsquo;une immense puissance statistique peut produire une p-valeur infinit\u00e9simale (par exemple <em>p<\/em> &lt; 0,0001) en pr\u00e9sence d&rsquo;un effet d&rsquo;une trivialit\u00e9 absolue sur le plan th\u00e9orique ou clinique. Inversement, une \u00e9tude men\u00e9e sur un \u00e9chantillon restreint peut mettre en \u00e9vidence un effet substantiel sans atteindre le seuil de significativit\u00e9 conventionnel en raison d&rsquo;une variance r\u00e9siduelle trop \u00e9lev\u00e9e.<\/p>\n<p>L&rsquo;utilisation de seuils critiques pr\u00e9\u00e9tablis, couramment fix\u00e9s \u00e0 &alpha; = 0,05, 0,01 ou 0,001, ob\u00e9it \u00e0 la logique de contr\u00f4le du taux de faux positifs d\u00e9velopp\u00e9e par <a rel=\"noopener noreferrer\" href=\"https:\/\/www.jstor.org\/stable\/91208\">Neyman et Pearson (1933)<\/a>. Si la p-valeur empirique se r\u00e9v\u00e8le strictement inf\u00e9rieure au seuil nominal &alpha;, le chercheur rejette formellement <em>H<\/em><sub>0<\/sub> au b\u00e9n\u00e9fice de <em>H<\/em><sub>1<\/sub>. Toutefois, l&rsquo;usage rigide et ritualis\u00e9 de ces seuils fait l&rsquo;objet de r\u00e9serves substantielles formul\u00e9es par la communaut\u00e9 m\u00e9thodologique internationale, illustr\u00e9es notamment par les d\u00e9clarations successives de l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/www.amstat.org\/\">American Statistical Association<\/a> (ASA) relatives \u00e0 la significativit\u00e9 statistique. La p-valeur doit \u00eatre appr\u00e9hend\u00e9e comme un indicateur continu de la compatibilit\u00e9 des donn\u00e9es empiriques avec le mod\u00e8le th\u00e9orique sous-jacent, et doit syst\u00e9matiquement \u00eatre interpr\u00e9t\u00e9e conjointement avec des m\u00e9triques de taille d&rsquo;effet standardis\u00e9es et leurs intervalles de confiance respectifs.<\/p>\n<h2>2. Compr\u00e9hension th\u00e9orique de la distribution de Fisher-Snedecor<\/h2>\n<h3>2.1 Propri\u00e9t\u00e9s math\u00e9matiques de la loi F<\/h3>\n<p>La distribution de Fisher-Snedecor, not\u00e9e usuellement <em>F<\/em>(<em>d<\/em><sub>1<\/sub>, <em>d<\/em><sub>2<\/sub>), est une distribution de probabilit\u00e9 continue d\u00e9finie exclusivement sur l&rsquo;intervalle des nombres r\u00e9els strictement positifs, avec un support formel born\u00e9 inf\u00e9rieurement par z\u00e9ro et s&rsquo;\u00e9tendant asymptotiquement jusqu&rsquo;\u00e0 l&rsquo;infini : [0, +&infin;[. Elle d\u00e9coule structurellement du comportement al\u00e9atoire de variables suivant des lois du Chi-deux ind\u00e9pendantes. Consid\u00e9rons deux variables al\u00e9atoires ind\u00e9pendantes, <em>U<\/em><sub>1<\/sub> et <em>U<\/em><sub>2<\/sub>, distribu\u00e9es respectivement selon des lois du Chi-deux \u00e0 <em>d<\/em><sub>1<\/sub> et <em>d<\/em><sub>2<\/sub> degr\u00e9s de libert\u00e9 :<\/p>\n<p><em>U<\/em><sub>1<\/sub> ~ &chi;<sup>2<\/sup>(<em>d<\/em><sub>1<\/sub>) et <em>U<\/em><sub>2<\/sub> ~ &chi;<sup>2<\/sup>(<em>d<\/em><sub>2<\/sub>)<\/p>\n<p>La variable al\u00e9atoire r\u00e9sultant du rapport de ces deux entit\u00e9s, chacune pr\u00e9alablement standardis\u00e9e par ses propres degr\u00e9s de libert\u00e9, suit par d\u00e9finition une distribution de Fisher-Snedecor :<\/p>\n<p><em>F<\/em> = (<em>U<\/em><sub>1<\/sub> \/ <em>d<\/em><sub>1<\/sub>) \/ (<em>U<\/em><sub>2<\/sub> \/ <em>d<\/em><sub>2<\/sub>) ~ <em>F<\/em>(<em>d<\/em><sub>1<\/sub>, <em>d<\/em><sub>2<\/sub>)<\/p>\n<p>Sur le plan topologique, la fonction de densit\u00e9 de probabilit\u00e9 de la loi <em>F<\/em> se caract\u00e9rise par une forte asym\u00e9trie positive (\u00e9talement prononc\u00e9 de la distribution vers la droite). La forme exacte de sa courbe d\u00e9pend int\u00e9gralement du couple de param\u00e8tres (<em>d<\/em><sub>1<\/sub>, <em>d<\/em><sub>2<\/sub>). Pour des valeurs modestes de <em>d<\/em><sub>1<\/sub> (en particulier lorsque <em>d<\/em><sub>1<\/sub> = 1 ou <em>d<\/em><sub>1<\/sub> = 2), la densit\u00e9 culmine rapidement pr\u00e8s de l&rsquo;origine avant de d\u00e9cliner de mani\u00e8re exponentielle. Lorsque <em>d<\/em><sub>1<\/sub> exc\u00e8de 2, la densit\u00e9 admet un mode unique situ\u00e9 strictement \u00e0 :<\/p>\n<p>Mode = ( (<em>d<\/em><sub>1<\/sub> &#8211; 2) \/ <em>d<\/em><sub>1<\/sub> ) &times; ( <em>d<\/em><sub>2<\/sub> \/ (<em>d<\/em><sub>2<\/sub> + 2) ) pour <em>d<\/em><sub>1<\/sub> &gt; 2<\/p>\n<p>L&rsquo;esp\u00e9rance math\u00e9matique d&rsquo;une variable distribu\u00e9e selon la loi <em>F<\/em> ne d\u00e9pend de mani\u00e8re remarquable que des degr\u00e9s de libert\u00e9 du d\u00e9nominateur <em>d<\/em><sub>2<\/sub>, \u00e0 condition expresse que ceux-ci soient strictement sup\u00e9rieurs \u00e0 2 :<\/p>\n<p><em>E<\/em>(<em>F<\/em>) = <em>d<\/em><sub>2<\/sub> \/ (<em>d<\/em><sub>2<\/sub> &#8211; 2) pour <em>d<\/em><sub>2<\/sub> &gt; 2<\/p>\n<p>On constate imm\u00e9diatement que lorsque la taille de l&rsquo;\u00e9chantillon s&rsquo;accro\u00eet, faisant tendre <em>d<\/em><sub>2<\/sub> vers l&rsquo;infini, l&rsquo;esp\u00e9rance math\u00e9matique converge asymptotiquement vers l&rsquo;unit\u00e9 exacte (1,00). Cette propri\u00e9t\u00e9 fournit un rep\u00e8re mn\u00e9monique indispensable : sous l&rsquo;hypoth\u00e8se nulle d&rsquo;absence d&rsquo;effet, la statistique <em>F<\/em> fluctue de fa\u00e7on typique autour de 1. Toute d\u00e9viation positive majeure au-del\u00e0 de cette valeur sugg\u00e8re une invalidation de la stochasticit\u00e9 pure.<\/p>\n<h3>2.2 R\u00f4le des degr\u00e9s de libert\u00e9 du num\u00e9rateur et du d\u00e9nominateur<\/h3>\n<p>Les deux param\u00e8tres fondamentaux r\u00e9gissant la morphologie de la distribution <em>F<\/em> sont les degr\u00e9s de libert\u00e9 du num\u00e9rateur (<em>df<\/em><sub>1<\/sub> ou <em>d<\/em><sub>1<\/sub>) et les degr\u00e9s de libert\u00e9 du d\u00e9nominateur (<em>df<\/em><sub>2<\/sub> ou <em>d<\/em><sub>2<\/sub>). Chacun d&rsquo;eux renvoie \u00e0 une dimension g\u00e9om\u00e9trique pr\u00e9cise de l&rsquo;espace vectoriel dans lequel s&rsquo;effectuent les projections des donn\u00e9es empiriques lors de l&rsquo;ajustement du mod\u00e8le statistique.<\/p>\n<p>Les degr\u00e9s de libert\u00e9 du num\u00e9rateur (<em>df<\/em><sub>1<\/sub>), \u00e9galement d\u00e9nomm\u00e9s degr\u00e9s de libert\u00e9 du mod\u00e8le ou du traitement, refl\u00e8tent le nombre de restrictions lin\u00e9aires ind\u00e9pendantes impos\u00e9es par l&rsquo;hypoth\u00e8se nulle, ou de mani\u00e8re \u00e9quivalente, le nombre de param\u00e8tres libres additionnels accord\u00e9s au mod\u00e8le complet par rapport au mod\u00e8le restreint. Dans une analyse de variance factorielle comportant <em>k<\/em> modalit\u00e9s de traitement, <em>df<\/em><sub>1<\/sub> = <em>k<\/em> &#8211; 1. Dans une r\u00e9gression lin\u00e9aire multiple incluant <em>p<\/em> covariables explicatives continues ou indicatrices, <em>df<\/em><sub>1<\/sub> = <em>p<\/em>. Une valeur de <em>df<\/em><sub>1<\/sub> \u00e9lev\u00e9e \u00e9largit la portion initiale de la densit\u00e9 et modifie la pente ascendante de la courbe aux abords de l&rsquo;origine.<\/p>\n<p>Les degr\u00e9s de libert\u00e9 du d\u00e9nominateur (<em>df<\/em><sub>2<\/sub>), syst\u00e9matiquement associ\u00e9s \u00e0 la somme des carr\u00e9s des erreurs r\u00e9siduelles, correspondent \u00e0 la quantit\u00e9 d&rsquo;informations ind\u00e9pendantes disponibles pour estimer la variance de l&rsquo;erreur pure (&sigma;<sup>2<\/sup>). Ils se calculent g\u00e9n\u00e9ralement selon la relation :<\/p>\n<p><em>df<\/em><sub>2<\/sub> = <em>N<\/em> &#8211; <em>p<\/em> &#8211; 1 (pour la r\u00e9gression) ou <em>df<\/em><sub>2<\/sub> = <em>N<\/em> &#8211; <em>k<\/em> (pour l&rsquo;ANOVA \u00e0 un facteur)<\/p>\n<p>o\u00f9 <em>N<\/em> symbolise la taille totale de l&rsquo;\u00e9chantillon d&rsquo;observations. L&rsquo;importance de <em>df<\/em><sub>2<\/sub> sur le comportement inf\u00e9rentiel est capitale : il d\u00e9termine la stabilit\u00e9 num\u00e9rique de l&rsquo;estimation de l&rsquo;erreur r\u00e9siduelle et dicte l&rsquo;\u00e9paisseur de la queue droite de distribution. Pour de faibles valeurs de <em>df<\/em><sub>2<\/sub>, l&rsquo;incertitude sur l&rsquo;estimation de l&rsquo;al\u00e9a r\u00e9siduel est pr\u00e9gnante, ce qui engendre une queue de distribution lourde et \u00e9tir\u00e9e, n\u00e9cessitant des valeurs empiriques de <em>F<\/em> extr\u00eamement \u00e9lev\u00e9es pour franchir les seuils de significativit\u00e9. \u00c0 mesure que <em>df<\/em><sub>2<\/sub> cro\u00eet vers de grands effectifs, la variance d&rsquo;\u00e9chantillonnage de l&rsquo;estimateur de l&rsquo;erreur diminue, la distribution s&rsquo;amincit et se stabilise, diminuant ainsi le seuil critique requis pour rejeter l&rsquo;hypoth\u00e8se nulle.<\/p>\n<h3>2.3 Nature unilat\u00e9rale droite du test F<\/h3>\n<p>Une particularit\u00e9 th\u00e9orique du test <em>F<\/em> r\u00e9side dans la configuration spatiale de sa r\u00e9gion critique : bien que testant souvent des hypoth\u00e8ses non directionnelles (comme l&rsquo;in\u00e9galit\u00e9 globale de multiples moyennes sans ordonnancement pr\u00e9alable), le test <em>F<\/em> s&rsquo;effectue syst\u00e9matiquement sous la forme d&rsquo;un test unilat\u00e9ral sup\u00e9rieur, localis\u00e9 exclusivement dans la queue droite de la distribution th\u00e9orique.<\/p>\n<p>Cette sp\u00e9cificit\u00e9 s&rsquo;explique par la construction quadratique des sommes de carr\u00e9s intervenant au num\u00e9rateur et au d\u00e9nominateur du ratio. La statistique <em>F<\/em> est structur\u00e9e de mani\u00e8re \u00e0 situer la variance expliqu\u00e9e par le mod\u00e8le au num\u00e9rateur et la variance r\u00e9siduelle au d\u00e9nominateur :<\/p>\n<p><em>F<\/em> = Carr\u00e9 Moyen du Mod\u00e8le \/ Carr\u00e9 Moyen R\u00e9siduel = <em>MS<\/em><sub>mod\u00e8le<\/sub> \/ <em>MS<\/em><sub>r\u00e9siduel<\/sub><\/p>\n<p>Lorsque l&rsquo;effet d&rsquo;un traitement s&rsquo;intensifie ou que la capacit\u00e9 pr\u00e9dictive des covariables s&rsquo;accro\u00eet, la dispersion inter-groupes s&rsquo;\u00e9l\u00e8ve de mani\u00e8re quadratique, ce qui augmente m\u00e9caniquement la valeur num\u00e9rique de <em>MS<\/em><sub>mod\u00e8le<\/sub>. L&rsquo;erreur r\u00e9siduelle <em>MS<\/em><sub>r\u00e9siduel<\/sub>, quant \u00e0 elle, demeure un estimateur sans biais de la variance intrins\u00e8que non expliqu\u00e9e. En cons\u00e9quence directe de cette formulation fractionnaire, toute pr\u00e9sence d&rsquo;un effet syst\u00e9matique se traduit obligatoirement par une translation du ratio <em>F<\/em> vers des valeurs largement sup\u00e9rieures \u00e0 1. Les valeurs proches de z\u00e9ro ne traduisent jamais un effet exp\u00e9rimental majeur, mais refl\u00e8tent au contraire un mod\u00e8le dont le pouvoir pr\u00e9dictif est inf\u00e9rieur aux attentes probabilistes moyennes (situation susceptible d&rsquo;indiquer une sur-sp\u00e9cification ou une violation m\u00e9thodologique des conditions d&rsquo;ind\u00e9pendance).<\/p>\n<p>Contrairement aux tests de Student (statistique <em>t<\/em>) ou aux tests normaux centr\u00e9s r\u00e9duits (statistique <em>z<\/em>), o\u00f9 les d\u00e9viations par rapport \u00e0 l&rsquo;hypoth\u00e8se nulle peuvent survenir dans les deux directions polaires de l&rsquo;espace r\u00e9el (justifiant une r\u00e9partition bilat\u00e9rale sym\u00e9trique du risque alpha en &alpha;\/2 de chaque c\u00f4t\u00e9), l&rsquo;espace de rejet de la distribution <em>F<\/em> est donc int\u00e9gralement concentr\u00e9 sur l&rsquo;asymptote droite. L&rsquo;\u00e9valuation informatique de la p-valeur consiste ainsi sans ambigu\u00eft\u00e9 \u00e0 calculer l&rsquo;int\u00e9grale d\u00e9finie de la fonction de densit\u00e9 entre la borne inf\u00e9rieure <em>F<\/em><sub>obs<\/sub> et la borne sup\u00e9rieure infinie +&infin;.<\/p>\n<h2>3. Syntaxe et fonctionnement de la fonction pf() dans l&rsquo;environnement R<\/h2>\n<h3>3.1 D\u00e9composition des arguments de la fonction pf()<\/h3>\n<p>L&rsquo;environnement statistique R int\u00e8gre nativement une famille coh\u00e9rente de quatre fonctions sp\u00e9cialis\u00e9es d\u00e9di\u00e9es \u00e0 la manipulation de la distribution de Fisher-Snedecor : <code>df()<\/code> pour la densit\u00e9 de probabilit\u00e9, <code>pf()<\/code> pour la fonction de r\u00e9partition cumulative (probabilit\u00e9), <code>qf()<\/code> pour la fonction quantile inverse, et <code>rf()<\/code> pour la g\u00e9n\u00e9ration de nombres pseudo-al\u00e9atoires distribu\u00e9s selon la loi <em>F<\/em>. Pour d\u00e9terminer avec une rigueur analytique la p-valeur associ\u00e9e \u00e0 une statistique de test empirique, c&rsquo;est la fonction de distribution cumulative <code>pf()<\/code> qui doit \u00eatre mobilis\u00e9e.<\/p>\n<p>La signature compl\u00e8te de cette fonction se structure formellement selon la syntaxe suivante :<\/p>\n<p><code>pf(q, df1, df2, ncp, lower.tail = TRUE, log.p = FALSE)<\/code><\/p>\n<p>La ma\u00eetrise op\u00e9rationnelle de cette commande implique l&rsquo;assimilation pr\u00e9cise de chacun de ses param\u00e8tres constitutifs :<\/p>\n<ul>\n<li><strong><code>q<\/code><\/strong> : Repr\u00e9sente le quantile ou la valeur num\u00e9rique de la statistique <em>F<\/em> empiriquement observ\u00e9e lors de l&rsquo;exp\u00e9rimentation ou extraite du mod\u00e8le ajust\u00e9. Cet argument vectoriel accepte soit une valeur scalaire isol\u00e9e, soit une structure vectorielle complexe permettant des calculs simultan\u00e9s.<\/li>\n<li><strong><code>df1<\/code><\/strong> : Sp\u00e9cifie les degr\u00e9s de libert\u00e9 associ\u00e9s au num\u00e9rateur du ratio de variance. Ce param\u00e8tre doit \u00eatre un nombre r\u00e9el strictement positif, correspondant g\u00e9n\u00e9ralement au nombre de pr\u00e9dicteurs ou de groupes factoriels diminu\u00e9 d&rsquo;une unit\u00e9.<\/li>\n<li><strong><code>df2<\/code><\/strong> : D\u00e9finit les degr\u00e9s de libert\u00e9 associ\u00e9s au d\u00e9nominateur du ratio de variance (erreur r\u00e9siduelle). Tout comme <code>df1<\/code>, cet argument requiert une valeur strictement positive.<\/li>\n<li><strong><code>ncp<\/code><\/strong> : Correspond au param\u00e8tre de non-centralit\u00e9 (<em>non-centrality parameter<\/em>). Par d\u00e9faut, cet argument est omis ou initialis\u00e9 \u00e0 z\u00e9ro, pla\u00e7ant de facto le calcul sous l&rsquo;\u00e9gide de la distribution <em>F<\/em> centrale standard, conforme aux r\u00e9quisitions du test de l&rsquo;hypoth\u00e8se nulle.<\/li>\n<li><strong><code>lower.tail<\/code><\/strong> : Param\u00e8tre bool\u00e9en logique. Lorsqu&rsquo;il est assign\u00e9 \u00e0 <code>TRUE<\/code> (comportement par d\u00e9faut dans R), la fonction renvoie la probabilit\u00e9 cumul\u00e9e inf\u00e9rieure, c&rsquo;est-\u00e0-dire <em>P<\/em>(<em>F<\/em> &le; <em>q<\/em>). Lorsqu&rsquo;il est configur\u00e9 sur <code>FALSE<\/code>, la fonction calcule directement la probabilit\u00e9 unilat\u00e9rale sup\u00e9rieure, soit <em>P<\/em>(<em>F<\/em> &gt; <em>q<\/em>).<\/li>\n<li><strong><code>log.p<\/code><\/strong> : Param\u00e8tre bool\u00e9en logique. S&rsquo;il est fix\u00e9 \u00e0 <code>TRUE<\/code>, les probabilit\u00e9s r\u00e9sultantes sont d\u00e9livr\u00e9es sous forme de logarithme naturel ln(<em>p<\/em>), un m\u00e9canisme crucial pour pr\u00e9server la pr\u00e9cision num\u00e9rique lors du traitement de p-valeurs infinit\u00e9simales.<\/li>\n<\/ul>\n<h3>3.2 Gestion du param\u00e8tre lower.tail pour obtenir la p-valeur<\/h3>\n<p>L&rsquo;obtention exacte de la p-valeur d\u00e9pend intimement de l&rsquo;assignation conf\u00e9r\u00e9e au param\u00e8tre logique <code>lower.tail<\/code>. Comme expos\u00e9 pr\u00e9alablement, la p-valeur associ\u00e9e au test <em>F<\/em> d\u00e9coule de l&rsquo;\u00e9valuation de la surface sous la courbe situ\u00e9e strictly \u00e0 droite de la valeur observ\u00e9e <em>F<\/em><sub>obs<\/sub>. Or, la d\u00e9finition classique de la fonction de r\u00e9partition cumulative <em>F<\/em>(<em>x<\/em>) adopt\u00e9e par la majorit\u00e9 des logiciels de calcul renvoie syst\u00e9matiquement l&rsquo;aire balay\u00e9e depuis l&rsquo;origine jusqu&rsquo;au point consid\u00e9r\u00e9 :<\/p>\n<p><em>F<\/em>(<em>q<\/em>) = <em>P<\/em>(<em>F<\/em> &le; <em>q<\/em>) = &int;<sub>0<\/sub><sup><em>q<\/em><\/sup> <em>f<\/em>(<em>u<\/em>) <em>d<\/em><em>u<\/em><\/p>\n<p>Si l&rsquo;analyste ex\u00e9cute l&rsquo;instruction <code>pf(q, df1, df2)<\/code> en conservant les valeurs initiales par d\u00e9faut (impliquant <code>lower.tail = TRUE<\/code>), la sortie produite correspondra \u00e0 la probabilit\u00e9 de non-d\u00e9passement, laquelle \u00e9quivaut \u00e0 1 &#8211; p-valeur. Pour isoler la p-valeur v\u00e9ritable, deux approches computationnelles s&rsquo;offrent conceptuellement \u00e0 l&rsquo;utilisateur :<\/p>\n<p>Approche arithm\u00e9tique compl\u00e9mentaire : <code>1 - pf(q, df1, df2)<\/code><br \/>\nApproche directe via param\u00e8tre natif : <code>pf(q, df1, df2, lower.tail = FALSE)<\/code><\/p>\n<p>Bien que ces deux expressions soient math\u00e9matiquement rigoureusement \u00e9quivalentes sur le plan formel, elles diff\u00e8rent radicalement au niveau de leur comportement computationnel au sein de l&rsquo;unit\u00e9 de calcul en virgule flottante. Selon la norme internationale <a rel=\"noopener noreferrer\" href=\"https:\/\/standards.ieee.org\/ieee\/754\/6210\/\">IEEE 754<\/a> r\u00e9gissant l&rsquo;arithm\u00e9tique binaire \u00e0 double pr\u00e9cision, les nombres r\u00e9els sont repr\u00e9sent\u00e9s avec une pr\u00e9cision limit\u00e9e \u00e0 environ 53 bits de mantisse, conf\u00e9rant un maximum de 15 \u00e0 17 chiffres d\u00e9cimaux significatifs. Lorsque la statistique <em>F<\/em> observ\u00e9e atteint des magnitudes tr\u00e8s \u00e9lev\u00e9es (par exemple <em>F<\/em> = 85 avec des degr\u00e9s de libert\u00e9 substantiels), la probabilit\u00e9 cumulative inf\u00e9rieure s&rsquo;approche de l&rsquo;unit\u00e9 avec une proximit\u00e9 telle que l&rsquo;ordinateur l&rsquo;arrondit \u00e0 1.0000000000000000. L&rsquo;ex\u00e9cution subs\u00e9quente de la soustraction arithm\u00e9tique <code>1 - 1<\/code> engendre un ph\u00e9nom\u00e8ne pernicieux d&rsquo;annulation catastrophique, g\u00e9n\u00e9rant en sortie un z\u00e9ro absolu (<code>0<\/code>) et provoquant une perte d&rsquo;information quant \u00e0 l&rsquo;ordre de grandeur r\u00e9el de la significativit\u00e9.<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;impl\u00e9mentation algorithmique interne r\u00e9gissant l&rsquo;instruction <code>lower.tail = FALSE<\/code> effectue une int\u00e9gration directe de la queue sup\u00e9rieure au moyen de d\u00e9veloppements asymptotiques et de fractions continues sophistiqu\u00e9es (notamment bas\u00e9es sur la fonction b\u00eata incompl\u00e8te r\u00e9gularis\u00e9e). Cette optimisation computationnelle permet \u00e0 R de renvoyer des p-valeurs d&rsquo;une finesse extr\u00eame, atteignant des \u00e9chelles de l&rsquo;ordre de 10<sup>-300<\/sup> sans encourir le moindre \u00e9cueil d&rsquo;arrondi pr\u00e9matur\u00e9. Ainsi, l&rsquo;usage d\u00e9lib\u00e9r\u00e9 et syst\u00e9matique de l&rsquo;argument <code>lower.tail = FALSE<\/code> constitue une r\u00e8gle de bonne pratique de programmation en R.<\/p>\n<h3>3.3 Param\u00e8tre ncp et tests de puissance statistique<\/h3>\n<p>L&rsquo;argument <code>ncp<\/code> (abr\u00e9viation de <em>non-centrality parameter<\/em>, ou param\u00e8tre de non-centralit\u00e9 &lambda;) \u00e9largit la distribution de Fisher-Snedecor au-del\u00e0 de son cadre central conventionnel. Alors que la distribution centrale mod\u00e9lise les donn\u00e9es lorsque l&rsquo;hypoth\u00e8se nulle <em>H<\/em><sub>0<\/sub> est rigoureusement exacte (&lambda; = 0), la distribution <em>F<\/em> non centrale d\u00e9crit le comportement de la statistique de test lorsque l&rsquo;hypoth\u00e8se alternative <em>H<\/em><sub>1<\/sub> pr\u00e9vaut dans la population sous-jacente (&lambda; &gt; 0). Le param\u00e8tre &lambda; est fonction directe de la v\u00e9ritable taille de l&rsquo;effet dans la population et de la dimension de l&rsquo;\u00e9chantillon d&rsquo;observation ; il traduit le d\u00e9calage syst\u00e9matique de la distribution vers la droite provoqu\u00e9 par l&rsquo;impact du traitement \u00e9tudi\u00e9.<\/p>\n<p>Dans le protocole usuel d&rsquo;\u00e9valuation de la significativit\u00e9 inf\u00e9rentielle, le chercheur cherche \u00e0 valider ou rejeter <em>H<\/em><sub>0<\/sub>. Il convient donc imp\u00e9rativement d&rsquo;omettre l&rsquo;argument <code>ncp<\/code> ou de le contraindre explicitement \u00e0 z\u00e9ro (<code>ncp = 0<\/code>). Toute assignation erron\u00e9e d&rsquo;une valeur positive \u00e0 ce param\u00e8tre fausserait int\u00e9gralement la p-valeur r\u00e9sultante en superposant une distribution th\u00e9orique d\u00e9centr\u00e9e sur le test de nullit\u00e9.<\/p>\n<p>Toutefois, la manipulation experte de l&rsquo;argument <code>ncp<\/code> s&rsquo;av\u00e8re indispensable lors de la planification exp\u00e9rimentale et de la mod\u00e9lisation de la puissance statistique (not\u00e9e 1 &#8211; &beta;). Si l&rsquo;on souhaite quantifier la probabilit\u00e9 th\u00e9orique de rejeter formellement <em>H<\/em><sub>0<\/sub> sachant qu&rsquo;un effet d&rsquo;une amplitude d\u00e9termin\u00e9e &lambda; existe r\u00e9ellement dans la population, la d\u00e9marche requiert une double \u00e9valuation :<\/p>\n<p>1. D\u00e9termination de la valeur critique <em>F<\/em><sub>crit<\/sub> sur la distribution centrale via la commande :<br \/>\n<code>F_crit &lt;- qf(p = 0.05, df1 = df1, df2 = df2, lower.tail = FALSE)<\/code><br \/>\n2. Calcul de la puissance statistique sur la distribution non centrale via l&rsquo;instruction :<br \/>\n<code>puissance &lt;- pf(q = F_crit, df1 = df1, df2 = df2, ncp = lambda, lower.tail = FALSE)<\/code><\/p>\n<p>Cette flexibilit\u00e9 op\u00e9rationnelle illustre la profondeur de l&rsquo;\u00e9cosyst\u00e8me algorithmique de R, capable de basculer de l&rsquo;\u00e9valuation inf\u00e9rentielle standard aux calculs prospectifs de sensibilit\u00e9 exp\u00e9rimentale par la modification raisonn\u00e9e d&rsquo;un param\u00e8tre math\u00e9matique unitaire.<\/p>\n<h2>4. Calcul direct pas \u00e0 pas de la p-valeur \u00e0 partir d&rsquo;une statistique F connue<\/h2>\n<h3>4.1 Application num\u00e9rique basique avec des param\u00e8tres d\u00e9termin\u00e9s<\/h3>\n<p>Afin d&rsquo;ancrer ces concepts dans une d\u00e9marche empirique imm\u00e9diatement transposable, examinons le cas concret d&rsquo;une \u00e9tude exp\u00e9rimentale o\u00f9 les calculs d&rsquo;analyse de variance pr\u00e9liminaires ont fourni les r\u00e9sultats suivants : une statistique <em>F<\/em> \u00e9gale \u00e0 5,00, calcul\u00e9e avec 3 degr\u00e9s de libert\u00e9 au num\u00e9rateur (indiquant une comparaison entre 4 groupes distincts de sujets, soit <em>k<\/em> &#8211; 1 = 3) et 14 degr\u00e9s de libert\u00e9 au d\u00e9nominateur (correspondant aux r\u00e9sidus d&rsquo;un collectif total de 18 participants).<\/p>\n<p>L&rsquo;ex\u00e9cution de la commande au sein de la console R s&rsquo;effectue par l&rsquo;assignation explicite des variables ou par le passage direct des scalaires au sein de l&rsquo;instruction <code>pf()<\/code> :<\/p>\n<p><code>statistique_F &lt;- 5.00<\/code><br \/>\n<code>degres_liberte_num &lt;- 3<\/code><br \/>\n<code>degres_liberte_den &lt;- 14<\/code><br \/>\n<code>p_valeur &lt;- pf(q = statistique_F, df1 = degres_liberte_num, df2 = degres_liberte_den, lower.tail = FALSE)<\/code><br \/>\n<code>print(p_valeur)<\/code><\/p>\n<p>\u00c0 l&rsquo;ex\u00e9cution, l&rsquo;interpr\u00e9teur de commandes R renvoie la valeur num\u00e9rique d\u00e9cimale brute suivante :<\/p>\n<p><code>[1] 0.01449174<\/code><\/p>\n<p>L&rsquo;analyse de ce r\u00e9sultat brut permet de formuler une conclusion inf\u00e9rentielle transparente. Sous l&rsquo;hypoth\u00e8se nulle stipulant que les 4 groupes proviennent de populations pr\u00e9sentant des moyennes rigoureusement identiques, la probabilit\u00e9 d&rsquo;obtenir par simple fluctuation d&rsquo;\u00e9chantillonnage al\u00e9atoire une statistique <em>F<\/em> sup\u00e9rieure ou \u00e9gale \u00e0 5,00 s&rsquo;\u00e9tablit \u00e0 environ 1,45 %. Si le seuil de d\u00e9cision conventionnel &alpha; a \u00e9t\u00e9 initialis\u00e9 en amont du protocole \u00e0 0,05, la p-valeur obtenue est strictement inf\u00e9rieure \u00e0 ce niveau de risque (0,0145 &lt; 0,05). En cons\u00e9quence, le chercheur rejette formellement l&rsquo;hypoth\u00e8se nulle d&rsquo;homog\u00e9n\u00e9it\u00e9 globale des moyennes et conclut \u00e0 l&rsquo;existence d&rsquo;une diff\u00e9rence statistiquement significative entre au moins deux conditions exp\u00e9rimentales.<\/p>\n<h3>4.2 Formatage et affichage scientifique du r\u00e9sultat num\u00e9rique<\/h3>\n<p>Dans le cadre de l&rsquo;\u00e9laboration de rapports techniques, de th\u00e8ses ou d&rsquo;articles soumis \u00e0 l&rsquo;\u00e9valuation par les pairs, la transcription des sorties d\u00e9cimales brutes g\u00e9n\u00e9r\u00e9es par R sous leur forme non travaill\u00e9e n&rsquo;est pas adapt\u00e9e. Deux imp\u00e9ratifs contradictoires doivent \u00eatre concili\u00e9s : \u00e9viter l&rsquo;accumulation superflue de d\u00e9cimales non significatives tout en pr\u00e9venant l&rsquo;apparition de repr\u00e9sentations exponentielles absconses (telles que <code>1.449174e-02<\/code>) qui nuisent \u00e0 la lisibilit\u00e9.<\/p>\n<p>Le langage R dispose d&rsquo;un jeu de primitives d\u00e9di\u00e9es \u00e0 la mise en forme textuelle des donn\u00e9es chiffr\u00e9es. La fonction <code>round()<\/code> permet d&rsquo;arrondir la p-valeur \u00e0 un nombre d\u00e9termin\u00e9 de chiffres apr\u00e8s la virgule d\u00e9cimale :<\/p>\n<p><code>p_arrondie &lt;- round(p_valeur, digits = 4)<\/code><br \/>\nCe qui g\u00e9n\u00e8re en sortie la cha\u00eene ou le scalaire : <code>0.0145<\/code>.<\/p>\n<p>Cependant, l&rsquo;utilisation exclusive de <code>round()<\/code> pr\u00e9sente une d\u00e9faillance majeure lorsque la p-valeur s&rsquo;av\u00e8re particuli\u00e8rement faible : l&rsquo;application de <code>round(0.000025, digits = 3)<\/code> produira un r\u00e9sultat affichant <code>0<\/code>, ce qui constitue une aberration sur le plan statistique (une p-valeur n&rsquo;\u00e9tant jamais strictement nulle au sein d&rsquo;une distribution asymptotique continue). Pour pallier cette difficult\u00e9, la fonction <code>format()<\/code> s&rsquo;impose par sa polyvalence :<\/p>\n<p><code>format(p_valeur, scientific = FALSE, digits = 4)<\/code><br \/>\nL&rsquo;argument <code>scientific = FALSE<\/code> neutralise formellement l&rsquo;encodage exponentiel, contraignant R \u00e0 maintenir une notation d\u00e9cimale standard jusqu&rsquo;\u00e0 ce que les limites mat\u00e9rielles du terminal d&rsquo;affichage soient atteintes.<\/p>\n<p>Pour automatiser int\u00e9gralement le respect des conventions acad\u00e9miques internationales prescrites par l&rsquo;APA, il est opportun d&rsquo;encapsuler ces r\u00e8gles dans une logique conditionnelle vectoris\u00e9e :<\/p>\n<p><code>p_formatee &lt;- ifelse(p_valeur &lt; 0.001, \"p &lt; .001\", paste0(\"p = \", sub(\"^0\\.\", \".\", sprintf(\"%.3f\", p_valeur))))<\/code><br \/>\nCette formulation avanc\u00e9e effectue deux op\u00e9rations d&rsquo;une grande \u00e9l\u00e9gance typographique : elle seuille automatiquement \u00e0 <code>p &lt; .001<\/code> toute probabilit\u00e9 inf\u00e9rieure au milli\u00e8me, et recourt \u00e0 l&rsquo;expression rationnelle <code>sub(\"^0\\.\", \".\", ...)<\/code> pour retrancher le z\u00e9ro pr\u00e9c\u00e9dant la virgule d\u00e9cimale, conform\u00e9ment aux exigences stylistiques anglophones et acad\u00e9miques qui pr\u00e9conisent d&rsquo;omettre ce chiffre pour toute grandeur math\u00e9matiquement born\u00e9e entre z\u00e9ro et un.<\/p>\n<h3>4.3 V\u00e9rification de la sensibilit\u00e9 du r\u00e9sultat selon les degr\u00e9s de libert\u00e9<\/h3>\n<p>L&rsquo;inf\u00e9rence adoss\u00e9e au test <em>F<\/em> ne se r\u00e9sume pas \u00e0 l&rsquo;amplitude intrins\u00e8que de la statistique calcul\u00e9e ; elle demeure tributaire des degr\u00e9s de libert\u00e9 engag\u00e9s. Pour saisir cette d\u00e9pendance structurelle, il est instructif d&rsquo;observer le comportement de la p-valeur lorsque l&rsquo;on fait varier m\u00e9thodiquement soit la complexit\u00e9 du mod\u00e8le (<em>df<\/em><sub>1<\/sub>), soit l&rsquo;effectif global d&rsquo;\u00e9chantillonnage (<em>df<\/em><sub>2<\/sub>), \u00e0 valeur de <em>F<\/em> rigoureusement constante.<\/p>\n<p>Consid\u00e9rons une s\u00e9quence de simulations ex\u00e9cut\u00e9e sous R, maintenant la statistique observ\u00e9e fix\u00e9e \u00e0 <em>F<\/em> = 3,50, tout en modulant graduellement les degr\u00e9s de libert\u00e9 r\u00e9siduels <em>df<\/em><sub>2<\/sub> (repr\u00e9sentatifs d&rsquo;un \u00e9chantillon s&rsquo;\u00e9tendant de 10 \u00e0 200 observations) pour un num\u00e9rateur fix\u00e9 \u00e0 <em>df<\/em><sub>1<\/sub> = 2 :<\/p>\n<p><code>df2_vecteur &lt;- c(8, 15, 30, 60, 120, 500)<\/code><br \/>\n<code>p_sensibilite &lt;- pf(q = 3.50, df1 = 2, df2 = df2_vecteur, lower.tail = FALSE)<\/code><br \/>\n<code>tableau_sensibilite &lt;- data.frame(df1 = 2, df2 = df2_vecteur, F_valeur = 3.50, p_valeur = round(p_sensibilite, 4))<\/code><br \/>\n<code>print(tableau_sensibilite)<\/code><\/p>\n<p>L&rsquo;affichage de ce tableau de sensibilit\u00e9 illustre une r\u00e9gularit\u00e9 math\u00e9matique fondamentale :<\/p>\n<ul>\n<li>Pour <em>df<\/em><sub>2<\/sub> = 8 : <em>p<\/em> = 0,0809 (l&rsquo;hypoth\u00e8se nulle ne peut \u00eatre rejet\u00e9e au seuil &alpha; = 0,05)<\/li>\n<li>Pour <em>df<\/em><sub>2<\/sub> = 15 : <em>p<\/em> = 0,0563 (la significativit\u00e9 conventionnelle n&rsquo;est toujours pas atteinte)<\/li>\n<li>Pour <em>df<\/em><sub>2<\/sub> = 30 : <em>p<\/em> = 0,0429 (l&rsquo;effet devient statistiquement significatif au seuil de 5 %)<\/li>\n<li>Pour <em>df<\/em><sub>2<\/sub> = 60 : <em>p<\/em> = 0,0365<\/li>\n<li>Pour <em>df<\/em><sub>2<\/sub> = 120 : <em>p<\/em> = 0,0333<\/li>\n<li>Pour <em>df<\/em><sub>2<\/sub> = 500 : <em>p<\/em> = 0,0309<\/li>\n<\/ul>\n<p>Ces donn\u00e9es d\u00e9montrent l&rsquo;impact de l&rsquo;augmentation des degr\u00e9s de libert\u00e9 du d\u00e9nominateur sur la d\u00e9cision statistique : \u00e0 mesure que <em>df<\/em><sub>2<\/sub> augmente, l&rsquo;incertitude sur la variance d&rsquo;\u00e9chantillonnage de l&rsquo;erreur diminue, ce qui affine la p-valeur. Cependant, cette relation n&rsquo;est pas lin\u00e9aire : le gain en pr\u00e9cision statistique est spectaculaire lors des incr\u00e9ments d&rsquo;\u00e9chantillon initiaux (entre 8 et 30), avant de se heurter \u00e0 un plateau asymptotique au-del\u00e0 d&rsquo;une certaine taille critique. Une telle analyse de sensibilit\u00e9 met en exergue l&rsquo;imp\u00e9rieuse n\u00e9cessit\u00e9 de dimensionner convenablement les cohortes d&rsquo;\u00e9tude avant la collecte effective des donn\u00e9es empiriques.<\/p>\n<h2>5. Application pratique : Le test F global dans la r\u00e9gression lin\u00e9aire simple et multiple<\/h2>\n<h3>5.1 Fondement du test F omnibus en r\u00e9gression<\/h3>\n<p>Dans l&rsquo;architecture de la mod\u00e9lisation lin\u00e9aire par les moindres carr\u00e9s ordinaires (MCO), le test <em>F<\/em> global, fr\u00e9quemment d\u00e9sign\u00e9 sous le terme de test <em>F<\/em> omnibus, remplit une mission d&rsquo;\u00e9valuation syst\u00e9mique. Alors que les tests <em>t<\/em> associ\u00e9s aux coefficients de r\u00e9gression individuels examinent l&rsquo;apport marginal de chaque pr\u00e9dicteur suppos\u00e9 ind\u00e9pendant, le test <em>F<\/em> global interroge la viabilit\u00e9 explicative de l&rsquo;ensemble de la combinaison lin\u00e9aire des variables int\u00e9gr\u00e9es au mod\u00e8le math\u00e9matique.<\/p>\n<p>Ce test repose th\u00e9oriquement sur le th\u00e9or\u00e8me fondamental de d\u00e9composition de la variabilit\u00e9, qui postule que la somme des carr\u00e9s totale (<em>SS<\/em><sub>totale<\/sub>), mesurant la dispersion globale des observations autour de la moyenne de la variable r\u00e9ponse <em>Y<\/em>, se scinde rigoureusement en deux composantes additives mutuellement orthogonales :<\/p>\n<p><em>SS<\/em><sub>totale<\/sub> = <em>SS<\/em><sub>r\u00e9gression<\/sub> + <em>SS<\/em><sub>r\u00e9siduelle<\/sub><\/p>\n<p>La premi\u00e8re composante, <em>SS<\/em><sub>r\u00e9gression<\/sub> (ou somme des carr\u00e9s expliqu\u00e9e), mesure la quantit\u00e9 d&rsquo;information captur\u00e9e par le faisceau de covariables. La seconde composante, <em>SS<\/em><sub>r\u00e9siduelle<\/sub>, quantifie la dispersion des \u00e9carts r\u00e9siduels entre les valeurs empiriques observ\u00e9es et les projections ajust\u00e9es par l&rsquo;hyperplan de r\u00e9gression. Pour transformer ces sommes brutes en estimateurs de variance non biais\u00e9s, chacune est rapport\u00e9e \u00e0 ses degr\u00e9s de libert\u00e9 respectifs :<\/p>\n<p><em>MS<\/em><sub>r\u00e9gression<\/sub> = <em>SS<\/em><sub>r\u00e9gression<\/sub> \/ <em>p<\/em> (o\u00f9 <em>p<\/em> est le nombre de covariables hors constante)<br \/>\n<em>MS<\/em><sub>r\u00e9siduelle<\/sub> = <em>SS<\/em><sub>r\u00e9siduelle<\/sub> \/ (<em>n<\/em> &#8211; <em>p<\/em> &#8211; 1)<\/p>\n<p>La statistique <em>F<\/em> omnibus d\u00e9coule alors directement de ce ratio : <em>F<\/em> = <em>MS<\/em><sub>r\u00e9gression<\/sub> \/ <em>MS<\/em><sub>r\u00e9siduelle<\/sub>. Ce ratio peut \u00e9galement s&rsquo;exprimer \u00e9l\u00e9gamment au moyen du coefficient de d\u00e9termination global (<em>R<\/em><sup>2<\/sup>) :<\/p>\n<p><em>F<\/em> = ( <em>R<\/em><sup>2<\/sup> \/ <em>p<\/em> ) \/ ( (1 &#8211; <em>R<\/em><sup>2<\/sup>) \/ (<em>n<\/em> &#8211; <em>p<\/em> &#8211; 1) )<\/p>\n<p>Sous l&rsquo;hypoth\u00e8se nulle <em>H<\/em><sub>0<\/sub> : &beta;<sub>1<\/sub> = &#8230; = &beta;<sub><em>p<\/em><\/sub> = 0, ce quotient suit rigoureusement la distribution <em>F<\/em>(<em>p<\/em>, <em>n<\/em> &#8211; <em>p<\/em> &#8211; 1). L&rsquo;obtention d&rsquo;une p-valeur significative pour ce test atteste formellement que le mod\u00e8le sp\u00e9cifi\u00e9 offre un ajustement sup\u00e9rieur aux donn\u00e9es empiriques par rapport au mod\u00e8le nul se cantonnant \u00e0 pr\u00e9dire chaque individu par la moyenne g\u00e9n\u00e9rale de l&rsquo;\u00e9chantillon.<\/p>\n<h3>5.2 Construction manuelle du jeu de donn\u00e9es sous R<\/h3>\n<p>Afin de contextualiser l&rsquo;application de ces concepts, g\u00e9n\u00e9rons un environnement de simulation contr\u00f4l\u00e9 sous R reproduisant un cadre de recherche en p\u00e9dagogie quantitative et psychom\u00e9trie universitaire. Nous cherchons \u00e0 pr\u00e9dire la note finale obtenue par des \u00e9tudiants \u00e0 un examen standardis\u00e9 de biostatistique (variable d\u00e9pendante continue not\u00e9e <code>performance<\/code>) \u00e0 partir de deux covariables explicatives continues : le volume horaire d&rsquo;\u00e9tude personnelle cumul\u00e9 au cours du semestre (<code>heures_etude<\/code>) et le nombre de sessions de r\u00e9vision m\u00e9thodologique en laboratoire suivies par l&rsquo;\u00e9tudiant (<code>sessions_revision<\/code>).<\/p>\n<p>Pour assurer la reproductibilit\u00e9 des r\u00e9sultats, nous initialisons le g\u00e9n\u00e9rateur de nombres pseudo-al\u00e9atoires au moyen de la primitive <code>set.seed()<\/code>, puis nous g\u00e9n\u00e9rons un \u00e9chantillon r\u00e9aliste comprenant <em>N<\/em> = 45 observations :<\/p>\n<p><code>set.seed(42)<\/code><br \/>\n<code>taille_echantillon &lt;- 45<\/code><br \/>\n<code>heures_etude &lt;- round(rnorm(n = taille_echantillon, mean = 25, sd = 6), 1)<\/code><br \/>\n<code>sessions_revision &lt;- rpois(n = taille_echantillon, lambda = 4)<\/code><br \/>\n<code>erreur_aleatoire &lt;- rnorm(n = taille_echantillon, mean = 0, sd = 4)<\/code><br \/>\n<code>performance &lt;- 40 + (1.2 * heures_etude) + (2.5 * sessions_revision) + erreur_aleatoire<\/code><br \/>\n<code>donnees_academiques &lt;- data.frame(performance, heures_etude, sessions_revision)<\/code><br \/>\n<code>head(donnees_academiques, n = 5)<\/code><\/p>\n<p>L&rsquo;inspection des premi\u00e8res lignes du tableau de donn\u00e9es confirme l&rsquo;int\u00e9grit\u00e9 structurelle des variables synth\u00e9tis\u00e9es. Les distributions marginales pr\u00e9sentent des ordres de grandeur compatibles avec une m\u00e9trique de notes acad\u00e9miques sur cent points. Les deux variables explicatives manifestent une corr\u00e9lation th\u00e9orique positive avec le crit\u00e8re d\u00e9pendant, tout en demeurant alt\u00e9r\u00e9es par un terme de bruit gaussien repr\u00e9sentatif de l&rsquo;erreur de mesure in\u00e9vitable en sciences sociales.<\/p>\n<h3>5.3 Ajustement du mod\u00e8le lin\u00e9aire via la fonction lm()<\/h3>\n<p>Dans l&rsquo;environnement de base de R, la mod\u00e9lisation de relations lin\u00e9aires continues s&rsquo;op\u00e8re au moyen de la fonction <a rel=\"noopener noreferrer\" href=\"https:\/\/stat.ethz.ch\/R-manual\/R-patched\/library\/stats\/html\/lm.html\"><code>lm()<\/code><\/a> (acronyme de <em>linear model<\/em>). Cette routine s&rsquo;appuie sur une d\u00e9composition QR hautement stabilis\u00e9e de la matrice des r\u00e9gresseurs, \u00e9cartant les risques d&rsquo;instabilit\u00e9 num\u00e9rique inh\u00e9rents \u00e0 l&rsquo;inversion directe des matrices de variance-covariance.<\/p>\n<p>La formulation syntaxique s&rsquo;articule autour de l&rsquo;op\u00e9rateur tilde (<code>~<\/code>), s\u00e9parant la variable r\u00e9ponse des termes pr\u00e9dictifs :<\/p>\n<p><code>modele_lineaire &lt;- lm(performance ~ heures_etude + sessions_revision, data = donnees_academiques)<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de cette commande instancie en m\u00e9moire vive un objet complexe de classe <code>lm<\/code>. Cet objet ne se r\u00e9sume pas \u00e0 un simple tableau r\u00e9capitulatif ; il contient l&rsquo;ensemble des \u00e9l\u00e9ments analytiques du mod\u00e8le : les coefficients estim\u00e9s du vecteur &beta;, les r\u00e9sidus ordinaires observ\u00e9s, les valeurs pr\u00e9dites ajust\u00e9es, ainsi que les structures internes de rang matriciel indispensables aux d\u00e9rivations inf\u00e9rentielles ult\u00e9rieures. L&rsquo;appel de l&rsquo;objet direct dans la console (<code>print(modele_lineaire)<\/code>) d\u00e9livre les coefficients d&rsquo;estimation ponctuels, mais n&rsquo;affiche ni les \u00e9carts-types des estimateurs, ni les valeurs de la statistique <em>F<\/em> omnibus, ni sa p-valeur associ\u00e9e, ces m\u00e9triques inf\u00e9rentielles \u00e9tant calcul\u00e9es \u00e0 la vol\u00e9e par des fonctions d&rsquo;inspection sp\u00e9cialis\u00e9es.<\/p>\n<h2>6. Extraction automatis\u00e9e de la statistique F et de sa p-valeur depuis un mod\u00e8le R<\/h2>\n<h3>6.1 Inspection de l&rsquo;objet summary.lm<\/h3>\n<p>L&rsquo;inspection analytique exhaustive d&rsquo;un mod\u00e8le estim\u00e9 par <code>lm()<\/code> s&rsquo;effectue usuellement via l&rsquo;invocation de la m\u00e9thode g\u00e9n\u00e9rique <code>summary()<\/code> :<\/p>\n<p><code>resume_modele &lt;- summary(modele_lineaire)<\/code><br \/>\n<code>print(resume_modele)<\/code><\/p>\n<p>La restitution visuelle g\u00e9n\u00e9r\u00e9e dans la console R structure les r\u00e9sultats en plusieurs blocs hi\u00e9rarchis\u00e9s : distribution des r\u00e9sidus, tableau des coefficients avec leurs statistiques <em>t<\/em> associ\u00e9es, \u00e9cart-type r\u00e9siduel (&sigma;), coefficient de d\u00e9termination (<em>R<\/em><sup>2<\/sup> multiple et ajust\u00e9), pour s&rsquo;achever syst\u00e9matiquement sur une ultime ligne r\u00e9capitulative pr\u00e9sentant l&rsquo;inf\u00e9rence omnibus sous la forme :<\/p>\n<p><code>F-statistic: 72.85 on 2 and 42 DF,  p-value: 2.134e-14<\/code><\/p>\n<p>Si cette pr\u00e9sentation textuelle satisfait pleinement l&rsquo;examen visuel interactif de l&rsquo;analyste humain, elle s&rsquo;av\u00e8re impropre \u00e0 l&rsquo;automatisation dans des scripts programmatiques ou des cha\u00eenes de traitement dynamique (telles que des pipelines R Markdown, Quarto ou Shiny). Pour extraire ces m\u00e9triques sous forme de scalaires r\u00e9utilisables, il est imp\u00e9ratif d&rsquo;interroger la structure sous-jacente de l&rsquo;objet produit par <code>summary()<\/code>. L&rsquo;ex\u00e9cution de la commande de diagnostic <code>names(resume_modele)<\/code> r\u00e9v\u00e8le que cet objet de classe <code>summary.lm<\/code> est une liste nomm\u00e9e constitu\u00e9e de onze composantes distinctes, parmi lesquelles figure l&rsquo;attribut vectoriel sp\u00e9cifique d\u00e9nomm\u00e9 <code>fstatistic<\/code>.<\/p>\n<h3>6.2 Extraction programmatique du vecteur fstatistic<\/h3>\n<p>L&rsquo;attribut <code>fstatistic<\/code> encapsul\u00e9 au sein de l&rsquo;objet <code>summary.lm<\/code> n&rsquo;est pas un scalaire isol\u00e9, mais un vecteur num\u00e9rique nomm\u00e9 \u00e0 trois dimensions, agenc\u00e9 selon l&rsquo;ordre canonique suivant : la valeur empirique de la statistique <em>F<\/em>, les degr\u00e9s de libert\u00e9 du mod\u00e8le (num\u00e9rateur), et les degr\u00e9s de libert\u00e9 des r\u00e9sidus (d\u00e9nominateur).<\/p>\n<p>Pour extraire ce vecteur et diss\u00e9quer programmatiquement ses composantes sans ambigu\u00eft\u00e9, la syntaxe d&rsquo;indexation par double crochet ou par l&rsquo;op\u00e9rateur dollar s&rsquo;av\u00e8re indispensable :<\/p>\n<p><code>vecteur_f &lt;- resume_modele$fstatistic<\/code><br \/>\n<code>print(vecteur_f)<\/code><br \/>\nL&rsquo;affichage console renvoie typiquement :<br \/>\n<code>    value     numdf     dendf <\/code><br \/>\n<code>72.846522  2.000000 42.000000 <\/code><\/p>\n<p>Pour d\u00e9river avec une rigueur absolue la p-valeur \u00e0 partir de ce vecteur dynamique sans recourir \u00e0 aucune saisie manuelle susceptible d&rsquo;introduire des erreurs de transcription, il suffit de combiner l&rsquo;assignation vectorielle avec la fonction <code>pf()<\/code> pr\u00e9c\u00e9demment \u00e9tudi\u00e9e :<\/p>\n<p><code>f_observe &lt;- vecteur_f[\"value\"]<\/code><br \/>\n<code>df_numerateur &lt;- vecteur_f[\"numdf\"]<\/code><br \/>\n<code>df_denominateur &lt;- vecteur_f[\"dendf\"]<\/code><br \/>\n<code>p_valeur_globale &lt;- pf(q = f_observe, df1 = df_numerateur, df2 = df_denominateur, lower.tail = FALSE)<\/code><br \/>\n<code>names(p_valeur_globale) &lt;- NULL<\/code><br \/>\n<code>print(p_valeur_globale)<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de cette s\u00e9quence assigne \u00e0 la variable <code>p_valeur_globale<\/code> le scalaire num\u00e9rique pur <code>2.134371e-14<\/code>. Il convient de souligner que l&rsquo;objet <code>summary.lm<\/code> ne stocke pas directement la p-valeur de l&rsquo;omnibus dans ses tiroirs internes, mais la recalcule dynamiquement \u00e0 l&rsquo;\u00e9cran lors de l&rsquo;ex\u00e9cution de sa fonction d&rsquo;affichage <code>print.summary.lm<\/code>. D\u00e8s lors, la m\u00e9thode d&rsquo;extraction analytique d\u00e9taill\u00e9e ci-dessus constitue l&rsquo;unique voie computationnelle pour r\u00e9cup\u00e9rer cette probabilit\u00e9 sous forme d&rsquo;un nombre \u00e0 virgule flottante au sein d&rsquo;un pipeline statistique automatis\u00e9.<\/p>\n<h3>6.3 Impl\u00e9mentation d&rsquo;une fonction R personnalis\u00e9e d&rsquo;extraction<\/h3>\n<p>Dans une perspective de g\u00e9nie logiciel et de d\u00e9veloppement de code robuste et r\u00e9utilisable, il est fortement d\u00e9conseill\u00e9 de r\u00e9\u00e9crire manuellement ces \u00e9tapes d&rsquo;indexation \u00e0 chaque ajustement de mod\u00e8le. La cr\u00e9ation d&rsquo;une fonction R modulaire, encapsulant la logique d&rsquo;extraction et int\u00e9grant des m\u00e9canismes de gestion des cas limites, repr\u00e9sente la solution m\u00e9thodologique optimale.<\/p>\n<p>Consid\u00e9rons l&rsquo;impl\u00e9mentation de la fonction d&rsquo;ing\u00e9nierie statistique <code>extraire_f_pvaleur()<\/code> suivante :<\/p>\n<p><code>extraire_f_pvaleur &lt;- function(modele) {<\/code><br \/>\n<code>  if (!inherits(modele, \"lm\")) {<\/code><br \/>\n<code>    stop(\"L'argument transmis doit obligatoirement \u00eatre un objet instanci\u00e9 de classe 'lm'.\")<\/code><br \/>\n<code>  }<\/code><br \/>\n<code>  resume &lt;- summary(modele)<\/code><br \/>\n<code>  fstat &lt;- resume$fstatistic<\/code><br \/>\n<code>  if (is.null(fstat)) {<\/code><br \/>\n<code>    warning(\"Le mod\u00e8le sp\u00e9cifi\u00e9 ne comporte aucun degr\u00e9 de libert\u00e9 au num\u00e9rateur (ex: mod\u00e8le vide r\u00e9duit \u00e0 l'intercept).\")<\/code><br \/>\n<code>    return(list(F_stat = NA, df1 = NA, df2 = NA, p_val = NA, texte_apa = NA))<\/code><br \/>\n<code>  }<\/code><br \/>\n<code>  val_f &lt;- as.numeric(fstat[\"value\"])<\/code><br \/>\n<code>  df_1  &lt;- as.numeric(fstat[\"numdf\"])<\/code><br \/>\n<code>  df_2  &lt;- as.numeric(fstat[\"dendf\"])<\/code><br \/>\n<code>  p_val &lt;- pf(q = val_f, df1 = df_1, df2 = df_2, lower.tail = FALSE)<\/code><br \/>\n<code>  chaine_p &lt;- ifelse(p_val &lt; 0.001, \"p &lt; .001\", paste0(\"p = \", sub(\"^0\\.\", \".\", sprintf(\"%.3f\", p_val))))<\/code><br \/>\n<code>  texte_apa &lt;- sprintf(\"F(%d, %d) = %.2f, %s\", df_1, df_2, val_f, chaine_p)<\/code><br \/>\n<code>  return(list(<\/code><br \/>\n<code>    statistique = val_f,<\/code><br \/>\n<code>    degres_liberte_num = df_1,<\/code><br \/>\n<code>    degres_liberte_den = df_2,<\/code><br \/>\n<code>    p_valeur = p_val,<\/code><br \/>\n<code>    format_apa = texte_apa<\/code><br \/>\n<code>  ))<\/code><br \/>\n<code>}<\/code><\/p>\n<p>L&rsquo;\u00e9valuation de notre mod\u00e8le au moyen de cette fonction personnalis\u00e9e s&rsquo;op\u00e8re avec une simplicit\u00e9 syntaxique exemplaire :<\/p>\n<p><code>resultat_extraction &lt;- extraire_f_pvaleur(modele_lineaire)<\/code><br \/>\n<code>print(resultat_extraction$format_apa)<\/code><br \/>\nCe qui restitue instantan\u00e9ment : <code>[1] \"F(2, 42) = 72.85, p &lt; .001\"<\/code>.<\/p>\n<p>Cette routine garantit non seulement l&rsquo;int\u00e9grit\u00e9 des types num\u00e9riques extraits, mais pr\u00e9vient l&rsquo;arr\u00eat inopin\u00e9 des scripts gr\u00e2ce au contr\u00f4le pr\u00e9liminaire de classe via <code>inherits()<\/code> et \u00e0 la gestion pr\u00e9ventive des mod\u00e8les nuls r\u00e9duits \u00e0 l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine (pour lesquels <code>resume$fstatistic<\/code> renvoie la valeur <code>NULL<\/code>).<\/p>\n<h2>7. Calcul de la p-valeur du F dans le cadre de l&rsquo;analyse de variance (ANOVA)<\/h2>\n<h3>7.1 Lien structurel entre r\u00e9gression lin\u00e9aire et ANOVA<\/h3>\n<p>Sur le plan fondamental de l&rsquo;alg\u00e8bre lin\u00e9aire, l&rsquo;analyse de variance (ANOVA) et la r\u00e9gression par les moindres carr\u00e9s ordinaires ne constituent en aucune fa\u00e7on des m\u00e9thodes statistiques divergentes, mais rel\u00e8vent rigoureusement du m\u00eame cadre d&rsquo;unification th\u00e9orique : le Mod\u00e8le Lin\u00e9aire G\u00e9n\u00e9ralis\u00e9 univari\u00e9. Tout protocole d&rsquo;ANOVA factorielle comportant des variables de regroupement nominales peut \u00eatre transcrit sous la forme d&rsquo;une r\u00e9gression multiple \u00e9quivalente par le recours au codage par variables indicatrices (commun\u00e9ment d\u00e9sign\u00e9es sous le vocable de <em>dummy variables<\/em> ou contrastes de traitement).<\/p>\n<p>La distinction entre ces deux approches r\u00e9side exclusivement dans la perspective de lecture privil\u00e9gi\u00e9e par l&rsquo;analyste. Alors que la r\u00e9gression met l&rsquo;accent sur les pentes marginales liant des dimensions sp\u00e9cifiques au crit\u00e8re continu, l&rsquo;ANOVA partitionne la somme des carr\u00e9s totale en blocs discrets de variance imputables aux diff\u00e9rents facteurs exp\u00e9rimentaux et \u00e0 leurs interactions potentielles. Dans une ANOVA \u00e0 un facteur \u00e0 effets fixes comportant <em>k<\/em> groupes, la variabilit\u00e9 inter-groupes refl\u00e8te la somme quadratique des \u00e9carts entre les moyennes de chaque sous-groupe et la grande moyenne de l&rsquo;\u00e9chantillon, pond\u00e9r\u00e9e par l&rsquo;effectif des cellules :<\/p>\n<p><em>SS<\/em><sub>traitements<\/sub> = &sum;<sub><em>i<\/em>=1<\/sub><sup><em>k<\/em><\/sup> <em>n<\/em><sub><em>i<\/em><\/sub> (<em>Y\u0304<\/em><sub><em>i<\/em>&bull;<\/sub> &#8211; <em>Y\u0304<\/em><sub>&bull;&bull;<\/sub>)<sup>2<\/sup><\/p>\n<p>La variance intra-groupe (somme des carr\u00e9s des r\u00e9sidus) quantifie la dispersion des individus autour du centre de gravit\u00e9 de leur propre condition d&rsquo;assignation :<\/p>\n<p><em>SS<\/em><sub>r\u00e9siduelle<\/sub> = &sum;<sub><em>i<\/em>=1<\/sub><sup><em>k<\/em><\/sup> &sum;<sub><em>j<\/em>=1<\/sub><sup><em>n<\/em><sub><em>i<\/em><\/sub><\/sup> (<em>Y<\/em><sub><em>ij<\/em><\/sub> &#8211; <em>Y\u0304<\/em><sub><em>i<\/em>&bull;<\/sub>)<sup>2<\/sup><\/p>\n<p>La statistique <em>F<\/em> issue du ratio des carr\u00e9s moyens associ\u00e9s d\u00e9coule de cette g\u00e9om\u00e9trie vectorielle. En cons\u00e9quence, les p-valeurs calcul\u00e9es dans une table d&rsquo;ANOVA r\u00e9sultent de la m\u00eame loi de Fisher-Snedecor que celles \u00e9mergeant du test global de r\u00e9gression, confirmant l&rsquo;identit\u00e9 math\u00e9matique sous-jacente aux diff\u00e9rents mod\u00e8les.<\/p>\n<h3>7.2 G\u00e9n\u00e9ration de la table d&rsquo;ANOVA avec anova() et aov()<\/h3>\n<p>Au sein de l&rsquo;environnement logiciel R, deux approches fonctionnelles compl\u00e9mentaires permettent de g\u00e9n\u00e9rer les tables de partitionnement de variance : la fonction <code>aov()<\/code>, qui encapsule l&rsquo;ajustement lin\u00e9aire dans une classe d\u00e9di\u00e9e optimis\u00e9e pour les comparaisons de moyennes exp\u00e9rimentales, et la fonction <code>anova()<\/code>, qui d\u00e9rive la table de d\u00e9composition s\u00e9quentielle (dite de Type I dans R) \u00e0 partir d&rsquo;un objet de mod\u00e9lisation pr\u00e9existant.<\/p>\n<p>Illustrons cette application en transformant une variable explicative en facteur discret \u00e0 trois modalit\u00e9s au sein d&rsquo;un nouveau mod\u00e8le :<\/p>\n<p><code>donnees_academiques$groupe_tutorat &lt;- factor(<\/code><br \/>\n<code>  sample(c(\"Aucun\", \"Standard\", \"Intensif\"), size = taille_echantillon, replace = TRUE),<\/code><br \/>\n<code>  levels = c(\"Aucun\", \"Standard\", \"Intensif\")<\/code><br \/>\n<code>)<\/code><br \/>\n<code>modele_aov &lt;- aov(performance ~ groupe_tutorat, data = donnees_academiques)<\/code><br \/>\n<code>table_aov &lt;- summary(modele_aov)<\/code><br \/>\n<code>print(table_aov)<\/code><\/p>\n<p>L&rsquo;affichage structur\u00e9 restitu\u00e9 par l&rsquo;environnement R d\u00e9ploie la table canonique sous la forme suivante :<\/p>\n<p><code>                Df Sum Sq Mean Sq F value Pr(&gt;F)<\/code><br \/>\n<code>groupe_tutorat   2   58.2   29.11   0.245  0.784<\/code><br \/>\n<code>Residuals       42 4991.6  118.85               <\/code><\/p>\n<p>La derni\u00e8re colonne de la table, intitul\u00e9e avec rigueur <code>Pr(&gt;F)<\/code>, mat\u00e9rialise l&rsquo;exacte p-valeur du test unilat\u00e9ral sup\u00e9rieur associ\u00e9 au ratio <em>F<\/em> calcul\u00e9 pour le facteur <code>groupe_tutorat<\/code>. Lorsque le mod\u00e8le est directement estim\u00e9 au moyen de la syntaxe <code>lm()<\/code>, l&rsquo;application de la commande <code>anova(modele_lm)<\/code> d\u00e9livre une table analytique strictement superposable, articulant les sommes de carr\u00e9s, les degr\u00e9s de libert\u00e9 sectoriels, les carr\u00e9s moyens, les valeurs empiriques du ratio <em>F<\/em> et les probabilit\u00e9s inf\u00e9rentielles associ\u00e9es.<\/p>\n<h3>7.3 Calcul manuel et v\u00e9rification des p-valeurs du tableau d&rsquo;ANOVA<\/h3>\n<p>Par souci de validation m\u00e9thodologique et de compr\u00e9hension algorithmique, il est salutaire d&rsquo;isoler les composantes num\u00e9riques de cette table d&rsquo;ANOVA afin de reconstituer manuellement le calcul de la p-valeur via la fonction <code>pf()<\/code>, attestant de l&rsquo;ad\u00e9quation absolue entre la table automatis\u00e9e et les \u00e9quations th\u00e9oriques.<\/p>\n<p>L&rsquo;objet produit par <code>anova(modele)<\/code> h\u00e9rite de la double classe <code>anova<\/code> et <code>data.frame<\/code>. Il peut par cons\u00e9quent \u00eatre interrog\u00e9 au moyen des m\u00e9canismes classiques d&rsquo;indexation matricielle :<\/p>\n<p><code>table_test &lt;- anova(lm(performance ~ groupe_tutorat, data = donnees_academiques))<\/code><br \/>\n<code>f_observe_aov &lt;- table_test[\"groupe_tutorat\", \"F value\"]<\/code><br \/>\n<code>df_facteur &lt;- table_test[\"groupe_tutorat\", \"Df\"]<\/code><br \/>\n<code>df_residus &lt;- table_test[\"Residuals\", \"Df\"]<\/code><br \/>\n<code>p_calculee_directement &lt;- pf(q = f_observe_aov, df1 = df_facteur, df2 = df_residus, lower.tail = FALSE)<\/code><br \/>\n<code>p_fournie_par_table &lt;- table_test[\"groupe_tutorat\", \"Pr(&gt;F)\"]<\/code><br \/>\n<code>print(all.equal(p_calculee_directement, p_fournie_par_table))<\/code><\/p>\n<p>L&rsquo;instruction <code>all.equal()<\/code> renvoie la confirmation logique <code>TRUE<\/code>, attestant que la valeur rapport\u00e9e dans la colonne <code>Pr(&gt;F)<\/code> d\u00e9rive rigoureusement de l&rsquo;\u00e9valuation de la fonction de r\u00e9partition cumulative compl\u00e9mentaire de Fisher-Snedecor. Cette \u00e9quivalence computationnelle d\u00e9montre que la ma\u00eetrise de la fonction <code>pf()<\/code> permet de recalculer, v\u00e9rifier ou ajuster n&rsquo;importe quelle composante d&rsquo;une table d&rsquo;ANOVA complexe, y compris dans le cadre de mod\u00e8les factoriels non orthogonaux ou \u00e0 mesures r\u00e9p\u00e9t\u00e9es.<\/p>\n<h2>8. Comparaison de mod\u00e8les embo\u00eet\u00e9s et test F incr\u00e9mental<\/h2>\n<h3>8.1 Principes du test F de changement de variance expliqu\u00e9e<\/h3>\n<p>Dans de nombreuses applications en recherche quantitative, la question pos\u00e9e ne porte pas tant sur la significativit\u00e9 isol\u00e9e d&rsquo;un mod\u00e8le unitaire, mais sur le gain explicatif incr\u00e9mental apport\u00e9 par l&rsquo;adjonction d&rsquo;un bloc de nouvelles variables th\u00e9oriques \u00e0 un mod\u00e8le de base pr\u00e9existant. Ce protocole rel\u00e8ve de la comparaison formelle de mod\u00e8les embo\u00eet\u00e9s (ou hi\u00e9rarchiques), op\u00e9rationnalis\u00e9e au moyen du test <em>F<\/em> de changement de variance expliqu\u00e9e, \u00e9galement d\u00e9sign\u00e9 sous l&rsquo;appellation de test <em>F<\/em> partiel ou test de restriction lin\u00e9aire.<\/p>\n<p>Consid\u00e9rons deux mod\u00e8les statistiques ajust\u00e9s sur le m\u00eame \u00e9chantillon d&rsquo;observations :<\/p>\n<p>1. Le mod\u00e8le restreint (<em>M<\/em><sub>0<\/sub>), comportant <em>k<\/em> pr\u00e9dicteurs de base :<br \/>\n<em>Y<\/em> = &beta;<sub>0<\/sub> + &beta;<sub>1<\/sub><em>X<\/em><sub>1<\/sub> + &#8230; + &beta;<sub><em>k<\/em><\/sub><em>X<\/em><sub><em>k<\/em><\/sub> + &epsilon;<\/p>\n<p>2. Le mod\u00e8le complet (<em>M<\/em><sub>1<\/sub>), enrichi par l&rsquo;int\u00e9gration de <em>m<\/em> variables explicatives additionnelles :<br \/>\n<em>Y<\/em> = &beta;<sub>0<\/sub> + &beta;<sub>1<\/sub><em>X<\/em><sub>1<\/sub> + &#8230; + &beta;<sub><em>k<\/em><\/sub><em>X<\/em><sub><em>k<\/em><\/sub> + &beta;<sub><em>k<\/em>+1<\/sub><em>X<\/em><sub><em>k<\/em>+1<\/sub> + &#8230; + &beta;<sub><em>k<\/em>+<em>m<\/em><\/sub><em>X<\/em><sub><em>k<\/em>+<em>m<\/em><\/sub> + &epsilon;<\/p>\n<p>L&rsquo;hypoth\u00e8se nulle soumise au test postule que l&rsquo;ensemble des coefficients associ\u00e9s aux variables additionnelles sont nuls dans la population :<\/p>\n<p><em>H<\/em><sub>0<\/sub> : &beta;<sub><em>k<\/em>+1<\/sub> = &beta;<sub><em>k<\/em>+2<\/sub> = &#8230; = &beta;<sub><em>k<\/em>+<em>m<\/em><\/sub> = 0<\/p>\n<p>Le test \u00e9value si la r\u00e9duction de la somme des carr\u00e9s r\u00e9siduels permise par le mod\u00e8le complet (not\u00e9e <em>SS<\/em><sub>res(M0)<\/sub> &#8211; <em>SS<\/em><sub>res(M1)<\/sub>) est statistiquement disproportionn\u00e9e par rapport au co\u00fbt en degr\u00e9s de libert\u00e9 (<em>m<\/em> param\u00e8tres suppl\u00e9mentaires estim\u00e9s). La statistique <em>F<\/em> partielle s&rsquo;exprime formellement ainsi :<\/p>\n<p><em>F<\/em><sub>partiel<\/sub> = [ ( <em>SS<\/em><sub>res(M0)<\/sub> &#8211; <em>SS<\/em><sub>res(M1)<\/sub> ) \/ <em>m<\/em> ] \/ [ <em>SS<\/em><sub>res(M1)<\/sub> \/ (<em>n<\/em> &#8211; <em>k<\/em> &#8211; <em>m<\/em> &#8211; 1) ]<\/p>\n<p>Ce ratio compare le gain quadratique moyen propre au bloc additionnel \u00e0 l&rsquo;estimation la plus raffin\u00e9e de la variance d&rsquo;erreur pure disponible, issue du mod\u00e8le complet non biais\u00e9.<\/p>\n<h3>8.2 Syntaxe de comparaison de deux mod\u00e8les avec anova()<\/h3>\n<p>L&rsquo;impl\u00e9mentation logicielle de cette comparaison d&rsquo;hypoth\u00e8ses embo\u00eet\u00e9es s&rsquo;effectue dans R avec une concision remarquable en transmettant conjointement les deux mod\u00e8les ajust\u00e9s comme arguments ordonn\u00e9s de la fonction g\u00e9n\u00e9rique <code>anova()<\/code>.<\/p>\n<p>D\u00e9veloppons un exemple empirique s&rsquo;appuyant sur les donn\u00e9es simul\u00e9es ant\u00e9rieures, en contrastant un mod\u00e8le restreint ne pr\u00e9disant la note que par les heures d&rsquo;\u00e9tude, et un mod\u00e8le complet int\u00e9grant conjointement les sessions de r\u00e9vision :<\/p>\n<p><code>modele_restreint &lt;- lm(performance ~ heures_etude, data = donnees_academiques)<\/code><br \/>\n<code>modele_complet   &lt;- lm(performance ~ heures_etude + sessions_revision, data = donnees_academiques)<\/code><br \/>\n<code>comparaison_modeles &lt;- anova(modele_restreint, modele_complet)<\/code><br \/>\n<code>print(comparaison_modeles)<\/code><\/p>\n<p>La sortie produite organise m\u00e9thodiquement les donn\u00e9es de la confrontation :<\/p>\n<p><code>Analysis of Variance Table<\/code><br \/>\n<code>Model 1: performance ~ heures_etude<\/code><br \/>\n<code>Model 2: performance ~ heures_etude + sessions_revision<\/code><br \/>\n<code>  Res.Df    RSS Df Sum of Sq      F    Pr(&gt;F)    <\/code><br \/>\n<code>1     43 1948.7                                  <\/code><br \/>\n<code>2     42  698.8  1    1249.9 75.122 1.341e-11 ***<\/code><\/p>\n<p>La seconde ligne du tableau isole la statistique <em>F<\/em> diff\u00e9rentielle (<code>F = 75.122<\/code>), associ\u00e9e \u00e0 1 degr\u00e9 de libert\u00e9 au num\u00e9rateur (adjonction d&rsquo;une variable unique) et 42 degr\u00e9s de libert\u00e9 au d\u00e9nominateur (degr\u00e9s de libert\u00e9 r\u00e9siduels du mod\u00e8le complet). La p-valeur associ\u00e9e (<code>Pr(&gt;F) = 1.341e-11<\/code>) confirme avec une certitude absolue que l&rsquo;adjonction de la variable mesurant les sessions de r\u00e9vision induit une am\u00e9lioration statistiquement hautement significative de la mod\u00e9lisation explicative des scores d&rsquo;examen.<\/p>\n<h3>8.3 Extraction manuelle et confirmation de la p-valeur incr\u00e9mentale<\/h3>\n<p>Pour parachever l&rsquo;appropriation conceptuelle du test <em>F<\/em> partiel, proc\u00e9dons \u00e0 la v\u00e9rification arithm\u00e9tique rigoureuse des grandeurs rapport\u00e9es dans cette table de comparaison embo\u00eet\u00e9e :<\/p>\n<p><code>rss_0 &lt;- deviance(modele_restreint)<\/code><br \/>\n<code>rss_1 &lt;- deviance(modele_complet)<\/code><br \/>\n<code>df_res_0 &lt;- df.residual(modele_restreint)<\/code><br \/>\n<code>df_res_1 &lt;- df.residual(modele_complet)<\/code><br \/>\n<code>delta_rss &lt;- rss_0 - rss_1<\/code><br \/>\n<code>delta_df  &lt;- df_res_0 - df_res_1<\/code><br \/>\n<code>ms_erreur_complet &lt;- rss_1 \/ df_res_1<\/code><br \/>\n<code>f_partiel_calcule &lt;- (delta_rss \/ delta_df) \/ ms_erreur_complet<\/code><br \/>\n<code>p_valeur_partielle &lt;- pf(q = f_partiel_calcule, df1 = delta_df, df2 = df_res_1, lower.tail = FALSE)<\/code><br \/>\n<code>cat(sprintf(\"Statistique F partielle: %.3fn\", f_partiel_calcule))<\/code><br \/>\n<code>cat(sprintf(\"P-valeur incr\u00e9mentale exacte: %.5en\", p_valeur_partielle))<\/code><\/p>\n<p>Ce script restitue des grandeurs strictement identiques \u00e0 la sortie g\u00e9n\u00e9r\u00e9e par <code>anova(modele_restreint, modele_complet)<\/code>. Cette d\u00e9marche d\u00e9montre le fonctionnement de l&rsquo;inf\u00e9rence hi\u00e9rarchique : le test <em>F<\/em> incr\u00e9mental mat\u00e9rialise l&rsquo;exact rapport d&rsquo;orthogonalit\u00e9 entre sous-espaces de projection, permettant de valider empiriquement des hypoth\u00e8ses th\u00e9oriques portant sur des blocs complexes de variables embo\u00eet\u00e9es.<\/p>\n<h2>9. Visualisation graphique de la statistique F et de la zone critique sous R<\/h2>\n<h3>9.1 Trac\u00e9 de la densit\u00e9 th\u00e9orique avec les fonctions graphiques de base<\/h3>\n<p>L&rsquo;appr\u00e9hension intuitive de la p-valeur b\u00e9n\u00e9ficie consid\u00e9rablement d&rsquo;une repr\u00e9sentation graphique spatialisant la densit\u00e9 th\u00e9orique de la loi de Fisher-Snedecor, la valeur empirique observ\u00e9e et l&rsquo;aire g\u00e9om\u00e9trique mat\u00e9rialisant la probabilit\u00e9 critique. Le syst\u00e8me graphique natif de base de R autorise la g\u00e9n\u00e9ration rapide et pr\u00e9cise de tels diagrammes au moyen d&rsquo;outils vectoriels \u00e9l\u00e9mentaires.<\/p>\n<p>Le protocole de visualisation requiert la g\u00e9n\u00e9ration d&rsquo;un maillage r\u00e9gulier de coordonn\u00e9es d&rsquo;abscisse via la fonction <code>seq()<\/code>, l&rsquo;\u00e9valuation de la densit\u00e9 th\u00e9orique associ\u00e9e via <code>df()<\/code>, et le coloriage du polygone d\u00e9limitant la queue critique sup\u00e9rieure au moyen de la primitive <code>polygon()<\/code> :<\/p>\n<p><code>stat_f &lt;- 2.85<\/code><br \/>\n<code>d1 &lt;- 3; d2 &lt;- 25<\/code><br \/>\n<code>abscisses &lt;- seq(from = 0, to = 6, length.out = 1000)<\/code><br \/>\n<code>densites  &lt;- df(x = abscisses, df1 = d1, df2 = d2)<\/code><br \/>\n<code>plot(abscisses, densites, type = \"l\", lwd = 2, col = \"black\",<\/code><br \/>\n<code>     main = \"Distribution de Fisher-Snedecor F(3, 25)\",<\/code><br \/>\n<code>     xlab = \"Valeur de la statistique F\", ylab = \"Densit\u00e9 de probabilit\u00e9\",<\/code><br \/>\n<code>     las = 1, bty = \"l\")<\/code><br \/>\n<code>abscisses_critiques &lt;- seq(from = stat_f, to = 6, length.out = 500)<\/code><br \/>\n<code>densites_critiques  &lt;- df(x = abscisses_critiques, df1 = d1, df2 = d2)<\/code><br \/>\n<code>polygon(x = c(stat_f, abscisses_critiques, 6),<\/code><br \/>\n<code>        y = c(0, densites_critiques, 0),<\/code><br \/>\n<code>        col = rgb(0.8, 0.1, 0.1, alpha = 0.5), border = NA)<\/code><br \/>\n<code>abline(v = stat_f, col = \"darkred\", lwd = 2, lty = 2)<\/code><\/p>\n<p>Ce script de base dresse un profil asym\u00e9trique caract\u00e9ristique. L&rsquo;aire teint\u00e9e de rouge s&rsquo;\u00e9tendant \u00e0 droite de la ligne pointill\u00e9e rouge mat\u00e9rialise directement la p-valeur empirique. Plus cette aire r\u00e9siduelle s&rsquo;amenuise, plus la statistique observ\u00e9e migre vers les zones de r\u00e9jection de la queue de distribution, rendant l&rsquo;hypoth\u00e8se de nullit\u00e9 de plus en plus invraisemblable.<\/p>\n<h3>9.2 Repr\u00e9sentation \u00e9l\u00e9gante de l&rsquo;aire de rejet avec ggplot2<\/h3>\n<p>Pour des standards de publication \u00e9ditoriale sup\u00e9rieure, le package sp\u00e9cialis\u00e9 <a rel=\"noopener noreferrer\" href=\"https:\/\/ggplot2.tidyverse.org\/\"><code>ggplot2<\/code><\/a>, articul\u00e9 autour des principes de la grammaire des graphiques de Leland Wilkinson, offre une expressivit\u00e9 visuelle et une pr\u00e9cision de composition incomparables.<\/p>\n<p>Afin de projeter la distribution <em>F<\/em> sous <code>ggplot2<\/code>, il convient d&rsquo;agencer pr\u00e9alablement les coordonn\u00e9es sous la forme d&rsquo;un tableau rectangulaire (<code>data.frame<\/code>), puis de superposer les couches g\u00e9om\u00e9triques de courbes (<code>geom_line<\/code>), de remplissage d&rsquo;aire conditionnelle (<code>geom_area<\/code>) et de rep\u00e8re vertical (<code>geom_vline<\/code>) :<\/p>\n<p><code>suppressPackageStartupMessages(library(ggplot2))<\/code><br \/>\n<code>df_graphique &lt;- data.frame(F_val = abscisses, Densite = densites)<\/code><br \/>\n<code>f_seuil_alpha &lt;- qf(p = 0.05, df1 = d1, df2 = d2, lower.tail = FALSE)<\/code><br \/>\n<code>p_val_num &lt;- pf(q = stat_f, df1 = d1, df2 = d2, lower.tail = FALSE)<\/code><br \/>\n<code>graphique_f &lt;- ggplot(df_graphique, aes(x = F_val, y = Densite)) +<\/code><br \/>\n<code>  geom_line(color = \"#2C3E50\", linewidth = 1.1) +<\/code><br \/>\n<code>  geom_area(data = subset(df_graphique, F_val &gt;= stat_f),<\/code><br \/>\n<code>            aes(y = Densite), fill = \"#E74C3C\", alpha = 0.5) +<\/code><br \/>\n<code>  geom_vline(xintercept = stat_f, color = \"#C0392B\", linetype = \"dashed\", linewidth = 0.9) +<\/code><br \/>\n<code>  geom_vline(xintercept = f_seuil_alpha, color = \"#27AE60\", linetype = \"dotted\", linewidth = 0.9) +<\/code><br \/>\n<code>  labs(<\/code><br \/>\n<code>    title = \"Densit\u00e9 th\u00e9orique de Fisher-Snedecor et r\u00e9gion critique\",<\/code><br \/>\n<code>    subtitle = sprintf(\"Mod\u00e8le F(%d, %d) | F observ\u00e9 = %.2f (p-valeur = %.4f)\", d1, d2, stat_f, p_val_num),<\/code><br \/>\n<code>    x = \"Statistique F\",<\/code><br \/>\n<code>    y = \"Densit\u00e9 f(F)\"<\/code><br \/>\n<code>  ) +<\/code><br \/>\n<code>  theme_minimal(base_size = 12) +<\/code><br \/>\n<code>  theme(<\/code><br \/>\n<code>    plot.title = element_text(face = \"bold\", size = 14, color = \"#2C3E50\"),<\/code><br \/>\n<code>    plot.subtitle = element_text(size = 11, color = \"#7F8C8D\"),<\/code><br \/>\n<code>    panel.grid.minor = element_blank()<\/code><br \/>\n<code>  )<\/code><br \/>\n<code>print(graphique_f)<\/code><\/p>\n<p>Cette visualisation met simultan\u00e9ment en regard deux rep\u00e8res d\u00e9cisionnels distincts : la ligne pointill\u00e9e verte, qui mat\u00e9rialise le seuil critique th\u00e9orique &alpha; = 0,05 (la valeur au-del\u00e0 de laquelle toute observation justifie un rejet formel de <em>H<\/em><sub>0<\/sub>), et la ligne bris\u00e9e rouge, correspondant \u00e0 la statistique <em>F<\/em> observ\u00e9e empiriquement. La surface ombr\u00e9e en rouge vif visualise directement la p-valeur. Ce type de repr\u00e9sentation permet de clarifier les arbitrages inf\u00e9rentiels lors de communications ou de soutenances acad\u00e9miques.<\/p>\n<h3>9.3 Annotation dynamique des param\u00e8tres et de la p-valeur sur le graphique<\/h3>\n<p>L&rsquo;adjonction d&rsquo;\u00e9l\u00e9ments typographiques explicatifs dynamiquement ancr\u00e9s au sein de l&rsquo;espace de coordonn\u00e9es graphiques enrichit l&rsquo;autonomie documentaire du visuel produit. La fonction <code>annotate()<\/code> de <code>ggplot2<\/code> permet de positionner des cartouches d&rsquo;information textuels calcul\u00e9s en temps r\u00e9el :<\/p>\n<p><code>etiquette_annotation &lt;- sprintf(<\/code><br \/>\n<code>  \"Zone de rejet observ\u00e9enp = %.4fnF(%d, %d) = %.2f\",<\/code><br \/>\n<code>  p_val_num, d1, d2, stat_f<\/code><br \/>\n<code>)<\/code><br \/>\n<code>graphique_annote &lt;- graphique_f +<\/code><br \/>\n<code>  annotate(<\/code><br \/>\n<code>    geom = \"text\",<\/code><br \/>\n<code>    x = stat_f + 1.2,<\/code><br \/>\n<code>    y = max(densites) * 0.5,<\/code><br \/>\n<code>    label = etiquette_annotation,<\/code><br \/>\n<code>    color = \"#C0392B\",<\/code><br \/>\n<code>    fontface = \"bold\",<\/code><br \/>\n<code>    size = 3.8,<\/code><br \/>\n<code>    hjust = 0<\/code><br \/>\n<code>  ) +<\/code><br \/>\n<code>  annotate(<\/code><br \/>\n<code>    geom = \"segment\",<\/code><br \/>\n<code>    x = stat_f + 1.1, xend = stat_f + 0.1,<\/code><br \/>\n<code>    y = max(densites) * 0.45, yend = df(stat_f + 0.1, d1, d2) * 0.8,<\/code><br \/>\n<code>    arrow = arrow(length = unit(0.2, \"cm\")),<\/code><br \/>\n<code>    color = \"#C0392B\"<\/code><br \/>\n<code>  )<\/code><br \/>\n<code>print(graphique_annote)<\/code><\/p>\n<p>Pour exporter le graphique produit vers un support \u00e9ditorial ou une revue scientifique, la fonction <code>ggsave()<\/code> s&rsquo;impose en privil\u00e9giant des formats vectoriels non compress\u00e9s pr\u00e9servant l&rsquo;int\u00e9grit\u00e9 de la r\u00e9solution typographique \u00e0 l&rsquo;impression :<\/p>\n<p><code>ggsave(filename = \"distribution_f_modele.pdf\", plot = graphique_annote, width = 8, height = 5, device = \"pdf\")<\/code><\/p>\n<h2>10. Erreurs courantes, pi\u00e8ges m\u00e9thodologiques et diagnostics<\/h2>\n<h3>10.1 Inversion de l&rsquo;argument lower.tail et interpr\u00e9tation erron\u00e9e<\/h3>\n<p>L&rsquo;omission de l&rsquo;argument <code>lower.tail = FALSE<\/code> lors de l&rsquo;appel de la commande <code>pf()<\/code> constitue l&rsquo;\u00e9cueil le plus fr\u00e9quent chez les praticiens manipulant l&rsquo;environnement R. Par coh\u00e9rence d&rsquo;architecture logicielle avec les distributions sym\u00e9triques centr\u00e9es ou les calculs de quantiles unilat\u00e9raux inf\u00e9rieurs, les concepteurs de R ont initialis\u00e9 le comportement par d\u00e9faut de la fonction sur <code>lower.tail = TRUE<\/code>.<\/p>\n<p>Examinons les cons\u00e9quences de cette m\u00e9prise informatique sur un test <em>F<\/em> empirique mesur\u00e9 \u00e0 <em>F<\/em>(2, 40) = 6,20 :<\/p>\n<p><code>erreur_inversion &lt;- pf(q = 6.20, df1 = 2, df2 = 40)<\/code><br \/>\n<code>print(erreur_inversion)<\/code><br \/>\nL&rsquo;interpr\u00e9teur renvoie la valeur : <code>0.9955746<\/code>.<\/p>\n<p>Un utilisateur non averti interpr\u00e9tant directement cette sortie num\u00e9rique brute pourrait en d\u00e9duire que le test affiche une probabilit\u00e9 avoisinant 0,996, et conclure de mani\u00e8re fallacieuse \u00e0 une absence totale de significativit\u00e9 statistique en consid\u00e9rant que <em>p<\/em> &gt;&gt; 0,05. Or, la p-valeur v\u00e9ritable correspond au compl\u00e9ment sup\u00e9rieur :<\/p>\n<p><code>vraie_pvaleur &lt;- pf(q = 6.20, df1 = 2, df2 = 40, lower.tail = FALSE)<\/code><br \/>\n<code>print(vraie_pvaleur)<\/code><br \/>\nLa valeur exacte s&rsquo;\u00e9tablit \u00e0 <code>0.004425<\/code>, soit un effet hautement significatif au seuil &alpha; = 0,01.<\/p>\n<p>Pour conjurer syst\u00e9matiquement cette source d&rsquo;erreur, une heuristique de contr\u00f4le de vraisemblance s&rsquo;impose : d\u00e8s lors qu&rsquo;une statistique <em>F<\/em> observ\u00e9e exc\u00e8de largement le seuil empirique de 2 ou 3 avec des degr\u00e9s de libert\u00e9 raisonnables, la p-valeur associ\u00e9e doit obligatoirement \u00eatre inf\u00e9rieure \u00e0 0,10. L&rsquo;obtention d&rsquo;une valeur sup\u00e9rieure \u00e0 0,50 pour un ratio <em>F<\/em> \u00e9lev\u00e9 constitue la signature univoque d&rsquo;une inversion de queue de distribution dans la syntaxe computationnelle.<\/p>\n<h3>10.2 Permutation des degr\u00e9s de libert\u00e9 du num\u00e9rateur et du d\u00e9nominateur<\/h3>\n<p>Une alt\u00e9ration sym\u00e9trique tout aussi d\u00e9vastatrice pour la validit\u00e9 de l&rsquo;inf\u00e9rence consiste \u00e0 permuter l&rsquo;attribution des degr\u00e9s de libert\u00e9 du num\u00e9rateur (<code>df1<\/code>) et du d\u00e9nominateur (<code>df2<\/code>). Contrairement \u00e0 d&rsquo;autres distributions sym\u00e9triques, la loi de Fisher-Snedecor est fondamentalement non commutative : <em>F<\/em>(<em>d<\/em><sub>1<\/sub>, <em>d<\/em><sub>2<\/sub>) ne pr\u00e9sente pas les m\u00eames propri\u00e9t\u00e9s que <em>F<\/em>(<em>d<\/em><sub>2<\/sub>, <em>d<\/em><sub>1<\/sub>).<\/p>\n<p>Consid\u00e9rons une situation exp\u00e9rimentale caract\u00e9ris\u00e9e par <em>F<\/em><sub>obs<\/sub> = 3,10, avec <em>df<\/em><sub>1<\/sub> = 3 et <em>df<\/em><sub>2<\/sub> = 80 (comparaison de 4 groupes sur 84 sujets) :<\/p>\n<p><code>p_correcte &lt;- pf(q = 3.10, df1 = 3, df2 = 80, lower.tail = FALSE)<\/code><br \/>\n<code>p_erronee  &lt;- pf(q = 3.10, df1 = 80, df2 = 3, lower.tail = FALSE)<\/code><br \/>\n<code>cat(sprintf(\"P-valeur param\u00e9tr\u00e9e correctement : %.4fn\", p_correcte))<\/code><br \/>\n<code>cat(sprintf(\"P-valeur avec permutation de df : %.4fn\", p_erronee))<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de ces instructions met en lumi\u00e8re une distorsion majeure :<\/p>\n<ul>\n<li>Calcul exact : <em>p<\/em> = 0,0313 (rejet de l&rsquo;hypoth\u00e8se nulle au seuil &alpha; = 0,05)<\/li>\n<li>Calcul biais\u00e9 par inversion : <em>p<\/em> = 0,1789 (\u00e9chec du rejet de l&rsquo;hypoth\u00e8se nulle)<\/li>\n<\/ul>\n<p>La r\u00e8gle mn\u00e9monique pour pr\u00e9venir cet \u00e9garement consiste \u00e0 se rem\u00e9morer la structure alg\u00e9brique du test : le mod\u00e8le explicatif pr\u00e9c\u00e8de l&rsquo;al\u00e9a r\u00e9siduel. Le premier degr\u00e9 de libert\u00e9 (<code>df1<\/code>) quantifie l&rsquo;effort d&rsquo;ajustement du mod\u00e8le (les pr\u00e9dicteurs, au num\u00e9rateur), tandis que le second degr\u00e9 de libert\u00e9 (<code>df2<\/code>) quantifie la taille de la cohorte d&rsquo;\u00e9valuation diminu\u00e9e des param\u00e8tres estim\u00e9s (les r\u00e9sidus, au d\u00e9nominateur). D\u00e8s lors, dans l&rsquo;immense majorit\u00e9 des applications empiriques, <code>df2<\/code> surpasse substantiellement <code>df1<\/code>.<\/p>\n<h3>10.3 Violation des postulats d&rsquo;application et fausses conclusions<\/h3>\n<p>La validit\u00e9 nominale de la p-valeur d\u00e9livr\u00e9e par la fonction <code>pf()<\/code> repose sur le strict respect d&rsquo;un ensemble de postulats probabilistes r\u00e9gissant le mod\u00e8le lin\u00e9aire gaussien :<\/p>\n<ul>\n<li><strong>Normalit\u00e9 des erreurs<\/strong> : Les r\u00e9sidus de population doivent suivre une distribution normale multivari\u00e9e centr\u00e9e sur z\u00e9ro.<\/li>\n<li><strong>Homosc\u00e9dasticit\u00e9<\/strong> : La variance conditionnelle des erreurs doit demeurer constante sur l&rsquo;ensemble du domaine de pr\u00e9diction (&sigma;<sub><em>i<\/em><\/sub><sup>2<\/sup> = &sigma;<sup>2<\/sup> pour tout <em>i<\/em>).<\/li>\n<li><strong>Ind\u00e9pendance stochastique<\/strong> : Les observations doivent \u00eatre rigoureusement ind\u00e9pendantes les unes des autres (absence d&rsquo;autocorr\u00e9lation temporelle, spatiale ou de structure de grappes).<\/li>\n<\/ul>\n<p>Lorsque ces postulats sont transgress\u00e9s, la statistique empirique calcul\u00e9e ne suit plus la distribution th\u00e9orique <em>F<\/em>(<em>d<\/em><sub>1<\/sub>, <em>d<\/em><sub>2<\/sub>). En particulier, la pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 s\u00e9v\u00e8re (variance des erreurs instable selon les groupes ou les valeurs ajust\u00e9es) gonfle artificiellement le taux de faux positifs, produisant des p-valeurs nominalement significatives (ex: <em>p<\/em> = 0,02) alors que le risque r\u00e9el de type I d\u00e9passe parfois 15 ou 20 %.<\/p>\n<p>En pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 confirm\u00e9e (par exemple via le test de Breusch-Pagan, <code>bptest()<\/code> du package <code>lmtest<\/code>), le recours direct \u00e0 la statistique <em>F<\/em> ordinaire est prohib\u00e9. Deux alternatives m\u00e9thodologiques s&rsquo;imposent :<\/p>\n<p>1. Dans le cadre de l&rsquo;ANOVA, mobiliser la proc\u00e9dure de Welch qui corrige simultan\u00e9ment les degr\u00e9s de libert\u00e9 r\u00e9siduels au moyen de la primitive native de R :<br \/>\n<code>oneway.test(performance ~ groupe_tutorat, data = donnees_academiques, var.equal = FALSE)<\/code><\/p>\n<p>2. Dans le cadre de la r\u00e9gression lin\u00e9aire multiple, utiliser un estimateur de matrice de variance-covariance robuste aux h\u00e9t\u00e9rosc\u00e9dasticit\u00e9s (dite correction de White ou HC3), accessible via le package sp\u00e9cialis\u00e9 <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=sandwich\"><code>sandwich<\/code><\/a> coupl\u00e9 au test de Wald robuste instrument\u00e9 par <code>car::linearHypothesis()<\/code> :<br \/>\n<code>library(sandwich)<\/code><br \/>\n<code>library(car)<\/code><br \/>\n<code>linearHypothesis(modele_lineaire, c(\"heures_etude = 0\", \"sessions_revision = 0\"), vcov = vcovHC(modele_lineaire, type = \"HC3\"))<\/code><\/p>\n<p>Cette d\u00e9marche ajuste directement la statistique de test et ses degr\u00e9s de libert\u00e9 r\u00e9siduels, restaurant la validit\u00e9 probabiliste de la p-valeur r\u00e9sultante.<\/p>\n<h2>11. R\u00e9daction et rapport acad\u00e9mique selon les normes APA<\/h2>\n<h3>11.1 Directives de pr\u00e9sentation textuelle des r\u00e9sultats du test F<\/h3>\n<p>La transcription des analyses statistiques dans les manuscrits scientifiques r\u00e9gis par le style de l&rsquo;American Psychological Association (notamment la 7e \u00e9dition du <em>Publication Manual of the APA<\/em>) ob\u00e9it \u00e0 un ensemble de r\u00e8gles typographiques standardis\u00e9es. L&rsquo;objectif est d&rsquo;assurer une lisibilit\u00e9 unifi\u00e9e et de fournir \u00e0 la communaut\u00e9 des pairs l&rsquo;ensemble des \u00e9l\u00e9ments informationnels n\u00e9cessaires \u00e0 l&rsquo;interpr\u00e9tation critique et \u00e0 la m\u00e9ta-analyse des donn\u00e9es publi\u00e9es.<\/p>\n<p>Pour tout test <em>F<\/em> rapport\u00e9 au sein du corps du texte, la formulation formelle impose le respect des crit\u00e8res suivants :<\/p>\n<ul>\n<li>La lettre majuscule <strong><em>F<\/em><\/strong> doit syst\u00e9matiquement \u00eatre compos\u00e9e en caract\u00e8res italiques.<\/li>\n<li>Les degr\u00e9s de libert\u00e9 du num\u00e9rateur et du d\u00e9nominateur doivent \u00eatre explicitement mentionn\u00e9s entre parenth\u00e8ses imm\u00e9diatement apr\u00e8s la lettre du test, s\u00e9par\u00e9s par une virgule suivie d&rsquo;un espace typographique, sans espace entre l&rsquo;ouverture de la parenth\u00e8se et le premier chiffre : <em>F<\/em>(<em>df<\/em><sub>1<\/sub>, <em>df<\/em><sub>2<\/sub>).<\/li>\n<li>La valeur empirique calcul\u00e9e pour la statistique <em>F<\/em> doit \u00eatre s\u00e9par\u00e9e des parenth\u00e8ses par une espace, suivie du signe \u00e9gal (=), et arrondie avec une rigueur constante \u00e0 deux d\u00e9cimales significatives.<\/li>\n<li>L&rsquo;\u00e9nonc\u00e9 de la statistique de test doit \u00eatre imm\u00e9diatement compl\u00e9t\u00e9 par la p-valeur exacte et par une m\u00e9trique standardis\u00e9e de la taille de l&rsquo;effet mesur\u00e9 (telle que le <em>R<\/em><sup>2<\/sup> ajust\u00e9 pour les r\u00e9gressions ou l&rsquo;\u00eata-carr\u00e9 partiel &eta;<sub>p<\/sub><sup>2<\/sup> pour les analyses de variance).<\/li>\n<\/ul>\n<h3>11.2 Formatage standardis\u00e9 de la p-valeur selon l&rsquo;APA 7e \u00e9dition<\/h3>\n<p>Les recommandations de l&rsquo;APA 7e \u00e9dition introduisent des contraintes \u00e9ditoriales pr\u00e9cises concernant la restitution des probabilit\u00e9s inf\u00e9rentielles, souvent n\u00e9glig\u00e9es lors du transfert des r\u00e9sultats depuis la console R vers le traitement de texte :<\/p>\n<p>1. <strong>Omission du z\u00e9ro initial<\/strong> : Dans la mesure o\u00f9 une p-valeur est math\u00e9matiquement born\u00e9e au sein du segment [0, 1] et ne peut th\u00e9oriquement jamais atteindre ni d\u00e9passer l&rsquo;unit\u00e9 absolue, les normes APA prohibent formellement l&rsquo;adjonction d&rsquo;un z\u00e9ro avant le point d\u00e9cimal. Il convient ainsi d&rsquo;\u00e9crire <em>p<\/em> = .024 et non pas <em>p<\/em> = 0.024.<\/p>\n<p>2. <strong>Pr\u00e9cision d\u00e9cimale<\/strong> : La p-valeur doit \u00eatre rapport\u00e9e avec une pr\u00e9cision fix\u00e9e \u00e0 trois chiffres d\u00e9cimaux significatifs, comme <em>p<\/em> = .042 ou <em>p<\/em> = .318.<\/p>\n<p>3. <strong>Gestion des valeurs infinit\u00e9simales<\/strong> : Lorsque la console R renvoie une p-valeur extr\u00eamement basse (affich\u00e9e sous une notation scientifique comme <code>2.13e-14<\/code>), il est formellement prohib\u00e9 d&rsquo;\u00e9crire <em>p<\/em> = .000, cette formulation constituant un non-sens inf\u00e9rentiel. Les directives imposent dans ce cas de tronquer la notation sous la forme d&rsquo;une in\u00e9galit\u00e9 stricte : <em>p<\/em> &lt; .001.<\/p>\n<p>4. <strong>Interdiction des mentions vagues<\/strong> : Les d\u00e9signations obsol\u00e8tes telles que \u00ab <em>p<\/em> = NS \u00bb (pour <em>non significatif<\/em>) doivent \u00eatre totalement proscrites au b\u00e9n\u00e9fice de la transcription exacte de la probabilit\u00e9 mesur\u00e9e (par exemple <em>p<\/em> = .412), garantissant ainsi une totale transparence m\u00e9trologique.<\/p>\n<h3>11.3 Exemples textuels complets contextualis\u00e9s en sciences du comportement<\/h3>\n<p>Afin d&rsquo;illustrer la mise en pratique de ces directives typographiques et conceptuelles, examinons deux canevas de r\u00e9daction repr\u00e9sentatifs des standards \u00e9ditoriaux en psychologie exp\u00e9rimentale et en neurosciences computationnelles.<\/p>\n<p><strong>Exemple 1 : Rapport d&rsquo;une r\u00e9gression multiple (\u00e9valuation pr\u00e9dictive de l&rsquo;anxi\u00e9t\u00e9)<\/strong><br \/>\n<em>\u00ab Une r\u00e9gression lin\u00e9aire multiple a \u00e9t\u00e9 conduite afin de d\u00e9terminer dans quelle mesure le niveau de stress per\u00e7u et la qualit\u00e9 du sommeil pr\u00e9disent les scores de symptomatologie anxieuse chez de jeunes adultes. Le mod\u00e8le explicatif global s&rsquo;est r\u00e9v\u00e9l\u00e9 hautement significatif, F(2, 87) = 24.38, p &lt; .001, attestant que la combinaison de ces pr\u00e9dicteurs rend compte d&rsquo;une part substantielle de la variance observ\u00e9e, R<sup>2<\/sup> = .359, R<sup>2<\/sup><sub>ajust\u00e9<\/sub> = .344. L&rsquo;examen des coefficients individuels indique que le stress per\u00e7u contribue positivement et significativement au mod\u00e8le (&beta; = 0.48, t = 5.12, p &lt; .001), tandis qu&rsquo;un sommeil de meilleure qualit\u00e9 exerce un effet protecteur ind\u00e9pendant (&beta; = -0.26, t = -2.78, p = .007). \u00bb<\/em><\/p>\n<p><strong>Exemple 2 : Rapport d&rsquo;une analyse de variance factorielle univari\u00e9e (intervention cognitive)<\/strong><br \/>\n<em>\u00ab Une analyse de variance univari\u00e9e \u00e0 un facteur a \u00e9t\u00e9 r\u00e9alis\u00e9e pour \u00e9valuer l&rsquo;efficacit\u00e9 de trois protocoles d&rsquo;entra\u00eenement de la m\u00e9moire de travail (Contr\u00f4le passif, Entra\u00eenement standardis\u00e9, Entra\u00eenement adaptatif) sur le temps de r\u00e9action \u00e0 une t\u00e2che d&rsquo;attention soutenue. Les r\u00e9sultats mettent en \u00e9vidence un effet principal significatif de la condition exp\u00e9rimentale, F(2, 57) = 6.84, p = .002, &eta;<sub>p<\/sub><sup>2<\/sup> = .193. Les comparaisons post-hoc ajust\u00e9es au moyen de la m\u00e9thode de Tukey indiquent que les participants assign\u00e9s au protocole adaptatif pr\u00e9sentent des temps de r\u00e9action significativement plus rapides que ceux du groupe contr\u00f4le (d = 0.82, p = .002), alors que la divergence entre le protocole standardis\u00e9 et le groupe contr\u00f4le n&rsquo;atteint pas le seuil de significativit\u00e9 nominal (d = 0.31, p = .341). \u00bb<\/em><\/p>\n<h2>12. Synth\u00e8se proc\u00e9durale et bonnes pratiques de programmation sous R<\/h2>\n<h3>12.1 Arbre de d\u00e9cision proc\u00e9dural pour le calcul de la p-valeur<\/h3>\n<p>L&rsquo;ex\u00e9cution optimale du calcul d&rsquo;une p-valeur adoss\u00e9e \u00e0 une statistique <em>F<\/em> sous R varie selon le contexte analytique rencontr\u00e9 par le praticien. Afin de structurer le choix de la m\u00e9thode computationnelle, l&rsquo;arbre de d\u00e9cision proc\u00e9dural suivant permet de guider l&rsquo;analyste vers la routine technique appropri\u00e9e :<\/p>\n<p><strong>Cas 1 : La statistique F et les degr\u00e9s de libert\u00e9 sont d\u00e9j\u00e0 connus isol\u00e9ment<\/strong> (ex: extraction d&rsquo;un article scientifique, donn\u00e9es tabul\u00e9es manuelles, sorties d&rsquo;un logiciel tiers).<br \/>\n&rarr; <em>Action recommand\u00e9e<\/em> : Recourir directement \u00e0 la fonction de distribution univari\u00e9e :<br \/>\n<code>pf(q = valeur_F, df1 = df_num, df2 = df_den, lower.tail = FALSE)<\/code><\/p>\n<p><strong>Cas 2 : Les donn\u00e9es brutes sont mod\u00e9lis\u00e9es via une r\u00e9gression lin\u00e9aire dans un script d&rsquo;analyse reproductible<\/strong>.<br \/>\n&rarr; <em>Action recommand\u00e9e<\/em> : Ajuster le mod\u00e8le avec <code>lm()<\/code>, extraire le vecteur des composantes inf\u00e9rentielles via <code>summary(modele)$fstatistic<\/code>, puis \u00e9valuer dynamiquement la probabilit\u00e9 unilat\u00e9rale sup\u00e9rieure via <code>pf()<\/code> pour \u00e9viter tout codage en dur de constantes chiffr\u00e9es dans le code source.<\/p>\n<p><strong>Cas 3 : L&rsquo;analyse rel\u00e8ve d&rsquo;une comparaison entre groupes factoriels exp\u00e9rimentaux (ANOVA)<\/strong>.<br \/>\n&rarr; <em>Action recommand\u00e9e<\/em> : Utiliser <code>anova(modele_lm)<\/code> ou <code>summary(aov(formule))<\/code> pour obtenir directement la table de d\u00e9composition int\u00e9grant la colonne standardis\u00e9e <code>Pr(&gt;F)<\/code>.<\/p>\n<p><strong>Cas 4 : L&rsquo;analyste cherche \u00e0 \u00e9valuer l&rsquo;apport incr\u00e9mental de nouvelles variables au sein de mod\u00e8les hi\u00e9rarchiques<\/strong>.<br \/>\n&rarr; <em>Action recommand\u00e9e<\/em> : Ajuster s\u00e9quentiellement le mod\u00e8le restreint et le mod\u00e8le enrichi, puis ex\u00e9cuter la comparaison globale via <code>anova(modele_restreint, modele_enrichi)<\/code>, qui isole automatiquement la statistique <em>F<\/em> incr\u00e9mentale et sa p-valeur partielle.<\/p>\n<h3>12.2 Cr\u00e9ation d&rsquo;un package ou script utilitaire r\u00e9utilisable<\/h3>\n<p>Pour assurer la p\u00e9rennit\u00e9 et la portabilit\u00e9 des d\u00e9veloppements analytiques au sein d&rsquo;un laboratoire de recherche ou d&rsquo;une direction technique, il est judicieux de regrouper ces routines de traitement statistique dans un script utilitaire autonome (par exemple nomm\u00e9 <code>outils_inference.R<\/code>) destin\u00e9 \u00e0 \u00eatre sourc\u00e9 \u00e0 l&rsquo;initialisation de chaque projet.<\/p>\n<p>Consid\u00e9rons l&rsquo;impl\u00e9mentation de la fonction modulaire avanc\u00e9e <code>calculer_f_pvaleur_robuste()<\/code>, dot\u00e9e de gardes d&rsquo;assertion garantissant la validit\u00e9 math\u00e9matique des param\u00e8tres d&rsquo;entr\u00e9e :<\/p>\n<p><code>calculer_f_pvaleur_robuste &lt;- function(f_stat, df1, df2, formater = TRUE) {<\/code><br \/>\n<code>  # Assertions d\u00e9fensives d'int\u00e9grit\u00e9 num\u00e9rique<\/code><br \/>\n<code>  stopifnot(is.numeric(f_stat), length(f_stat) == 1, !is.na(f_stat))<\/code><br \/>\n<code>  stopifnot(is.numeric(df1), length(df1) == 1, df1 &gt; 0)<\/code><br \/>\n<code>  stopifnot(is.numeric(df2), length(df2) == 1, df2 &gt; 0)<\/code><br \/>\n<code>  if (f_stat &lt; 0) {<\/code><br \/>\n<code>    warning(\"Une statistique F ne peut \u00eatre n\u00e9gative. Valeur contrainte \u00e0 0.\")<\/code><br \/>\n<code>    f_stat &lt;- 0<\/code><br \/>\n<code>  }<\/code><br \/>\n<code>  # Calcul haute pr\u00e9cision sur queue droite<\/code><br \/>\n<code>  p_brute &lt;- pf(q = f_stat, df1 = df1, df2 = df2, lower.tail = FALSE)<\/code><br \/>\n<code>  if (!formater) {<\/code><br \/>\n<code>    return(p_brute)<\/code><br \/>\n<code>  }<\/code><br \/>\n<code>  chaine_p &lt;- if (p_brute &lt; 0.001) {<\/code><br \/>\n<code>    \"p &lt; .001\"<\/code><br \/>\n<code>  } else {<\/code><br \/>\n<code>    paste0(\"p = \", sub(\"^0\\.\", \".\", sprintf(\"%.3f\", p_brute)))<\/code><br \/>\n<code>  }<\/code><br \/>\n<code>  return(list(<\/code><br \/>\n<code>    statistique_f = f_stat,<\/code><br \/>\n<code>    df_num = df1,<\/code><br \/>\n<code>    df_den = df2,<\/code><br \/>\n<code>    p_valeur_exacte = p_brute,<\/code><br \/>\n<code>    rapport_apa = sprintf(\"F(%s, %s) = %.2f, %s\", df1, df2, f_stat, chaine_p)<\/code><br \/>\n<code>  ))<\/code><br \/>\n<code>}<\/code><\/p>\n<p>L&rsquo;int\u00e9gration de tests unitaires formels adoss\u00e9s au package <a rel=\"noopener noreferrer\" href=\"https:\/\/testthat.r-lib.org\/\"><code>testthat<\/code><\/a> permet de certifier que les fonctions utilitaires r\u00e9sistent aux valeurs aberrantes, aux types inattendus et aux structures matricielles incompatibles, \u00e9levant ainsi les scripts d&rsquo;analyse de donn\u00e9es aux standards de l&rsquo;ing\u00e9nierie logicielle contemporaine.<\/p>\n<h3>12.3 Perspectives et transition vers l&rsquo;inf\u00e9rence bay\u00e9sienne<\/h3>\n<p>Bien que le calcul de la p-valeur d&rsquo;une statistique <em>F<\/em> demeure l&rsquo;une des d\u00e9marches les plus fr\u00e9quemment sollicit\u00e9es dans la recherche empirique contemporaine, il convient d&rsquo;en reconna\u00eetre lucidement les limites inh\u00e9rentes \u00e0 la m\u00e9canique fr\u00e9quentiste. Une p-valeur infinit\u00e9simale n&rsquo;autorise en aucun cas \u00e0 quantifier la plausibilit\u00e9 de l&rsquo;hypoth\u00e8se de recherche au regard des th\u00e9ories concurrentes ; elle indique simplement que les donn\u00e9es observ\u00e9es s&rsquo;\u00e9cartent sensiblement des pr\u00e9dictions d&rsquo;un mod\u00e8le nul d&rsquo;invariance souvent irr\u00e9aliste.<\/p>\n<p>Pour d\u00e9passer les controverses associ\u00e9es \u00e0 la dichotomisation rigide des d\u00e9cisions inf\u00e9rentielles (le fameux \u00e9cueil du seuil <em>p<\/em> &lt; .05), un nombre croissant de m\u00e9thodologistes pr\u00e9conisent de compl\u00e9ter ou de substituer le test <em>F<\/em> classique par des d\u00e9marches issues de l&rsquo;inf\u00e9rence bay\u00e9sienne. Au c\u0153ur de ce paradigme, le facteur de Bayes (not\u00e9 <em>BF<\/em><sub>10<\/sub>) quantifie directement la force probante relative accord\u00e9e par les donn\u00e9es empiriques \u00e0 l&rsquo;hypoth\u00e8se alternative par rapport \u00e0 l&rsquo;hypoth\u00e8se nulle :<\/p>\n<p><em>BF<\/em><sub>10<\/sub> = <em>P<\/em>(Donn\u00e9es | <em>H<\/em><sub>1<\/sub>) \/ <em>P<\/em>(Donn\u00e9es | <em>H<\/em><sub>0<\/sub>)<\/p>\n<p>L&rsquo;environnement R offre des biblioth\u00e8ques logicielles de pointe d\u00e9di\u00e9es \u00e0 cette transition computationnelle, au premier rang desquelles figurent le package <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=BayesFactor\"><code>BayesFactor<\/code><\/a> de Richard Morey et Jeffrey Rouder, ainsi que le package de mod\u00e9lisation bay\u00e9sienne par cha\u00eenes de Markov de Monte-Carlo (MCMC) <a rel=\"noopener noreferrer\" href=\"https:\/\/paul-buerkner.github.io\/brms\/\"><code>brms<\/code><\/a>. L&rsquo;\u00e9valuation bay\u00e9sienne d&rsquo;un mod\u00e8le lin\u00e9aire s&rsquo;op\u00e8re par exemple au moyen de l&rsquo;instruction :<\/p>\n<p><code>suppressPackageStartupMessages(library(BayesFactor))<\/code><br \/>\n<code>bf_modele &lt;- regressionBF(performance ~ heures_etude + sessions_revision, data = donnees_academiques)<\/code><br \/>\n<code>print(bf_modele)<\/code><\/p>\n<p>La sortie g\u00e9n\u00e9r\u00e9e exprime directement le facteur de Bayes attestant que les donn\u00e9es observ\u00e9es sont, par exemple, plusieurs millions de fois plus probables sous le mod\u00e8le incluant les covariables explicatives que sous le mod\u00e8le nul de pure stochasticit\u00e9. La ma\u00eetrise op\u00e9rationnelle conjointe de la p-valeur fr\u00e9quentiste issue de <code>pf()<\/code> et des facteurs de Bayes permet ainsi au chercheur contemporain de d\u00e9ployer une vision int\u00e9gr\u00e9e, nuanc\u00e9e et transparente de la preuve statistique.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7e \u00e9d.). 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>Fisher, R. A. (1925). <em>Statistical methods for research workers<\/em>. Oliver and Boyd.<\/li>\n<li>Fox, J., &amp; Weisberg, S. (2019). <em>An R companion to applied regression<\/em> (3e \u00e9d.). SAGE Publications. <a rel=\"noopener noreferrer\" href=\"https:\/\/socialsciences.mcmaster.ca\/jfox\/Books\/Companion\/\">https:\/\/socialsciences.mcmaster.ca\/jfox\/Books\/Companion\/<\/a><\/li>\n<li>Institute of Electrical and Electronics Engineers. (2019). <em>IEEE standard for floating-point arithmetic<\/em> (IEEE Std 754-2019). IEEE. <a rel=\"noopener noreferrer\" href=\"https:\/\/standards.ieee.org\/ieee\/754\/6210\/\">https:\/\/standards.ieee.org\/ieee\/754\/6210\/<\/a><\/li>\n<li>Neyman, J., &amp; Pearson, E. S. (1933). On the problem of the most efficient tests of statistical hypotheses. <em>Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character<\/em>, 231(694-706), 289-337. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1098\/rsta.1933.0009\">https:\/\/doi.org\/10.1098\/rsta.1933.0009<\/a><\/li>\n<li>R Core Team. (2023). <em>R: A language and environment for statistical computing<\/em>. R Foundation for Statistical Computing, Vienna, Austria. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.r-project.org\/\">https:\/\/www.r-project.org\/<\/a><\/li>\n<li>Snedecor, G. W., &amp; Cochran, W. G. (1989). <em>Statistical methods<\/em> (8e \u00e9d.). Iowa State University Press.<\/li>\n<li>Wasserstein, R. L., &amp; Lazar, N. A. (2016). The ASA statement on p-values: Context, process, and purpose. <em>The American Statistician<\/em>, 70(2), 129-133. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00031305.2016.1154108\">https:\/\/doi.org\/10.1080\/00031305.2016.1154108<\/a><\/li>\n<li>Wickham, H. (2016). <em>ggplot2: Elegant graphics for data analysis<\/em>. Springer-Verlag New York. <a rel=\"noopener noreferrer\" href=\"https:\/\/ggplot2.tidyverse.org\/\">https:\/\/ggplot2.tidyverse.org\/<\/a><\/li>\n<li>Zeileis, A. (2004). Econometric computing with HC and HAC covariance matrix estimators. <em>Journal of Statistical Software<\/em>, 11(10), 1-17. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.18637\/jss.v011.i10\">https:\/\/doi.org\/10.18637\/jss.v011.i10<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour calculer et interpr\u00e9ter la p-valeur d&rsquo;une statistique F dans R via la fonction pf() et les mod\u00e8les de r\u00e9gression lin\u00e9aire.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2193,1838],"tags":[2224,3406,2658,1910,3405],"class_list":["post-1627","psych_stat","type-psych_stat","status-publish","hentry","category-programmation-r","category-statistiques","tag-anova","tag-p-valeur","tag-r","tag-regression-lineaire","tag-statistique-f"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1627","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=1627"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1627\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1627"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1627"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1627"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}