{"id":6309,"date":"2026-09-16T07:54:52","date_gmt":"2026-09-16T07:54:52","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-les-residus-analyse-regression\/"},"modified":"2026-09-16T07:54:58","modified_gmt":"2026-09-16T07:54:58","slug":"comment-calculer-les-residus-analyse-regression","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-les-residus-analyse-regression","title":{"rendered":"Comment calculer les r\u00e9sidus dans l&rsquo;analyse de r\u00e9gression"},"content":{"rendered":"<p>Dans le champ des sciences humaines et comportementales, la mod\u00e9lisation statistique ne constitue pas une simple formalit\u00e9 calculatoire, mais un instrument \u00e9pist\u00e9mologique d&rsquo;investigation des r\u00e9gularit\u00e9s psychologiques. Qu&rsquo;il s&rsquo;agisse d&rsquo;\u00e9valuer l&rsquo;impact du stress sur la m\u00e9moire de travail, de pr\u00e9dire l&rsquo;adaptation sociale \u00e0 partir de traits de personnalit\u00e9, ou de sonder l&rsquo;efficacit\u00e9 d&rsquo;un protocole psychoth\u00e9rapeutique, les chercheurs recourent fr\u00e9quemment aux mod\u00e8les pr\u00e9dictifs lin\u00e9aires. Toutefois, la valeur heuristique d&rsquo;un mod\u00e8le ne se mesure pas uniquement \u00e0 l&rsquo;aune de sa capacit\u00e9 \u00e0 tracer une tendance centrale \u00e9l\u00e9gante \u00e0 travers un nuage d&rsquo;observations empiriques. Elle repose de mani\u00e8re tout aussi critique sur l&rsquo;exploration m\u00e9thodique de ce que l&rsquo;\u00e9quation laisse de c\u00f4t\u00e9, c&rsquo;est-\u00e0-dire les d\u00e9viations individuelles que l&rsquo;on nomme formellement les r\u00e9sidus.<\/p>\n<p>Le calcul et le diagnostic des r\u00e9sidus incarnent la pierre angulaire de l&rsquo;inf\u00e9rence statistique moderne. Trop souvent rel\u00e9gu\u00e9e au rang d&rsquo;artefact secondaire par les utilisateurs novices de logiciels de traitement de donn\u00e9es, l&rsquo;analyse r\u00e9siduelle offre un miroir direct de la validit\u00e9 interne d&rsquo;une recherche. En confrontant les pr\u00e9dictions th\u00e9oriques issues de la droite de r\u00e9gression aux donn\u00e9es effectivement recueillies sur le terrain, le r\u00e9sidu mat\u00e9rialise la part d&rsquo;inconnu, de bruit de mesure ou de complexit\u00e9 psychologique non capt\u00e9e par les variables explicatives s\u00e9lectionn\u00e9es. Loin d&rsquo;\u00eatre un simple d\u00e9chet algorithmique, il constitue la mati\u00e8re premi\u00e8re \u00e0 partir de laquelle s&rsquo;\u00e9value l&rsquo;ad\u00e9quation d&rsquo;un postulat th\u00e9orique face \u00e0 la r\u00e9alit\u00e9 empirique.<\/p>\n<p>Ce guide exhaustif a pour vocation d&rsquo;accompagner les chercheurs, psychom\u00e9triciens, \u00e9tudiants de cycles sup\u00e9rieurs et analystes de donn\u00e9es dans la ma\u00eetrise int\u00e9grale des r\u00e9sidus en analyse de r\u00e9gression. De la d\u00e9rivation alg\u00e9brique fondamentale aux proc\u00e9dures computationnelles avanc\u00e9es sous R, Python, SPSS et JASP, nous examinerons en d\u00e9tail les fondements math\u00e9matiques, les typologies de standardisation, les tests de conformit\u00e9 des postulats de Gauss-Markov et les strat\u00e9gies de rem\u00e9diation face aux donn\u00e9es atypiques. \u00c0 travers une d\u00e9marche r\u00e9solument ancr\u00e9e dans la rigueur m\u00e9thodologique, cet ouvrage propose une immersion compl\u00e8te au c\u0153ur de l&rsquo;analyse des \u00e9carts, transformant l&rsquo;erreur de pr\u00e9diction en un puissant vecteur de d\u00e9couverte scientifique.<\/p>\n<h2>1. Introduction aux r\u00e9sidus dans l&rsquo;analyse de r\u00e9gression lin\u00e9aire en psychologie<\/h2>\n<h3>1.1 D\u00e9finition conceptuelle et r\u00f4le \u00e9pist\u00e9mologique du r\u00e9sidu<\/h3>\n<p>D&rsquo;un point de vue \u00e9pist\u00e9mologique, la psychologie scientifique tente de formaliser des lois g\u00e9n\u00e9rales du comportement tout en composant avec l&rsquo;in\u00e9vitable singularit\u00e9 des individus observ\u00e9s. Lorsqu&rsquo;un chercheur applique un mod\u00e8le de r\u00e9gression lin\u00e9aire, il postule qu&rsquo;une variable ind\u00e9pendante exerce une influence structurelle, syst\u00e9matique et mesurable sur une variable d\u00e9pendante. N\u00e9anmoins, l&rsquo;observation d&rsquo;un participant particulier ne co\u00efncide que tr\u00e8s rarement avec la trajectoire moyenne pr\u00e9dite par la th\u00e9orie. Cet \u00e9cart quantitatif pr\u00e9cis, qui s\u00e9pare la r\u00e9ponse empirique manifest\u00e9e par le sujet de la r\u00e9ponse attendue selon la fonction math\u00e9matique, d\u00e9finit formellement le r\u00e9sidu.<\/p>\n<p>Dans le cadre de la mod\u00e9lisation comportementale, l&rsquo;erreur d&rsquo;estimation ne doit pas \u00eatre appr\u00e9hend\u00e9e comme une banale maladresse exp\u00e9rimentale ou une imperfection n\u00e9gligeable de l&rsquo;instrumentation. Elle refl\u00e8te la part intrins\u00e8que de variance psychologique non expliqu\u00e9e par les d\u00e9terminants int\u00e9gr\u00e9s au mod\u00e8le. Cette variance r\u00e9siduelle agr\u00e8ge de multiples sources : fluctuations attentionnelles passag\u00e8res, ant\u00e9c\u00e9dents biographiques idiosyncrasiques, variations neurobiologiques momentan\u00e9es, erreurs al\u00e9atoires d&rsquo;\u00e9chantillonnage psychom\u00e9trique ou omissions d\u00e9lib\u00e9r\u00e9es d&rsquo;autres variables explicatives pertinentes au sein de l&rsquo;architecture du protocole exp\u00e9rimental.<\/p>\n<p>Il est imp\u00e9ratif d&rsquo;op\u00e9rer une distinction s\u00e9mantique et math\u00e9matique irr\u00e9ductible entre l&rsquo;erreur th\u00e9orique de population, traditionnellement not\u00e9e par la lettre grecque epsilon (&epsilon;), et le r\u00e9sidu empirique \u00e9chantillonnal, d\u00e9sign\u00e9 par la lettre latine <em>e<\/em>. L&rsquo;erreur de population repr\u00e9sente un param\u00e8tre latent et fondamentalement inobservable ; elle renvoie \u00e0 la distance exacte entre une observation individuelle et la droite de r\u00e9gression vraie, issue de la population globale th\u00e9orique. \u00c0 l&rsquo;inverse, le r\u00e9sidu empirique est une quantit\u00e9 calculable, estim\u00e9e directement \u00e0 partir d&rsquo;un \u00e9chantillon fini de donn\u00e9es via l&rsquo;ajustement des param\u00e8tres d&rsquo;\u00e9chantillonnage. Ignorer cette divergence conceptuelle reviendrait \u00e0 confondre les propri\u00e9t\u00e9s ontologiques de l&rsquo;univers avec les approximations statistiques permises par un \u00e9chantillon circonscrit.<\/p>\n<p>L&rsquo;analyse rigoureuse des r\u00e9sidus constitue ainsi le m\u00e9canisme privil\u00e9gi\u00e9 par lequel la psychologie quantitative \u00e9prouve la robustesse de ses \u00e9difices th\u00e9oriques. Si un mod\u00e8le conceptuel pr\u00e9tend rendre compte d&rsquo;un ph\u00e9nom\u00e8ne psychologique pr\u00e9cis, les \u00e9carts \u00e0 ce mod\u00e8le ne doivent manifester aucune structure logique, aucune directionnalit\u00e9 sous-jacente ni aucune r\u00e9gularit\u00e9 non ma\u00eetris\u00e9e. L&rsquo;examen des r\u00e9sidus offre alors un test de falsification popp\u00e9rien par excellence : la pr\u00e9sence de r\u00e9gularit\u00e9s dans l&rsquo;erreur d\u00e9montre sans \u00e9quivoque que la th\u00e9orie sous-jacente est incompl\u00e8te ou fondamentalement mal sp\u00e9cifi\u00e9e.<\/p>\n<h3>1.2 Le mod\u00e8le de r\u00e9gression lin\u00e9aire simple : variables pr\u00e9dictive et r\u00e9ponse<\/h3>\n<p>La formulation alg\u00e9brique canonique de la <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/R%C3%A9gression_lin%C3%A9aire_simple\">r\u00e9gression lin\u00e9aire simple<\/a> s&rsquo;\u00e9nonce selon l&rsquo;\u00e9quation fondamentale : <em>Y = b<sub>0<\/sub> + b<sub>1<\/sub>X + e<\/em>. Dans ce dispositif, <em>Y<\/em> repr\u00e9sente la variable r\u00e9ponse, \u00e9galement qualifi\u00e9e de crit\u00e8re ou d&rsquo;observation d\u00e9pendante, tandis que <em>X<\/em> d\u00e9signe le pr\u00e9dicteur, commun\u00e9ment identifi\u00e9 comme la variable explicative ou ind\u00e9pendante. Les coefficients <em>b<sub>0<\/sub><\/em> et <em>b<sub>1<\/sub><\/em> constituent respectivement l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine (ou constante du mod\u00e8le) et le coefficient de pente empirique, traduisant l&rsquo;amplitude de la variation attendue de <em>Y<\/em> cons\u00e9cutivement \u00e0 une augmentation unitaire de <em>X<\/em>.<\/p>\n<p>Consid\u00e9rons une investigation psychom\u00e9trique classique visant \u00e0 mod\u00e9liser le niveau d&rsquo;anxi\u00e9t\u00e9 situationnelle (<em>Y<\/em>), op\u00e9rationnalis\u00e9 par un score standardis\u00e9 sur une \u00e9chelle d&rsquo;\u00e9valuation clinique, en fonction de la charge cognitive de travail (<em>X<\/em>) impos\u00e9e lors d&rsquo;une t\u00e2che d&rsquo;attention divis\u00e9e informatis\u00e9e. L&rsquo;\u00e9quation de r\u00e9gression cherche \u00e0 formaliser comment l&rsquo;accumulation des contraintes cognitives amplifie le sentiment d&rsquo;anxi\u00e9t\u00e9. D\u00e8s lors que le protocole recueille des donn\u00e9es r\u00e9elles, chaque participant pr\u00e9sente un score singulier d&rsquo;anxi\u00e9t\u00e9 qui s&rsquo;\u00e9carte de la valeur g\u00e9n\u00e9rale mod\u00e9lis\u00e9e pour son niveau de charge cognitive sp\u00e9cifique. L&rsquo;\u00e9cart entre ce score clinique recueilli et le score th\u00e9orique g\u00e9n\u00e9r\u00e9 par la combinaison lin\u00e9aire de la constante et de la pente constitue le r\u00e9sidu sp\u00e9cifique attribu\u00e9 \u00e0 cet individu.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/weight_height1.jpg\" alt=\"Simple linear regression\" \/><figcaption>Simple linear regression<\/figcaption><\/figure>\n<p>Sur le plan de l&rsquo;interpr\u00e9tation g\u00e9om\u00e9trique, la r\u00e9gression lin\u00e9aire simple par la m\u00e9thode des moindres carr\u00e9s ordinaires (MCO) vise \u00e0 tracer une trajectoire rectiligne au sein du plan cart\u00e9sien bidimensionnel d\u00e9fini par les axes des abscisses (<em>X<\/em>) et des ordonn\u00e9es (<em>Y<\/em>). Cette droite optimale ne passe pas n\u00e9cessairement par l&rsquo;ensemble des points empiriques ; au contraire, elle traverse le nuage de points de telle fa\u00e7on que la somme des distances verticales au carr\u00e9 s\u00e9parant chaque point de la droite soit rendue rigoureusement minimale. Le r\u00e9sidu s&rsquo;incarne donc physiquement dans cette distance verticale mesur\u00e9e parall\u00e8lement \u00e0 l&rsquo;axe des ordonn\u00e9es, et non pas selon une projection perpendiculaire ou orthogonale directe \u00e0 la droite elle-m\u00eame.<\/p>\n<h3>1.3 Pourquoi le calcul exact des r\u00e9sidus est indispensable au chercheur<\/h3>\n<p>L&rsquo;estimation num\u00e9rique pr\u00e9cise de chaque r\u00e9sidu individuel n&rsquo;est pas un luxe computationnel : elle conditionne l&rsquo;\u00e9valuation de l&rsquo;ad\u00e9quation globale du mod\u00e8le statistique, concept d\u00e9sign\u00e9 dans la litt\u00e9rature internationale sous le terme d&rsquo;<a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Goodness_of_fit\">ad\u00e9quation du mod\u00e8le (goodness-of-fit)<\/a>. Il ne suffit pas d&rsquo;obtenir un coefficient de corr\u00e9lation \u00e9lev\u00e9 ou une valeur de significativit\u00e9 statistique infinit\u00e9simale (<em>p<\/em> &lt; .001) pour conclure \u00e0 la validit\u00e9 d&rsquo;une r\u00e9gression. Un mod\u00e8le peut afficher un coefficient de d\u00e9termination s\u00e9duisant tout en \u00e9chouant dramatiquement \u00e0 capturer la g\u00e9om\u00e9trie r\u00e9elle du ph\u00e9nom\u00e8ne \u00e9tudi\u00e9 en raison d&rsquo;une mauvaise r\u00e9partition interne des erreurs d&rsquo;estimation.<\/p>\n<p>Le calcul exhaustif des r\u00e9sidus repr\u00e9sente la d\u00e9marche incontournable pour v\u00e9rifier si les postulats fondamentaux du th\u00e9or\u00e8me de Gauss-Markov sont respect\u00e9s. Ces postulats &ndash; incluant la lin\u00e9arit\u00e9 intrins\u00e8que de la fonction, l&rsquo;homog\u00e9n\u00e9it\u00e9 de la variance des erreurs et l&rsquo;ind\u00e9pendance mutuelle des observations &ndash; reposent enti\u00e8rement sur les propri\u00e9t\u00e9s structurelles de la distribution r\u00e9siduelle. Si un chercheur omet d&rsquo;extraire et d&rsquo;analyser ses r\u00e9sidus, il court le risque grave d&rsquo;\u00e9mettre des inf\u00e9rences statistiques erron\u00e9es, de sous-estimer de fa\u00e7on trompeuse les erreurs types associ\u00e9es \u00e0 ses coefficients de r\u00e9gression et, in fine, de publier des conclusions scientifiques fallacieuses bas\u00e9es sur des faux positifs m\u00e9thodologiques.<\/p>\n<p>Au-del\u00e0 de la validation technique des hypoth\u00e8ses distributionnelles, l&rsquo;exploration des grandeurs r\u00e9siduelles s&rsquo;av\u00e8re f\u00e9conde pour le diagnostic clinique et psychom\u00e9trique appliqu\u00e9. L&rsquo;analyse des \u00e9carts met en lumi\u00e8re les profils comportementaux atypiques : des sujets dont les scores s&rsquo;av\u00e8rent radicalement divergents de la pr\u00e9diction commune. Ces individus aux r\u00e9sidus disproportionn\u00e9s peuvent r\u00e9v\u00e9ler des strat\u00e9gies cognitives compensatoires non r\u00e9pertori\u00e9es, des biais de d\u00e9sirabilit\u00e9 sociale dans la passation des auto-questionnaires, ou signaler l&rsquo;existence de sous-groupes cliniques m\u00e9connus dont la prise en compte n\u00e9cessitera une r\u00e9vision th\u00e9orique de grande envergure.<\/p>\n<h2>2. Fondements math\u00e9matiques et formulation alg\u00e9brique des r\u00e9sidus<\/h2>\n<h3>2.1 L&rsquo;\u00e9quation formelle du r\u00e9sidu ordinaire<\/h3>\n<p>Dans le formalisme de l&rsquo;analyse univari\u00e9e et multivari\u00e9e, l&rsquo;expression alg\u00e9brique canonique d\u00e9finissant le r\u00e9sidu ordinaire pour la <em>i<\/em>-i\u00e8me observation d&rsquo;un \u00e9chantillon de taille <em>n<\/em> s&rsquo;\u00e9crit avec une sobri\u00e9t\u00e9 math\u00e9matique universelle :<\/p>\n<p><em>e<sub>i<\/sub> = y<sub>i<\/sub> &minus; \u0177<sub>i<\/sub><\/em><\/p>\n<p>Dans cette formulation \u00e9l\u00e9mentaire, <em>y<sub>i<\/sub><\/em> symbolise le score empirique r\u00e9el recueilli aupr\u00e8s du participant exp\u00e9rimental <em>i<\/em>. Il s&rsquo;agit de la valeur brute, concr\u00e8te et non modifi\u00e9e, observ\u00e9e sur l&rsquo;instrument psychom\u00e9trique ou le capteur psychophysiologique. En contrepartie, <em>\u0177<sub>i<\/sub><\/em> (prononc\u00e9 <em>y chapeau<\/em>) caract\u00e9rise la valeur ajust\u00e9e ou pr\u00e9dite par l&rsquo;\u00e9quation de r\u00e9gression lin\u00e9aire pour ce participant particulier, calcul\u00e9e en injectant sa valeur observ\u00e9e sur le pr\u00e9dicteur <em>x<sub>i<\/sub><\/em> dans le syst\u00e8me : <em>\u0177<sub>i<\/sub> = b<sub>0<\/sub> + b<sub>1<\/sub>x<sub>i<\/sub><\/em>.<\/p>\n<p>L&rsquo;analyse du signe arithm\u00e9tique issu de cette soustraction rev\u00eat une importance analytique directe. Lorsque le r\u00e9sidu est strictement positif (<em>e<sub>i<\/sub> &gt; 0<\/em>), la valeur observ\u00e9e <em>y<sub>i<\/sub><\/em> est sup\u00e9rieure \u00e0 la pr\u00e9diction th\u00e9orique <em>\u0177<sub>i<\/sub><\/em> : le mod\u00e8le a donc commis une sous-estimation du score du sujet. R\u00e9ciproquement, d\u00e8s lors que le r\u00e9sidu affiche une valeur strictement n\u00e9gative (<em>e<sub>i<\/sub> &lt; 0<\/em>), la valeur empirique se situe en de\u00e7\u00e0 de la coordonn\u00e9e projet\u00e9e sur la droite : le mod\u00e8le math\u00e9matique a produit une surestimation manifeste de la r\u00e9ponse psychologique r\u00e9elle.<\/p>\n<p>Un r\u00e9sidu d&rsquo;une magnitude rigoureusement nulle (<em>e<sub>i<\/sub> = 0<\/em>) caract\u00e9rise une co\u00efncidence parfaite, indiquant que le point empirique repr\u00e9sentatif du sujet repose exactement sur l&rsquo;axe rectiligne de l&rsquo;\u00e9quation pr\u00e9dictive. Dans la pratique empirique des sciences du comportement, cette situation demeure rarissime en raison de l&rsquo;in\u00e9vitable variance d&rsquo;erreur inh\u00e9rente \u00e0 la mesure du comportement humain.<\/p>\n<h3>2.2 La propri\u00e9t\u00e9 fondamentale de la somme nulle des moindres carr\u00e9s<\/h3>\n<p>L&rsquo;une des propri\u00e9t\u00e9s fondamentales de l&rsquo;estimation par les moindres carr\u00e9s ordinaires r\u00e9side dans la nullit\u00e9 obligatoire de la somme alg\u00e9brique des r\u00e9sidus ordinaires, d\u00e8s lors que le mod\u00e8le int\u00e8gre une ordonn\u00e9e \u00e0 l&rsquo;origine (<em>b<sub>0<\/sub><\/em>). Cette caract\u00e9ristique d\u00e9coule de la d\u00e9rivation analytique des conditions du premier ordre de minimisation math\u00e9matique.<\/p>\n<p>Pour optimiser le mod\u00e8le, l&rsquo;algorithme des MCO cherche \u00e0 minimiser la fonctionnelle de perte repr\u00e9sent\u00e9e par la somme des carr\u00e9s des erreurs r\u00e9siduelles (not\u00e9e <em>SSE<\/em> pour <em>Sum of Squared Errors<\/em>) :<\/p>\n<p><em>S(b<sub>0<\/sub>, b<sub>1<\/sub>) = &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &minus; b<sub>0<\/sub> &minus; b<sub>1<\/sub>x<sub>i<\/sub>)<sup>2<\/sup><\/em><\/p>\n<p>En calculant la d\u00e9riv\u00e9e partielle de cette fonction d&rsquo;erreur par rapport au param\u00e8tre constant <em>b<sub>0<\/sub><\/em> et en posant cette d\u00e9riv\u00e9e \u00e9gale \u00e0 z\u00e9ro pour atteindre l&rsquo;extremum local, nous obtenons l&rsquo;\u00e9galit\u00e9 suivante :<\/p>\n<p><em>&part;S \/ &part;b<sub>0<\/sub> = &minus;2 &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &minus; b<sub>0<\/sub> &minus; b<sub>1<\/sub>x<sub>i<\/sub>) = 0<\/em><\/p>\n<p>En \u00e9liminant le facteur multiplicatif &minus;2 et en reconnaissant que l&rsquo;expression entre parenth\u00e8ses \u00e9quivaut rigoureusement au r\u00e9sidu individuel <em>e<sub>i<\/sub><\/em>, la relation d\u00e9bouche in\u00e9vitablement sur la loi fondamentale :<\/p>\n<p><em>&sum;<sub>i=1<\/sub><sup>n<\/sup> e<sub>i<\/sub> = 0<\/em><\/p>\n<p>Cette contrainte structurelle engendre une cons\u00e9quence directe : la moyenne arithm\u00e9tique des r\u00e9sidus bruts au sein de l&rsquo;\u00e9chantillon d&rsquo;ajustement est invariablement et exactement \u00e9gale \u00e0 z\u00e9ro (<em>\u0113 = 0<\/em>). Sur le plan de la g\u00e9om\u00e9trie vectorielle hilbertienne, cela implique que le vecteur des r\u00e9sidus se d\u00e9ploie orthogonalement \u00e0 l&rsquo;hyperplan sous-tendu par les pr\u00e9dicteurs du mod\u00e8le et le vecteur unitaire associ\u00e9 \u00e0 la constante. Le produit scalaire entre le vecteur des valeurs pr\u00e9dites <em>\u0176<\/em> et le vecteur des r\u00e9sidus <em>e<\/em> est donc nul (<em>\u0176<sup>T<\/sup>e = 0<\/em>), ce qui formalise l&rsquo;ind\u00e9pendance lin\u00e9aire stricte entre la pr\u00e9diction et l&rsquo;erreur d&rsquo;ajustement globale.<\/p>\n<h3>2.3 D\u00e9composition de la variance totale et r\u00f4le des r\u00e9sidus<\/h3>\n<p>L&rsquo;orthogonalit\u00e9 math\u00e9matique entre les pr\u00e9dictions du mod\u00e8le et les erreurs calcul\u00e9es constitue le pilier du principe de d\u00e9composition de la variance, central dans l&rsquo;analyse de r\u00e9gression et l&rsquo;analyse de variance (<a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Analyse_de_la_variance\">ANOVA<\/a>). La dispersion totale des scores observ\u00e9s par rapport \u00e0 leur moyenne g\u00e9n\u00e9rale s&rsquo;articule autour de la Somme des Carr\u00e9s Totale (<em>SST<\/em>, <em>Total Sum of Squares<\/em>) :<\/p>\n<p><em>SST = &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &minus; \u0233)<sup>2<\/sup><\/em><\/p>\n<p>Cette variabilit\u00e9 globale se partitionne de fa\u00e7on exacte en deux composantes mutuellement exclusives et exhaustives : la Somme des Carr\u00e9s de la R\u00e9gression (<em>SSR<\/em>, <em>Regression Sum of Squares<\/em>), qui capture la dispersion des pr\u00e9dictions th\u00e9oriques autour de la moyenne, et la Somme des Carr\u00e9s des Erreurs (<em>SSE<\/em>, <em>Residual Sum of Squares<\/em>), qui agr\u00e8ge l&rsquo;amplitude quadratique des r\u00e9sidus bruts :<\/p>\n<p><em>SST = SSR + SSE<\/em>, soit : <em>&sum; (y<sub>i<\/sub> &minus; \u0233)<sup>2<\/sup> = &sum; (\u0177<sub>i<\/sub> &minus; \u0233)<sup>2<\/sup> + &sum; (y<sub>i<\/sub> &minus; \u0177<sub>i<\/sub>)<sup>2<\/sup><\/em><\/p>\n<p>Cette relation met en \u00e9vidence la nature du <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Coefficient_de_d%C3%A9termination\">coefficient de d\u00e9termination<\/a> <em>R<sup>2<\/sup><\/em>, qui mesure la proportion de la variance de la variable crit\u00e8re directement attribuable aux variations de la variable explicative. Alg\u00e9briquement, le <em>R<sup>2<\/sup><\/em> peut \u00eatre d\u00e9riv\u00e9 directement de la quantit\u00e9 r\u00e9siduelle calcul\u00e9e :<\/p>\n<p><em>R<sup>2<\/sup> = 1 &minus; (SSE \/ SST)<\/em><\/p>\n<p>Enfin, la variance r\u00e9siduelle de l&rsquo;\u00e9chantillon, appel\u00e9e \u00e9galement le Carr\u00e9 Moyen de l&rsquo;Erreur (<em>MSE<\/em> pour <em>Mean Squared Error<\/em>), s&rsquo;obtient en divisant la somme des carr\u00e9s des r\u00e9sidus par les degr\u00e9s de libert\u00e9 r\u00e9siduels associ\u00e9s au syst\u00e8me, \u00e0 savoir la taille de l&rsquo;\u00e9chantillon <em>n<\/em> amput\u00e9e du nombre de pr\u00e9dicteurs <em>k<\/em> et de la constante unitaire :<\/p>\n<p><em>s<sub>e<\/sub><sup>2<\/sup> = MSE = SSE \/ (n &minus; k &minus; 1)<\/em><\/p>\n<p>L&rsquo;\u00e9cart-type r\u00e9siduel (<em>s<sub>e<\/sub> = &radic;MSE<\/em>), souvent d\u00e9nomm\u00e9 erreur type de l&rsquo;estimation, quantifie la dispersion moyenne des observations empiriques autour de la droite de r\u00e9gression finale.<\/p>\n<h2>3. Proc\u00e9dure manuelle \u00e9tape par \u00e9tape pour calculer les r\u00e9sidus<\/h2>\n<h3>3.1 \u00c9tape 1 : D\u00e9termination des param\u00e8tres de la droite de r\u00e9gression<\/h3>\n<p>L&rsquo;extraction manuelle des r\u00e9sidus exige en premier lieu l&rsquo;estimation rigoureuse des coefficients de la fonction de r\u00e9gression par les MCO. Pour une r\u00e9gression lin\u00e9aire simple reliant un pr\u00e9dicteur <em>X<\/em> \u00e0 une variable d\u00e9pendante <em>Y<\/em>, le calcul commence par la d\u00e9termination de la moyenne empirique de chaque variable, not\u00e9es respectivement <em>x\u0304<\/em> et <em>\u0233<\/em>. L&rsquo;analyste doit ensuite calculer la somme des produits des \u00e9carts \u00e0 la moyenne &ndash; \u00e9quivalente \u00e0 la covariance non normalis\u00e9e &ndash; ainsi que la somme des carr\u00e9s des \u00e9carts de la variable explicative <em>X<\/em> (sa variance non normalis\u00e9e).<\/p>\n<p>Le coefficient de r\u00e9gression de la pente (<em>b<sub>1<\/sub><\/em>) est d\u00e9riv\u00e9 en \u00e9tablissant le quotient suivant :<\/p>\n<p><em>b<sub>1<\/sub> = [&sum;<sub>i=1<\/sub><sup>n<\/sup> (x<sub>i<\/sub> &minus; x\u0304)(y<sub>i<\/sub> &minus; \u0233)] \/ [&sum;<sub>i=1<\/sub><sup>n<\/sup> (x<sub>i<\/sub> &minus; x\u0304)<sup>2<\/sup>]<\/em><\/p>\n<p>Une fois le coefficient directeur obtenu, l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine (<em>b<sub>0<\/sub><\/em>) d\u00e9coule imm\u00e9diatement du fait obligatoire que la droite de r\u00e9gression des moindres carr\u00e9s passe invariablement par le centre de gravit\u00e9 bivari\u00e9 de l&rsquo;\u00e9chantillon, c&rsquo;est-\u00e0-dire le point de coordonn\u00e9es moyennes <em>(x\u0304, \u0233)<\/em>. L&rsquo;\u00e9quation se r\u00e9sout alors de fa\u00e7on lin\u00e9aire :<\/p>\n<p><em>b<sub>0<\/sub> = \u0233 &minus; b<sub>1<\/sub>x\u0304<\/em><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/statology_residuals1.jpg\" alt=\"Scatterplot example\" \/><figcaption>Scatterplot example<\/figcaption><\/figure>\n<p>L&rsquo;\u00e9quation pr\u00e9dictive finale s&rsquo;\u00e9crit d\u00e8s lors sous sa forme op\u00e9rationnelle : <em>\u0177 = b<sub>0<\/sub> + b<sub>1<\/sub>x<\/em>. Cette \u00e9quation fixe le standard pr\u00e9dictif \u00e0 partir duquel l&rsquo;ensemble des \u00e9carts individuels va \u00eatre \u00e9valu\u00e9.<\/p>\n<h3>3.2 \u00c9tape 2 : Calcul syst\u00e9matique des valeurs pr\u00e9dites (\u0177)<\/h3>\n<p>La deuxi\u00e8me s\u00e9quence op\u00e9ratoire consiste \u00e0 calculer, pour chacun des participants composant la cohorte d&rsquo;\u00e9tude, la valeur sp\u00e9cifiquement attendue selon le mod\u00e8le. Pour chaque sujet <em>i<\/em>, la valeur mesur\u00e9e sur le pr\u00e9dicteur <em>x<sub>i<\/sub><\/em> est introduite directement dans l&rsquo;\u00e9quation param\u00e9tr\u00e9e lors de l&rsquo;\u00e9tape pr\u00e9c\u00e9dente.<\/p>\n<p>Il est fondamental de maintenir une rigueur dimensionnelle et arithm\u00e9tique totale lors de cette substitution. Si l&rsquo;\u00e9valuation de la charge de travail est mesur\u00e9e sur une \u00e9chelle de 0 \u00e0 100, la coordonn\u00e9e <em>x<sub>i<\/sub><\/em> doit \u00eatre trait\u00e9e avec pr\u00e9cision sans arrondis h\u00e2tifs, sous peine d&rsquo;induire des erreurs d&rsquo;arrondi cumulatives lors de l&rsquo;estimation subs\u00e9quente des composantes r\u00e9siduelles.<\/p>\n<p>Dans le cadre d&rsquo;un traitement psychom\u00e9trique rigoureux, il est recommand\u00e9 de compiler ces valeurs pr\u00e9dites dans une colonne structur\u00e9e au sein d&rsquo;une grille ou d&rsquo;un tableau d&rsquo;analyse. Chaque ligne correspond \u00e0 un sujet unique, mat\u00e9rialisant le passage continu du niveau empirique brut au niveau th\u00e9orique pr\u00e9dit. Cette table servira d&rsquo;assise aux contr\u00f4les de coh\u00e9rence ult\u00e9rieurs.<\/p>\n<h3>3.3 \u00c9tape 3 : Soustraction pour obtenir les r\u00e9sidus bruts<\/h3>\n<p>L&rsquo;\u00e9tape charni\u00e8re r\u00e9side dans l&rsquo;application rigoureuse et s\u00e9quentielle de la formule de soustraction canonique <em>e<sub>i<\/sub> = y<sub>i<\/sub> &minus; \u0177<sub>i<\/sub><\/em> pour chacune des lignes du jeu de donn\u00e9es. La valeur th\u00e9orique calcul\u00e9e <em>\u0177<sub>i<\/sub><\/em> est retranch\u00e9e de la mesure observ\u00e9e <em>y<sub>i<\/sub><\/em>.<\/p>\n<p>Cette \u00e9tape requiert une attention soutenue \u00e0 la pr\u00e9servation rigoureuse des signes math\u00e9matiques. Une erreur fr\u00e9quente chez les chercheurs en d\u00e9but de parcours consiste \u00e0 inverser les termes de l&rsquo;op\u00e9rateur (en calculant <em>\u0177<sub>i<\/sub> &minus; y<sub>i<\/sub><\/em>) ou \u00e0 convertir d&rsquo;autorit\u00e9 l&rsquo;ensemble des r\u00e9sultats en valeurs absolues positives. Une telle inversion fausse compl\u00e8tement l&rsquo;analyse ult\u00e9rieure de la direction des biais de pr\u00e9diction.<\/p>\n<p>D\u00e8s l&rsquo;ach\u00e8vement de la colonne des r\u00e9sidus, une proc\u00e9dure de v\u00e9rification interne doit \u00eatre ex\u00e9cut\u00e9e : la sommation cumulative de l&rsquo;int\u00e9gralit\u00e9 des termes <em>e<sub>i<\/sub><\/em> doit converger vers z\u00e9ro (aux micro-impr\u00e9cisions d&rsquo;arrondis machine pr\u00e8s). Si la somme d\u00e9vie substantiellement de z\u00e9ro, cela signale une erreur syst\u00e9matique dans le calcul pr\u00e9alable des coefficients <em>b<sub>0<\/sub><\/em> ou <em>b<sub>1<\/sub><\/em>. Enfin, chaque r\u00e9sidu individuel est \u00e9lev\u00e9 au carr\u00e9 (<em>e<sub>i<\/sub><sup>2<\/sup><\/em>) dans une colonne adjacente ; leur agr\u00e9gation finale g\u00e9n\u00e8re la <em>SSE<\/em> indispensable au diagnostic d&rsquo;erreur global du protocole.<\/p>\n<h3>3.4 Exemple num\u00e9rique guid\u00e9 sur un jeu de donn\u00e9es psychologiques<\/h3>\n<p>Afin d&rsquo;illustrer concr\u00e8tement la s\u00e9quence des op\u00e9rations math\u00e9matiques d\u00e9crites, examinons un protocole exp\u00e9rimental men\u00e9 sur un \u00e9chantillon restreint de 10 participants adultes (<em>n = 10<\/em>). Cette \u00e9tude explore l&rsquo;influence du score de stress per\u00e7u (variable explicative <em>X<\/em>, mesur\u00e9e de 1 \u00e0 10 sur une \u00e9chelle de tension psychologique) sur le temps de r\u00e9action moyen \u00e0 une t\u00e2che d&rsquo;inhibition cognitive de type Stroop (variable r\u00e9ponse <em>Y<\/em>, mesur\u00e9e en centi\u00e8mes de seconde).<\/p>\n<p>Les donn\u00e9es brutes recueillies sont les suivantes :<\/p>\n<ul>\n<li>Sujet 1 : X = 2, Y = 35<\/li>\n<li>Sujet 2 : X = 3, Y = 42<\/li>\n<li>Sujet 3 : X = 4, Y = 40<\/li>\n<li>Sujet 4 : X = 5, Y = 55<\/li>\n<li>Sujet 5 : X = 6, Y = 50<\/li>\n<li>Sujet 6 : X = 7, Y = 68<\/li>\n<li>Sujet 7 : X = 8, Y = 62<\/li>\n<li>Sujet 8 : X = 9, Y = 75<\/li>\n<li>Sujet 9 : X = 10, Y = 82<\/li>\n<li>Sujet 10 : X = 11, Y = 81<\/li>\n<\/ul>\n<p>Calculons en premier lieu les statistiques descriptives univari\u00e9es de base :<\/p>\n<p><em>x\u0304 = (2 + 3 + 4 + 5 + 6 + 7 + 8 + 9 + 10 + 11) \/ 10 = 65 \/ 10 = 6,5<\/em><\/p>\n<p><em>\u0233 = (35 + 42 + 40 + 55 + 50 + 68 + 62 + 75 + 82 + 81) \/ 10 = 590 \/ 10 = 59,0<\/em><\/p>\n<p>Calculons ensuite la somme des carr\u00e9s des \u00e9carts de <em>X<\/em> et la somme des produits crois\u00e9s :<\/p>\n<p><em>&sum; (x<sub>i<\/sub> &minus; x\u0304)<sup>2<\/sup> = (-4,5)<sup>2<\/sup> + (-3,5)<sup>2<\/sup> + (-2,5)<sup>2<\/sup> + (-1,5)<sup>2<\/sup> + (-0,5)<sup>2<\/sup> + (0,5)<sup>2<\/sup> + (1,5)<sup>2<\/sup> + (2,5)<sup>2<\/sup> + (3,5)<sup>2<\/sup> + (4,5)<sup>2<\/sup> = 20,25 + 12,25 + 6,25 + 2,25 + 0,25 + 0,25 + 2,25 + 6,25 + 12,25 + 20,25 = 82,5<\/em><\/p>\n<p><em>&sum; (x<sub>i<\/sub> &minus; x\u0304)(y<sub>i<\/sub> &minus; \u0233) = [(-4,5)(-24)] + [(-3,5)(-17)] + [(-2,5)(-19)] + [(-1,5)(-4)] + [(-0,5)(-9)] + [(0,5)(9)] + [(1,5)(3)] + [(2,5)(16)] + [(3,5)(23)] + [(4,5)(22)] = 108 + 59,5 + 47,5 + 6 + 4,5 + 4,5 + 4,5 + 40 + 80,5 + 99 = 454,0<\/em><\/p>\n<p>D\u00e9terminons d\u00e9sormais les coefficients de la droite de r\u00e9gression :<\/p>\n<p><em>b<sub>1<\/sub> = 454,0 \/ 82,5 &asymp; 5,50303<\/em><\/p>\n<p><em>b<sub>0<\/sub> = 59,0 &minus; (5,50303 &times; 6,5) = 59,0 &minus; 35,7697 = 23,2303<\/em><\/p>\n<p>L&rsquo;\u00e9quation pr\u00e9dictive ajust\u00e9e s&rsquo;\u00e9tablit par cons\u00e9quent selon la formule : <em>\u0177<sub>i<\/sub> = 23,2303 + 5,50303 &times; x<sub>i<\/sub><\/em>. Nous pouvons d\u00e9sormais calculer rigoureusement pour chaque participant sa pr\u00e9diction, son r\u00e9sidu ordinaire ainsi que le carr\u00e9 de ce dernier.<\/p>\n<table border=\"1\" cellpadding=\"6\" cellspacing=\"0\" style=\"border-collapse: collapse;width: 100%;text-align: right\">\n<thead style=\"background-color: #f2f2f2\">\n<tr>\n<th style=\"text-align: center\">Participant (<em>i<\/em>)<\/th>\n<th style=\"text-align: center\">Stress (<em>x<sub>i<\/sub><\/em>)<\/th>\n<th style=\"text-align: center\">Temps obs. (<em>y<sub>i<\/sub><\/em>)<\/th>\n<th style=\"text-align: center\">Temps pr\u00e9dit (<em>\u0177<sub>i<\/sub><\/em>)<\/th>\n<th style=\"text-align: center\">R\u00e9sidu brut (<em>e<sub>i<\/sub><\/em>)<\/th>\n<th style=\"text-align: center\">R\u00e9sidu au carr\u00e9 (<em>e<sub>i<\/sub><sup>2<\/sup><\/em>)<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td style=\"text-align: center\">1<\/td>\n<td>2<\/td>\n<td>35<\/td>\n<td>34,236<\/td>\n<td>+0,764<\/td>\n<td>0,584<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">2<\/td>\n<td>3<\/td>\n<td>42<\/td>\n<td>39,739<\/td>\n<td>+2,261<\/td>\n<td>5,112<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">3<\/td>\n<td>4<\/td>\n<td>40<\/td>\n<td>45,242<\/td>\n<td>-5,242<\/td>\n<td>27,479<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">4<\/td>\n<td>5<\/td>\n<td>55<\/td>\n<td>50,745<\/td>\n<td>+4,255<\/td>\n<td>18,105<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">5<\/td>\n<td>6<\/td>\n<td>50<\/td>\n<td>56,248<\/td>\n<td>-6,248<\/td>\n<td>39,038<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">6<\/td>\n<td>7<\/td>\n<td>68<\/td>\n<td>61,752<\/td>\n<td>+6,248<\/td>\n<td>39,038<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">7<\/td>\n<td>8<\/td>\n<td>62<\/td>\n<td>67,255<\/td>\n<td>-5,255<\/td>\n<td>27,615<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">8<\/td>\n<td>9<\/td>\n<td>75<\/td>\n<td>72,758<\/td>\n<td>+2,242<\/td>\n<td>5,027<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">9<\/td>\n<td>10<\/td>\n<td>82<\/td>\n<td>78,261<\/td>\n<td>+3,739<\/td>\n<td>13,980<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center\">10<\/td>\n<td>11<\/td>\n<td>81<\/td>\n<td>83,764<\/td>\n<td>-2,764<\/td>\n<td>7,640<\/td>\n<\/tr>\n<\/tbody>\n<tfoot style=\"font-weight: bold;background-color: #fafafa\">\n<tr>\n<td style=\"text-align: center\">Total \/ Moyenne<\/td>\n<td>x\u0304 = 6,5<\/td>\n<td>\u0233 = 59,0<\/td>\n<td>59,000<\/td>\n<td>&sum; e<sub>i<\/sub> = 0,000<\/td>\n<td>SSE = 183,818<\/td>\n<\/tr>\n<\/tfoot>\n<\/table>\n<p>Ce tableau d\u00e9montre empiriquement la confirmation de la propri\u00e9t\u00e9 th\u00e9orique \u00e9nonc\u00e9e : la sommation alg\u00e9brique pr\u00e9cise de l&rsquo;ensemble des r\u00e9sidus ordinaires donne exactement 0,000. Par ailleurs, la somme des carr\u00e9s r\u00e9siduels (<em>SSE<\/em>) s&rsquo;\u00e9l\u00e8ve \u00e0 183,818. Consid\u00e9rant que les degr\u00e9s de libert\u00e9 r\u00e9siduels s&rsquo;\u00e9l\u00e8vent ici \u00e0 <em>n &minus; k &minus; 1 = 10 &minus; 1 &minus; 1 = 8<\/em>, nous pouvons d\u00e9duire la variance d&rsquo;erreur r\u00e9siduelle non biais\u00e9e : <em>s<sub>e<\/sub><sup>2<\/sup> = 183,818 \/ 8 = 22,977<\/em>, ce qui produit un \u00e9cart-type r\u00e9siduel (erreur type de l&rsquo;estimation) de <em>s<sub>e<\/sub> = &radic;22,977 &asymp; 4,793<\/em> centi\u00e8mes de seconde.<\/p>\n<h2>4. Typologie avanc\u00e9e des r\u00e9sidus : au-del\u00e0 des r\u00e9sidus bruts<\/h2>\n<h3>4.1 Les r\u00e9sidus standardis\u00e9s (ou r\u00e9sidus de Pearson)<\/h3>\n<p>Si les r\u00e9sidus ordinaires constituent l&rsquo;expression la plus imm\u00e9diate de l&rsquo;erreur d&rsquo;ajustement, ils pr\u00e9sentent une limite m\u00e9thodologique majeure : ils sont intrins\u00e8quement d\u00e9pendants des unit\u00e9s de mesure d&rsquo;origine associ\u00e9es \u00e0 la variable d\u00e9pendante. Comparer l&rsquo;amplitude brute d&rsquo;un r\u00e9sidu issu d&rsquo;un temps de r\u00e9action en millisecondes avec celle d&rsquo;un r\u00e9sidu issu d&rsquo;un score sur une \u00e9chelle de d\u00e9pression en 20 points s&rsquo;av\u00e8re rigoureusement impossible. Pour surmonter cette entrave m\u00e9trique, la litt\u00e9rature statistique recourt aux <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Standardized_residual\">r\u00e9sidus standardis\u00e9s<\/a>, souvent d\u00e9sign\u00e9s sous l&rsquo;appellation de r\u00e9sidus de Pearson.<\/p>\n<p>La formulation alg\u00e9brique d&rsquo;un r\u00e9sidu standardis\u00e9 (not\u00e9 <em>z<sub>i<\/sub><\/em>) consiste \u00e0 diviser le r\u00e9sidu brut par l&rsquo;estimation globale de l&rsquo;\u00e9cart-type r\u00e9siduel de l&rsquo;\u00e9chantillon :<\/p>\n<p><em>z<sub>i<\/sub> = e<sub>i<\/sub> \/ s<sub>e<\/sub> = (y<sub>i<\/sub> &minus; \u0177<sub>i<\/sub>) \/ &radic;MSE<\/em><\/p>\n<p>Cette op\u00e9ration r\u00e9alise une normalisation m\u00e9trique de la distribution d&rsquo;erreur, conf\u00e9rant aux r\u00e9sidus une \u00e9chelle adimensionnelle caract\u00e9ris\u00e9e approximativement par une moyenne de z\u00e9ro et un \u00e9cart-type unitaire. Gr\u00e2ce \u00e0 cette transformation, le statisticien peut \u00e9valuer l&rsquo;\u00e9cart relatif d&rsquo;une observation \u00e0 la droite ind\u00e9pendamment de l&rsquo;\u00e9chelle d&rsquo;origine. En psychom\u00e9trie computationnelle, les seuils d\u00e9cisionnels conventionnels appliqu\u00e9s aux r\u00e9sidus standardis\u00e9s s&rsquo;alignent sur les propri\u00e9t\u00e9s de la loi normale centr\u00e9e r\u00e9duite : un r\u00e9sidu standardis\u00e9 dont la valeur absolue d\u00e9passe 2 (|<em>z<sub>i<\/sub><\/em>| &gt; 2) signale une d\u00e9viation inhabituelle pr\u00e9sente chez moins de 5 % des individus th\u00e9oriques, tandis qu&rsquo;une valeur absolue exc\u00e9dant 3 (|<em>z<sub>i<\/sub><\/em>| &gt; 3) identifie une observation suspecte m\u00e9ritant un audit minutieux.<\/p>\n<h3>4.2 Les r\u00e9sidus studentis\u00e9s internes et externes<\/h3>\n<p>La standardisation simple par <em>s<sub>e<\/sub><\/em> souffre d&rsquo;un biais th\u00e9orique sous-jacent : elle postule implicitement que la variance de chaque r\u00e9sidu individuel est constante sur l&rsquo;ensemble du domaine de pr\u00e9diction. Or, les d\u00e9veloppements math\u00e9matiques d\u00e9montrent que la pr\u00e9cision pr\u00e9dictive d&rsquo;une droite de r\u00e9gression des moindres carr\u00e9s est maximale au niveau du centre de gravit\u00e9 des donn\u00e9es (la moyenne des pr\u00e9dicteurs) et se d\u00e9grade au fur et \u00e0 mesure que l&rsquo;on s&rsquo;\u00e9carte vers les extr\u00e9mit\u00e9s de la distribution du pr\u00e9dicteur. En cons\u00e9quence, la variance r\u00e9elle d&rsquo;un r\u00e9sidu <em>e<sub>i<\/sub><\/em> d\u00e9pend formellement de la position relative de son observation sur l&rsquo;axe des <em>X<\/em>, grandeur formalis\u00e9e par le levier (<em>leverage<\/em>, not\u00e9 <em>h<sub>ii<\/sub><\/em>).<\/p>\n<p>La variance exacte du <em>i<\/em>-i\u00e8me r\u00e9sidu brut s&rsquo;\u00e9crit de mani\u00e8re formelle : <em>Var(e<sub>i<\/sub>) = &sigma;<sup>2<\/sup>(1 &minus; h<sub>ii<\/sub>)<\/em>. Pour tenir compte de cette h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 structurelle, le r\u00e9sidu studentis\u00e9 interne (souvent appel\u00e9 <em>internally studentized residual<\/em>) ajuste chaque \u00e9cart en le divisant par son erreur type locale :<\/p>\n<p><em>r<sub>i<\/sub> = e<sub>i<\/sub> \/ [s<sub>e<\/sub> &radic;(1 &minus; h<sub>ii<\/sub>)]<\/em><\/p>\n<p>Bien que cette m\u00e9trique apporte un raffinement notable, elle demeure sensible au fait qu&rsquo;une observation atypique extr\u00eame exerce une influence d\u00e9mesur\u00e9e sur l&rsquo;estimation globale de <em>s<sub>e<\/sub><\/em>, tendant \u00e0 gonfler artificiellement ce d\u00e9nominateur et \u00e0 masquer le caract\u00e8re divergent de ce point pr\u00e9cis (ph\u00e9nom\u00e8ne de masquage ou <em>masking<\/em>). Pour pallier ce biais, les statisticiens ont d\u00e9velopp\u00e9 le concept de r\u00e9sidu studentis\u00e9 externe, fr\u00e9quemment appel\u00e9 r\u00e9sidu supprim\u00e9 (<em>jackknife<\/em> ou <em>externally studentized residual<\/em>, not\u00e9 <em>t<sub>i<\/sub><\/em>). Dans cette formulation, l&rsquo;erreur type r\u00e9siduelle (not\u00e9e <em>s<sub>(i)<\/sub><\/em>) est recalcul\u00e9e en excluant explicitement la <em>i<\/em>-i\u00e8me observation de l&rsquo;\u00e9chantillon d&rsquo;estimation :<\/p>\n<p><em>t<sub>i<\/sub> = e<sub>i<\/sub> \/ [s<sub>(i)<\/sub> &radic;(1 &minus; h<sub>ii<\/sub>)]<\/em><\/p>\n<p>Le r\u00e9sidu studentis\u00e9 externe suit rigoureusement une loi de Student \u00e0 <em>n &minus; k &minus; 2<\/em> degr\u00e9s de libert\u00e9. Cette propri\u00e9t\u00e9 offre un avantage m\u00e9thodologique d\u00e9cisif pour d\u00e9tecter les cas aberrants : elle autorise l&rsquo;application de tests d&rsquo;hypoth\u00e8ses formels avec ajustement pour comparaisons multiples (par exemple via la correction de Bonferroni), fournissant un crit\u00e8re d&rsquo;exclusion statistique robuste et th\u00e9oriquement fond\u00e9.<\/p>\n<h3>4.3 Les r\u00e9sidus partiels et r\u00e9sidus de r\u00e9gression partielle<\/h3>\n<p>Dans les configurations d&rsquo;analyse impliquant plusieurs r\u00e9gresseurs simultan\u00e9s (r\u00e9gression multiple), l&rsquo;examen des r\u00e9sidus ordinaires globaux devient insuffisant pour d\u00e9m\u00ealer les relations sp\u00e9cifiques unissant la variable crit\u00e8re \u00e0 un pr\u00e9dicteur unique isol\u00e9. C&rsquo;est ici qu&rsquo;interviennent les concepts de r\u00e9sidus partiels et de graphiques de r\u00e9gression partielle (souvent d\u00e9nomm\u00e9s <em>Added Variable Plots<\/em> ou graphiques de variables ajout\u00e9es).<\/p>\n<p>Le r\u00e9sidu partiel d\u00e9di\u00e9 au <em>j<\/em>-i\u00e8me pr\u00e9dicteur (not\u00e9 <em>e<sub>i(j)<\/sub><\/em>) est construit en additionnant la composante de l&rsquo;effet lin\u00e9aire mod\u00e9lis\u00e9 de ce pr\u00e9dicteur sp\u00e9cifique au r\u00e9sidu global du mod\u00e8le complet :<\/p>\n<p><em>e<sub>i(j)<\/sub> = e<sub>i<\/sub> + b<sub>j<\/sub> x<sub>ij<\/sub><\/em><\/p>\n<p>Sur le plan g\u00e9om\u00e9trique, tracer le nuage de points reliant la variable <em>x<sub>j<\/sub><\/em> \u00e0 ces r\u00e9sidus partiels permet de visualiser la trajectoire non fauss\u00e9e du lien entre ce r\u00e9gresseur et le crit\u00e8re, une fois le signal des autres variables explicatives retir\u00e9. Cette approche graphique s&rsquo;av\u00e8re indispensable pour d\u00e9masquer d&rsquo;\u00e9ventuelles relations curvilin\u00e9aires ou des seuils de rupture masqu\u00e9s par la colin\u00e9arit\u00e9 statistique dans les investigations psychologiques complexes portant sur des covariables entrem\u00eal\u00e9es (par exemple, isoler l&rsquo;impact net de l&rsquo;anxi\u00e9t\u00e9 sur les performances de rappel mn\u00e9sique apr\u00e8s neutralisation math\u00e9matique de l&rsquo;\u00e2ge et de la vitesse de traitement de base).<\/p>\n<h2>5. V\u00e9rification des postulats de la r\u00e9gression par l&rsquo;examen des r\u00e9sidus<\/h2>\n<h3>5.1 Postulat de lin\u00e9arit\u00e9 de la relation<\/h3>\n<p>L&rsquo;inf\u00e9rence statistique classique sous-jacente \u00e0 la r\u00e9gression des moindres carr\u00e9s postule de mani\u00e8re irr\u00e9ductible que la relation reliant les variables explicatives \u00e0 l&rsquo;esp\u00e9rance math\u00e9matique de la variable crit\u00e8re est de nature strictement lin\u00e9aire. Une sp\u00e9cification erron\u00e9e de cette morphologie fonctionnelle invalide d&#8217;embl\u00e9e l&rsquo;int\u00e9gralit\u00e9 des estimations des pentes. Pour auditer la validit\u00e9 de ce postulat, le statisticien s&rsquo;appuie sur le graphique fondamental des r\u00e9sidus en fonction des valeurs pr\u00e9dites (g\u00e9n\u00e9ralement d\u00e9sign\u00e9 dans la litt\u00e9rature sous le vocable anglophone de diagramme <em>Residuals vs Fitted<\/em>).<\/p>\n<p>Sous l&rsquo;hypoth\u00e8se de lin\u00e9arit\u00e9 parfaite, le nuage de points de ce graphique doit pr\u00e9senter l&rsquo;apparence d&rsquo;une bande horizontale homog\u00e8ne, exempte de toute organisation spatiale pr\u00e9visible. Les points doivent se distribuer de fa\u00e7on al\u00e9atoire et sym\u00e9trique de part et d&rsquo;autre de la ligne horizontale de r\u00e9f\u00e9rence situ\u00e9e \u00e0 l&rsquo;ordonn\u00e9e z\u00e9ro.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/statology_residuals2.jpg\" alt=\"Trend line on scatterplot in Excel\" \/><figcaption>Trend line on scatterplot in Excel<\/figcaption><\/figure>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;\u00e9mergence d&rsquo;un motif curviligne discernable &ndash; telle qu&rsquo;une parabole en forme de U, une configuration en U invers\u00e9 ou une oscillation sinuso\u00efdale &ndash; r\u00e9v\u00e8le que le mod\u00e8le a \u00e9chou\u00e9 \u00e0 capter une courbure intrins\u00e8que aux donn\u00e9es comportementales. Par exemple, la c\u00e9l\u00e8bre loi de Yerkes-Dodson reliant l&rsquo;\u00e9veil motivationnel \u00e0 la performance cognitive d\u00e9crit typiquement une relation non lin\u00e9aire en U invers\u00e9. Tenter de capturer ce ph\u00e9nom\u00e8ne par une \u00e9quation lin\u00e9aire simple engendre une concentration massive de r\u00e9sidus positifs aux valeurs interm\u00e9diaires et de r\u00e9sidus n\u00e9gatifs aux extr\u00eames. La rem\u00e9diation m\u00e9thodique impose alors l&rsquo;introduction de termes polynomiaux (tels que <em>X<sup>2<\/sup><\/em>) ou l&rsquo;application de transformations non lin\u00e9aires (logarithmiques, inverses) sur les pr\u00e9dicteurs.<\/p>\n<h3>5.2 Postulat d&rsquo;homosc\u00e9dasticit\u00e9 (constance de la variance de l&rsquo;erreur)<\/h3>\n<p>Le postulat d&rsquo;<a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Homosc%C3%A9dasticit%C3%A9\">homosc\u00e9dasticit\u00e9<\/a> stipule que la variance conditionnelle des termes d&rsquo;erreur r\u00e9siduelle demeure rigoureusement constante \u00e0 travers tous les niveaux de pr\u00e9diction du mod\u00e8le. Formellement, cette condition s&rsquo;\u00e9nonce : <em>Var(e<sub>i<\/sub> | X) = &sigma;<sup>2<\/sup><\/em> pour tout <em>i<\/em>. Lorsque cette hypoth\u00e8se vacille, les donn\u00e9es sont frapp\u00e9es d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9. Bien que les estimateurs des pentes MCO demeurent sans biais en pr\u00e9sence de variances in\u00e9gales, leurs erreurs types associ\u00e9es deviennent gravement inexactes, faussant les calculs de significativit\u00e9 des tests <em>t<\/em> et les intervalles de confiance correspondants.<\/p>\n<p>L&rsquo;inspection visuelle du graphique des r\u00e9sidus studentis\u00e9s en fonction des valeurs pr\u00e9dites (ou du graphique dit <em>Scale-Location<\/em>, exploitant la racine carr\u00e9e des r\u00e9sidus standardis\u00e9s) fournit la premi\u00e8re ligne diagnostique. Une manifestation embl\u00e9matique de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 r\u00e9side dans l&rsquo;apparition d&rsquo;un patron d&rsquo;\u00e9vasement g\u00e9om\u00e9trique en c\u00f4ne, en m\u00e9gaphone ou en entonnoir : les r\u00e9sidus, initialement tr\u00e8s resserr\u00e9s pour les pr\u00e9dictions faibles, voient leur dispersion spatiale s&rsquo;amplifier consid\u00e9rablement \u00e0 mesure que les scores pr\u00e9dits progressent.<\/p>\n<p>Cette observation subjective doit imp\u00e9rativement \u00eatre corrobor\u00e9e par des tests statistiques inf\u00e9rentiels rigoureux :<\/p>\n<ul>\n<li><strong>Le test de Breusch-Pagan :<\/strong> particuli\u00e8rement performant pour d\u00e9tecter les h\u00e9t\u00e9rosc\u00e9dasticit\u00e9s de nature lin\u00e9aire en mod\u00e9lisant les carr\u00e9s des r\u00e9sidus en fonction des r\u00e9gresseurs d&rsquo;origine.<\/li>\n<li><strong>Le test de White :<\/strong> un test conservateur g\u00e9n\u00e9ral ne n\u00e9cessitant aucun a priori sur la structure de la variance et int\u00e9grant les carr\u00e9s et produits crois\u00e9s des pr\u00e9dicteurs.<\/li>\n<li><strong>Le test de Koenker :<\/strong> une variante robuste du test de Breusch-Pagan, particuli\u00e8rement recommand\u00e9e en psychom\u00e9trie lorsque le postulat sous-jacent de normalit\u00e9 bivari\u00e9e ne peut \u00eatre pleinement garanti.<\/li>\n<\/ul>\n<h3>5.3 Postulat de normalit\u00e9 de la distribution des erreurs<\/h3>\n<p>Pour assurer la validit\u00e9 de l&rsquo;inf\u00e9rence statistique aux petits et moyens \u00e9chantillons, les termes d&rsquo;erreur r\u00e9siduelle doivent se distribuer de mani\u00e8re congruente avec une distribution gaussienne th\u00e9orique, condition formalis\u00e9e par <em>e<sub>i<\/sub> ~ N(0, &sigma;<sup>2<\/sup>)<\/em>. Si le th\u00e9or\u00e8me central limite garantit une convergence asymptotique de la distribution d&rsquo;\u00e9chantillonnage des coefficients vers la normalit\u00e9 dans les cohortes volumineuses, les recherches psychologiques men\u00e9es sur des \u00e9chantillons restreints requi\u00e8rent un examen approfondi de la normalit\u00e9 des r\u00e9sidus.<\/p>\n<p>L&rsquo;outil s\u00e9miologique privil\u00e9gi\u00e9 pour cette \u00e9valuation r\u00e9side dans le diagramme quantile-quantile normal (couramment baptis\u00e9 <em>Normal Q-Q Plot<\/em>). Ce graphique projette les quantiles observ\u00e9s de la distribution des r\u00e9sidus standardis\u00e9s contre les quantiles th\u00e9oriques attendus sous une loi normale id\u00e9ale. Lorsque le postulat est v\u00e9rifi\u00e9, les observations empiriques s&rsquo;alignent \u00e9troitement sur la bissectrice diagonale \u00e0 45 degr\u00e9s.<\/p>\n<p>Les d\u00e9viations caract\u00e9ristiques \u00e0 cette trajectoire id\u00e9ale permettent de poser un diagnostic morphologique pr\u00e9cis :<\/p>\n<ul>\n<li><strong>Distributions \u00e0 queues lourdes (leptokurtiques) :<\/strong> les points s&rsquo;\u00e9cartent vers le haut \u00e0 l&rsquo;extr\u00e9mit\u00e9 droite et vers le bas \u00e0 l&rsquo;extr\u00e9mit\u00e9 gauche, formant une silhouette en &lsquo;S&rsquo; prononc\u00e9e et signalant une abondance anormale de valeurs atypiques.<\/li>\n<li><strong>Distributions asym\u00e9triques (skewness positive ou n\u00e9gative) :<\/strong> le nuage pr\u00e9sente une convexit\u00e9 unilat\u00e9rale continue par rapport \u00e0 la diagonale de r\u00e9f\u00e9rence, r\u00e9v\u00e9lant un plancher ou un plafond de mesure dans l&rsquo;\u00e9valuation psychologique.<\/li>\n<\/ul>\n<p>Ces \u00e9l\u00e9ments qualitatifs s&rsquo;accompagnent de v\u00e9rifications formelles via le test de normalit\u00e9 de Shapiro-Wilk (adapt\u00e9 aux \u00e9chantillons o\u00f9 <em>n &lt; 50<\/em>) ou le test de Kolmogorov-Smirnov assorti de la correction syst\u00e9matique de Lilliefors pour les \u00e9chantillons plus amples.<\/p>\n<h3>5.4 Postulat d&rsquo;ind\u00e9pendance des r\u00e9sidus (absence d&rsquo;autocorr\u00e9lation)<\/h3>\n<p>Le th\u00e9or\u00e8me de Gauss-Markov exige l&rsquo;ind\u00e9pendance statistique mutuelle de toutes les paires d&rsquo;erreurs r\u00e9siduelles : la covariance entre le r\u00e9sidu de l&rsquo;individu <em>i<\/em> et celui de l&rsquo;individu <em>j<\/em> doit \u00eatre strictement nulle (<em>Cov(e<sub>i<\/sub>, e<sub>j<\/sub>) = 0<\/em> pour tout <em>i &ne; j<\/em>). Cette condition est fr\u00e9quemment compromise d\u00e8s lors que le protocole de recueil introduit une s\u00e9quentialit\u00e9 temporelle (mesures r\u00e9p\u00e9t\u00e9es longitudinales, protocoles d&rsquo;\u00e9valuation \u00e9cologique momentan\u00e9e &#8211; EMA) ou une proximit\u00e9 spatiale et groupale (\u00e9l\u00e8ves issus des m\u00eames classes scolaires, patients pris en charge par le m\u00eame th\u00e9rapeute).<\/p>\n<p>La persistance d&rsquo;une d\u00e9pendance s\u00e9rielle d&rsquo;ordre un &ndash; o\u00f9 le r\u00e9sidu d&rsquo;une observation au temps <em>t<\/em> pr\u00e9dit partiellement le r\u00e9sidu au temps <em>t+1<\/em> &ndash; fausse gravement l&rsquo;\u00e9valuation de la variance d&rsquo;\u00e9chantillonnage des estimateurs. Une autocorr\u00e9lation positive comprime artificiellement les erreurs types calcul\u00e9es, conf\u00e9rant aux tests <em>t<\/em> de Student une confiance fallacieuse et conduisant \u00e0 l&rsquo;affirmation indue d&rsquo;effets th\u00e9rapeutiques ou comportementaux non significatifs dans la population cible.<\/p>\n<p>L&rsquo;indicateur classique d&rsquo;\u00e9valuation de cette propri\u00e9t\u00e9 est le <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Test_de_Durbin-Watson\">test de Durbin-Watson<\/a>. La statistique <em>d<\/em> qui en d\u00e9coule est calcul\u00e9e selon l&rsquo;algorithme :<\/p>\n<p><em>d = [&sum;<sub>t=2<\/sub><sup>n<\/sup> (e<sub>t<\/sub> &minus; e<sub>t&minus;1<\/sub>)<sup>2<\/sup>] \/ [&sum;<sub>t=1<\/sub><sup>n<\/sup> e<sub>t<\/sub><sup>2<\/sup>]<\/em><\/p>\n<p>La m\u00e9trique <em>d<\/em> oscille sur un continuum born\u00e9 entre 0 et 4. Une valeur proche de 2,0 indique une ind\u00e9pendance s\u00e9rielle satisfaisante des r\u00e9sidus. \u00c0 l&rsquo;oppos\u00e9, une valeur descendant nettement en de\u00e7\u00e0 de 1,5 d\u00e9nonce une autocorr\u00e9lation positive critique, tandis qu&rsquo;une valeur exc\u00e9dant 2,5 traduit une autocorr\u00e9lation s\u00e9rielle n\u00e9gative. Face \u00e0 de telles configurations, le cadre de la r\u00e9gression lin\u00e9aire simple MCO doit \u00eatre abandonn\u00e9 au profit de mod\u00e8les autor\u00e9gressifs int\u00e9gr\u00e9s (ARIMA) ou de mod\u00e9lisations mixtes lin\u00e9aires int\u00e9grant des matrices de covariance r\u00e9siduelle structur\u00e9es (autor\u00e9gressives d&rsquo;ordre 1, toeplitz ou non structur\u00e9es).<\/p>\n<h2>6. D\u00e9tection des donn\u00e9es aberrantes et des observations influentes<\/h2>\n<h3>6.1 Distinction th\u00e9orique entre valeur aberrante, levier et influence<\/h3>\n<p>Dans la litt\u00e9rature m\u00e9thodologique consacr\u00e9e au diagnostic de r\u00e9gression, il convient de distinguer avec rigueur trois dimensions topologiques distinctes : la valeur aberrante (<em>outlier<\/em>), le point de levier (<em>leverage<\/em>) et l&rsquo;observation influente (<em>influential observation<\/em>). Confondre ces terminologies expose le chercheur \u00e0 des interpr\u00e9tations partielles et \u00e0 des d\u00e9cisions d&rsquo;exclusion d&rsquo;observations injustifi\u00e9es.<\/p>\n<p>Une valeur aberrante r\u00e9siduelle qualifie une observation qui pr\u00e9sente un \u00e9cart vertical majeur par rapport \u00e0 la trajectoire mod\u00e9lis\u00e9e. Il s&rsquo;agit d&rsquo;un point dont la valeur empirique <em>y<sub>i<\/sub><\/em> s&rsquo;av\u00e8re profond\u00e9ment discordante par rapport \u00e0 sa pr\u00e9diction <em>\u0177<sub>i<\/sub><\/em>, g\u00e9n\u00e9rant un r\u00e9sidu standardis\u00e9 ou studentis\u00e9 de magnitude extr\u00eame (|<em>t<sub>i<\/sub><\/em>| &gt; 3), bien que ses valeurs sur les variables pr\u00e9dictives <em>X<\/em> puissent se situer dans la moyenne de la population.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/statology_residuals2.jpg\" alt=\"Trend line on scatterplot in Excel\" \/><figcaption>Trend line on scatterplot in Excel<\/figcaption><\/figure>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;effet de levier mesure l&rsquo;\u00e9loignement d&rsquo;une observation sur l&rsquo;axe horizontal constitu\u00e9 par l&rsquo;espace multidimensionnel des pr\u00e9dicteurs, ind\u00e9pendamment de sa coordonn\u00e9e sur la variable crit\u00e8re <em>Y<\/em>. Formalis\u00e9 par l&rsquo;\u00e9l\u00e9ment diagonal de la matrice de projection chapeau <em>h<sub>ii<\/sub><\/em>, le levier traduit la distance de Mahalanobis d&rsquo;un individu par rapport au centro\u00efde des variables ind\u00e9pendantes. Une observation peut afficher un levier consid\u00e9rable sans pour autant fausser le mod\u00e8le si elle s&rsquo;aligne fid\u00e8lement sur l&rsquo;axe th\u00e9orique naturel de la r\u00e9gression.<\/p>\n<p>L&rsquo;observation influente, enfin, repr\u00e9sente l&rsquo;intersection dynamique de ces deux propri\u00e9t\u00e9s. Un cas est qualifi\u00e9 d&rsquo;influent d\u00e8s lors que sa pr\u00e9sence &ndash; ou son omission d\u00e9lib\u00e9r\u00e9e &ndash; modifie substantiellement les coefficients de r\u00e9gression estim\u00e9s, la significativit\u00e9 statistique des tests d&rsquo;hypoth\u00e8ses ou la qualit\u00e9 de pr\u00e9diction globale du syst\u00e8me. Une observation influente conjugue presque toujours un effet de levier mod\u00e9r\u00e9 \u00e0 substantiel avec un r\u00e9sidu de grande ampleur, agissant tel un point d&rsquo;appui physique faisant pivoter artificiellement la droite de r\u00e9gression dans sa direction.<\/p>\n<h3>6.2 La distance de Cook : formulation et crit\u00e8res d&rsquo;application<\/h3>\n<p>L&rsquo;indicateur le plus r\u00e9put\u00e9 pour \u00e9valuer l&rsquo;influence exerc\u00e9e par une observation singuli\u00e8re est la <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Cook%27s_distance\">distance de Cook<\/a> (not\u00e9e <em>D<sub>i<\/sub><\/em>). Con\u00e7ue par le statisticien R. Dennis Cook, cette m\u00e9trique synth\u00e9tique quantifie le d\u00e9placement multidimensionnel global op\u00e9r\u00e9 sur l&rsquo;ensemble du vecteur des pr\u00e9dictions lorsque le mod\u00e8le est r\u00e9estim\u00e9 en soustrayant l&rsquo;observation d&rsquo;int\u00e9r\u00eat de l&rsquo;ensemble de donn\u00e9es.<\/p>\n<p>Sur le plan computationnel, la distance de Cook se formule sans imposer le r\u00e9ajustement s\u00e9quentiel de <em>n<\/em> r\u00e9gressions distinctes, gr\u00e2ce \u00e0 une relation combinant \u00e9l\u00e9gamment le r\u00e9sidu studentis\u00e9 interne <em>r<sub>i<\/sub><\/em> et le coefficient de levier <em>h<sub>ii<\/sub><\/em> :<\/p>\n<p><em>D<sub>i<\/sub> = [e<sub>i<\/sub><sup>2<\/sup> \/ ((k + 1) &times; MSE)] &times; [h<sub>ii<\/sub> \/ (1 &minus; h<sub>ii<\/sub>)<sup>2<\/sup>] = [r<sub>i<\/sub><sup>2<\/sup> \/ (k + 1)] &times; [h<sub>ii<\/sub> \/ (1 &minus; h<sub>ii<\/sub>)]<\/em><\/p>\n<p>Cette \u00e9criture d\u00e9montre que la distance de Cook est le produit direct de deux facteurs : la taille relative du r\u00e9sidu d&rsquo;ajustement (premi\u00e8re parenth\u00e8se) et la susceptibilit\u00e9 g\u00e9om\u00e9trique du mod\u00e8le aux perturbations structurelles li\u00e9e au levier (seconde parenth\u00e8se).<\/p>\n<p>Quant aux seuils d\u00e9cisionnels r\u00e9gissant l&rsquo;interpr\u00e9tation de <em>D<sub>i<\/sub><\/em>, deux approches coexistent dans la communaut\u00e9 scientifique :<\/p>\n<ul>\n<li><strong>Le seuil absolu classique :<\/strong> une valeur de <em>D<sub>i<\/sub> &gt; 1,0<\/em> signale une observation influente majeure, susceptible de justifier une r\u00e9\u00e9valuation de l&rsquo;estimation. Certains auteurs \u00e9tendent cette alerte d\u00e8s <em>D<sub>i<\/sub> &gt; 0,5<\/em>.<\/li>\n<li><strong>Le seuil relatif aux cohortes comportementales :<\/strong> dans les \u00e9tudes o\u00f9 la taille de l&rsquo;\u00e9chantillon est plus modeste, la borne critique ajust\u00e9e <em>D<sub>i<\/sub> &gt; 4 \/ n<\/em> (ou <em>4 \/ (n &minus; k &minus; 1)<\/em>) est pr\u00e9conis\u00e9e. Elle permet d&rsquo;identifier les individus dont l&rsquo;influence est disproportionn\u00e9e par rapport \u00e0 la moyenne du groupe \u00e9tudi\u00e9.<\/li>\n<\/ul>\n<h3>6.3 Les indices DFFITS et DFBETAS<\/h3>\n<p>Pour affiner l&rsquo;\u00e9valuation diagnostique au-del\u00e0 de la mesure globale fournie par la distance de Cook, le psychom\u00e9tricien dispose d&rsquo;outils analytiques individualis\u00e9s : les indices DFBETAS et DFFITS. Ces indices d\u00e9composent l&rsquo;impact de chaque observation cas par cas.<\/p>\n<p>L&rsquo;indice DFBETAS \u00e9value la d\u00e9viation sp\u00e9cifique imprim\u00e9e par l&rsquo;observation <em>i<\/em> sur un coefficient de pente <em>b<sub>j<\/sub><\/em> pris isol\u00e9ment. Exprim\u00e9 en unit\u00e9s d&rsquo;erreur type, il s&rsquo;\u00e9nonce selon le rapport :<\/p>\n<p><em>DFBETAS<sub>j(i)<\/sub> = (b<sub>j<\/sub> &minus; b<sub>j(i)<\/sub>) \/ [s<sub>(i)<\/sub> &radic;(&sum; (x<sub>ij<\/sub> &minus; x\u0304<sub>j<\/sub>)<sup>2<\/sup>)]<\/em><\/p>\n<p>Dans cette formule, <em>b<sub>j(i)<\/sub><\/em> repr\u00e9sente la valeur r\u00e9estim\u00e9e du param\u00e8tre de r\u00e9gression apr\u00e8s \u00e9viction compl\u00e8te du <em>i<\/em>-i\u00e8me cas. Un seuil d\u00e9cisionnel conventionnel stipule qu&rsquo;une valeur absolue |DFBETAS<sub>j(i)<\/sub>| exc\u00e9dant l&rsquo;indice <em>2 \/ &radic;n<\/em> d\u00e9signe un individu modifiant de mani\u00e8re disproportionn\u00e9e la pente de la variable explicative consid\u00e9r\u00e9e, risquant d&rsquo;en alt\u00e9rer la signification clinique.<\/p>\n<p>De fa\u00e7on analogue, l&rsquo;indice DFFITS mesure le nombre d&rsquo;\u00e9carts-types dont la valeur pr\u00e9dite pour le sujet <em>i<\/em> lui-m\u00eame oscille selon que ce cas particulier est int\u00e9gr\u00e9 ou retir\u00e9 de la matrice d&rsquo;apprentissage :<\/p>\n<p><em>DFFITS<sub>i<\/sub> = (\u0177<sub>i<\/sub> &minus; \u0177<sub>i(i)<\/sub>) \/ [s<sub>(i)<\/sub> &radic;h<sub>ii<\/sub>] = t<sub>i<\/sub> &radic;[h<sub>ii<\/sub> \/ (1 &minus; h<sub>ii<\/sub>)]<\/em><\/p>\n<p>Un seuil de coupure empirique \u00e9tabli \u00e0 |DFFITS<sub>i<\/sub>| &gt; 2 &radic;[(k + 1) \/ n] signale un cas alt\u00e9rant sa propre pr\u00e9diction. Face \u00e0 ces diagnostics, la r\u00e8gle \u00e9thique proscrit toute suppression machinale des observations divergentes : chaque cas rep\u00e9r\u00e9 doit d&rsquo;abord faire l&rsquo;objet d&rsquo;une tra\u00e7abilit\u00e9 archivistique (erreur d&rsquo;encodage, non-respect av\u00e9r\u00e9 des consignes du test), puis \u00eatre analys\u00e9 dans le cadre d&rsquo;analyses de sensibilit\u00e9 avec et sans l&rsquo;observation litigieuse.<\/p>\n<h2>7. Calcul et diagnostic des r\u00e9sidus sous R pour la recherche empirique<\/h2>\n<h3>7.1 Extraction des r\u00e9sidus avec les fonctions de base<\/h3>\n<p>Le langage et environnement statistique open-source R propose une architecture native pour l&rsquo;estimation et l&rsquo;extraction des composantes r\u00e9siduelles. L&rsquo;ajustement initial d&rsquo;une r\u00e9gression s&rsquo;effectue via la fonction g\u00e9n\u00e9rique <code>lm()<\/code> (<em>linear model<\/em>), qui configure la structure du syst\u00e8me \u00e0 partir d&rsquo;une formule symbolique et d&rsquo;un bloc de donn\u00e9es d\u00e9sign\u00e9.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/weight_height1.jpg\" alt=\"Simple linear regression\" \/><figcaption>Simple linear regression<\/figcaption><\/figure>\n<p>Une fois l&rsquo;objet de r\u00e9gression instanci\u00e9 dans l&rsquo;espace de travail (par exemple sous la d\u00e9nomination <code>modele_psy<\/code>), le chercheur extrait les r\u00e9sidus ordinaires bruts au moyen des fonctions \u00e9quivalentes <code>residuals(modele_psy)<\/code> ou <code>resid(modele_psy)<\/code>. Simultan\u00e9ment, la commande <code>fitted(modele_psy)<\/code> fournit le vecteur exhaustif des pr\u00e9dictions th\u00e9oriques <em>\u0177<sub>i<\/sub><\/em> associ\u00e9es \u00e0 chaque ligne de la matrice exp\u00e9rimentale.<\/p>\n<p>Pour acc\u00e9der aux formes standardis\u00e9es et studentis\u00e9es sans recourir \u00e0 des formules manuelles, R dispose des fonctions de base <code>rstandard(modele_psy)<\/code> et <code>rstudent(modele_psy)<\/code>. L&rsquo;appel \u00e0 <code>rstudent()<\/code> applique automatiquement la formulation jackknife en calculant l&rsquo;erreur type r\u00e9siduelle externe par suppression it\u00e9rative d&rsquo;observation, fournissant des vecteurs pr\u00eats pour l&rsquo;inf\u00e9rence.<\/p>\n<h3>7.2 G\u00e9n\u00e9ration des quatre graphiques diagnostiques automatis\u00e9s<\/h3>\n<p>L&rsquo;une des fonctionnalit\u00e9s les plus puissantes du moteur statistique de R r\u00e9side dans sa routine graphique int\u00e9gr\u00e9e. L&rsquo;application de la fonction g\u00e9n\u00e9rique <code>plot()<\/code> sur un objet de mod\u00e9lisation lin\u00e9aire &ndash; formul\u00e9e sous la syntaxe <code>plot(modele_psy)<\/code> &ndash; engendre de mani\u00e8re s\u00e9quentielle un tableau d&rsquo;\u00e9valuation clinique compos\u00e9 de quatre cadrans diagnostiques fondamentaux :<\/p>\n<p>Le premier panneau (<em>Residuals vs Fitted<\/em>) trace la configuration des r\u00e9sidus ordinaires en relation avec les niveaux de pr\u00e9diction th\u00e9orique, assortie d&rsquo;un ajustement liss\u00e9 non param\u00e9trique local (courbe rouge de r\u00e9gression LOESS). Cette visualisation permet d&rsquo;identifier d&rsquo;\u00e9ventuelles non-lin\u00e9arit\u00e9s structurelles lorsque la courbe rouge d\u00e9vie d&rsquo;une rectitude horizontale.<\/p>\n<p>Le deuxi\u00e8me panneau (<em>Normal Q-Q<\/em>) confronte les quantiles empiriques des r\u00e9sidus standardis\u00e9s \u00e0 la distribution normale id\u00e9ale, mettant en relief les \u00e9carts de sym\u00e9trie et l&rsquo;\u00e9paisseur anormale des queues de distribution.<\/p>\n<p>Le troisi\u00e8me panneau (<em>Scale-Location<\/em>, ou <em>Spread-Location<\/em>) projette la racine carr\u00e9e des r\u00e9sidus standardis\u00e9s absolus (&radic;|<em>z<sub>i<\/sub><\/em>|) en fonction des valeurs ajust\u00e9es. Cette transformation stabilise la variance de l&rsquo;indicateur graphique, facilitant le rep\u00e9rage de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 par la surveillance de la pente du lissage LOESS.<\/p>\n<p>Le quatri\u00e8me panneau (<em>Residuals vs Leverage<\/em>), enfin, cartographie les r\u00e9sidus standardis\u00e9s contre les coefficients de levier individuel <em>h<sub>ii<\/sub><\/em>. Ce graphe superpose des courbes d&rsquo;isocontours en pointill\u00e9s d\u00e9limitant les zones critiques de la distance de Cook (typiquement pour <em>D = 0.5<\/em> et <em>D = 1.0<\/em>), permettant d&rsquo;isoler au premier regard les observations qui conjuguent un fort levier horizontal et un r\u00e9sidu vertical marqu\u00e9.<\/p>\n<h3>7.3 Visualisation avanc\u00e9e avec ggplot2 et packages sp\u00e9cialis\u00e9s<\/h3>\n<p>Pour r\u00e9pondre aux exigences \u00e9ditoriales modernes, l&rsquo;\u00e9cosyst\u00e8me graphique <code>ggplot2<\/code> offre une flexibilit\u00e9 de composition sup\u00e9rieure. En exploitant la fonction <code>geom_segment()<\/code>, l&rsquo;analyste peut g\u00e9n\u00e9rer un diagramme de dispersion enrichi visualisant directement la projection verticale des r\u00e9sidus sous forme de trac\u00e9s reliant chaque point empirique <em>(x<sub>i<\/sub>, y<sub>i<\/sub>)<\/em> \u00e0 sa projection pr\u00e9dite sur la droite <em>(x<sub>i<\/sub>, \u0177<sub>i<\/sub>)<\/em> :<\/p>\n<p>Cette mise en sc\u00e8ne graphique met en \u00e9vidence l&rsquo;\u00e9cart vertical individuel et facilite l&rsquo;explication visuelle des notions d&rsquo;erreurs quadratiques aupr\u00e8s d&rsquo;un public non sp\u00e9cialiste.<\/p>\n<p>Parall\u00e8lement, la biblioth\u00e8que sp\u00e9cialis\u00e9e <code>car<\/code> (<em>Companion to Applied Regression<\/em>), d\u00e9velopp\u00e9e par Fox et Weisberg, fournit une suite d&rsquo;outils analytiques d&rsquo;une pr\u00e9cision \u00e9prouv\u00e9e. La fonction <code>residualPlots(modele_psy)<\/code> calcule des tests d&rsquo;ad\u00e9quation non lin\u00e9aires pour chaque r\u00e9gresseur individuel, tandis que la fonction <code>influenceIndexPlot(modele_psy)<\/code> g\u00e9n\u00e8re une matrice synoptique comprenant les graphiques d&rsquo;alignement des r\u00e9sidus studentis\u00e9s, des distances de Cook, des indices de levier et des modifications structurelles DFBETAS.<\/p>\n<p>Enfin, l&rsquo;extension <code>broom<\/code> transforme ces flux analytiques en donn\u00e9es structur\u00e9es exploitables via la commande <code>augment(modele_psy)<\/code>. Cette fonction ajoute au jeu de donn\u00e9es initial l&rsquo;int\u00e9gralit\u00e9 des indicateurs r\u00e9siduels calcul\u00e9s (<code>.fitted<\/code>, <code>.resid<\/code>, <code>.std.resid<\/code>, <code>.cooksd<\/code>, <code>.hat<\/code>), facilitant les op\u00e9rations d&rsquo;exportation vers des tableaux de synth\u00e8se aux normes de publication de l&rsquo;American Psychological Association (<a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\/\">APA Style<\/a>).<\/p>\n<h2>8. Calcul et analyse des r\u00e9sidus avec Python (Statsmodels et Scikit-Learn)<\/h2>\n<h3>8.1 Mise en \u0153uvre compl\u00e8te avec Statsmodels<\/h3>\n<p>Dans l&rsquo;\u00e9cosyst\u00e8me open-source Python d\u00e9di\u00e9 \u00e0 la science des donn\u00e9es, la biblioth\u00e8que <code>statsmodels<\/code> constitue la r\u00e9f\u00e9rence acad\u00e9mique privil\u00e9gi\u00e9e pour conduire des mod\u00e9lisations fond\u00e9es sur la rigueur inf\u00e9rentielle. Contrairement aux biblioth\u00e8ques orient\u00e9es exclusivement vers l&rsquo;apprentissage automatique pr\u00e9dictif, <code>statsmodels<\/code> conserve l&rsquo;int\u00e9gralit\u00e9 des m\u00e9triques diagnostiques ch\u00e8res aux biostatisticiens et psychom\u00e9triciens.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/weight_height1.jpg\" alt=\"Simple linear regression\" \/><figcaption>Simple linear regression<\/figcaption><\/figure>\n<p>La configuration d&rsquo;une r\u00e9gression lin\u00e9aire s&rsquo;op\u00e8re via la classe <code>OLS<\/code> (<em>Ordinary Least Squares<\/em>), typiquement aliment\u00e9e par la m\u00e9thode matricielle ou via le sous-module d\u00e9claratif <code>statsmodels.formula.api<\/code>. Une fois la m\u00e9thode <code>.fit()<\/code> ex\u00e9cut\u00e9e sur l&rsquo;instance sp\u00e9cifi\u00e9e, le mod\u00e8le produit une structure d&rsquo;attributs d\u00e9di\u00e9e \u00e0 l&rsquo;extraction imm\u00e9diate des termes r\u00e9siduels. L&rsquo;attribut <code>results.resid<\/code> d\u00e9livre le vecteur des r\u00e9sidus ordinaires <em>e<sub>i<\/sub><\/em>, tandis que <code>results.fittedvalues<\/code> rassemble les projections th\u00e9oriques <em>\u0177<sub>i<\/sub><\/em>. La variance r\u00e9siduelle moyenne non biais\u00e9e (<em>MSE<\/em>) est accessible via l&rsquo;attribut <code>results.mse_resid<\/code>.<\/p>\n<p>Pour d\u00e9ployer des analyses d&rsquo;influence exhaustives, le module <code>statsmodels.stats.outliers_influence<\/code> offre la classe sp\u00e9cialis\u00e9e <code>OLSInfluence<\/code>. En instanciant cet objet \u00e0 partir des r\u00e9sultats de l&rsquo;ajustement, le chercheur extrait directement une structure de donn\u00e9es contenant les r\u00e9sidus studentis\u00e9s externes (via <code>resid_studentized_external<\/code>), les valeurs de levier diagonal <em>h<sub>ii<\/sub><\/em> (via <code>hat_matrix_diag<\/code>), et les distances de Cook (via la m\u00e9thode <code>cooks_distance<\/code>), simplifiant les op\u00e9rations d&rsquo;audit statistique.<\/p>\n<h3>8.2 Calcul manuel et vectoris\u00e9 avec Scikit-Learn et NumPy<\/h3>\n<p>Lorsque les volum\u00e9tries de donn\u00e9es s&rsquo;accroissent ou que l&rsquo;architecture d&rsquo;analyse s&rsquo;int\u00e8gre au sein d&rsquo;un pipeline de machine learning, la biblioth\u00e8que <code>scikit-learn<\/code> s&rsquo;impose par son efficience computationnelle. La classe <code>LinearRegression<\/code> issue du module <code>sklearn.linear_model<\/code> optimise l&rsquo;ajustement des param\u00e8tres sur les jeux d&rsquo;apprentissage.<\/p>\n<p>Toutefois, <code>scikit-learn<\/code> n&rsquo;encapsule pas directement les r\u00e9sidus sous forme d&rsquo;attributs d&rsquo;objets statistiques. Il incombe \u00e0 l&rsquo;analyste d&rsquo;ex\u00e9cuter l&rsquo;\u00e9valuation vectorielle par programmation matricielle avec <code>numpy<\/code>. Apr\u00e8s avoir extrait le vecteur des pr\u00e9dictions via l&rsquo;application de <code>model.predict(X)<\/code>, le calcul des r\u00e9sidus s&rsquo;op\u00e8re par une soustraction vectoris\u00e9e :<\/p>\n<p><em>residus = y &minus; model.predict(X)<\/em><\/p>\n<p>\u00c0 partir de ce vecteur r\u00e9siduel, l&rsquo;analyste d\u00e9rive facilement les m\u00e9triques synth\u00e9tiques fondamentales via le module <code>sklearn.metrics<\/code> : l&rsquo;erreur absolue moyenne (<em>MAE<\/em>, <em>Mean Absolute Error<\/em>), l&rsquo;erreur quadratique moyenne (<em>MSE<\/em>, <em>Mean Squared Error<\/em>), et sa racine carr\u00e9e (<em>RMSE<\/em>, <em>Root Mean Squared Error<\/em>).<\/p>\n<p>L&rsquo;exploration distributionnelle s&rsquo;effectue ensuite par le biais des biblioth\u00e8ques de rendu graphique <code>matplotlib<\/code> et <code>seaborn<\/code>. Un histogramme enrichi d&rsquo;une estimation de densit\u00e9 par noyau (KDE) appliqu\u00e9 aux r\u00e9sidus ordinaires, combin\u00e9 \u00e0 un nuage de dispersion <code>scatterplot(x=y_pred, y=residus)<\/code> assorti d&rsquo;une ligne de r\u00e9f\u00e9rence horizontale trac\u00e9e \u00e0 z\u00e9ro, permet de valider rapidement l&rsquo;absence de biais flagrant dans le syst\u00e8me pr\u00e9dictif.<\/p>\n<h3>8.3 Automatisation des tests d&rsquo;hypoth\u00e8ses sur les r\u00e9sidus<\/h3>\n<p>L&rsquo;automatisation des diagnostics de r\u00e9sidus au sein de scripts Python robustes s&rsquo;appuie sur la combinaison des fonctions analytiques de <code>scipy.stats<\/code> et de <code>statsmodels.stats.diagnostic<\/code>. Ce couplage permet d&rsquo;auditer les hypoth\u00e8ses de Gauss-Markov sans d\u00e9pendre d&rsquo;inspections visuelles manuelles.<\/p>\n<p>L&rsquo;\u00e9valuation formelle de la normalit\u00e9 de la distribution des erreurs s&rsquo;ex\u00e9cute de fa\u00e7on programmatique via le test de Shapiro-Wilk accessible dans la m\u00e9thode <code>scipy.stats.shapiro(residus)<\/code>, ou par le test d&rsquo;asym\u00e9trie et d&rsquo;aplatissement combin\u00e9 de Jarque-Bera au moyen de <code>statsmodels.stats.stattools.jarque_bera(residus)<\/code>. L&rsquo;obtention d&rsquo;une valeur <em>p<\/em> sup\u00e9rieure au seuil &alpha; conventionnel de 0,05 conclut au non-rejet de l&rsquo;hypoth\u00e8se de normalit\u00e9.<\/p>\n<p>La v\u00e9rification de l&rsquo;homosc\u00e9dasticit\u00e9 s&rsquo;automatise quant \u00e0 elle par l&rsquo;appel \u00e0 la fonction <code>het_breuschpagan(residus, X_avec_constante)<\/code>, laquelle retourne la statistique du multiplicateur de Lagrange ainsi que sa <em>p<\/em>-valeur associ\u00e9e. Enfin, le contr\u00f4le de l&rsquo;ind\u00e9pendance mutuelle des r\u00e9sidus s\u00e9quentiels est confi\u00e9 \u00e0 la fonction <code>stat_durbin_watson(residus)<\/code> issue de <code>statsmodels.stats.stattools<\/code>, qui confirme l&rsquo;absence d&rsquo;autocorr\u00e9lation lin\u00e9aire d&rsquo;ordre un lorsque la statistique converge autour de la valeur de r\u00e9f\u00e9rence 2,0.<\/p>\n<h2>9. Calcul et interpr\u00e9tation des r\u00e9sidus dans SPSS et JASP<\/h2>\n<h3>9.1 Param\u00e9trage de la r\u00e9gression lin\u00e9aire sous IBM SPSS<\/h3>\n<p>Le logiciel commercial <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/www.ibm.com\/fr-fr\/products\/spss-statistics\">IBM SPSS Statistics<\/a> demeure une plateforme tr\u00e8s r\u00e9pandue pour le traitement des donn\u00e9es en psychologie clinique, en sciences de l&rsquo;\u00e9ducation et en sociologie quantitative. L&rsquo;analyse de r\u00e9gression s&rsquo;initialise via l&rsquo;arborescence des menus : <em>Analyser &gt; R\u00e9gression &gt; Lin\u00e9aire&#8230;<\/em> Apr\u00e8s avoir positionn\u00e9 la variable crit\u00e8re et le ou les r\u00e9gresseurs dans leurs champs respectifs, l&rsquo;acc\u00e8s aux composantes r\u00e9siduelles requiert la configuration fine de plusieurs fen\u00eatres de dialogue associ\u00e9es.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/statology_residuals3.jpg\" alt=\"\" \/><\/figure>\n<p>Dans le sous-menu intitul\u00e9 <em>Enregistrer&#8230;<\/em> (<em>Save<\/em>), le chercheur coche les indicateurs qu&rsquo;il souhaite injecter directement dans sa matrice de donn\u00e9es active. Il est recommand\u00e9 de s\u00e9lectionner :<\/p>\n<ul>\n<li>Dans la section <em>R\u00e9sidus<\/em> : les r\u00e9sidus non standardis\u00e9s (<em>Unstandardized<\/em>), les r\u00e9sidus standardis\u00e9s (<em>Standardized<\/em>) et les r\u00e9sidus studentis\u00e9s supprim\u00e9s (<em>Studentized deleted<\/em>).<\/li>\n<li>Dans la section <em>Distances<\/em> : les distances de Cook (<em>Cook&rsquo;s<\/em>) ainsi que les valeurs de levier (<em>Leverage values<\/em>).<\/li>\n<li>Dans la section <em>Statistiques d&rsquo;influence<\/em> : les indices DfBeta(s) et DfFit.<\/li>\n<\/ul>\n<p>Parall\u00e8lement, la bo\u00eete de dialogue <em>Trac\u00e9s&#8230;<\/em> (<em>Plots<\/em>) permet de configurer la g\u00e9n\u00e9ration automatis\u00e9e des graphiques de contr\u00f4le. L&rsquo;introduction de la variable <code>*ZPRED<\/code> (valeurs pr\u00e9dites standardis\u00e9es) sur l&rsquo;axe des <em>X<\/em> et de la variable <code>*ZRESID<\/code> (r\u00e9sidus standardis\u00e9s) sur l&rsquo;axe des <em>Y<\/em> produit le diagramme de dispersion indispensable \u00e0 l&rsquo;\u00e9valuation de l&rsquo;homosc\u00e9dasticit\u00e9. Le cochage syst\u00e9matique de la case <em>Trac\u00e9 du graphique de probabilit\u00e9 normale<\/em> (<em>Normal probability plot<\/em>) commande quant \u00e0 lui la sortie du diagramme Q-Q r\u00e9siduel.<\/p>\n<h3>9.2 Proc\u00e9dure d&rsquo;analyse intuitive dans JASP<\/h3>\n<p>Proposant une alternative libre, moderne et ergonomique aux logiciels d&rsquo;analyse traditionnels, la plateforme logicielle open-source <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/jasp-stats.org\/\">JASP<\/a>, d\u00e9velopp\u00e9e \u00e0 l&rsquo;Universit\u00e9 d&rsquo;Amsterdam, simplifie l&rsquo;ex\u00e9cution des diagnostics d&rsquo;ad\u00e9quation statistique. L&rsquo;analyse s&rsquo;amorce en s\u00e9lectionnant l&rsquo;onglet <em>Regression<\/em>, suivi de la commande <em>Linear Regression<\/em>.<\/p>\n<p>L&rsquo;affectation des variables s&rsquo;effectue par glisser-d\u00e9poser interactif au sein des panneaux <em>Dependent Variable<\/em> et <em>Covariates<\/em>. La singularit\u00e9 de JASP r\u00e9side dans la mise \u00e0 jour dynamique et imm\u00e9diate des r\u00e9sultats au sein de la fen\u00eatre de visualisation lat\u00e9rale d\u00e8s qu&rsquo;une option est coch\u00e9e, favorisant un flux de travail transparent.<\/p>\n<p>En d\u00e9ployant l&rsquo;onglet d\u00e9roulant intitul\u00e9 <em>Diagnostics<\/em>, l&rsquo;utilisateur active un \u00e9ventail complet d&rsquo;options d&rsquo;\u00e9valuation :<\/p>\n<ul>\n<li>La s\u00e9lection de <em>Residuals vs. fitted<\/em> g\u00e9n\u00e8re instantan\u00e9ment le graphique de dispersion des r\u00e9sidus centr\u00e9 sur l&rsquo;axe nul.<\/li>\n<li>L&rsquo;option <em>Q-Q plot residuals<\/em> trace le graphique quantile-quantile des r\u00e9sidus, assorti de bandes de confiance empiriques.<\/li>\n<li>L&rsquo;activation de l&rsquo;option <em>Casewise diagnostics<\/em> g\u00e9n\u00e8re un tableau identifiant automatiquement les observations pr\u00e9sentant des r\u00e9sidus standardis\u00e9s exc\u00e9dant deux ou trois \u00e9carts-types.<\/li>\n<li>L&rsquo;affichage des coefficients d&rsquo;asym\u00e9trie (<em>skewness<\/em>) et d&rsquo;aplatissement (<em>kurtosis<\/em>) des r\u00e9sidus permet d&rsquo;\u00e9valuer la forme distributionnelle sans recourir \u00e0 des calculs manuels.<\/li>\n<\/ul>\n<h3>9.3 Lecture critique des tableaux de sortie logiciels<\/h3>\n<p>L&rsquo;interpr\u00e9tation des r\u00e9sultats g\u00e9n\u00e9r\u00e9s par ces environnements logiciels exige une lecture rigoureuse des tableaux de synth\u00e8se standards, \u00e9vitant de se focaliser pr\u00e9matur\u00e9ment sur les seules valeurs <em>p<\/em> des coefficients de r\u00e9gression.<\/p>\n<p>Le premier tableau \u00e0 inspecter attentivement est le <em>R\u00e9capitulatif des mod\u00e8les<\/em> (<em>Model Summary<\/em>). Ce tableau fournit le <em>R<\/em>, le <em>R-deux<\/em> (coefficient de d\u00e9termination), ainsi que l&rsquo;<em>Erreur standard de l&rsquo;estimation<\/em>. Cette derni\u00e8re correspond formellement \u00e0 <em>s<sub>e<\/sub> = &radic;MSE<\/em> ; elle quantifie la dispersion moyenne absolue des erreurs r\u00e9siduelles exprim\u00e9e dans l&rsquo;unit\u00e9 de mesure de la variable r\u00e9ponse. Si un mod\u00e8le affiche un <em>R<sup>2<\/sup><\/em> \u00e9lev\u00e9 mais une erreur standard disproportionn\u00e9e au regard de l&rsquo;\u00e9chelle d&rsquo;\u00e9valuation clinique, sa valeur pr\u00e9dictive unitaire demeure modeste.<\/p>\n<p>Le tableau de l&rsquo;<em>Analyse de variance<\/em> (<em>ANOVA<\/em>) explicite ensuite la d\u00e9composition formelle des sommes de carr\u00e9s (<em>Regression<\/em>, <em>Residual<\/em>, <em>Total<\/em>). Le rapport entre le carr\u00e9 moyen de la r\u00e9gression (<em>MSR = SSR \/ k<\/em>) et le carr\u00e9 moyen de l&rsquo;erreur (<em>MSE = SSE \/ (n &minus; k &minus; 1)<\/em>) produit la statistique du test <em>F<\/em> de Fisher. Un test <em>F<\/em> statistiquement significatif indique que le mod\u00e8le lin\u00e9aire extrait une variance substantiellement sup\u00e9rieure \u00e0 celle g\u00e9n\u00e9r\u00e9e par un processus r\u00e9siduel al\u00e9atoire.<\/p>\n<p>Enfin, le tableau du <em>Diagnostic des cas extr\u00eames<\/em> (<em>Casewise Diagnostics<\/em>) dresse un inventaire des observations atypiques. L&rsquo;analyste doit y v\u00e9rifier si certains cas isol\u00e9s pr\u00e9sentent des r\u00e9sidus standardis\u00e9s extr\u00eames (|<em>z<\/em>| &gt; 3,0). La pr\u00e9sence conjointe de ces r\u00e9sidus \u00e9lev\u00e9s avec des distances de Cook critiques &ndash; report\u00e9es dans les colonnes ajout\u00e9es au fichier de donn\u00e9es &ndash; impose un examen d\u00e9taill\u00e9 des donn\u00e9es de ces participants avant toute conclusion d\u00e9finitive.<\/p>\n<h2>10. Les r\u00e9sidus dans la r\u00e9gression multiple et les mod\u00e8les non lin\u00e9aires<\/h2>\n<h3>10.1 Calcul matriciel des r\u00e9sidus en r\u00e9gression multiple<\/h3>\n<p>L&rsquo;extension du mod\u00e8le de r\u00e9gression \u00e0 un ensemble multidimensionnel regroupant <em>k<\/em> pr\u00e9dicteurs explicatifs requiert l&rsquo;utilisation de l&rsquo;<a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Alg%C3%A8bre_lin%C3%A9aire\">alg\u00e8bre matricielle<\/a>. L&rsquo;\u00e9quation globale du syst\u00e8me prend la forme condens\u00e9e : <em>Y = X&beta; + e<\/em>, o\u00f9 <em>Y<\/em> est un vecteur colonne d&rsquo;observations de dimension <em>n &times; 1<\/em>, <em>X<\/em> repr\u00e9sente la matrice du plan exp\u00e9rimental de dimension <em>n &times; (k + 1)<\/em> &ndash; int\u00e9grant une colonne unitaire initiale assign\u00e9e \u00e0 la constante &ndash;, <em>&beta;<\/em> symbolise le vecteur des param\u00e8tres inconnus de dimension <em>(k + 1) &times; 1<\/em>, et <em>e<\/em> d\u00e9signe le vecteur des r\u00e9sidus de dimension <em>n &times; 1<\/em>.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/statology_residuals4.jpg\" alt=\"\" \/><\/figure>\n<p>L&rsquo;estimateur des moindres carr\u00e9s ordinaires s&rsquo;obtient par la r\u00e9solution de l&rsquo;\u00e9quation matricielle classique : <em>&beta;&#770; = (X<sup>T<\/sup>X)<sup>&minus;1<\/sup>X<sup>T<\/sup>Y<\/em>. Le vecteur des valeurs th\u00e9oriques pr\u00e9dites d\u00e9coule de l&rsquo;application : <em>\u0176 = X&beta;&#770; = X(X<sup>T<\/sup>X)<sup>&minus;1<\/sup>X<sup>T<\/sup>Y<\/em>. L&rsquo;op\u00e9rateur matriciel central <em>H = X(X<sup>T<\/sup>X)<sup>&minus;1<\/sup>X<sup>T<\/sup><\/em> est formellement d\u00e9sign\u00e9 sous le terme de matrice chapeau (ou <em>Hat Matrix<\/em>), pr\u00e9cis\u00e9ment parce qu&rsquo;il projette le vecteur des valeurs observ\u00e9es <em>Y<\/em> vers son approximation vectorielle ajust\u00e9e <em>\u0176<\/em> : <em>\u0176 = HY<\/em>.<\/p>\n<p>Le calcul du vecteur r\u00e9siduel s&rsquo;exprime alors \u00e0 travers une relation matricielle directe :<\/p>\n<p><em>e = Y &minus; \u0176 = Y &minus; HY = (I &minus; H)Y<\/em><\/p>\n<p>Dans cette formulation, <em>I<\/em> caract\u00e9rise la matrice identit\u00e9 de dimension <em>n &times; n<\/em>. La matrice <em>(I &minus; H)<\/em> poss\u00e8de des propri\u00e9t\u00e9s g\u00e9om\u00e9triques fondamentales : elle est \u00e0 la fois sym\u00e9trique (<em>(I &minus; H)<sup>T<\/sup> = (I &minus; H)<\/em>) et idempotente (<em>(I &minus; H)(I &minus; H) = (I &minus; H)<\/em>). Cela traduit le fait qu&rsquo;elle op\u00e8re une projection orthogonale du vecteur observ\u00e9 <em>Y<\/em> sur le sous-espace vectoriel r\u00e9siduel d&rsquo;ordre <em>n &minus; k &minus; 1<\/em>, math\u00e9matiquement orthogonal \u00e0 l&rsquo;espace vectoriel engendr\u00e9 par les colonnes de pr\u00e9dicteurs de <em>X<\/em>.<\/p>\n<h3>10.2 R\u00e9sidus partiels et graphiques de r\u00e9gression r\u00e9siduelle partielle<\/h3>\n<p>Lorsque plusieurs pr\u00e9dicteurs corr\u00e9l\u00e9s interviennent conjointement dans l&rsquo;explication d&rsquo;un comportement, le simple examen du lien bivari\u00e9 brut entre un pr\u00e9dicteur <em>X<sub>j<\/sub><\/em> et la r\u00e9ponse <em>Y<\/em> peut v\u00e9hiculer des conclusions trompeuses en raison du ph\u00e9nom\u00e8ne de confusion statistique ou de suppression. Pour appr\u00e9hender l&rsquo;impact net associ\u00e9 \u00e0 une covariable particuli\u00e8re, la construction de r\u00e9sidus de r\u00e9gression partielle s&rsquo;av\u00e8re indispensable.<\/p>\n<p>Cette approche implique la g\u00e9n\u00e9ration conjointe de deux vecteurs r\u00e9siduels ind\u00e9pendants :<\/p>\n<ul>\n<li>On mod\u00e9lise en premier lieu <em>Y<\/em> par l&rsquo;ensemble des pr\u00e9dicteurs en excluant d\u00e9lib\u00e9r\u00e9ment <em>X<sub>j<\/sub><\/em>, et l&rsquo;on extrait le vecteur des r\u00e9sidus r\u00e9sultant <em>e<sub>Y | X(-j)<\/sub><\/em>. Ce vecteur repr\u00e9sente la variance de la r\u00e9ponse qui ne peut pas \u00eatre expliqu\u00e9e par les autres covariables.<\/li>\n<li>On mod\u00e9lise en second lieu le pr\u00e9dicteur <em>X<sub>j<\/sub><\/em> lui-m\u00eame en fonction de l&rsquo;ensemble de ces m\u00eames pr\u00e9dicteurs restants, et l&rsquo;on isole le second vecteur r\u00e9siduel <em>e<sub>Xj | X(-j)<\/sub><\/em>. Ce vecteur caract\u00e9rise la part de variabilit\u00e9 unique du pr\u00e9dicteur <em>X<sub>j<\/sub><\/em>, ind\u00e9pendante de ses relations avec les autres variables explicatives.<\/li>\n<\/ul>\n<p>Le trac\u00e9 du diagramme de dispersion reliant ces deux grandeurs r\u00e9siduelles &ndash; <em>e<sub>Y | X(-j)<\/sub><\/em> en ordonn\u00e9e et <em>e<sub>Xj | X(-j)<\/sub><\/em> en abscisse &ndash; constitue le graphique de r\u00e9gression partielle (ou <em>Added-Variable Plot<\/em>). La pente de la droite des moindres carr\u00e9s traversant ce graphique \u00e9quivaut exactement au coefficient de r\u00e9gression partielle <em>b<sub>j<\/sub><\/em> calcul\u00e9 au sein du mod\u00e8le multiple complet. Ce dispositif offre ainsi une fen\u00eatre visuelle directe pour inspecter d&rsquo;\u00e9ventuelles non-lin\u00e9arit\u00e9s localis\u00e9es sur un pr\u00e9dicteur sp\u00e9cifique, une fois contr\u00f4l\u00e9e l&rsquo;influence des autres variables de contr\u00f4le.<\/p>\n<h3>10.3 R\u00e9sidus dans les mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s (GLM) et logistiques<\/h3>\n<p>L&rsquo;\u00e9difice des r\u00e9sidus ordinaires <em>e<sub>i<\/sub> = y<sub>i<\/sub> &minus; \u0177<sub>i<\/sub><\/em> atteint ses limites m\u00e9thodologiques d\u00e8s lors que la variable r\u00e9ponse s&rsquo;\u00e9carte du cadre des distributions continues et sym\u00e9triques. En psychom\u00e9trie exp\u00e9rimentale, les chercheurs recueillent r\u00e9guli\u00e8rement des variables dichotomiques (ex. : r\u00e9ussite ou \u00e9chec \u00e0 un item de test, pr\u00e9sence ou absence d&rsquo;un trouble psychopathologique) ou des donn\u00e9es de comptage (ex. : nombre d&rsquo;erreurs commises lors d&rsquo;une t\u00e2che d&rsquo;attention soutenue). Ces configurations rel\u00e8vent des <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Mod%C3%A8le_lin%C3%A9aire_g%C3%A9n%C3%A9ralis%C3%A9\">mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s<\/a> (GLM).<\/p>\n<p>Dans le cas d&rsquo;une r\u00e9gression logistique o\u00f9 la r\u00e9ponse <em>y<sub>i<\/sub> &isin; {0, 1}<\/em> et o\u00f9 la pr\u00e9diction correspond \u00e0 une probabilit\u00e9 conditionnelle <em>p&#770;<sub>i<\/sub><\/em>, le r\u00e9sidu brut ne peut adopter que deux valeurs possibles pour une pr\u00e9diction donn\u00e9e : <em>1 &minus; p&#770;<sub>i<\/sub><\/em> ou <em>0 &minus; p&#770;<sub>i<\/sub><\/em>. Ce caract\u00e8re discret invalide les postulats traditionnels d&rsquo;homosc\u00e9dasticit\u00e9 et de normalit\u00e9. Deux formulations r\u00e9siduelles adapt\u00e9es s&rsquo;imposent alors :<\/p>\n<p>Le <strong>r\u00e9sidu de Pearson<\/strong> pour GLM, qui divise l&rsquo;\u00e9cart brut par l&rsquo;\u00e9cart-type th\u00e9orique d\u00e9coulant de la loi de distribution sous-jacente (loi binomiale dans ce contexte) :<\/p>\n<p><em>r<sub>P, i<\/sub> = (y<sub>i<\/sub> &minus; p&#770;<sub>i<\/sub>) \/ &radic;[p&#770;<sub>i<\/sub>(1 &minus; p&#770;<sub>i<\/sub>)]<\/em><\/p>\n<p>Le <strong>r\u00e9sidu de d\u00e9viance<\/strong> (<em>Deviance residual<\/em>, not\u00e9 <em>d<sub>i<\/sub><\/em>), qui d\u00e9compose la d\u00e9viance totale du mod\u00e8le (\u00e9quivalent fonctionnel de la somme des carr\u00e9s r\u00e9siduels dans le paradigme de la r\u00e9gression logistique) en contributions individuelles sign\u00e9es :<\/p>\n<p><em>d<sub>i<\/sub> = signe(y<sub>i<\/sub> &minus; p&#770;<sub>i<\/sub>) &times; &radic;{2 [y<sub>i<\/sub> ln(y<sub>i<\/sub> \/ p&#770;<sub>i<\/sub>) + (1 &minus; y<sub>i<\/sub>) ln((1 &minus; y<sub>i<\/sub>) \/ (1 &minus; p&#770;<sub>i<\/sub>))]}<\/em><\/p>\n<p>Pour les donn\u00e9es de comptage psychom\u00e9triques (mod\u00e9lis\u00e9es via une r\u00e9gression de Poisson ou binomiale n\u00e9gative), ces r\u00e9sidus de d\u00e9viance jouent le m\u00eame r\u00f4le que les r\u00e9sidus ordinaires en r\u00e9gression classique : leur convergence vers une distribution approximativement gaussienne autorise l&rsquo;identification des profils atypiques et la d\u00e9tection des ph\u00e9nom\u00e8nes de surdispersion.<\/p>\n<h2>11. Strat\u00e9gies m\u00e9thodologiques face \u00e0 des anomalies de r\u00e9sidus<\/h2>\n<h3>11.1 Transformations de variables pour corriger la non-lin\u00e9arit\u00e9 et l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9<\/h3>\n<p>La d\u00e9couverte d&rsquo;anomalies structurelles lors du diagnostic des r\u00e9sidus &ndash; telles qu&rsquo;une courbure non lin\u00e9aire ou un \u00e9largissement h\u00e9t\u00e9rosc\u00e9dastique en entonnoir &ndash; ne commande pas n\u00e9cessairement l&rsquo;abandon imm\u00e9diat du paradigme lin\u00e9aire. La premi\u00e8re ligne d&rsquo;intervention m\u00e9thodologique consiste \u00e0 appliquer une transformation math\u00e9matique appropri\u00e9e sur la variable d\u00e9pendante <em>Y<\/em>, sur les pr\u00e9dicteurs <em>X<\/em>, ou sur les deux composantes simultan\u00e9ment.<\/p>\n<p>L&rsquo;application d&rsquo;une fonction concave &ndash; telle que le logarithme n\u00e9p\u00e9rien <em>ln(Y)<\/em> ou la racine carr\u00e9e <em>&radic;Y<\/em> &ndash; est efficace pour corriger une asym\u00e9trie r\u00e9siduelle positive persistante combin\u00e9e \u00e0 une h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 o\u00f9 la variance augmente proportionnellement \u00e0 la pr\u00e9diction. Dans le contexte de l&rsquo;\u00e9valuation de latences cognitives ou de seuils perceptifs, la transformation inverse (<em>1 \/ Y<\/em>) permet souvent de stabiliser les distributions d&rsquo;erreurs.<\/p>\n<p>Pour d\u00e9passer les t\u00e2tonnements empiriques informels, les statisticiens s&rsquo;appuient sur la famille param\u00e9trique des <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Box%E2%80%93Cox_transformation\">transformations de Box-Cox<\/a>. Cette m\u00e9thode d&rsquo;optimisation par le maximum de vraisemblance estime un exposant de puissance &lambda; adapt\u00e9 aux donn\u00e9es :<\/p>\n<p><em>y<sup>(&lambda;)<\/sup> = (y<sup>&lambda;<\/sup> &minus; 1) \/ &lambda;<\/em> si &lambda; &ne; 0, et <em>y<sup>(&lambda;)<\/sup> = ln(y)<\/em> si &lambda; = 0<\/p>\n<p>Si la transformation de Box-Cox r\u00e9gularise souvent la forme des r\u00e9sidus et restaure l&rsquo;homog\u00e9n\u00e9it\u00e9 de leur variance, elle impose une contrepartie interpr\u00e9tative : les coefficients de r\u00e9gression estim\u00e9s ne s&rsquo;interpr\u00e8tent plus sur l&rsquo;\u00e9chelle m\u00e9trique d&rsquo;origine, n\u00e9cessitant des conversions adapt\u00e9es pour communiquer les conclusions cliniques aux non-statisticiens.<\/p>\n<h3>11.2 R\u00e9gression robuste et moindres carr\u00e9s pond\u00e9r\u00e9s (WLS)<\/h3>\n<p>Lorsque la structure de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 est explicitement connue &ndash; par exemple lorsque la variance de l&rsquo;erreur r\u00e9siduelle est proportionnelle \u00e0 la valeur d&rsquo;une variable ind\u00e9pendante connue <em>x<sub>i<\/sub><\/em> ou \u00e0 la taille des groupes d&rsquo;agr\u00e9gation empiriques &ndash;, l&rsquo;approche d&rsquo;estimation par les Moindres Carr\u00e9s Pond\u00e9r\u00e9s (<em>Weighted Least Squares<\/em> ou WLS) offre une solution performante. La m\u00e9thode consiste \u00e0 attribuer \u00e0 chaque observation un coefficient de pond\u00e9ration inversement proportionnel \u00e0 sa variance r\u00e9siduelle propre : <em>w<sub>i<\/sub> = 1 \/ &sigma;<sub>i<\/sub><sup>2<\/sup><\/em>. Ce m\u00e9canisme temp\u00e8re l&rsquo;impact des observations tr\u00e8s dispers\u00e9es au profit de celles issues de segments stables.<\/p>\n<p>Dans d&rsquo;autres configurations de recherche, les r\u00e9sidus r\u00e9v\u00e8lent des violations de la normalit\u00e9 caus\u00e9es par des observations extr\u00eames et atypiques r\u00e9currentes, sans motif m\u00e9thodologique justifiant leur exclusion unilat\u00e9rale. Dans ce cadre, le recours aux r\u00e9gressions robustes &ndash; et singuli\u00e8rement aux estimateurs de type M (dont l&rsquo;arch\u00e9type est l&rsquo;<a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Estimateur_de_Huber\">estimateur de Huber<\/a>) &ndash; s&rsquo;av\u00e8re particuli\u00e8rement indiqu\u00e9.<\/p>\n<p>L&rsquo;estimation de Huber modifie la fonction de minimisation des \u00e9carts : elle conserve un comportement quadratique classique (type MCO) pour les r\u00e9sidus de faible amplitude situ\u00e9s \u00e0 l&rsquo;int\u00e9rieur d&rsquo;un seuil critique (g\u00e9n\u00e9ralement fix\u00e9 \u00e0 <em>|e<sub>i<\/sub>| &le; 1,345 s<sub>e<\/sub><\/em>), mais commute vers une p\u00e9nalisation lin\u00e9aire absolue pour les r\u00e9sidus plus \u00e9loign\u00e9s. Cette transformation limite le pouvoir d&rsquo;attraction des r\u00e9sidus \u00e9lev\u00e9s sur le calcul de la pente :<\/p>\n<p>Pour des situations de contamination plus s\u00e9v\u00e8re, la fonction de perte bisquarr\u00e9e (<em>biweight<\/em>) de Tukey neutralise compl\u00e8tement l&rsquo;influence des points dont le r\u00e9sidu d\u00e9passe un seuil de dispersion fix\u00e9, garantissant une estimation robuste des coefficients structuraux.<\/p>\n<h3>11.3 R\u00e9vision de la sp\u00e9cification th\u00e9orique du mod\u00e8le psychologique<\/h3>\n<p>L&rsquo;observation d&rsquo;une anomalie r\u00e9siduelle majeure constitue souvent le sympt\u00f4me direct d&rsquo;une inad\u00e9quation th\u00e9orique sous-jacente, soulignant la n\u00e9cessit\u00e9 d&rsquo;une r\u00e9vision du mod\u00e8le conceptuel. L&rsquo;\u00e9mergence d&rsquo;une structure non al\u00e9atoire dans les r\u00e9sidus signale fr\u00e9quemment l&rsquo;omission d&rsquo;un pr\u00e9dicteur pertinent ou d&rsquo;un terme d&rsquo;interaction mod\u00e9ratrice.<\/p>\n<p>Si la trajectoire r\u00e9siduelle prend une allure concave ou convexe selon les niveaux d&rsquo;un pr\u00e9dicteur continu, l&rsquo;introduction formelle d&rsquo;un terme polynomial d&rsquo;ordre sup\u00e9rieur (<em>X<sup>2<\/sup><\/em>) dans l&rsquo;\u00e9quation de r\u00e9gression peut redresser le mod\u00e8le :<\/p>\n<p><em>Y = b<sub>0<\/sub> + b<sub>1<\/sub>X + b<sub>2<\/sub>X<sup>2<\/sup> + e<\/em><\/p>\n<p>De m\u00eame, si la variance r\u00e9siduelle oscille en fonction des niveaux d&rsquo;un facteur qualitatif cat\u00e9goriel (comme le groupe d&rsquo;appartenance ou le genre), cette h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 r\u00e9siduelle traduit souvent l&rsquo;omission d&rsquo;un effet d&rsquo;interaction multiplicatif (<em>X &times; Groupe<\/em>). L&rsquo;int\u00e9gration de ce terme d&rsquo;interaction permet de mod\u00e9liser une pente distincte pour chaque sous-groupe, r\u00e9sorbant ainsi la dispersion apparente des r\u00e9sidus.<\/p>\n<p>Enfin, lorsque les diagnostics mettent en \u00e9vidence une violation flagrante du postulat d&rsquo;ind\u00e9pendance des erreurs due \u00e0 un regroupement naturel des participants (classes, services hospitaliers, fratries), le paradigme MCO classique s&rsquo;av\u00e8re inadapt\u00e9. Le chercheur doit alors mobiliser les mod\u00e8les lin\u00e9aires mixtes (<em>Linear Mixed Models<\/em>, LMM) ou la mod\u00e9lisation hi\u00e9rarchique \u00e0 effets al\u00e9atoires. Ces architectures avanc\u00e9es d\u00e9composent la variance r\u00e9siduelle en deux niveaux : une variance r\u00e9siduelle inter-groupes (entre clusters) et une variance r\u00e9siduelle intra-groupe (individuelle), restaurant la rigueur de l&rsquo;inf\u00e9rence statistique sur les donn\u00e9es imbriqu\u00e9es.<\/p>\n<h2>12. \u00c9tude de cas int\u00e9grative en psychom\u00e9trie : \u00c9puisement professionnel et empathie<\/h2>\n<h3>12.1 Cadre th\u00e9orique, recueil des donn\u00e9es et mod\u00e9lisation initiale<\/h3>\n<p>Afin de consolider les principes th\u00e9oriques et op\u00e9ratoires d\u00e9velopp\u00e9s tout au long de cet ouvrage, examinons une \u00e9tude empirique men\u00e9e aupr\u00e8s d&rsquo;une cohorte de <em>n = 120<\/em> personnels soignants hospitaliers en exercice dans des unit\u00e9s de soins critiques. L&rsquo;objectif th\u00e9orique de cette investigation en psychologie de la sant\u00e9 au travail consiste \u00e0 tester l&rsquo;hypoth\u00e8se de l&rsquo;\u00e9rosion de l&#8217;empathie face \u00e0 l&rsquo;usure professionnelle continue. Le score d&rsquo;\u00e9puisement \u00e9motionnel (\u00e9valu\u00e9 via la sous-\u00e9chelle d\u00e9di\u00e9e du <em>Maslach Burnout Inventory<\/em>, not\u00e9 <em>Burnout<\/em>, variant de 0 \u00e0 54) fait office de variable pr\u00e9dictive (<em>X<\/em>), tandis que l&rsquo;\u00e9valuation du niveau d&#8217;empathie clinique (mesur\u00e9e via la <em>Jefferson Scale of Empathy<\/em>, not\u00e9e <em>Empathie<\/em>, s&rsquo;\u00e9chelonnant sur un continuum de 20 \u00e0 140 points) constitue la variable r\u00e9ponse crit\u00e8re (<em>Y<\/em>).<\/p>\n<p>Le chercheur ajuste initialement un mod\u00e8le de r\u00e9gression lin\u00e9aire simple classique par les moindres carr\u00e9s ordinaires. Les r\u00e9sultats param\u00e9triques initiaux se formulent selon l&rsquo;\u00e9quation de r\u00e9gression suivante :<\/p>\n<p><em>Empathie&#770; = 118,45 &minus; 0,62 &times; Burnout<\/em><\/p>\n<p>Le test de conformit\u00e9 globale affiche un <em>F(1, 118) = 28,34<\/em>, associ\u00e9 \u00e0 une significativit\u00e9 statistique marqu\u00e9e (<em>p<\/em> &lt; .001), et le coefficient de d\u00e9termination s&rsquo;\u00e9tablit \u00e0 <em>R<sup>2<\/sup> = 0,194<\/em>. L&rsquo;interpr\u00e9tation h\u00e2tive de cette mod\u00e9lisation pr\u00e9liminaire sugg\u00e8re qu&rsquo;une \u00e9l\u00e9vation d&rsquo;une unit\u00e9 sur l&rsquo;\u00e9chelle du burnout provoque une r\u00e9duction attendue moyenne de 0,62 point d&#8217;empathie clinique chez ces soignants, ce pr\u00e9dicteur expliquant pr\u00e8s de 19,4 % de la variance de la r\u00e9ponse. N\u00e9anmoins, valider ce mod\u00e8le sans conduire un audit approfondi des r\u00e9sidus reviendrait \u00e0 ignorer de potentielles violations m\u00e9thodologiques fragilisant ces conclusions.<\/p>\n<h3>12.2 Audit complet des r\u00e9sidus et identification des violations<\/h3>\n<p>L&rsquo;investigateur engage alors l&rsquo;audit m\u00e9thodique des r\u00e9sidus d&rsquo;estimation g\u00e9n\u00e9r\u00e9s par le mod\u00e8le initial. L&rsquo;extraction syst\u00e9matique du vecteur r\u00e9siduel <em>e<sub>i<\/sub><\/em> et l&rsquo;\u00e9tablissement des valeurs pr\u00e9dites <em>\u0177<sub>i<\/sub><\/em> permettent de construire la batterie de contr\u00f4les diagnostiques fondamentaux.<\/p>\n<p>L&rsquo;\u00e9valuation commence par l&rsquo;examen du diagramme de dispersion des r\u00e9sidus standardis\u00e9s en relation avec les valeurs ajust\u00e9es (<em>Residuals vs Fitted<\/em>). Loin d&rsquo;afficher une bande horizontale d&rsquo;\u00e9paisseur homog\u00e8ne, le nuage de points dessine une forme en m\u00e9gaphone : pour les valeurs pr\u00e9dites \u00e9lev\u00e9es d&#8217;empathie (associ\u00e9es \u00e0 des scores d&rsquo;\u00e9puisement faibles), les r\u00e9sidus standardis\u00e9s sont fortement resserr\u00e9s autour de l&rsquo;axe z\u00e9ro (intervalle compris entre &minus;1,2 et +1,1). En revanche, \u00e0 mesure que la valeur pr\u00e9dite d\u00e9cline (associ\u00e9e \u00e0 un \u00e9puisement professionnel \u00e9lev\u00e9), la dispersion verticale des r\u00e9sidus s&rsquo;\u00e9largit substantiellement, oscillant entre &minus;3,4 et +2,8. Cette configuration h\u00e9t\u00e9rosc\u00e9dastique est confirm\u00e9e par le test de Breusch-Pagan, qui affiche une significativit\u00e9 marqu\u00e9e : <em>BP = 14,82, df = 1, p &lt; .001<\/em>.<\/p>\n<p>Simultan\u00e9ment, l&rsquo;analyse du diagramme quantile-quantile (<em>Normal Q-Q Plot<\/em>) montre un d\u00e9tachement manifeste des r\u00e9sidus dans les quantiles inf\u00e9rieurs de la distribution, caract\u00e9ris\u00e9 par un fl\u00e9chissement sous la diagonale th\u00e9orique \u00e0 45 degr\u00e9s. Le test de Shapiro-Wilk confirme cette violation de la normalit\u00e9 distributionnelle (<em>W = 0,948, p = .0002<\/em>).<\/p>\n<p>Enfin, le calcul des distances de Cook met en lumi\u00e8re deux observations influentes majeures. En appliquant le seuil critique adapt\u00e9 \u00e0 la cohorte <em>D<sub>crit<\/sub> = 4 \/ n = 4 \/ 120 &asymp; 0,033<\/em>, le sujet n&deg; 47 (<em>D<sub>47<\/sub> = 0,112<\/em>, combinant un r\u00e9sidu standardis\u00e9 de <em>z = &minus;3,41<\/em> et un levier de <em>h = 0,058<\/em>) et le sujet n&deg; 89 (<em>D<sub>89<\/sub> = 0,089<\/em>, pr\u00e9sentant un r\u00e9sidu de <em>z = &minus;2,95<\/em> et un levier de <em>h = 0,061<\/em>) p\u00e8sent de mani\u00e8re disproportionn\u00e9e sur le calcul de la pente. Un audit des dossiers d&rsquo;encodage confirme que ces deux cas correspondent \u00e0 des personnels soignants soumis \u00e0 des protocoles de travail de nuit atypiques, non contr\u00f4l\u00e9s dans le devis d&rsquo;\u00e9chantillonnage initial.<\/p>\n<h3>12.3 Ajustement du mod\u00e8le et interpr\u00e9tation scientifique finale<\/h3>\n<p>Face \u00e0 cette double violation des postulats d&rsquo;homosc\u00e9dasticit\u00e9 et de normalit\u00e9, assortie de la pr\u00e9sence d&rsquo;observations influentes document\u00e9es, l&rsquo;ajustement du mod\u00e8le devient m\u00e9thodologiquement n\u00e9cessaire. Le protocole d&rsquo;assainissement s&rsquo;articule en deux d\u00e9marches compl\u00e9mentaires :<\/p>\n<p>Dans un premier temps, une r\u00e9gression robuste utilisant un estimateur M de Huber est mise en \u0153uvre afin de limiter l&rsquo;influence disproportionn\u00e9e des observations n&deg; 47 et n&deg; 89, tout en conservant l&rsquo;int\u00e9gralit\u00e9 de l&rsquo;\u00e9chantillon au sein de l&rsquo;analyse. Dans un second temps, pour stabiliser la variance r\u00e9siduelle croissante sans distordre l&rsquo;\u00e9chelle psychom\u00e9trique, des erreurs types robustes \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 (estimateur sandwich HC3 de MacKinnon et White) sont mobilis\u00e9es lors du calcul des inf\u00e9rences statistiques.<\/p>\n<p>Le recalcul des coefficients sous ce cadre analytique robuste produit la fonction d&rsquo;estimation r\u00e9vis\u00e9e suivante :<\/p>\n<p><em>Empathie&#770; = 114,20 &minus; 0,48 &times; Burnout<\/em><\/p>\n<p>L&rsquo;\u00e9valuation des nouveaux param\u00e8tres met en \u00e9vidence une modification notable de la pente : celle-ci passe de &minus;0,62 \u00e0 &minus;0,48 point par unit\u00e9 d&rsquo;\u00e9puisement. L&rsquo;erreur type du coefficient de pente, initialement sous-estim\u00e9e sous les MCO ordinaires (<em>SE<\/em> = 0,116), s&rsquo;\u00e9tablit \u00e0 <em>SE<sub>HC3<\/sub> = 0,094<\/em> apr\u00e8s ajustement robuste, confirmant n\u00e9anmoins la significativit\u00e9 du lien empirique (<em>t = &minus;5,10, p &lt; .001, IC \u00e0 95% [&minus;0,67 ; &minus;0,29]<\/em>).<\/p>\n<p>L&rsquo;audit des r\u00e9sidus recalcul\u00e9s sur le mod\u00e8le robuste d\u00e9montre une r\u00e9sorption satisfaisante des anomalies distributionnelles. Le trac\u00e9 <em>Residuals vs Fitted<\/em> stabilis\u00e9 ne manifeste plus d&rsquo;\u00e9vasement h\u00e9t\u00e9rosc\u00e9dastique marqu\u00e9. De m\u00eame, les r\u00e9sidus studentis\u00e9s externes issus du recalibrage convergent vers l&rsquo;alignement sur la diagonale th\u00e9orique du diagramme Q-Q plot.<\/p>\n<p>La r\u00e9daction acad\u00e9mique finale peut d\u00e8s lors consigner ces r\u00e9sultats selon les normes \u00e9ditoriales les plus strictes :<\/p>\n<blockquote>\n<p>&laquo; Une analyse de r\u00e9gression lin\u00e9aire a \u00e9t\u00e9 conduite afin de pr\u00e9dire le niveau d&#8217;empathie clinique \u00e0 partir du score d&rsquo;\u00e9puisement professionnel chez 120 soignants en soins critiques. L&rsquo;audit initial des r\u00e9sidus par la m\u00e9thode des moindres carr\u00e9s ordinaires a r\u00e9v\u00e9l\u00e9 une h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 significative (test de Breusch-Pagan : <em>BP = 14,82, p &lt; .001<\/em>) ainsi qu&rsquo;une influence disproportionn\u00e9e exerc\u00e9e par deux cas singuliers identifi\u00e9s par la distance de Cook (<em>D<sub>i<\/sub> &gt; 4\/n<\/em>). L&rsquo;ajustement a \u00e9t\u00e9 reconduit via une r\u00e9gression robuste utilisant l&rsquo;estimateur de Huber, assortie d&rsquo;une correction d&rsquo;erreur type h\u00e9t\u00e9rosc\u00e9dastique HC3. Le mod\u00e8le final indique une diminution statistiquement significative de l&#8217;empathie clinique en fonction de l&rsquo;\u00e9puisement \u00e9motionnel (<em>b = &minus;0,48, SE = 0,094, t = &minus;5,10, p &lt; .001, IC \u00e0 95 % [&minus;0,67 ; &minus;0,29]<\/em>). Cette d\u00e9marche garantit la validit\u00e9 de l&rsquo;inf\u00e9rence th\u00e9orique en neutralisant les biais d&rsquo;estimation document\u00e9s lors du diagnostic des r\u00e9sidus. &raquo;<\/p>\n<\/blockquote>\n<p>Cette \u00e9tude de cas illustre ainsi concr\u00e8tement la d\u00e9marche d&rsquo;analyse r\u00e9siduelle : le calcul m\u00e9ticuleux des r\u00e9sidus n&rsquo;a pas constitu\u00e9 une simple formalit\u00e9 technique d&rsquo;ajustement, mais un filtre analytique indispensable pour pr\u00e9server la recherche d&rsquo;une mauvaise sp\u00e9cification statistique et d\u00e9livrer des conclusions solides \u00e0 la communaut\u00e9 scientifique.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Belsley, D. A., Kuh, E., &amp; Welsch, R. E. (1980). <em>Regression diagnostics: Identifying influential data and sources of collinearity<\/em>. John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471725153\">https:\/\/doi.org\/10.1002\/0471725153<\/a><\/li>\n<li>Breusch, T. S., &amp; Pagan, A. R. (1979). A simple test for heteroscedasticity and random coefficient variation. <em>Econometrica<\/em>, 47(5), 1287-1294. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1911963\">https:\/\/doi.org\/10.2307\/1911963<\/a><\/li>\n<li>Cook, R. D. (1977). Detection of influential observation in linear regression. <em>Technometrics<\/em>, 19(1), 15-18. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00401706.1977.10489493\">https:\/\/doi.org\/10.1080\/00401706.1977.10489493<\/a><\/li>\n<li>Draper, N. R., &amp; Smith, H. (1998). <em>Applied regression analysis<\/em> (3rd ed.). John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/9781118625590\">https:\/\/doi.org\/10.1002\/9781118625590<\/a><\/li>\n<li>Durbin, J., &amp; Watson, G. S. (1951). Testing for serial correlation in least squares regression. II. <em>Biometrika<\/em>, 38(1\/2), 159-177. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2332325\">https:\/\/doi.org\/10.2307\/2332325<\/a><\/li>\n<li>Fox, J., &amp; Weisberg, S. (2019). <em>An R companion to applied regression<\/em> (3rd ed.). SAGE Publications. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/socialsciences.mcmaster.ca\/jfox\/Books\/Companion\/\">https:\/\/socialsciences.mcmaster.ca\/jfox\/Books\/Companion\/<\/a><\/li>\n<li>Huber, P. J. (1981). <em>Robust statistics<\/em>. John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471725250\">https:\/\/doi.org\/10.1002\/0471725250<\/a><\/li>\n<li>MacKinnon, J. G., &amp; White, H. (1985). Some heteroskedasticity-consistent covariance matrix estimators with improved finite sample properties. <em>Journal of Econometrics<\/em>, 29(3), 305-325. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/0304-4076(85)90158-7\">https:\/\/doi.org\/10.1016\/0304-4076(85)90158-7<\/a><\/li>\n<li>Shapiro, S. S., &amp; Wilk, M. B. (1965). An analysis of variance test for normality (complete samples). <em>Biometrika<\/em>, 52(3\/4), 591-611. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2333709\">https:\/\/doi.org\/10.2307\/2333709<\/a><\/li>\n<li>White, H. (1980). A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity. <em>Econometrica<\/em>, 48(4), 817-838. <a target=\"_blank\" rel=\"noopener noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1912934\">https:\/\/doi.org\/10.2307\/1912934<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide m\u00e9thodologique complet sur le calcul math\u00e9matique, l&rsquo;analyse diagnostique et l&rsquo;interpr\u00e9tation des r\u00e9sidus de r\u00e9gression en sciences psychologiques.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1881,1835],"tags":[3896,3126,2230,16,8012],"class_list":["post-6309","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-de-la-recherche","category-statistiques-en-psychologie","tag-analyse-de-regression","tag-diagnostic-statistique","tag-moindres-carres-ordinaires","tag-psychometrie","tag-residus-de-regression"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6309","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=6309"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6309\/revisions"}],"predecessor-version":[{"id":6318,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6309\/revisions\/6318"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=6309"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=6309"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=6309"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}