{"id":8251,"date":"2026-09-17T05:01:40","date_gmt":"2026-09-17T05:01:40","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-erreur-absolue-moyenne-dans-r\/"},"modified":"2026-09-17T05:01:40","modified_gmt":"2026-09-17T05:01:40","slug":"comment-calculer-erreur-absolue-moyenne-dans-r","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-erreur-absolue-moyenne-dans-r","title":{"rendered":"Comment calculer l&rsquo;erreur absolue moyenne dans R"},"content":{"rendered":"<p>L&rsquo;\u00e9valuation de la performance pr\u00e9dictive d&rsquo;un mod\u00e8le statistique constitue le pivot fondamental sur lequel repose toute inf\u00e9rence empirique rigoureuse. Au sein de la mod\u00e9lisation statistique contemporaine et de l&rsquo;apprentissage automatique, la quantification de la discordance entre les observations empiriques et les estimations mod\u00e9lis\u00e9es n\u00e9cessite le recours \u00e0 des m\u00e9triques d&rsquo;erreur robustes, interpr\u00e9tables et th\u00e9oriquement fond\u00e9es. Parmi l&rsquo;arsenal d&rsquo;indicateurs d\u00e9velopp\u00e9s par la th\u00e9orie de l&rsquo;estimation, l&rsquo;erreur absolue moyenne (couramment d\u00e9sign\u00e9e sous l&rsquo;acronyme anglo-saxon <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Mean_absolute_error\">MAE pour <em>Mean Absolute Error<\/em><\/a>) occupe une position privil\u00e9gi\u00e9e. Elle offre une mesure lin\u00e9aire et directe de la magnitude moyenne des \u00e9carts r\u00e9siduels, sans alt\u00e9ration quadratique pr\u00e9alable de l&rsquo;\u00e9chelle des donn\u00e9es observ\u00e9es.<\/p>\n<p>Dans l&rsquo;\u00e9cosyst\u00e8me du langage de programmation statistique <a rel=\"noopener noreferrer\" href=\"https:\/\/www.r-project.org\/\">R<\/a>, plateforme de r\u00e9f\u00e9rence pour les biostatisticiens, \u00e9conom\u00e8tres, psychom\u00e9triciens et scientifiques des donn\u00e9es, le calcul de la MAE s&rsquo;articule autour de multiples paradigmes d&rsquo;impl\u00e9mentation. Qu&rsquo;il s&rsquo;agisse de routines vectorielles primitives en R base, de structures modulaires issues de biblioth\u00e8ques sp\u00e9cialis\u00e9es telles que <code>Metrics<\/code>, ou encore d&rsquo;architectures unifi\u00e9es int\u00e9gr\u00e9es dans des m\u00e9ta-frameworks de pointe comme <code>caret<\/code> et <code>tidymodels<\/code>, R conf\u00e8re aux analystes une flexibilit\u00e9 computationnelle in\u00e9gal\u00e9e. N\u00e9anmoins, l&rsquo;obtention d&rsquo;une valeur scalaire brute ne saurait suffire : la ma\u00eetrise de la MAE requiert une compr\u00e9hension profonde de ses fondements math\u00e9matiques, de son comportement face aux distributions asym\u00e9triques et aux valeurs atypiques, ainsi que de ses r\u00e9percussions pratiques sur l&rsquo;arbitrage entre la m\u00e9diane et la moyenne conditionnelle.<\/p>\n<p>Ce guide exhaustif a pour vocation d&rsquo;analyser en profondeur les m\u00e9canismes th\u00e9oriques, fonctionnels et algorithmiques de l&rsquo;erreur absolue moyenne dans R. \u00c0 travers un continuum didactique reliant la d\u00e9composition analytique des pertes de type L1 jusqu&rsquo;\u00e0 l&rsquo;orchestration de protocoles de validation crois\u00e9e complexes et de r\u00e9\u00e9chantillonnage non param\u00e9trique par <em>bootstrap<\/em>, cet article explore les dimensions conceptuelles, computationnelles et appliqu\u00e9es de la m\u00e9trique. Que vous mod\u00e9lisiez des temps de r\u00e9action en chronom\u00e9trie cognitive, des scores composites d&rsquo;inventaires psychologiques ou des trajectoires longitudinales, l&rsquo;analyse suivante fournit les fondations statistiques et les m\u00e9thodologies concr\u00e8tes indispensables pour une mise en \u0153uvre sans faille de la MAE.<\/p>\n<h2>1. Fondements th\u00e9oriques et d\u00e9finition conceptuelle de l&rsquo;erreur absolue moyenne<\/h2>\n<h3>1.1 D\u00e9finition \u00e9pist\u00e9mologique et r\u00f4le de la MAE dans l&rsquo;\u00e9valuation des mod\u00e8les<\/h3>\n<p>L&rsquo;\u00e9mergence des mesures de pr\u00e9cision dans la th\u00e9orie de l&rsquo;estimation statistique r\u00e9pond \u00e0 une probl\u00e9matique \u00e9pist\u00e9mologique fondamentale : comment quantifier, de mani\u00e8re univoque et op\u00e9rationnelle, la fid\u00e9lit\u00e9 d&rsquo;une abstraction formelle vis-\u00e0-vis d&rsquo;une r\u00e9alit\u00e9 empirique fluctuante ? Historiquement formalis\u00e9e dans les travaux pionniers d&rsquo;astronomes et de math\u00e9maticiens tels qu&rsquo;Adrien-Marie Legendre et Carl Friedrich Gauss lors de la r\u00e9solution du probl\u00e8me de l&rsquo;ajustement des orbites c\u00e9lestes, l&rsquo;\u00e9valuation de l&rsquo;erreur s&rsquo;est scind\u00e9e entre deux grandes approches : la p\u00e9nalisation quadratique des \u00e9carts et la quantification lin\u00e9aire des d\u00e9viations absolues.<\/p>\n<p>L&rsquo;erreur absolue moyenne s&rsquo;ancre dans le concept d&rsquo;\u00e9cart absolu univari\u00e9 appliqu\u00e9 aux variables continues. Contrairement aux approches visant la variance r\u00e9siduelle pure, la MAE traite l&rsquo;erreur comme une distance physique directe entre deux points dans un espace unidimensionnel. Il est \u00e0 cet \u00e9gard imp\u00e9ratif d&rsquo;\u00e9tablir une diff\u00e9renciation fondamentale entre le biais directionnel et la magnitude de l&rsquo;erreur. Une erreur arithm\u00e9tique simple non absolue permet des compensations mutuelles : une sous-estimation syst\u00e9matique de dix unit\u00e9s sur un groupe de sujets peut \u00eatre math\u00e9matiquement neutralis\u00e9e par une surestimation sym\u00e9trique de dix unit\u00e9s sur un second groupe, g\u00e9n\u00e9rant une moyenne des erreurs nulle qui dissimulerait une invalidit\u00e9 pr\u00e9dictive compl\u00e8te. La MAE pr\u00e9vient ce ph\u00e9nom\u00e8ne d&rsquo;annulation en appliquant l&rsquo;op\u00e9rateur valeur absolue \u00e0 chaque \u00e9cart r\u00e9siduel avant l&rsquo;agr\u00e9gation.<\/p>\n<p>Dans l&rsquo;\u00e9valuation empirique des mod\u00e8les psychologiques et pr\u00e9dictifs, cette caract\u00e9ristique conf\u00e8re \u00e0 la m\u00e9trique une pertinence heuristique majeure. Qu&rsquo;il s&rsquo;agisse de pr\u00e9dire le niveau d&rsquo;\u00e9puisement professionnel sur une \u00e9chelle psychom\u00e9trique standardis\u00e9e ou de quantifier la r\u00e9activit\u00e9 cardiovasculaire sous induction de stress exp\u00e9rimental, la MAE fournit au chercheur la magnitude moyenne attendue de l&rsquo;erreur individuelle, pr\u00e9servant l&rsquo;int\u00e9grit\u00e9 du jugement clinique et statistique sans induire de distorsions computationnelles arbitraires.<\/p>\n<h3>1.2 Propri\u00e9t\u00e9s math\u00e9matiques et comportement face aux distributions<\/h3>\n<p>Sur le plan math\u00e9matique, l&rsquo;erreur absolue moyenne se distingue par la lin\u00e9arit\u00e9 parfaite de son \u00e9chelle de p\u00e9nalit\u00e9s. Alors que les fonctions de co\u00fbt quadratiques imposent une p\u00e9nalisation proportionnelle au carr\u00e9 de l&rsquo;\u00e9cart \u2014 accordant un poids exponentiellement sup\u00e9rieur aux r\u00e9sidus de grande amplitude \u2014, la MAE conserve une fonction de perte homog\u00e8ne de degr\u00e9 un. Chaque incr\u00e9ment unitaire d&rsquo;erreur r\u00e9siduelle, qu&rsquo;il s&rsquo;agisse du passage d&rsquo;un r\u00e9sidu de 1 \u00e0 2 ou de 100 \u00e0 101, g\u00e9n\u00e8re rigoureusement la m\u00eame augmentation marginale de la p\u00e9nalit\u00e9 globale.<\/p>\n<p>Cette lin\u00e9arit\u00e9 engendre toutefois une particularit\u00e9 d&rsquo;ordre analytique : la non-diff\u00e9rentiabilit\u00e9 de la fonction de perte au point z\u00e9ro. D&rsquo;un point de vue de l&rsquo;optimisation math\u00e9matique, la d\u00e9riv\u00e9e de la fonction valeur absolue n&rsquo;est pas continue en z\u00e9ro, basculant brutalement d&rsquo;une pente de -1 \u00e0 une pente de +1. Cette singularit\u00e9 interdit l&rsquo;usage imm\u00e9diat d&rsquo;algorithmes de descente de gradient standard reposant sur la diff\u00e9rentiabilit\u00e9 continue de classe C1, imposant le recours aux sous-gradients, \u00e0 la programmation lin\u00e9aire ou \u00e0 des approximations liss\u00e9es telles que la <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Huber_loss\">fonction de perte de Huber<\/a>.<\/p>\n<p>D&rsquo;un point de vue g\u00e9om\u00e9trique, la MAE mat\u00e9rialise la distance de Manhattan, correspondant \u00e0 la norme L1 dans les espaces vectoriels fonctionnels norm\u00e9s. Cette relation g\u00e9om\u00e9trique explique son insensibilit\u00e9 relative aux valeurs aberrantes (<em>outliers<\/em>) par rapport aux normes euclidiennes (norme L2). Les distributions empiriques pr\u00e9sentant des queues lourdes (leptokurtiques), courantes dans les sciences du comportement et les donn\u00e9es financi\u00e8res, n&rsquo;exercent pas sur la MAE l&rsquo;attraction gravitationnelle disproportionn\u00e9e qu&rsquo;elles exercent sur la moyenne des erreurs quadratiques, assurant ainsi une stabilit\u00e9 structurelle remarquable du diagnostic statistique.<\/p>\n<h3>1.3 Interpr\u00e9tabilit\u00e9 dimensionnelle et s\u00e9mantique des r\u00e9sultats<\/h3>\n<p>Le principal avantage op\u00e9rationnel de l&rsquo;erreur absolue moyenne r\u00e9side dans la conservation stricte des unit\u00e9s de mesure de la variable d\u00e9pendante d&rsquo;origine. Si un mod\u00e8le statistique est \u00e9labor\u00e9 pour pr\u00e9dire des scores d&rsquo;anxi\u00e9t\u00e9 mesur\u00e9s sur l&rsquo;\u00e9chelle STAI (<em>State-Trait Anxiety Inventory<\/em>), oscillant entre 20 et 80 points, une MAE de 4,25 s&rsquo;interpr\u00e8te directement comme une marge d&rsquo;erreur moyenne de 4,25 points sur ladite \u00e9chelle. Cette propri\u00e9t\u00e9 dimensionnelle \u00e9vite le recours \u00e0 des op\u00e9rations mentales de transformation inverse, indispensables lors de l&rsquo;usage de la variance r\u00e9siduelle (exprim\u00e9e en unit\u00e9s au carr\u00e9) ou du coefficient de d\u00e9termination qui fournit un ratio sans dimension physique directe.<\/p>\n<p>Pour les chercheurs en sciences quantitatives, cette facilit\u00e9 d&rsquo;interpr\u00e9tation fluidifie la communication interdisciplinaire, permettant \u00e0 des praticiens, cliniciens ou d\u00e9cideurs institutionnels de saisir imm\u00e9diatement l&rsquo;impr\u00e9cision concr\u00e8te d&rsquo;un outil diagnostique ou d\u00e9cisionnel. Il importe n\u00e9anmoins de souligner avec rigueur la distinction entre l&rsquo;erreur moyenne arithm\u00e9tique et l&rsquo;esp\u00e9rance conditionnelle : la minimisation de la MAE identifie la <em>m\u00e9diane conditionnelle<\/em> de la distribution conditionnelle des observations sachant les pr\u00e9dicteurs, et non sa <em>moyenne conditionnelle<\/em>.<\/p>\n<p>Enfin, l&rsquo;interpr\u00e9tation purement unilat\u00e9rale des r\u00e9sidus absolus comporte des limites inh\u00e9rentes qu&rsquo;il convient de ne pas occulter. En amalgamant sous une m\u00eame valeur scalaire les sous-estimations et les surestimations, la MAE dissimule l&rsquo;asym\u00e9trie potentielle du comportement du mod\u00e8le. Deux algorithmes pr\u00e9dictifs affichant une MAE rigoureusement identique de 5 points peuvent pr\u00e9senter des profils cliniques diam\u00e9tralement oppos\u00e9s : l&rsquo;un sous-estimant syst\u00e9matiquement les manifestations pathologiques aigu\u00ebs, l&rsquo;autre sur\u00e9valuant de mani\u00e8re uniforme les \u00e9tats stables. L&rsquo;examen de la MAE doit par cons\u00e9quent \u00eatre articul\u00e9 avec une analyse conjointe du biais r\u00e9siduel moyen.<\/p>\n<h2>2. Formalisation math\u00e9matique et d\u00e9composition analytique de la formule<\/h2>\n<h3>2.1 Anatomie de l&rsquo;\u00e9quation standard de la MAE<\/h3>\n<p>La formalisation math\u00e9matique de l&rsquo;erreur absolue moyenne repose sur une structure alg\u00e9brique limpide, caract\u00e9ris\u00e9e par une d\u00e9composition indicielle explicite. Soit un \u00e9chantillon de <em>n<\/em> observations appari\u00e9es, o\u00f9 <em>y<sub>i<\/sub><\/em> repr\u00e9sente la valeur empirique observ\u00e9e pour le <em>i<\/em>-\u00e8me sujet ou unit\u00e9 statistique, et <em>x<sub>i<\/sub><\/em> (ou <em>\u0177<sub>i<\/sub><\/em>) d\u00e9signe l&rsquo;estimation mod\u00e9lis\u00e9e correspondante produite par le mod\u00e8le statistique. La formule standard s&rsquo;\u00e9nonce comme suit :<\/p>\n<p><strong>MAE = (1 \/ n) * &Sigma; |y<sub>i<\/sub> &minus; x<sub>i<\/sub>|<\/strong> (la sommation s&rsquo;\u00e9tendant de <em>i<\/em> = 1 jusqu&rsquo;\u00e0 <em>n<\/em>).<\/p>\n<p>Dans cette architecture analytique, chaque composante remplit une fonction fondamentale :<\/p>\n<ul>\n<li>Le terme <strong>n<\/strong> repr\u00e9sente la taille \u00e9chantillonnale effective ou le nombre de degr\u00e9s d&rsquo;observation ind\u00e9pendants. Il agit comme un facteur normalisateur scalaire, transformant la sommation totale des \u00e9carts en une grandeur intensive ind\u00e9pendante du volume de l&rsquo;\u00e9chantillon analys\u00e9.<\/li>\n<li>L&rsquo;op\u00e9rateur <strong>sommatoire (&Sigma;)<\/strong> assure l&rsquo;agr\u00e9gation s\u00e9quentielle des divergences individuelles sur l&rsquo;ensemble de l&rsquo;espace \u00e9chantillonnal.<\/li>\n<li>L&rsquo;op\u00e9rateur de <strong>valeur absolue |&middot;|<\/strong> supprime le signe alg\u00e9brique de la d\u00e9viation univari\u00e9e, garantissant que les vecteurs d&rsquo;erreurs r\u00e9siduelles de directions oppos\u00e9es s&rsquo;additionnent au lieu de se compenser.<\/li>\n<li>Les composantes <strong>y<sub>i<\/sub><\/strong> et <strong>x<sub>i<\/sub><\/strong> d\u00e9finissent respectivement la r\u00e9alit\u00e9 ontologique capt\u00e9e par la mesure instrumentale et la projection th\u00e9orique issue des postulats param\u00e9triques ou algorithmiques du mod\u00e8le.<\/li>\n<\/ul>\n<h3>2.2 Comparaison fonctionnelle entre perte L1 et perte L2<\/h3>\n<p>L&rsquo;analyse comparative des fonctions de perte met en lumi\u00e8re la divergence structurelle entre l&rsquo;optimisation sous contrainte de norme L1 (fondement de la MAE) et sous norme L2 (fondement de l&rsquo;erreur quadratique moyenne ou MSE). La perte L1, d\u00e9finie par <em>L<sub>1<\/sub>(y, x) = |y &minus; x|<\/em>, se caract\u00e9rise par une d\u00e9riv\u00e9e premi\u00e8re constante par morceaux : \u00e9gale \u00e0 &minus;1 lorsque l&rsquo;estimation exc\u00e8de l&rsquo;observation, et \u00e0 +1 dans le cas inverse. En corollaire, le gradient n&rsquo;augmente pas lorsque l&rsquo;erreur s&rsquo;accro\u00eet, privant l&rsquo;optimiseur d&rsquo;une force de rappel disproportionn\u00e9e envers les \u00e9carts majeurs.<\/p>\n<p>\u00c0 l&rsquo;inverse, la perte quadratique L2, formalis\u00e9e par <em>L<sub>2<\/sub>(y, x) = (y &minus; x)<sup>2<\/sup><\/em>, g\u00e9n\u00e8re un gradient lin\u00e9aire <em>2(x &minus; y)<\/em>. Cette structure conf\u00e8re une p\u00e9nalit\u00e9 exponentielle aux r\u00e9sidus \u00e9loign\u00e9s. Sur le plan de la th\u00e9orie asymptotique de l&rsquo;estimation, ce distinguo poss\u00e8de des r\u00e9percussions consid\u00e9rables : minimiser la somme des pertes quadratiques conduit asymptotiquement \u00e0 l&rsquo;estimation de la moyenne conditionnelle de la variable cible (<em>E[Y|X]<\/em>), tandis que la minimisation de la perte L1 converge vers la m\u00e9diane conditionnelle (<em>Med[Y|X]<\/em>).<\/p>\n<p>Dans les contextes exp\u00e9rimentaux marqu\u00e9s par des asym\u00e9tries de distribution prononc\u00e9es ou par la pr\u00e9sence de contaminations stochastiques \u00e0 queue lourde, l&rsquo;estimateur de la m\u00e9diane conditionnelle d\u00e9riv\u00e9 de la norme L1 pr\u00e9sente des propri\u00e9t\u00e9s de robustesse bien document\u00e9es par la statistique non param\u00e9trique. La m\u00e9diane \u00e9tant caract\u00e9ris\u00e9e par un point de rupture de 50 %, l&rsquo;estimation de la MAE r\u00e9siste \u00e0 des contaminations locales de donn\u00e9es qui d\u00e9stabiliseraient irr\u00e9m\u00e9diablement toute estimation quadratique.<\/p>\n<h3>2.3 D\u00e9rivation de la MAE pond\u00e9r\u00e9e et variantes th\u00e9oriques<\/h3>\n<p>Dans de multiples contextes empiriques, l&rsquo;hypoth\u00e8se sous-jacente d&rsquo;homosc\u00e9dasticit\u00e9 et d&rsquo;\u00e9quiprobabilit\u00e9 des observations ne peut \u00eatre soutenue. Lorsque la fiabilit\u00e9 des instruments de mesure varie \u00e0 travers l&rsquo;\u00e9chantillon, ou lorsque certaines strates d\u00e9mographiques sont sous-repr\u00e9sent\u00e9es au sein d&rsquo;un protocole d&rsquo;\u00e9chantillonnage complexe, le recours \u00e0 l&rsquo;erreur absolue moyenne pond\u00e9r\u00e9e (WMAE, pour <em>Weighted Mean Absolute Error<\/em>) devient math\u00e9matiquement indispensable. Sa formulation analytique int\u00e8gre un vecteur de poids normalis\u00e9s <em>w<sub>i<\/sub><\/em> :<\/p>\n<p><strong>WMAE = (&Sigma; w<sub>i<\/sub> * |y<sub>i<\/sub> &minus; x<sub>i<\/sub>|) \/ (&Sigma; w<sub>i<\/sub>)<\/strong><\/p>\n<p>Les poids <em>w<sub>i<\/sub><\/em> peuvent \u00eatre inversement proportionnels \u00e0 la variance d&rsquo;erreur de mesure du sous-groupe consid\u00e9r\u00e9 (pond\u00e9ration par la pr\u00e9cision) ou proportionnels \u00e0 l&rsquo;inverse de la probabilit\u00e9 d&rsquo;inclusion dans l&rsquo;\u00e9chantillon (redressement par calage sur marges en \u00e9pid\u00e9miologie ou psychologie sociale).<\/p>\n<p>Une seconde d\u00e9clinaison th\u00e9orique concerne la normalisation scalaire de la m\u00e9trique, aboutissant \u00e0 la MAE normalis\u00e9e (NMAE). Cette variante vise \u00e0 neutraliser l&rsquo;impact de l&rsquo;amplitude de la variable d&rsquo;int\u00e9r\u00eat afin de rendre possibles des comparaisons inter-domaines. La normalisation s&rsquo;op\u00e8re usuellement soit en divisant la MAE par l&rsquo;\u00e9tendue des valeurs observ\u00e9es (<em>y<sub>max<\/sub> &minus; y<sub>min<\/sub><\/em>), soit par la moyenne empirique de l&rsquo;\u00e9chantillon (<em>y\u0304<\/em>), produisant une quantit\u00e9 adimensionnelle. Enfin, l&rsquo;int\u00e9gration de fonctions de co\u00fbt asym\u00e9triques (par exemple, la r\u00e9gression quantile et la perte de type <em>pinball<\/em>) permet de surpond\u00e9rer d\u00e9lib\u00e9r\u00e9ment les sous-estimations par rapport aux surestimations, r\u00e9pondant aux exigences d\u00e9ontologiques ou pratiques o\u00f9 un faux n\u00e9gatif entra\u00eene des cons\u00e9quences cliniques infiniment plus dommageables qu&rsquo;un faux positif.<\/p>\n<h2>3. Analyse comparative : MAE, RMSE, MSE et MAPE dans la mod\u00e9lisation statistique<\/h2>\n<h3>3.1 MAE versus RMSE (Root Mean Squared Error) : dynamique de sensibilit\u00e9<\/h3>\n<p>Le d\u00e9bat m\u00e9thodologique opposant la MAE \u00e0 la racine de l&rsquo;erreur quadratique moyenne (<a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Root-mean-square_deviation\">RMSE<\/a>) a suscit\u00e9 d&rsquo;abondantes controverses dans la litt\u00e9rature statistique contemporaine, notamment stimul\u00e9es par les travaux de Willmott et Matsuura d\u00e9montrant les distorsions induites par la RMSE dans les \u00e9valuations g\u00e9ophysiques et environnementales. La divergence essentielle r\u00e9side dans l&rsquo;amplification quadratique des \u00e9carts op\u00e9r\u00e9e par la RMSE : en \u00e9levant chaque r\u00e9sidu au carr\u00e9 avant d&rsquo;en calculer la moyenne puis d&rsquo;en extraire la racine, la RMSE p\u00e9nalise s\u00e9v\u00e8rement les erreurs de grande amplitude.<\/p>\n<p>Cette divergence fonctionnelle autorise une exploitation diagnostique de grande valeur : le ratio <strong>RMSE \/ MAE<\/strong>. Sur le plan th\u00e9orique, la RMSE est math\u00e9matiquement sup\u00e9rieure ou \u00e9gale \u00e0 la MAE (en vertu de l&rsquo;in\u00e9galit\u00e9 de Cauchy-Schwarz). Le ratio atteint sa borne inf\u00e9rieure rigoureuse de 1 lorsque toutes les erreurs individuelles ont exactement la m\u00eame amplitude absolue (<em>|e<sub>1<\/sub>| = |e<sub>2<\/sub>| = &#8230; = |e<sub>n<\/sub>|<\/em>). \u00c0 l&rsquo;inverse, plus la dispersion et la variabilit\u00e9 de la magnitude des erreurs augmentent \u2014 notamment en pr\u00e9sence de quelques d\u00e9viations monumentales \u2014, plus le ratio RMSE\/MAE diverge vers des valeurs \u00e9lev\u00e9es.<\/p>\n<p>Dans les jeux de donn\u00e9es psychom\u00e9triques o\u00f9 des profils neuroatypiques ou des r\u00e9ponses d&rsquo;inattention peuvent g\u00e9n\u00e9rer des discordances ponctuelles consid\u00e9rables, la RMSE tend \u00e0 surdiagnostiquer une mauvaise ad\u00e9quation globale du mod\u00e8le, alors que la MAE rend fid\u00e8lement compte du niveau de pr\u00e9cision pr\u00e9dictive typique pour la vaste majorit\u00e9 des sujets. Le choix entre ces deux indicateurs doit imp\u00e9rativement d\u00e9couler des objectifs de validit\u00e9 interne du chercheur : s&rsquo;agit-il de minimiser les d\u00e9faillances catastrophiques (privil\u00e9gier la RMSE) ou de garantir la meilleure exactitude m\u00e9diane globale (privil\u00e9gier la MAE) ?<\/p>\n<h3>3.2 MAE versus MSE (Mean Squared Error) : dimensions et interpr\u00e9tation<\/h3>\n<p>La confrontation entre l&rsquo;erreur absolue moyenne et l&rsquo;erreur quadratique moyenne (MSE) s&rsquo;articule autour d&rsquo;un arbitrage classique entre tractabilit\u00e9 computationnelle et lisibilit\u00e9 op\u00e9rationnelle. L&rsquo;\u00e9l\u00e9vation au carr\u00e9 inh\u00e9rente \u00e0 la MSE alt\u00e8re irr\u00e9versiblement l&rsquo;unit\u00e9 de mesure de la variable d\u00e9pendante : une estimation d&rsquo;un temps de r\u00e9action mesur\u00e9 en millisecondes produira une MSE exprim\u00e9e en <em>millisecondes au carr\u00e9<\/em>, une entit\u00e9 dimensionnelle d\u00e9pourvue de toute r\u00e9sonance clinique intuitive.<\/p>\n<p>N\u00e9anmoins, la MSE demeure la m\u00e9trique de pr\u00e9dilection au sein des algorithmes d&rsquo;optimisation num\u00e9rique, sp\u00e9cifiquement les m\u00e9thodes de descente de gradient standard, les mod\u00e8les de r\u00e9gression lin\u00e9aire par moindres carr\u00e9s ordinaires (MCO) et les r\u00e9seaux de neurones profonds. La diff\u00e9rentiabilit\u00e9 continue de la fonction quadratique permet d&rsquo;obtenir des gradients analytiques \u00e9l\u00e9gants et directeurs, facilitant une convergence num\u00e9rique rapide et sans heurts.<\/p>\n<p>\u00c0 l&rsquo;oppos\u00e9, la MAE s&rsquo;impose comme l&rsquo;outil privil\u00e9gi\u00e9 pour la restitution descriptive et la communication des performances statistiques aupr\u00e8s de panels d&rsquo;experts non statisticiens. Elle pr\u00e9serve l&rsquo;\u00e9chelle naturelle du ph\u00e9nom\u00e8ne observ\u00e9, \u00e9liminant tout artefact d&rsquo;interpr\u00e9tation. En pratique de mod\u00e9lisation, il est courant d&rsquo;optimiser les param\u00e8tres internes d&rsquo;un mod\u00e8le sous contrainte de MSE (en raison de sa convexit\u00e9 et diff\u00e9rentiabilit\u00e9 fluide) tout en \u00e9valuant et en communiquant la pr\u00e9cision op\u00e9rationnelle finale via la MAE.<\/p>\n<h3>3.3 MAE versus MAPE (Mean Absolute Percentage Error) : limites d&rsquo;\u00e9chelle<\/h3>\n<p>L&rsquo;erreur absolue moyenne en pourcentage (<a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Mean_absolute_percentage_error\">MAPE<\/a>) tente de r\u00e9soudre la question de l&rsquo;interpr\u00e9tabilit\u00e9 en exprimant la divergence r\u00e9siduelle sous la forme d&rsquo;un pourcentage relatif de la valeur observ\u00e9e :<\/p>\n<p><strong>MAPE = (100% \/ n) * &Sigma; |(y<sub>i<\/sub> &minus; x<sub>i<\/sub>) \/ y<sub>i<\/sub>|<\/strong><\/p>\n<p>Si la formulation en pourcentage pr\u00e9sente une s\u00e9duction \u00e9vidente pour les gestionnaires et \u00e9conomistes en qu\u00eate d&rsquo;un indicateur universellement comparable, elle introduit des pathologies statistiques s\u00e9v\u00e8res qui rendent son usage extr\u00eamement risqu\u00e9 dans de nombreuses disciplines quantitatives. La premi\u00e8re pathologie r\u00e9side dans son instabilit\u00e9 asymptotique \u00e0 l&rsquo;approche de z\u00e9ro : d\u00e8s lors qu&rsquo;une valeur observ\u00e9e <em>y<sub>i<\/sub><\/em> tend vers z\u00e9ro, le d\u00e9nominateur de l&rsquo;\u00e9quation s&rsquo;annule, provoquant l&rsquo;explosion math\u00e9matique du ratio vers l&rsquo;infini. Si un jeu de donn\u00e9es comprend ne serait-ce qu&rsquo;une seule observation \u00e9gale \u00e0 z\u00e9ro, le calcul de la MAPE devient purement ind\u00e9fini.<\/p>\n<p>La seconde distorsion majeure concerne l&rsquo;asym\u00e9trie de p\u00e9nalisation inh\u00e9rente \u00e0 la MAPE. En raison de la pr\u00e9sence de la valeur observ\u00e9e <em>y<sub>i<\/sub><\/em> au d\u00e9nominateur, la m\u00e9trique inflige une p\u00e9nalit\u00e9 consid\u00e9rablement plus lourde lorsque la valeur pr\u00e9dite est sup\u00e9rieure \u00e0 la valeur observ\u00e9e que lorsqu&rsquo;elle lui est inf\u00e9rieure. Une surestimation d&rsquo;un facteur 2 produit une erreur de 100 %, tandis qu&rsquo;une sous-estimation d&rsquo;un facteur 2 (pr\u00e9diction \u00e9gale \u00e0 la moiti\u00e9 de l&rsquo;observation) ne produit qu&rsquo;une erreur de 50 %. Cette distorsion syst\u00e9matique favorise artificiellement les mod\u00e8les produisant des pr\u00e9dictions structurellement biais\u00e9es vers le bas. Pour les scores psychologiques standardis\u00e9s (tels que les scores de QI ou les indices de d\u00e9pression standardis\u00e9s en scores Z ou T), o\u00f9 le point z\u00e9ro est purement arbitraire et non absolu, la MAPE est conceptuellement caduque, consacrant la MAE comme la seule m\u00e9trique robuste et universellement l\u00e9gitime.<\/p>\n<h2>4. Configuration de l&rsquo;environnement de programmation statistique R<\/h2>\n<h3>4.1 Initialisation de la session et gestion des d\u00e9pendances<\/h3>\n<p>La mise en \u0153uvre d&rsquo;analyses quantitatives rigoureuses dans le langage R requiert la configuration m\u00e9thodique d&rsquo;un environnement de travail reproductible. Avant tout traitement matriciel ou mod\u00e9lisation statistique, l&rsquo;initialisation de la session implique l&rsquo;inspection des versions d&rsquo;ex\u00e9cution du moteur R et la standardisation des options num\u00e9riques globales. L&rsquo;instruction R <code>sessionInfo()<\/code> permet de consigner formellement la version exacte du syst\u00e8me d&rsquo;exploitation et des compilateurs sous-jacents, gage indispensable de tra\u00e7abilit\u00e9 scientifique.<\/p>\n<p>La gestion de la pr\u00e9cision num\u00e9rique s&rsquo;effectue via l&rsquo;ajustement des options d&rsquo;environnement. Il est recommand\u00e9 de param\u00e9trer explicitement le nombre de chiffres significatifs restitu\u00e9s dans la console via l&rsquo;instruction <code>options(digits = 7, scipen = 999)<\/code>, \u00e9liminant ainsi l&rsquo;affichage intempestif en notation scientifique pour les r\u00e9sidus de faible amplitude et harmonisant la pr\u00e9cision des nombres \u00e0 virgule flottante r\u00e9pondant \u00e0 la norme IEEE 754. De surcro\u00eet, la d\u00e9finition d&rsquo;un r\u00e9pertoire de travail explicite via <code>setwd()<\/code> ou l&rsquo;utilisation structur\u00e9e de projets RStudio coupl\u00e9s au package <code>here<\/code> garantit l&rsquo;ind\u00e9pendance des scripts \u00e0 l&rsquo;\u00e9gard de l&rsquo;arborescence des syst\u00e8mes de fichiers locaux.<\/p>\n<p>Enfin, la reproductibilit\u00e9 parfaite des simulations stochastiques et des d\u00e9coupages d&rsquo;\u00e9chantillons exige l&rsquo;initialisation d\u00e9lib\u00e9r\u00e9e du g\u00e9n\u00e9rateur congruentiel de nombres pseudo-al\u00e9atoires au moyen de la fonction primitive <code>set.seed()<\/code>. Le recours \u00e0 des sessions sandbox\u00e9es, orchestr\u00e9es via des gestionnaires de d\u00e9pendances herm\u00e9tiques tels que <code>renv<\/code>, pr\u00e9vient toute alt\u00e9ration des r\u00e9sultats statistiques induite par des mises \u00e0 jour r\u00e9tro-incompatibles des packages du CRAN.<\/p>\n<h3>4.2 Pr\u00e9sentation et installation du package Metrics<\/h3>\n<p>Parmi les extensions sp\u00e9cialis\u00e9es disponibles sur le Comprehensive R Archive Network (<a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/\">CRAN<\/a>), le package <code>Metrics<\/code> constitue une r\u00e9f\u00e9rence historique pour l&rsquo;\u00e9valuation empirique des algorithmes d&rsquo;apprentissage supervis\u00e9 et de r\u00e9gression statistique. D\u00e9velopp\u00e9 initialement pour fournir des routines de m\u00e9triques standardis\u00e9es, ce package regroupe un ensemble exhaustif de fonctions vectoris\u00e9es calculant les pertes de pr\u00e9diction courantes, incluant l&rsquo;exactitude de classification, la d\u00e9viance logistique, la RMSE et la MAE.<\/p>\n<p>L&rsquo;installation du package depuis les miroirs officiels du CRAN s&rsquo;op\u00e8re au moyen de la commande suivante ex\u00e9cut\u00e9e dans l&rsquo;interpr\u00e9teur interactif :<\/p>\n<p><code>install.packages(\"Metrics\")<\/code><\/p>\n<p>Cette proc\u00e9dure compile et d\u00e9ploie les d\u00e9pendances n\u00e9cessaires au fonctionnement du package. Son chargement en m\u00e9moire active s&rsquo;effectue ensuite par l&rsquo;instruction :<\/p>\n<p><code>library(Metrics)<\/code><\/p>\n<p>Il convient de souligner que le package <code>Metrics<\/code> se caract\u00e9rise par une empreinte m\u00e9moire extr\u00eamement l\u00e9g\u00e8re et une absence quasi totale de d\u00e9pendances externes lourdes, ce qui en fait un composant de choix pour les environnements de production l\u00e9gers ou les scripts d&rsquo;analyse rapide. L&rsquo;exploration de son espace de noms r\u00e9v\u00e8le que sa fonction d\u00e9di\u00e9e <code>mae()<\/code> accepte comme arguments minimaux deux vecteurs num\u00e9riques de dimensions identiques, repr\u00e9sentant respectivement les observations empiriques et les estimations mod\u00e9lis\u00e9es.<\/p>\n<h3>4.3 Alternatives logicielles contemporaines au sein du CRAN<\/h3>\n<p>Bien que le package <code>Metrics<\/code> r\u00e9ponde parfaitement aux besoins \u00e9l\u00e9mentaires d&rsquo;\u00e9valuation m\u00e9trique, le paysage logiciel de R s&rsquo;est consid\u00e9rablement enrichi au cours de la derni\u00e8re d\u00e9cennie, donnant naissance \u00e0 des architectures d&rsquo;analyse unifi\u00e9es adapt\u00e9es aux flux de donn\u00e9es massifs et complexes :<\/p>\n<ul>\n<li>Le m\u00e9ta-framework <strong>tidymodels<\/strong>, et plus sp\u00e9cifiquement son composant d&rsquo;\u00e9valuation modulaire <strong>yardstick<\/strong>, s&rsquo;impose aujourd&rsquo;hui comme le standard moderne de la science des donn\u00e9es en R. Con\u00e7u pour s&rsquo;int\u00e9grer harmonieusement avec les tables de donn\u00e9es <code>tibble<\/code> et la syntaxe fluide du <code>tidyverse<\/code>, il permet des calculs m\u00e9triques structur\u00e9s, respectueux du groupement fonctionnel des donn\u00e9es et directement compatibles avec l&rsquo;\u00e9valuation distribu\u00e9e.<\/li>\n<li>L&rsquo;environnement historique <strong>caret<\/strong> (<em>Classification And REgression Training<\/em>), con\u00e7u par Max Kuhn, int\u00e8gre nativement le calcul de la MAE pour orchestrer les grilles d&rsquo;optimisation d&rsquo;hyperparam\u00e8tres et la s\u00e9lection automatique de mod\u00e8les pr\u00e9dictifs.<\/li>\n<li>Des biblioth\u00e8ques ultra-l\u00e9g\u00e8res et hautement optimis\u00e9es telles que <strong>MLmetrics<\/strong> ou <strong>ModelMetrics<\/strong> proposent des impl\u00e9mentations de m\u00e9triques pr\u00e9dictives \u00e9crites directement en C++ via <code>Rcpp<\/code>. <code>ModelMetrics<\/code>, en particulier, excelle dans la vitesse d&rsquo;ex\u00e9cution lors du traitement de vecteurs g\u00e9ants comptant des dizaines de millions d&rsquo;observations.<\/li>\n<\/ul>\n<p>La s\u00e9lection du package ad\u00e9quat d\u00e9pend donc intrins\u00e8quement de l&rsquo;architecture logicielle globale du projet : une fonction personnalis\u00e9e en R base pour une autonomie totale sans d\u00e9pendance externe, <code>Metrics<\/code> pour la concision scripturale, <code>yardstick<\/code> pour des pipelines statistiques \u00e9l\u00e9gants et modernes, ou <code>ModelMetrics<\/code> pour des contraintes d&rsquo;optimisation computationnelle critique.<\/p>\n<h2>5. Calcul de la MAE sans package externe : impl\u00e9mentation en R Base<\/h2>\n<h3>5.1 Vectorisation native et programmation d&rsquo;une fonction personnalis\u00e9e<\/h3>\n<p>L&rsquo;un des atouts les plus remarquables du langage R r\u00e9side dans son paradigme vectoriel natif, h\u00e9rit\u00e9 de la logique matricielle de S. Dans ce mod\u00e8le d&rsquo;ex\u00e9cution, les op\u00e9rations math\u00e9matiques appliqu\u00e9es \u00e0 des structures vectorielles ne n\u00e9cessitent aucun recours explicite \u00e0 des boucles it\u00e9ratives au niveau de l&rsquo;utilisateur. Les fonctions primitives op\u00e8rent directement sur l&rsquo;ensemble des \u00e9l\u00e9ments via des routines sous-jacentes hautement optimis\u00e9es et compil\u00e9es en C et Fortran.<\/p>\n<p>L&rsquo;impl\u00e9mentation algorithmique de la MAE en R base d\u00e9coule directement de la combinaison de la fonction de soustraction vectorielle, de la primitive <code>abs()<\/code> calculant la valeur absolue terme \u00e0 terme, et de la fonction primitive <code>mean()<\/code> calculant l&rsquo;esp\u00e9rance arithm\u00e9tique. Consid\u00e9rons deux vecteurs simul\u00e9s :<\/p>\n<p><code>valeurs_observees &lt;- c(14.2, 16.8, 12.1, 15.5, 18.9)<\/code><br \/>\n<code>valeurs_predites &lt;- c(13.9, 17.2, 11.5, 16.1, 18.0)<\/code><\/p>\n<p>Le calcul vectoriel s&rsquo;ex\u00e9cute de mani\u00e8re limpide : l&rsquo;op\u00e9ration <code>valeurs_observees - valeurs_predites<\/code> g\u00e9n\u00e8re un vecteur r\u00e9siduel <em>(0.3, -0.4, 0.6, -0.6, 0.9)<\/em>. L&rsquo;application de <code>abs()<\/code> transforme imm\u00e9diatement ce r\u00e9sultat en d\u00e9viations absolues non n\u00e9gatives <em>(0.3, 0.4, 0.6, 0.6, 0.9)<\/em>. L&rsquo;invocation terminale de <code>mean()<\/code> agr\u00e8ge ces amplitudes pour retourner l&rsquo;erreur absolue moyenne finale, soit 0.56.<\/p>\n<p>Afin de conf\u00e9rer \u00e0 cette s\u00e9quence un caract\u00e8re modulaire et p\u00e9renne, il est m\u00e9thodologiquement recommand\u00e9 d&rsquo;encapsuler cette logique au sein d&rsquo;une fonction personnalis\u00e9e rigoureusement param\u00e9tr\u00e9e :<\/p>\n<p><code>calcul_mae_base &lt;- function(reel, pred, na.rm = TRUE) {<br \/>\n&nbsp;&nbsp;return(mean(abs(reel - pred), na.rm = na.rm))<br \/>\n}<\/code><\/p>\n<p>Cette fonction \u00e9l\u00e9mentaire constitue une alternative universelle, dispensant int\u00e9gralement l&rsquo;analyste de l&rsquo;importation de biblioth\u00e8ques tierces au sein de son environnement de calcul.<\/p>\n<h3>5.2 Gestion des valeurs manquantes et contraintes vectorielles<\/h3>\n<p>Dans les contextes d&rsquo;exp\u00e9rimentation r\u00e9elle et de collecte de donn\u00e9es observationnelles, l&rsquo;int\u00e9grit\u00e9 des vecteurs de mesures est fr\u00e9quemment alt\u00e9r\u00e9e par la survenue d&rsquo;\u00e9l\u00e9ments manquants de type <code>NA<\/code> (<em>Not Available<\/em>), r\u00e9sultant d&rsquo;abandons de participants, de d\u00e9faillances de capteurs ou de donn\u00e9es aberrantes pr\u00e9alablement expurg\u00e9es. En R, la propagation des valeurs manquantes ob\u00e9it \u00e0 une logique stricte : toute op\u00e9ration arithm\u00e9tique impliquant un \u00e9l\u00e9ment <code>NA<\/code> renvoie invariablement <code>NA<\/code>.<\/p>\n<p>Si un vecteur d&rsquo;observations ou de pr\u00e9dictions contient ne serait-ce qu&rsquo;une seule occurrence non assign\u00e9e, l&rsquo;appel non param\u00e9tr\u00e9 de la fonction <code>mean(abs(reel - pred))<\/code> \u00e9chouera irr\u00e9m\u00e9diablement en renvoyant <code>NA<\/code>. Pour pr\u00e9server la validit\u00e9 du calcul sans fausser la taille d&rsquo;\u00e9chantillon restante, le param\u00e9trage explicite de l&rsquo;argument <code>na.rm = TRUE<\/code> au sein de la primitive <code>mean()<\/code> est indispensable. Cette directive ordonne au moteur d&rsquo;ex\u00e9cution d&rsquo;ignorer les paires d&rsquo;observations incompl\u00e8tes lors de la sommation et de recalculer le d\u00e9nominateur <em>n<\/em> sur la base exclusive des cas complets.<\/p>\n<p>Par ailleurs, la robustesse logicielle exige une v\u00e9rification formelle des contraintes dimensionnelles des structures vectorielles d&rsquo;entr\u00e9e. En vertu du m\u00e9canisme de recyclage automatique des vecteurs propre \u00e0 R, si l&rsquo;on tente de soustraire un vecteur de longueur 5 \u00e0 un vecteur de longueur 10, R ne d\u00e9clenchera pas n\u00e9cessairement d&rsquo;erreur fatale mais recyclera de fa\u00e7on cyclique les \u00e9l\u00e9ments du vecteur court, produisant des calculs de r\u00e9sidus silencieusement corrompus et d\u00e9nu\u00e9s de sens scientifique. L&rsquo;int\u00e9gration de tests d&rsquo;assertion formels \u00e0 l&rsquo;aide de la fonction <code>stopifnot()<\/code> au sein de la fonction d&rsquo;estimation s&rsquo;av\u00e8re donc imp\u00e9rative :<\/p>\n<p><code>calcul_mae_robuste &lt;- function(reel, pred, na.rm = TRUE) {<br \/>\n&nbsp;&nbsp;stopifnot(is.numeric(reel), is.numeric(pred))<br \/>\n&nbsp;&nbsp;stopifnot(length(reel) == length(pred))<br \/>\n&nbsp;&nbsp;return(mean(abs(reel - pred), na.rm = na.rm))<br \/>\n}<\/code><\/p>\n<h3>5.3 Analyse de performance et efficacit\u00e9 computationnelle<\/h3>\n<p>L&rsquo;efficacit\u00e9 computationnelle du calcul de la MAE en R base repose enti\u00e8rement sur l&rsquo;\u00e9vitement d\u00e9lib\u00e9r\u00e9 des structures it\u00e9ratives explicites du type boucle <code>for<\/code>. Dans les langages interpr\u00e9t\u00e9s dynamiquement, chaque it\u00e9ration d&rsquo;une boucle <code>for<\/code> entra\u00eene une surcharge d&rsquo;interpr\u00e9tation substantielle, induisant des co\u00fbts cognitifs et temporels prohibitifs lorsque le nombre d&rsquo;observations se chiffre en millions.<\/p>\n<p>Le benchmarking empirique \u00e0 l&rsquo;aide du package sp\u00e9cialis\u00e9 <code>microbenchmark<\/code> d\u00e9montre de mani\u00e8re spectaculaire cette divergence. Lorsque l&rsquo;on confronte une boucle it\u00e9rative classique \u2014 initialisant un accumulateur scalaire et incr\u00e9mentant la somme absolue \u00e0 travers chaque \u00e9l\u00e9ment \u2014 \u00e0 la formulation vectoris\u00e9e <code>mean(abs(reel - pred))<\/code> sur un \u00e9chantillon synth\u00e9tique de 10 millions d&rsquo;observations g\u00e9n\u00e9r\u00e9 al\u00e9atoirement, l&rsquo;impl\u00e9mentation vectoris\u00e9e native se r\u00e9v\u00e8le g\u00e9n\u00e9ralement entre 15 et 40 fois plus rapide. Ce diff\u00e9rentiel spectaculaire s&rsquo;explique par l&rsquo;ex\u00e9cution de la boucle interne au niveau du code C sous-jacent, exploitant de mani\u00e8re optimale les registres de m\u00e9moire vectorielle du processeur central (instructions SIMD).<\/p>\n<p>Sur le plan de l&#8217;empreinte m\u00e9moire, l&rsquo;impl\u00e9mentation directe <code>mean(abs(reel - pred))<\/code> requiert n\u00e9anmoins l&rsquo;allocation transitoire d&rsquo;un vecteur num\u00e9rique interm\u00e9diaire pour stocker les r\u00e9sultats de la diff\u00e9rence vectorielle, ainsi que d&rsquo;un second vecteur pour les valeurs absolues correspondantes. Pour des volum\u00e9tries massives d\u00e9passant la capacit\u00e9 de la m\u00e9moire vive (RAM), l&rsquo;analyste se tournera vers des fonctions compil\u00e9es via <code>Rcpp<\/code> op\u00e9rant un calcul d&rsquo;accumulation en passe unique, \u00e9liminant ainsi toute duplication de vecteurs temporaires au sein de l&rsquo;environnement R.<\/p>\n<h2>6. Impl\u00e9mentation pratique avec le package Metrics sur des vecteurs<\/h2>\n<h3>6.1 Application \u00e9l\u00e9mentaire de la fonction mae() sur deux vecteurs<\/h3>\n<p>L&rsquo;exploitation op\u00e9rationnelle du package <code>Metrics<\/code> vise \u00e0 standardiser les proc\u00e9dures d&rsquo;\u00e9valuation statistique \u00e0 travers une syntaxe concise, directement inspir\u00e9e de la nomenclature math\u00e9matique. Au c\u0153ur de cette biblioth\u00e8que figure la fonction <code>mae()<\/code>, con\u00e7ue pour op\u00e9rer sur des structures vectorielles \u00e0 une dimension.<\/p>\n<p>La syntaxe canonique d&rsquo;appel de cette fonction se structure formellement autour de deux arguments positionnels obligatoires :<\/p>\n<p><code>mae(actual, predicted)<\/code><\/p>\n<p>Le premier argument, <code>actual<\/code>, re\u00e7oit le vecteur des valeurs empiriques r\u00e9elles constat\u00e9es exp\u00e9rimentalement. Le second argument, <code>predicted<\/code>, accueille le vecteur des estimations mod\u00e9lis\u00e9es issues du mod\u00e8le statistique test\u00e9. Sous le capot, la fonction applique une s\u00e9quence d&rsquo;instructions similaire \u00e0 l&rsquo;impl\u00e9mentation R base vectoris\u00e9e, tout en int\u00e9grant des v\u00e9rifications internes de typage et de structure.<\/p>\n<p>L&rsquo;ex\u00e9cution de cette commande retourne un scalaire num\u00e9rique unique repr\u00e9sentant l&rsquo;erreur absolue moyenne globale. L&rsquo;interpr\u00e9tation de ce r\u00e9sultat doit imp\u00e9rativement s&rsquo;accompagner d&rsquo;une contextualisation : une MAE \u00e9gale \u00e0 2,45 n&rsquo;a de sens statistique et pratique qu&rsquo;au regard de la dispersion intrins\u00e8que de la variable cible, appr\u00e9ci\u00e9e par exemple au travers de son \u00e9cart-type ou de son intervalle interquartile.<\/p>\n<h3>6.2 Analyse approfondie de l&rsquo;exemple num\u00e9rique classique<\/h3>\n<p>Afin de d\u00e9composer pas \u00e0 pas la m\u00e9canique arithm\u00e9tique sous-jacente \u00e0 la fonction <code>mae()<\/code> du package <code>Metrics<\/code>, examinons en d\u00e9tail une s\u00e9rie canonique d&rsquo;observations et de pr\u00e9dictions fr\u00e9quemment mobilis\u00e9e dans la litt\u00e9rature didactique :<\/p>\n<p>Soit le vecteur des observations r\u00e9elles :<br \/>\n<strong>y<\/strong> = (12, 13, 14, 15, 15, 22, 27, 29, 29, 30, 32)<br \/>\nSoit le vecteur des estimations mod\u00e9lis\u00e9es correspondantes :<br \/>\n<strong>x<\/strong> = (11, 13, 14, 14, 16, 19, 24, 30, 32, 36, 30)<\/p>\n<p>L&rsquo;\u00e9chantillon comporte <em>n<\/em> = 11 observations appari\u00e9es. La d\u00e9composition s\u00e9quentielle des r\u00e9sidus directionnels (<em>y<sub>i<\/sub> &minus; x<sub>i<\/sub><\/em>) et de leurs amplitudes absolues respectives (<em>|y<sub>i<\/sub> &minus; x<sub>i<\/sub>|<\/em>) s&rsquo;\u00e9tablit comme suit :<\/p>\n<ul>\n<li>Observation 1 : |12 &minus; 11| = |1| = 1<\/li>\n<li>Observation 2 : |13 &minus; 13| = |0| = 0<\/li>\n<li>Observation 3 : |14 &minus; 14| = |0| = 0<\/li>\n<li>Observation 4 : |15 &minus; 14| = |1| = 1<\/li>\n<li>Observation 5 : |15 &minus; 16| = |-1| = 1<\/li>\n<li>Observation 6 : |22 &minus; 19| = |3| = 3<\/li>\n<li>Observation 7 : |27 &minus; 24| = |3| = 3<\/li>\n<li>Observation 8 : |29 &minus; 30| = |-1| = 1<\/li>\n<li>Observation 9 : |29 &minus; 32| = |-3| = 3<\/li>\n<li>Observation 10 : |30 &minus; 36| = |-6| = 6<\/li>\n<li>Observation 11 : |32 &minus; 30| = |2| = 2<\/li>\n<\/ul>\n<p>L&rsquo;agr\u00e9gation par sommation de ces d\u00e9viations absolues individuelles produit la somme totale des \u00e9carts :<br \/>\n&Sigma; |y<sub>i<\/sub> &minus; x<sub>i<\/sub>| = 1 + 0 + 0 + 1 + 1 + 3 + 3 + 1 + 3 + 6 + 2 = 21.<\/p>\n<p>La division normalisatrice par le nombre total d&rsquo;observations (<em>n<\/em> = 11) fournit le r\u00e9sultat analytique :<br \/>\nMAE = 21 \/ 11 &approx; <strong>1.909091<\/strong>.<\/p>\n<p>Lorsque ces vecteurs sont d\u00e9clar\u00e9s sous R et transmis \u00e0 l&rsquo;instruction <code>Metrics::mae(y, x)<\/code>, l&rsquo;interpr\u00e9teur renvoie exactement la valeur flottante <code>1.909091<\/code>, validant empiriquement la parfaite concordance entre la r\u00e9solution math\u00e9matique formelle et le calcul logiciel automatis\u00e9.<\/p>\n<h3>6.3 D\u00e9pannage syntaxique et erreurs courantes lors de l&rsquo;usage de Metrics<\/h3>\n<p>Bien que la manipulation du package <code>Metrics<\/code> se distingue par son accessibilit\u00e9 conceptuelle, plusieurs \u00e9cueils syntaxiques et structurels peuvent alt\u00e9rer le d\u00e9roulement des calculs m\u00e9triques au sein de flux d&rsquo;analyses statistiques :<\/p>\n<p>L&rsquo;erreur la plus fr\u00e9quente concerne l&rsquo;interversion des arguments au sein de fonctions de perte asym\u00e9triques. Si, pour la MAE stricte, l&rsquo;inversion des param\u00e8tres <code>actual<\/code> et <code>predicted<\/code> n&rsquo;affecte pas le r\u00e9sultat en raison de la propri\u00e9t\u00e9 de commutativit\u00e9 de la valeur absolue (|y &minus; x| = |x &minus; y|), cette mauvaise habitude de programmation s&rsquo;av\u00e8re d\u00e9sastreuse d\u00e8s lors que l&rsquo;analyste transpose son code vers d&rsquo;autres fonctions du package (telles que <code>mape()<\/code> ou les d\u00e9viances logistiques) o\u00f9 l&rsquo;argument <code>actual<\/code> figure au d\u00e9nominateur.<\/p>\n<p>Un second \u00e9cueil critique r\u00e9side dans les conflits de types de donn\u00e9es. La fonction <code>mae()<\/code> exige des vecteurs strictement num\u00e9riques de classe <code>integer<\/code> ou <code>numeric<\/code>. La pr\u00e9sence accidentelle d&rsquo;une variable cat\u00e9gorielle encod\u00e9e sous la forme de facteurs (<code>factor<\/code>) ou de cha\u00eenes de caract\u00e8res (<code>character<\/code>) \u2014 situation fr\u00e9quente lors de l&rsquo;importation de fichiers textuels ou CSV comportant des marqueurs d&rsquo;absence non conventionnels \u2014 provoquera une rupture imm\u00e9diate de l&rsquo;ex\u00e9cution avec un message d&rsquo;erreur explicite signalant l&rsquo;impossibilit\u00e9 d&rsquo;appliquer des op\u00e9rateurs binaires \u00e0 des donn\u00e9es non num\u00e9riques.<\/p>\n<p>Enfin, les risques de masquage d&rsquo;espace de noms (<em>namespace masking<\/em>) ne doivent pas \u00eatre sous-estim\u00e9s. Si un environnement de travail charge simultan\u00e9ment plusieurs biblioth\u00e8ques disposant de fonctions homonymes, l&rsquo;appel non qualifi\u00e9 <code>mae()<\/code> peut pointer vers une m\u00e9thode concurrente au comportement inattendu. Il demeure par cons\u00e9quent une bonne pratique d&rsquo;encapsulation de pr\u00e9fixer formellement l&rsquo;instruction au moyen de l&rsquo;op\u00e9rateur de double deux-points : <code>Metrics::mae(actual = reel, predicted = pred)<\/code>.<\/p>\n<h2>7. Calcul de la MAE appliqu\u00e9e aux mod\u00e8les de r\u00e9gression lin\u00e9aire sous R<\/h2>\n<h3>7.1 Ajustement d&rsquo;un mod\u00e8le lin\u00e9aire standard via la fonction lm()<\/h3>\n<p>Dans la pratique statistique quotidienne, l&rsquo;\u00e9valuation de la MAE ne s&rsquo;effectue pas sur des vecteurs synth\u00e9tiques isol\u00e9s, mais intervient en aval de la mod\u00e9lisation formelle de relations structurelles entre une variable d\u00e9pendante continue et une constellation de r\u00e9gresseurs. L&rsquo;outil fondamental de l&rsquo;analyse lin\u00e9aire sous R est incarn\u00e9 par la fonction primitive <a rel=\"noopener noreferrer\" href=\"https:\/\/stat.ethz.ch\/R-manual\/R-patched\/library\/stats\/html\/lm.html\"><code>lm()<\/code><\/a>, qui estime les coefficients optimaux par la m\u00e9thode des moindres carr\u00e9s ordinaires.<\/p>\n<p>Consid\u00e9rons l&rsquo;ajustement d&rsquo;un mod\u00e8le de r\u00e9gression multivari\u00e9e visant \u00e0 mod\u00e9liser la performance cognitive \u00e0 partir de pr\u00e9dicteurs tels que la qualit\u00e9 du sommeil, l&rsquo;\u00e2ge chronologique et l&rsquo;indice de r\u00e9serve cognitive. La syntaxe standard s&rsquo;articule comme suit :<\/p>\n<p><code>modele_lineaire &lt;- lm(score_cognitif ~ heures_sommeil + age + reserve_cognitive, data = donnees_etude)<\/code><\/p>\n<p>Une fois le mod\u00e8le converg\u00e9 et les coefficients stabilis\u00e9s, la g\u00e9n\u00e9ration des estimations mod\u00e9lis\u00e9es sur l&rsquo;\u00e9chantillon de calibration s&rsquo;op\u00e8re au moyen de la m\u00e9thode g\u00e9n\u00e9rique <code>predict()<\/code> :<\/p>\n<p><code>predictions_calib &lt;- predict(modele_lineaire, newdata = donnees_etude)<\/code><\/p>\n<p>L&rsquo;extraction des valeurs observ\u00e9es r\u00e9elles depuis le tableau de donn\u00e9es originel permet de calculer l&rsquo;erreur absolue moyenne d&rsquo;ajustement interne (<em>in-sample MAE<\/em>) :<\/p>\n<p><code>mae_ajustement &lt;- Metrics::mae(donnees_etude$score_cognitif, predictions_calib)<\/code><\/p>\n<p>Cet indicateur quantifie la distance absolue moyenne r\u00e9siduelle entre les observations des sujets ayant servi \u00e0 b\u00e2tir l&rsquo;\u00e9quation et la droite (ou l&rsquo;hyperplan) de r\u00e9gression calibr\u00e9e.<\/p>\n<h3>7.2 \u00c9valuation de la g\u00e9n\u00e9ralisation : partitionnement Apprentissage et Test<\/h3>\n<p>L&rsquo;erreur absolue moyenne calcul\u00e9e sur l&rsquo;\u00e9chantillon d&rsquo;ajustement interne (<em>in-sample error<\/em>) offre une vision syst\u00e9matiquement optimiste des capacit\u00e9s r\u00e9elles du mod\u00e8le, en raison du ph\u00e9nom\u00e8ne de sur-apprentissage (<em>overfitting<\/em>). Les param\u00e8tres estim\u00e9s absorbent non seulement le signal structurel de la population, mais \u00e9galement le bruit stochastique propre \u00e0 l&rsquo;\u00e9chantillon de recueil. L&rsquo;appr\u00e9ciation rigoureuse du pouvoir de g\u00e9n\u00e9ralisation requiert la scission des donn\u00e9es en deux sous-ensembles ind\u00e9pendants : une partition d&rsquo;apprentissage (<em>training set<\/em>) et une partition de test herm\u00e9tique (<em>testing set<\/em>).<\/p>\n<p>Ce protocole de partitionnement al\u00e9atoire (couramment articul\u00e9 autour d&rsquo;un ratio 80\/20 ou 70\/30) se programme avec pr\u00e9cision en R base :<\/p>\n<p><code>set.seed(42)<br \/>\nindices_train &lt;- sample(seq_len(nrow(donnees_etude)), size = floor(0.80 * nrow(donnees_etude)))<br \/>\ntrain_set &lt;- donnees_etude[indices_train, ]<br \/>\ntest_set  &lt;- donnees_etude[-indices_train, ]<\/code><\/p>\n<p>Le mod\u00e8le est alors ajust\u00e9 exclusivement sur la partition <code>train_set<\/code>. Les coefficients ainsi fix\u00e9s sont ensuite appliqu\u00e9s aux donn\u00e9es in\u00e9dites de la partition <code>test_set<\/code> via l&rsquo;argument <code>newdata<\/code> de la fonction de pr\u00e9diction :<\/p>\n<p><code>modele_split &lt;- lm(score_cognitif ~ heures_sommeil + age + reserve_cognitive, data = train_set)<br \/>\npred_test &lt;- predict(modele_split, newdata = test_set)<br \/>\nmae_generalisation &lt;- Metrics::mae(test_set$score_cognitif, pred_test)<\/code><\/p>\n<p>La comparaison diff\u00e9rentielle entre la MAE d&rsquo;apprentissage et la MAE de g\u00e9n\u00e9ralisation constitue un indicateur diagnostique fondamental : une discordance mod\u00e9r\u00e9e t\u00e9moigne d&rsquo;un mod\u00e8le robuste et stable, tandis qu&rsquo;une explosion de la MAE sur l&rsquo;ensemble de test r\u00e9v\u00e8le un sur-ajustement pathologique de l&rsquo;\u00e9quation pr\u00e9dictive.<\/p>\n<h3>7.3 Diagnostic des r\u00e9sidus et corr\u00e9lation avec la m\u00e9trique MAE<\/h3>\n<p>L&rsquo;exploration approfondie des r\u00e9sidus d&rsquo;ajustement apporte un \u00e9clairage indispensable \u00e0 la bonne compr\u00e9hension de la MAE calcul\u00e9e. En vertu des propri\u00e9t\u00e9s des mod\u00e8les lin\u00e9aires classiques sous R, les r\u00e9sidus bruts d&rsquo;ajustement peuvent \u00eatre directement extraits de l&rsquo;objet de classe <code>lm<\/code> au moyen de la m\u00e9thode d\u00e9di\u00e9e <code>residuals(modele_lineaire)<\/code>.<\/p>\n<p>Il existe une \u00e9quivalence analytique parfaite, sur l&rsquo;ensemble d&rsquo;apprentissage, entre la MAE et l&rsquo;esp\u00e9rance arithm\u00e9tique de la valeur absolue des r\u00e9sidus du mod\u00e8le :<\/p>\n<p><code>mean(abs(residuals(modele_lineaire))) == Metrics::mae(donnees_etude$score_cognitif, fitted(modele_lineaire))<\/code><\/p>\n<p>Cette relation permet d&rsquo;\u00e9valuer directement l&rsquo;impact de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 \u2014 c&rsquo;est-\u00e0-dire la non-constance de la variance r\u00e9siduelle \u00e0 travers les diff\u00e9rentes valeurs pr\u00e9dites \u2014 sur la dispersion spatiale de l&rsquo;erreur absolue. Si la variance des erreurs s&rsquo;accro\u00eet pour les scores extr\u00eames de la variable cible, la MAE moyenne globale aura tendance \u00e0 sous-estimer la marge d&rsquo;incertitude pesant sur les hauts scores tout en la surestimant pour les valeurs m\u00e9dianes. L&rsquo;analyse des quantiles conditionnels des valeurs absolues de r\u00e9sidus permet alors de rep\u00e9rer avec pr\u00e9cision les sous-groupes de sujets ou d&rsquo;unit\u00e9s exp\u00e9rimentales qui concentrent les magnitudes d&rsquo;erreurs pr\u00e9dictives les plus pr\u00e9occupantes.<\/p>\n<h2>8. Applications de la MAE en psychologie quantitative et sciences du comportement<\/h2>\n<h3>8.1 Mod\u00e9lisation des \u00e9chelles psychom\u00e9triques et inventaires de personnalit\u00e9<\/h3>\n<p>Dans le domaine de la psychologie quantitative et de l&rsquo;\u00e9valuation psychom\u00e9trique, la quantification de la pr\u00e9cision pr\u00e9dictive rev\u00eat des implications directes sur la validit\u00e9 du diagnostic individuel. Les scores obtenus sur des inventaires cliniques valid\u00e9s \u2014 \u00e0 l&rsquo;instar de l&rsquo;inventaire de d\u00e9pression de Beck (BDI) ou de l&rsquo;inventaire de personnalit\u00e9 NEO-PI-R \u2014 constituent des variables continues born\u00e9es, r\u00e9sultant de l&rsquo;agr\u00e9gation d&rsquo;\u00e9chelles de Likert ordinales.<\/p>\n<p>L&rsquo;utilisation de l&rsquo;erreur absolue moyenne s&rsquo;av\u00e8re particuli\u00e8rement appropri\u00e9e dans ce champ disciplinaire en raison de la notion clinique de \u00ab\u00a0diff\u00e9rence minimale cliniquement significative\u00a0\u00bb (MCID, pour <em>Minimal Clinically Important Difference<\/em>). Si la litt\u00e9rature psychom\u00e9trique \u00e9tablit qu&rsquo;une variation de 3 points sur le score BDI refl\u00e8te un changement perceptible de l&rsquo;\u00e9tat clinique d&rsquo;un patient, une MAE de 2,1 points obtenue par un algorithme pr\u00e9dictif fond\u00e9 sur des biomarqueurs vocaux ou num\u00e9riques indique imm\u00e9diatement au clinicien que l&rsquo;erreur typique du mod\u00e8le demeure circonscrite sous le seuil d&rsquo;alt\u00e9ration clinique observable.<\/p>\n<p>En outre, l&rsquo;adoption de la MAE pr\u00e9serve l&rsquo;\u00e9thique de l&rsquo;\u00e9valuation psychologique. L&rsquo;usage de la RMSE tend \u00e0 disqualifier des mod\u00e8les pr\u00e9dictifs qui \u00e9chouent lourdement sur une poign\u00e9e de profils psychologiques atypiques ou de r\u00e9ponses extr\u00eames non pathologiques. La lin\u00e9arit\u00e9 de la MAE \u00e9vite cette surpond\u00e9ration excessive des scores aberrants, permettant de valoriser la pr\u00e9cision r\u00e9guli\u00e8re et stable de l&rsquo;outil pour la grande majorit\u00e9 de la patient\u00e8le.<\/p>\n<h3>8.2 Temps de r\u00e9action et chronom\u00e9trie cognitive<\/h3>\n<p>L&rsquo;\u00e9tude des temps de r\u00e9action (TR) dans les t\u00e2ches de chronom\u00e9trie cognitive (telles que les paradigmes de Stroop, de Simon ou de d\u00e9tection de signaux) confronte syst\u00e9matiquement le statisticien \u00e0 des distributions de donn\u00e9es fortement asym\u00e9triques, traditionnellement mod\u00e9lis\u00e9es par des lois ex-gaussiennes ou des distributions de Wald. Ces structures temporelles se singularisent par une queue \u00e9tir\u00e9e vers les latences \u00e9lev\u00e9es, mat\u00e9rialisant des d\u00e9crochages attentionnels passagers de la part des participants.<\/p>\n<p>Dans ce contexte exp\u00e9rimental, l&rsquo;\u00e9valuation des architectures computationnelles de la prise de d\u00e9cision \u2014 tels que les mod\u00e8les de diffusion neuronale (DDM pour <em>Drift-Diffusion Models<\/em>) \u2014 par des m\u00e9triques d&rsquo;erreur quadratiques comme la RMSE aboutit fr\u00e9quemment \u00e0 des conclusions statistiques erron\u00e9es. La p\u00e9nalisation quadratique contraint les algorithmes d&rsquo;optimisation \u00e0 d\u00e9placer artificiellement leur estimation centrale vers la queue de distribution afin d&rsquo;amortir l&rsquo;impact de quelques temps de r\u00e9action exceptionnellement longs (par exemple, un temps de r\u00e9ponse de 2500 ms dans une t\u00e2che calibr\u00e9e sur une m\u00e9diane de 450 ms).<\/p>\n<p>La MAE, en mesurant l&rsquo;\u00e9cart absolu moyen sans amplification quadratique, offre une \u00e9valuation remarquablement fid\u00e8le de l&rsquo;ad\u00e9quation du mod\u00e8le cognitif \u00e0 la zone de densit\u00e9 d\u00e9cisionnelle principale. Elle assure une mod\u00e9lisation stable et reproductible des constantes de temps cognitives, autorisant des comparaisons inter-sujets fiables m\u00eame au sein d&rsquo;exp\u00e9riences comportementales administr\u00e9es en ligne, soumises \u00e0 d&rsquo;importantes fluctuations environnementales non contr\u00f4l\u00e9es.<\/p>\n<h3>8.3 Mod\u00e9lisation longitudinale de comportements th\u00e9rapeutiques<\/h3>\n<p>La recherche sur l&rsquo;efficacit\u00e9 des interventions psychoth\u00e9rapeutiques et comportementales s&rsquo;appuie d\u00e9sormais massivement sur des devis longitudinaux intensifs, mat\u00e9rialis\u00e9s par les m\u00e9thodes d&rsquo;\u00e9chantillonnage des exp\u00e9riences (ESM pour <em>Experience Sampling Methodology<\/em>) ou d&rsquo;\u00e9valuation \u00e9cologique instantan\u00e9e (EMA). Ces protocoles recueillent des mesures r\u00e9p\u00e9t\u00e9es de l&rsquo;affect, de l&rsquo;alliance th\u00e9rapeutique et de l&rsquo;adh\u00e9sion au traitement plusieurs fois par jour durant des semaines.<\/p>\n<p>L&rsquo;application de mod\u00e8les lin\u00e9aires mixtes (multiniveaux) \u00e0 ces s\u00e9ries chronologiques denses permet de d\u00e9composer l&rsquo;erreur absolue moyenne \u00e0 deux \u00e9chelons analytiques distincts : la MAE intra-individuelle et la MAE inter-individuelle. Le calcul de la MAE intra-individuelle renseigne le chercheur sur la capacit\u00e9 du mod\u00e8le \u00e0 anticiper avec pr\u00e9cision les oscillations dynamiques quotidiennes d&rsquo;un patient donn\u00e9, tandis que la MAE inter-individuelle \u00e9value l&rsquo;exactitude de la pr\u00e9diction du profil symptomatique moyen \u00e0 travers l&rsquo;\u00e9chantillon global.<\/p>\n<p>Cette granularit\u00e9 analytique trouve une r\u00e9sonance directe dans le d\u00e9veloppement d&rsquo;algorithmes de m\u00e9decine personnalis\u00e9e et de psychiatrie de pr\u00e9cision. En exprimant la d\u00e9rive pr\u00e9dictive en unit\u00e9s r\u00e9elles de l&rsquo;\u00e9chelle d&rsquo;auto-\u00e9valuation (par exemple, des unit\u00e9s de s\u00e9v\u00e9rit\u00e9 sur une \u00e9chelle visuelle analogique gradu\u00e9e de 0 \u00e0 100), l&rsquo;\u00e9quipe clinique dispose d&rsquo;une m\u00e9trique d&rsquo;\u00e9valuation directement traduisible en protocoles d&rsquo;alerte algorithmique et d&rsquo;ajustement th\u00e9rapeutique en temps r\u00e9el.<\/p>\n<h2>9. Automatisation et calcul avanc\u00e9 avec caret et l&rsquo;\u00e9cosyst\u00e8me tidymodels<\/h2>\n<h3>9.1 Validation crois\u00e9e k-fold orchestr\u00e9e avec caret<\/h3>\n<p>L&rsquo;automatisation des proc\u00e9dures d&rsquo;\u00e9valuation statistique \u00e9vite les \u00e9cueils de variabilit\u00e9 inh\u00e9rents \u00e0 un d\u00e9coupage unique en ensembles d&rsquo;apprentissage et de test. Le package <code>caret<\/code> offre une structure unifi\u00e9e permettant d&rsquo;orchestrer des protocoles rigoureux de validation crois\u00e9e \u00e0 <em>k<\/em> blocs (<em>k-fold cross-validation<\/em>) r\u00e9p\u00e9t\u00e9s, tout en imposant la MAE comme crit\u00e8re r\u00e9gulateur du tuning algorithmique.<\/p>\n<p>La configuration de la strat\u00e9gie d&rsquo;\u00e9chantillonnage s&rsquo;op\u00e8re via la fonction <code>trainControl()<\/code>, au sein de laquelle la m\u00e9thode de validation crois\u00e9e est stipul\u00e9e :<\/p>\n<p><code>controle_cv &lt;- caret::trainControl(<br \/>\n&nbsp;&nbsp;method = \"repeatedcv\",<br \/>\n&nbsp;&nbsp;number = 10,<br \/>\n&nbsp;&nbsp;repeats = 5,<br \/>\n&nbsp;&nbsp;verboseIter = FALSE<br \/>\n)<\/code><\/p>\n<p>L&rsquo;entra\u00eenement du mod\u00e8le pr\u00e9dictif et l&rsquo;exploration des hyperparam\u00e8tres sont ensuite d\u00e9clench\u00e9s par l&rsquo;appel de la fonction principale <code>train()<\/code>. L&rsquo;argument explicite <code>metric = \"MAE\"<\/code> ordonne au moteur de recherche de s\u00e9lectionner la combinaison d&rsquo;hyperparam\u00e8tres qui minimise strictement l&rsquo;erreur absolue moyenne \u00e0 travers l&rsquo;ensemble des plis de validation :<\/p>\n<p><code>modele_optimise &lt;- caret::train(<br \/>\n&nbsp;&nbsp;score_clinique ~ .,<br \/>\n&nbsp;&nbsp;data = donnees_psychometrie,<br \/>\n&nbsp;&nbsp;method = \"glmnet\",<br \/>\n&nbsp;&nbsp;trControl = controle_cv,<br \/>\n&nbsp;&nbsp;metric = \"MAE\"<br \/>\n)<\/code><\/p>\n<p>L&rsquo;extraction de l&rsquo;objet final via <code>modele_optimise$results<\/code> permet de consulter la moyenne et l&rsquo;\u00e9cart-type de la MAE calcul\u00e9e \u00e0 travers les 50 plis de validation crois\u00e9e (10 blocs multipli\u00e9s par 5 r\u00e9p\u00e9titions), garantissant une estimation asymptotiquement non biais\u00e9e du pouvoir g\u00e9n\u00e9ralisateur de l&rsquo;architecture pr\u00e9dictive retenue.<\/p>\n<h3>9.2 Mise en \u0153uvre moderne avec le package yardstick de tidymodels<\/h3>\n<p>L&rsquo;int\u00e9gration de la mod\u00e9lisation statistique au sein de l&rsquo;\u00e9cosyst\u00e8me contemporain <code>tidymodels<\/code> consacre le package <a rel=\"noopener noreferrer\" href=\"https:\/\/yardstick.tidymodels.org\/\"><code>yardstick<\/code><\/a> comme la r\u00e9f\u00e9rence pour l&rsquo;\u00e9valuation empirique des performances. S&rsquo;inscrivant dans la philosophie du <em>tidy data<\/em>, <code>yardstick<\/code> traite les sorties pr\u00e9dictives sous la forme de structures <code>tibble<\/code> pr\u00e9servant l&rsquo;alignement strict entre v\u00e9rit\u00e9 terrain et estimations.<\/p>\n<p>Consid\u00e9rons un tableau de pr\u00e9dictions <code>resultats_tibble<\/code> contenant la colonne r\u00e9elle <code>anxiete_observee<\/code> et la colonne estim\u00e9e <code>anxiete_predite<\/code>. L&rsquo;\u00e9valuation de l&rsquo;erreur absolue moyenne s&rsquo;effectue via une syntaxe standardis\u00e9e et \u00e9l\u00e9gante :<\/p>\n<p><code>metrique_mae &lt;- yardstick::mae(<br \/>\n&nbsp;&nbsp;data = resultats_tibble,<br \/>\n&nbsp;&nbsp;truth = anxiete_observee,<br \/>\n&nbsp;&nbsp;estimate = anxiete_predite<br \/>\n)<\/code><\/p>\n<p>Cette fonction retourne une table format\u00e9e r\u00e9pertoriant le nom de la m\u00e9trique (<code>.metric<\/code>), l&rsquo;estimateur employ\u00e9 (<code>.estimator<\/code>, g\u00e9n\u00e9ralement <code>standard<\/code>) et la valeur calcul\u00e9e (<code>.estimate<\/code>). L&rsquo;un des avantages conceptuels majeurs de <code>yardstick<\/code> r\u00e9side dans sa capacit\u00e9 \u00e0 construire des ensembles m\u00e9triques composites via <code>metric_set()<\/code> :<\/p>\n<p><code>ensemble_metriques &lt;- yardstick::metric_set(yardstick::mae, yardstick::rmse, yardstick::rsq)<br \/>\nevaluation_globale &lt;- ensemble_metriques(resultats_tibble, truth = anxiete_observee, estimate = anxiete_predite)<\/code><\/p>\n<p>De surcro\u00eet, la compatibilit\u00e9 native avec l&rsquo;op\u00e9rateur de groupement <code>dplyr::group_by()<\/code> permet d&rsquo;\u00e9valuer instantan\u00e9ment la MAE \u00e0 travers diff\u00e9rentes strates de participants (par exemple, selon le genre, la cohorte d&rsquo;\u00e2ge ou le site clinique de collecte) en une unique instruction vectoris\u00e9e, sans n\u00e9cessiter la moindre boucle explicite :<\/p>\n<p><code>mae_par_cohorte &lt;- resultats_tibble %&gt;%<br \/>\n&nbsp;&nbsp;dplyr::group_by(cohorte_clinique) %&gt;%<br \/>\n&nbsp;&nbsp;yardstick::mae(truth = anxiete_observee, estimate = anxiete_predite)<\/code><\/p>\n<h3>9.3 Calcul de la MAE pour les algorithmes non lin\u00e9aires (Random Forest, SVM)<\/h3>\n<p>L&rsquo;\u00e9valuation comparative de la performance statistique atteint sa pleine complexit\u00e9 lors de la confrontation d&rsquo;approches non lin\u00e9aires hautement flexibles, telles que les for\u00eats al\u00e9atoires (<em>Random Forest<\/em>) ou les machines \u00e0 vecteurs de support (<em>Support Vector Machines<\/em> &#8211; SVM), sur des donn\u00e9es multidimensionnelles.<\/p>\n<p>Ces architectures non param\u00e9triques disposent d&rsquo;une capacit\u00e9 d&rsquo;approximation universelle, susceptible de distordre l&rsquo;erreur r\u00e9siduelle de mani\u00e8re disparate \u00e0 travers l&rsquo;espace des pr\u00e9dicteurs. Par exemple, une for\u00eat al\u00e9atoire optimis\u00e9e par d\u00e9faut sous une perte quadratique tend \u00e0 lisser les pr\u00e9dictions vers la moyenne locale de chaque feuille terminale. L&rsquo;application terminale de la MAE sur les donn\u00e9es de test ind\u00e9pendantes constitue alors le juge de paix impartial :<\/p>\n<p><code>pred_rf &lt;- predict(modele_randomforest, newdata = test_data)<br \/>\npred_svm &lt;- predict(modele_svm, newdata = test_data)<\/p>\n<p>mae_rf  &lt;- Metrics::mae(test_data$score, pred_rf)<br \/>\nmae_svm &lt;- Metrics::mae(test_data$score, pred_svm)<\/code><\/p>\n<p>Cette confrontation directe permet de d\u00e9terminer si la complexit\u00e9 structurelle additionnelle introduite par des algorithmes \u00e0 base de noyaux (SVM) ou d&rsquo;assemblage d&rsquo;arbres d\u00e9cisionnels compense r\u00e9ellement le manque d&rsquo;interpr\u00e9tabilit\u00e9 des coefficients par une r\u00e9duction cliniquement substantielle de l&rsquo;erreur absolue moyenne globale, guidant le choix ultime du mod\u00e8le destin\u00e9 au d\u00e9ploiement en recherche translationnelle.<\/p>\n<h2>10. Gestion des valeurs aberrantes et situations limites de donn\u00e9es<\/h2>\n<h3>10.1 Influence des valeurs extr\u00eames sur la MAE comparativement aux m\u00e9triques quadratiques<\/h3>\n<p>La sensibilit\u00e9 diff\u00e9rentielle des m\u00e9triques statistiques face aux donn\u00e9es extr\u00eames (<em>outliers<\/em>) constitue l&rsquo;un des crit\u00e8res de s\u00e9lection les plus cruciaux dans la conception d&rsquo;un pipeline analytique. Afin de formaliser empiriquement cette dynamique, consid\u00e9rons une simulation num\u00e9rique o\u00f9 l&rsquo;on injecte d\u00e9lib\u00e9r\u00e9ment une contamination de magnitude croissante au sein d&rsquo;une seule observation d&rsquo;un \u00e9chantillon psychologique de taille <em>N<\/em> = 100.<\/p>\n<p>Soit un vecteur r\u00e9siduel homog\u00e8ne o\u00f9 l&rsquo;ensemble des d\u00e9viations individuelles oscille initialement autour de 2 unit\u00e9s, produisant une MAE de 2,0 et une RMSE de 2,0. Si une unique observation enregistre accidentellement une divergence erratique de 100 unit\u00e9s (induite par une d\u00e9faillance de saisie informatique ou un choc physiologique transitoire), l&rsquo;impact marginal sur les deux m\u00e9triques diverge radicalement :<\/p>\n<ul>\n<li>La <strong>MAE<\/strong> augmente de la valeur scalaire (100 &minus; 2) \/ 100 = 0,98 unit\u00e9, s&rsquo;\u00e9tablissant \u00e0 environ 2,98. La d\u00e9gradation de la performance est proportionnelle \u00e0 la taille d&rsquo;\u00e9chantillon.<\/li>\n<li>La <strong>RMSE<\/strong> int\u00e8gre la variation quadratique [(99 * 2<sup>2<\/sup> + 100<sup>2<\/sup>) \/ 100]<sup>1\/2<\/sup> &approx; [10396 \/ 100]<sup>1\/2<\/sup> &approx; 10,20 unit\u00e9s. L&rsquo;indicateur quadratique est multipli\u00e9 par plus de 5, sous l&rsquo;influence exclusive d&rsquo;une seule mesure discordante.<\/li>\n<\/ul>\n<p>Cet exemple formalise le concept de robustesse en th\u00e9orie de l&rsquo;estimation. L&rsquo;estimateur L1 poss\u00e8de une r\u00e9sistance intrins\u00e8que qui \u00e9pargne au chercheur l&rsquo;obligation d&rsquo;\u00e9liminer cavali\u00e8rement des donn\u00e9es atypiques authentiques. Alors que la RMSE force fr\u00e9quemment l&rsquo;analyste \u00e0 un rejet d&rsquo;outliers discutable pour restaurer la stabilit\u00e9 apparente de son mod\u00e8le, la MAE absorbe la perturbation avec une mod\u00e9ration math\u00e9matique exemplaire.<\/p>\n<h3>10.2 Gestion des valeurs manquantes, infinies et biais d&rsquo;\u00e9chantillonnage<\/h3>\n<p>Les situations limites de donn\u00e9es incluent fr\u00e9quemment l&rsquo;\u00e9mergence inattendue de valeurs non finies (<code>Inf<\/code>, <code>-Inf<\/code>) ou de formes ind\u00e9termin\u00e9es (<code>NaN<\/code> pour <em>Not a Number<\/em>). Ces entit\u00e9s math\u00e9matiques apparaissent typiquement lors de phases de pr\u00e9-traitement impliquant des transformations non lin\u00e9aires (telles que le logarithme de z\u00e9ro ou la division par une variance nulle).<\/p>\n<p>Lors de l&rsquo;\u00e9valuation de la MAE, la pr\u00e9sence d&rsquo;une seule valeur <code>NaN<\/code> ou <code>Inf<\/code> au sein du vecteur d&rsquo;estimations mod\u00e9lis\u00e9es corrompt irr\u00e9m\u00e9diablement le calcul scalaire en renvoyant une divergence ind\u00e9termin\u00e9e. Pour assainir rigoureusement les pipelines statistiques, un filtrage formel en amont au moyen de fonctions R base telles que <code>complete.cases()<\/code> ou de la routine <code>tidyr::drop_na()<\/code> est formellement requis :<\/p>\n<p><code>indices_valides &lt;- is.finite(valeurs_observees) &amp; is.finite(valeurs_predites)<br \/>\nmae_nettoyee &lt;- Metrics::mae(valeurs_observees[indices_valides], valeurs_predites[indices_valides])<\/code><\/p>\n<p>Par ailleurs, les donn\u00e9es issues d&rsquo;\u00e9chelles psychologiques standardis\u00e9es se heurtent couramment \u00e0 des effets de plancher (<em>floor effects<\/em>) ou de plafond (<em>ceiling effects<\/em>). Ces troncatures structurelles induisent une distorsion s\u00e9v\u00e8re : le mod\u00e8le pr\u00e9dit des valeurs au-del\u00e0 des bornes licites de l&rsquo;instrument (par exemple un score pr\u00e9dit de 85 sur une \u00e9chelle plafonn\u00e9e \u00e0 80). Dans ce cadre, la MAE brute peut refl\u00e9ter une erreur algorithmique artificielle li\u00e9e \u00e0 l&rsquo;absence de contraintes aux limites. Une correction m\u00e9thodologique consiste soit \u00e0 appliquer un seuillage (<em>clipping<\/em>) des pr\u00e9dictions aux bornes de l&rsquo;instrument avant le calcul de la m\u00e9trique, soit \u00e0 mod\u00e9liser le ph\u00e9nom\u00e8ne via des r\u00e9gressions censur\u00e9es (mod\u00e8les Tobit) adapt\u00e9es \u00e0 l&rsquo;\u00e9valuation des erreurs absolues sous troncature.<\/p>\n<h3>10.3 Invariance d&rsquo;\u00e9chelle et transformations de variables<\/h3>\n<p>Dans de nombreuses applications d&rsquo;inf\u00e9rence statistique, les variables cibles asym\u00e9triques subissent des transformations stabilisatrices de variance, au premier rang desquelles figurent la transformation logarithmique <em>log(y)<\/em> ou la transformation param\u00e9trique de <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Box%E2%80%93Cox_transformation\">Box-Cox<\/a>. Le calcul de la MAE au sein de l&rsquo;espace transform\u00e9 introduit toutefois un biais d&rsquo;interpr\u00e9tation dimensionnelle majeur.<\/p>\n<p>Si un mod\u00e8le pr\u00e9dictif est calibr\u00e9 sur une \u00e9chelle log-transform\u00e9e <em>ln(Y)<\/em>, l&rsquo;\u00e9valuation de la m\u00e9trique <code>mean(abs(log(y) - log(y_pred)))<\/code> quantifie en r\u00e9alit\u00e9 l&rsquo;erreur absolue moyenne des logarithmes, et non l&rsquo;erreur sur l&rsquo;\u00e9chelle physiologique ou comportementale d&rsquo;origine. Cette grandeur dans l&rsquo;espace transform\u00e9 poss\u00e8de une signification relative (apparent\u00e9e \u00e0 un ratio multiplicatif moyen d&rsquo;erreur), mais perd toute correspondance directe avec les unit\u00e9s physiques initiales.<\/p>\n<p>Pour r\u00e9tablir l&rsquo;int\u00e9grit\u00e9 de la communication scientifique, l&rsquo;analyste doit imp\u00e9rativement r\u00e9-exprimer les pr\u00e9dictions dans l&rsquo;espace m\u00e9trique originel via la fonction inverse correspondante (l&rsquo;exponentielle <code>exp()<\/code> pour le logarithme) avant d&rsquo;ex\u00e9cuter l&rsquo;\u00e9valuation de la MAE. Il convient cependant de garder \u00e0 l&rsquo;esprit un th\u00e9or\u00e8me fondamental de la th\u00e9orie statistique : l&rsquo;inversion non lin\u00e9aire d&rsquo;une esp\u00e9rance conditionnelle n&rsquo;est pas \u00e9gale \u00e0 l&rsquo;esp\u00e9rance conditionnelle inverse (en vertu de l&rsquo;in\u00e9galit\u00e9 de Jensen). N\u00e9anmoins, dans le cas strict de la MAE \u2014 qui optimise la m\u00e9diane conditionnelle et non la moyenne \u2014, l&rsquo;application d&rsquo;une fonction strictement monotone poss\u00e8de la propri\u00e9t\u00e9 remarquable de commuter parfaitement avec l&rsquo;op\u00e9rateur m\u00e9diane :<\/p>\n<p><strong>Med[f(Y)|X] = f(Med[Y|X])<\/strong><\/p>\n<p>Cette invariance de transformation monotone constitue un argument \u00e9pist\u00e9mologique et pratique massif en faveur de l&rsquo;\u00e9valuation par la MAE lors de l&rsquo;usage de donn\u00e9es transform\u00e9es.<\/p>\n<h2>11. Visualisation graphique de la MAE et analyse diagnostique des r\u00e9sidus<\/h2>\n<h3>11.1 Repr\u00e9sentation cart\u00e9sienne des \u00e9carts absolus avec ggplot2<\/h3>\n<p>L&rsquo;exploration visuelle constitue l&rsquo;adjuvant indispensable de toute analyse num\u00e9rique quantitative. Au sein de l&rsquo;environnement R, le package <a rel=\"noopener noreferrer\" href=\"https:\/\/ggplot2.tidyverse.org\/\"><code>ggplot2<\/code><\/a> offre un cadre grammatical puissant pour mat\u00e9rialiser spatialement l&rsquo;erreur absolue moyenne et diagnostiquer les distorsions pr\u00e9dictives locales. La repr\u00e9sentation cart\u00e9sienne de r\u00e9f\u00e9rence repose sur la confrontation directe des observations r\u00e9elles en abscisse et des estimations mod\u00e9lis\u00e9es en ordonn\u00e9e.<\/p>\n<p>Le fondement de ce trac\u00e9 diagnostique r\u00e9side dans l&rsquo;adjonction de la bissectrice parfaite (la droite identit\u00e9 <em>y = x<\/em>, caract\u00e9ris\u00e9e par une pente de 1 et une ordonn\u00e9e \u00e0 l&rsquo;origine de 0). Tout point s&rsquo;\u00e9cartant de cette ligne diagonale mat\u00e9rialise une d\u00e9faillance pr\u00e9dictive. L&rsquo;originalit\u00e9 p\u00e9dagogique et analytique d&rsquo;une telle visualisation consiste \u00e0 tracer explicitement les segments de r\u00e9sidus verticaux reliant chaque point coordonn\u00e9 (<em>x<sub>i<\/sub>, y<sub>i<\/sub><\/em>) \u00e0 la droite d&rsquo;exactitude th\u00e9orique.<\/p>\n<p>Le script d\u00e9claratif sous <code>ggplot2<\/code> s&rsquo;articule comme suit :<\/p>\n<p><code>library(ggplot2)<\/p>\n<p>graphique_dispersion &lt;- ggplot(donnees_eval, aes(x = observe, y = predit)) +<br \/>\n&nbsp;&nbsp;geom_segment(aes(xend = observe, yend = observe), color = \"firebrick\", alpha = 0.4) +<br \/>\n&nbsp;&nbsp;geom_point(color = \"steelblue\", size = 2.5, alpha = 0.8) +<br \/>\n&nbsp;&nbsp;geom_abline(slope = 1, intercept = 0, linetype = \"dashed\", color = \"black\", linewidth = 0.8) +<br \/>\n&nbsp;&nbsp;annotate(\"text\", x = min(donnees_eval<span class=\"arabpsy-math-inline\" data-latex=\"observe) * 1.1, y = \\max(donnees_eval\">$observe) * 1.1, y = \\max(donnees_eval$<\/span>predit) * 0.9,<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;label = paste(\"MAE =\", round(Metrics::mae(donnees_eval<span class=\"arabpsy-math-inline\" data-latex=\"observe, donnees_eval\">$observe, donnees_eval$<\/span>predit), 3)),<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;size = 5, fontface = \"bold\") +<br \/>\n&nbsp;&nbsp;labs(title = \"Diagnostic de pr\u00e9cision : Observations versus Pr\u00e9dictions\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;subtitle = \"Les segments rouges repr\u00e9sentent la magnitude |yi - xi| pour chaque sujet\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x = \"Valeurs Observ\u00e9es (V\u00e9rit\u00e9 Terrain)\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y = \"Valeurs Pr\u00e9dites par le Mod\u00e8le\") +<br \/>\n&nbsp;&nbsp;theme_minimal()<\/code><\/p>\n<p>La pr\u00e9sence d&rsquo;une zone ombr\u00e9e mat\u00e9rialisant un ruban d&rsquo;incertitude de largeur \u00e9gale \u00e0 &plusmn; 1 MAE de part et d&rsquo;autre de la bissectrice permet de visualiser instantan\u00e9ment la proportion de la cohorte dont l&rsquo;estimation s&rsquo;inscrit au sein de la marge d&rsquo;erreur moyenne attendue.<\/p>\n<h3>11.2 Distribution empirique des r\u00e9sidus absolus<\/h3>\n<p>L&rsquo;inspection de la valeur scalaire de la MAE dissimule n\u00e9cessairement la forme g\u00e9om\u00e9trique de la distribution des d\u00e9viations. Deux jeux de donn\u00e9es peuvent g\u00e9n\u00e9rer une MAE rigoureusement identique de 3,0 unit\u00e9s tout en arborant des structures r\u00e9siduelles totalement discordantes : le premier pr\u00e9sentant une distribution uniforme des \u00e9carts entre 0 et 6, le second combinant 90 % d&rsquo;erreurs infimes de 0,5 unit\u00e9 et 10 % d&rsquo;erreurs massives de 25,5 unit\u00e9s.<\/p>\n<p>Pour lever cette ambigu\u00eft\u00e9 diagnostique, la mod\u00e9lisation de la densit\u00e9 empirique des r\u00e9sidus absolus s&rsquo;av\u00e8re indispensable. Le calcul de la variable d\u00e9riv\u00e9e <code>residu_absolu = abs(observe - predit)<\/code> permet de construire un histogramme de fr\u00e9quences surmont\u00e9 d&rsquo;une estimation non param\u00e9trique de densit\u00e9 par noyau (<em>kernel density estimation<\/em>) :<\/p>\n<p><code>donnees_eval<span class=\"arabpsy-math-inline\" data-latex=\"residu_\\abs &lt;- \\abs(donnees_eval\">$residu_\\abs &lt;- \\abs(donnees_eval$<\/span>observe - donnees_eval$predit)<\/p>\n<p>graphique_densite &lt;- ggplot(donnees_eval, aes(x = residu_abs)) +<br \/>\n&nbsp;&nbsp;geom_histogram(aes(y = after_stat(density)), bins = 30, fill = \"gray80\", color = \"white\") +<br \/>\n&nbsp;&nbsp;geom_density(color = \"darkblue\", linewidth = 1) +<br \/>\n&nbsp;&nbsp;geom_vline(aes(xintercept = mean(residu_abs)), color = \"firebrick\", linetype = \"solid\", linewidth = 1) +<br \/>\n&nbsp;&nbsp;geom_vline(aes(xintercept = median(residu_abs)), color = \"forestgreen\", linetype = \"dotted\", linewidth = 1) +<br \/>\n&nbsp;&nbsp;labs(title = \"Distribution empirique de l'amplitude des erreurs\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;subtitle = \"Ligne continue rouge = MAE ; Ligne pointill\u00e9e verte = M\u00e9diane de l'erreur absolue\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x = \"Erreur Absolue |yi - xi|\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y = \"Densit\u00e9 de probabilit\u00e9\") +<br \/>\n&nbsp;&nbsp;theme_light()<\/code><\/p>\n<p>L&rsquo;\u00e9cart spatial entre la ligne marquant la MAE (la moyenne arithm\u00e9tique de l&rsquo;\u00e9cart absolu) et la m\u00e9diane de l&rsquo;erreur absolue met imm\u00e9diatement en \u00e9vidence le degr\u00e9 d&rsquo;asym\u00e9trie positive de la distribution des d\u00e9faillances. Si la MAE s&rsquo;\u00e9carte substantiellement vers la droite par rapport \u00e0 la m\u00e9diane des erreurs, l&rsquo;analyste identifie imm\u00e9diatement l&rsquo;existence d&rsquo;une grappe de participants concentrant des erreurs atypiques majeures m\u00e9ritant une investigation clinique isol\u00e9e.<\/p>\n<h3>11.3 Graphiques de diagnostic multi-mod\u00e8les et comparaisons crois\u00e9es<\/h3>\n<p>Dans la phase conclusive de s\u00e9lection d&rsquo;une mod\u00e9lisation, le chercheur est invariablement amen\u00e9 \u00e0 comparer les distributions de MAE g\u00e9n\u00e9r\u00e9es par une pluralit\u00e9 d&rsquo;architectures concurrentes (par exemple : r\u00e9gression lin\u00e9aire multiple standard, r\u00e9gression r\u00e9gularis\u00e9e ElasticNet, mod\u00e8les additifs g\u00e9n\u00e9ralis\u00e9s GAM et for\u00eats al\u00e9atoires) soumises \u00e0 une validation crois\u00e9e r\u00e9p\u00e9t\u00e9e.<\/p>\n<p>La repr\u00e9sentation sous forme de diagrammes en bo\u00eete \u00e0 moustaches (<em>boxplots<\/em>) combin\u00e9s \u00e0 la superposition des points individuels de chaque pli (<em>jittered points<\/em>) offre une lecture diagnostique de premier ordre pour \u00e9valuer la stabilit\u00e9 pr\u00e9dictive :<\/p>\n<p><code>graphique_comparatif &lt;- ggplot(resultats_cv_modeles, aes(x = algorithme, y = mae_pli, fill = algorithme)) +<br \/>\n&nbsp;&nbsp;geom_boxplot(alpha = 0.5, outlier.shape = NA) +<br \/>\n&nbsp;&nbsp;geom_jitter(width = 0.2, alpha = 0.7, size = 1.8) +<br \/>\n&nbsp;&nbsp;stat_summary(fun = mean, geom = \"point\", shape = 18, size = 4, color = \"darkred\") +<br \/>\n&nbsp;&nbsp;labs(title = \"Comparaison de la MAE \u00e0 travers les plis de validation crois\u00e9e\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;subtitle = \"Les losanges rouges indiquent la MAE moyenne globale de chaque mod\u00e8le\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x = \"Architecture algorithmique\",<br \/>\n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y = \"Erreur Absolue Moyenne (MAE)\") +<br \/>\n&nbsp;&nbsp;theme_classic() +<br \/>\n&nbsp;&nbsp;theme(legend.position = \"none\")<\/code><\/p>\n<p>L&rsquo;exploitation conjointe des facettes graphiques via <code>facet_wrap(~ sous_groupe_clinique)<\/code> permet d&rsquo;\u00e9tendre cette comparaison \u00e0 travers diff\u00e9rentes strates diagnostiques ou d\u00e9mographiques, r\u00e9v\u00e9lant si la sup\u00e9riorit\u00e9 pr\u00e9dictive globale d&rsquo;un algorithme s&rsquo;exerce de mani\u00e8re homog\u00e8ne sur tous les profils de sujets ou si elle dissimule des d\u00e9gradations locales de performance sur des sous-populations sp\u00e9cifiques.<\/p>\n<h2>12. Rigueur m\u00e9thodologique, communication des r\u00e9sultats et bonnes pratiques<\/h2>\n<h3>12.1 Normes de reporting acad\u00e9mique et publications scientifiques<\/h3>\n<p>La communication des indices d&rsquo;ad\u00e9quation pr\u00e9dictive au sein des revues scientifiques internationales \u00e0 comit\u00e9 de lecture est r\u00e9gie par des standards m\u00e9thodologiques stricts, formalis\u00e9s notamment par l&rsquo;American Psychological Association (<a rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\/\">APA<\/a>) et les lignes directrices TRIPOD (pour les mod\u00e8les d&rsquo;apprentissage pr\u00e9dictifs en sant\u00e9). La pr\u00e9sentation isol\u00e9e d&rsquo;un coefficient de d\u00e9termination R<sup>2<\/sup> est aujourd&rsquo;hui consid\u00e9r\u00e9e comme notoirement insuffisante, ce dernier \u00e9tant insensible aux translations syst\u00e9matiques du biais directionnel.<\/p>\n<p>Le reporting scientifique rigoureux impose la communication conjointe de trois grandeurs m\u00e9triques fondamentales :<\/p>\n<ul>\n<li>Une mesure d&rsquo;ad\u00e9quation globale robuste et dimensionnellement pr\u00e9serv\u00e9e : la <strong>MAE<\/strong>.<\/li>\n<li>Une m\u00e9trique quadratique compl\u00e9mentaire : la <strong>RMSE<\/strong>, dont le rapport avec la MAE renseigne directement la communaut\u00e9 scientifique sur la variance de la dispersion des r\u00e9sidus.<\/li>\n<li>Le <strong>biais directionnel moyen<\/strong> (<em>Mean Error<\/em> ou ME = mean(pred &minus; obs)), indispensable pour attester de l&rsquo;absence de sur-estimation ou de sous-estimation structurelle globale de l&rsquo;outil mod\u00e9lis\u00e9.<\/li>\n<\/ul>\n<p>De surcro\u00eet, les normes acad\u00e9miques contemporaines exigent une transparence computationnelle totale. Le manuscrit soumis doit sp\u00e9cifier sans ambigu\u00eft\u00e9 : le syst\u00e8me d&rsquo;exploitation utilis\u00e9, la version de distribution du langage R, l&rsquo;int\u00e9gralit\u00e9 des packages d\u00e9ploy\u00e9s avec leurs num\u00e9ros de version respectifs, ainsi que les valeurs initiales des g\u00e9n\u00e9rateurs pseudo-al\u00e9atoires ayant permis d&rsquo;obtenir les r\u00e9sultats pr\u00e9sent\u00e9s. Les scripts R document\u00e9s doivent \u00eatre mis \u00e0 disposition dans des d\u00e9p\u00f4ts ouverts p\u00e9rennes (tels que l&rsquo;Open Science Framework ou Zenodo), garantissant la reproductibilit\u00e9 empirique sans restriction.<\/p>\n<h3>12.2 Calcul d&rsquo;intervalles de confiance pour la MAE par r\u00e9\u00e9chantillonnage bootstrap<\/h3>\n<p>Une erreur absolue moyenne issue d&rsquo;un ensemble de test ou d&rsquo;un \u00e9chantillon observationnel constitue un estimateur ponctuel d&rsquo;un param\u00e8tre populationnel inconnu. En tant que telle, elle est soumise \u00e0 la variabilit\u00e9 d&rsquo;\u00e9chantillonnage stochastique. Communiquer une MAE sans son intervalle de confiance associ\u00e9 prive le lecteur de toute estimation de la pr\u00e9cision de la m\u00e9trique elle-m\u00eame. La distribution th\u00e9orique de la moyenne des \u00e9carts absolus n&rsquo;ob\u00e9issant pas n\u00e9cessairement \u00e0 une loi normale connue \u2014 en particulier lorsque les r\u00e9sidus sont asym\u00e9triques \u2014, le recours au r\u00e9\u00e9chantillonnage non param\u00e9trique par <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Bootstrapping_(statistics)\">bootstrap<\/a> s&rsquo;impose comme la m\u00e9thodologie de r\u00e9f\u00e9rence.<\/p>\n<p>L&rsquo;impl\u00e9mentation formelle sous R s&rsquo;articule \u00e9l\u00e9gamment autour du package standard <code>boot<\/code>. La premi\u00e8re \u00e9tape consiste \u00e0 programmer la fonction statistique qui sera \u00e9valu\u00e9e sur chaque it\u00e9ration r\u00e9\u00e9chantillonn\u00e9e avec remise :<\/p>\n<p><code>library(boot)<\/p>\n<p>statistique_mae_boot &lt;- function(data, indices) {<br \/>\n&nbsp;&nbsp;d &lt;- data[indices, ]<br \/>\n&nbsp;&nbsp;return(Metrics::mae(d<span class=\"arabpsy-math-inline\" data-latex=\"observe, d\">$observe, d$<\/span>predit))<br \/>\n}<\/code><\/p>\n<p>L&rsquo;orchestration des it\u00e9rations de r\u00e9\u00e9chantillonnage s&rsquo;ex\u00e9cute ensuite via la fonction principale <code>boot()<\/code>, configur\u00e9e usuellement pour g\u00e9n\u00e9rer au minimum 2000 r\u00e9plications empiriques :<\/p>\n<p><code>set.seed(1234)<br \/>\nresultats_bootstrap &lt;- boot(<br \/>\n&nbsp;&nbsp;data = data.frame(observe = valeurs_reelles, predit = valeurs_estimees),<br \/>\n&nbsp;&nbsp;statistic = statistique_mae_boot,<br \/>\n&nbsp;&nbsp;R = 2000<br \/>\n)<\/code><\/p>\n<p>L&rsquo;extraction de l&rsquo;intervalle de confiance non param\u00e9trique s&rsquo;op\u00e8re enfin via la m\u00e9thode <code>boot.ci()<\/code>. Il est fortement pr\u00e9conis\u00e9 de retenir l&rsquo;intervalle corrig\u00e9 du biais et acc\u00e9l\u00e9r\u00e9 (BCa pour <em>bias-corrected and accelerated<\/em>), qui ajuste de fa\u00e7on rigoureuse l&rsquo;estimation en pr\u00e9sence d&rsquo;asym\u00e9trie dans la distribution bootstrap des erreurs absolues :<\/p>\n<p><code>intervalle_bca &lt;- boot.ci(resultats_bootstrap, type = \"bca\")<br \/>\nprint(intervalle_bca)<\/code><\/p>\n<p>La formulation acad\u00e9mique finale s&rsquo;\u00e9noncera alors selon la typographie APA recommand\u00e9e : <em>\u00ab\u00a0Le mod\u00e8le pr\u00e9dictif affiche une erreur absolue moyenne globale MAE = 3,42 points, IC BCa \u00e0 95 % [3,12 ; 3,78], attestant d&rsquo;une fid\u00e9lit\u00e9 pr\u00e9dictive statistiquement stable au sein de la population cible.\u00a0\u00bb<\/em><\/p>\n<h3>12.3 Synth\u00e8se et arbre de d\u00e9cision m\u00e9thodologique<\/h3>\n<p>Pour clore cet examen exhaustif de l&rsquo;erreur absolue moyenne dans R, il convient de synth\u00e9tiser la d\u00e9marche sous la forme d&rsquo;un arbre d\u00e9cisionnel guidant le choix m\u00e9thodologique de l&rsquo;analyste face \u00e0 des donn\u00e9es continues :<\/p>\n<ul>\n<li><strong>Question 1 : Quelle est la nature de l&rsquo;\u00e9chelle de la variable cible ?<\/strong>\n<ul>\n<li>Si l&rsquo;\u00e9chelle poss\u00e8de des bornes arbitraires sans z\u00e9ro absolu (ex. : \u00e9chelles d&rsquo;attitudes, scores Z, inventaires psychom\u00e9triques), <em>exclure imp\u00e9rativement la MAPE<\/em> et privil\u00e9gier la <strong>MAE<\/strong> ou la <strong>RMSE<\/strong>.<\/li>\n<li>Si l&rsquo;\u00e9chelle poss\u00e8de un z\u00e9ro physique absolu et strictement positif avec des besoins de comparaison d&rsquo;\u00e9chelles disparates (ex. : volumes mon\u00e9taires relatifs), la MAPE peut \u00eatre consult\u00e9e avec r\u00e9serve.<\/li>\n<\/ul>\n<\/li>\n<li><strong>Question 2 : Quelle est la politique de p\u00e9nalisation des d\u00e9viations extr\u00eames ?<\/strong>\n<ul>\n<li>Si quelques erreurs massives ont des cons\u00e9quences sanitaires ou s\u00e9curitaires critiques et doivent \u00eatre p\u00e9nalis\u00e9es avec une s\u00e9v\u00e9rit\u00e9 disproportionn\u00e9e : retenir la <strong>RMSE<\/strong> comme crit\u00e8re d&rsquo;optimisation central.<\/li>\n<li>Si l&rsquo;objectif est d&rsquo;estimer l&rsquo;erreur typique, robuste, attendue pour le sujet m\u00e9dian en limitant la vuln\u00e9rabilit\u00e9 aux donn\u00e9es contamin\u00e9es : opter r\u00e9solument pour la <strong>MAE<\/strong>.<\/li>\n<\/ul>\n<\/li>\n<li><strong>Question 3 : Quelle architecture logicielle d\u00e9ployer dans R ?<\/strong>\n<ul>\n<li>Pour une autonomie totale et des calculs vectoriels ultra-rapides sans d\u00e9pendance : <code>mean(abs(reel - pred), na.rm = TRUE)<\/code>.<\/li>\n<li>Pour des scripts d&rsquo;\u00e9valuation classiques et normalis\u00e9s : <code>Metrics::mae(actual, predicted)<\/code>.<\/li>\n<li>Pour des pipelines d&rsquo;apprentissage supervis\u00e9 modernes int\u00e9grant la validation crois\u00e9e et le calcul distribu\u00e9 par strates : <code>yardstick::mae()<\/code> au sein de <code>tidymodels<\/code>.<\/li>\n<\/ul>\n<\/li>\n<\/ul>\n<p>En int\u00e9grant ces principes de rigueur th\u00e9orique, d&rsquo;int\u00e9grit\u00e9 num\u00e9rique et de pertinence op\u00e9rationnelle, l&rsquo;analyste quantitative conf\u00e8re \u00e0 ses mod\u00e9lisations sous R le niveau d&rsquo;exactitude et de robustesse exig\u00e9 par la science des donn\u00e9es contemporaine.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7th ed.). American Psychological Association. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/li>\n<li>Box, G. E., &amp; Cox, D. R. (1964). An analysis of transformations. <em>Journal of the Royal Statistical Society: Series B (Methodological)<\/em>, 26(2), 211\u2013243. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x\">https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x<\/a><\/li>\n<li>Canty, A., &amp; Ripley, B. D. (2021). <em>boot: Bootstrap R (S-Plus) Functions<\/em> (R package version 1.3-28) [Computer software]. Comprehensive R Archive Network. <a rel=\"noopener noreferrer\" href=\"https:\/\/CRAN.R-project.org\/package=boot\">https:\/\/CRAN.R-project.org\/package=boot<\/a><\/li>\n<li>Davison, A. C., &amp; Hinkley, D. V. (1997). <em>Bootstrap methods and their application<\/em>. Cambridge University Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1017\/CBO9780511802843\">https:\/\/doi.org\/10.1017\/CBO9780511802843<\/a><\/li>\n<li>Frane, A. V. (2015). Power and type I error rate of new variants of the Mann-Whitney-Wilcoxon test. <em>Behavior Research Methods<\/em>, 47(4), 1146\u20131156. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.3758\/s13428-014-0545-2\">https:\/\/doi.org\/10.3758\/s13428-014-0545-2<\/a><\/li>\n<li>Hamner, B., &amp; Frasco, M. (2018). <em>Metrics: Evaluation metrics for machine learning<\/em> (R package version 0.1.4) [Computer software]. Comprehensive R Archive Network. <a rel=\"noopener noreferrer\" href=\"https:\/\/CRAN.R-project.org\/package=Metrics\">https:\/\/CRAN.R-project.org\/package=Metrics<\/a><\/li>\n<li>Huber, P. J. (1981). <em>Robust statistics<\/em>. John Wiley &amp; Sons. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471725250\">https:\/\/doi.org\/10.1002\/0471725250<\/a><\/li>\n<li>Hyndman, R. J., &amp; Koehler, A. B. (2006). Another look at measures of forecast accuracy. <em>International Journal of Forecasting<\/em>, 22(4), 679\u2013688. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/j.ijforecast.2006.03.001\">https:\/\/doi.org\/10.1016\/j.ijforecast.2006.03.001<\/a><\/li>\n<li>Kuhn, M. (2008). Building predictive models in R using the caret package. <em>Journal of Statistical Software<\/em>, 28(5), 1\u201326. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.18637\/jss.v028.i05\">https:\/\/doi.org\/10.18637\/jss.v028.i05<\/a><\/li>\n<li>Kuhn, M., &amp; Wickham, H. (2020). <em>Tidymodels: A collection of packages for modeling and machine learning using tidyverse principles<\/em> [Computer software]. Comprehensive R Archive Network. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.tidymodels.org\">https:\/\/www.tidymodels.org<\/a><\/li>\n<li>R Core Team. (2023). <em>R: A language and environment for statistical computing<\/em> [Computer software]. R Foundation for Statistical Computing. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.R-project.org\/\">https:\/\/www.R-project.org\/<\/a><\/li>\n<li>Wickham, H. (2016). <em>ggplot2: Elegant graphics for data analysis<\/em>. Springer-Verlag. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-3-319-24277-4\">https:\/\/doi.org\/10.1007\/978-3-319-24277-4<\/a><\/li>\n<li>Willmott, C. J., &amp; Matsuura, K. (2005). Advantages of the mean absolute error (MAE) over the root mean squared error (RMSE) in assessing average model performance. <em>Climate Research<\/em>, 30(1), 79\u201382. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.3354\/cr030079\">https:\/\/doi.org\/10.3354\/cr030079<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet d\u00e9taillant le calcul de l&rsquo;erreur absolue moyenne (MAE) dans R : th\u00e9orie statistique, syntaxe, packages et applications empiriques.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2193],"tags":[9394,2196,9395,16,2658],"class_list":["post-8251","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-programmation-r","tag-erreur-absolue-moyenne","tag-modelisation-statistique","tag-package-metrics","tag-psychometrie","tag-r"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8251","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=8251"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/8251\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=8251"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=8251"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=8251"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}