{"id":1299,"date":"2026-10-14T18:00:38","date_gmt":"2026-10-14T18:00:38","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-le-r-carre-ajuste-en-python\/"},"modified":"2026-10-14T18:00:38","modified_gmt":"2026-10-14T18:00:38","slug":"comment-calculer-le-r-carre-ajuste-en-python","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-le-r-carre-ajuste-en-python","title":{"rendered":"Comment calculer le R-carr\u00e9 ajust\u00e9 en Python"},"content":{"rendered":"<p>Dans le champ de l&rsquo;analyse statistique moderne et de la mod\u00e9lisation pr\u00e9dictive, l&rsquo;\u00e9valuation de la qualit\u00e9 d&rsquo;ajustement d&rsquo;un mod\u00e8le de r\u00e9gression lin\u00e9aire constitue une \u00e9tape pivotale. Les praticiens des donn\u00e9es, qu&rsquo;ils exercent en \u00e9conom\u00e9trie quantitative, en psychom\u00e9trie comportementale ou en apprentissage automatique, sont syst\u00e9matiquement confront\u00e9s au d\u00e9fi d&rsquo;arbitrer entre la repr\u00e9sentativit\u00e9 empirique d&rsquo;un mod\u00e8le et sa n\u00e9cessaire parcimonie. L&rsquo;indicateur classique d&rsquo;ajustement, le coefficient de d\u00e9termination usuel not\u00e9 <em>R\u00b2<\/em>, a longtemps servi de crit\u00e8re cardinal pour mesurer la proportion de variance de la variable cible expliqu\u00e9e par un faisceau de pr\u00e9dicteurs. N\u00e9anmoins, son utilisation exclusive dans des contextes multivari\u00e9s s&rsquo;av\u00e8re trompeuse en raison d&rsquo;une faille fondamentale : son incapacit\u00e9 intrins\u00e8que \u00e0 p\u00e9naliser la complexit\u00e9 param\u00e9trique superflue.<\/p>\n<p>Pour pallier ce biais m\u00e9canique d&rsquo;inflation qui fausse les inf\u00e9rences scientifiques et favorise le surapprentissage statistique, la litt\u00e9rature \u00e9conom\u00e9trique a vu \u00e9merger une m\u00e9trique corrig\u00e9e : le coefficient de d\u00e9termination ajust\u00e9, ou <strong>R-carr\u00e9 ajust\u00e9<\/strong> (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Coefficient_of_determination#Adjusted_R2\">Adjusted R-squared<\/a>). Formalis\u00e9 originellement par l&rsquo;\u00e9conom\u00e8tre Henri Theil au d\u00e9but des ann\u00e9es 1960, cet indicateur int\u00e8gre explicitement les degr\u00e9s de libert\u00e9 r\u00e9siduels et totaux, r\u00e9\u00e9quilibrant ainsi l&rsquo;apport informatif marginal de chaque nouveau r\u00e9gresseur face au co\u00fbt structurel de son inclusion. D\u00e8s lors, le calcul rigoureux et l&rsquo;interpr\u00e9tation critique du R-carr\u00e9 ajust\u00e9 ne repr\u00e9sentent pas une simple formalit\u00e9 arithm\u00e9tique, mais une exigence m\u00e9thodologique majeure pour garantir la validit\u00e9 interne et la r\u00e9plicabilit\u00e9 externe des conclusions d&rsquo;une \u00e9tude.<\/p>\n<p>L&rsquo;\u00e9cosyst\u00e8me Python s&rsquo;est impos\u00e9 comme l&rsquo;environnement computationnel pr\u00e9dominant pour mener de telles investigations empiriques, gr\u00e2ce \u00e0 la compl\u00e9mentarit\u00e9 de ses biblioth\u00e8ques sp\u00e9cialis\u00e9es telles que <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/numpy.org\">NumPy<\/a>, <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pandas.pydata.org\">Pandas<\/a>, <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/scikit-learn.org\">Scikit-Learn<\/a> et <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.statsmodels.org\">Statsmodels<\/a>. Cet article propose une exploration exhaustive, conceptuelle et pratique du calcul du R-carr\u00e9 ajust\u00e9 en Python. \u00c0 travers un continuum articulant th\u00e9orie des probabilit\u00e9s, d\u00e9rivations d&rsquo;alg\u00e8bre matricielle, codage algorithmique optimis\u00e9 et directives de diffusion acad\u00e9mique, ce guide offre une r\u00e9f\u00e9rence exhaustive destin\u00e9e aux chercheurs, ing\u00e9nieurs et analystes exigeant une rigueur m\u00e9thodologique sans compromis.<\/p>\n<h2>1. Introduction th\u00e9orique : Du coefficient de d\u00e9termination au R-carr\u00e9 ajust\u00e9<\/h2>\n<h3>1.1 D\u00e9finition et r\u00f4le du R-carr\u00e9 classique dans la mod\u00e9lisation lin\u00e9aire<\/h3>\n<p>Le coefficient de d\u00e9termination ordinaire, commun\u00e9ment d\u00e9sign\u00e9 sous le symbole <em>R\u00b2<\/em>, repr\u00e9sente l&rsquo;une des statistiques descriptives les plus r\u00e9pandues dans l&rsquo;analyse de r\u00e9gression lin\u00e9aire univari\u00e9e et multivari\u00e9e. Son r\u00f4le fondamental r\u00e9side dans la quantification de la proportion de la variance totale observ\u00e9e au sein de la variable d\u00e9pendante (ou variable r\u00e9ponse, not\u00e9e conventionnellement <em>y<\/em>) qui se trouve fid\u00e8lement restitu\u00e9e par la combinaison lin\u00e9aire des variables ind\u00e9pendantes ou pr\u00e9dictives (not\u00e9es dans la matrice de conception <em>X<\/em>). Pour appr\u00e9hender formellement cette m\u00e9trique, il convient de d\u00e9composer la variabilit\u00e9 de la s\u00e9rie observ\u00e9e en sommes de carr\u00e9s distinctes.<\/p>\n<p>La variabilit\u00e9 globale des observations autour de leur moyenne arithm\u00e9tique empirique est captur\u00e9e par la somme totale des carr\u00e9s, d\u00e9sign\u00e9e par l&rsquo;acronyme <strong>SST<\/strong> (<em>Total Sum of Squares<\/em>). Alg\u00e9briquement, pour un \u00e9chantillon constitu\u00e9 de <em>n<\/em> observations, cette grandeur s&rsquo;exprime par la sommation pour l&rsquo;indice <em>i<\/em> allant de 1 \u00e0 <em>n<\/em> des \u00e9carts quadratiques entre chaque valeur observ\u00e9e <em>y<sub>i<\/sub><\/em> et la moyenne g\u00e9n\u00e9rale empirique <em>y\u0304<\/em>. Parall\u00e8lement, l&rsquo;\u00e9cart r\u00e9siduel inh\u00e9rent \u00e0 l&rsquo;ajustement du mod\u00e8le lin\u00e9aire, traduisant la fraction d&rsquo;information que les pr\u00e9dicteurs s&rsquo;av\u00e8rent incapables d&rsquo;expliquer, est formalis\u00e9 par la somme des carr\u00e9s des r\u00e9sidus, not\u00e9e <strong>SSR<\/strong> (<em>Sum of Squared Residuals<\/em>) ou parfois <strong>SSE<\/strong> (<em>Sum of Squared Errors<\/em>). Cette derni\u00e8re \u00e9quivaut \u00e0 la sommation sur l&rsquo;\u00e9chantillon des carr\u00e9s des diff\u00e9rences entre les valeurs r\u00e9elles <em>y<sub>i<\/sub><\/em> et les pr\u00e9dictions g\u00e9n\u00e9r\u00e9es par l&rsquo;hyperplan de r\u00e9gression, <em>\u0177<sub>i<\/sub><\/em>.<\/p>\n<p>Dans le cadre strict des moindres carr\u00e9s ordinaires (MCO) int\u00e9grant une constante d&rsquo;ajustement (l&rsquo;intercept), la somme totale des carr\u00e9s se partitionne orthogonalement selon l&rsquo;identit\u00e9 fondamentale de la d\u00e9composition de la variance : SST = SSR + SSM, o\u00f9 SSM (<em>Model Sum of Squares<\/em>) refl\u00e8te la dispersion imputable au mod\u00e8le pr\u00e9dictif. Le coefficient de d\u00e9termination s&rsquo;\u00e9tablit d\u00e8s lors sous la forme du ratio standardis\u00e9 : R\u00b2 = 1 &#8211; (SSR \/ SST). Normalis\u00e9 par construction au sein de l&rsquo;intervalle ferm\u00e9 [0, 1] lorsque l&rsquo;estimation s&rsquo;effectue avec un terme constant, le R-carr\u00e9 propose une \u00e9chelle de mesure intuitive o\u00f9 une valeur de 0 stipule une incapacit\u00e9 totale du mod\u00e8le \u00e0 faire mieux qu&rsquo;une pr\u00e9diction constante \u00e9gale \u00e0 la moyenne, tandis qu&rsquo;une valeur de 1 d\u00e9note un ajustement d\u00e9terministe parfait, l&rsquo;ensemble des r\u00e9sidus \u00e9tant rigoureusement nuls.<\/p>\n<p>N\u00e9anmoins, en d\u00e9pit de son \u00e9l\u00e9gance p\u00e9dagogique, le coefficient de d\u00e9termination classique souffre de limites conceptuelles r\u00e9dhibitoires d\u00e8s lors que l&rsquo;analyste proc\u00e8de \u00e0 des r\u00e9gressions lin\u00e9aires multiples. En effet, sa formulation ne comporte aucun m\u00e9canisme intrins\u00e8que de r\u00e9gulation prenant en compte le nombre de covariables introduites. Cette carence structurelle fait du R-carr\u00e9 une mesure na\u00efve de la performance globale, vuln\u00e9rable \u00e0 la prolif\u00e9ration de variables descriptives dans la matrice explicative.<\/p>\n<h3>1.2 Le biais d&rsquo;inflation syst\u00e9matique du R-carr\u00e9 ordinaire<\/h3>\n<p>Le principal \u00e9cueil m\u00e9thodologique associ\u00e9 au R-carr\u00e9 ordinaire r\u00e9side dans sa propri\u00e9t\u00e9 math\u00e9matique de monotonie non d\u00e9croissante face \u00e0 l&rsquo;adjonction de variables explicatives. Sur le plan g\u00e9om\u00e9trique, l&rsquo;estimation des param\u00e8tres par la m\u00e9thode des moindres carr\u00e9s ordinaires consiste \u00e0 projeter orthogonalement le vecteur d&rsquo;observations <em>y<\/em> sur le sous-espace vectoriel vectoriellement engendr\u00e9 par les colonnes de la matrice des r\u00e9gresseurs <em>X<\/em>. Lorsqu&rsquo;un pr\u00e9dicteur additionnel est int\u00e9gr\u00e9 au mod\u00e8le, la dimension de ce sous-espace vectoriel augmente n\u00e9cessairement, ou demeure au minimum identique si le vecteur ajout\u00e9 s&rsquo;av\u00e8re parfaitement colin\u00e9aire aux pr\u00e9c\u00e9dents.<\/p>\n<p>Par cons\u00e9quent, la distance euclidienne s\u00e9parant le vecteur des valeurs observ\u00e9es du sous-espace projet\u00e9 \u2014 distance dont le carr\u00e9 repr\u00e9sente pr\u00e9cis\u00e9ment la somme des r\u00e9sidus au carr\u00e9 (SSR) \u2014 ne peut math\u00e9matiquement jamais cro\u00eetre lors de l&rsquo;ajout d&rsquo;une nouvelle variable. Dans le cas le plus d\u00e9favorable o\u00f9 le pr\u00e9dicteur introduit ne poss\u00e9derait strictement aucune relation th\u00e9orique ou stochastique avec la variable d&rsquo;int\u00e9r\u00eat, les moindres carr\u00e9s lui affecteront simplement un coefficient proche de z\u00e9ro, profitant des fluctuations \u00e9chantillonnales opportunistes pour r\u00e9duire de mani\u00e8re marginale mais syst\u00e9matique la SSR. En d\u00e9coule un m\u00e9canisme d&rsquo;inflation artificielle continue du R-carr\u00e9, qui augmente irr\u00e9m\u00e9diablement ou stagne, mais ne diminue jamais, ind\u00e9pendamment de la pertinence causale ou statistique des descripteurs ajout\u00e9s.<\/p>\n<p>Ce ph\u00e9nom\u00e8ne d&rsquo;optimisme statistique engendre des cons\u00e9quences s\u00e9v\u00e8res sur le plan inf\u00e9rentiel, d\u00e9sign\u00e9es sous le vocable de <strong>surajustement<\/strong> (<em>overfitting<\/em>). En cherchant \u00e0 maximiser aveugl\u00e9ment le R-carr\u00e9 ordinaire par l&rsquo;accumulation de variables, le mod\u00e9lisateur contraint l&rsquo;\u00e9quation de r\u00e9gression \u00e0 capturer le bruit al\u00e9atoire inh\u00e9rent \u00e0 l&rsquo;\u00e9chantillon d&rsquo;entra\u00eenement plut\u00f4t que les structures g\u00e9n\u00e9ratrices sous-jacentes. La validit\u00e9 externe du mod\u00e8le s&rsquo;effondre alors drastiquement : une \u00e9quation affichant un R\u00b2 spectaculairement \u00e9lev\u00e9 sur ses donn\u00e9es d&rsquo;estimation produira des erreurs de pr\u00e9diction consid\u00e9rables lorsqu&rsquo;elle sera confront\u00e9e \u00e0 de nouvelles observations ind\u00e9pendantes.<\/p>\n<p>Dans les sciences du comportement, en \u00e9pid\u00e9miologie ou en psychom\u00e9trie quantitative, ce biais engendre de faux espoirs explicatifs, poussant les chercheurs \u00e0 multiplier les covariables sans justification \u00e9pist\u00e9mologique sous pr\u00e9texte d&rsquo;accro\u00eetre la variance expliqu\u00e9e. Il devient par cons\u00e9quent indispensable de faire appel \u00e0 une proc\u00e9dure de p\u00e9nalisation th\u00e9orique qui indexe la performance d&rsquo;ajustement sur les degr\u00e9s de libert\u00e9 consomm\u00e9s lors de la param\u00e9trisation.<\/p>\n<h3>1.3 \u00c9mergence et pertinence du R-carr\u00e9 ajust\u00e9<\/h3>\n<p>Face aux insuffisances criantes du coefficient de d\u00e9termination non corrig\u00e9, l&rsquo;\u00e9conom\u00e8tre n\u00e9erlandais Henri Theil a propos\u00e9, dans son ouvrage s\u00e9minal de 1961 intitul\u00e9 <em>Economic Forecasts and Policy<\/em>, une correction syst\u00e9matique calibr\u00e9e sur les degr\u00e9s de libert\u00e9 statistiques. L&rsquo;objectif explicite de cette innovation m\u00e9thodologique \u00e9tait de substituer aux sommes de carr\u00e9s brutes des estimateurs convergents et sans biais des variances de population correspondantes, introduisant ainsi le concept de R-carr\u00e9 ajust\u00e9 (souvent not\u00e9 <em>R\u0304\u00b2<\/em> ou <em>R\u00b2<sub>adj<\/sub><\/em>).<\/p>\n<p>L&rsquo;apport novateur de Theil repose sur la distinction nette entre la variance d&rsquo;\u00e9chantillon et la variance de l&rsquo;al\u00e9a structurel sous-jacent. En ajustant le ratio SSR\/SST par le rapport de leurs degr\u00e9s de libert\u00e9 respectifs, la m\u00e9trique impose un tribut statistique proportionnel \u00e0 la complexit\u00e9 structurelle impos\u00e9e au mod\u00e8le. Chaque param\u00e8tre suppl\u00e9mentaire inclus dans l&rsquo;architecture de r\u00e9gression absorbe un degr\u00e9 de libert\u00e9 r\u00e9siduel, augmentant m\u00e9caniquement le d\u00e9nominateur de p\u00e9nalisation. Ainsi, le R-carr\u00e9 ajust\u00e9 ne progresse que si la r\u00e9duction de la dispersion r\u00e9siduelle induite par le nouveau pr\u00e9dicteur surpasse le co\u00fbt associ\u00e9 \u00e0 la perte de ce degr\u00e9 de libert\u00e9.<\/p>\n<p>Cette formulation conf\u00e8re au R-carr\u00e9 ajust\u00e9 des propri\u00e9t\u00e9s asymptotiques remarquables. Lorsque la taille de l&rsquo;\u00e9chantillon <em>n<\/em> tend vers l&rsquo;infini pour un nombre fixe de r\u00e9gresseurs <em>k<\/em>, le terme de correction des degr\u00e9s de libert\u00e9 tend vers l&rsquo;unit\u00e9, de sorte que l&rsquo;\u00e9cart entre le R\u00b2 ordinaire et le R\u00b2 ajust\u00e9 s&rsquo;annule asymptotiquement. En revanche, pour des \u00e9chantillons finis, particuli\u00e8rement dans les \u00e9tudes cliniques ou les designs quasi-exp\u00e9rimentaux o\u00f9 <em>n<\/em> est modeste et <em>k<\/em> relativement \u00e9lev\u00e9, la sanction impos\u00e9e par le R-carr\u00e9 ajust\u00e9 se fait impitoyable.<\/p>\n<p>Une particularit\u00e9 th\u00e9orique et empirique majeure de cette m\u00e9trique r\u00e9side dans le fait qu&rsquo;elle n&rsquo;est plus strictement confin\u00e9e \u00e0 l&rsquo;intervalle [0, 1]. Si le pouvoir explicatif r\u00e9el des variables est excessivement faible au point que le mod\u00e8le sous-jacent soit moins performant qu&rsquo;une simple assignation de la moyenne, le ratio p\u00e9nalis\u00e9 exc\u00e8de l&rsquo;unit\u00e9, provoquant l&rsquo;apparition d&rsquo;une valeur n\u00e9gative. Loin de constituer une aberration num\u00e9rique, l&rsquo;obtention d&rsquo;un R-carr\u00e9 ajust\u00e9 n\u00e9gatif transmet une information diagnostique essentielle : le mod\u00e8le sp\u00e9cifi\u00e9 est profond\u00e9ment inad\u00e9quat, surajust\u00e9, et ses coefficients n&rsquo;apportent aucun gain informatif cr\u00e9dible au regard des degr\u00e9s de libert\u00e9 gaspill\u00e9s.<\/p>\n<h2>2. Fondements math\u00e9matiques et \u00e9quations de r\u00e9gression<\/h2>\n<h3>2.1 Analyse formelle de la formule du R-carr\u00e9 ajust\u00e9<\/h3>\n<p>La formulation alg\u00e9brique canonique du coefficient de d\u00e9termination ajust\u00e9 proc\u00e8de d&rsquo;une rectification rigoureuse du terme d&rsquo;erreur r\u00e9siduelle standardis\u00e9. Pour un mod\u00e8le de r\u00e9gression lin\u00e9aire multiple comportant <em>n<\/em> observations ind\u00e9pendantes et <em>k<\/em> variables explicatives distinctes (en excluant formellement la constante d&rsquo;ajustement du d\u00e9compte de <em>k<\/em>), la relation math\u00e9matique s&rsquo;\u00e9tablit traditionnellement comme suit :<\/p>\n<p><strong>R\u00b2<sub>adj<\/sub> = 1 &#8211; [ (1 &#8211; R\u00b2) \u00d7 (n &#8211; 1) \/ (n &#8211; k &#8211; 1) ]<\/strong><\/p>\n<p>Une analyse diss\u00e9qu\u00e9e de cette expression r\u00e9v\u00e8le la logique de son architecture. Le terme <em>(1 &#8211; R\u00b2)<\/em> repr\u00e9sente la proportion brute de variance non expliqu\u00e9e par l&rsquo;ajustement des moindres carr\u00e9s, soit exactement le quotient SSR \/ SST. Ce terme scalaire se trouve ensuite multipli\u00e9 par le facteur d&rsquo;\u00e9chelle <em>(n &#8211; 1) \/ (n &#8211; k &#8211; 1)<\/em>, lequel constitue le c\u0153ur du m\u00e9canisme de correction.<\/p>\n<p>Le num\u00e9rateur de ce facteur d&rsquo;\u00e9chelle, <em>(n &#8211; 1)<\/em>, correspond pr\u00e9cis\u00e9ment aux degr\u00e9s de libert\u00e9 associ\u00e9s \u00e0 la somme totale des carr\u00e9s (SST). En effet, l&rsquo;estimation pr\u00e9alable de la moyenne arithm\u00e9tique globale <em>y\u0304<\/em> impose une contrainte lin\u00e9aire unique sur les <em>n<\/em> observations, r\u00e9duisant l&rsquo;espace des \u00e9carts ind\u00e9pendants \u00e0 une dimension de <em>n &#8211; 1<\/em>. Inversement, le d\u00e9nominateur <em>(n &#8211; k &#8211; 1)<\/em> capture les degr\u00e9s de libert\u00e9 de la somme des r\u00e9sidus au carr\u00e9 (SSR). L&rsquo;estimation simultan\u00e9e des <em>k<\/em> pentes partielles (coefficients de r\u00e9gression) ainsi que du terme constant (l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine) pr\u00e9l\u00e8ve au total <em>k + 1<\/em> contraintes lin\u00e9aires distinctes sur le syst\u00e8me d&rsquo;\u00e9quations normales. Par cons\u00e9quent, les r\u00e9sidus sont contraints de r\u00e9sider dans un sous-espace orthogonal de dimension <em>n &#8211; k &#8211; 1<\/em>.<\/p>\n<p>\u00c9tant donn\u00e9 que <em>k \u2265 1<\/em>, le d\u00e9nominateur <em>(n &#8211; k &#8211; 1)<\/em> est rigoureusement inf\u00e9rieur au num\u00e9rateur <em>(n &#8211; 1)<\/em>, ce qui signifie que le ratio multiplicateur <em>(n &#8211; 1) \/ (n &#8211; k &#8211; 1)<\/em> est structurellement sup\u00e9rieur \u00e0 1. Ce ratio amplifie la fraction non expliqu\u00e9e <em>(1 &#8211; R\u00b2)<\/em>, venant par soustraction d\u00e9pr\u00e9cier la valeur finale du R\u00b2 ajust\u00e9 par rapport \u00e0 sa version ordinaire. L&rsquo;ampleur de cette d\u00e9cote est inversement proportionnelle au ratio d&rsquo;\u00e9chantillonnage <em>n \/ k<\/em> : plus le nombre d&rsquo;observations par variable est restreint, plus la p\u00e9nalit\u00e9 s&rsquo;alourdit de fa\u00e7on exponentielle.<\/p>\n<h3>2.2 Comparaison entre variance empirique et variance non biais\u00e9e<\/h3>\n<p>Pour saisir l&rsquo;\u00e9l\u00e9gance th\u00e9orique de la m\u00e9trique de Theil, il est imp\u00e9ratif d&rsquo;examiner le passage des moments empiriques de l&rsquo;\u00e9chantillon aux estimateurs sans biais des param\u00e8tres de la population parente. Dans le formalisme classique de Gauss-Markov, le mod\u00e8le de r\u00e9gression lin\u00e9aire sous-jacent s&rsquo;\u00e9crit matriciellement <em>y = X\u03b2 + \u03b5<\/em>, o\u00f9 le vecteur des perturbations stochastiques inobservables <em>\u03b5<\/em> est suppos\u00e9 homosc\u00e9dastique et non autocorr\u00e9l\u00e9, de variance scalaire <em>\u03c3\u00b2<\/em>.<\/p>\n<p>Si l&rsquo;on \u00e9value na\u00efvement la variance r\u00e9siduelle de l&rsquo;\u00e9chantillon \u00e0 l&rsquo;aide de la formule de la variance empirique brute, on obtient la quantit\u00e9 SSR \/ n. Cependant, il est d\u00e9montr\u00e9 en th\u00e9orie de l&rsquo;estimation que cet estimateur pr\u00e9sente un biais n\u00e9gatif syst\u00e9matique, car les r\u00e9sidus empiriques minimisent par construction la somme quadratique autour des estimations \u00e9chantillonnales plut\u00f4t qu&rsquo;autour de la v\u00e9ritable esp\u00e9rance math\u00e9matique conditionnelle. L&rsquo;estimateur non biais\u00e9 de la variance des perturbations \u2014 qui co\u00efncide avec l&rsquo;<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Mean_squared_error\">erreur quadratique moyenne<\/a> des r\u00e9sidus, not\u00e9e MSE ou <em>s\u00b2<\/em> \u2014 n\u00e9cessite de diviser la somme des carr\u00e9s par ses degr\u00e9s de libert\u00e9 r\u00e9siduels :<\/p>\n<p><strong>s\u00b2 = MSE = SSR \/ (n &#8211; k &#8211; 1)<\/strong><\/p>\n<p>Simultan\u00e9ment, l&rsquo;estimateur sans biais de la variance marginale de la variable d\u00e9pendante <em>y<\/em> n&rsquo;est pas SST \/ n, mais correspond \u00e0 la variance \u00e9chantillonnale corrig\u00e9e :<\/p>\n<p><strong>s\u00b2<sub>y<\/sub> = SST \/ (n &#8211; 1)<\/strong><\/p>\n<p>En reformulant le coefficient de d\u00e9termination comme le compl\u00e9ment \u00e0 l&rsquo;unit\u00e9 du rapport de ces estimateurs non biais\u00e9s de variance, la relation s&rsquo;\u00e9tablit avec clart\u00e9 :<\/p>\n<p><strong>R\u00b2<sub>adj<\/sub> = 1 &#8211; (s\u00b2 \/ s\u00b2<sub>y<\/sub>) = 1 &#8211; [ (SSR \/ (n &#8211; k &#8211; 1)) \/ (SST \/ (n &#8211; 1)) ]<\/strong><\/p>\n<p>Cette d\u00e9rivation d\u00e9montre avec \u00e9clat que le R-carr\u00e9 ajust\u00e9 n&rsquo;est pas un artifice empirique ad hoc, mais l&rsquo;analogue rigoureux du R-carr\u00e9 exprim\u00e9 en termes de ratios de variances non biais\u00e9es de la population. D\u00e8s lors, le comportement de la m\u00e9trique est intrins\u00e8quement sensible au rapport de complexit\u00e9 structurelle <em>n \/ k<\/em>. Si ce ratio est faible, traduisant une saturation param\u00e9trique de l&rsquo;\u00e9chantillon, la variance r\u00e9siduelle non biais\u00e9e s&rsquo;envole, sanctionnant instantan\u00e9ment le pouvoir explicatif revendiqu\u00e9.<\/p>\n<h3>2.3 Crit\u00e8res de d\u00e9cision statistiques bas\u00e9s sur le R-carr\u00e9 ajust\u00e9<\/h3>\n<p>L&rsquo;utilisation du R-carr\u00e9 ajust\u00e9 comme crit\u00e8re de s\u00e9lection et d&rsquo;\u00e9lagage de mod\u00e8les (<em>model selection<\/em>) repose sur des propri\u00e9t\u00e9s math\u00e9matiques \u00e9l\u00e9gantes reliant sa variation marginale aux tests d&rsquo;hypoth\u00e8ses classiques. Consid\u00e9rons un mod\u00e8le de base embo\u00eet\u00e9 contenant <em>k<\/em> r\u00e9gresseurs, auquel un chercheur envisage d&rsquo;adjoindre un pr\u00e9dicteur suppl\u00e9mentaire, portant la dimension \u00e0 <em>k + 1<\/em>.<\/p>\n<p>Contrairement au R\u00b2 ordinaire qui augmentera in\u00e9vitablement, le R\u00b2 ajust\u00e9 n&rsquo;augmentera \u00e0 la suite de cette adjonction que si, et seulement si, l&rsquo;apport d&rsquo;information nouvelle surpasse la perte du degr\u00e9 de libert\u00e9 induit. Une d\u00e9monstration alg\u00e9brique \u00e9tablie de longue date d\u00e9montre que l&rsquo;inclusion d&rsquo;une variable unique am\u00e9liore le R-carr\u00e9 ajust\u00e9 si et seulement si la valeur absolue de la statistique <em>t<\/em> de Student associ\u00e9e au coefficient de cette variable est strictement sup\u00e9rieure \u00e0 l&rsquo;unit\u00e9 : <strong>|t| &gt; 1<\/strong> (ou, de mani\u00e8re \u00e9quivalente, si la statistique de Fisher partielle associ\u00e9e \u00e0 l&rsquo;adjonction v\u00e9rifie <strong>F &gt; 1<\/strong>).<\/p>\n<p>Ce r\u00e9sultat fondamental m\u00e9rite d&rsquo;\u00eatre mis en perspective avec le seuil traditionnel de signification statistique. Pour un seuil de confiance bilat\u00e9ral usuel de 5 % (\u03b1 = 0,05) et pour un \u00e9chantillon de taille raisonnable, la valeur critique de la statistique <em>t<\/em> oscille approximativement autour de 1,96. Il appara\u00eet ainsi que la condition requise pour maximiser le R-carr\u00e9 ajust\u00e9 (|t| &gt; 1, ce qui correspond g\u00e9n\u00e9ralement \u00e0 une valeur de p inf\u00e9rieure \u00e0 environ 0,32) est nettement plus tol\u00e9rante que l&rsquo;exigence formelle de rejet de l&rsquo;hypoth\u00e8se nulle \u00e0 p &lt; 0,05.<\/p>\n<p>Par cons\u00e9quent, retenir un pr\u00e9dicteur sous le seul pr\u00e9texte qu&rsquo;il rehausse marginalement le R\u00b2 ajust\u00e9 ne garantit nullement sa significativit\u00e9 statistique aux normes acad\u00e9miques strictes. N\u00e9anmoins, ce crit\u00e8re constitue un garde-fou imp\u00e9rieux : toute variable dont l&rsquo;introduction engendre une statistique <em>t<\/em> inf\u00e9rieure \u00e0 1 provoque une d\u00e9gradation irr\u00e9versible du R\u00b2 ajust\u00e9, signalant math\u00e9matiquement que son inclusion d\u00e9grade la pr\u00e9cision globale de l&rsquo;estimation de la variance r\u00e9siduelle. L&rsquo;arbitrage entre parcimonie th\u00e9orique et pouvoir descriptif trouve ainsi une assise quantifiable et reproductible.<\/p>\n<h2>3. Pertinence du R-carr\u00e9 ajust\u00e9 dans les recherches quantitatives et psychologiques<\/h2>\n<h3>3.1 Application aux \u00e9chelles psychom\u00e9triques et \u00e9tudes comportementales<\/h3>\n<p>Dans le domaine des sciences psychologiques, de la neuropsychologie et du comportement organisationnel, la mod\u00e9lisation statistique se heurte fr\u00e9quemment \u00e0 des singularit\u00e9s m\u00e9thodologiques complexes. Les variables explicatives employ\u00e9es ne sont pas des mesures physiques parfaites, mais des scores agr\u00e9g\u00e9s issus de questionnaires psychom\u00e9triques (\u00e9chelles de Likert), des temps de r\u00e9action cognitifs ou des indicateurs dimensionnels de traits latents. La tentation est souvent grande d&rsquo;introduire de multiples facettes d&rsquo;un m\u00eame construit afin d&rsquo;expliquer une variable d\u00e9pendante telle que la performance cognitive, la r\u00e9activit\u00e9 \u00e9motionnelle ou l&rsquo;observance th\u00e9rapeutique.<\/p>\n<p>Lorsqu&rsquo;un chercheur int\u00e8gre simultan\u00e9ment plusieurs sous-\u00e9chelles fortement interconnect\u00e9es au sein d&rsquo;une \u00e9quation de r\u00e9gression multiple, le R\u00b2 ordinaire affiche fr\u00e9quemment une progression trompeuse, sugg\u00e9rant \u00e0 tort un gain explicatif substantiel. En r\u00e9alit\u00e9, cette hausse r\u00e9sulte souvent de la redondance informationnelle et du captage opportuniste de la variance d&rsquo;erreur de mesure. Le recours syst\u00e9matique au R-carr\u00e9 ajust\u00e9 permet d&rsquo;assainir m\u00e9thodologiquement ces d\u00e9marches. Si l&rsquo;adjonction d&rsquo;une cinqui\u00e8me facette de la personnalit\u00e9 n&rsquo;apporte aucune validit\u00e9 incr\u00e9mentale r\u00e9elle par rapport aux quatre pr\u00e9c\u00e9dentes, le R\u00b2 ajust\u00e9 stagne ou r\u00e9gresse, informant imm\u00e9diatement le psychom\u00e9tricien de la superfluit\u00e9 de son instrument.<\/p>\n<p>Cette pr\u00e9caution s&rsquo;av\u00e8re particuli\u00e8rement cruciale dans les \u00e9tudes comportementales portant sur des populations cliniques sp\u00e9cifiques \u2014 telles que des patients atteints de pathologies neurod\u00e9g\u00e9n\u00e9ratives rares ou des enfants pr\u00e9sentant des troubles neurod\u00e9veloppementaux s\u00e9v\u00e8res. Dans ces contextes exp\u00e9rimentaux, le recrutement de cohortes massives se heurte \u00e0 des contraintes \u00e9thiques et logistiques insurmontables, limitant fr\u00e9quemment la taille d&rsquo;\u00e9chantillon <em>n<\/em> \u00e0 quelques dizaines de sujets. Appliquer un mod\u00e8le lin\u00e9aire complexe comportant 8 ou 10 pr\u00e9dicteurs comportementaux sur un \u00e9chantillon de 40 participants expose l&rsquo;analyse \u00e0 une perte drastique de degr\u00e9s de libert\u00e9. Seul le R-carr\u00e9 ajust\u00e9 rend compte fid\u00e8lement de la fragilit\u00e9 structurelle du mod\u00e8le, immunisant le clinicien contre des d\u00e9ductions th\u00e9rapeutiques erron\u00e9es bas\u00e9es sur des corr\u00e9lations fortuites.<\/p>\n<h3>3.2 Pr\u00e9vention des faux positifs dans l&rsquo;exploration de donn\u00e9es observationnelles<\/h3>\n<p>La reproductibilit\u00e9 des r\u00e9sultats scientifiques traverse une crise document\u00e9e, exacerb\u00e9e par des pratiques de recherche douteuses connues sous le terme g\u00e9n\u00e9rique de <em>p-hacking<\/em> ou de manipulation s\u00e9lective des sp\u00e9cifications statistiques. Dans les vastes ensembles de donn\u00e9es observationnelles \u2014 tels que les registres \u00e9pid\u00e9miologiques ou les enqu\u00eates sociologiques longitudinales \u2014 les analystes disposent de centaines de covariables d\u00e9mographiques, socio-\u00e9conomiques et comportementales susceptibles de faire office de variables de contr\u00f4le.<\/p>\n<p>La pratique consistant \u00e0 tester une multitude de combinaisons de covariables pour ne rapporter publiquement que celle qui maximise le R-carr\u00e9 ordinaire g\u00e9n\u00e8re une inflation dramatique du taux de faux positifs (erreur de type I). Un mod\u00e8le gonfl\u00e9 par dix variables de contr\u00f4le purement al\u00e9atoires verra invariablement son R\u00b2 progresser, donnant l&rsquo;illusion d&rsquo;une meilleure ma\u00eetrise des facteurs de confusion. \u00c0 l&rsquo;inverse, l&rsquo;examen scrupuleux du R-carr\u00e9 ajust\u00e9 agit comme un r\u00e9v\u00e9lateur d&rsquo;artifices : d\u00e8s lors que les variables de contr\u00f4le ajout\u00e9es n&rsquo;exercent aucune m\u00e9diation th\u00e9orique authentique, le coefficient ajust\u00e9 subit une d\u00e9cote visible, d\u00e9masquant l&rsquo;absence d&rsquo;assise empirique r\u00e9elle du mod\u00e8le exploratoire.<\/p>\n<p>Les comit\u00e9s \u00e9ditoriaux des revues acad\u00e9miques de premier rang en sciences humaines et sociales exigent d\u00e9sormais une transparence totale quant \u00e0 la dynamique s\u00e9quentielle des r\u00e9gressions hi\u00e9rarchiques. Les auteurs doivent pr\u00e9senter la trajectoire conjointe du R\u00b2 ordinaire, du R\u00b2 ajust\u00e9 et de la variation de F (<em>\u0394F<\/em>) \u00e0 chaque \u00e9tape d&rsquo;incorporation de blocs de variables. Cette standardisation garantit que l&rsquo;\u00e9largissement d&rsquo;un paradigme th\u00e9orique n&rsquo;est pas valid\u00e9 sur la base d&rsquo;une simple inflation arithm\u00e9tique de la variance r\u00e9siduelle, consolidant la robustesse des th\u00e9ories psychologiques transmises \u00e0 la communaut\u00e9 scientifique.<\/p>\n<h3>3.3 R\u00e8gles empiriques d&rsquo;\u00e9valuation des seuils de variance en sciences humaines<\/h3>\n<p>L&rsquo;\u00e9valuation quantitative des mod\u00e8les statistiques en sciences humaines soul\u00e8ve la question r\u00e9currente des seuils d&rsquo;acceptabilit\u00e9. Dans les sciences de l&rsquo;ing\u00e9nieur ou en m\u00e9canique c\u00e9leste, des valeurs de R\u00b2 ajust\u00e9 inf\u00e9rieures \u00e0 0,90 sont souvent jug\u00e9es indicatrices d&rsquo;un mod\u00e8le m\u00e9diocre. En revanche, dans l&rsquo;\u00e9tude des comportements humains, caract\u00e9ris\u00e9s par une stochasticit\u00e9 intrins\u00e8que \u00e9lev\u00e9e, une multiplicit\u00e9 de causes non mesur\u00e9es et des bruits d&rsquo;\u00e9chantillonnage incompressibles, les standards d&rsquo;interpr\u00e9tation diff\u00e8rent radicalement.<\/p>\n<p>Les travaux fondateurs de Jacob Cohen (1988) sur la taille d&rsquo;effet statistique, compl\u00e9t\u00e9s par les recommandations de Ferguson (2009) pour les sciences sociales, ont \u00e9tabli des rep\u00e8res d&rsquo;\u00e9valuation empirique qui font autorit\u00e9 :<\/p>\n<ul>\n<li><strong>Effet faible :<\/strong> R\u00b2 ajust\u00e9 avoisinant 0,02 (environ 2 % de variance expliqu\u00e9e). Ce seuil, bien que modeste, peut s&rsquo;av\u00e9rer substantiellement pertinent dans les \u00e9tudes \u00e9pid\u00e9miologiques \u00e0 large \u00e9chelle o\u00f9 une intervention pr\u00e9ventive minimale influence une fraction minime de la population g\u00e9n\u00e9rale.<\/li>\n<li><strong>Effet mod\u00e9r\u00e9 :<\/strong> R\u00b2 ajust\u00e9 compris entre 0,13 et 0,15 (environ 13 % \u00e0 15 % de variance expliqu\u00e9e). Ce palier repr\u00e9sente le niveau d&rsquo;explication typique des mod\u00e8les psychologiques multidimensionnels int\u00e9grant des facteurs cognitifs et motivationnels.<\/li>\n<li><strong>Effet fort :<\/strong> R\u00b2 ajust\u00e9 atteignant ou d\u00e9passant 0,26 \u00e0 0,35. En recherche comportementale, expliquer plus d&rsquo;un tiers de la variance totale constitue une performance explicative de premier plan, rarement d\u00e9pass\u00e9e sans suspecter une redondance tautologique entre pr\u00e9dicteurs et crit\u00e8res.<\/li>\n<\/ul>\n<p>Ces conventions doivent imp\u00e9rativement \u00eatre temp\u00e9r\u00e9es par l&rsquo;erreur de mesure propre aux instruments de recueil. Selon la th\u00e9orie classique des tests, la fid\u00e9lit\u00e9 (not\u00e9e <em>r<sub>xx<\/sub><\/em>) d&rsquo;un questionnaire plafonne math\u00e9matiquement la corr\u00e9lation observable maximale. Par cons\u00e9quent, un R-carr\u00e9 ajust\u00e9 de 0,25 obtenu \u00e0 partir d&rsquo;\u00e9chelles pr\u00e9sentant une consistance interne modeste (alpha de Cronbach d&rsquo;environ 0,70) peut en r\u00e9alit\u00e9 refl\u00e9ter une relation latente quasi-parfaite, dont l&rsquo;expression statistique apparente est att\u00e9nu\u00e9e par l&rsquo;imperfection des m\u00e9triques. Le chercheur doit donc imp\u00e9rativement discriminer la signification purement arithm\u00e9tique du coefficient de son retentissement th\u00e9orique et clinique concret.<\/p>\n<h2>4. Configuration de l&rsquo;environnement Python et d\u00e9pendances logicielles<\/h2>\n<h3>4.1 S\u00e9lection des biblioth\u00e8ques pour l&rsquo;analyse statistique avanc\u00e9e<\/h3>\n<p>L&rsquo;impl\u00e9mentation op\u00e9rationnelle du calcul du R-carr\u00e9 ajust\u00e9 sous Python s&rsquo;appuie sur un socle technologique hautement optimis\u00e9, compos\u00e9 de quatre biblioth\u00e8ques sp\u00e9cialis\u00e9es interconnect\u00e9es, chacune r\u00e9pondant \u00e0 des paradigmes computationnels distincts :<\/p>\n<ul>\n<li><strong>NumPy (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/numpy.org\">Numerical Python<\/a>) :<\/strong> Constitue la fondation du calcul scientifique. Elle procure des structures de tableaux multidimensionnels homog\u00e8nes (<em>ndarray<\/em>) hautement performantes, impl\u00e9ment\u00e9es en C. Son utilisation est primordiale pour vectoriser les calculs d&rsquo;alg\u00e8bre lin\u00e9aire, \u00e9vitant les surco\u00fbts d&rsquo;interpr\u00e9tation des boucles natives de Python lors de l&rsquo;\u00e9valuation des sommes quadratiques matricielles.<\/li>\n<li><strong>Pandas (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/pandas.pydata.org\">Python Data Analysis Library<\/a>) :<\/strong> Fournit les structures <em>DataFrame<\/em> et <em>Series<\/em>, essentielles \u00e0 l&rsquo;ingestion, au nettoyage, \u00e0 l&rsquo;alignement temporel et \u00e0 la manipulation de matrices de donn\u00e9es h\u00e9t\u00e9rog\u00e8nes \u00e9tiquet\u00e9es. Elle simplifie grandement l&rsquo;extraction des dimensions \u00e9chantillonnales et la gestion transparente des m\u00e9tadonn\u00e9es statistiques.<\/li>\n<li><strong>Scikit-Learn (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/scikit-learn.org\">Machine Learning in Python<\/a>) :<\/strong> S&rsquo;impose comme la r\u00e9f\u00e9rence industrielle pour la construction de pipelines de mod\u00e9lisation pr\u00e9dictive. Sa classe <code>LinearRegression<\/code> et son module de m\u00e9triques permettent un entra\u00eenement rapide et une int\u00e9gration standardis\u00e9e au sein de protocoles de validation crois\u00e9e.<\/li>\n<li><strong>Statsmodels (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.statsmodels.org\">Econometric and Statistical Modeling with Python<\/a>) :<\/strong> Indispensable pour l&rsquo;inf\u00e9rence statistique formelle. Contrairement \u00e0 Scikit-Learn qui est orient\u00e9 vers l&rsquo;optimisation pr\u00e9dictive hors-\u00e9chantillon, Statsmodels \u00e9mule l&rsquo;environnement de mod\u00e9lisation du langage R, g\u00e9n\u00e9rant des diagnostics d&rsquo;ad\u00e9quation complets, incluant le calcul natif automatis\u00e9 du R\u00b2 ajust\u00e9, des erreurs-types h\u00e9t\u00e9rosc\u00e9dastiques et des tests d&rsquo;hypoth\u00e8ses sur les param\u00e8tres.<\/li>\n<\/ul>\n<h3>4.2 Gestion des environnements virtuels et compatibilit\u00e9 des versions<\/h3>\n<p>La reproductibilit\u00e9 computationnelle exige une isolation stricte des d\u00e9pendances logicielles via des environnements virtuels d\u00e9di\u00e9s. Les variations de versions sous-jacentes des biblioth\u00e8ques scientifiques (notamment l&rsquo;interface binaire entre SciPy, NumPy et Statsmodels) peuvent modifier le comportement des solveurs matriciels ou la gestion des valeurs singuli\u00e8res. L&rsquo;initialisation d&rsquo;un espace de travail herm\u00e9tique sous Linux, macOS ou Windows s&rsquo;op\u00e8re classiquement \u00e0 l&rsquo;aide de l&rsquo;utilitaire standard <code>venv<\/code> :<\/p>\n<p>L&rsquo;environnement \u00e9tant activ\u00e9, l&rsquo;installation conjointe des d\u00e9pendances cibles dans leurs r\u00e9visions stables s&rsquo;effectue via le gestionnaire de paquets officiel. L&rsquo;ex\u00e9cution de la commande d&rsquo;installation garantit l&rsquo;alignement des liaisons dynamiques avec les biblioth\u00e8ques d&rsquo;acc\u00e9l\u00e9ration mat\u00e9rielle BLAS et LAPACK :<\/p>\n<p><code>pip install numpy pandas scipy scikit-learn statsmodels matplotlib seaborn jupyter<\/code><\/p>\n<p>Il est fortement recommand\u00e9 d&rsquo;utiliser des carnets interactifs (<em>Jupyter Notebooks<\/em>) ou des scripts Python dot\u00e9s d&rsquo;un typage statique rigoureux (via le module standard <code>typing<\/code>). La documentation formelle des signatures de fonctions \u2014 pr\u00e9cisant explicitement les types scalaires et matriciels accept\u00e9s (tels que <code>numpy.ndarray<\/code> ou <code>pandas.DataFrame<\/code>) \u2014 pr\u00e9vient les erreurs silencieuses d&rsquo;alignement dimensionnel qui faussent fr\u00e9quemment le d\u00e9compte des degr\u00e9s de libert\u00e9.<\/p>\n<h3>4.3 Structure modulaire d&rsquo;un script d&rsquo;\u00e9valuation de r\u00e9gression<\/h3>\n<p>Pour s&rsquo;inscrire dans les standards du g\u00e9nie logiciel scientifique, l&rsquo;\u00e9valuation \u00e9conom\u00e9trique ne doit pas \u00eatre ex\u00e9cut\u00e9e au moyen de fragments de scripts monolithiques. Une architecture logicielle p\u00e9renne requiert une s\u00e9gr\u00e9gation nette des responsabilit\u00e9s \u00e0 travers des modules sp\u00e9cialis\u00e9s :<\/p>\n<ol>\n<li><strong>Module d&rsquo;ingestion et validation :<\/strong> Responsable du chargement des tenseurs de donn\u00e9es, de la v\u00e9rification de l&rsquo;int\u00e9grit\u00e9 des types num\u00e9riques, et de la validation conditionnelle stricte des contraintes d&rsquo;identifiabilit\u00e9 du mod\u00e8le (notamment la v\u00e9rification que le nombre d&rsquo;observations <em>n<\/em> est strictement sup\u00e9rieur au nombre de param\u00e8tres <em>k + 1<\/em>).<\/li>\n<li><strong>Module d&rsquo;ajustement et d&rsquo;estimation :<\/strong> Isole l&rsquo;instanciation des algorithmes d&rsquo;estimation matricielle (MCO, moindres carr\u00e9s pond\u00e9r\u00e9s ou r\u00e9gularis\u00e9s) de la couche de pr\u00e9sentation.<\/li>\n<li><strong>Module m\u00e9trique et diagnostique :<\/strong> Encapsule les fonctions pures calculant le R-carr\u00e9 ordinaire, le R-carr\u00e9 ajust\u00e9, l&rsquo;erreur quadratique moyenne et les diagnostics de r\u00e9sidus. Ces fonctions doivent accepter indiff\u00e9remment des structures de donn\u00e9es NumPy ou Pandas sans cr\u00e9er d&rsquo;effets de bord.<\/li>\n<li><strong>Module de journalisation et reporting :<\/strong> G\u00e9n\u00e8re des sorties structur\u00e9es sous formats tabulaires standardis\u00e9s (LaTeX, CSV, JSON), garantissant la tra\u00e7abilit\u00e9 des param\u00e8tres exp\u00e9rimentaux et l&rsquo;archivage p\u00e9renne des r\u00e9sultats d&rsquo;analyse.<\/li>\n<\/ol>\n<h2>5. Pr\u00e9paration et structuration de donn\u00e9es quantitatives pour la mod\u00e9lisation<\/h2>\n<h3>5.1 Chargement et inspection pr\u00e9liminaire des matrices de donn\u00e9es<\/h3>\n<p>La mise en \u0153uvre pratique de toute d\u00e9marche de r\u00e9gression d\u00e9bute par la conversion de fichiers de donn\u00e9es bruts en matrices math\u00e9matiques exploitables. Dans le cadre d&rsquo;une recherche neuropsychologique typique examinant le d\u00e9clin cognitif li\u00e9 \u00e0 l&rsquo;\u00e2ge, consid\u00e9rons l&rsquo;ingestion d&rsquo;un fichier tabulaire mesurant un score d&rsquo;efficience cognitive globale (not\u00e9 <em>y<\/em>) en fonction de pr\u00e9dicteurs tels que l&rsquo;\u00e2ge chronologique, le niveau de r\u00e9serve cognitive (ann\u00e9es d&rsquo;\u00e9ducation formelle), la charge vasculaire syst\u00e9mique et la qualit\u00e9 moyenne du sommeil.<\/p>\n<p>Le chargement s&rsquo;op\u00e8re instantan\u00e9ment via la m\u00e9thode <code>pandas.read_csv()<\/code>. L&rsquo;analyste doit imm\u00e9diatement proc\u00e9der \u00e0 l&rsquo;extraction programmatique des dimensions du syst\u00e8me lin\u00e9aire. L&rsquo;attribut <code>df.shape<\/code> fournit un tuple fondamental <em>(n, m)<\/em>. Si la variable cible <em>y<\/em> est isol\u00e9e sous la forme d&rsquo;un vecteur unidimensionnel de longueur <em>n<\/em>, la matrice explicative <em>X<\/em> comportera une dimension <em>n \u00d7 k<\/em>, o\u00f9 <em>n<\/em> repr\u00e9sente le volume de participants et <em>k<\/em> le nombre exact de pr\u00e9dicteurs initiaux.<\/p>\n<p>Une inspection descriptive pr\u00e9liminaire via <code>df.describe()<\/code> et <code>df.info()<\/code> s&rsquo;impose pour valider les moments empiriques des distributions : moyennes, \u00e9carts-types, asym\u00e9tries (<em>skewness<\/em>) et aplatissements (<em>kurtosis<\/em>). Cette analyse exploratoire bivari\u00e9e permet de s&rsquo;assurer de l&rsquo;absence de valeurs hors normes caus\u00e9es par des erreurs de saisie num\u00e9rique, lesquelles d\u00e9formeraient substantiellement le calcul des sommes quadratiques SST et SSR.<\/p>\n<h3>5.2 Nettoyage, traitement des donn\u00e9es manquantes et imputation<\/h3>\n<p>L&rsquo;omnipr\u00e9sence de donn\u00e9es manquantes (<em>missing data<\/em>) constitue l&rsquo;un des d\u00e9fis les plus redoutables pour l&rsquo;estimation sans biais du R-carr\u00e9 ajust\u00e9. Il est imp\u00e9ratif de diagnostiquer la typologie stochastique r\u00e9gissant ces manques, selon la classification de Rubin : manques compl\u00e8tement al\u00e9atoires (<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Missing_data\">MCAR<\/a>), manques al\u00e9atoires conditionnels (MAR) ou manques non al\u00e9atoires (MNAR).<\/p>\n<p>L&rsquo;approche simpliste consistant \u00e0 appliquer syst\u00e9matiquement une suppression par liste compl\u00e8te (<em>listwise deletion<\/em> via <code>df.dropna()<\/code>) engendre une contraction brutale de la taille d&rsquo;\u00e9chantillon <em>n<\/em>. Or, comme explicit\u00e9 par la formule du R-carr\u00e9 ajust\u00e9, toute amputation de <em>n<\/em> majore directement le ratio de p\u00e9nalit\u00e9 <em>(n &#8211; 1) \/ (n &#8211; k &#8211; 1)<\/em>, d\u00e9gradant m\u00e9caniquement la m\u00e9trique m\u00eame si le pouvoir d&rsquo;ajustement structurel reste constant. De surcro\u00eet, si le m\u00e9canisme de perte n&rsquo;est pas strictement MCAR, cette suppression induit des biais d&rsquo;\u00e9chantillonnage s\u00e9v\u00e8res sur l&rsquo;estimation des coefficients de r\u00e9gression.<\/p>\n<p>Dans les cadres m\u00e9thodologiques contemporains, il est recommand\u00e9 de privil\u00e9gier des techniques d&rsquo;imputation robustes, telles que l&rsquo;imputation it\u00e9rative par cha\u00eenes d&rsquo;\u00e9quations multivari\u00e9es (impl\u00e9ment\u00e9e dans Scikit-Learn via la classe <code>IterativeImputer<\/code>), ou \u00e0 d\u00e9faut une imputation m\u00e9diane conditionnelle pr\u00e9servant les relations structurelles. Dans tous les cas, le nombre final d&rsquo;observations effectives <em>n<\/em> entrant dans la matrice d&rsquo;estimation doit \u00eatre scrupuleusement consign\u00e9, car il conditionne l&rsquo;exactitude arithm\u00e9tique du d\u00e9nominateur des degr\u00e9s de libert\u00e9.<\/p>\n<h3>5.3 Normalisation, centrage et encodage des covariables<\/h3>\n<p>Bien que le calcul matriciel du R-carr\u00e9 au sens strict soit invariant par transformation affine univari\u00e9e des r\u00e9gresseurs (une mise \u00e0 l&rsquo;\u00e9chelle lin\u00e9aire ne modifiant pas le sous-espace vectoriel projet\u00e9), la standardisation pr\u00e9alable des variables continues s&rsquo;av\u00e8re indispensable pour comparer l&rsquo;ampleur relative des coefficients et stabiliser num\u00e9riquement l&rsquo;inversion des matrices mal conditionn\u00e9es. La standardisation en scores Z \u2014 soustrayant la moyenne empirique et divisant par l&rsquo;\u00e9cart-type d&rsquo;\u00e9chantillon \u2014 s&rsquo;effectue ais\u00e9ment \u00e0 l&rsquo;aide de la classe <code>StandardScaler<\/code> de Scikit-Learn.<\/p>\n<p>Lorsque le mod\u00e8le statistique incorpore des covariables cat\u00e9gorielles (par exemple, le statut tabagique ou le g\u00e9notype d&rsquo;un patient), celles-ci doivent faire l&rsquo;objet d&rsquo;un encodage disjonctif complet en variables muettes (<em>dummy encoding<\/em> via <code>pandas.get_dummies(drop_first=True)<\/code> ou <code>OneHotEncoder(drop='first')<\/code> de Scikit-Learn). L&rsquo;omission d\u00e9lib\u00e9r\u00e9e de la modalit\u00e9 de r\u00e9f\u00e9rence (l&rsquo;option <code>drop_first=True<\/code>) est une exigence absolue : son maintien engendrerait une colin\u00e9arit\u00e9 stricte entre les colonnes encod\u00e9es et le terme constant (probl\u00e8me de la \u00ab trappe des variables indicatrices \u00bb), rendant la matrice <em>X&rsquo;X<\/em> non inversible.<\/p>\n<p>Sur le plan m\u00e9thodologique du R-carr\u00e9 ajust\u00e9, il importe de souligner qu&rsquo;une variable qualitative poss\u00e9dant <em>C<\/em> modalit\u00e9s distinctes consomme exactement <em>C &#8211; 1<\/em> degr\u00e9s de libert\u00e9 dans le mod\u00e8le, puisqu&rsquo;elle requiert l&rsquo;estimation de <em>C &#8211; 1<\/em> coefficients ind\u00e9pendants. Par cons\u00e9quent, l&rsquo;inclusion d&rsquo;une variable nominale \u00e0 forte cardinalit\u00e9 \u00e9l\u00e8ve consid\u00e9rablement le param\u00e8tre <em>k<\/em>, imposant une sanction imm\u00e9diate sur le R\u00b2 ajust\u00e9 si ces subdivisions cat\u00e9gorielles n&rsquo;expliquent pas substantiellement la variance de <em>y<\/em>.<\/p>\n<h2>6. Calcul vectoriel direct du R-carr\u00e9 ajust\u00e9 avec NumPy<\/h2>\n<h3>6.1 Extraction des param\u00e8tres fondamentaux de dimensionnalit\u00e9<\/h3>\n<p>Le calcul pur et direct du R-carr\u00e9 ajust\u00e9 via la biblioth\u00e8que NumPy permet d&rsquo;appr\u00e9hender le fonctionnement intime de la m\u00e9trique, sans d\u00e9pendre d&rsquo;abstractions de haut niveau. La premi\u00e8re \u00e9tape algorithmique requiert l&rsquo;extraction sans ambigu\u00eft\u00e9 des grandeurs scalaires fondamentales : le nombre d&rsquo;observations <em>n<\/em> et le nombre de variables explicatives effectives <em>k<\/em>.<\/p>\n<p>Consid\u00e9rons un tableau NumPy \u00e0 deux dimensions <code>X<\/code> repr\u00e9sentant les pr\u00e9dicteurs, et un vecteur unidimensionnel <code>y_true<\/code> contenant les observations cibles. La dimensionnalit\u00e9 s&rsquo;obtient instantan\u00e9ment \u00e0 partir de l&rsquo;attribut <code>shape<\/code> de la matrice :<\/p>\n<p>Le nombre total d&rsquo;observations <em>n<\/em> correspond \u00e0 <code>X.shape[0]<\/code>, tandis que le nombre de r\u00e9gresseurs <em>k<\/em> est fourni par <code>X.shape[1]<\/code>, \u00e0 condition expresse que la matrice <code>X<\/code> ne contienne pas d\u00e9j\u00e0 une colonne unitaire d\u00e9di\u00e9e \u00e0 la constante (l&rsquo;intercept). Si un vecteur colonne de 1 a \u00e9t\u00e9 pr\u00e9alablement concat\u00e9n\u00e9 \u00e0 <code>X<\/code> pour permettre l&rsquo;estimation de l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine, le nombre r\u00e9el de pr\u00e9dicteurs statistiques <em>k<\/em> doit \u00eatre imp\u00e9rativement r\u00e9ajust\u00e9 \u00e0 <code>X.shape[1] - 1<\/code>.<\/p>\n<p>Une validation conditionnelle d&rsquo;int\u00e9grit\u00e9 computationnelle doit imp\u00e9rativement \u00eatre ex\u00e9cut\u00e9e avant tout calcul : il convient de v\u00e9rifier par assertion programmatique que l&rsquo;in\u00e9galit\u00e9 stricte <code>n &gt; k + 1<\/code> est respect\u00e9e. Dans le cas pathologique inverse o\u00f9 <code>n = k + 1<\/code>, le d\u00e9nominateur des degr\u00e9s de libert\u00e9 <em>(n &#8211; k &#8211; 1)<\/em> s&rsquo;annule, provoquant une division par z\u00e9ro qui g\u00e9n\u00e8re un scalaire infini (<code>numpy.inf<\/code>). Si <code>n &lt; k + 1<\/code>, le d\u00e9nominateur devient n\u00e9gatif, inversant absurdeusement le sens de la correction de p\u00e9nalit\u00e9.<\/p>\n<h3>6.2 Impl\u00e9mentation algorithmique de la formule alg\u00e9brique<\/h3>\n<p>L&rsquo;estimation des pr\u00e9dictions <code>y_pred<\/code> \u00e9tant disponible (issue de l&rsquo;application des estimateurs des moindres carr\u00e9s ou de tout algorithme d&rsquo;ajustement lin\u00e9aire), l&rsquo;\u00e9valuation des sommes quadratiques s&rsquo;op\u00e8re par des op\u00e9rations vectoris\u00e9es tirant profit des instructions SIMD des microprocesseurs modernes.<\/p>\n<p>Le vecteur des r\u00e9sidus empiriques est calcul\u00e9 par soustraction \u00e9l\u00e9ment par \u00e9l\u00e9ment : <code>residuals = y_true - y_pred<\/code>. La somme des carr\u00e9s des r\u00e9sidus (SSR) s&rsquo;obtient ensuite par le produit scalaire vectoris\u00e9 ou la fonction de sommation quadratique : <code>ssr = np.sum(residuals ** 2)<\/code>. Parall\u00e8lement, la somme totale des carr\u00e9s (SST) exige le calcul de l&rsquo;\u00e9cart par rapport \u00e0 la moyenne empirique globale : <code>sst = np.sum((y_true - np.mean(y_true)) ** 2)<\/code>.<\/p>\n<p>Dans un premier temps, le R-carr\u00e9 ordinaire est calcul\u00e9 selon son expression canonique : <code>r2 = 1.0 - (ssr \/ sst)<\/code>. Dans un second temps, l&rsquo;application de la correction de Theil s&rsquo;ex\u00e9cute en double pr\u00e9cision (format <code>float64<\/code> natif de NumPy) :<\/p>\n<p><code>adj_r2 = 1.0 - (1.0 - r2) * ((n - 1) \/ (n - k - 1))<\/code><\/p>\n<p>L&rsquo;\u00e9limination totale de structures it\u00e9ratives manuelles (boucles <code>for<\/code> ou compr\u00e9hensions de listes) garantit un temps d&rsquo;ex\u00e9cution quasi-instantan\u00e9, m\u00eame pour des matrices comportant des centaines de milliers de lignes, tout en \u00e9vitant les erreurs d&rsquo;accumulation de troncature num\u00e9rique.<\/p>\n<h3>6.3 Encapsulation dans une fonction Python robuste et typ\u00e9e<\/h3>\n<p>Pour assurer une r\u00e9utilisabilit\u00e9 optimale au sein d&rsquo;environnements de production ou de recherche, ce calcul vectoriel doit \u00eatre encapsul\u00e9 au sein d&rsquo;une fonction logicielle rigoureusement document\u00e9e, int\u00e9grant des annotations de type (<em>type hints<\/em>) et des contr\u00f4les d&rsquo;exceptions stricts :<\/p>\n<p>La fonction peut \u00eatre formalis\u00e9e comme suit : elle re\u00e7oit les param\u00e8tres <code>y_true<\/code> et <code>y_pred<\/code> sous forme de structures <code>numpy.ndarray<\/code> ou d&rsquo;it\u00e9rables convertibles, ainsi que l&rsquo;entier scalaire <code>k<\/code> d\u00e9signant le d\u00e9compte des r\u00e9gresseurs. La fonction v\u00e9rifie au pr\u00e9alable que les vecteurs <code>y_true<\/code> et <code>y_pred<\/code> pr\u00e9sentent un alignement dimensionnel parfait (<code>y_true.shape == y_pred.shape<\/code>) et une dimensionalit\u00e9 unidimensionnelle stricte.<\/p>\n<p>Si la variance empirique de <code>y_true<\/code> est identifi\u00e9e comme rigoureusement nulle (cas d\u00e9g\u00e9n\u00e9r\u00e9 o\u00f9 toutes les valeurs observ\u00e9es sont identiques, rendant SST \u00e9gal \u00e0 z\u00e9ro), la fonction intercepte l&rsquo;erreur et l\u00e8ve une exception statistique explicite (<code>ValueError<\/code>) pour pr\u00e9venir toute division ind\u00e9finie. De m\u00eame, la gestion des cas limites \u2014 lorsque SSR s&rsquo;av\u00e8re sup\u00e9rieur \u00e0 SST, conduisant \u00e0 un R\u00b2 ordinaire n\u00e9gatif \u2014 est trait\u00e9e sans interruption brutale, retournant fid\u00e8lement la valeur d\u00e9cot\u00e9e n\u00e9gative calcul\u00e9e. Cette rigueur algorithmique constitue le fondement indispensable pour int\u00e9grer le calcul du R\u00b2 ajust\u00e9 dans des architectures logicielles plus complexes.<\/p>\n<h2>7. Impl\u00e9mentation avec scikit-learn : Mod\u00e9lisation et calcul personnalis\u00e9<\/h2>\n<h3>7.1 Ajustement du mod\u00e8le LinearRegression et extraction de R-carr\u00e9<\/h3>\n<p>La biblioth\u00e8que Scikit-Learn s&rsquo;est impos\u00e9e comme le standard international pour l&rsquo;ing\u00e9nierie pr\u00e9dictive en Python. Cependant, sa conception philosophique est d\u00e9lib\u00e9r\u00e9ment orient\u00e9e vers l&rsquo;\u00e9valuation de la capacit\u00e9 de g\u00e9n\u00e9ralisation pr\u00e9dictive hors-\u00e9chantillon (<em>out-of-sample<\/em>) plut\u00f4t que vers l&rsquo;inf\u00e9rence \u00e9conom\u00e9trique descriptive in-sample.<\/p>\n<p>L&rsquo;ajustement d&rsquo;un mod\u00e8le lin\u00e9aire s&rsquo;op\u00e8re via la classe <code>sklearn.linear_model.LinearRegression<\/code>. L&rsquo;instanciation de l&rsquo;estimateur suivie de l&rsquo;appel \u00e0 la m\u00e9thode <code>fit(X, y)<\/code> d\u00e9clenche la r\u00e9solution du syst\u00e8me des \u00e9quations normales \u00e0 l&rsquo;aide de la d\u00e9composition en valeurs singuli\u00e8res (SVD) de LAPACK, r\u00e9put\u00e9e pour sa stabilit\u00e9 num\u00e9rique en pr\u00e9sence de quasi-colin\u00e9arit\u00e9 :<\/p>\n<p>Une fois le mod\u00e8le calibr\u00e9, la pr\u00e9diction s&rsquo;obtient via <code>y_pred = model.predict(X)<\/code>. Scikit-Learn fournit la fonction de calcul de score <code>sklearn.metrics.r2_score(y_true, y_pred)<\/code>, laquelle peut \u00e9galement \u00eatre invoqu\u00e9e directement via la m\u00e9thode d&rsquo;\u00e9valuation de l&rsquo;estimateur : <code>model.score(X, y)<\/code>. Cependant, une carence majeure subsiste : Scikit-Learn ne propose nativement aucune fonction ni param\u00e8tre optionnel permettant de calculer directement le R-carr\u00e9 ajust\u00e9.<\/p>\n<h3>7.2 Cr\u00e9ation d&rsquo;une fonction d&rsquo;ajustement compatible avec les API scikit-learn<\/h3>\n<p>Pour combler cette lacune tout en respectant l&rsquo;idiome de conception de Scikit-Learn, il est opportun de d\u00e9velopper une fonction utilitaire wrapper ou une m\u00e9trique personnalis\u00e9e capable d&rsquo;extraire automatiquement la dimensionnalit\u00e9 depuis la matrice de conception <em>X<\/em> :<\/p>\n<p>Cette fonction sur mesure re\u00e7oit en arguments le mod\u00e8le ajust\u00e9 (ou directement les tenseurs <code>y_true<\/code> et <code>y_pred<\/code>) ainsi que la matrice des caract\u00e9ristiques <code>X<\/code>. La fonction extrait automatiquement la taille d&rsquo;\u00e9chantillon <code>n = X.shape[0]<\/code> et le nombre de r\u00e9gresseurs <code>k = X.shape[1]<\/code>. Elle prend \u00e9galement en consid\u00e9ration la pr\u00e9sence effective de l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine en interrogeant l&rsquo;attribut bool\u00e9en <code>model.fit_intercept<\/code>. Si ce dernier est configur\u00e9 sur <code>False<\/code>, le d\u00e9compte des degr\u00e9s de libert\u00e9 de la r\u00e9gression doit \u00eatre adapt\u00e9 : le mod\u00e8le ne consommant pas de degr\u00e9 de libert\u00e9 pour la moyenne, le d\u00e9nominateur de la formule de Theil devient <em>n &#8211; k<\/em> au lieu de <em>n &#8211; k &#8211; 1<\/em>, et le num\u00e9rateur devient <em>n<\/em> au lieu de <em>n &#8211; 1<\/em>.<\/p>\n<p>En int\u00e9grant ce wrapper au sein de l&rsquo;environnement de travail, l&rsquo;analyste r\u00e9concilie la fluidit\u00e9 de l&rsquo;API de Scikit-Learn avec la rigueur diagnostique propre \u00e0 l&rsquo;\u00e9conom\u00e9trie multivari\u00e9e.<\/p>\n<h3>7.3 Int\u00e9gration dans un pipeline de validation de mod\u00e8les multiples<\/h3>\n<p>La v\u00e9ritable puissance de cette impl\u00e9mentation personnalis\u00e9e se manifeste lorsqu&rsquo;elle est combin\u00e9e avec les modules de s\u00e9lection s\u00e9quentielle de caract\u00e9ristiques (<em>feature selection<\/em>) propos\u00e9s par Scikit-Learn, \u00e0 l&rsquo;instar de l&rsquo;\u00e9limination r\u00e9cursive de variables (<strong>RFE<\/strong>, <em>Recursive Feature Elimination<\/em>) ou de la s\u00e9lection univari\u00e9e <code>SelectKBest<\/code>.<\/p>\n<p>Consid\u00e9rons un sc\u00e9nario de recherche dans lequel l&rsquo;exp\u00e9rimentateur dispose de 100 covariables potentielles et souhaite identifier la configuration optimale maximisant le pouvoir explicatif net sans induire de surapprentissage. En encapsulant le mod\u00e8le lin\u00e9aire et le wrapper de R-carr\u00e9 ajust\u00e9 au sein d&rsquo;une boucle it\u00e9rative \u00e9valuant des sous-ensembles croissants de variables (de 1 \u00e0 100), il devient ais\u00e9 de tracer la trajectoire comparative du R\u00b2 classique face au R\u00b2 ajust\u00e9.<\/p>\n<p>Tandis que le R\u00b2 classique cro\u00eet de mani\u00e8re monotone jusqu&rsquo;\u00e0 saturer artificiellement \u00e0 l&rsquo;approche de 100 r\u00e9gresseurs, la courbe du R\u00b2 ajust\u00e9 affiche une dynamique concave typique : elle s&rsquo;\u00e9l\u00e8ve initialement tant que les pr\u00e9dicteurs introduits apportent de l&rsquo;information pr\u00e9dictive authentique, atteint un sommet global (l&rsquo;optimum de parcimonie param\u00e9trique), puis s&rsquo;infl\u00e9chit nettement vers le bas lorsque les variables superflues d\u00e9gradent les degr\u00e9s de libert\u00e9 r\u00e9siduels. Le point culminant de cette courbe fournit ainsi un crit\u00e8re objectif et pragmatique pour stopper l&rsquo;inclusion de covariables.<\/p>\n<h2>8. Automatisation compl\u00e8te avec statsmodels : Extraction native et analyse approfondie<\/h2>\n<h3>8.1 Sp\u00e9cification de mod\u00e8les avec l&rsquo;API OLS de statsmodels<\/h3>\n<p>Pour les chercheurs exigeant un reporting statistique complet conforme aux traditions acad\u00e9miques, la biblioth\u00e8que Statsmodels repr\u00e9sente l&rsquo;outil de pr\u00e9dilection en Python. Contrairement \u00e0 Scikit-Learn, Statsmodels place l&rsquo;inf\u00e9rence th\u00e9orique et la restitution des degr\u00e9s de libert\u00e9 au centre de son architecture.<\/p>\n<p>L&rsquo;utilisation standard de Statsmodels requiert une pr\u00e9caution fondamentale : par d\u00e9faut, la classe <code>statsmodels.api.OLS<\/code> n&rsquo;inclut pas de terme constant dans la matrice de conception. L&rsquo;analyste doit formellement adjoindre une colonne de 1 en invoquant la fonction <code>statsmodels.api.add_constant(X)<\/code>. L&rsquo;ajustement s&rsquo;effectue ensuite par l&rsquo;appel \u00e0 la m\u00e9thode <code>fit()<\/code> :<\/p>\n<p><code>import statsmodels.api as sm<\/code><br \/>\n<code>X_with_const = sm.add_constant(X)<\/code><br \/>\n<code>ols_model = sm.OLS(y, X_with_const).fit()<\/code><\/p>\n<p>D\u00e8s l&rsquo;ach\u00e8vement de l&rsquo;estimation, l&rsquo;objet de r\u00e9sultats <code>RegressionResultsWrapper<\/code> instanci\u00e9 calcule et stocke en m\u00e9moire l&rsquo;ensemble des grandeurs inf\u00e9rentielles. Le R-carr\u00e9 ajust\u00e9 est extrait de mani\u00e8re totalement native et transparente via l&rsquo;attribut d\u00e9di\u00e9 : <code>ols_model.rsquared_adj<\/code>.<\/p>\n<p>Statsmodels propose \u00e9galement une syntaxe d\u00e9clarative par formule calqu\u00e9e sur le standard S\/R, via le module <code>statsmodels.formula.api<\/code>. En \u00e9crivant une \u00e9quation sous la forme <code>'score_cognitif ~ age + education + volume_hippocampe'<\/code>, l&rsquo;ajout du terme constant et l&rsquo;encodage des variables qualitatives sont enti\u00e8rement automatis\u00e9s, et le R\u00b2 ajust\u00e9 demeure imm\u00e9diatement accessible via le m\u00eame attribut <code>rsquared_adj<\/code>.<\/p>\n<h3>8.2 Exploitation du tableau r\u00e9capitulatif summary()<\/h3>\n<p>L&rsquo;un des avantages comp\u00e9titifs majeurs de Statsmodels r\u00e9side dans la g\u00e9n\u00e9ration instantan\u00e9e d&rsquo;un tableau synth\u00e9tique universellement reconnu dans la communaut\u00e9 scientifique via l&rsquo;invocation de <code>ols_model.summary()<\/code>. Ce rapport textuel condense l&rsquo;int\u00e9gralit\u00e9 des diagnostics d&rsquo;ad\u00e9quation en trois blocs distincts :<\/p>\n<p>Le premier bloc r\u00e9capitulatif pr\u00e9sente conjointement la m\u00e9trique <strong>R-squared<\/strong> et la m\u00e9trique <strong>Adj. R-squared<\/strong>. Cette mise en regard directe permet d&rsquo;\u00e9valuer d&rsquo;un simple coup d&rsquo;\u0153il l&rsquo;\u00e9cart de p\u00e9nalit\u00e9 impos\u00e9 au mod\u00e8le. Si l&rsquo;\u00e9cart entre ces deux grandeurs d\u00e9passe quelques centi\u00e8mes, cela signale imm\u00e9diatement une saturation inutile de l&rsquo;espace des variables explicatives face \u00e0 la taille d&rsquo;\u00e9chantillon disponible.<\/p>\n<p>Ce bloc sup\u00e9rieur pr\u00e9sente \u00e9galement la <strong>F-statistic<\/strong> globale et sa significativit\u00e9 associ\u00e9e (<em>Prob (F-statistic)<\/em>). Ce test d&rsquo;hypoth\u00e8se globale \u00e9value l&rsquo;hypoth\u00e8se nulle selon laquelle l&rsquo;ensemble des coefficients de pente <em>\u03b2<sub>1<\/sub>, &#8230;, \u03b2<sub>k<\/sub><\/em> seraient conjointement \u00e9gaux \u00e0 z\u00e9ro. Un mod\u00e8le pr\u00e9sentant un R\u00b2 ajust\u00e9 modeste peut n\u00e9anmoins afficher une statistique F hautement significative si la taille de l&rsquo;\u00e9chantillon est substantielle. Enfin, les crit\u00e8res d&rsquo;information d&rsquo;Akaike (<strong>AIC<\/strong>) et de Schwarz-Bayes (<strong>BIC<\/strong>) compl\u00e8tent utilement l&rsquo;\u00e9valuation : tout comme le R\u00b2 ajust\u00e9, ils p\u00e9nalisent la complexit\u00e9 param\u00e9trique, le BIC imposant une sanction logarithmique proportionnelle \u00e0 <em>ln(n)<\/em> encore plus s\u00e9v\u00e8re pour la s\u00e9lection de mod\u00e8les parcimonieux.<\/p>\n<h3>8.3 Extraction programmatique pour les m\u00e9ta-analyses computationnelles<\/h3>\n<p>Lors de recherches m\u00e9ta-analytiques ou de simulations computationnelles intensives (telles que des \u00e9tudes de Monte-Carlo), l&rsquo;examen visuel des r\u00e9sum\u00e9s textuels devient inapplicable. Il est indispensable d&rsquo;extraire programmatiquement les grandeurs d&rsquo;int\u00e9r\u00eat pour les structurer au sein de bases de donn\u00e9es r\u00e9capitulatives :<\/p>\n<p>L&rsquo;acc\u00e8s direct aux attributs <code>ols_model.rsquared<\/code>, <code>ols_model.rsquared_adj<\/code>, <code>ols_model.fvalue<\/code>, <code>ols_model.f_pvalue<\/code> ainsi qu&rsquo;au vecteur des degr\u00e9s de libert\u00e9 r\u00e9siduels <code>ols_model.df_resid<\/code> permet d&rsquo;it\u00e9rer efficacement sur des centaines de configurations de mod\u00e8les embo\u00eet\u00e9s. Dans le cadre de r\u00e9gressions hi\u00e9rarchiques pas-\u00e0-pas en psychom\u00e9trie, un script Python peut ainsi automatiquement compiler un tableau de comparaison structur\u00e9 sous forme de <code>pandas.DataFrame<\/code>, calculant les variations incr\u00e9mentales de R\u00b2 ajust\u00e9 (<em>\u0394R\u00b2<sub>adj<\/sub><\/em>) entre blocs successifs, puis exporter directement ce tableau au format LaTeX via <code>df.to_latex()<\/code> pour une int\u00e9gration imm\u00e9diate dans un manuscrit acad\u00e9mique.<\/p>\n<h2>9. Comparaison m\u00e9thodique entre R-carr\u00e9 ordinaire et R-carr\u00e9 ajust\u00e9<\/h2>\n<h3>9.1 Simulation exp\u00e9rimentale de l&rsquo;adjonction de bruit blanc<\/h3>\n<p>Pour illustrer sans \u00e9quivoque le comportement divergent des deux coefficients, rien n&rsquo;\u00e9gale une simulation num\u00e9rique contr\u00f4l\u00e9e. Construisons un paradigme exp\u00e9rimental synth\u00e9tique : g\u00e9n\u00e9rons un \u00e9chantillon de taille modeste (<em>n = 100<\/em> observations) d\u00e9pendant causalement de seulement deux variables explicatives authentiques (<em>x<sub>1<\/sub><\/em> et <em>x<sub>2<\/sub><\/em>), reli\u00e9es \u00e0 la variable <em>y<\/em> par l&rsquo;\u00e9quation g\u00e9n\u00e9ratrice <em>y = 3.0 \u00d7 x<sub>1<\/sub> &#8211; 2.0 \u00d7 x<sub>2<\/sub> + \u03b5<\/em>, o\u00f9 <em>\u03b5<\/em> repr\u00e9sente un bruit gaussien de variance unitaire.<\/p>\n<p>Ex\u00e9cutons ensuite un protocole it\u00e9ratif en 50 \u00e9tapes : \u00e0 chaque \u00e9tape, nous adjoignons artificiellement au mod\u00e8le une nouvelle variable pr\u00e9dictive purement al\u00e9atoire, tir\u00e9e ind\u00e9pendamment selon une loi normale centr\u00e9e r\u00e9duite <em>N(0, 1)<\/em>, n&rsquo;ayant strictement aucun lien causal ni stochastique avec <em>y<\/em>. \u00c0 chaque incorporation de ce bruit blanc, nous r\u00e9ajustons le mod\u00e8le lin\u00e9aire par moindres carr\u00e9s ordinaires et enregistrons simultan\u00e9ment la valeur du R\u00b2 ordinaire et celle du R\u00b2 ajust\u00e9.<\/p>\n<p>Les r\u00e9sultats d&rsquo;une telle simulation d\u00e9montrent de mani\u00e8re \u00e9clatante la faille du coefficient classique : alors que les variables ajout\u00e9es ne constituent que du bruit al\u00e9atoire, le R-carr\u00e9 ordinaire s&rsquo;\u00e9l\u00e8ve continuellement, passant par exemple de 0,55 avec les deux vrais pr\u00e9dicteurs \u00e0 plus de 0,75 lorsque 50 bruits blancs ont \u00e9t\u00e9 agr\u00e9g\u00e9s. \u00c0 l&rsquo;oppos\u00e9, le R-carr\u00e9 ajust\u00e9, apr\u00e8s avoir atteint son niveau v\u00e9ridique autour de 0,54, d\u00e9cline progressivement et m\u00e9thodiquement au fur et \u00e0 mesure de l&rsquo;adjonction des pr\u00e9dicteurs parasites, sanctionnant la d\u00e9perdition des degr\u00e9s de libert\u00e9 r\u00e9siduels (qui passent de 97 \u00e0 47) et r\u00e9v\u00e9lant sans fard l&rsquo;ineptie du mod\u00e8le surajust\u00e9.<\/p>\n<h3>9.2 Visualisation graphique sous Matplotlib et Seaborn<\/h3>\n<p>La transposition visuelle de cette divergence constitue un instrument p\u00e9dagogique et scientifique de premier ordre. \u00c0 l&rsquo;aide des biblioth\u00e8ques <code>matplotlib.pyplot<\/code> et <code>seaborn<\/code>, l&rsquo;analyste peut \u00e9laborer des graphiques d&rsquo;une grande rigueur esth\u00e9tique, conformes aux standards graphiques de l&rsquo;American Psychological Association (<strong>APA<\/strong>) :<\/p>\n<p>Le trac\u00e9 superpose deux courbes en fonction du nombre de variables parasites introduites sur l&rsquo;axe des abscisses : une courbe ascendante continue illustrant la hausse illusoire du R\u00b2, et une courbe descendante en pointill\u00e9s mat\u00e9rialisant l&rsquo;\u00e9rosion du R\u00b2 ajust\u00e9. L&rsquo;adjonction de bandes d&rsquo;incertitude obtenues par r\u00e9\u00e9chantillonnage bootstrap ou par r\u00e9p\u00e9tition de 1 000 simulations de Monte-Carlo met en exergue l&rsquo;instabilit\u00e9 grandissante des estimateurs au fur et \u00e0 mesure que l&rsquo;espace param\u00e9trique se sature.<\/p>\n<p>L&rsquo;\u00e9l\u00e9gance du graphique est consolid\u00e9e par l&rsquo;application d&rsquo;un th\u00e8me sobre (<code>seaborn-v0_8-whitegrid<\/code>), des polices avec empattement adapt\u00e9es aux publications universitaires, des \u00e9tiquettes d&rsquo;axes explicitant clairement les notations math\u00e9matiques (<em>R\u00b2<\/em> vs <em>Adjusted R\u00b2<\/em>) et l&rsquo;insertion d&rsquo;une ligne horizontale marquant le niveau th\u00e9orique optimal de variance explicable. Une telle figure d\u00e9montre visuellement l&rsquo;imp\u00e9rieuse n\u00e9cessit\u00e9 du coefficient ajust\u00e9 pour tout arbitrage structurel.<\/p>\n<h3>9.3 Analyse du comportement asymptotique selon la taille d&rsquo;\u00e9chantillon<\/h3>\n<p>L&rsquo;\u00e9cart observ\u00e9 entre le R\u00b2 conventionnel et sa variante ajust\u00e9e est fondamentalement modul\u00e9 par la taille absolue de l&rsquo;\u00e9chantillon <em>n<\/em>. L&rsquo;examen attentif du terme de correction <em>(n &#8211; 1) \/ (n &#8211; k &#8211; 1)<\/em> met en lumi\u00e8re les conditions r\u00e9gissant la convergence asymptotique des deux m\u00e9triques.<\/p>\n<p>Consid\u00e9rons le contexte contemporain des m\u00e9gadonn\u00e9es (<em>Big Data<\/em>), o\u00f9 les jeux de donn\u00e9es \u00e9pid\u00e9miologiques ou g\u00e9nomiques agr\u00e8gent des centaines de milliers d&rsquo;observations (par exemple, <em>n = 500 000<\/em> participants dans la UK Biobank) pour un ensemble mod\u00e9r\u00e9 de covariables d\u00e9mographiques (<em>k = 20<\/em>). Dans cette configuration, le ratio multiplicateur \u00e9quivaut \u00e0 <em>499 999 \/ 499 979 \u2248 1,00004<\/em>. La p\u00e9nalit\u00e9 devient alors infinit\u00e9simale : le R\u00b2 ordinaire et le R\u00b2 ajust\u00e9 co\u00efncident jusqu&rsquo;\u00e0 la quatri\u00e8me ou cinqui\u00e8me d\u00e9cimale. Dans ces conditions asymptotiques, la m\u00e9trique ajust\u00e9e perd de son utilit\u00e9 diagnostique imm\u00e9diate, car les degr\u00e9s de libert\u00e9 r\u00e9siduels sont pl\u00e9thoriques.<\/p>\n<p>\u00c0 l&rsquo;extr\u00eame inverse, dans les recherches neuropsychologiques ou psychiatriques sp\u00e9cialis\u00e9es impliquant des cohortes cliniques rares (par exemple, <em>n = 25<\/em> patients victimes d&rsquo;un accident vasculaire c\u00e9r\u00e9bral focal) test\u00e9es sur <em>k = 6<\/em> \u00e9preuves cognitives, le facteur d&rsquo;ajustement devient <em>(24) \/ (18) = 1,333<\/em>. La fraction non expliqu\u00e9e <em>(1 &#8211; R\u00b2)<\/em> se trouve amplifi\u00e9e d&rsquo;un tiers. Dans ce sc\u00e9nario de faible puissance statistique, l&rsquo;\u00e9cart entre les deux indices est b\u00e9ant : un mod\u00e8le affichant un R\u00b2 ordinaire faussement flatteur de 0,40 s&rsquo;effondrera \u00e0 un R\u00b2 ajust\u00e9 r\u00e9el de 0,20. Les directives m\u00e9thodologiques imposent donc d&rsquo;accorder une vigilance d\u00e9cupl\u00e9e au R-carr\u00e9 ajust\u00e9 d\u00e8s lors que le ratio d&rsquo;\u00e9chantillonnage <em>n \/ k<\/em> descend en de\u00e7\u00e0 d&rsquo;un seuil critique estim\u00e9 empiriquement \u00e0 20 ou 30 observations par r\u00e9gresseur.<\/p>\n<h2>10. Diagnostics de r\u00e9gression, multicolin\u00e9arit\u00e9 et impact sur le R-carr\u00e9 ajust\u00e9<\/h2>\n<h3>10.1 Influence de la multicolin\u00e9arit\u00e9 sur l&rsquo;instabilit\u00e9 des estimations<\/h3>\n<p>La multicolin\u00e9arit\u00e9 d\u00e9signe la pr\u00e9sence de fortes interd\u00e9pendances lin\u00e9aires entre deux ou plusieurs variables ind\u00e9pendantes au sein de la matrice de conception <em>X<\/em>. Si la colin\u00e9arit\u00e9 exacte rend l&rsquo;inversion de la matrice <em>X&rsquo;X<\/em> rigoureusement impossible, la quasi-colin\u00e9arit\u00e9 \u2014 fr\u00e9quente dans les batteries de tests psychom\u00e9triques o\u00f9 plusieurs scores \u00e9valuent des composantes similaires de l&rsquo;anxi\u00e9t\u00e9 ou de la m\u00e9moire de travail \u2014 engendre des distorsions diagnostiques majeures.<\/p>\n<p>Sur le plan computationnel, une multicolin\u00e9arit\u00e9 s\u00e9v\u00e8re n&rsquo;alt\u00e8re pas directement la capacit\u00e9 globale du mod\u00e8le \u00e0 s&rsquo;ajuster aux donn\u00e9es d&rsquo;\u00e9chantillon : le R\u00b2 et le R\u00b2 ajust\u00e9 conservent leur validit\u00e9 arithm\u00e9tique globale. Cependant, elle provoque une inflation colossale de la variance des estimateurs des coefficients de r\u00e9gression individuels (les pentes <em>\u03b2<sub>j<\/sub><\/em>), \u00e9largissant consid\u00e9rablement leurs intervalles de confiance et provoquant l&rsquo;effondrement des statistiques <em>t<\/em> associ\u00e9es. Le mod\u00e9lisateur peut ainsi \u00eatre confront\u00e9 \u00e0 un paradoxe d\u00e9routant : un mod\u00e8le affichant un R-carr\u00e9 ajust\u00e9 tr\u00e8s \u00e9lev\u00e9 (par exemple 0,65) accompagn\u00e9 d&rsquo;une statistique F globale hautement significative, tandis qu&rsquo;aucun des pr\u00e9dicteurs individuels ne franchit le seuil de significativit\u00e9 marginale <em>p &lt; 0,05<\/em>.<\/p>\n<p>L&rsquo;\u00e9valuation pratique de cette pathologie statistique sous Python s&rsquo;op\u00e8re par le calcul du Facteur d&rsquo;Inflation de la Variance (<strong>VIF<\/strong>, <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Variance_inflation_factor\">Variance Inflation Factor<\/a>). Le module <code>statsmodels.stats.outliers_influence<\/code> fournit la fonction <code>variance_inflation_factor<\/code>. Un VIF calcul\u00e9 pour un pr\u00e9dicteur <em>j<\/em> \u00e9value le R\u00b2 issu de la r\u00e9gression de ce pr\u00e9dicteur sur tous les autres r\u00e9gresseurs : <em>VIF<sub>j<\/sub> = 1 \/ (1 &#8211; R\u00b2<sub>j<\/sub>)<\/em>. Les r\u00e8gles empiriques stipulent qu&rsquo;un VIF exc\u00e9dant 5 ou 10 traduit une redondance critique. L&rsquo;\u00e9limination des variables affichant les VIF les plus \u00e9lev\u00e9s assainit le mod\u00e8le : bien que cette suppression puisse abaisser marginalement le R\u00b2 ordinaire, elle stabilise les erreurs-types et se traduit fr\u00e9quemment par un maintien, voire une progression du R\u00b2 ajust\u00e9 gr\u00e2ce \u00e0 l&rsquo;\u00e9conomie de degr\u00e9s de libert\u00e9 r\u00e9alis\u00e9e.<\/p>\n<h3>10.2 Impact des points aberrants et leviers statistiques<\/h3>\n<p>Le calcul du R-carr\u00e9 ajust\u00e9 repose int\u00e9gralement sur la m\u00e9thode des moindres carr\u00e9s, laquelle est notoire pour son hypersensibilit\u00e9 aux observations extr\u00eames. Puisque la fonction de perte quadratique \u00e9l\u00e8ve au carr\u00e9 l&rsquo;ensemble des \u00e9carts r\u00e9siduels, une observation aberrante isol\u00e9e (<em>outlier<\/em>) dans la variable d\u00e9pendante ou un point pr\u00e9sentant un fort levier dans l&rsquo;espace des pr\u00e9dicteurs (<em>leverage point<\/em>) peut d\u00e9former disproportionnellement la surface de r\u00e9gression.<\/p>\n<p>L&rsquo;impact d&rsquo;une telle observation atypique sur le R\u00b2 ajust\u00e9 peut s&rsquo;av\u00e9rer destructeur ou, au contraire, trompeusement amplificateur :<\/p>\n<ul>\n<li>Si une observation aberrante se situe \u00e0 l&rsquo;\u00e9cart de la trajectoire g\u00e9n\u00e9rale du nuage de points, son r\u00e9sidu quadratique sera colossal, faisant exploser la somme des r\u00e9sidus (SSR). Le R\u00b2 ajust\u00e9 s&rsquo;effondrera instantan\u00e9ment, masquant une relation lin\u00e9aire pourtant hautement fid\u00e8le sur les 99 % restants de l&rsquo;\u00e9chantillon.<\/li>\n<li>Inversement, si un point extr\u00eame cumule un fort levier et s&rsquo;aligne opportun\u00e9ment dans la direction de la pente estim\u00e9e, il tirera l&rsquo;hyperplan vers lui, gonflant artificiellement la somme totale des carr\u00e9s (SST) et comprimant les r\u00e9sidus relatifs, ce qui propulsera le R\u00b2 ajust\u00e9 \u00e0 des niveaux fictifs.<\/li>\n<\/ul>\n<p>Sous Python, la d\u00e9tection rigoureuse de ces points s&rsquo;effectue via l&rsquo;extraction de la distance de Cook et des valeurs de levier (\u00e9l\u00e9ments diagonaux de la matrice chapeau <em>H = X(X&rsquo;X)\u207b\u00b9X&rsquo;<\/em>) accessibles directement sur l&rsquo;objet d&rsquo;ajustement Statsmodels : <code>influence = ols_model.get_influence()<\/code>, puis <code>cooks_d = influence.cooks_distance[0]<\/code>. Face \u00e0 des distributions asym\u00e9triques ou contamin\u00e9es, l&rsquo;analyste se doit de comparer le R\u00b2 ajust\u00e9 des moindres carr\u00e9s \u00e0 des m\u00e9triques robustes issues de mod\u00e8les de r\u00e9gression robuste (impl\u00e9ment\u00e9s dans <code>statsmodels.formula.api.rlm<\/code> \u00e0 l&rsquo;aide de fonctions de perte de Huber ou de Tukey), garantissant une \u00e9valuation immunis\u00e9e contre les artefacts d&rsquo;\u00e9chantillonnage.<\/p>\n<h3>10.3 Validation des hypoth\u00e8ses classiques de Gauss-Markov<\/h3>\n<p>L&rsquo;interpr\u00e9tation inf\u00e9rentielle du R-carr\u00e9 ajust\u00e9 s&rsquo;ancre dans le respect strict des hypoth\u00e8ses du th\u00e9or\u00e8me de Gauss-Markov. Si la formule arithm\u00e9tique de Theil peut toujours \u00eatre math\u00e9matiquement calcul\u00e9e sur n&rsquo;importe quel tableau de nombres, sa validit\u00e9 en tant qu&rsquo;estimateur sans biais du ratio de variance sous-jacent s&rsquo;effondre lorsque ces postulats sont transgress\u00e9s.<\/p>\n<p>La premi\u00e8re hypoth\u00e8se cardinale concerne l&rsquo;<strong>homosc\u00e9dasticit\u00e9<\/strong> : la variance des r\u00e9sidus doit demeurer constante le long de toutes les valeurs pr\u00e9dictives. En pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 (fr\u00e9quente dans les \u00e9tudes \u00e9conomiques o\u00f9 la variance des d\u00e9penses s&rsquo;accro\u00eet proportionnellement au revenu), l&rsquo;estimateur MCO de la variance r\u00e9siduelle devient inefficace et les erreurs-types sont biais\u00e9es. Sous Python, la d\u00e9tection de cette faille s&rsquo;ex\u00e9cute via les tests formels de Breusch-Pagan ou de White disponibles dans <code>statsmodels.stats.diagnostic.het_breuschpagan<\/code>. En cas de violation av\u00e9r\u00e9e, le recours \u00e0 des estimateurs de covariance robustes \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 (matrices sandwich HC0 \u00e0 HC3) s&rsquo;impose, alt\u00e9rant la signification inf\u00e9rentielle des statistiques rattach\u00e9es au R\u00b2 ajust\u00e9.<\/p>\n<p>La seconde exigence majeure implique l&rsquo;<strong>ind\u00e9pendance des erreurs<\/strong> stochastiques. Dans les donn\u00e9es longitudinales, chronologiques ou spatialement agr\u00e9g\u00e9es, les r\u00e9sidus pr\u00e9sentent fr\u00e9quemment une autocorr\u00e9lation s\u00e9rielle. Le test de Durbin-Watson \u2014 dont la statistique est directement restitu\u00e9e dans le tableau de synth\u00e8se <code>summary()<\/code> de Statsmodels \u2014 permet de diagnostiquer ce biais (une valeur proche de 2 traduisant l&rsquo;absence d&rsquo;autocorr\u00e9lation). L&rsquo;autocorr\u00e9lation positive conduit \u00e0 une sous-estimation syst\u00e9matique de la v\u00e9ritable variance r\u00e9siduelle, provoquant une sur\u00e9valuation spectaculaire et artificielle du R-carr\u00e9 ajust\u00e9. Enfin, la normalit\u00e9 des r\u00e9sidus \u2014 indispensable pour valider la distribution des statistiques F associ\u00e9es au R\u00b2 ajust\u00e9 en petits \u00e9chantillons \u2014 doit \u00eatre attest\u00e9e par les tests de Shapiro-Wilk (<code>scipy.stats.shapiro<\/code>) ou le test de Jarque-Bera rapport\u00e9 par Statsmodels.<\/p>\n<h2>11. Validation crois\u00e9e et estimation out-of-sample du R-carr\u00e9<\/h2>\n<h3>11.1 Le R-carr\u00e9 ajust\u00e9 face aux m\u00e9triques de g\u00e9n\u00e9ralisation hors-\u00e9chantillon<\/h3>\n<p>Une distinction \u00e9pist\u00e9mologique fondamentale doit \u00eatre op\u00e9r\u00e9e entre l&rsquo;ad\u00e9quation d&rsquo;un mod\u00e8le sur ses propres donn\u00e9es d&rsquo;entra\u00eenement (<em>in-sample fit<\/em>) et son aptitude \u00e0 formuler des pr\u00e9dictions exactes sur des observations futures non mobilis\u00e9es lors de l&rsquo;estimation (<em>out-of-sample generalization<\/em>). Bien que le R-carr\u00e9 ajust\u00e9 introduise une p\u00e9nalit\u00e9 th\u00e9orique fond\u00e9e sur les degr\u00e9s de libert\u00e9 pour mitiger l&rsquo;optimisme des moindres carr\u00e9s, il demeure par d\u00e9finition une m\u00e9trique <strong>in-sample<\/strong>.<\/p>\n<p>Le R\u00b2 ajust\u00e9 suppose implicitement que les conditions d&rsquo;\u00e9chantillonnage sont repr\u00e9sentatives et que la sp\u00e9cification lin\u00e9aire est structurellement exacte. Or, dans les contextes r\u00e9els, des non-lin\u00e9arit\u00e9s sous-jacentes m\u00e9connues, des ruptures structurelles ou des biais de s\u00e9lection de covariables peuvent conf\u00e9rer \u00e0 un mod\u00e8le un R-carr\u00e9 ajust\u00e9 in-sample remarquable, alors que ses pr\u00e9dictions appliqu\u00e9es \u00e0 une cohorte de validation externe r\u00e9v\u00e9leront une erreur quadratique moyenne de pr\u00e9diction (<strong>RMSEP<\/strong>) d\u00e9sastreuse.<\/p>\n<p>Dans la perspective de l&rsquo;apprentissage automatique, le dilemme biais-variance ne se r\u00e9sout pas exclusivement par des ajustements alg\u00e9briques des degr\u00e9s de libert\u00e9, mais par une confrontation directe avec des donn\u00e9es non vues. Il est donc m\u00e9thodologiquement imprudent d&rsquo;\u00e9riger le R\u00b2 ajust\u00e9 en indicateur unique de performance pr\u00e9dictive sans l&rsquo;adosser \u00e0 des protocoles de partitionnement rigoureux des donn\u00e9es.<\/p>\n<h3>11.2 Calcul du R-carr\u00e9 par validation crois\u00e9e K-Fold avec Python<\/h3>\n<p>Pour mesurer la v\u00e9ritable capacit\u00e9 d&rsquo;un mod\u00e8le \u00e0 expliquer la variance hors-\u00e9chantillon, la proc\u00e9dure \u00e9talon r\u00e9side dans la validation crois\u00e9e \u00e0 <em>k<\/em> blocs (<strong>K-Fold Cross-Validation<\/strong>). Dans ce paradigme, le jeu de donn\u00e9es initial est scind\u00e9 en <em>K<\/em> partitions mutuellement exclusives (par exemple <em>K = 5<\/em> ou <em>K = 10<\/em>). \u00c0 tour de r\u00f4le, chaque pli est r\u00e9serv\u00e9 pour le test tandis que les <em>K &#8211; 1<\/em> plis restants servent \u00e0 ajuster les param\u00e8tres de r\u00e9gression.<\/p>\n<p>L&rsquo;impl\u00e9mentation sous Python est grandement facilit\u00e9e par les classes <code>sklearn.model_selection.KFold<\/code> et <code>sklearn.model_selection.cross_val_score<\/code>. Sur chaque ensemble de test ind\u00e9pendant, le score R-carr\u00e9 calcul\u00e9 refl\u00e8te fid\u00e8lement la proportion de variance r\u00e9ellement captur\u00e9e hors-\u00e9chantillon :<\/p>\n<p>Contrairement au R-carr\u00e9 ordinaire in-sample qui est born\u00e9 inf\u00e9rieurement par 0, le R-carr\u00e9 calcul\u00e9 par validation crois\u00e9e hors-\u00e9chantillon (souvent d\u00e9sign\u00e9 sous l&rsquo;appellation de <em>Q\u00b2<\/em> ou R\u00b2 pr\u00e9dictif) peut ais\u00e9ment plonger dans des valeurs n\u00e9gatives si l&rsquo;hyperplan ajust\u00e9 g\u00e9n\u00e9ralise moins bien que la simple moyenne arithm\u00e9tique de la cohorte d&rsquo;entra\u00eenement. La comparaison m\u00e9thodique entre le R\u00b2 ajust\u00e9 in-sample et le score moyen de validation crois\u00e9e hors-\u00e9chantillon constitue le test de robustesse ultime : un mod\u00e8le authentiquement parcimonieux et stable doit afficher une convergence \u00e9troite entre ces deux m\u00e9triques.<\/p>\n<h3>11.3 Techniques de r\u00e9\u00e9chantillonnage bootstrap pour des intervalles de confiance<\/h3>\n<p>L&rsquo;une des limites historiques du R-carr\u00e9 ajust\u00e9 dans les logiciels statistiques conventionnels est l&rsquo;absence de fourniture syst\u00e9matique d&rsquo;un intervalle de confiance. En raison de la complexit\u00e9 de la distribution d&rsquo;\u00e9chantillonnage analytique d&rsquo;un ratio de formes quadratiques sous des distributions non parfaitement gaussiennes, d\u00e9river math\u00e9matiquement l&rsquo;erreur-type exacte du R\u00b2 ajust\u00e9 est une t\u00e2che analytique ardue.<\/p>\n<p>La puissance computationnelle contemporaine permet de contourner cette contrainte gr\u00e2ce au r\u00e9\u00e9chantillonnage stochastique par <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Bootstrapping_(statistics)\">bootstrap<\/a> non param\u00e9trique. En tirant de mani\u00e8re r\u00e9p\u00e9t\u00e9e (par exemple <em>B = 2 000<\/em> it\u00e9rations) des sous-\u00e9chantillons de taille <em>n<\/em> avec remise au sein du jeu de donn\u00e9es original, on ajuste le mod\u00e8le de r\u00e9gression et l&rsquo;on calcule le R-carr\u00e9 ajust\u00e9 sur chaque r\u00e9plique.<\/p>\n<p>La distribution empirique r\u00e9sultante des 2 000 valeurs de R\u00b2 ajust\u00e9 permet d&rsquo;appr\u00e9hender directement son asym\u00e9trie d&rsquo;\u00e9chantillonnage et d&rsquo;en d\u00e9duire des intervalles de confiance non param\u00e9triques par la m\u00e9thode des percentiles (par exemple, en retenant les 2,5<sup>e<\/sup> et 97,5<sup>e<\/sup> percentiles pour un intervalle \u00e0 95 %). En neuropsychologie ou en m\u00e9decine personnalis\u00e9e, le reporting de cet intervalle de confiance bootstrap\u00e9 apporte une plus-value scientifique d\u00e9terminante : un R\u00b2 ajust\u00e9 ponctuel de 0,30 dont l&rsquo;intervalle \u00e0 95 % s&rsquo;\u00e9tend de 0,04 \u00e0 0,55 informe imm\u00e9diatement la communaut\u00e9 acad\u00e9mique de la forte volatilit\u00e9 de l&rsquo;estimation, incitant \u00e0 une mod\u00e9ration l\u00e9gitime des conclusions cliniques.<\/p>\n<h2>12. Bonnes pratiques de reporting acad\u00e9mique et diffusion des r\u00e9sultats<\/h2>\n<h3>12.1 Mise en forme des r\u00e9sultats selon les normes APA (7e \u00e9dition)<\/h3>\n<p>La diffusion publique de mod\u00e8les de r\u00e9gression au sein de revues scientifiques \u00e0 comit\u00e9 de lecture est assujettie \u00e0 des r\u00e8gles typographiques et m\u00e9thodologiques strictes, notamment codifi\u00e9es par le manuel de publication de l&rsquo;<a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\">American Psychological Association<\/a> (APA, 7<sup>e<\/sup> \u00e9dition). Les prescriptions exigent une transparence sans \u00e9quivoque, proscrivant l&rsquo;annonce isol\u00e9e d&rsquo;un coefficient d&rsquo;ajustement sans son armature inf\u00e9rentielle.<\/p>\n<p>Dans le corps du texte, la formulation canonique requiert le signalement conjoint du mod\u00e8le global, des degr\u00e9s de libert\u00e9, de la valeur du test F, de sa significativit\u00e9 probabiliste exacte, ainsi que des deux indicateurs de d\u00e9termination :<\/p>\n<p><em>\u00ab Une r\u00e9gression lin\u00e9aire multiple a \u00e9t\u00e9 r\u00e9alis\u00e9e pour pr\u00e9dire le score d&rsquo;efficience cognitive en fonction de l&rsquo;\u00e2ge, du niveau d&rsquo;\u00e9ducation et du volume hippocampique. L&rsquo;ensemble des pr\u00e9dicteurs explique une part statistiquement significative de la variance, F(3, 96) = 14,82, p &lt; ,001. Le mod\u00e8le affiche un coefficient de d\u00e9termination R\u00b2 = ,316, correspondant \u00e0 un R\u00b2 ajust\u00e9 de ,295, ce qui traduit une taille d&rsquo;effet mod\u00e9r\u00e9e \u00e0 forte selon les crit\u00e8res conventionnels. \u00bb<\/em><\/p>\n<p>Notez qu&rsquo;en syntaxe typographique APA, les symboles statistiques s&rsquo;\u00e9crivent obligatoirement en italique (<em>F<\/em>, <em>p<\/em>, <em>R\u00b2<\/em>). De plus, lorsque la grandeur statistique ne peut math\u00e9matiquement jamais d\u00e9passer 1 (comme c&rsquo;est le cas des probabilit\u00e9s <em>p<\/em> ou du R\u00b2 ordinaire), le z\u00e9ro pr\u00e9c\u00e9dant la virgule (ou le point d\u00e9cimal anglo-saxon) est conventionnellement omis dans les manuscrits anglophones.<\/p>\n<p>Dans les tableaux synth\u00e9tiques de r\u00e9gressions hi\u00e9rarchiques par blocs, chaque \u00e9tape d&rsquo;incorporation de variables doit comporter les colonnes suivantes : les coefficients non standardis\u00e9s <em>B<\/em> avec leurs erreurs-types (<em>SE<\/em>), les coefficients standardis\u00e9s <em>\u03b2<\/em>, la statistique de Student <em>t<\/em>, le <em>R\u00b2<\/em> ordinaire, le <em>R\u00b2 ajust\u00e9<\/em>, ainsi que la statistique de variation incr\u00e9mentale <em>\u0394R\u00b2<\/em> et sa significativit\u00e9 propre (<em>F change<\/em>). Ce niveau de d\u00e9tail permet aux relecteurs d&rsquo;\u00e9valuer imm\u00e9diatement si l&rsquo;adjonction d&rsquo;un bloc de variables s&rsquo;est traduite par un gain de variance explicative authentique ou par une simple d\u00e9gradation p\u00e9nalis\u00e9e des degr\u00e9s de libert\u00e9.<\/p>\n<h3>12.2 Automatisation du rapport scientifique avec Python et Markdown<\/h3>\n<p>L&rsquo;int\u00e9gration de la programmation lettr\u00e9e et des pipelines d&rsquo;analyse reproductibles via des environnements tels que <strong>Quarto<\/strong> ou <strong>JupyterLab<\/strong> permet de supprimer d\u00e9finitivement le recopiage manuel des r\u00e9sultats num\u00e9riques \u2014 source notoire d&rsquo;erreurs typographiques dans les publications scientifiques.<\/p>\n<p>En tirant parti des m\u00e9thodes d&rsquo;exportation de Pandas et Statsmodels, un script d&rsquo;analyse peut formater directement les matrices de coefficients et les statistiques d&rsquo;ad\u00e9quation (incluant <code>rsquared_adj<\/code>) au format LaTeX (<code>df.to_latex(index=False, escape=False)<\/code>) ou au format HTML\/Markdown pour les documents dynamiques. L&rsquo;incorporation d&rsquo;instructions Python imbriqu\u00e9es directement dans le texte d&rsquo;un document Quarto (ex\u00e9cutant par exemple <code>`python f\"{model.rsquared_adj:.3f}\"`<\/code>) garantit que toute r\u00e9vision du jeu de donn\u00e9es brut ou tout filtrage d&rsquo;observations aberrantes se r\u00e9percute instantan\u00e9ment sur l&rsquo;ensemble des valeurs num\u00e9riques et des tableaux du manuscrit final.<\/p>\n<p>Cette automatisation rigoureuse consolide l&rsquo;auditabilit\u00e9 de la recherche scientifique, permettant \u00e0 des pairs examinateurs de r\u00e9ex\u00e9cuter l&rsquo;int\u00e9gralit\u00e9 du pipeline computationnel depuis l&rsquo;ingestion brute jusqu&rsquo;au rendu textuel du R-carr\u00e9 ajust\u00e9.<\/p>\n<h3>12.3 Synth\u00e8se critique et recommandations m\u00e9thodologiques finales<\/h3>\n<p>Au terme de cette exploration th\u00e9orique et computationnelle, il convient de synth\u00e9tiser les directives d\u00e9cisionnelles que tout mod\u00e9lisateur doit garder \u00e0 l&rsquo;esprit lors de l&rsquo;\u00e9valuation de mod\u00e8les lin\u00e9aires :<\/p>\n<ul>\n<li><strong>Le R-carr\u00e9 ordinaire ne doit jamais \u00eatre utilis\u00e9 isol\u00e9ment<\/strong> pour d\u00e9partager des sp\u00e9cifications multivari\u00e9es embo\u00eet\u00e9es ou non embo\u00eet\u00e9es, en raison de son incapacit\u00e9 structurelle \u00e0 p\u00e9naliser la prolif\u00e9ration param\u00e9trique.<\/li>\n<li><strong>Le R-carr\u00e9 ajust\u00e9 est un indicateur de parcimonie in-sample indispensable<\/strong> : son utilisation est particuli\u00e8rement critique lorsque le ratio d&rsquo;observations par pr\u00e9dicteur <em>n \/ k<\/em> est restreint (inf\u00e9rieur \u00e0 30:1), contexte dans lequel il immunise l&rsquo;analyste contre le surapprentissage fallacieux.<\/li>\n<li><strong>L&rsquo;adjonction d&rsquo;une variable n&rsquo;accro\u00eet le R\u00b2 ajust\u00e9 que si sa statistique |t| exc\u00e8de 1<\/strong> : cette condition math\u00e9matique, bien que plus souple que le seuil formel de rejet <em>p &lt; 0,05<\/em>, offre un rep\u00e8re analytique imm\u00e9diat pour proscrire les r\u00e9gresseurs d\u00e9pourvus de tout apport informatif.<\/li>\n<li><strong>L&rsquo;arsenal logiciel Python doit \u00eatre choisi selon l&rsquo;objectif poursuivi<\/strong> : NumPy pour l&rsquo;optimisation algorithmique pure et la compr\u00e9hension vectorielle matricielle ; Scikit-Learn pour l&rsquo;int\u00e9gration dans des pipelines d&rsquo;apprentissage supervis\u00e9 et de s\u00e9lection de caract\u00e9ristiques ; Statsmodels pour l&rsquo;inf\u00e9rence formelle, l&rsquo;acc\u00e8s natif imm\u00e9diat \u00e0 <code>rsquared_adj<\/code> et les diagnostics \u00e9conom\u00e9triques exhaustifs.<\/li>\n<li><strong>La m\u00e9trique statistique ne remplace jamais l&rsquo;assise th\u00e9orique<\/strong> : un R-carr\u00e9 ajust\u00e9 \u00e9lev\u00e9 n&rsquo;atteste en rien de la nature causale des associations observ\u00e9es, et un mod\u00e8le parcimonieux fond\u00e9 sur des principes th\u00e9oriques substantiels pr\u00e9vaudra toujours sur une \u00e9quation empiriquement optimis\u00e9e mais \u00e9pist\u00e9mologiquement st\u00e9rile.<\/li>\n<\/ul>\n<p>Pour les horizons m\u00e9thodologiques avanc\u00e9s, l&rsquo;analyste sera amen\u00e9 \u00e0 explorer les extensions contemporaines de cette m\u00e9trique : les coefficients de d\u00e9termination conditionnels et marginaux dans les mod\u00e8les lin\u00e9aires mixtes hi\u00e9rarchiques (<em>LMM<\/em>), ainsi que les m\u00e9thodes de r\u00e9gularisation p\u00e9nalis\u00e9e (<strong>Ridge<\/strong>, <strong>Lasso<\/strong>, <strong>ElasticNet<\/strong>), o\u00f9 la complexit\u00e9 n&rsquo;est plus contr\u00f4l\u00e9e par le simple d\u00e9compte discret des degr\u00e9s de libert\u00e9 <em>k<\/em>, mais par la contraction continue de la norme des coefficients r\u00e9gul\u00e9s.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Cohen, J. (1988). <em>Statistical Power Analysis for the Behavioral Sciences<\/em> (2nd ed.). Lawrence Erlbaum Associates.<\/li>\n<li>Ferguson, C. J. (2009). An effect size primer: A guide for clinicians and researchers. <em>Professional Psychology: Research and Practice<\/em>, 40(5), 532\u2013538. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/a0015808\">https:\/\/doi.org\/10.1037\/a0015808<\/a><\/li>\n<li>Hair, J. F., Black, W. C., Babin, B. J., &amp; Anderson, R. E. (2019). <em>Multivariate Data Analysis<\/em> (8th ed.). Cengage Learning.<\/li>\n<li>McKinney, W. (2022). <em>Python for Data Analysis: Data Wrangling with pandas, NumPy, and Jupyter<\/em> (3rd ed.). O&rsquo;Reilly Media.<\/li>\n<li>Pedhazur, E. J. (1997). <em>Multiple Regression in Behavioral Research: Explanation and Prediction<\/em> (3rd ed.). Wadsworth Publishing.<\/li>\n<li>Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., Blondel, M., Prettenhofer, P., Weiss, R., Dubourg, V., Vanderplas, J., Passos, A., Cournapeau, D., Brucher, M., Perrot, M., &amp; Duchesnay, E. (2011). Scikit-learn: Machine learning in Python. <em>Journal of Machine Learning Research<\/em>, 12, 2825\u20132830. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/jmlr.org\/papers\/v12\/pedregosa11a.html\">https:\/\/jmlr.org\/papers\/v12\/pedregosa11a.html<\/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> (pp. 92\u201396). <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.25080\/Majora-ebaa42b7-011\">https:\/\/doi.org\/10.25080\/Majora-ebaa42b7-011<\/a><\/li>\n<li>Theil, H. (1961). <em>Economic Forecasts and Policy<\/em> (Contributions to Economic Analysis, Vol. 7). North-Holland Publishing Company.<\/li>\n<li>van Rossum, G., &amp; Drake, F. L. (2009). <em>Python 3 Reference Manual<\/em>. CreateSpace.<\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour calculer et interpr\u00e9ter le R-carr\u00e9 ajust\u00e9 en Python avec scikit-learn, statsmodels et NumPy dans l&rsquo;analyse de r\u00e9gression.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2465],"tags":[2468,2902,1910,2904,2903],"class_list":["post-1299","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-programmation-python","tag-python","tag-r-carre-ajuste","tag-regression-lineaire","tag-scikit-learn","tag-statsmodels"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1299","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=1299"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1299\/revisions"}],"predecessor-version":[{"id":16217,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1299\/revisions\/16217"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1299"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1299"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1299"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}