{"id":5541,"date":"2026-09-11T22:39:24","date_gmt":"2026-09-11T22:39:24","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-residus-studentises-python-2\/"},"modified":"2026-09-11T22:39:25","modified_gmt":"2026-09-11T22:39:25","slug":"comment-calculer-residus-studentises-python-2","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-residus-studentises-python-2","title":{"rendered":"Comment calculer les r\u00e9sidus studentis\u00e9s en Python"},"content":{"rendered":"<p>L&rsquo;\u00e9valuation de la validit\u00e9 d&rsquo;un mod\u00e8le de r\u00e9gression lin\u00e9aire classique repose sur un ensemble d&rsquo;hypoth\u00e8ses distributionnelles et structurelles qui, lorsqu&rsquo;elles sont viol\u00e9es, compromettent gravement l&rsquo;inf\u00e9rence statistique. Parmi ces postulats figurent la lin\u00e9arit\u00e9 de l&rsquo;esp\u00e9rance conditionnelle, l&rsquo;homosc\u00e9dasticit\u00e9 de la variance d&rsquo;erreur, l&rsquo;ind\u00e9pendance des perturbations stochastiques et leur normalit\u00e9 asymptotique. Dans la pratique de la mod\u00e9lisation statistique appliqu\u00e9e aux sciences empiriques, l&rsquo;analyste se trouve confront\u00e9 \u00e0 des jeux de donn\u00e9es r\u00e9els o\u00f9 certaines observations s&rsquo;\u00e9cartent substantiellement de la tendance g\u00e9n\u00e9rale. Si les r\u00e9sidus bruts constituent le point de d\u00e9part intuitif pour quantifier ces \u00e9carts, leur utilisation brute pour tester des hypoth\u00e8ses de valeurs aberrantes s&rsquo;av\u00e8re m\u00e9thodologiquement d\u00e9ficiente en raison de leur h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 intrins\u00e8que de variance.<\/p>\n<p>Cette limite fondamentale a motiv\u00e9 le d\u00e9veloppement de r\u00e9sidus standardis\u00e9s et plus particuli\u00e8rement des r\u00e9sidus studentis\u00e9s. En divisant l&rsquo;\u00e9cart d&rsquo;ajustement par un estimateur rigoureux de son \u00e9cart-type conditionnel \u00e0 la position de l&rsquo;observation dans l&rsquo;espace multidimensionnel des r\u00e9gresseurs, la studentisation conf\u00e8re au r\u00e9sidu une m\u00e9trique adimensionnelle et une distribution th\u00e9orique connue sous l&rsquo;hypoth\u00e8se nulle. La distinction cardinale entre studentisation interne et studentisation externe (ou jackknife) repr\u00e9sente une avanc\u00e9e majeure dans la d\u00e9tection des points atypiques, permettant de neutraliser l&rsquo;effet d&rsquo;att\u00e9nuation qu&rsquo;une observation extr\u00eame exerce sur l&rsquo;estimation de sa propre variance d&rsquo;erreur.<\/p>\n<p>L&rsquo;\u00e9cosyst\u00e8me scientifique du langage Python, structur\u00e9 autour de biblioth\u00e8ques performantes telles que NumPy, Pandas et Statsmodels, fournit un cadre computationnel id\u00e9al pour op\u00e9rationnaliser ces diagnostics avanc\u00e9s. Ce guide exhaustif explore l&rsquo;ensemble des fondements th\u00e9oriques, des formulations math\u00e9matiques et des impl\u00e9mentations algorithmiques n\u00e9cessaires au calcul, \u00e0 l&rsquo;interpr\u00e9tation et \u00e0 la visualisation des r\u00e9sidus studentis\u00e9s en Python. \u00c0 travers des cas d&rsquo;application concrets issus des sciences du comportement et de la mod\u00e9lisation quantitative, nous d\u00e9taillerons les protocoles m\u00e9thodologiques assurant une inf\u00e9rence robuste et reproductible.<\/p>\n<h2>1. Introduction th\u00e9orique aux r\u00e9sidus studentis\u00e9s dans l&rsquo;analyse de r\u00e9gression<\/h2>\n<p>### 1.1 D\u00e9finition fondamentale et formalisation math\u00e9matique<\/p>\n<p>Dans le cadre d&rsquo;un mod\u00e8le de r\u00e9gression lin\u00e9aire multiple exprim\u00e9 sous la forme matricielle canonique <span class=\"arabpsy-math-inline\" data-latex=\"Y = X\\beta + \\varepsilon\">$Y = X\\beta + \\varepsilon$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> d\u00e9signe le vecteur des r\u00e9ponses de dimension <span class=\"arabpsy-math-inline\" data-latex=\"n \\times 1\">$n \\times 1$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> la matrice de conception de dimension <span class=\"arabpsy-math-inline\" data-latex=\"n \\times p\">$n \\times p$<\/span> incluant la constante, <span class=\"arabpsy-math-inline\" data-latex=\"\\beta\">$\\beta$<\/span> le vecteur des param\u00e8tres inconnus de dimension <span class=\"arabpsy-math-inline\" data-latex=\"p \\times 1\">$p \\times 1$<\/span>, et <span class=\"arabpsy-math-inline\" data-latex=\"varepsilon\">$varepsilon$<\/span> le vecteur des erreurs stochastiques non observables suivant une distribution gaussienne <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{N}(0, \\sigma^2 I_n)\">$\\mathcal{N}(0, \\sigma^2 I_n)$<\/span>, le vecteur des r\u00e9sidus ordinaires (ou bruts) est d\u00e9fini par la diff\u00e9rence alg\u00e9brique entre les valeurs observ\u00e9es et les valeurs ajust\u00e9es par la m\u00e9thode des moindres carr\u00e9s ordinaires (MCO) :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"e = Y - \\hat{Y} = Y - X\\hat{\\beta} = (I_n - H)Y\">$e = Y &#8211; \\hat{Y} = Y &#8211; X\\hat{\\beta} = (I_n &#8211; H)Y$<\/span><\/p>\n<p>L&rsquo;op\u00e9rateur <span class=\"arabpsy-math-inline\" data-latex=\"H = X(X^T X)^{-1} X^T\">$H = X(X^T X)^{-1} X^T$<\/span> repr\u00e9sente la matrice de projection orthogonale sur l&rsquo;espace vectoriel engendr\u00e9 par les colonnes de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span>, commun\u00e9ment appel\u00e9e matrice chapeau (*hat matrix*). Bien que l&rsquo;hypoth\u00e8se classique postule que les erreurs th\u00e9oriques <span class=\"arabpsy-math-inline\" data-latex=\"\\varepsilon_i\">$\\varepsilon_i$<\/span> partagent une variance commune et constante <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span>, les r\u00e9sidus empiriques <span class=\"arabpsy-math-inline\" data-latex=\"e_i\">$e_i$<\/span> ne poss\u00e8dent pas cette propri\u00e9t\u00e9 d&rsquo;homosc\u00e9dasticit\u00e9. En d\u00e9veloppant la matrice de covariance du vecteur r\u00e9siduel <span class=\"arabpsy-math-inline\" data-latex=\"e\">$e$<\/span>, nous obtenons :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\text{Var}(e) = \\text{Var}((I_n - H)Y) = (I_n - H) \\text{Var}(Y) (I_n - H)^T = \\sigma^2 (I_n - H)\">$\\text{Var}(e) = \\text{Var}((I_n &#8211; H)Y) = (I_n &#8211; H) \\text{Var}(Y) (I_n &#8211; H)^T = \\sigma^2 (I_n &#8211; H)$<\/span><\/p>\n<p>Il en d\u00e9coule directement que la variance th\u00e9orique du <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me r\u00e9sidu brut d\u00e9pend explicitement de la position occup\u00e9e par l&rsquo;observation correspondante dans l&rsquo;espace des covariables, mat\u00e9rialis\u00e9e par le <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me \u00e9l\u00e9ment diagonal <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> de la matrice <span class=\"arabpsy-math-inline\" data-latex=\"H\">$H$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\text{Var}(e_i) = \\sigma^2 (1 - h_{ii})\">$\\text{Var}(e_i) = \\sigma^2 (1 &#8211; h_{ii})$<\/span><\/p>\n<p>Cette d\u00e9pendance met en lumi\u00e8re les limites intrins\u00e8ques des r\u00e9sidus bruts comme outil de diagnostic comparatif. Les observations situ\u00e9es aux confins de la distribution des variables explicatives pr\u00e9sentent une valeur de levier <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> \u00e9lev\u00e9e, ce qui comprime m\u00e9caniquement la variance de leur r\u00e9sidu brut vers z\u00e9ro. Par cons\u00e9quent, un point poss\u00e9dant un fort effet de levier peut forcer la droite de r\u00e9gression \u00e0 passer excessivement pr\u00e8s de lui, affichant un r\u00e9sidu brut artificiellement modeste alors m\u00eame qu&rsquo;il d\u00e9stabilise profond\u00e9ment l&rsquo;estimation des coefficients <span class=\"arabpsy-math-inline\" data-latex=\"\\beta\">$\\beta$<\/span>. Pour rem\u00e9dier \u00e0 cette inhomog\u00e9n\u00e9it\u00e9, la proc\u00e9dure de studentisation consiste \u00e0 diviser le r\u00e9sidu brut par l&rsquo;estimation de son \u00e9cart-type exact, g\u00e9n\u00e9rant ainsi une statistique d&rsquo;ajustement normalis\u00e9e et adimensionnelle.<\/p>\n<p>Sous l&rsquo;hypoth\u00e8se de normalit\u00e9 des perturbations stochastiques <span class=\"arabpsy-math-inline\" data-latex=\"\\varepsilon \\sim \\mathcal{N}(0, \\sigma^2 I_n)\">$\\varepsilon \\sim \\mathcal{N}(0, \\sigma^2 I_n)$<\/span>, la transformation normalis\u00e9e par l&rsquo;\u00e9cart-type r\u00e9siduel ne suit pas une loi normale standard en raison de l&rsquo;estimation de la variance d&rsquo;erreur par un carr\u00e9 moyen r\u00e9siduel comportant un nombre fini de degr\u00e9s de libert\u00e9. Selon les travaux fondamentaux de William Sealy Gosset (publi\u00e9s sous le pseudonyme de <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2331554\">Student en 1908<\/a>), ce ratio d&rsquo;une variable gaussienne par la racine carr\u00e9e d&rsquo;une variable du khi-deux ind\u00e9pendante divis\u00e9e par ses degr\u00e9s de libert\u00e9 engendre une distribution de Student. Cette propri\u00e9t\u00e9 formelle offre une base probabiliste rigoureuse pour construire des tests de signification d&rsquo;atypicit\u00e9 statistique.<\/p>\n<p>### 1.2 R\u00e9sidus studentis\u00e9s internes versus externes<\/p>\n<p>La litt\u00e9rature \u00e9conom\u00e9trique et statistique \u00e9tablit une distinction m\u00e9thodologique imp\u00e9rative entre les r\u00e9sidus studentis\u00e9s internes (parfois qualifi\u00e9s de r\u00e9sidus standardis\u00e9s) et les r\u00e9sidus studentis\u00e9s externes (\u00e9galement appel\u00e9s r\u00e9sidus jackknife, r\u00e9sidus d&rsquo;exclusion ou r\u00e9sidus *leave-one-out*). Cette divergence repose fondamentalement sur le sous-\u00e9chantillon utilis\u00e9 pour estimer la variance r\u00e9siduelle inconnue <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span>.<\/p>\n<p>Le r\u00e9sidu studentis\u00e9 interne, couramment not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"r_i\">$r_i$<\/span>, s&rsquo;obtient en utilisant l&rsquo;estimateur sans biais global de la variance d&rsquo;erreur calcul\u00e9 sur l&rsquo;int\u00e9gralit\u00e9 des <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> observations de l&rsquo;\u00e9chantillon, soit le carr\u00e9 moyen de l&rsquo;erreur <span class=\"arabpsy-math-inline\" data-latex=\"s^2 = \\frac{e^T e}{n - p} = \\frac{SSE}{n - p}\">$s^2 = \\frac{e^T e}{n &#8211; p} = \\frac{SSE}{n &#8211; p}$<\/span>. Sa formulation analytique s&rsquo;\u00e9crit :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"r_i = \\frac{e_i}{s \\sqrt{1 - h_{ii}}}\">$r_i = \\frac{e_i}{s \\sqrt{1 &#8211; h_{ii}}}$<\/span><\/p>\n<p>Bien que cette m\u00e9trique corrige l&rsquo;effet de structure imputable \u00e0 la matrice chapeau via le facteur <span class=\"arabpsy-math-inline\" data-latex=\"\\sqrt{1 - h_{ii}}\">$\\sqrt{1 &#8211; h_{ii}}$<\/span>, elle souffre d&rsquo;un biais d&rsquo;auto-contamination. Si la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me observation est une valeur aberrante extr\u00eame (un *outlier* majeur), son propre r\u00e9sidu brut <span class=\"arabpsy-math-inline\" data-latex=\"e_i\">$e_i$<\/span> de grande amplitude contribue de mani\u00e8re disproportionn\u00e9e \u00e0 la somme des carr\u00e9s des erreurs r\u00e9siduelles (<span class=\"arabpsy-math-inline\" data-latex=\"SSE\">$SSE$<\/span>). Cette inflation m\u00e9canique de <span class=\"arabpsy-math-inline\" data-latex=\"SSE\">$SSE$<\/span> gonfle artificiellement l&rsquo;estimateur global <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>, ce qui a pour effet pervers de r\u00e9duire le d\u00e9nominateur de <span class=\"arabpsy-math-inline\" data-latex=\"r_i\">$r_i$<\/span>. Ainsi, le r\u00e9sidu studentis\u00e9 interne d&rsquo;une observation aberrante se trouve auto-masqu\u00e9 par sa propre pr\u00e9sence dans le calcul de la dispersion globale, diminuant la puissance diagnostique du test.<\/p>\n<p>Pour contourner ce ph\u00e9nom\u00e8ne d&rsquo;att\u00e9nuation, le r\u00e9sidu studentis\u00e9 externe, not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>, substitue \u00e0 l&rsquo;estimateur global <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span> un estimateur d&rsquo;exclusion not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}\">$s_{(-i)}$<\/span>, calcul\u00e9 sur le jeu de donn\u00e9es duquel la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me observation a \u00e9t\u00e9 pr\u00e9alablement retranch\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"t_i = \\frac{e_i}{s_{(-i)} \\sqrt{1 - h_{ii}}}\">$t_i = \\frac{e_i}{s_{(-i)} \\sqrt{1 &#8211; h_{ii}}}$<\/span><\/p>\n<p>Dans cette configuration, l&rsquo;estimateur <span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2\">$s_{(-i)}^2$<\/span> est calcul\u00e9 sur <span class=\"arabpsy-math-inline\" data-latex=\"n - 1\">$n &#8211; 1$<\/span> observations et <span class=\"arabpsy-math-inline\" data-latex=\"n - p - 1\">$n &#8211; p &#8211; 1$<\/span> degr\u00e9s de libert\u00e9. Par construction, si l&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> est fortement discordante, son exclusion \u00e9limine son influence contaminatrice sur l&rsquo;estimation de la variance de fond du mod\u00e8le. L&rsquo;\u00e9cart-type d&rsquo;erreur <span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}\">$s_{(-i)}$<\/span> refl\u00e8te fid\u00e8lement la variabilit\u00e9 des donn\u00e9es non contamin\u00e9es, ce qui propulse le ratio <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span> vers des valeurs d&rsquo;amplitude consid\u00e9rablement plus \u00e9lev\u00e9es.<\/p>\n<p>Sur le plan strictement distributionnel, tandis que le r\u00e9sidu interne <span class=\"arabpsy-math-inline\" data-latex=\"r_i^2 \/ (n - p)\">$r_i^2 \/ (n &#8211; p)$<\/span> suit une loi B\u00eata en raison de la d\u00e9pendance stochastique entre son num\u00e9rateur <span class=\"arabpsy-math-inline\" data-latex=\"e_i\">$e_i$<\/span> et le d\u00e9nominateur <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>, le r\u00e9sidu externe <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span> est rigoureusement distribu\u00e9 selon une loi de Student \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"n - p - 1\">$n &#8211; p &#8211; 1$<\/span> degr\u00e9s de libert\u00e9 sous l&rsquo;hypoth\u00e8se nulle stipulant que l&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> est issue du m\u00eame processus g\u00e9n\u00e9rateur que le reste de la cohorte. Cette propri\u00e9t\u00e9 statistique exacte justifie la sup\u00e9riorit\u00e9 des r\u00e9sidus studentis\u00e9s externes pour les tests d&rsquo;hypoth\u00e8se formels relatifs aux valeurs aberrantes dans les pipelines d&rsquo;audit diagnostique.<\/p>\n<p>### 1.3 R\u00f4le des r\u00e9sidus studentis\u00e9s dans l&rsquo;inf\u00e9rence statistique<\/p>\n<p>L&rsquo;analyse des r\u00e9sidus studentis\u00e9s remplit une fonction heuristique et confirmatoire centrale dans l&rsquo;\u00e9valuation de la robustesse des mod\u00e8les de r\u00e9gression appliqu\u00e9s aux sciences empiriques. Le postulat fondamental d&rsquo;homosc\u00e9dasticit\u00e9 exige que l&rsquo;erreur stochastique poss\u00e8de une variance constante <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2\">$\\sigma^2$<\/span> sur l&rsquo;ensemble du domaine de pr\u00e9diction. Lorsque l&rsquo;analyste examine un diagramme de dispersion opposant les r\u00e9sidus bruts aux valeurs pr\u00e9dites, la contraction structurelle de la variance induite par les points de fort levier peut mimer ou au contraire occulter une structure h\u00e9t\u00e9rosc\u00e9dastique r\u00e9elle. Les r\u00e9sidus studentis\u00e9s, en \u00e9liminant cette h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 artefactuelle li\u00e9e \u00e0 la configuration spatiale des covariables, permettent d&rsquo;isoler avec une fid\u00e9lit\u00e9 optimale les v\u00e9ritables distorsions de la variance r\u00e9siduelle, telles que les patrons en entonnoir repr\u00e9sentatifs d&rsquo;une variance proportionnelle \u00e0 l&rsquo;amplitude des variables d\u00e9pendantes.<\/p>\n<p>Sur le plan de la g\u00e9om\u00e9trie multidimensionnelle, l&rsquo;inf\u00e9rence moderne distingue formellement les points \u00e0 fort r\u00e9sidu des points \u00e0 fort levier. Une observation peut afficher une position excentrique dans l&rsquo;espace factoriel des variables explicatives (levier \u00e9lev\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span>) tout en s&rsquo;alignant scrupuleusement sur la trajectoire lin\u00e9aire du mod\u00e8le ; un tel point poss\u00e8de un r\u00e9sidu studentis\u00e9 proche de z\u00e9ro et ne perturbe nullement l&rsquo;estimation des coefficients directeurs. \u00c0 l&rsquo;inverse, une observation peut se situer \u00e0 proximit\u00e9 imm\u00e9diate du centro\u00efde des r\u00e9gresseurs (faible levier) mais d\u00e9vier drastiquement le long de l&rsquo;axe de la variable d\u00e9pendante <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>. Les r\u00e9sidus studentis\u00e9s externes constituent l&rsquo;unique sonde statistique capable de quantifier cette discordance verticale ind\u00e9pendamment du positionnement horizontal de l&rsquo;observation.<\/p>\n<p>Dans les protocoles de recherche quantitative en sciences m\u00e9dicales, \u00e9conomiques et sociales, la non-d\u00e9tection d&rsquo;observations aberrantes menace directement la validit\u00e9 interne des inf\u00e9rences causales et corr\u00e9lationnelles. La m\u00e9thode des moindres carr\u00e9s ordinaires minimise la somme des \u00e9carts quadratiques, ce qui octroie un poids quadratiquement proportionnel aux d\u00e9viations extr\u00eames. Une seule observation erron\u00e9e non identifi\u00e9e peut infl\u00e9chir la pente d&rsquo;un coefficient <span class=\"arabpsy-math-inline\" data-latex=\"\\beta\">$\\beta$<\/span>, inverser son signe alg\u00e9brique, ou faire basculer artificiellement une valeur <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> en de\u00e7\u00e0 ou au-del\u00e0 du seuil de significativit\u00e9 conventionnel <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,05\">$\\alpha = 0,05$<\/span>. L&rsquo;\u00e9valuation syst\u00e9matique des r\u00e9sidus studentis\u00e9s pr\u00e9serve la stabilit\u00e9 structurelle du mod\u00e8le en garantissant que les conclusions th\u00e9oriques tir\u00e9es des donn\u00e9es ne sont pas l&rsquo;art\u00e9fact d&rsquo;une fraction infinit\u00e9simale d&rsquo;observations pathologiques.<\/p>\n<p>## 2. Pertinence diagnostique des r\u00e9sidus studentis\u00e9s en psychom\u00e9trie et sciences du comportement<\/p>\n<p>### 2.1 Vuln\u00e9rabilit\u00e9 des donn\u00e9es psychologiques aux observations atypiques<\/p>\n<p>Les sciences du comportement, la neuropsychologie et la psychom\u00e9trie contemporaine op\u00e8rent par nature sur des donn\u00e9es caract\u00e9ris\u00e9es par un ratio signal-sur-bruit intrins\u00e8quement modeste et une vuln\u00e9rabilit\u00e9 aigu\u00eb aux perturbations m\u00e9thodologiques. Lors de l&rsquo;administration de protocoles exp\u00e9rimentaux chronom\u00e9tr\u00e9s ou de batteries d&rsquo;\u00e9valuation psychologique standardis\u00e9es (telles que les \u00e9chelles de Wechsler ou les inventaires de personnalit\u00e9 NEO-PI), de multiples facteurs parasites peuvent corrompre l&rsquo;acquisition des donn\u00e9es. Les biais d&rsquo;acquiescement, la fatigue attentionnelle, les comportements de r\u00e9ponse al\u00e9atoire (*careless responding*), la falsification d\u00e9lib\u00e9r\u00e9e, ou de simples d\u00e9faillances logicielles d&rsquo;enregistrement des temps de r\u00e9action engendrent des profils de r\u00e9ponse aberrants qui s&rsquo;\u00e9cartent radicalement des pr\u00e9dictions th\u00e9oriques des mod\u00e8les lin\u00e9aires.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/studres1-1.png\" alt=\"Studentized residuals in Python\" \/><figcaption>Studentized residuals in Python<\/figcaption><\/figure>\n<p>Il est imp\u00e9ratif pour le chercheur d&rsquo;\u00e9tablir une distinction \u00e9pist\u00e9mologique stricte entre une valeur atypique issue d&rsquo;une erreur d&rsquo;\u00e9chantillonnage ou de mesure, et une valeur extr\u00eame l\u00e9gitime refl\u00e9tant la variabilit\u00e9 ph\u00e9notypique r\u00e9elle de la population cible. Dans une cohorte clinique examinant le d\u00e9clin cognitif li\u00e9 au vieillissement, un score anormalement bas aux \u00e9preuves de m\u00e9moire de travail peut constituer un signal neurophysiologique valide de prodrome neurod\u00e9g\u00e9n\u00e9ratif, ou d\u00e9couler d&rsquo;une incompr\u00e9hension des consignes par le participant. L&rsquo;utilisation non critique des tests de rejet peut induire une inflation dramatique de l&rsquo;erreur de type I (faux positifs) par suppression injustifi\u00e9e de la variance \u00e9cologique, ou de l&rsquo;erreur de type II (faux n\u00e9gatifs) par conservation de points aberrants qui \u00e9crasent la puissance statistique des tests d&rsquo;interaction.<\/p>\n<p>Les comit\u00e9s de r\u00e9daction des revues scientifiques de premier rang exigent d\u00e9sormais une transparence m\u00e9thodologique absolue quant au traitement des donn\u00e9es divergentes, conform\u00e9ment aux lignes directrices des protocoles de r\u00e9plication et aux normes de l&rsquo;American Psychological Association (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\">APA Reporting Standards<\/a>). La studentisation externe des r\u00e9sidus fournit un ancrage statistique probabiliste objectif qui affranchit l&rsquo;investigateur de l&rsquo;arbitraire des r\u00e8gles empiriques subjectives, garantissant une d\u00e9tection standardis\u00e9e, quantifiable et auditable des d\u00e9viations majeures au mod\u00e8le psychom\u00e9trique ajust\u00e9.<\/p>\n<p>### 2.2 Levier, distance et influence dans les cohortes comportementales<\/p>\n<p>L&rsquo;analyse de r\u00e9gression appliqu\u00e9e aux cohortes comportementales requiert une articulation conceptuelle pr\u00e9cise entre les notions de levier (*leverage*), d&rsquo;atypicit\u00e9 r\u00e9siduelle (*outlyingness*) et d&rsquo;influence globale (*influence*). Comme expos\u00e9 pr\u00e9c\u00e9demment, la mesure de levier, codifi\u00e9e par la valeur <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> extraite de la diagonale de la matrice de projection, quantifie l&rsquo;\u00e9loignement d&rsquo;un individu statistique par rapport au barycentre des covariables explicatives de l&rsquo;\u00e9chantillon. Un participant pr\u00e9sentant un profil sociod\u00e9mographique ou neuropsychologique exceptionnel \u2014 par exemple un individu affichant un \u00e2ge tr\u00e8s avanc\u00e9 coupl\u00e9 \u00e0 un niveau de scolarit\u00e9 atypiquement \u00e9lev\u00e9 au sein d&rsquo;une \u00e9tude sur l&rsquo;anxi\u00e9t\u00e9 \u2014 obtiendra une valeur <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> substantielle, sans que cela n&rsquo;implique n\u00e9cessairement une discordance comportementale sur la variable pr\u00e9dite.<\/p>\n<p>L&rsquo;\u00e9valuation de l&rsquo;influence d\u00e9terminante exerc\u00e9e par une observation sur l&rsquo;ensemble de la topologie du mod\u00e8le n\u00e9cessite de combiner l&rsquo;information port\u00e9e par le r\u00e9sidu studentis\u00e9 et celle port\u00e9e par le levier. Cette synth\u00e8se math\u00e9matique s&rsquo;incarne notamment dans la distance de Cook (<span class=\"arabpsy-math-inline\" data-latex=\"D_i\">$D_i$<\/span>), dont l&rsquo;expression analytique peut \u00eatre directement formul\u00e9e en fonction du r\u00e9sidu studentis\u00e9 interne <span class=\"arabpsy-math-inline\" data-latex=\"r_i\">$r_i$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"D_i = \\frac{r_i^2}{p} \\left( \\frac{h_{ii}}{1 - h_{ii}} \\right)\">$D_i = \\frac{r_i^2}{p} \\left( \\frac{h_{ii}}{1 &#8211; h_{ii}} \\right)$<\/span><\/p>\n<p>Cette formulation met en exergue la dynamique multiplicative de l&rsquo;influence statistique : une observation n&rsquo;exerce une modification substantielle des estimations param\u00e9triques globales que si elle cumule simultan\u00e9ment un r\u00e9sidu studentis\u00e9 d&rsquo;amplitude marqu\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"r_i^2\">$r_i^2$<\/span>) et une position d&rsquo;ancrage excentrique dans l&rsquo;espace factoriel des pr\u00e9dicteurs (<span class=\"arabpsy-math-inline\" data-latex=\"h_{ii} \/ (1 - h_{ii})\">$h_{ii} \/ (1 &#8211; h_{ii})$<\/span>). <\/p>\n<p>Consid\u00e9rons \u00e0 titre d&rsquo;illustration une \u00e9tude psychom\u00e9trique analysant la relation entre le niveau de stress per\u00e7u et la performance \u00e0 une t\u00e2che de flexibilit\u00e9 cognitive chez des adultes neurotypiques et des individus pr\u00e9sentant un profil neurodivergent atypique (par exemple, un trouble du spectre de l&rsquo;autisme \u00e0 haut niveau de fonctionnement associ\u00e9 \u00e0 une synesth\u00e9sie). Si cet individu atypique obtient des scores de stress mod\u00e9r\u00e9s mais r\u00e9ussit la t\u00e2che cognitive avec un score stratosph\u00e9rique d\u00e9fiant les patterns habituels, son r\u00e9sidu studentis\u00e9 externe <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span> atteindra des grandeurs extr\u00eames. Si son profil de covariables est par ailleurs standard, son levier restera mod\u00e9r\u00e9, limitant son impact d\u00e9stabilisateur sur les coefficients directeurs g\u00e9n\u00e9raux, tout en signalant sans \u00e9quivoque son statut de singularit\u00e9 clinique n\u00e9cessitant une analyse qualitative distincte.<\/p>\n<p>L&rsquo;arbitrage entre l&rsquo;exclusion d&rsquo;un profil et son maintien dans le mod\u00e8le final conditionne la repr\u00e9sentativit\u00e9 \u00e9cologique des r\u00e9sultats. L&rsquo;analyste comportemental doit s&rsquo;appuyer sur un faisceau convergent d&rsquo;indices diagnostiques \u2014 combinant seuils critiques sur <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>, examen des variations param\u00e9triques par DFBETAS et calcul de la distance de Cook \u2014 afin de prendre des d\u00e9cisions motiv\u00e9es qui ne biaisent pas artificiellement la composition de l&rsquo;\u00e9chantillon \u00e9tudi\u00e9.<\/p>\n<p>## 3. Formulation math\u00e9matique et distribution des r\u00e9sidus studentis\u00e9s<\/p>\n<p>### 3.1 Matrice de projection et calcul des valeurs de levier<\/p>\n<p>La d\u00e9rivation alg\u00e9brique rigoureuse des r\u00e9sidus studentis\u00e9s repose sur l&rsquo;alg\u00e8bre lin\u00e9aire appliqu\u00e9e \u00e0 l&rsquo;espace de Hilbert des variables al\u00e9atoires de carr\u00e9 int\u00e9grable. Soit <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> la matrice de conception de rang plein <span class=\"arabpsy-math-inline\" data-latex=\"p \\leq n\">$p \\leq n$<\/span>. L&rsquo;op\u00e9rateur de projection orthogonale <span class=\"arabpsy-math-inline\" data-latex=\"H\">$H$<\/span>, qui projette le vecteur des observations <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> sur le sous-espace vectoriel de r\u00e9gression <span class=\"arabpsy-math-inline\" data-latex=\"\\text{Im}(X)\">$\\text{Im}(X)$<\/span>, est formellement d\u00e9fini par l&rsquo;expression :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"H = X (X^T X)^{-1} X^T\">$H = X (X^T X)^{-1} X^T$<\/span><\/p>\n<p>Cette matrice carr\u00e9e de dimension <span class=\"arabpsy-math-inline\" data-latex=\"n \\times n\">$n \\times n$<\/span> poss\u00e8de des propri\u00e9t\u00e9s g\u00e9om\u00e9triques fondamentales qui dictent le comportement des r\u00e9sidus :<\/p>\n<ul>\n<li>Elle est <strong>sym\u00e9trique<\/strong>, v\u00e9rifiant rigoureusement <span class=\"arabpsy-math-inline\" data-latex=\"H^T = H\">$H^T = H$<\/span>.<\/li>\n<li>Elle est <strong>idempotente<\/strong>, c&rsquo;est-\u00e0-dire que <span class=\"arabpsy-math-inline\" data-latex=\"H^2 = H\">$H^2 = H$<\/span>, ce qui implique g\u00e9om\u00e9triquement qu&rsquo;une projection successive r\u00e9it\u00e9r\u00e9e n&rsquo;alt\u00e8re plus le vecteur projet\u00e9.<\/li>\n<li>Sa trace, repr\u00e9sentant la somme de ses \u00e9l\u00e9ments diagonaux, est invariante et \u00e9gale \u00e0 la dimension du sous-espace d&rsquo;ajustement : <span class=\"arabpsy-math-inline\" data-latex=\"\\text{Tr}(H) = \\sum_{i=1}^n h_{ii} = \\text{rg}(X) = p\">$\\text{Tr}(H) = \\sum_{i=1}^n h_{ii} = \\text{rg}(X) = p$<\/span>.<\/li>\n<\/ul>\n<p>Les \u00e9l\u00e9ments diagonaux <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span>, commun\u00e9ment d\u00e9sign\u00e9s sous l&rsquo;appellation de valeurs de levier (*leverage values*), poss\u00e8dent un encadrement th\u00e9orique strict. En exploitant l&rsquo;idempotence et la sym\u00e9trie de la matrice <span class=\"arabpsy-math-inline\" data-latex=\"H\">$H$<\/span>, nous pouvons expliciter <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> sous forme d&rsquo;un produit scalaire interne :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"h_{ii} = \\sum_{j=1}^n h_{ij}^2 = h_{ii}^2 + \\sum_{j \\neq i} h_{ij}^2\">$h_{ii} = \\sum_{j=1}^n h_{ij}^2 = h_{ii}^2 + \\sum_{j \\neq i} h_{ij}^2$<\/span><\/p>\n<p>\u00c9tant donn\u00e9 que la somme quadratique <span class=\"arabpsy-math-inline\" data-latex=\"\\sum_{j \\neq i} h_{ij}^2\">$\\sum_{j \\neq i} h_{ij}^2$<\/span> est obligatoirement positive ou nulle, il s&rsquo;ensuit imm\u00e9diatement que <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii} \\geq h_{ii}^2\">$h_{ii} \\geq h_{ii}^2$<\/span>, ce qui impose la double in\u00e9galit\u00e9 fondamentale :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"0 \\leq h_{ii} \\leq 1\">$0 \\leq h_{ii} \\leq 1$<\/span><\/p>\n<p>Lorsque la matrice de r\u00e9gression <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> int\u00e8gre un terme constant (colonne de uns), la borne inf\u00e9rieure se resserre \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"1\/n \\leq h_{ii} \\leq 1\">$1\/n \\leq h_{ii} \\leq 1$<\/span>. Sur le plan g\u00e9om\u00e9trique, l&rsquo;\u00e9l\u00e9ment <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> mesure le carr\u00e9 de la distance de Mahalanobis standardis\u00e9e s\u00e9parant le point <span class=\"arabpsy-math-inline\" data-latex=\"x_i = (x_{i1}, x_{i2}, dots, x_{ip})^T\">$x_i = (x_{i1}, x_{i2}, dots, x_{ip})^T$<\/span> du barycentre vectoriel des variables explicatives <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{x}\">$\\bar{x}$<\/span>. Plus une observation s&rsquo;isole aux p\u00e9riph\u00e9ries de la distribution conjointe des pr\u00e9dicteurs, plus sa valeur <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> converge vers l&rsquo;unit\u00e9. <\/p>\n<p>En rappelant que la matrice de projection r\u00e9siduelle <span class=\"arabpsy-math-inline\" data-latex=\"M = I_n - H\">$M = I_n &#8211; H$<\/span> partage les m\u00eames propri\u00e9t\u00e9s de sym\u00e9trie et d&rsquo;idempotence, la variance du vecteur des r\u00e9sidus ordinaires <span class=\"arabpsy-math-inline\" data-latex=\"e = MY\">$e = MY$<\/span> se d\u00e9ploie comme suit :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\text{Var}(e) = \\text{Var}(M(X\\beta + \\varepsilon)) = M \\text{Var}(\\varepsilon) M^T = \\sigma^2 M M = \\sigma^2 (I_n - H)\">$\\text{Var}(e) = \\text{Var}(M(X\\beta + \\varepsilon)) = M \\text{Var}(\\varepsilon) M^T = \\sigma^2 M M = \\sigma^2 (I_n &#8211; H)$<\/span><\/p>\n<p>Cette \u00e9quation \u00e9tablit sans ambigu\u00eft\u00e9 que l&rsquo;\u00e9cart-type th\u00e9orique exact du r\u00e9sidu associ\u00e9 au <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me individu est donn\u00e9 par la fonction analytique <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma \\sqrt{1 - h_{ii}}\">$\\sigma \\sqrt{1 &#8211; h_{ii}}$<\/span>. Lorsque <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> approche de 1, la variance du r\u00e9sidu tend vers z\u00e9ro, confirmant math\u00e9matiquement que les points de fort levier contraignent l&rsquo;ajustement du mod\u00e8le \u00e0 s&rsquo;annuler sur leurs propres coordonn\u00e9es observ\u00e9es.<\/p>\n<p>### 3.2 Algorithme d&rsquo;exclusion et loi de Student associ\u00e9e<\/p>\n<p>Pour d\u00e9river la statistique du r\u00e9sidu studentis\u00e9 externe sans compromettre l&rsquo;ind\u00e9pendance stochastique requise pour un test d&rsquo;hypoth\u00e8se rigoureux, il convient d&rsquo;expliciter le mod\u00e8le lin\u00e9aire auxiliaire d&rsquo;exclusion. Consid\u00e9rons le mod\u00e8le de r\u00e9gression ajust\u00e9 sur l&rsquo;\u00e9chantillon tronqu\u00e9 d&rsquo;effectif <span class=\"arabpsy-math-inline\" data-latex=\"n - 1\">$n &#8211; 1$<\/span>, dans lequel la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me observation <span class=\"arabpsy-math-inline\" data-latex=\"(x_i, y_i)\">$(x_i, y_i)$<\/span> a \u00e9t\u00e9 omise. Soit <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}_{(-i)}\">$\\hat{\\beta}_{(-i)}$<\/span> le vecteur des coefficients MCO obtenu sur ce sous-\u00e9chantillon. L&rsquo;erreur de pr\u00e9diction hors-\u00e9chantillon (*deleted residual*) pour cette observation s&rsquo;\u00e9nonce :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"d_i = y_i - x_i^T \\hat{\\beta}_{(-i)}\">$d_i = y_i &#8211; x_i^T \\hat{\\beta}_{(-i)}$<\/span><\/p>\n<p>Gr\u00e2ce \u00e0 l&rsquo;identit\u00e9 alg\u00e9brique de Sherman-Morrison-Woodbury relative aux inversions de rang un des matrices modifi\u00e9es, il est possible de d\u00e9montrer que ce r\u00e9sidu de pr\u00e9diction d&rsquo;exclusion est directement proportionnel au r\u00e9sidu brut ordinaire <span class=\"arabpsy-math-inline\" data-latex=\"e_i\">$e_i$<\/span> selon la relation :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"d_i = \\frac{e_i}{1 - h_{ii}}\">$d_i = \\frac{e_i}{1 &#8211; h_{ii}}$<\/span><\/p>\n<p>La variance th\u00e9orique de cette erreur de pr\u00e9diction hors-\u00e9chantillon, conditionn\u00e9e aux r\u00e9gresseurs, int\u00e8gre l&rsquo;incertitude sur la nouvelle observation combin\u00e9e \u00e0 l&rsquo;impr\u00e9cision d&rsquo;estimation de <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}_{(-i)}\">$\\hat{\\beta}_{(-i)}$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\text{Var}(d_i) = \\text{Var}(y_i - x_i^T \\hat{\\beta}_{(-i)}) = \\sigma^2 \\left( 1 + x_i^T (X_{(-i)}^T X_{(-i)})^{-1} x_i \\right) = \\frac{\\sigma^2}{1 - h_{ii}}\">$\\text{Var}(d_i) = \\text{Var}(y_i &#8211; x_i^T \\hat{\\beta}_{(-i)}) = \\sigma^2 \\left( 1 + x_i^T (X_{(-i)}^T X_{(-i)})^{-1} x_i \\right) = \\frac{\\sigma^2}{1 &#8211; h_{ii}}$<\/span><\/p>\n<p>Dans ce cadre, la variance d&rsquo;erreur r\u00e9siduelle jackknife <span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2\">$s_{(-i)}^2$<\/span> est calcul\u00e9e \u00e0 partir de la somme des carr\u00e9s des erreurs r\u00e9siduelles du mod\u00e8le tronqu\u00e9, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"SSE_{(-i)}\">$SSE_{(-i)}$<\/span>, divis\u00e9e par ses degr\u00e9s de libert\u00e9 r\u00e9siduels <span class=\"arabpsy-math-inline\" data-latex=\"n - p - 1\">$n &#8211; p &#8211; 1$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2 = \\frac{SSE_{(-i)}}{n - p - 1}\">$s_{(-i)}^2 = \\frac{SSE_{(-i)}}{n &#8211; p &#8211; 1}$<\/span><\/p>\n<p>L&rsquo;application directe de la mise \u00e0 jour alg\u00e9brique \u00e9vite d&rsquo;avoir \u00e0 r\u00e9ajuster informatiquement <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> r\u00e9gressions distinctes. La relation ferm\u00e9e reliant <span class=\"arabpsy-math-inline\" data-latex=\"SSE_{(-i)}\">$SSE_{(-i)}$<\/span> \u00e0 la somme r\u00e9siduelle totale <span class=\"arabpsy-math-inline\" data-latex=\"SSE\">$SSE$<\/span> s&rsquo;\u00e9tablit de mani\u00e8re exacte par :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"SSE_{(-i)} = SSE - \\frac{e_i^2}{1 - h_{ii}}\">$SSE_{(-i)} = SSE &#8211; \\frac{e_i^2}{1 &#8211; h_{ii}}$<\/span><\/p>\n<p>En substituant cette \u00e9quivalence dans l&rsquo;estimation de la variance d&rsquo;erreur hors-\u00e9chantillon, nous obtenons la formulation de calcul instantan\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2 = \\frac{(n - p)s^2 - \\frac{e_i^2}{1 - h_{ii}}}{n - p - 1}\">$s_{(-i)}^2 = \\frac{(n &#8211; p)s^2 &#8211; \\frac{e_i^2}{1 &#8211; h_{ii}}}{n &#8211; p &#8211; 1}$<\/span><\/p>\n<p>Le r\u00e9sidu studentis\u00e9 externe <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span> se d\u00e9finit d\u00e8s lors comme le quotient du r\u00e9sidu d&rsquo;exclusion <span class=\"arabpsy-math-inline\" data-latex=\"d_i\">$d_i$<\/span> par son \u00e9cart-type estim\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"t_i = \\frac{d_i}{\\sqrt{\\widehat{\\text{Var}}(d_i)}} = \\frac{\\frac{e_i}{1 - h_{ii}}}{\\sqrt{\\frac{s_{(-i)}^2}{1 - h_{ii}}}} = \\frac{e_i}{s_{(-i)} \\sqrt{1 - h_{ii}}}\">$t_i = \\frac{d_i}{\\sqrt{\\widehat{\\text{Var}}(d_i)}} = \\frac{\\frac{e_i}{1 &#8211; h_{ii}}}{\\sqrt{\\frac{s_{(-i)}^2}{1 &#8211; h_{ii}}}} = \\frac{e_i}{s_{(-i)} \\sqrt{1 &#8211; h_{ii}}}$<\/span><\/p>\n<p>D\u00e9monstration de la loi de distribution : Par construction sous l&rsquo;hypoth\u00e8se nulle <span class=\"arabpsy-math-inline\" data-latex=\"H_0 : \\mathbb{E}[y_i] = x_i^T \\beta\">$H_0 : \\mathbb{E}[y_i] = x_i^T \\beta$<\/span>, le num\u00e9rateur transform\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"e_i \/ (\\sigma \\sqrt{1 - h_{ii}})\">$e_i \/ (\\sigma \\sqrt{1 &#8211; h_{ii}})$<\/span> suit une loi normale standard <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{N}(0, 1)\">$\\mathcal{N}(0, 1)$<\/span>. Parall\u00e8lement, la statistique d&rsquo;erreur r\u00e9siduelle jackknife normalis\u00e9e suit une distribution du khi-deux ind\u00e9pendante :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\frac{(n - p - 1)s_{(-i)}^2}{\\sigma^2} \\sim \\chi^2(n - p - 1)\">$\\frac{(n &#8211; p &#8211; 1)s_{(-i)}^2}{\\sigma^2} \\sim \\chi^2(n &#8211; p &#8211; 1)$<\/span><\/p>\n<p>Puisque l&rsquo;estimation <span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2\">$s_{(-i)}^2$<\/span> est calcul\u00e9e sans la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me observation, elle est math\u00e9matiquement ind\u00e9pendante de <span class=\"arabpsy-math-inline\" data-latex=\"y_i\">$y_i$<\/span>, et par cons\u00e9quent stochastiquement ind\u00e9pendante de <span class=\"arabpsy-math-inline\" data-latex=\"e_i\">$e_i$<\/span>. Le rapport <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span> correspond donc sans \u00e9quivoque au quotient d&rsquo;une variable gaussienne centr\u00e9e r\u00e9duite par la racine carr\u00e9e d&rsquo;une variable ind\u00e9pendante du khi-deux divis\u00e9e par ses degr\u00e9s de libert\u00e9. En vertu du th\u00e9or\u00e8me fondamental de Gosset, la statistique r\u00e9sultante converge exactement vers une distribution de Student :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"t_i \\sim \\mathcal{T}(n - p - 1)\">$t_i \\sim \\mathcal{T}(n &#8211; p &#8211; 1)$<\/span><\/p>\n<p>Cette propri\u00e9t\u00e9 distributionnelle exacte conf\u00e8re aux r\u00e9sidus studentis\u00e9s externes une puissance statistique nettement sup\u00e9rieure \u00e0 celle des <span class=\"arabpsy-math-inline\" data-latex=\"z\">$z$<\/span>-scores classiques bas\u00e9s sur une normalit\u00e9 asymptotique souvent viol\u00e9e dans les \u00e9chantillons de taille finie.<\/p>\n<p>## 4. Environnement Python pour l&rsquo;analyse diagnostique de r\u00e9gression<\/p>\n<p>### 4.1 \u00c9cosyst\u00e8me scientifique et d\u00e9pendances requises<\/p>\n<p>L&rsquo;\u00e9cosyst\u00e8me open-source de Python s&rsquo;est impos\u00e9 comme la r\u00e9f\u00e9rence contemporaine pour l&rsquo;analyse statistique avanc\u00e9e, l&rsquo;\u00e9conom\u00e9trie et la mod\u00e9lisation en sciences quantitatives. Pour mettre en \u0153uvre le calcul et l&rsquo;\u00e9valuation rigoureuse des r\u00e9sidus studentis\u00e9s, une synergie entre plusieurs biblioth\u00e8ques fondamentales est requise :<\/p>\n<ul>\n<li><strong>NumPy<\/strong> : Socle computationnel fournissant des structures de tableaux multidimensionnels \u00e0 hautes performances (ndarrays) et les primitives de bas niveau pour l&rsquo;alg\u00e8bre matricielle via LAPACK et BLAS.<\/li>\n<li><strong>Pandas<\/strong> : Outil de r\u00e9f\u00e9rence pour l&rsquo;ingestion, le nettoyage, la manipulation indic\u00e9e et la structuration tabulaire des jeux de donn\u00e9es psychom\u00e9triques ou exp\u00e9rimentaux.<\/li>\n<li><strong>Statsmodels<\/strong> : Biblioth\u00e8que fondamentale d&rsquo;\u00e9conom\u00e9trie et de statistique math\u00e9matique, d\u00e9velopp\u00e9e initialement pour pallier les orientations purement pr\u00e9dictives de Scikit-Learn. Statsmodels fournit des sorties d&rsquo;estimation exhaustives (coefficients, erreurs-types asymptotiques, tests diagnostiques, m\u00e9triques d&rsquo;influence) comparables \u00e0 celles de logiciels historiques comme R ou Stata.<\/li>\n<li><strong>SciPy<\/strong> : Sp\u00e9cifiquement son sous-module <code>scipy.stats<\/code>, essentiel pour \u00e9valuer de mani\u00e8re vectoris\u00e9e les fonctions de r\u00e9partition, les densit\u00e9s probabilistes et les fonctions de survie associ\u00e9es aux lois de Student et de Fisher.<\/li>\n<li><strong>Matplotlib et Seaborn<\/strong> : Biblioth\u00e8ques graphiques permettant de g\u00e9n\u00e9rer des planches visuelles diagnostiques conformes aux standards de publication acad\u00e9mique.<\/li>\n<\/ul>\n<p>Pour garantir la p\u00e9rennit\u00e9 et la reproductibilit\u00e9 computationnelle de vos protocoles diagnostiques, l&rsquo;isolation des d\u00e9pendances logicielles au sein d&rsquo;environnements virtuels (via `venv`, `conda` ou `poetry`) demeure une exigence m\u00e9thodologique stricte. Les variations mineures de version dans l&rsquo;impl\u00e9mentation de LAPACK ou de Statsmodels pouvant impacter la pr\u00e9cision en virgule flottante des inversions matricielles, le figeage des versions via un fichier de configuration assure la stricte reproductibilit\u00e9 des inf\u00e9rences.<\/p>\n<p>### 4.2 Importation des biblioth\u00e8ques et structure de travail<\/p>\n<p>L&rsquo;initialisation d&rsquo;un pipeline diagnostique requiert le chargement des alias conventionnels garantissant la lisibilit\u00e9 et l&rsquo;interop\u00e9rabilit\u00e9 du code source. Deux points d&rsquo;acc\u00e8s principaux existent dans Statsmodels : l&rsquo;API classique (`statsmodels.api`), qui attend une matrice de r\u00e9gresseurs explicite requ\u00e9rant l&rsquo;adjonction manuelle d&rsquo;une constante, et l&rsquo;API de formules (`statsmodels.formula.api`), qui \u00e9mule la syntaxe d\u00e9clarative de la biblioth\u00e8que `Patsy` inspir\u00e9e du syst\u00e8me de formules de R.<\/p>\n<p>Pour \u00e9viter la troncature pr\u00e9judiciable des sorties diagnostiques tabulaires lors du traitement de cohortes comptant plusieurs centaines de participants, il convient d&rsquo;ajuster les options globales de formatage du moteur d&rsquo;affichage de Pandas. L&rsquo;analyste configure notamment la pr\u00e9cision d&rsquo;affichage des nombres r\u00e9els, le nombre maximal de lignes affich\u00e9es et la largeur d&rsquo;\u00e9cran. <\/p>\n<p>Enfin, la reproductibilit\u00e9 parfaite des simulations computationnelles et des v\u00e9rifications empiriques exige l&rsquo;initialisation explicite des g\u00e9n\u00e9rateurs de nombres pseudo-al\u00e9atoires. L&rsquo;utilisation conjointe des g\u00e9n\u00e9rateurs de NumPy (`numpy.random.seed` ou de la classe moderne `default_rng`) garantit que chaque \u00e9tape d&rsquo;injection de perturbations stochastiques ou de valeurs aberrantes demeure rigoureusement stable d&rsquo;une ex\u00e9cution \u00e0 l&rsquo;autre.<\/p>\n<p>## 5. Ajustement d&rsquo;un mod\u00e8le de r\u00e9gression lin\u00e9aire avec statsmodels<\/p>\n<p>### 5.1 Cr\u00e9ation et structuration d&rsquo;un jeu de donn\u00e9es psychom\u00e9triques<\/p>\n<p>Pour contextualiser l&rsquo;application pratique du diagnostic de r\u00e9gression, nous concevons un environnement de simulation contr\u00f4l\u00e9 refl\u00e9tant les dynamiques observ\u00e9es en psychologie de l&rsquo;\u00e9ducation et en psychom\u00e9trie clinique. Consid\u00e9rons un protocole \u00e9valuant la performance acad\u00e9mique finale (not\u00e9e sur une \u00e9chelle continue de 0 \u00e0 100) en fonction de trois pr\u00e9dicteurs comportementaux : l&rsquo;anxi\u00e9t\u00e9 cognitive d&rsquo;\u00e9valuation (mesur\u00e9e par une \u00e9chelle standardis\u00e9e de type Spielberger), le volume hebdomadaire d&rsquo;heures d&rsquo;\u00e9tude structur\u00e9e, et un score de conscience professionnelle extrait du mod\u00e8le \u00e0 cinq facteurs (Big Five).<\/p>\n<p>Dans un \u00e9chantillon repr\u00e9sentatif de <span class=\"arabpsy-math-inline\" data-latex=\"n = 120\">$n = 120$<\/span> \u00e9tudiants, la relation sous-jacente suit un processus lin\u00e9aire bruit\u00e9 par une composante gaussienne classique. Afin de tester la sensibilit\u00e9 des r\u00e9sidus studentis\u00e9s, nous introduisons d\u00e9lib\u00e9r\u00e9ment dans ce jeu de donn\u00e9es des observations aberrantes synth\u00e9tiques, imitant fid\u00e8lement les art\u00e9facts m\u00e9thodologiques r\u00e9currents en milieu universitaire : un individu pr\u00e9sentant un niveau d&rsquo;anxi\u00e9t\u00e9 normal mais un \u00e9chec catastrophique li\u00e9 \u00e0 une panne mat\u00e9rielle durant la passation informatis\u00e9e, ainsi qu&rsquo;un profil affichant un engagement d&rsquo;\u00e9tude marginal coupl\u00e9 \u00e0 une note anormalement brillante imputable \u00e0 un comportement de tricherie ou \u00e0 un effet de plafond intellectuel non mesur\u00e9.<\/p>\n<p>Les donn\u00e9es simul\u00e9es sont consolid\u00e9es dans un DataFrame Pandas. Une inspection pr\u00e9liminaire syst\u00e9matique permet de v\u00e9rifier l&rsquo;absence de valeurs manquantes fortuites, de confirmer le typage num\u00e9rique rigoureux des variables continues et d&rsquo;analyser les statistiques descriptives initiales (moyenne, \u00e9cart-type, \u00e9tendue, asym\u00e9trie et kurtosis) avant toute soumission aux algorithmes d&rsquo;estimation matricielle.<\/p>\n<p>### 5.2 Sp\u00e9cification et entra\u00eenement du mod\u00e8le OLS<\/p>\n<p>L&rsquo;ajustement du mod\u00e8le de r\u00e9gression par les moindres carr\u00e9s ordinaires s&rsquo;effectue au moyen de la m\u00e9thode `ols()` issue du module `statsmodels.formula.api`. Cette approche permet de d\u00e9clarer la structure structurelle du mod\u00e8le par une cha\u00eene syntaxique claire :<\/p>\n<p>`performance ~ anxiete + heures_etude + conscience`<\/p>\n<p>L&rsquo;ex\u00e9cution de la m\u00e9thode `.fit()` d\u00e9clenche l&rsquo;estimation des moindres carr\u00e9s ordinaires en proc\u00e9dant \u00e0 la r\u00e9solution du syst\u00e8me d&rsquo;\u00e9quations normales par d\u00e9composition QR ou d\u00e9composition en valeurs singuli\u00e8res (SVD) de la matrice de conception. Cette op\u00e9ration instancie un objet complexe de classe `RegressionResults` (ou plus sp\u00e9cifiquement `OLSResults`), qui encapsulate l&rsquo;int\u00e9gralit\u00e9 des statistiques d&rsquo;ajustement, des param\u00e8tres vectoriels estim\u00e9s, de la matrice de covariance variance-covariance des coefficients et des structures r\u00e9siduelles n\u00e9cessaires au diagnostic.<\/p>\n<p>L&rsquo;interrogation synth\u00e9tique du mod\u00e8le s&rsquo;op\u00e8re usuellement par l&rsquo;appel de la m\u00e9thode `.summary()`, qui d\u00e9livre un rapport exhaustif divis\u00e9 en trois volets analytiques :<\/p>\n<ol>\n<li>Le volet sup\u00e9rieur d\u00e9taillant les caract\u00e9ristiques de l&rsquo;ajustement global : coefficient de d\u00e9termination <span class=\"arabpsy-math-inline\" data-latex=\"R^2\">$R^2$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"R^2\">$R^2$<\/span> ajust\u00e9, statistique globale de Fisher (F-statistic) testant l&rsquo;hypoth\u00e8se conjointe de nullit\u00e9 des coefficients non constants, et crit\u00e8res d&rsquo;information d&rsquo;Akaike (AIC) et bay\u00e9sien (BIC).<\/li>\n<li>Le volet interm\u00e9diaire exposant les estimations ponctuelles des coefficients <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}_j\">$\\hat{\\beta}_j$<\/span>, leurs erreurs-types standard corrig\u00e9es, les statistiques <span class=\"arabpsy-math-inline\" data-latex=\"t\">$t$<\/span> de Wald associ\u00e9es et les valeurs <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> bilat\u00e9rales permettant de trancher la significativit\u00e9 marginale de chaque covariable.<\/li>\n<li>Le volet inf\u00e9rieur pr\u00e9sentant les premiers indicateurs de validit\u00e9 r\u00e9siduelle : test d&rsquo;asym\u00e9trie et d&rsquo;aplatissement d&rsquo;Omnibus, statistique d&rsquo;autocorr\u00e9lation de Durbin-Watson, test de Jarque-Bera et conditionnement num\u00e9rique de la matrice de conception (Condition Number), alertant sur la pr\u00e9sence de multicolin\u00e9arit\u00e9 s\u00e9v\u00e8re.<\/li>\n<\/ol>\n<p>Ces m\u00e9triques globales, bien qu&rsquo;essentielles, demeurent des agr\u00e9gats d&rsquo;ensemble incapables d&rsquo;isoler les points pathologiques au niveau individuel. C&rsquo;est pr\u00e9cis\u00e9ment \u00e0 cette jonction m\u00e9thodologique que s&rsquo;impose le d\u00e9ploiement des r\u00e9sidus studentis\u00e9s.<\/p>\n<p>## 6. Calcul des r\u00e9sidus studentis\u00e9s via la m\u00e9thode outlier_test()<\/p>\n<p>### 6.1 Syntaxe op\u00e9rationnelle et m\u00e9canisme interne de OLSResults.outlier_test()<\/p>\n<p>L&rsquo;acc\u00e8s le plus direct, robuste et optimis\u00e9 pour le calcul des r\u00e9sidus studentis\u00e9s externes au sein de Statsmodels r\u00e9side dans la m\u00e9thode native `outlier_test()`, accessible directement depuis l&rsquo;instance r\u00e9sultante de l&rsquo;ajustement `OLSResults`. Cette routine ex\u00e9cute de mani\u00e8re enti\u00e8rement vectoris\u00e9e l&rsquo;algorithme d&rsquo;exclusion pr\u00e9sent\u00e9 \u00e0 la section 3.2, dispensant l&rsquo;utilisateur de l&rsquo;\u00e9criture fastidieuse et inefficace d&rsquo;une boucle it\u00e9rative proc\u00e9dant \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> r\u00e9estimations individuelles du mod\u00e8le.<\/p>\n<p>L&rsquo;invocation basique de la fonction s&rsquo;ex\u00e9cute selon l&rsquo;instruction :<\/p>\n<p>`resultats_outliers = modele_ajuste.outlier_test(method=&rsquo;bonferroni&rsquo;, alpha=0.05)`<\/p>\n<p>La m\u00e9thode `outlier_test()` accepte plusieurs param\u00e8tres fondamentaux qui pilotent le traitement inf\u00e9rentiel :<\/p>\n<ul>\n<li><code>method<\/code> : Sp\u00e9cifie la technique de correction pour comparaisons multiples appliqu\u00e9e aux probabilit\u00e9s critiques. Les arguments recevables incluent <code>'bonferroni'<\/code> (ajustement conservateur par d\u00e9faut), <code>'fdr_bh'<\/code> (contr\u00f4le du taux de fausses d\u00e9couvertes selon Benjamini-Hochberg), <code>'sidak'<\/code> ou <code>'none'<\/code> pour obtenir les probabilit\u00e9s brutes non ajust\u00e9es.<\/li>\n<li><code>alpha<\/code> : Fixe le seuil de tol\u00e9rance nominal pour le taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce global appliqu\u00e9 au test d&rsquo;atypicit\u00e9 conjointe.<\/li>\n<li><code>cutoff<\/code> : Seuil scalaire optionnel permettant de filtrer automatiquement les sorties du tableau pour ne retourner que les observations dont la valeur <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> ajust\u00e9e est strictement inf\u00e9rieure \u00e0 la borne sp\u00e9cifi\u00e9e. Par d\u00e9faut, la m\u00e9thode retourne la table compl\u00e8te des <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> observations.<\/li>\n<\/ul>\n<p>Le fonctionnement sous-jacent de cette m\u00e9thode mobilise les \u00e9quations de mise \u00e0 jour alg\u00e9brique pour d\u00e9duire instantan\u00e9ment les variances r\u00e9siduelles d&rsquo;exclusion <span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2\">$s_{(-i)}^2$<\/span> \u00e0 partir des r\u00e9sidus ordinaires et de la diagonale de projection <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span>. La routine convertit ensuite ces r\u00e9sidus externes en statistiques de Student pour interroger la fonction de survie de la distribution <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{T}(n - p - 1)\">$\\mathcal{T}(n &#8211; p &#8211; 1)$<\/span>.<\/p>\n<p>### 6.2 D\u00e9composition des composantes du tableau de sortie<\/p>\n<p>L&rsquo;objet retourn\u00e9 par l&rsquo;appel de `outlier_test()` est un DataFrame Pandas index\u00e9 scrupuleusement sur l&rsquo;index d&rsquo;origine du jeu de donn\u00e9es ayant servi \u00e0 l&rsquo;estimation. Ce tableau structure trois colonnes fondamentales pour l&rsquo;investigation empirique :<\/p>\n<ol>\n<li><code>student_resid<\/code> : Cette colonne stocke la valeur num\u00e9rique exacte du r\u00e9sidu studentis\u00e9 externe <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>. Elle repr\u00e9sente l&rsquo;ampleur standardis\u00e9e de la d\u00e9viation de l&rsquo;observation mesur\u00e9e en unit\u00e9s d&rsquo;\u00e9cart-type jackknife. Une valeur positive indique une performance observ\u00e9e sup\u00e9rieure aux pr\u00e9dictions conditionnelles du mod\u00e8le, tandis qu&rsquo;une valeur n\u00e9gative d\u00e9note une sous-performance majeure. En raison de sa conformit\u00e9 \u00e0 une loi de Student \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"n - p - 1\">$n &#8211; p &#8211; 1$<\/span> degr\u00e9s de libert\u00e9, cette m\u00e9trique est imm\u00e9diatement assimilable \u00e0 un score de d\u00e9viation critique standardis\u00e9.<\/li>\n<li><code>unadj_p<\/code> : Cette colonne contient la valeur <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> brute non ajust\u00e9e issue du test d&rsquo;hypoth\u00e8se bilat\u00e9ral <span class=\"arabpsy-math-inline\" data-latex=\"H_0 : \\mathbb{E}[\\varepsilon_i] = 0\">$H_0 : \\mathbb{E}[\\varepsilon_i] = 0$<\/span> contre <span class=\"arabpsy-math-inline\" data-latex=\"H_1 : \\mathbb{E}[\\varepsilon_i] \\neq 0\">$H_1 : \\mathbb{E}[\\varepsilon_i] \\neq 0$<\/span>. Elle est calcul\u00e9e formellement par l&rsquo;expression <span class=\"arabpsy-math-inline\" data-latex=\"2 \\times (1 - F_{\\mathcal{T}}(|t_i|))\">$2 \\times (1 &#8211; F_{\\mathcal{T}}(|t_i|))$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"F_{\\mathcal{T}}\">$F_{\\mathcal{T}}$<\/span> d\u00e9signe la fonction de r\u00e9partition d&rsquo;une variable al\u00e9atoire suivant la loi de Student \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"n - p - 1\">$n &#8211; p &#8211; 1$<\/span> degr\u00e9s de libert\u00e9. Cette probabilit\u00e9 quantifie le risque d&rsquo;observer une d\u00e9viation au moins aussi extr\u00eame que <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span> sous la conjecture stricte que le point <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> est issu du mod\u00e8le nominal.<\/li>\n<li><code>bonf(p)<\/code> : Cette colonne d\u00e9livre la probabilit\u00e9 critique corrig\u00e9e pour tenir compte du risque d&rsquo;inflation d&rsquo;erreur globale li\u00e9 \u00e0 l&rsquo;\u00e9valuation simultan\u00e9e de l&rsquo;ensemble des points de l&rsquo;\u00e9chantillon. Par d\u00e9faut, cette valeur correspond \u00e0 l&rsquo;application rigoureuse de la borne de Bonferroni, calcul\u00e9e par la multiplication <span class=\"arabpsy-math-inline\" data-latex=\"\\min(1, n \\times \\text{unadj_p})\">$\\min(1, n \\times \\text{unadj_p})$<\/span>. Une valeur <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> corrig\u00e9e inf\u00e9rieure au seuil de significativit\u00e9 nominal <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,05\">$\\alpha = 0,05$<\/span> atteste que l&rsquo;observation s&rsquo;\u00e9carte de la surface de r\u00e9gression avec une intensit\u00e9 statistiquement incompatible avec les simples fluctuations d&rsquo;\u00e9chantillonnage de l&rsquo;\u00e9chantillon global.<\/li>\n<\/ol>\n<p>La conservation int\u00e9grale de l&rsquo;index du DataFrame source garantit un alignement d\u00e9terministe imm\u00e9diat avec les m\u00e9tadonn\u00e9es de l&rsquo;\u00e9tude (identifiants anonymis\u00e9s des participants, variables descriptives compl\u00e9mentaires ou centres d&rsquo;acquisition clinique), autorisant un audit cibl\u00e9 imm\u00e9diat des cas isol\u00e9s.<\/p>\n<p>### 6.3 Acc\u00e8s alternatif via la classe OLSInfluence<\/p>\n<p>Bien que la fonction `outlier_test()` constitue la voie la plus rapide pour le d\u00e9pistage d&rsquo;observations discordantes, Statsmodels propose une architecture objet plus approfondie via la classe `OLSInfluence`. Cette interface s&rsquo;instancie par la m\u00e9thode `.get_influence()` appel\u00e9e sur les r\u00e9sultats de la r\u00e9gression :<\/p>\n<p>`influence_diagnostics = modele_ajuste.get_influence()`<\/p>\n<p>L&rsquo;instance d&rsquo;objet ainsi g\u00e9n\u00e9r\u00e9e offre un acc\u00e8s unifi\u00e9 \u00e0 l&rsquo;int\u00e9gralit\u00e9 du spectre diagnostique de l&rsquo;\u00e9cole de r\u00e9gression classique de Belsley, Kuh, Welsch et Cook. Elle expose notamment deux propri\u00e9t\u00e9s sp\u00e9cifiques relatives aux r\u00e9sidus normalis\u00e9s :<\/p>\n<ul>\n<li><code>influence_diagnostics.resid_studentized_internal<\/code> : Retourne sous la forme d&rsquo;un tableau NumPy unidimensionnel les r\u00e9sidus studentis\u00e9s internes <span class=\"arabpsy-math-inline\" data-latex=\"r_i\">$r_i$<\/span>. Ces valeurs s&rsquo;appuient sur l&rsquo;estimateur sans biais global <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>, repr\u00e9sentant le vecteur des r\u00e9sidus standardis\u00e9s sans exclusion de l&rsquo;observation.<\/li>\n<li><code>influence_diagnostics.resid_studentized_external<\/code> : Retourne le tableau vectoris\u00e9 des r\u00e9sidus studentis\u00e9s externes <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>, \u00e9quivalant strictement aux donn\u00e9es num\u00e9riques de la colonne <code>student_resid<\/code> produite par <code>outlier_test()<\/code>.<\/li>\n<\/ul>\n<p>Le tableau comparatif suivant synth\u00e9tise les propri\u00e9t\u00e9s op\u00e9rationnelles de ces deux approches d&rsquo;acc\u00e8s programmatique :<\/p>\n<p>| Propri\u00e9t\u00e9 \/ Crit\u00e8re | M\u00e9thode `outlier_test()` | Attributs de `get_influence()` |<br \/>\n| :&#8212; | :&#8212; | :&#8212; |<br \/>\n| **Type de r\u00e9sidu principal** | R\u00e9sidu studentis\u00e9 externe (<span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>) | Internes (<span class=\"arabpsy-math-inline\" data-latex=\"r_i\">$r_i$<\/span>) et Externes (<span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>) |<br \/>\n| **Calcul des valeurs <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>** | Automatique (brutes et corrig\u00e9es) | Non inclus directement |<br \/>\n| **Correction pour tests multiples**| Bonferroni, FDR de Benjamini-Hochberg | Non applicable sans calcul tiers |<br \/>\n| **Format de retour** | DataFrame Pandas structur\u00e9 avec index | Tableaux NumPy ou s\u00e9ries Pandas |<br \/>\n| **M\u00e9triques suppl\u00e9mentaires** | Aucune (focalisation sur l&rsquo;atypicit\u00e9) | Distance de Cook, DFBETAS, DFFITS, Levier |<br \/>\n| **Empreinte computationnelle** | Faible \u00e0 mod\u00e9r\u00e9e | Plus \u00e9lev\u00e9e (calcul de tout le spectre) |<\/p>\n<p>L&rsquo;exploitation conjointe de ces structures programmatiques s&rsquo;av\u00e8re particuli\u00e8rement pertinente : `outlier_test()` permet d&rsquo;ex\u00e9cuter le test de d\u00e9tection formel, tandis que `get_influence()` apporte les m\u00e9triques compl\u00e9mentaires pour contextualiser si l&rsquo;observation atypique identifi\u00e9e d\u00e9stabilise significativement l&rsquo;architecture param\u00e9trique globale.<\/p>\n<p>## 7. Seuils critiques, tests d&rsquo;hypoth\u00e8se et ajustement pour comparaisons multiples<\/p>\n<p>### 7.1 La r\u00e8gle heuristique absolue du seuil de 3<\/p>\n<p>Dans la pratique de la mod\u00e9lisation statistique appliqu\u00e9e, de nombreux analystes recourent \u00e0 des heuristiques op\u00e9rationnelles pour isoler rapidement les observations probl\u00e9matiques. La r\u00e8gle empirique la plus largement diss\u00e9min\u00e9e postule qu&rsquo;un r\u00e9sidu studentis\u00e9 dont la valeur absolue exc\u00e8de un seuil critique arbitraire fix\u00e9 \u00e0 2 ou 3 doit \u00eatre \u00e9tiquet\u00e9 comme suspect ou aberrant :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"|t_i| &gt; 2 \\quad \\text{(zone d&#039;avertissement)} \\qquad \\text{ou} \\qquad |t_i| &gt; 3 \\quad \\text{(valeur aberrante suspecte)}\">$|t_i| &gt; 2 \\quad \\text{(zone d&rsquo;avertissement)} \\qquad \\text{ou} \\qquad |t_i| &gt; 3 \\quad \\text{(valeur aberrante suspecte)}$<\/span><\/p>\n<p>Ce crit\u00e8re s&rsquo;enracine th\u00e9oriquement dans la r\u00e8gle empirique de la loi normale standard (la r\u00e8gle des <span class=\"arabpsy-math-inline\" data-latex=\"3\\sigma\">$3\\sigma$<\/span>). Dans une distribution gaussienne univari\u00e9e id\u00e9ale, la probabilit\u00e9 qu&rsquo;une observation s&rsquo;\u00e9carte de son esp\u00e9rance de plus de trois \u00e9carts-types est d&rsquo;approximativement <span class=\"arabpsy-math-inline\" data-latex=\"P(|Z| &gt; 3) \\approx 0,0027\">$P(|Z| &gt; 3) \\approx 0,0027$<\/span>, soit un peu plus de deux observations pour mille. Lorsque les degr\u00e9s de libert\u00e9 r\u00e9siduels <span class=\"arabpsy-math-inline\" data-latex=\"n - p - 1\">$n &#8211; p &#8211; 1$<\/span> sont mod\u00e9r\u00e9s \u00e0 \u00e9lev\u00e9s (sup\u00e9rieurs \u00e0 30 ou 40), la loi de Student <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{T}(n - p - 1)\">$\\mathcal{T}(n &#8211; p &#8211; 1)$<\/span> pr\u00e9sente une allure morphologique tr\u00e8s proche de la gaussienne, justifiant intuitivement la transposition de ce cut-off de 3.<\/p>\n<p>Toutefois, l&rsquo;application aveugle et acritique de cette r\u00e8gle d\u00e9terministe pr\u00e9sente des failles m\u00e9thodologiques majeures conditionn\u00e9es par la taille totale de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span>. Consid\u00e9rons un \u00e9chantillon de psychom\u00e9trie comportementale comprenant <span class=\"arabpsy-math-inline\" data-latex=\"n = 1,000\">$n = 1,000$<\/span> participants. M\u00eame dans l&rsquo;hypoth\u00e8se o\u00f9 le mod\u00e8le lin\u00e9aire serait parfaitement sp\u00e9cifi\u00e9 et o\u00f9 l&rsquo;ensemble des perturbations stochastiques d\u00e9riverait d&rsquo;une distribution normale pure et irr\u00e9prochable, l&rsquo;esp\u00e9rance math\u00e9matique du nombre d&rsquo;observations d\u00e9passant le seuil critique <span class=\"arabpsy-math-inline\" data-latex=\"|t_i| &gt; 3\">$|t_i| &gt; 3$<\/span> par pure fluctuation al\u00e9atoire n&rsquo;est pas nulle, mais s&rsquo;\u00e9l\u00e8ve \u00e0 :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{E}[N_{|t_i| &gt; 3}] = n \\times P(|t_i| &gt; 3) \\approx 1,000 \\times 0,0027 = 2,7 \\text{ observations}\">$\\mathbb{E}[N_{|t_i| &gt; 3}] = n \\times P(|t_i| &gt; 3) \\approx 1,000 \\times 0,0027 = 2,7 \\text{ observations}$<\/span><\/p>\n<p>Dans une cohorte de grande \u00e9chelle (<span class=\"arabpsy-math-inline\" data-latex=\"n \\geq 10,000\">$n \\geq 10,000$<\/span>), plusieurs dizaines d&rsquo;individus parfaitement r\u00e9guliers afficheront m\u00e9caniquement des r\u00e9sidus studentis\u00e9s externes exc\u00e9dant 3 sans refl\u00e9ter la moindre anomalie m\u00e9thodologique. \u00c0 l&rsquo;inverse, dans un protocole clinique ou neuropsychologique restreint typique des \u00e9tudes en imagerie c\u00e9r\u00e9brale fonctionnelle ou sur des pathologies rares (<span class=\"arabpsy-math-inline\" data-latex=\"n = 15\">$n = 15$<\/span> \u00e0 $25$), un seuil d&rsquo;alerte de 3 s&rsquo;av\u00e8re excessivement conservateur. Dans un mod\u00e8le \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"n = 20\">$n = 20$<\/span> avec <span class=\"arabpsy-math-inline\" data-latex=\"p = 4\">$p = 4$<\/span>, la valeur critique bilat\u00e9rale exacte d&rsquo;une loi de Student \u00e0 15 degr\u00e9s de libert\u00e9 pour un seuil nominal <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,05\">$\\alpha = 0,05$<\/span> s&rsquo;\u00e9tablit \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"t_{0,025; 15} \\approx 2,131\">$t_{0,025; 15} \\approx 2,131$<\/span>. Fixer un seuil fixe \u00e0 3 conduirait \u00e0 n\u00e9gliger des d\u00e9viations majeures susceptibles de bouleverser enti\u00e8rement l&rsquo;estimation des pentes du mod\u00e8le.<\/p>\n<p>### 7.2 Ajustement de Bonferroni et contr\u00f4le du taux d&rsquo;erreur par famille<\/p>\n<p>Pour conf\u00e9rer une rigueur probabiliste authentique \u00e0 l&rsquo;\u00e9valuation diagnostique, la recherche d&rsquo;observations aberrantes doit \u00eatre formellement mod\u00e9lis\u00e9e comme un probl\u00e8me de tests d&rsquo;hypoth\u00e8ses multiples simultan\u00e9s. Lorsque le chercheur examine les r\u00e9sidus studentis\u00e9s de l&rsquo;ensemble de son \u00e9chantillon, il ne conduit pas un test d&rsquo;hypoth\u00e8se isol\u00e9 sur un participant choisi a priori. Il teste conjointement <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> hypoth\u00e8ses nulles distinctes :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"H_{0, i} : \\mathbb{E}[\\varepsilon_i] = 0 \\quad \\text{pour tout } i in {1, 2, dots, n}\">$H_{0, i} : \\mathbb{E}[\\varepsilon_i] = 0 \\quad \\text{pour tout } i in {1, 2, dots, n}$<\/span><\/p>\n<p>Si l&rsquo;analyste applique le seuil critique standard <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,05\">$\\alpha = 0,05$<\/span> ind\u00e9pendamment \u00e0 chaque observation, le taux d&rsquo;erreur par famille (FWER, pour *Family-Wise Error Rate*), c&rsquo;est-\u00e0-dire la probabilit\u00e9 d&rsquo;identifier \u00e0 tort au moins une observation comme aberrante par simple fluctuation stochastique, subit une inflation drastique r\u00e9gie par l&rsquo;in\u00e9galit\u00e9 de Boole-Fr\u00e9chet :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\text{FWER} = P\\left( \\big\\cup_{i=1}^n { \\text{Rejet de } H_{0, i} mid H_{0, i} \\text{ est vraie} } \\right) \\leq 1 - (1 - \\alpha)^n\">$\\text{FWER} = P\\left( \\big\\cup_{i=1}^n { \\text{Rejet de } H_{0, i} mid H_{0, i} \\text{ est vraie} } \\right) \\leq 1 &#8211; (1 &#8211; \\alpha)^n$<\/span><\/p>\n<p>Pour un \u00e9chantillon modeste de <span class=\"arabpsy-math-inline\" data-latex=\"n = 100\">$n = 100$<\/span> individus avec <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,05\">$\\alpha = 0,05$<\/span>, ce risque global th\u00e9orique atteint <span class=\"arabpsy-math-inline\" data-latex=\"1 - (1 - 0,05)^{100} \\approx 0,994\">$1 &#8211; (1 &#8211; 0,05)^{100} \\approx 0,994$<\/span>. Il est donc virtuellement certain (plus de 99 % de chances) que le chercheur identifiera \u00e0 tort au moins un pr\u00e9tendu \u00ab\u00a0outlier\u00a0\u00bb au sein d&rsquo;un \u00e9chantillon parfait s&rsquo;il omet d&rsquo;ajuster son seuil d\u00e9cisionnel.<\/p>\n<p>La proc\u00e9dure impl\u00e9ment\u00e9e par d\u00e9faut dans `OLSResults.outlier_test()` r\u00e9pond \u00e0 cet imp\u00e9ratif en appliquant la correction de Bonferroni. Cette m\u00e9thode garantit le contr\u00f4le strict du FWER \u00e0 un seuil global <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha_{\\text{global}}\">$\\alpha_{\\text{global}}$<\/span> (typiquement 0,05) en exigeant que chaque test individuel soit \u00e9valu\u00e9 \u00e0 un niveau de significativit\u00e9 nominal drastiquement r\u00e9duit :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\alpha_{\\text{local}} = \\frac{\\alpha_{\\text{global}}}{n}\">$\\alpha_{\\text{local}} = \\frac{\\alpha_{\\text{global}}}{n}$<\/span><\/p>\n<p>Sous ce prisme, une observation n&rsquo;est d\u00e9clar\u00e9e statistiquement discordante de mani\u00e8re confirmatoire que si sa probabilit\u00e9 critique bilat\u00e9rale brute v\u00e9rifie <span class=\"arabpsy-math-inline\" data-latex=\"\\text{unadj_p} &lt; \\frac{0,05}{n}\">$\\text{unadj_p} &lt; \\frac{0,05}{n}$<\/span>, ce qui \u00e9quivaut rigoureusement \u00e0 constater que la colonne `bonf(p)` de la sortie Statsmodels affiche une valeur strictement inf\u00e9rieure \u00e0 0,05. <\/p>\n<p>Bien que la correction de Bonferroni assure une protection robuste contre les fausses d\u00e9couvertes, elle est notoirement conservatrice lorsque la taille d&#039;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> cro\u00eet substantiellement. Dans les jeux de donn\u00e9es volumineux, l&#039;analyste peut avantageusement substituer au contr\u00f4le du FWER le contr\u00f4le du taux de fausses d\u00e9couvertes (FDR, *False Discovery Rate*) th\u00e9oris\u00e9 par <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2517-6161.1995.tb02031.x\">Benjamini et Hochberg (1995)<\/a>. En configurant l&rsquo;argument `method=&rsquo;fdr_bh&rsquo;`, l&rsquo;algorithme ordonne les valeurs <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> croissantes et optimise le compromis entre la d\u00e9tection des v\u00e9ritables anomalies et la pr\u00e9servation de la puissance statistique globale.<\/p>\n<p>### 7.3 Prise de d\u00e9cision empirique face \u00e0 une valeur p significative<\/p>\n<p>L&rsquo;obtention d&rsquo;une valeur <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> corrig\u00e9e statistiquement significative (`bonf(p) &lt; 0.05`) constitue un signal d&#039;alerte algorithmique, mais ne doit en aucun cas \u00eatre interpr\u00e9t\u00e9e comme un mandat automatique prescrivant la suppression imm\u00e9diate et m\u00e9canique de l&#039;observation concern\u00e9e. Une telle pratique de purge automatique rel\u00e8ve d&#039;une faute m\u00e9thodologique grave qui vicie l&#039;int\u00e9grit\u00e9 de la d\u00e9marche scientifique et s&#039;apparente \u00e0 une forme insidieuse de dragage de donn\u00e9es (*data dredging* ou *p-hacking*).<\/p>\n<p>La confrontation \u00e0 une observation affichant un r\u00e9sidu studentis\u00e9 significatif impose une d\u00e9marche d&#039;investigation s\u00e9quentielle articul\u00e9e autour de trois axes cardinaux :<\/p>\n<ol>\n<li><strong>V\u00e9rification de l&rsquo;int\u00e9grit\u00e9 factuelle<\/strong> : L&rsquo;investigateur doit imm\u00e9diatement inspecter la cha\u00eene de traitement num\u00e9rique pour traquer d&rsquo;\u00e9ventuelles erreurs d&rsquo;encodage (par exemple, la saisie d&rsquo;un score de 999 codant initialement une donn\u00e9e manquante comme une valeur num\u00e9rique r\u00e9elle), des d\u00e9faillances mat\u00e9rielles de capteurs psychophysiologiques ou des violations explicites des crit\u00e8res d&rsquo;inclusion du protocole clinique.<\/li>\n<li><strong>\u00c9valuation de la pertinence \u00e9cologique et clinique<\/strong> : Si la mesure est factuellement authentique, l&rsquo;analyste doit \u00e9valuer si le point refl\u00e8te une sous-population non anticip\u00e9e mais \u00e9cologiquement valide. Dans une r\u00e9gression mod\u00e9lisant la r\u00e9ponse th\u00e9rapeutique \u00e0 une psychoth\u00e9rapie, un patient d\u00e9montrant une d\u00e9t\u00e9rioration s\u00e9v\u00e8re et inattendue n&rsquo;est pas un \u00ab\u00a0bruit\u00a0\u00bb statistique \u00e0 \u00e9liminer pour am\u00e9liorer le <span class=\"arabpsy-math-inline\" data-latex=\"R^2\">$R^2$<\/span>, mais un profil de non-r\u00e9pondeur d&rsquo;un int\u00e9r\u00eat clinique fondamental dont l&rsquo;exclusion masquerait les limites r\u00e9elles de l&rsquo;intervention test\u00e9e.<\/li>\n<li><strong>\u00c9valuation de la sp\u00e9cification structurelle du mod\u00e8le<\/strong> : Un r\u00e9sidu studentis\u00e9 d&rsquo;amplitude anormale peut r\u00e9v\u00e9ler non pas une anomalie inh\u00e9rente \u00e0 l&rsquo;individu, mais une mauvaise sp\u00e9cification fonctionnelle du mod\u00e8le lin\u00e9aire. L&rsquo;omission d&rsquo;un terme polynomial d&rsquo;ordre sup\u00e9rieur (ph\u00e9nom\u00e8ne non lin\u00e9aire quadratique), l&rsquo;absence d&rsquo;un terme d&rsquo;interaction critique entre deux variables, ou le recours inappropri\u00e9 \u00e0 un lien lin\u00e9aire pour des donn\u00e9es born\u00e9es ou asym\u00e9triques suffisent \u00e0 g\u00e9n\u00e9rer des r\u00e9sidus studentis\u00e9s externes hautement significatifs aux marges de la distribution.<\/li>\n<\/ol>\n<p>Toute d\u00e9cision finale \u2014 qu&rsquo;elle consiste \u00e0 conserver l&rsquo;observation, \u00e0 proc\u00e9der \u00e0 une mod\u00e9lisation alternative robuste ou, dans des cas strictement circonscrits et document\u00e9s, \u00e0 \u00e9carter l&rsquo;observation de l&rsquo;analyse primaire \u2014 doit \u00eatre consign\u00e9e de mani\u00e8re exhaustive et transparente dans les rapports de recherche finaux.<\/p>\n<p>## 8. Impl\u00e9mentation pas \u00e0 pas et vectoris\u00e9e avec NumPy et Pandas<\/p>\n<p>### 8.1 Calcul vectoriel des valeurs diagonales de la matrice chapeau<\/p>\n<p>Pour appr\u00e9hender intimement la m\u00e9canique algorithmique des diagnostics de r\u00e9gression et acqu\u00e9rir une autonomie totale ind\u00e9pendamment des modules de haut niveau, il est hautement formateur de d\u00e9velopper une impl\u00e9mentation vectoris\u00e9e native exploitant la puissance du calcul matriciel dans NumPy et Pandas. La premi\u00e8re \u00e9tape computationnelle consiste \u00e0 d\u00e9river de mani\u00e8re analytique et efficiente les valeurs diagonales de la matrice chapeau <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span>.<\/p>\n<p>Soit <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> la matrice de r\u00e9gression compl\u00e8te de dimension <span class=\"arabpsy-math-inline\" data-latex=\"n \\times p\">$n \\times p$<\/span>, incorporant explicitement la premi\u00e8re colonne constante de uns pour l&rsquo;interception. Si un mod\u00e8le a \u00e9t\u00e9 pr\u00e9alablement ajust\u00e9 avec Statsmodels, cette matrice de conception \u00e9pur\u00e9e est accessible directement via l&rsquo;attribut `modele_ajuste.model.exog`. L&rsquo;\u00e9valuation informatique directe de la matrice de projection globale <span class=\"arabpsy-math-inline\" data-latex=\"H = X(X^T X)^{-1} X^T\">$H = X(X^T X)^{-1} X^T$<\/span> repr\u00e9sente une op\u00e9ration inefficace sur le plan de l&rsquo;allocation m\u00e9moire d\u00e8s lors que <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> grandit, puisqu&rsquo;elle n\u00e9cessite de mat\u00e9rialiser en m\u00e9moire vive un tableau bidimensionnel de dimension <span class=\"arabpsy-math-inline\" data-latex=\"n \\times n\">$n \\times n$<\/span> (occupant pr\u00e8s de 800 m\u00e9gaoctets pour <span class=\"arabpsy-math-inline\" data-latex=\"n = 10,000\">$n = 10,000$<\/span>).<\/p>\n<p>Puisque nous n&rsquo;avons besoin que des seuls \u00e9l\u00e9ments situ\u00e9s sur la diagonale principale de <span class=\"arabpsy-math-inline\" data-latex=\"H\">$H$<\/span>, nous pouvons exploiter les propri\u00e9t\u00e9s du produit matriciel pour d\u00e9river une formulation vectoris\u00e9e hautement optimis\u00e9e. Notons que :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"h_{ii} = [X (X^T X)^{-1} X^T]_{ii} = \\sum_{j=1}^p [X]_{ij} [(X^T X)^{-1} X^T]_{ji} = \\sum_{j=1}^p [X]_{ij} [X (X^T X)^{-1}]_{ij}\">$h_{ii} = [X (X^T X)^{-1} X^T]_{ii} = \\sum_{j=1}^p [X]_{ij} [(X^T X)^{-1} X^T]_{ji} = \\sum_{j=1}^p [X]_{ij} [X (X^T X)^{-1}]_{ij}$<\/span><\/p>\n<p>Cette identit\u00e9 alg\u00e9brique permet de calculer le vecteur colonne des <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> sans jamais instancier la matrice compl\u00e8te de dimension <span class=\"arabpsy-math-inline\" data-latex=\"n \\times n\">$n \\times n$<\/span>. L&rsquo;algorithme se d\u00e9roule computationnellement comme suit :<\/p>\n<ol>\n<li>Extraction du tableau NumPy des variables exog\u00e8nes : <code>X = modele_ajuste.model.exog<\/code>.<\/li>\n<li>Calcul du produit matriciel de Gram et de son inverse par les routines d&rsquo;alg\u00e8bre lin\u00e9aire de NumPy : <code>XtX_inv = np.linalg.inv(X.T @ X)<\/code>.<\/li>\n<li>Calcul du produit interm\u00e9diaire de dimension <span class=\"arabpsy-math-inline\" data-latex=\"n \\times p\">$n \\times p$<\/span> : <code>X_hat = X @ XtX_inv<\/code>.<\/li>\n<li>Sommation vectoris\u00e9e du produit terme \u00e0 terme (produit de Hadamard) le long des colonnes : <code>leverage = np.sum(X_hat * X, axis=1)<\/code>.<\/li>\n<\/ol>\n<p>Cette approche vectoris\u00e9e r\u00e9duit la complexit\u00e9 spatiale d&rsquo;un facteur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{O}(n^2)\">$\\mathcal{O}(n^2)$<\/span> \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{O}(np)\">$\\mathcal{O}(np)$<\/span>, permettant le traitement instantan\u00e9 d&rsquo;\u00e9chantillons de plusieurs dizaines de milliers d&rsquo;observations sur une station de travail standard.<\/p>\n<p>### 8.2 Estimation de l&rsquo;erreur standard d&rsquo;exclusion et finalisation du r\u00e9sidu<\/p>\n<p>Une fois le vecteur des valeurs de levier <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> correctement d\u00e9riv\u00e9, la seconde phase algorithmique op\u00e8re l&rsquo;assemblage matriciel des composantes r\u00e9siduelles pour finaliser l&rsquo;estimation des r\u00e9sidus studentis\u00e9s externes. L&rsquo;extraction des r\u00e9sidus ordinaires bruts s&rsquo;effectue via l&rsquo;attribut `modele_ajuste.resid`, dont nous extrayons le tableau NumPy sous-jacent.<\/p>\n<p>Nous proc\u00e9dons ensuite \u00e0 la d\u00e9rivation de la somme des carr\u00e9s des erreurs r\u00e9siduelles totales :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"SSE = \\sum_{i=1}^n e_i^2 = e^T e\">$SSE = \\sum_{i=1}^n e_i^2 = e^T e$<\/span><\/p>\n<p>Puis de l&rsquo;estimateur sans biais de la variance globale du mod\u00e8le d&rsquo;origine :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"s^2 = \\frac{SSE}{n - p}\">$s^2 = \\frac{SSE}{n &#8211; p}$<\/span><\/p>\n<p>En vertu de l&rsquo;\u00e9quation de mise \u00e0 jour rapide formalis\u00e9e dans la section 3.2, la variance r\u00e9siduelle jackknife <span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2\">$s_{(-i)}^2$<\/span> s&rsquo;obtient pour chaque observation par une instruction enti\u00e8rement vectoris\u00e9e, sans aucune boucle s\u00e9quentielle :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"s_{(-i)}^2 = \\frac{(n - p) s^2 - \\frac{e_i^2}{1 - h_{ii}}}{n - p - 1}\">$s_{(-i)}^2 = \\frac{(n &#8211; p) s^2 &#8211; \\frac{e_i^2}{1 &#8211; h_{ii}}}{n &#8211; p &#8211; 1}$<\/span><\/p>\n<p>Il convient d&rsquo;introduire une garde de s\u00e9curit\u00e9 num\u00e9rique dans le code pour pr\u00e9venir toute instabilit\u00e9 computationnelle (par exemple, si une observation poss\u00e8de un levier th\u00e9orique identique \u00e0 1 ou si une valeur r\u00e9siduelle anormale g\u00e9n\u00e8re un terme r\u00e9siduel d&rsquo;exclusion infinit\u00e9simalement n\u00e9gatif par erreur d&rsquo;arrondi en virgule flottante). Nous appliquons \u00e0 cette fin la fonction `np.maximum(s_ext_sq, 1e-15)` avant l&rsquo;extraction de la racine carr\u00e9e.<\/p>\n<p>Le calcul final du vecteur des r\u00e9sidus studentis\u00e9s externes s&rsquo;\u00e9crit d\u00e8s lors sous la forme compacte :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"t_i = \\frac{e_i}{\\sqrt{s_{(-i)}^2 \\times (1 - h_{ii})}}\">$t_i = \\frac{e_i}{\\sqrt{s_{(-i)}^2 \\times (1 &#8211; h_{ii})}}$<\/span><\/p>\n<p>Pour d\u00e9duire les probabilit\u00e9s critiques bilat\u00e9rales associ\u00e9es \u00e0 chaque statistique <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>, nous mobilisons la fonction de survie (survival function, not\u00e9e `sf`) issue du module `scipy.stats.t`. La fonction de survie, \u00e9quivalente \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"1 - F_{\\mathcal{T}}(|t_i|)\">$1 &#8211; F_{\\mathcal{T}}(|t_i|)$<\/span>, offre une pr\u00e9cision en virgule flottante consid\u00e9rablement sup\u00e9rieure \u00e0 l&rsquo;\u00e9valuation de <span class=\"arabpsy-math-inline\" data-latex=\"1 - \\text{cdf}(x)\">$1 &#8211; \\text{cdf}(x)$<\/span> lorsque <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span> atteint des queues de distribution lointaines o\u00f9 les nombres r\u00e9els subissent une sous-verse computationnelle (*underflow*). Le calcul s&rsquo;ex\u00e9cute ainsi :<\/p>\n<p>`p_valeurs_brutes = 2 * stats.t.sf(np.abs(t_resid), df=n &#8211; p &#8211; 1)`<\/p>\n<p>L&rsquo;ajustement conservateur de Bonferroni s&rsquo;obtient simplement par l&rsquo;op\u00e9ration vectoris\u00e9e :<\/p>\n<p>`p_valeurs_bonf = np.clip(p_valeurs_brutes * n, 0.0, 1.0)`<\/p>\n<p>L&rsquo;ensemble de ces colonnes d\u00e9riv\u00e9es est ensuite encapsul\u00e9 dans un DataFrame Pandas parfaitement align\u00e9 sur l&rsquo;index d&rsquo;origine, reproduisant \u00e0 la perfection les sorties fournies par les routines internes de haut niveau.<\/p>\n<p>### 8.3 Benchmarking et validation crois\u00e9e des calculs manuels<\/p>\n<p>Pour certifier la validit\u00e9 scientifique de l&rsquo;impl\u00e9mentation vectoris\u00e9e personnalis\u00e9e, il est imp\u00e9ratif de proc\u00e9der \u00e0 une validation crois\u00e9e syst\u00e9matique en comparant num\u00e9riquement ses sorties avec celles issues de la fonction de r\u00e9f\u00e9rence `OLSResults.outlier_test()`.<\/p>\n<p>Cette \u00e9tape de test unitaire s&rsquo;appuie sur la fonction `numpy.allclose()`, qui v\u00e9rifie si deux tableaux de r\u00e9els sont num\u00e9riquement \u00e9quivalents dans les limites d&rsquo;une tol\u00e9rance absolue et relative sp\u00e9cifi\u00e9e :<\/p>\n<p>`np.allclose(resultat_manuel[&lsquo;t_student&rsquo;], resultat_statsmodels[&lsquo;student_resid&rsquo;], rtol=1e-10, atol=1e-12)`<\/p>\n<p>Une telle v\u00e9rification garantit que l&rsquo;algorithme manuel ne d\u00e9rive pas sur le plan arithm\u00e9tique et int\u00e8gre fid\u00e8lement toutes les subtilit\u00e9s d&rsquo;arrondis et de degr\u00e9s de libert\u00e9 de l&rsquo;alg\u00e8bre de r\u00e9gression.<\/p>\n<p>Sur le plan de l&rsquo;efficience algorithmique et du benchmarking de performance, l&rsquo;impl\u00e9mentation manuelle vectoris\u00e9e pr\u00e9sente des gains d&rsquo;ex\u00e9cution substantiels sur des jeux de donn\u00e9es volumineux. Lorsque l&rsquo;on soumet un tableau de <span class=\"arabpsy-math-inline\" data-latex=\"n = 100,000\">$n = 100,000$<\/span> observations avec <span class=\"arabpsy-math-inline\" data-latex=\"p = 10\">$p = 10$<\/span> r\u00e9gresseurs \u00e0 un protocole de test de vitesse (via le module `timeit`), l&rsquo;approche vectoris\u00e9e directe s&rsquo;ex\u00e9cute en une fraction de seconde, contournant l&rsquo;instanciation des structures d&rsquo;objets interm\u00e9diaires associ\u00e9es \u00e0 la classe `OLSInfluence`. <\/p>\n<p>Au-del\u00e0 de ces consid\u00e9rations de vitesse brute, la valeur didactique de ce d\u00e9ploiement algorithmique est majeure : elle d\u00e9mythifie la notion de \u00ab\u00a0bo\u00eete noire\u00a0\u00bb diagnostique et ancre d\u00e9finitivement la compr\u00e9hension de l&rsquo;analyste dans les r\u00e9alit\u00e9s math\u00e9matiques et g\u00e9om\u00e9triques sous-tendant la th\u00e9orie de la r\u00e9gression lin\u00e9aire.<\/p>\n<p>## 9. Visualisation diagnostique des r\u00e9sidus studentis\u00e9s pour l&rsquo;\u00e9valuation de mod\u00e8les<\/p>\n<p>### 9.1 Graphique des r\u00e9sidus studentis\u00e9s contre les valeurs ajust\u00e9es<\/p>\n<p>La visualisation graphique des r\u00e9sidus studentis\u00e9s constitue la premi\u00e8re ligne de d\u00e9fense de l&rsquo;analyste pour d\u00e9celer les violations des postulats fondamentaux du mod\u00e8le lin\u00e9aire. Le diagramme fondamental oppose les valeurs pr\u00e9dites par le mod\u00e8le ajust\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"\\hat{Y}_i\">$\\hat{Y}_i$<\/span>) sur l&rsquo;axe des abscisses aux r\u00e9sidus studentis\u00e9s externes (<span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>) sur l&rsquo;axe des ordonn\u00e9es.<\/p>\n<p>Pour conf\u00e9rer \u00e0 ce graphique une port\u00e9e diagnostique imm\u00e9diate, il est indispensable de tracer plusieurs lignes de r\u00e9f\u00e9rence horizontales conventionnelles :<\/p>\n<ul>\n<li>Une ligne m\u00e9diane solide positionn\u00e9e \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"y = 0\">$y = 0$<\/span>, mat\u00e9rialisant l&rsquo;esp\u00e9rance th\u00e9orique des r\u00e9sidus sous l&rsquo;hypoth\u00e8se d&rsquo;une sp\u00e9cification lin\u00e9aire correcte.<\/li>\n<li>Deux lignes en tirets \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"y = +2\">$y = +2$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"y = -2\">$y = -2$<\/span>, d\u00e9limitant la bande de dispersion normale attendue pour environ 95 % des observations dans un cadre gaussien.<\/li>\n<li>Deux lignes d&rsquo;alerte critique en pointill\u00e9s \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"y = +3\">$y = +3$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"y = -3\">$y = -3$<\/span>, au-del\u00e0 desquelles les observations sont formellement isol\u00e9es pour un audit approfondi.<\/li>\n<\/ul>\n<p>L&rsquo;inspection morphologique de ce nuage de points permet de porter deux diagnostics critiques :<\/p>\n<ol>\n<li><strong>D\u00e9tection des non-lin\u00e9arit\u00e9s structurales<\/strong> : Dans un mod\u00e8le correctement calibr\u00e9, les r\u00e9sidus studentis\u00e9s doivent se distribuer de mani\u00e8re al\u00e9atoire et uniforme, formant une bande de dispersion horizontale homog\u00e8ne sans structure discernable. Si le nuage r\u00e9siduel adopte un profil curviligne distinct \u2014 tel qu&rsquo;une trajectoire parabolique en U ou en U invers\u00e9 \u2014, cela r\u00e9v\u00e8le formellement que la fonction de r\u00e9gression omet une composante quadratique ou une transformation non lin\u00e9aire essentielle d&rsquo;un des pr\u00e9dicteurs.<\/li>\n<li><strong>D\u00e9tection de l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9<\/strong> : L&rsquo;apparition d&rsquo;un patron morphologique en entonnoir (*funnel shape*), caract\u00e9ris\u00e9 par un \u00e9largissement ou un resserrement progressif de la dispersion des r\u00e9sidus au fur et \u00e0 mesure que la valeur pr\u00e9dite <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{Y}\">$\\hat{Y}$<\/span> cro\u00eet, signale la violation de l&rsquo;hypoth\u00e8se de variance constante. Les r\u00e9sidus studentis\u00e9s externes, \u00e9tant pr\u00e9alablement corrig\u00e9s de la structure de levier horizontal, refl\u00e8tent une v\u00e9ritable h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 du processus stochastique, invitant l&rsquo;analyste \u00e0 appliquer des estimateurs de variance-covariance robustes (tels que les estimateurs sandwich de White-Huber HC3) ou \u00e0 recourir aux moindres carr\u00e9s pond\u00e9r\u00e9s (WLS).<\/li>\n<\/ol>\n<p>### 9.2 Diagramme quantile-quantile (Q-Q plot) avec enveloppes de confiance<\/p>\n<p>L&rsquo;\u00e9valuation de la normalit\u00e9 des perturbations stochastiques s&rsquo;op\u00e8re de mani\u00e8re privil\u00e9gi\u00e9e via la construction d&rsquo;un diagramme quantile-quantile (Q-Q plot). Cette repr\u00e9sentation graphique projette les quantiles empiriques ordonn\u00e9s des r\u00e9sidus studentis\u00e9s calcul\u00e9s sur l&rsquo;\u00e9chantillon contre les quantiles th\u00e9oriques attendus d&rsquo;une loi de Student \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"n - p - 1\">$n &#8211; p &#8211; 1$<\/span> degr\u00e9s de libert\u00e9 (ou d&rsquo;une distribution normale standard lorsque les degr\u00e9s de libert\u00e9 r\u00e9siduels sont \u00e9lev\u00e9s).<\/p>\n<p>L&rsquo;\u00e9cosyst\u00e8me Statsmodels facilite la g\u00e9n\u00e9ration de ces visualisations gr\u00e2ce \u00e0 la classe `ProbPlot` issue du sous-module `statsmodels.graphics.gofplots`. Lorsque les postulats distributionnels sont rigoureusement satisfaits, les points empiriques s&rsquo;alignent \u00e9troitement le long de la bissectrice th\u00e9orique diagonale rouge \u00e0 45 degr\u00e9s.<\/p>\n<p>Les d\u00e9viations syst\u00e9matiques par rapport \u00e0 cette droite de r\u00e9f\u00e9rence fournissent des indications morphologiques pr\u00e9cises quant \u00e0 la nature des pathologies de distribution affectant les donn\u00e9es :<\/p>\n<ul>\n<li><strong>Distributions \u00e0 queues lourdes (leptokurtose)<\/strong> : Si le trac\u00e9 r\u00e9siduel prend une forme caract\u00e9ristique en \u00ab S \u00bb, s&rsquo;\u00e9cartant vers le haut de la droite \u00e0 l&rsquo;extr\u00e9mit\u00e9 droite et plongeant sous la droite \u00e0 l&rsquo;extr\u00e9mit\u00e9 gauche, cela d\u00e9montre une concentration de valeurs extr\u00eames sup\u00e9rieure \u00e0 celle pr\u00e9dite par la distribution normale. Ce patron d&rsquo;aplatissement excessif (kurtosis \u00e9lev\u00e9) avertit l&rsquo;analyste que les tests d&rsquo;inf\u00e9rence standards fond\u00e9s sur la normalit\u00e9 sous-estiment syst\u00e9matiquement la probabilit\u00e9 des \u00e9v\u00e9nements extr\u00eames.<\/li>\n<li><strong>Asym\u00e9trie distributionnelle (skewness)<\/strong> : Si les r\u00e9sidus s&rsquo;\u00e9cartent unilat\u00e9ralement de la droite de r\u00e9f\u00e9rence \u2014 par exemple en s&rsquo;arquant de fa\u00e7on convexe vers le haut \u2014, cela indique une asym\u00e9trie positive persistante des donn\u00e9es r\u00e9siduelles, sugg\u00e9rant l&rsquo;opportunit\u00e9 d&rsquo;appliquer une transformation stabilisatrice (telle qu&rsquo;une transformation logarithmique ou une transformation de Box-Cox) sur la variable d\u00e9pendante <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>.<\/li>\n<\/ul>\n<p>Pour consolider la d\u00e9cision visuelle, l&rsquo;adjonction d&rsquo;enveloppes de confiance par r\u00e9\u00e9chantillonnage bootstrap ou simulations de Monte-Carlo permet de mat\u00e9rialiser les bornes de variabilit\u00e9 acceptable au niveau des queues de distribution, pr\u00e9venant ainsi la surinterpr\u00e9tation d&rsquo;ondulations r\u00e9siduelles mineures aux extr\u00e9mit\u00e9s du diagramme.<\/p>\n<p>### 9.3 Graphique combin\u00e9 influence-levier et seuils d&rsquo;alerte<\/p>\n<p>L&rsquo;int\u00e9gration conclusive des m\u00e9triques diagnostiques s&rsquo;incarne dans la construction d&rsquo;un graphique combin\u00e9 croisant les valeurs de levier <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> sur l&rsquo;axe des abscisses et les r\u00e9sidus studentis\u00e9s (id\u00e9alement port\u00e9s au carr\u00e9 pour accentuer les \u00e9carts, ou conserv\u00e9s avec leur signe) sur l&rsquo;axe des ordonn\u00e9es. Cette typologie de visualisation \u00e0 bulles permet de d\u00e9composer instantan\u00e9ment l&rsquo;espace g\u00e9om\u00e9trique en quatre quadrants distincts :<\/p>\n<ol>\n<li><strong>Quadrant inf\u00e9rieur gauche (Faible levier, Faible r\u00e9sidu)<\/strong> : Zone de conformit\u00e9 nominale o\u00f9 se concentre la vaste majorit\u00e9 des participants respectant fid\u00e8lement le mod\u00e8le.<\/li>\n<li><strong>Quadrant sup\u00e9rieur gauche (Faible levier, Fort r\u00e9sidu)<\/strong> : Zone des valeurs aberrantes verticales (*outliers* purs). Ces observations s&rsquo;\u00e9cartent radicalement de la surface de pr\u00e9diction mais, du fait de leur proximit\u00e9 avec le barycentre des r\u00e9gresseurs, n&rsquo;exercent qu&rsquo;une force de levier limit\u00e9e sur la rotation du plan d&rsquo;ajustement.<\/li>\n<li><strong>Quadrant inf\u00e9rieur droit (Fort levier, Faible r\u00e9sidu)<\/strong> : Zone des points d&rsquo;ancrage concordants (*good leverage points*). Ces individus pr\u00e9sentent des configurations de covariables exceptionnelles mais se positionnent harmonieusement sur la trajectoire naturelle du mod\u00e8le. Loin de d\u00e9grader l&rsquo;estimation, ils renforcent la pr\u00e9cision des coefficients en \u00e9tendant artificiellement la dispersion du design exp\u00e9rimental.<\/li>\n<li><strong>Quadrant sup\u00e9rieur droit (Fort levier, Fort r\u00e9sidu)<\/strong> : Zone critique des points influents destructeurs (*bad leverage points*). Ces observations combinent une excentricit\u00e9 spatiale aigu\u00eb et un d\u00e9saccord violent avec la droite de r\u00e9gression. Ce sont ces points qui basculent l&rsquo;estimation MCO et modifient radicalement les conclusions de l&rsquo;\u00e9tude.<\/li>\n<\/ol>\n<p>La fonction native `statsmodels.graphics.regressionplots.plot_leverage_resid2` fournit une impl\u00e9mentation standardis\u00e9e de ce concept, tra\u00e7ant les leviers contre le carr\u00e9 des r\u00e9sidus studentis\u00e9s. L&rsquo;analyste averti enrichit fr\u00e9quemment ce graphique en proportionnant le diam\u00e8tre des marqueurs \u00e0 la distance de Cook (<span class=\"arabpsy-math-inline\" data-latex=\"D_i\">$D_i$<\/span>) et en incrustant des annotations textuelles automatiques ciblant les identifiants uniques des sujets qui violent simultan\u00e9ment les seuils critiques de levier (<span class=\"arabpsy-math-inline\" data-latex=\"h_{ii} &gt; 2p\/n\">$h_{ii} &gt; 2p\/n$<\/span>) et de r\u00e9sidu studentis\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"|t_i| &gt; 3\">$|t_i| &gt; 3$<\/span>).<\/p>\n<p>## 10. Protocoles de traitement des observations atypiques en recherche quantitative<\/p>\n<p>### 10.1 Strat\u00e9gies d&rsquo;investigation pr\u00e9-d\u00e9cisionnelles<\/p>\n<p>La d\u00e9tection objective d&rsquo;un r\u00e9sidu studentis\u00e9 externe significatif inaugure obligatoirement une phase d&rsquo;audit pr\u00e9-d\u00e9cisionnel rigoureusement encadr\u00e9e. La premi\u00e8re d\u00e9marche consiste en une v\u00e9rification r\u00e9trospective scrupuleuse des journaux d&rsquo;acquisition de donn\u00e9es (*data logs*). Dans les environnements exp\u00e9rimentaux contemporains coupl\u00e9s \u00e0 des interfaces informatis\u00e9es ou des biomarqueurs physiologiques (eyetracking, EEG, actim\u00e9trie), il est fr\u00e9quent qu&rsquo;un pic d&rsquo;atypicit\u00e9 co\u00efncide avec un artefact technique document\u00e9 : perte temporaire de signal, interf\u00e9rence \u00e9lectromagn\u00e9tique, recalibration d\u00e9faillante d&rsquo;un capteur, ou interruption intempestive de la session par un \u00e9v\u00e9nement ext\u00e9rieur. Dans une telle hypoth\u00e8se, l&rsquo;observation est entach\u00e9e d&rsquo;une invalidit\u00e9 instrumentale indiscutable, justifiant pleinement sa neutralisation.<\/p>\n<p>Lorsque l&rsquo;int\u00e9grit\u00e9 de la collecte est mat\u00e9riellement attest\u00e9e, l&rsquo;analyste explore le profil psychom\u00e9trique ou comportemental holistique du participant. L&rsquo;\u00e9valuation d&rsquo;indices de r\u00e9ponse inattentive \u2014 tels que la longueur maximale de r\u00e9ponses identiques cons\u00e9cutives (*long-string index*), la vitesse anormale d&rsquo;ex\u00e9cution d\u00e9fiant les temps incompressibles de traitement cognitif (temps de r\u00e9action inf\u00e9rieurs \u00e0 150 millisecondes), ou une discordance aigu\u00eb mesur\u00e9e par la distance de Mahalanobis intra-individuelle sur des \u00e9chelles invers\u00e9es \u2014 permet d&rsquo;\u00e9tayer empiriquement l&rsquo;hypoth\u00e8se d&rsquo;une r\u00e9ponse non investie.<\/p>\n<p>Enfin, l&rsquo;ensemble de ces d\u00e9marches doit \u00eatre confront\u00e9 aux crit\u00e8res d&rsquo;exclusion sp\u00e9cifi\u00e9s a priori dans le protocole de pr\u00e9-enregistrement de la recherche (sur des plateformes comme l&rsquo;Open Science Framework ou AsPredicted). La science ouverte exige que les crit\u00e8res d\u00e9finissant une valeur aberrante excluable soient formalis\u00e9s avant la phase d&rsquo;acquisition des donn\u00e9es. Toute exclusion post-hoc non stipul\u00e9e en amont du protocole doit \u00eatre explicitement d\u00e9sign\u00e9e comme une d\u00e9marche exploratoire et faire l&rsquo;objet d&rsquo;une justification th\u00e9orique exhaustive.<\/p>\n<p>### 10.2 Analyse de sensibilit\u00e9 et mod\u00e9lisation robuste<\/p>\n<p>Face \u00e0 la persistance d&rsquo;observations atypiques dont l&rsquo;\u00e9viction d\u00e9finitive ne peut \u00eatre tranch\u00e9e sans \u00e9quivoque, le standard m\u00e9thodologique absolu r\u00e9side dans la r\u00e9alisation syst\u00e9matique d&rsquo;une analyse de sensibilit\u00e9 coupl\u00e9e. Cette proc\u00e9dure consiste \u00e0 ajuster et rapporter en miroir deux mod\u00e8les de r\u00e9gression distincts : le mod\u00e8le primaire calibr\u00e9 sur l&rsquo;int\u00e9gralit\u00e9 de l&rsquo;\u00e9chantillon collect\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span>), et le mod\u00e8le secondaire r\u00e9ajust\u00e9 sur le sous-\u00e9chantillon \u00e9pur\u00e9 des observations isol\u00e9es par le test des r\u00e9sidus studentis\u00e9s (<span class=\"arabpsy-math-inline\" data-latex=\"n - k\">$n &#8211; k$<\/span>).<\/p>\n<p>L&rsquo;analyste examine ensuite la stabilit\u00e9 topologique des estimations en \u00e9valuant la variation relative des coefficients de r\u00e9gression :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\Delta \\beta_j = \\frac{\\hat{\\beta}_j^{(\\text{complet})} - \\hat{\\beta}_j^{(\\text{\u00e9pur\u00e9})}}{\\hat{\\beta}_j^{(\\text{complet})}} \\times 100\">$\\Delta \\beta_j = \\frac{\\hat{\\beta}_j^{(\\text{complet})} &#8211; \\hat{\\beta}_j^{(\\text{\u00e9pur\u00e9})}}{\\hat{\\beta}_j^{(\\text{complet})}} \\times 100$<\/span><\/p>\n<p>Si les conclusions substantielles de la recherche (signe des relations, ampleur clinique des effets, niveaux de significativit\u00e9 des tests d&rsquo;hypoth\u00e8se) demeurent invariantes entre les deux versions, le mod\u00e8le est qualifi\u00e9 de robuste \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 d&rsquo;\u00e9chantillonnage, d\u00e9samor\u00e7ant les doutes quant \u00e0 un potentiel biais d&rsquo;inf\u00e9rence. Si, \u00e0 l&rsquo;inverse, l&rsquo;exclusion d&rsquo;une seule observation atypique suffit \u00e0 faire basculer un pr\u00e9dicteur de la significativit\u00e9 statistique \u00e0 la non-significativit\u00e9, l&rsquo;effet d\u00e9couvert doit \u00eatre consid\u00e9r\u00e9 comme structurellement fragile et interpr\u00e9t\u00e9 avec une r\u00e9serve critique imp\u00e9rative.<\/p>\n<p>Une alternative analytique \u00e9l\u00e9gante \u00e0 l&rsquo;exclusion binaire consiste \u00e0 abandonner les moindres carr\u00e9s ordinaires au profit des m\u00e9thodes d&rsquo;estimation lin\u00e9aire robuste (RLM, *Robust Linear Models*). L&rsquo;estimateur <span class=\"arabpsy-math-inline\" data-latex=\"M\">$M$<\/span> de Huber, directement accessible via la classe `statsmodels.api.RLM`, substitue \u00e0 la fonction de perte quadratique <span class=\"arabpsy-math-inline\" data-latex=\"\\frac{1}{2}e^2\">$\\frac{1}{2}e^2$<\/span> une fonction de co\u00fbt hybride qui p\u00e9nalise lin\u00e9airement et non plus quadratiquement les r\u00e9sidus au-del\u00e0 d&rsquo;un seuil de d\u00e9viation param\u00e9trable (usuellement <span class=\"arabpsy-math-inline\" data-latex=\"c = 1,345\">$c = 1,345$<\/span>). Pour les jeux de donn\u00e9es lourdement contamin\u00e9s, les pond\u00e9rations par bivariance de Tukey (*Tukey&rsquo;s biweight*) permettent d&rsquo;annihiler compl\u00e8tement le poids accord\u00e9 aux valeurs extr\u00eames tout en conservant l&rsquo;int\u00e9gralit\u00e9 de l&rsquo;\u00e9chantillon au sein d&rsquo;un formalisme it\u00e9ratif \u00e0 moindres carr\u00e9s repond\u00e9r\u00e9s (IRLS).<\/p>\n<p>Enfin, lorsque les d\u00e9viations r\u00e9siduelles r\u00e9sultent d&rsquo;une distribution sous-jacente structurellement asym\u00e9trique ou \u00e0 queues lourdes r\u00e9fractaire \u00e0 toute normalisation, l&rsquo;analyste peut d\u00e9ployer la r\u00e9gression quantile (via `statsmodels.formula.api.quantreg`). En mod\u00e9lisant la m\u00e9diane conditionnelle (quantile 0,50) par minimisation des d\u00e9viations absolues (m\u00e9thode <span class=\"arabpsy-math-inline\" data-latex=\"L_1\">$L_1$<\/span>), cette approche offre une insensibilit\u00e9 math\u00e9matique totale aux r\u00e9sidus studentis\u00e9s extr\u00eames sur la variable d\u00e9pendante.<\/p>\n<p>### 10.3 Normes de transparence et directives de r\u00e9daction acad\u00e9mique<\/p>\n<p>Les directives \u00e9ditoriales acad\u00e9miques contemporaines convergent vers une exigence intransigeante de transparence m\u00e9thodologique. Les recommandations formelles de l&rsquo;APA (Publication Manual of the American Psychological Association, 7e \u00e9dition, section Journal Article Reporting Standards ou JARS) prescrivent explicitement de documenter :<\/p>\n<ul>\n<li>La liste exhaustive des tests diagnostiques r\u00e9siduels conduits lors des phases pr\u00e9liminaires de l&rsquo;analyse.<\/li>\n<li>Le nombre pr\u00e9cis d&rsquo;observations identifi\u00e9es comme discordantes, en sp\u00e9cifiant les m\u00e9triques utilis\u00e9es (valeurs exactes des r\u00e9sidus studentis\u00e9s externes <span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>, leviers <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> et distances de Cook <span class=\"arabpsy-math-inline\" data-latex=\"D_i\">$D_i$<\/span>).<\/li>\n<li>La justification explicite et motiv\u00e9e sous-tendant le maintien ou le rejet de chaque cas atypique.<\/li>\n<li>La restitution exhaustive des estimations comparatives issues des analyses de sensibilit\u00e9 afin de pr\u00e9munir les pairs contre tout embellissement artificiel des donn\u00e9es.<\/li>\n<\/ul>\n<p>Sur le plan de l&rsquo;\u00e9thique computationnelle et de la science ouverte, le chercheur contemporain consolide son manuscrit en publiant l&rsquo;int\u00e9gralit\u00e9 de ses scripts de nettoyage et d&rsquo;analyse Python sur des plateformes d&rsquo;archivage p\u00e9rennes et certifi\u00e9es (telles que le Curation and Preservation Repository Zenodo ou l&rsquo;Open Science Framework). Le d\u00e9p\u00f4t public d&rsquo;un notebook Jupyter document\u00e9, assorti de son environnement d&rsquo;ex\u00e9cution explicite (fichier `requirements.txt` ou `environment.yml`), garantit la r\u00e9plicabilit\u00e9 parfaite de la cha\u00eene de d\u00e9tection des r\u00e9sidus studentis\u00e9s, scellant ainsi l&rsquo;int\u00e9grit\u00e9 d\u00e9ontologique de la recherche publi\u00e9e.<\/p>\n<p>## 11. Extensions analytiques : Mod\u00e8les multivari\u00e9s et structures hi\u00e9rarchiques<\/p>\n<p>### 11.1 R\u00e9gression lin\u00e9aire multiple et colin\u00e9arit\u00e9 des pr\u00e9dicteurs<\/p>\n<p>La complexit\u00e9 math\u00e9matique et diagnostique des r\u00e9sidus studentis\u00e9s s&rsquo;accro\u00eet substantiellement lorsque le mod\u00e8le de r\u00e9gression incorpore un vecteur multidimensionnel de covariables corr\u00e9l\u00e9es entre elles. Dans un cadre de r\u00e9gression lin\u00e9aire multiple comportant <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> pr\u00e9dicteurs, la pr\u00e9sence de multicolin\u00e9arit\u00e9 s\u00e9v\u00e8re \u2014 quantifi\u00e9e par des facteurs d&rsquo;inflation de la variance \u00e9lev\u00e9s (<span class=\"arabpsy-math-inline\" data-latex=\"VIF_j &gt; 5\">$VIF_j &gt; 5$<\/span> ou $10$) \u2014 alt\u00e8re profond\u00e9ment la structure de la matrice de projection <span class=\"arabpsy-math-inline\" data-latex=\"H\">$H$<\/span> et la r\u00e9partition spatiale des valeurs de levier <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span>.<\/p>\n<p>Lorsque deux ou plusieurs pr\u00e9dicteurs partagent une intercorr\u00e9lation quasi-parfaite, le d\u00e9terminant de la matrice de Gram <span class=\"arabpsy-math-inline\" data-latex=\"X^T X\">$X^T X$<\/span> tend vers z\u00e9ro, provoquant une instabilit\u00e9 num\u00e9rique majeure dans le calcul de son inverse <span class=\"arabpsy-math-inline\" data-latex=\"(X^T X)^{-1}\">$(X^T X)^{-1}$<\/span>. Les erreurs-types des coefficients <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}_j\">$\\hat{\\beta}_j$<\/span> subissent alors une inflation dramatique, ce qui d\u00e9forme indirectement l&rsquo;\u00e9valuation des r\u00e9sidus studentis\u00e9s. Dans cette g\u00e9om\u00e9trie contract\u00e9e, une observation peut pr\u00e9senter une valeur parfaitement standard sur chacune des covariables examin\u00e9es isol\u00e9ment (dans des analyses descriptives univari\u00e9es classiques), mais constituer une singularit\u00e9 aberrante dans l&rsquo;espace multidimensionnel conjoint des r\u00e9gresseurs. <\/p>\n<p>Ce ph\u00e9nom\u00e8ne d&rsquo;atypicit\u00e9 cach\u00e9e n&rsquo;est d\u00e9tectable que par l&rsquo;interm\u00e9diaire de la matrice chapeau globale, qui capture l&rsquo;\u00e9loignement d&rsquo;une observation par rapport au plan factoriel commun. Les r\u00e9sidus studentis\u00e9s externes calcul\u00e9s dans ces contextes de forte colin\u00e9arit\u00e9 deviennent extr\u00eamement sensibles aux variations marginales d&rsquo;\u00e9chantillonnage. <\/p>\n<p>Par ailleurs, l&rsquo;inclusion dans le mod\u00e8le de termes d&rsquo;interactions d&rsquo;ordre sup\u00e9rieur (mod\u00e9rations multiplicatives de type <span class=\"arabpsy-math-inline\" data-latex=\"X_1 \\times X_2\">$X_1 \\times X_2$<\/span>) ou de composantes polynomiales g\u00e9n\u00e8re m\u00e9caniquement des points de fort levier structurels aux confins de la distribution. L&rsquo;analyste doit imp\u00e9rativement centrer ou standardiser les variables continues en amont du calcul des produits afin de minimiser la colin\u00e9arit\u00e9 non essentielle entre les termes simples et leurs interactions, pr\u00e9servant ainsi la validit\u00e9 distributionnelle du calcul des r\u00e9sidus studentis\u00e9s externes.<\/p>\n<p>### 11.2 Adaptation aux mod\u00e8les lin\u00e9aires \u00e0 effets mixtes (LMM)<\/p>\n<p>Dans les devis exp\u00e9rimentaux comportant des mesures r\u00e9p\u00e9t\u00e9es, des protocoles longitudinaux ou des cohortes comportementales hi\u00e9rarchis\u00e9es (\u00e9l\u00e8ves nich\u00e9s au sein d&rsquo;\u00e9coles, patients assign\u00e9s \u00e0 diff\u00e9rents centres hospitaliers), l&rsquo;ind\u00e9pendance stochastique des perturbations est structurellement viol\u00e9e. Le recours aux mod\u00e8les lin\u00e9aires \u00e0 effets mixtes (LMM, *Linear Mixed Models*) devient alors obligatoire :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"Y = X\\beta + Z u + \\varepsilon\">$Y = X\\beta + Z u + \\varepsilon$<\/span><\/p>\n<p>O\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> correspond \u00e0 la matrice de conception des effets fixes, <span class=\"arabpsy-math-inline\" data-latex=\"Z\">$Z$<\/span> d\u00e9signe la matrice des effets al\u00e9atoires associ\u00e9e aux grappes, <span class=\"arabpsy-math-inline\" data-latex=\"u \\sim \\mathcal{N}(0, G)\">$u \\sim \\mathcal{N}(0, G)$<\/span> repr\u00e9sente le vecteur des \u00e9carts al\u00e9atoires de groupe, et <span class=\"arabpsy-math-inline\" data-latex=\"\\varepsilon \\sim \\mathcal{N}(0, R)\">$\\varepsilon \\sim \\mathcal{N}(0, R)$<\/span> incarne le vecteur des erreurs conditionnelles r\u00e9siduelles.<\/p>\n<p>Dans cette architecture \u00e0 deux niveaux de variance, la d\u00e9finition et le calcul des r\u00e9sidus studentis\u00e9s se scindent en deux paradigmes analytiques distincts :<\/p>\n<ol>\n<li><strong>Les r\u00e9sidus marginaux (ou de population)<\/strong> : Calcul\u00e9s \u00e0 partir des seuls effets fixes globaux, <span class=\"arabpsy-math-inline\" data-latex=\"e_m = Y - X\\hat{\\beta}\">$e_m = Y &#8211; X\\hat{\\beta}$<\/span>. Ces r\u00e9sidus \u00e9valuent la d\u00e9viation d&rsquo;un individu par rapport \u00e0 la trajectoire moyenne de la population globale, ignorant son appartenance \u00e0 une grappe sp\u00e9cifique.<\/li>\n<li><strong>Les r\u00e9sidus conditionnels<\/strong> : Calcul\u00e9s en retranchant les effets fixes et la contribution des effets al\u00e9atoires sp\u00e9cifiques pr\u00e9dits au niveau du groupe (BLUPs, *Best Linear Unbiased Predictors*) : <span class=\"arabpsy-math-inline\" data-latex=\"e_c = Y - X\\hat{\\beta} - Z\\hat{u}\">$e_c = Y &#8211; X\\hat{\\beta} &#8211; Z\\hat{u}$<\/span>. Ces r\u00e9sidus \u00e9valuent l&rsquo;\u00e9cart entre la performance d&rsquo;un individu et la trajectoire sp\u00e9cifique de son groupe d&rsquo;appartenance.<\/li>\n<\/ol>\n<p>La studentisation de ces r\u00e9sidus au sein de l&rsquo;environnement Statsmodels (via le module `statsmodels.regression.mixed_linear_model.MixedLM`) requiert une d\u00e9composition matricielle complexe de la matrice de variance-covariance totale <span class=\"arabpsy-math-inline\" data-latex=\"V = Z G Z^T + R\">$V = Z G Z^T + R$<\/span>. La matrice de projection conditionnelle incorpore \u00e0 la fois la projection des effets fixes et la p\u00e9nalisation de cr\u00eate li\u00e9e \u00e0 la dispersion des effets al\u00e9atoires. La d\u00e9tection de r\u00e9sidus conditionnels studentis\u00e9s significatifs permet d&rsquo;isoler des trajectoires individuelles pathologiques au sein d&rsquo;une cohorte clinique donn\u00e9e, tandis que l&rsquo;examen des r\u00e9sidus de niveau sup\u00e9rieur (niveau grappe) permet d&rsquo;identifier des centres de recherche entiers dont les donn\u00e9es d\u00e9rivent de mani\u00e8re suspecte par rapport au protocole d&rsquo;essai multicentrique.<\/p>\n<p>## 12. \u00c9tude de cas compl\u00e8te et reproductible de bout en bout<\/p>\n<p>### 12.1 Mise en situation empirique et simulation de la cohorte<\/p>\n<p>Afin de cristalliser l&rsquo;ensemble des concepts th\u00e9oriques et des m\u00e9thodes algorithmiques abord\u00e9s, nous d\u00e9ployons une \u00e9tude de cas clinique compl\u00e8te, rigoureuse et enti\u00e8rement reproductible. Consid\u00e9rons un contexte de recherche en psychologie de la sant\u00e9 au travail analysant l&rsquo;\u00e9puisement professionnel (*burnout*) chez les personnels soignants en milieu hospitalier. Le protocole \u00e9value la pr\u00e9diction de l&rsquo;indice de bien-\u00eatre psychologique subjectif (variable d\u00e9pendante <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>, \u00e9talonn\u00e9e sur une \u00e9chelle continue de 0 \u00e0 100) en fonction de deux pr\u00e9dicteurs fondamentaux : le niveau de charge \u00e9motionnelle per\u00e7ue (score d&rsquo;\u00e9puisement continu) et le score global sur l&rsquo;\u00e9chelle de r\u00e9silience psychologique individuelle de Connor-Davidson.<\/p>\n<p>L&rsquo;\u00e9quation th\u00e9orique structurelle du processus de mod\u00e9lisation s&rsquo;\u00e9nonce :<\/p>\n<p><span class=\"arabpsy-math-inline\" data-latex=\"\\text{BienEtre}_i = \\beta_0 + \\beta_1 , \\text{Epuisement}_i + \\beta_2 , \\text{Resilience}_i + \\varepsilon_i\">$\\text{BienEtre}_i = \\beta_0 + \\beta_1 , \\text{Epuisement}_i + \\beta_2 , \\text{Resilience}_i + \\varepsilon_i$<\/span><\/p>\n<p>Nous mod\u00e9lisons un \u00e9chantillon de <span class=\"arabpsy-math-inline\" data-latex=\"n = 150\">$n = 150$<\/span> professionnels de sant\u00e9. Pour \u00e9prouver l&rsquo;efficacit\u00e9 diagnostique de notre pipeline, nous ins\u00e9rons d\u00e9lib\u00e9r\u00e9ment dans le processus g\u00e9n\u00e9rateur de donn\u00e9es deux anomalies cliniquement plausibles :<\/p>\n<ul>\n<li><strong>Le cas n\u00b0 42<\/strong> : Une valeur aberrante verticale pure. Ce soignant affiche des niveaux d&rsquo;\u00e9puisement et de r\u00e9silience rigoureusement conformes \u00e0 la moyenne de la cohorte, mais enregistre un score de bien-\u00eatre catastrophique proche de z\u00e9ro en raison d&rsquo;un \u00e9v\u00e9nement traumatique personnel survenu la veille de l&rsquo;\u00e9valuation (donn\u00e9e valide \u00e9cologiquement mais discordante vis-\u00e0-vis du mod\u00e8le).<\/li>\n<li><strong>Le cas n\u00b0 118<\/strong> : Un point de fort levier discordant (*bad leverage point*). Ce soignant cumule un niveau d&rsquo;\u00e9puisement exceptionnellement \u00e9lev\u00e9 coupl\u00e9 \u00e0 une r\u00e9silience minimale, mais affiche paradoxalement un bien-\u00eatre subjectif maximal imputable \u00e0 un biais de d\u00e9sirabilit\u00e9 sociale extr\u00eame ou \u00e0 un comportement de d\u00e9ni document\u00e9.<\/li>\n<\/ul>\n<p>### 12.2 Pipeline complet d&rsquo;analyse et d&rsquo;extraction diagnostique<\/p>\n<p>Le protocole analytique se d\u00e9ploie \u00e0 travers un pipeline s\u00e9quentiel continu assurant l&rsquo;ajustement du mod\u00e8le par les moindres carr\u00e9s ordinaires, l&rsquo;extraction matricielle des diagnostics, le calcul formel des r\u00e9sidus studentis\u00e9s externes via `outlier_test()`, et la consolidation des m\u00e9triques de levier et de distance de Cook.<\/p>\n<p>Le tableau ci-dessous illustre la structure des donn\u00e9es diagnostiques obtenue pour les observations cl\u00e9s de l&rsquo;\u00e9chantillon apr\u00e8s ex\u00e9cution du pipeline d&rsquo;analyse :<\/p>\n<p>| Identifiant (#) | R\u00e9sidu brut (<span class=\"arabpsy-math-inline\" data-latex=\"e_i\">$e_i$<\/span>) | Levier (<span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span>) | R\u00e9sidu studentis\u00e9 externe (<span class=\"arabpsy-math-inline\" data-latex=\"t_i\">$t_i$<\/span>) | Probabilit\u00e9 brute (`unadj_p`) | Probabilit\u00e9 ajust\u00e9e (`bonf(p)`) | Distance de Cook (<span class=\"arabpsy-math-inline\" data-latex=\"D_i\">$D_i$<\/span>) | Statut Diagnostique |<br \/>\n| :&#8212;: | :&#8212;: | :&#8212;: | :&#8212;: | :&#8212;: | :&#8212;: | :&#8212;: | :&#8212;: |<br \/>\n| **Observation 12** | <span class=\"arabpsy-math-inline\" data-latex=\"+1,42\">$+1,42$<\/span> | $0,018$ | <span class=\"arabpsy-math-inline\" data-latex=\"+0,48\">$+0,48$<\/span> | $0,631$ | $1,000$ | $0,001$ | Conforme |<br \/>\n| **Observation 42** | <span class=\"arabpsy-math-inline\" data-latex=\"-28,75\">$-28,75$<\/span> | $0,019$ | <span class=\"arabpsy-math-inline\" data-latex=\"-4,12\">$-4,12$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"6,2 \\times 10^{-5}\">$6,2 \\times 10^{-5}$<\/span> | $0,009$ | $0,112$ | Outlier vertical pur |<br \/>\n| **Observation 75** | <span class=\"arabpsy-math-inline\" data-latex=\"+2,15\">$+2,15$<\/span> | $0,022$ | <span class=\"arabpsy-math-inline\" data-latex=\"+0,73\">$+0,73$<\/span> | $0,466$ | $1,000$ | $0,004$ | Conforme |<br \/>\n| **Observation 118** | <span class=\"arabpsy-math-inline\" data-latex=\"+34,12\">$+34,12$<\/span> | $0,085$ | <span class=\"arabpsy-math-inline\" data-latex=\"+5,31\">$+5,31$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"3,8 \\times 10^{-7}\">$3,8 \\times 10^{-7}$<\/span> | $0,00005$ | $0,782$ | Point influent critique |<br \/>\n| **Observation 149** | <span class=\"arabpsy-math-inline\" data-latex=\"-0,88\">$-0,88$<\/span> | $0,015$ | <span class=\"arabpsy-math-inline\" data-latex=\"-0,29\">$-0,29$<\/span> | $0,772$ | $1,000$ | $0,000$ | Conforme |<\/p>\n<p>L&rsquo;examen du tableau permet de constater imm\u00e9diatement l&rsquo;efficience de la studentisation externe coupl\u00e9e \u00e0 la correction de Bonferroni. Les observations nominales (comme les sujets 12, 75 et 149) pr\u00e9sentent des r\u00e9sidus studentis\u00e9s compris dans l&rsquo;intervalle canonique <span class=\"arabpsy-math-inline\" data-latex=\"[-1, +1]\">$[-1, +1]$<\/span> et des probabilit\u00e9s critiques corrig\u00e9es \u00e9gales \u00e0 1, attestant de leur parfaite concordance avec le mod\u00e8le lin\u00e9aire.<\/p>\n<p>\u00c0 l&rsquo;inverse, les deux observations inject\u00e9es sont formellement isol\u00e9es :<\/p>\n<ol>\n<li>L&rsquo;observation 42 affiche un r\u00e9sidu studentis\u00e9 externe de <span class=\"arabpsy-math-inline\" data-latex=\"t_{42} = -4,12\">$t_{42} = -4,12$<\/span>, associ\u00e9 \u00e0 une probabilit\u00e9 ajust\u00e9e de Bonferroni hautement significative (`bonf(p) = 0,009 &lt; 0,05`). Son levier modeste (<span class=\"arabpsy-math-inline\" data-latex=\"h_{42} = 0,019\">$h_{42} = 0,019$<\/span>, proche de la moyenne th\u00e9orique <span class=\"arabpsy-math-inline\" data-latex=\"p\/n = 3\/150 = 0,020\">$p\/n = 3\/150 = 0,020$<\/span>) confirme qu&#039;il s&#039;agit d&#039;une d\u00e9viation verticale pure n&#039;exer\u00e7ant qu&#039;une influence globale limit\u00e9e (distance de Cook <span class=\"arabpsy-math-inline\" data-latex=\"D_{42} = 0,112\">$D_{42} = 0,112$<\/span>).<\/li>\n<li>L&rsquo;observation 118 pr\u00e9sente une d\u00e9viation encore plus spectaculaire avec un r\u00e9sidu studentis\u00e9 externe de <span class=\"arabpsy-math-inline\" data-latex=\"t_{118} = +5,31\">$t_{118} = +5,31$<\/span> (`bonf(p) &lt; 0,0001`). Conjugu\u00e9 \u00e0 une valeur de levier de <span class=\"arabpsy-math-inline\" data-latex=\"h_{118} = 0,085\">$h_{118} = 0,085$<\/span> (d\u00e9passant le seuil d&#039;alerte classique fix\u00e9 \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"2p\/n = 0,040\">$2p\/n = 0,040$<\/span>), ce r\u00e9sidu extr\u00eame propulse la distance de Cook \u00e0 une valeur massive de <span class=\"arabpsy-math-inline\" data-latex=\"D_{118} = 0,782\">$D_{118} = 0,782$<\/span>, signalant un point influent majeur capable de distordre unilat\u00e9ralement la pente des r\u00e9gresseurs.<\/li>\n<\/ol>\n<p>La consolidation graphique de ces m\u00e9triques via une planche diagnostique \u00e0 quatre quadrants (diagramme r\u00e9sidus studentis\u00e9s contre valeurs ajust\u00e9es, diagramme quantile-quantile de Student, diagramme de dispersion des leviers et cartographie d&rsquo;influence globale) fournit \u00e0 l&rsquo;analyste un panorama exhaustif corroborant de mani\u00e8re univoque les conclusions tabulaires.<\/p>\n<p>### 12.3 R\u00e9daction du paragraphe de r\u00e9sultats aux normes scientifiques<\/p>\n<p>Pour conclure la d\u00e9marche empirique, il convient de formaliser la restitution des analyses conform\u00e9ment aux standards r\u00e9dactionnels et typographiques de la psychom\u00e9trie contemporaine et de la 7e \u00e9dition du manuel de style de l&rsquo;APA.<\/p>\n<p>Le tableau comparatif suivant illustre formellement l&rsquo;impact substantiel de la prise en compte diagnostique sur les estimations param\u00e9triques du mod\u00e8le de r\u00e9gression :<\/p>\n<p>| Param\u00e8tre structurel | Mod\u00e8le 1 : \u00c9chantillon complet (<span class=\"arabpsy-math-inline\" data-latex=\"n = 150\">$n = 150$<\/span>) | Mod\u00e8le 2 : \u00c9chantillon \u00e9pur\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"n = 148\">$n = 148$<\/span>) | Variation relative (<span class=\"arabpsy-math-inline\" data-latex=\"\\Delta %\">$\\Delta %$<\/span>) |<br \/>\n| :&#8212; | :&#8212;: | :&#8212;: | :&#8212;: |<br \/>\n| **Constante (<span class=\"arabpsy-math-inline\" data-latex=\"\\beta_0\">$\\beta_0$<\/span>)** | <span class=\"arabpsy-math-inline\" data-latex=\"52,14 \\quad [46,80,;,57,48]\">$52,14 \\quad [46,80,;,57,48]$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"54,82 \\quad [50,31,;,59,33]\">$54,82 \\quad [50,31,;,59,33]$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"+5,1 %\">$+5,1 %$<\/span> |<br \/>\n| **\u00c9puisement (<span class=\"arabpsy-math-inline\" data-latex=\"\\beta_1\">$\\beta_1$<\/span>)** | <span class=\"arabpsy-math-inline\" data-latex=\"-0,42 \\quad [-0,54,;,-0,30] ; (p &lt; 0,001)\">$-0,42 \\quad [-0,54,;,-0,30] ; (p &lt; 0,001)$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"-0,51 \\quad [-0,61,;,-0,41] ; (p &lt; 0,001)\">$-0,51 \\quad [-0,61,;,-0,41] ; (p &lt; 0,001)$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"+21,4 %\">$+21,4 %$<\/span> |<br \/>\n| **R\u00e9silience (<span class=\"arabpsy-math-inline\" data-latex=\"\\beta_2\">$\\beta_2$<\/span>)** | <span class=\"arabpsy-math-inline\" data-latex=\"+0,38 \\quad [+0,25,;,+0,51] ; (p &lt; 0,001)\">$+0,38 \\quad [+0,25,;,+0,51] ; (p &lt; 0,001)$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"+0,46 \\quad [+0,35,;,+0,57] ; (p &lt; 0,001)\">$+0,46 \\quad [+0,35,;,+0,57] ; (p &lt; 0,001)$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"+21,1 %\">$+21,1 %$<\/span> |<br \/>\n| **Erreur r\u00e9siduelle (<span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>)** | $7,84$ | $6,12$ | <span class=\"arabpsy-math-inline\" data-latex=\"-21,9 %\">$-21,9 %$<\/span> |<br \/>\n| **Coefficient <span class=\"arabpsy-math-inline\" data-latex=\"R^2\">$R^2$<\/span>** | $0,482$ | $0,614$ | <span class=\"arabpsy-math-inline\" data-latex=\"+27,4 %\">$+27,4 %$<\/span> |<br \/>\n| **F-statistique globale** | <span class=\"arabpsy-math-inline\" data-latex=\"68,41 ; (p &lt; 0,001)\">$68,41 ; (p &lt; 0,001)$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"116,84 ; (p &lt; 0,001)\">$116,84 ; (p &lt; 0,001)$<\/span> | <span class=\"arabpsy-math-inline\" data-latex=\"+70,8 %\">$+70,8 %$<\/span> |<\/p>\n<p>Ce comparatif met en lumi\u00e8re l&#039;ampleur de la distorsion caus\u00e9e par la pr\u00e9sence non corrig\u00e9e des observations atypiques : l&#039;estimation de l&#039;ampleur des effets protecteurs de la r\u00e9silience et des impacts d\u00e9l\u00e9t\u00e8res de l&#039;\u00e9puisement se trouvait att\u00e9nu\u00e9e d&#039;environ 21 % dans le mod\u00e8le non \u00e9pur\u00e9, tandis que la variance expliqu\u00e9e globale du mod\u00e8le (<span class=\"arabpsy-math-inline\" data-latex=\"R^2\">$R^2$<\/span>) subissait une d\u00e9gradation artificielle de pr\u00e8s de 13 points de pourcentage en raison du gonflement de la variance r\u00e9siduelle globale <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>.<\/p>\n<p>Le mod\u00e8le r\u00e9dactionnel acad\u00e9mique standard s&#039;\u00e9nonce comme suit :<\/p>\n<p><em>\u00ab Une r\u00e9gression lin\u00e9aire multiple par les moindres carr\u00e9s ordinaires a \u00e9t\u00e9 conduite pour \u00e9valuer la capacit\u00e9 pr\u00e9dictive du niveau d&rsquo;\u00e9puisement \u00e9motionnel et de la r\u00e9silience psychologique sur le bien-\u00eatre subjectif au sein d&rsquo;une cohorte de professionnels de sant\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"n = 150\">$n = 150$<\/span>). L&rsquo;\u00e9valuation pr\u00e9liminaire de la validit\u00e9 interne du mod\u00e8le a repos\u00e9 sur l&rsquo;audit syst\u00e9matique des r\u00e9sidus studentis\u00e9s externes (r\u00e9sidus jackknife) calcul\u00e9s via la routine <code>outlier_test<\/code> de la biblioth\u00e8que <code>statsmodels<\/code> en environnement Python. Le contr\u00f4le du taux d&rsquo;erreur par famille a \u00e9t\u00e9 op\u00e9r\u00e9 au moyen d&rsquo;un ajustement conservateur de Bonferroni au seuil nominal global de <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0,05\">$\\alpha = 0,05$<\/span>.<\/em><\/p>\n<p><em>Cette proc\u00e9dure diagnostique a permis d&rsquo;isoler formellement deux observations pr\u00e9sentant des d\u00e9viations r\u00e9siduelles statistiquement significatives : le participant #42 (<span class=\"arabpsy-math-inline\" data-latex=\"t_{146} = -4,12,; p_{\\text{ajust\u00e9e}} = 0,009\">$t_{146} = -4,12,; p_{\\text{ajust\u00e9e}} = 0,009$<\/span>) et le participant #118 (<span class=\"arabpsy-math-inline\" data-latex=\"t_{146} = +5,31,; p_{\\text{ajust\u00e9e}}  2p\/n,; D_{118} = 0,782\">$t_{146} = +5,31,; p_{\\text{ajust\u00e9e}}  2p\/n,; D_{118} = 0,782$<\/span>). L&rsquo;audit qualitatif des journaux de collecte a corrobor\u00e9 l&rsquo;existence d&rsquo;une r\u00e9ponse non valide li\u00e9e \u00e0 un biais de d\u00e9sirabilit\u00e9 extr\u00eame pour ce participant, tandis que l&rsquo;observation #42 s&rsquo;expliquait par un choc \u00e9motionnel aigu exog\u00e8ne document\u00e9.<\/em><\/p>\n<p><em>Une analyse de sensibilit\u00e9 comparative conduite entre le mod\u00e8le initial (<span class=\"arabpsy-math-inline\" data-latex=\"n = 150,; R^2 = 0,482\">$n = 150,; R^2 = 0,482$<\/span>) et le mod\u00e8le ajust\u00e9 sans les deux observations aberrantes (<span class=\"arabpsy-math-inline\" data-latex=\"n = 148,; R^2 = 0,614\">$n = 148,; R^2 = 0,614$<\/span>) a confirm\u00e9 la robustesse directionnelle des relations observ\u00e9es, tout en r\u00e9v\u00e9lant une sous-estimation substantielle de l&rsquo;intensit\u00e9 de l&rsquo;effet de l&rsquo;\u00e9puisement (<span class=\"arabpsy-math-inline\" data-latex=\"\\Delta \\beta_1 = +21,4%\">$\\Delta \\beta_1 = +21,4%$<\/span>) et de la r\u00e9silience (<span class=\"arabpsy-math-inline\" data-latex=\"\\Delta \\beta_2 = +21,1%\">$\\Delta \\beta_2 = +21,1%$<\/span>) en pr\u00e9sence des cas contaminateurs. Les estimations d\u00e9finitives issues du mod\u00e8le robuste d\u00e9montrent qu&rsquo;une charge d&rsquo;\u00e9puisement \u00e9lev\u00e9e est significativement associ\u00e9e \u00e0 une diminution du bien-\u00eatre subjectif (<span class=\"arabpsy-math-inline\" data-latex=\"\\beta = -0,51,; IC_{95%} = [-0,61,; -0,41],; t(145) = -10,02,; p &lt; 0,001\">$\\beta = -0,51,; IC_{95%} = [-0,61,; -0,41],; t(145) = -10,02,; p &lt; 0,001$<\/span>), tandis qu&#039;une r\u00e9silience \u00e9lev\u00e9e exerce un effet tampon pr\u00e9dictif positif substantiel (<span class=\"arabpsy-math-inline\" data-latex=\"\\beta = +0,46,; IC_{95%} = [+0,35,; +0,57],; t(145) = +8,38,; p &lt; 0,001\">$\\beta = +0,46,; IC_{95%} = [+0,35,; +0,57],; t(145) = +8,38,; p &lt; 0,001$<\/span>). L&#039;ensemble des scripts de simulation et de diagnostic Python a \u00e9t\u00e9 d\u00e9pos\u00e9 en acc\u00e8s ouvert sur le r\u00e9f\u00e9rentiel p\u00e9renne Zenodo conform\u00e9ment aux pr\u00e9conisations de la science ouverte. \u00bb<\/em><\/p>\n<p>En appliquant rigoureusement ce protocole matriciel, algorithmique et r\u00e9dactionnel, l&rsquo;analyste de donn\u00e9es s&rsquo;assure d&rsquo;une mod\u00e9lisation lin\u00e9aire exempte de distorsions artefactuelles, conf\u00e9rant \u00e0 ses inf\u00e9rences statistiques le plus haut standard de rigueur et de cr\u00e9dibilit\u00e9 scientifique.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7th ed.). American Psychological Association. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/li>\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 noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471725153\">https:\/\/doi.org\/10.1002\/0471725153<\/a><\/li>\n<li>Benjamini, Y., &amp; Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. <em>Journal of the Royal Statistical Society: Series B (Methodological)<\/em>, 57(1), 289\u2013300. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2517-6161.1995.tb02031.x\">https:\/\/doi.org\/10.1111\/j.2517-6161.1995.tb02031.x<\/a><\/li>\n<li>Cook, R. D. (1977). Detection of influential observation in linear regression. <em>Technometrics<\/em>, 19(1), 15\u201318. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00401706.1977.10489493\">https:\/\/doi.org\/10.1080\/00401706.1977.10489493<\/a><\/li>\n<li>Cook, R. D., &amp; Weisberg, S. (1982). <em>Residuals and influence in regression<\/em>. Chapman and Hall.<\/li>\n<li>Fox, J. (2015). <em>Applied regression analysis and generalized linear models<\/em> (3rd ed.). SAGE Publications.<\/li>\n<li>Huber, P. J. (1981). <em>Robust statistics<\/em>. John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471725250\">https:\/\/doi.org\/10.1002\/0471725250<\/a><\/li>\n<li>Rousseeuw, P. J., &amp; Leroy, A. M. (1987). <em>Robust regression and outlier detection<\/em>. John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471852333\">https:\/\/doi.org\/10.1002\/0471852333<\/a><\/li>\n<li>Seabold, S., &amp; Perktold, J. (2010). statsmodels: Econometric and statistical modeling with python. In <em>Proceedings of the 9th Python in Science Conference<\/em> (Vol. 57, pp. 92\u201396). <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.25080\/Majora-92bf1922-011\">https:\/\/doi.org\/10.25080\/Majora-92bf1922-011<\/a><\/li>\n<li>Student. (1908). The probable error of a mean. <em>Biometrika<\/em>, 6(1), 1\u201325. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2331554\">https:\/\/doi.org\/10.2307\/2331554<\/a><\/li>\n<li>Weisberg, S. (2014). <em>Applied linear regression<\/em> (4th ed.). John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/9781118594858\">https:\/\/doi.org\/10.1002\/9781118594858<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour calculer et interpr\u00e9ter les r\u00e9sidus studentis\u00e9s en Python avec statsmodels, optimis\u00e9 pour l&rsquo;analyse des donn\u00e9es comportementales.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2465,2992],"tags":[4229,3034,1910,7446,7450],"class_list":["post-5541","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-programmation-python","category-science-des-donnees-en-psychologie","tag-detection-des-valeurs-aberrantes","tag-psychometrie-quantitative","tag-regression-lineaire","tag-residus-studentises","tag-statsmodels-python"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/5541","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=5541"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/5541\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=5541"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=5541"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=5541"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}