{"id":6088,"date":"2026-09-16T04:41:34","date_gmt":"2026-09-16T04:41:34","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-somme-carres-residus-python\/"},"modified":"2026-09-16T04:41:34","modified_gmt":"2026-09-16T04:41:34","slug":"comment-calculer-somme-carres-residus-python","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-somme-carres-residus-python","title":{"rendered":"Comment calculer la somme des carr\u00e9s des r\u00e9sidus en Python"},"content":{"rendered":"<p>La mod\u00e9lisation statistique contemporaine repose sur une dialectique constante entre les structures th\u00e9oriques id\u00e9alis\u00e9es et la contingence irr\u00e9ductible des donn\u00e9es empiriques. Au c\u0153ur de cette confrontation se situe l&rsquo;analyse des r\u00e9sidus, composante fondamentale de l&rsquo;inf\u00e9rence statistique et de l&rsquo;apprentissage automatique supervis\u00e9. Dans le cadre de la r\u00e9gression lin\u00e9aire classique, l&rsquo;\u00e9valuation de la qualit\u00e9 de l&rsquo;ajustement d&rsquo;un mod\u00e8le math\u00e9matique aux ph\u00e9nom\u00e8nes observ\u00e9s n\u00e9cessite des indicateurs quantitatifs rigoureux. Parmi ces m\u00e9triques, la somme des carr\u00e9s des r\u00e9sidus, commun\u00e9ment d\u00e9sign\u00e9e sous l&rsquo;acronyme RSS pour <em>Residual Sum of Squares<\/em> ou SSR pour <em>Sum of Squared Residuals<\/em>, occupe une place pr\u00e9pond\u00e9rante. Elle constitue la fonction de co\u00fbt canonique dont la minimisation sous-tend le paradigme des moindres carr\u00e9s ordinaires, tout en servant d&rsquo;ancrage m\u00e9thodologique \u00e0 de multiples tests d&rsquo;hypoth\u00e8ses et crit\u00e8res de s\u00e9lection de mod\u00e8les.<\/p>\n<p>L&rsquo;av\u00e8nement de l&rsquo;\u00e9cosyst\u00e8me scientifique Python a profond\u00e9ment transform\u00e9 la mise en \u0153uvre pratique de ces concepts th\u00e9oriques. Gr\u00e2ce \u00e0 des environnements de calcul matriciel de haute performance tels que <a rel=\"noopener noreferrer\" href=\"https:\/\/numpy.org\/\">NumPy<\/a>, des biblioth\u00e8ques de manipulation de donn\u00e9es expressives \u00e0 l&rsquo;instar de <a rel=\"noopener noreferrer\" href=\"https:\/\/pandas.pydata.org\/\">Pandas<\/a>, ainsi que des cadres sp\u00e9cialis\u00e9s pour la mod\u00e9lisation statistique et pr\u00e9dictive comme <a rel=\"noopener noreferrer\" href=\"https:\/\/www.statsmodels.org\/\">Statsmodels<\/a> et <a rel=\"noopener noreferrer\" href=\"https:\/\/scikit-learn.org\/\">Scikit-Learn<\/a>, les chercheurs et analystes disposent d&rsquo;un arsenal logiciel sans pr\u00e9c\u00e9dent. N\u00e9anmoins, l&rsquo;accessibilit\u00e9 technique de ces biblioth\u00e8ques masque parfois des subtilit\u00e9s alg\u00e9briques, num\u00e9riques et distributionnelles cruciales, risquant de rel\u00e9guer l&rsquo;\u00e9valuation des mod\u00e8les \u00e0 une ex\u00e9cution automatis\u00e9e et acritique.<\/p>\n<p>Cet article propose une exploration exhaustive, rigoureuse et multidimensionnelle du calcul et de l&rsquo;interpr\u00e9tation de la somme des carr\u00e9s des r\u00e9sidus en Python. En ancrant notre r\u00e9flexion dans les fondements math\u00e9matiques de l&rsquo;alg\u00e8bre lin\u00e9aire et du calcul diff\u00e9rentiel, nous \u00e9tudierons comment cette grandeur s&rsquo;articule avec la d\u00e9composition globale de la variance. Nous analyserons ensuite ses applications pratiques, notamment dans les sciences comportementales et psychom\u00e9triques o\u00f9 la variabilit\u00e9 interindividuelle et les erreurs de mesure soul\u00e8vent des d\u00e9fis m\u00e9thodologiques consid\u00e9rables. Enfin, nous d\u00e9taillerons trois m\u00e9thodologies d&rsquo;impl\u00e9mentation distinctes au sein de l&rsquo;environnement Python, compl\u00e9t\u00e9es par des protocoles avanc\u00e9s de diagnostic des r\u00e9sidus et des recommandations pour pr\u00e9venir les \u00e9cueils de l&rsquo;inf\u00e9rence empirique.<\/p>\n<h2>1. Fondements th\u00e9oriques des r\u00e9sidus en r\u00e9gression lin\u00e9aire<\/h2>\n<h3>1.1 D\u00e9finition conceptuelle et math\u00e9matique de l&rsquo;erreur r\u00e9siduelle<\/h3>\n<p>Dans le cadre formel du mod\u00e8le de r\u00e9gression lin\u00e9aire multiple, la relation postul\u00e9e entre une variable d\u00e9pendante continue, not\u00e9e <em>y<\/em>, et un ensemble de <em>p<\/em> variables explicatives ou r\u00e9gresseurs, regroup\u00e9s au sein d&rsquo;un vecteur de caract\u00e9ristiques <em>x<\/em>, est soumise \u00e0 une composante d\u00e9terministe et \u00e0 une composante stochastique. L&rsquo;expression math\u00e9matique g\u00e9n\u00e9rale r\u00e9gissant la <em>i<\/em>-\u00e8me observation, issue d&rsquo;un \u00e9chantillon al\u00e9atoire de taille <em>n<\/em>, s&rsquo;\u00e9tablit selon l&rsquo;\u00e9quation :<\/p>\n<p><em>y<sub>i<\/sub> = &beta;<sub>0<\/sub> + &beta;<sub>1<\/sub>x<sub>i1<\/sub> + &beta;<sub>2<\/sub>x<sub>i2<\/sub> + &#8230; + &beta;<sub>p<\/sub>x<sub>ip<\/sub> + &epsilon;<sub>i<\/sub><\/em><\/p>\n<p>Dans cette formulation, les coefficients &beta; repr\u00e9sentent les param\u00e8tres r\u00e9els et inconnus de la population, tandis que le terme &epsilon;<sub>i<\/sub> d\u00e9signe la perturbation stochastique ou l&rsquo;erreur th\u00e9orique inh\u00e9rente au processus g\u00e9n\u00e9rateur de donn\u00e9es. Cette erreur th\u00e9orique &epsilon;<sub>i<\/sub> demeure par nature inobservable, puisqu&rsquo;elle appartient \u00e0 la structure latente de la population sous-jacente. \u00c0 l&rsquo;inverse, d\u00e8s lors que les param\u00e8tres &beta; sont estim\u00e9s empiriquement \u00e0 partir des donn\u00e9es de l&rsquo;\u00e9chantillon par des estimateurs statistiques, le mod\u00e8le produit une valeur pr\u00e9dite ou ajust\u00e9e pour chaque individu statistique, traditionnellement not\u00e9e <em>&ycirc;<sub>i<\/sub><\/em>. L&rsquo;\u00e9cart alg\u00e9brique individuel entre la valeur empirique r\u00e9ellement observ\u00e9e <em>y<sub>i<\/sub><\/em> et cette projection d\u00e9terministe <em>&ycirc;<sub>i<\/sub><\/em> constitue le r\u00e9sidu empirique, formellement d\u00e9fini par :<\/p>\n<p><em>e<sub>i<\/sub> = y<sub>i<\/sub> &#8211; &ycirc;<sub>i<\/sub><\/em><\/p>\n<p>Une distinction \u00e9pist\u00e9mologique stricte doit \u00eatre \u00e9tablie entre l&rsquo;erreur stochastique &epsilon;<sub>i<\/sub> et le r\u00e9sidu empirique <em>e<sub>i<\/sub><\/em>. Alors que le terme d&rsquo;erreur th\u00e9orique rel\u00e8ve de l&rsquo;ontologie du ph\u00e9nom\u00e8ne physique ou psychologique mod\u00e9lis\u00e9 \u2014 englobant l&rsquo;ensemble des d\u00e9terminants infinit\u00e9simaux non captur\u00e9s par les covariables \u2014, le r\u00e9sidu empirique est un artefact calculatoire conditionn\u00e9 par l&rsquo;\u00e9chantillon sp\u00e9cifique pr\u00e9lev\u00e9 et par l&rsquo;estimateur math\u00e9matique retenu. Le <a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/Gauss-MarkovTheorem.html\">th\u00e9or\u00e8me de Gauss-Markov<\/a> formalise les hypoth\u00e8ses structurales requises pour que les estimateurs obtenus par la m\u00e9thode des moindres carr\u00e9s soient sans biais et de variance minimale (propri\u00e9t\u00e9 BLUE : <em>Best Linear Unbiased Estimator<\/em>).<\/p>\n<p>Ces hypoth\u00e8ses imposent que l&rsquo;esp\u00e9rance conditionnelle des perturbations stochastiques soit strictement nulle, que leur variance soit homosc\u00e9dastique (constante pour toute configuration des r\u00e9gresseurs) et qu&rsquo;il n&rsquo;existe aucune corr\u00e9lation s\u00e9rielle ou spatiale entre les diff\u00e9rentes erreurs d&rsquo;observation. L&rsquo;implication imm\u00e9diate de l&rsquo;esp\u00e9rance math\u00e9matique nulle des r\u00e9sidus dans le mod\u00e8le empirique avec constante d&rsquo;interception est que la moyenne arithm\u00e9tique des r\u00e9sidus calcul\u00e9s est rigoureusement \u00e9gale \u00e0 z\u00e9ro, ce qui conditionne la coh\u00e9rence g\u00e9om\u00e9trique de l&rsquo;espace de projection des variables.<\/p>\n<h3>1.2 R\u00f4le des r\u00e9sidus dans l&rsquo;\u00e9valuation empirique des mod\u00e8les<\/h3>\n<p>Le r\u00e9sidu empirique constitue l&rsquo;unit\u00e9 fondamentale de mesure de la discordance entre le formalisme math\u00e9matique postul\u00e9 et la ph\u00e9nom\u00e9nologie observ\u00e9e. Dans toute investigation empirique, un mod\u00e8le statistique ne pr\u00e9tend pas reproduire l&rsquo;int\u00e9gralit\u00e9 de la r\u00e9alit\u00e9 physique ou sociale, mais aspire \u00e0 en capturer les tendances structurelles majeures. Les r\u00e9sidus repr\u00e9sentent pr\u00e9cis\u00e9ment la concr\u00e9tisation de l&rsquo;\u00e9chec partiel du mod\u00e8le \u00e0 appr\u00e9hender la totalit\u00e9 de l&rsquo;information contenue dans les donn\u00e9es. Par cons\u00e9quent, l&rsquo;analyse d\u00e9taill\u00e9e de leur magnitude globale, de leur distribution spatiale et de leurs corr\u00e9lations constitue la premi\u00e8re \u00e9tape d&rsquo;une \u00e9valuation critique de la sp\u00e9cification adopt\u00e9e.<\/p>\n<p>L&rsquo;agr\u00e9gation des r\u00e9sidus individuels permet de quantifier la part des variations de la variable r\u00e9ponse qui demeure inexpliqu\u00e9e par les covariables s\u00e9lectionn\u00e9es. Si les r\u00e9sidus d&rsquo;un mod\u00e8le pr\u00e9sentent une amplitude d\u00e9mesur\u00e9e par rapport \u00e0 la dispersion de la variable d\u00e9pendante d&rsquo;origine, le pouvoir explicatif de l&rsquo;architecture choisie doit \u00eatre r\u00e9\u00e9valu\u00e9. De surcro\u00eet, la pr\u00e9sence de structures persistantes, telles que des courbures syst\u00e9matiques ou des variations d&rsquo;amplitude conditionnelles dans les graphiques de r\u00e9sidus, d\u00e9nonce imm\u00e9diatement une mauvaise sp\u00e9cification fonctionnelle, comme l&rsquo;omission d&rsquo;effets d&rsquo;interaction ou de termes non lin\u00e9aires d&rsquo;ordre sup\u00e9rieur.<\/p>\n<p>C&rsquo;est sur cette mesure de l&rsquo;\u00e9cart que repose le principe fondateur de l&rsquo;optimisation par la m\u00e9thode des <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Ordinary_least_squares\">moindres carr\u00e9s ordinaires (MCO)<\/a>. Plut\u00f4t que de rechercher une solution subjective ou heuristique pour r\u00e9concilier les divergences empiriques, la m\u00e9thode des moindres carr\u00e9s \u00e9rige la minimisation math\u00e9matique d&rsquo;une fonction agr\u00e9g\u00e9e des r\u00e9sidus en crit\u00e8re unique d&rsquo;admissibilit\u00e9 des param\u00e8tres. Les r\u00e9sidus cessent ainsi d&rsquo;\u00eatre de simples sous-produits du calcul pour devenir la mati\u00e8re premi\u00e8re du processus d&rsquo;estimation param\u00e9trique.<\/p>\n<h3>1.3 Application des r\u00e9sidus aux sciences psychologiques et comportementales<\/h3>\n<p>Dans les disciplines psychologiques, psychom\u00e9triques et neurocognitives, la mod\u00e9lisation statistique se heurte \u00e0 une complexit\u00e9 idiosyncrasique li\u00e9e \u00e0 la variabilit\u00e9 intrins\u00e8que des comportements humains. Contrairement aux sciences physiques o\u00f9 les conditions exp\u00e9rimentales peuvent \u00eatre isol\u00e9es avec une pr\u00e9cision quasi absolue, l&rsquo;\u00e9valuation des dimensions psychologiques implique une part substantielle de variance d&rsquo;erreur. Les tests psychom\u00e9triques mesurant des construits latents tels que l&rsquo;intelligence fluide, l&rsquo;anxi\u00e9t\u00e9 situationnelle ou l&rsquo;efficacit\u00e9 de la m\u00e9moire de travail sont sujets \u00e0 des bruits de mesure consid\u00e9rables, attribuables \u00e0 la fatigue passag\u00e8re, aux fluctuations attentionnelles ou \u00e0 des biais de d\u00e9sirabilit\u00e9 sociale.<\/p>\n<p>Dans ce contexte, la variabilit\u00e9 r\u00e9siduelle acquiert une signification th\u00e9orique profonde. Loin de repr\u00e9senter uniquement du bruit parasitaire sans int\u00e9r\u00eat, le r\u00e9sidu quantifie la variabilit\u00e9 comportementale interindividuelle qui r\u00e9siste aux taxonomies psychologiques conventionnelles. Lorsqu&rsquo;un mod\u00e8le de r\u00e9gression pr\u00e9dit la vitesse de traitement cognitif en fonction de l&rsquo;\u00e2ge chronologique et du niveau d&rsquo;\u00e9tudes, un r\u00e9sidu positif substantiel chez un sujet d\u00e9signe un fonctionnement neurocognitif sup\u00e9rieur aux attentes th\u00e9oriques de son groupe de r\u00e9f\u00e9rence, signalant potentiellement des facteurs de r\u00e9serve cognitive ou des strat\u00e9gies compensatoires non r\u00e9pertori\u00e9es.<\/p>\n<p>De plus, dans le diagnostic comportemental clinique et l&rsquo;\u00e9valuation neurocognitive informatis\u00e9e, l&rsquo;\u00e9valuation de la pr\u00e9cision des scores pr\u00e9dits s&rsquo;av\u00e8re vitale. L&rsquo;estimation d&rsquo;un score attendu \u00e0 une t\u00e2che de rappel diff\u00e9r\u00e9 permet d&rsquo;identifier des dissociations pathologiques lorsque le r\u00e9sidu observ\u00e9 s&rsquo;\u00e9carte significativement de la distribution normale des erreurs d&rsquo;ajustement. Ainsi, la quantification syst\u00e9matique des r\u00e9sidus fournit aux praticiens et chercheurs un barom\u00e8tre objectif pour isoler l&rsquo;influence des covariables environnementales et \u00e9valuer la validit\u00e9 incr\u00e9mentale de nouveaux protocoles th\u00e9rapeutiques ou psychop\u00e9dagogiques.<\/p>\n<h2>2. D\u00e9finition et propri\u00e9t\u00e9s de la somme des carr\u00e9s des r\u00e9sidus (RSS)<\/h2>\n<h3>2.1 Formulation alg\u00e9brique et calcul diff\u00e9rentiel<\/h3>\n<p>Pour d\u00e9passer l&rsquo;\u00e9chelle individuelle des observations et synth\u00e9tiser l&rsquo;impr\u00e9cision d&rsquo;ajustement globale d&rsquo;un mod\u00e8le lin\u00e9aire sur l&rsquo;ensemble d&rsquo;un \u00e9chantillon statistique de taille <em>n<\/em>, il est n\u00e9cessaire de recourir \u00e0 une fonction de co\u00fbt scalaire. L&rsquo;addition directe des r\u00e9sidus simples s&rsquo;av\u00e8re inop\u00e9rante, car en pr\u00e9sence d&rsquo;une constante dans le mod\u00e8le, les \u00e9carts positifs et n\u00e9gatifs se compensent exactement, produisant une somme rigoureusement nulle :<\/p>\n<p>&sum;<sub>i=1<\/sub><sup>n<\/sup> <em>e<sub>i<\/sub><\/em> = 0<\/p>\n<p>Pour conf\u00e9rer un poids positif \u00e0 chaque d\u00e9viation par rapport \u00e0 la surface de r\u00e9gression, l&rsquo;approche math\u00e9matique standard consiste \u00e0 \u00e9lever chaque r\u00e9sidu au carr\u00e9 avant de proc\u00e9der \u00e0 la sommation. La somme des carr\u00e9s des r\u00e9sidus, not\u00e9e RSS (<em>Residual Sum of Squares<\/em>), s&rsquo;\u00e9nonce formellement sous la forme :<\/p>\n<p><em>RSS<\/em> = &sum;<sub>i=1<\/sub><sup>n<\/sup> <em>e<sub>i<\/sub><\/em><sup>2<\/sup> = &sum;<sub>i=1<\/sub><sup>n<\/sup> (<em>y<sub>i<\/sub> &#8211; &ycirc;<sub>i<\/sub><\/em>)<sup>2<\/sup> = &sum;<sub>i=1<\/sub><sup>n<\/sup> (<em>y<sub>i<\/sub> &#8211; &beta;<sub>0<\/sub> &#8211; &sum;<sub>j=1<\/sub><sup>p<\/sup> &beta;<sub>j<\/sub>x<sub>ij<\/sub><\/em>)<sup>2<\/sup><\/p>\n<p>Le choix m\u00e9thodologique d&rsquo;une fonction quadratique plut\u00f4t que d&rsquo;une somme de valeurs absolues (telle qu&#8217;employ\u00e9e dans la r\u00e9gression de la m\u00e9diane ou r\u00e9gression L1) repose sur des imp\u00e9ratifs d&rsquo;analyse diff\u00e9rentielle et de g\u00e9om\u00e9trie riemannienne. La fonction de perte quadratique est une fonction de classe C<sup>&infin;<\/sup>, c&rsquo;est-\u00e0-dire infiniment diff\u00e9rentiable par rapport \u00e0 chaque coefficient &beta;<sub>j<\/sub>. Cette d\u00e9rivabilit\u00e9 continue garantit l&rsquo;existence d&rsquo;un gradient analytique unique et continu sur l&rsquo;ensemble de l&rsquo;espace des param\u00e8tres r\u00e9els &reals;<sup>p+1<\/sup>. En d\u00e9rivant partiellement la fonction <em>RSS<\/em> par rapport \u00e0 chaque param\u00e8tre &beta;<sub>k<\/sub> et en \u00e9galisant ces d\u00e9riv\u00e9es partielles \u00e0 z\u00e9ro, on engendre le syst\u00e8me d&rsquo;\u00e9quations normales qui admet une solution ferm\u00e9e exacte, sans n\u00e9cessiter d&rsquo;approximations par sous-gradients caract\u00e9ristiques des fonctions non d\u00e9rivables en z\u00e9ro.<\/p>\n<h3>2.2 D\u00e9composition de la variance : lien entre TSS, ESS et RSS<\/h3>\n<p>L&rsquo;importance fondamentale de la RSS d\u00e9coule de son insertion canonique dans l&rsquo;\u00e9quation de d\u00e9composition de la variabilit\u00e9 totale de la variable d\u00e9pendante. La somme totale des carr\u00e9s, d\u00e9sign\u00e9e par <em>TSS<\/em> pour <em>Total Sum of Squares<\/em>, quantifie la dispersion globale des observations empiriques <em>y<sub>i<\/sub><\/em> autour de leur moyenne arithm\u00e9tique <em>&ymacr;<\/em> :<\/p>\n<p><em>TSS<\/em> = &sum;<sub>i=1<\/sub><sup>n<\/sup> (<em>y<sub>i<\/sub> &#8211; &ymacr;<\/em>)<sup>2<\/sup><\/p>\n<p>Sous r\u00e9serve que le mod\u00e8le de r\u00e9gression int\u00e8gre un terme d&rsquo;interception (constante), l&rsquo;identit\u00e9 alg\u00e9brique fondamentale \u00e9tablit que cette dispersion totale se partitionne de fa\u00e7on univoque et exhaustive en deux composantes orthogonales :<\/p>\n<p><em>TSS = ESS + RSS<\/em><\/p>\n<p>o\u00f9 <em>ESS<\/em> d\u00e9signe la somme des carr\u00e9s expliqu\u00e9e (<em>Explained Sum of Squares<\/em>, parfois appel\u00e9e <em>Model Sum of Squares<\/em>), formul\u00e9e par :<\/p>\n<p><em>ESS<\/em> = &sum;<sub>i=1<\/sub><sup>n<\/sup> (<em>&ycirc;<sub>i<\/sub> &#8211; &ymacr;<\/em>)<sup>2<\/sup><\/p>\n<p>Cette d\u00e9composition d\u00e9coule directement du <a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/PythagoreanTheorem.html\">th\u00e9or\u00e8me de Pythagore<\/a> appliqu\u00e9 \u00e0 l&rsquo;espace euclidien \u00e0 <em>n<\/em> dimensions &reals;<sup>n<\/sup>. Le vecteur des observations centr\u00e9es <em>(y &#8211; &ymacr;1)<\/em> se projette orthogonalement sur le sous-espace vectoriel engendr\u00e9 par les colonnes de la matrice de design. Le vecteur des valeurs ajust\u00e9es centr\u00e9es <em>(&ycirc; &#8211; &ymacr;1)<\/em> et le vecteur des r\u00e9sidus empiriques <em>e<\/em> forment ainsi deux vecteurs rigoureusement orthogonaux, garantissant que le produit scalaire <em>(&ycirc; &#8211; &ymacr;1)<sup>T<\/sup>e<\/em> soit nul. Par cons\u00e9quent, le ratio <em>RSS \/ TSS<\/em> exprime la proportion exacte de variance inexpliqu\u00e9e par la r\u00e9gression, constituant le compl\u00e9ment direct \u00e0 l&rsquo;unit\u00e9 de la proportion de variance expliqu\u00e9e.<\/p>\n<h3>2.3 Interpr\u00e9tation de la magnitude de la RSS<\/h3>\n<p>L&rsquo;\u00e9valuation brute de la magnitude de la RSS soul\u00e8ve d&rsquo;importantes questions d&rsquo;\u00e9chelle et d&rsquo;interpr\u00e9tabilit\u00e9 m\u00e9trologique. Sur le plan heuristique, une RSS tendant asymptotiquement vers z\u00e9ro indique que les points de donn\u00e9es co\u00efncident parfaitement avec l&rsquo;hyperplan de r\u00e9gression calcul\u00e9. Dans cette configuration limite, l&rsquo;impr\u00e9cision d\u00e9terministe est enti\u00e8rement r\u00e9sorb\u00e9e et les valeurs observ\u00e9es correspondent exactement aux valeurs pr\u00e9dites par les variables s\u00e9lectionn\u00e9es.<\/p>\n<p>Cependant, la valeur absolue de la RSS est une grandeur dimensionn\u00e9e, intrins\u00e8quement assujettie \u00e0 l&rsquo;unit\u00e9 de mesure dans laquelle la variable d\u00e9pendante <em>y<\/em> est exprim\u00e9e. Si une variable de temps de r\u00e9action cognitive est convertie de secondes en millisecondes, chaque r\u00e9sidu individuel se trouve multipli\u00e9 par un facteur 1000, ce qui entra\u00eene ipso facto la multiplication de la RSS par un facteur 1 000 000 (10<sup>6<\/sup>), sans que la structure de d\u00e9pendance lin\u00e9aire sous-jacente ou la qualit\u00e9 pr\u00e9dictive intrins\u00e8que du mod\u00e8le n&rsquo;aient subi la moindre alt\u00e9ration.<\/p>\n<p>Cette d\u00e9pendance d&rsquo;\u00e9chelle interdit toute comparaison directe de la RSS brute entre des mod\u00e8les portant sur des variables crit\u00e8res distinctes ou mesur\u00e9es sur des \u00e9chelles h\u00e9t\u00e9rog\u00e8nes. Pour comparer la performance d&rsquo;ajustement entre \u00e9tudes h\u00e9t\u00e9roclites, les statisticiens doivent imp\u00e9rativement recourir \u00e0 des transformations adimensionnelles, telles que la standardisation pr\u00e9alable des variables (scores Z), ou d\u00e9river de la RSS des indices normalis\u00e9s comme le coefficient de d\u00e9termination <em>R<sup>2<\/sup><\/em> ou l&rsquo;erreur quadratique moyenne standardis\u00e9e.<\/p>\n<h2>3. Importance de la RSS dans l&rsquo;analyse quantitative et psychom\u00e9trique<\/h2>\n<h3>3.1 Validation de l&rsquo;ad\u00e9quation d&rsquo;un mod\u00e8le en psychologie exp\u00e9rimentale<\/h3>\n<p>Dans le domaine de la psychologie exp\u00e9rimentale, o\u00f9 les protocoles testent des hypoth\u00e8ses relatives \u00e0 l&rsquo;architecture des fonctions mentales, l&rsquo;\u00e9valuation de la validit\u00e9 d&rsquo;ajustement d&rsquo;un mod\u00e8le via la RSS constitue une \u00e9tape d\u00e9terminante de la d\u00e9marche inf\u00e9rentielle. Lors d&rsquo;exp\u00e9riences portant sur l&rsquo;apprentissage perceptif ou la consolidation mn\u00e9sique sous diff\u00e9rentes conditions d&rsquo;interf\u00e9rence, plusieurs mod\u00e8les fonctionnels concurrents s&rsquo;opposent fr\u00e9quemment dans la litt\u00e9rature : mod\u00e8les lin\u00e9aires additifs, fonctions logarithmiques simulant la loi de Weber-Fechner, ou cin\u00e9tiques exponentielles d&rsquo;oubli formalis\u00e9es par Ebbinghaus.<\/p>\n<p>La confrontation de la RSS issue de ces diff\u00e9rentes sp\u00e9cifications architecturales permet de trancher empiriquement entre des th\u00e9ories rivales du traitement de l&rsquo;information. Un mod\u00e8le qui parvient \u00e0 r\u00e9duire drastiquement la RSS sans recourir \u00e0 une prolif\u00e9ration injustifi\u00e9e de degr\u00e9s de libert\u00e9 d\u00e9montre une capacit\u00e9 sup\u00e9rieure \u00e0 refl\u00e9ter les contraintes neurobiologiques de la t\u00e2che exp\u00e9rimentale.<\/p>\n<p>Parall\u00e8lement, la quantification de la RSS s&rsquo;av\u00e8re indispensable pour estimer la part de variance attribuable aux traits de personnalit\u00e9 latents lors de l&rsquo;\u00e9tude de r\u00e9ponses comportementales en milieu social. En comparant la RSS d&rsquo;un mod\u00e8le n&rsquo;incluant que des facteurs contextuels imm\u00e9diats \u00e0 celle d&rsquo;un mod\u00e8le enrichi par les dimensions du mod\u00e8le \u00e0 cinq facteurs (Big Five), le chercheur quantifie rigoureusement l&rsquo;impact structurel des dispositions individuelles stables sur l&rsquo;ajustement global.<\/p>\n<h3>3.2 Comparaison de mod\u00e8les hi\u00e9rarchiques et embo\u00eet\u00e9s<\/h3>\n<p>L&rsquo;une des fonctions fondamentales de la RSS dans l&rsquo;analyse quantitative moderne r\u00e9side dans son utilisation comme pivot math\u00e9matique pour la comparaison formelle de mod\u00e8les statistiques embo\u00eet\u00e9s. Dans de nombreux programmes de recherche, la mod\u00e9lisation proc\u00e8de par \u00e9tapes hi\u00e9rarchiques : un mod\u00e8le restreint <em>M<sub>1<\/sub><\/em> contenant <em>p<sub>1<\/sub><\/em> pr\u00e9dicteurs de base est compar\u00e9 \u00e0 un mod\u00e8le \u00e9tendu ou complet <em>M<sub>2<\/sub><\/em> incluant <em>p<sub>2<\/sub><\/em> r\u00e9gresseurs additionnels (avec <em>p<sub>2<\/sub> &gt; p<sub>1<\/sub><\/em>).<\/p>\n<p>Sur le plan g\u00e9om\u00e9trique, le sous-espace vectoriel de pr\u00e9diction engendr\u00e9 par le mod\u00e8le restreint est strictement inclus dans celui du mod\u00e8le complet. En vertu des lois de la projection orthogonale dans les espaces de Hilbert, l&rsquo;adjonction de variables explicatives suppl\u00e9mentaires entra\u00eene syst\u00e9matiquement une diminution m\u00e9canique (ou au minimum une stagnation) de la somme des carr\u00e9s des r\u00e9sidus, de telle sorte que <em>RSS<sub>2<\/sub> &le; RSS<sub>1<\/sub><\/em>. L&rsquo;enjeu m\u00e9thodologique consiste \u00e0 d\u00e9terminer si cette r\u00e9duction constat\u00e9e, not\u00e9e <em>&Delta;RSS = RSS<sub>1<\/sub> &#8211; RSS<sub>2<\/sub><\/em>, refl\u00e8te un apport pr\u00e9dictif substantiel et statistiquement fiable, ou si elle d\u00e9coule d&rsquo;une simple exploitation probabiliste du bruit d&rsquo;\u00e9chantillonnage.<\/p>\n<p>Ce questionnement est r\u00e9solu par le calcul de la statistique du test <em>F<\/em> de changement (ou test de Fisher-Snedecor pour mod\u00e8les embo\u00eet\u00e9s), dont la formule repose explicitement sur la variation relative des RSS :<\/p>\n<p><em>F = [ (RSS<sub>1<\/sub> &#8211; RSS<sub>2<\/sub>) \/ (df<sub>1<\/sub> &#8211; df<sub>2<\/sub>) ] \/ [ RSS<sub>2<\/sub> \/ df<sub>2<\/sub> ]<\/em><\/p>\n<p>o\u00f9 <em>df<sub>1<\/sub><\/em> et <em>df<sub>2<\/sub><\/em> d\u00e9signent respectivement les degr\u00e9s de libert\u00e9 r\u00e9siduels des mod\u00e8les restreint et complet (avec <em>df = n &#8211; p &#8211; 1<\/em>). Ce ratio probabiliste met en \u00e9quilibre le principe de parcimonie scientifique, h\u00e9rit\u00e9 du rasoir d&rsquo;Ockham, et le gain d&rsquo;ajustement empirique mesur\u00e9 par la d\u00e9croissance de la RSS.<\/p>\n<h3>3.3 Lien fondamental avec l&rsquo;estimation du maximum de vraisemblance<\/h3>\n<p>La port\u00e9e statistique de la RSS d\u00e9passe le seul cadre g\u00e9om\u00e9trique des moindres carr\u00e9s ordinaires gr\u00e2ce \u00e0 son articulation th\u00e9orique avec le paradigme du <a rel=\"noopener noreferrer\" href=\"https:\/\/plato.stanford.edu\/entries\/statistics\/\">maximum de vraisemblance (Maximum Likelihood Estimation &#8211; MLE)<\/a>. Sous l&rsquo;hypoth\u00e8se classique selon laquelle les perturbations stochastiques &epsilon;<sub>i<\/sub> sont ind\u00e9pendantes et distribu\u00e9es selon une loi normale centr\u00e9e d&rsquo;\u00e9cart-type constant &sigma; (loi <em>N(0, &sigma;<sup>2<\/sup>)<\/em>), la fonction de log-vraisemblance de l&rsquo;\u00e9chantillon conditionnelle aux param\u00e8tres du vecteur &beta; et de la variance r\u00e9siduelle s&rsquo;\u00e9crit :<\/p>\n<p><em>ln L(&beta;, &sigma;<sup>2<\/sup> | y, X) = &#8211; (n \/ 2) ln(2&pi;) &#8211; (n \/ 2) ln(&sigma;<sup>2<\/sup>) &#8211; [ 1 \/ (2&sigma;<sup>2<\/sup>) ] &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &#8211; x<sub>i<\/sub><sup>T<\/sup>&beta;)<sup>2<\/sup><\/em><\/p>\n<p>En observant cette \u00e9quation, on constate imm\u00e9diatement que la somme apparaissant dans le terme de droite est pr\u00e9cis\u00e9ment l&rsquo;expression analytique de la RSS. D\u00e8s lors, pour toute variance r\u00e9siduelle &sigma;<sup>2<\/sup> fix\u00e9e positive, maximiser la log-vraisemblance par rapport aux coefficients de r\u00e9gression &beta; \u00e9quivaut strictement \u00e0 minimiser le terme soustrait, c&rsquo;est-\u00e0-dire \u00e0 minimiser la RSS. Les estimateurs des MCO et les estimateurs du maximum de vraisemblance des coefficients &beta; sont donc rigoureusement identiques sous condition de normalit\u00e9 des perturbations.<\/p>\n<p>De plus, cette connexion th\u00e9orique permet d&rsquo;estimer directement la variance du bruit al\u00e9atoire de l&rsquo;\u00e9chantillon \u00e0 partir de la RSS. L&rsquo;estimateur sans biais de la variance r\u00e9siduelle de la population, not\u00e9 <em>s<sup>2<\/sup><\/em> ou <em>&sigma;&circ;<sup>2<\/sup><\/em>, s&rsquo;obtient en divisant la RSS par les degr\u00e9s de libert\u00e9 r\u00e9siduels de l&rsquo;\u00e9chantillon :<\/p>\n<p><em>s<sup>2<\/sup> = RSS \/ (n &#8211; p &#8211; 1)<\/em><\/p>\n<p>La racine carr\u00e9e de cette quantit\u00e9, <em>s<\/em>, d\u00e9finit l&rsquo;erreur standard d&rsquo;estimation (<em>Standard Error of the Estimate<\/em>). Dans les contextes psychom\u00e9triques d&rsquo;\u00e9valuation standardis\u00e9e du quotient intellectuel ou des aptitudes scolaires, cette grandeur quantifie la marge d&rsquo;incertitude type encadrant le diagnostic individuel et r\u00e9git le calcul des intervalles de confiance pr\u00e9dictifs.<\/p>\n<h2>4. Configuration de l&rsquo;environnement scientifique Python<\/h2>\n<h3>4.1 Choix et installation des biblioth\u00e8ques analytiques<\/h3>\n<p>La r\u00e9alisation d&rsquo;analyses statistiques avanc\u00e9es et le calcul computationnel de la RSS requi\u00e8rent un \u00e9cosyst\u00e8me logiciel alliant performance algorithmique de bas niveau et souplesse syntaxique de haut niveau. Le langage de programmation Python s&rsquo;est impos\u00e9 comme la r\u00e9f\u00e9rence internationale pour ces op\u00e9rations gr\u00e2ce \u00e0 un ensemble de biblioth\u00e8ques hautement sp\u00e9cialis\u00e9es et compl\u00e9mentaires.<\/p>\n<p>Au c\u0153ur de cet outillage figure <a rel=\"noopener noreferrer\" href=\"https:\/\/numpy.org\/\">NumPy<\/a>, qui constitue le socle du calcul num\u00e9rique vectoriel. D\u00e9velopp\u00e9e en langage C et interfac\u00e9e avec des biblioth\u00e8ques de sous-programmes d&rsquo;alg\u00e8bre lin\u00e9aire hautement optimis\u00e9es telles que BLAS (<em>Basic Linear Algebra Subprograms<\/em>) et LAPACK (<em>Linear Algebra Package<\/em>), la biblioth\u00e8que NumPy introduit la structure d&rsquo;objet <code>ndarray<\/code>. Cette architecture en m\u00e9moire contigu\u00eb permet d&rsquo;ex\u00e9cuter des op\u00e9rations matricielles complexes \u2014 inversions, factorisations de Cholesky, d\u00e9compositions QR et produits tensoriels \u2014 \u00e0 des vitesses proches du silicium mat\u00e9riel, sans le surco\u00fbt d&rsquo;interpr\u00e9tation dynamique propre \u00e0 la boucle native Python.<\/p>\n<p>Pour la manipulation de structures tabulaires h\u00e9t\u00e9rog\u00e8nes repr\u00e9sentatives de collectes exp\u00e9rimentales r\u00e9elles, <a rel=\"noopener noreferrer\" href=\"https:\/\/pandas.pydata.org\/\">Pandas<\/a> offre l&rsquo;objet <code>DataFrame<\/code>. Cette structure facilite l&rsquo;alignement index\u00e9 des donn\u00e9es, la gestion robuste des m\u00e9tadonn\u00e9es, le filtrage bool\u00e9en et le traitement transparent des valeurs manquantes. Pour la phase de mod\u00e9lisation math\u00e9matique, deux philosophies distinctes mais convergentes sont mobilis\u00e9es :<\/p>\n<ul>\n<li><strong>Statsmodels :<\/strong> Biblioth\u00e8que orient\u00e9e vers l&rsquo;\u00e9conom\u00e9trie formelle et la mod\u00e9lisation statistique pure. Elle produit des r\u00e9sum\u00e9s exhaustifs int\u00e9grant les valeurs de p, les intervalles de confiance param\u00e9triques, les tests omnibus et les m\u00e9triques r\u00e9siduelles canoniques comme la RSS.<\/li>\n<li><strong>Scikit-Learn :<\/strong> Biblioth\u00e8que orient\u00e9e vers l&rsquo;apprentissage automatique pr\u00e9dictif, l&rsquo;ing\u00e9nierie des caract\u00e9ristiques et l&rsquo;optimisation par validation crois\u00e9e \u00e0 grande \u00e9chelle.<\/li>\n<\/ul>\n<p>L&rsquo;installation recommand\u00e9e pour garantir l&rsquo;absence de conflits de d\u00e9pendances binaires s&rsquo;op\u00e8re dans un environnement virtuel d\u00e9di\u00e9 (par exemple via l&rsquo;outil <code>venv<\/code> ou <code>conda<\/code>). L&rsquo;ex\u00e9cution en ligne de commande standard se formule ainsi :<\/p>\n<p><code>pip install numpy pandas statsmodels scikit-learn scipy matplotlib seaborn<\/code><\/p>\n<h3>4.2 Structure de projet et reproductibilit\u00e9 computationnelle<\/h3>\n<p>La reproductibilit\u00e9 computationnelle est un imp\u00e9ratif \u00e9thique et scientifique fondamental de la recherche universitaire contemporaine. L&rsquo;obtention d&rsquo;une valeur num\u00e9rique sp\u00e9cifique pour la RSS au sein d&rsquo;une simulation statistique ou d&rsquo;un protocole empirique ne poss\u00e8de aucune valeur probante si l&rsquo;environnement d&rsquo;ex\u00e9cution et les conditions de g\u00e9n\u00e9ration des nombres pseudo-al\u00e9atoires ne sont pas formellement verrouill\u00e9s.<\/p>\n<p>Une organisation de projet rigoureuse impose la s\u00e9gr\u00e9gation des dossiers de travail : un sous-r\u00e9pertoire <code>data\/<\/code> contenant les donn\u00e9es brutes immutables prot\u00e9g\u00e9es en \u00e9criture, un dossier <code>src\/<\/code> ou <code>scripts\/<\/code> rassemblant les modules fonctionnels test\u00e9s unitairement, et un espace <code>notebooks\/<\/code> r\u00e9serv\u00e9 \u00e0 l&rsquo;exploration interactive document\u00e9e. Les versions exactes des paquets logiciels doivent \u00eatre fig\u00e9es dans un fichier d\u00e9claratif <code>requirements.txt<\/code> ou <code>environment.yml<\/code> g\u00e9n\u00e9r\u00e9 automatiquement.<\/p>\n<p>De surcro\u00eet, la variabilit\u00e9 stochastique induite par les proc\u00e9dures d&rsquo;\u00e9chantillonnage al\u00e9atoire, les fractionnements de validation crois\u00e9e ou la g\u00e9n\u00e9ration de donn\u00e9es synth\u00e9tiques requiert la fixation explicite des graines al\u00e9atoires (<em>random seeds<\/em>). En Python, cela n\u00e9cessite d&rsquo;initialiser les g\u00e9n\u00e9rateurs d&rsquo;\u00e9tats \u00e0 la fois dans la biblioth\u00e8que standard et au niveau de l&rsquo;infrastructure NumPy :<\/p>\n<p><code>import numpy as np<\/code><br \/>\n<code>import random<\/code><br \/>\n<code>SEED = 42<\/code><br \/>\n<code>np.random.seed(SEED)<\/code><br \/>\n<code>random.seed(SEED)<\/code><\/p>\n<p>Dans les versions modernes de NumPy (ant\u00e9rieures et post\u00e9rieures \u00e0 la version 1.17), la bonne pratique acad\u00e9mique pr\u00e9conise l&rsquo;instanciation directe d&rsquo;un g\u00e9n\u00e9rateur ind\u00e9pendant via l&rsquo;instruction <code>rng = np.random.default_rng(SEED)<\/code>, qui isole l&rsquo;\u00e9tat stochastique du script de toute interf\u00e9rence globale caus\u00e9e par des modules tiers.<\/p>\n<h2>5. Pr\u00e9paration et structuration des donn\u00e9es exp\u00e9rimentales avec Pandas<\/h2>\n<h3>5.1 Mod\u00e9lisation du jeu de donn\u00e9es d&rsquo;apprentissage et de performance<\/h3>\n<p>Pour illustrer de mani\u00e8re concr\u00e8te les calculs math\u00e9matiques et la mise en \u0153uvre informatique de la RSS, construisons un jeu de donn\u00e9es repr\u00e9sentatif d&rsquo;une cohorte acad\u00e9mique dans un d\u00e9partement de sciences cognitives. L&rsquo;objectif consiste \u00e0 mod\u00e9liser la performance finale obtenue \u00e0 une \u00e9preuve d&rsquo;\u00e9valuation cognitive formelle (variable d\u00e9pendante <em>y<\/em>, not\u00e9e <code>score_examen<\/code>, exprim\u00e9e sur une \u00e9chelle continue standardis\u00e9e allant de 0 \u00e0 100) en fonction de deux pr\u00e9dicteurs psychop\u00e9dagogiques continus :<\/p>\n<ul>\n<li><code>heures_etude<\/code> : Le volume cumul\u00e9 d&rsquo;heures d&rsquo;apprentissage individuel investi par l&rsquo;\u00e9tudiant au cours des quatre semaines pr\u00e9c\u00e9dant l&rsquo;\u00e9valuation.<\/li>\n<li><code>examens_blancs<\/code> : Le nombre d&rsquo;\u00e9preuves d&rsquo;entra\u00eenement formatives pass\u00e9es sous conditions chronom\u00e9tr\u00e9es avant l&rsquo;\u00e9valuation finale.<\/li>\n<\/ul>\n<p>Sous Pandas, la cr\u00e9ation et la configuration initiale de ce DataFrame s&rsquo;op\u00e8rent par le biais de dictionnaires de vecteurs typ\u00e9s ou par la lecture de fichiers tabulaires d\u00e9limit\u00e9s. L&rsquo;inspection pr\u00e9alable des types de variables via l&rsquo;attribut <code>.dtypes<\/code> et la production des statistiques descriptives de position et de dispersion au moyen de la m\u00e9thode <code>.describe()<\/code> constituent des \u00e9tapes obligatoires pour s&rsquo;assurer que les donn\u00e9es respectent les hypoth\u00e8ses de continuit\u00e9 num\u00e9rique (variables de type <code>float64<\/code>) et que leurs distributions empiriques s&rsquo;inscrivent dans des plages physiologiquement et psychom\u00e9triquement r\u00e9alistes.<\/p>\n<h3>5.2 Nettoyage statistique et traitement des valeurs manquantes<\/h3>\n<p>Les collectes de donn\u00e9es empiriques en milieu universitaire ou clinique souffrent fr\u00e9quemment de d\u00e9faillances m\u00e9trologiques : omissions de saisie de la part des participants, d\u00e9connexions intempestives des syst\u00e8mes informatis\u00e9s de recueil, ou refus explicite de r\u00e9ponse \u00e0 certains items sp\u00e9cifiques. Face \u00e0 des valeurs manquantes (repr\u00e9sent\u00e9es par l&rsquo;objet sentinelle <code>numpy.nan<\/code>), l&rsquo;analyste doit choisir une strat\u00e9gie d&rsquo;assainissement statistique adapt\u00e9e \u00e0 la taille de son \u00e9chantillon.<\/p>\n<p>La suppression pure et simple des lignes incompl\u00e8tes, connue sous le terme d&rsquo;analyse des cas complets (<em>listwise deletion<\/em> via <code>pandas.DataFrame.dropna()<\/code>), n&rsquo;est m\u00e9thodologiquement d\u00e9fendable que si le m\u00e9canisme g\u00e9n\u00e9rateur des donn\u00e9es manquantes satisfait \u00e0 l&rsquo;hypoth\u00e8se stricte de manque compl\u00e8tement al\u00e9atoire (MCAR : <em>Missing Completely at Random<\/em>). Lorsque cette hypoth\u00e8se est viol\u00e9e ou lorsque l&rsquo;\u00e9chantillon initial est modeste \u2014 ce qui est fr\u00e9quent lors d&rsquo;exp\u00e9rimentations psychologiques co\u00fbteuses \u2014, la suppression r\u00e9duit drastiquement la puissance statistique du test et risque d&rsquo;introduire des biais d&rsquo;\u00e9chantillonnage s\u00e9v\u00e8res.<\/p>\n<p>Des strat\u00e9gies d&rsquo;imputation univari\u00e9e (m\u00e9diane, moyenne conditionnelle) ou multivari\u00e9e par cha\u00eenes d&rsquo;\u00e9quations it\u00e9ratives (MICE via le module <code>sklearn.impute.IterativeImputer<\/code>) permettent de pr\u00e9server la compl\u00e9tude du tableau tout en conservant la structure de covariance inter-items. De m\u00eame, l&rsquo;identification des aberrations num\u00e9riques par l&rsquo;examen des \u00e9carts interquartiles (crit\u00e8re de Tukey) ou par le calcul des distances de Mahalanobis avant l&rsquo;ajustement du mod\u00e8le garantit que le calcul ult\u00e9rieur de la RSS ne sera pas contamin\u00e9 par des erreurs de saisie instrumentale.<\/p>\n<h2>6. M\u00e9thode 1 : Calcul de la RSS \u00e0 partir des principes fondamentaux avec NumPy<\/h2>\n<h3>6.1 Formulation matricielle et calcul des moindres carr\u00e9s manuels<\/h3>\n<p>L&rsquo;impl\u00e9mentation de la r\u00e9gression lin\u00e9aire sous sa forme matricielle originelle permet d&rsquo;appr\u00e9hender le fonctionnement intime du calcul statistique sans l&rsquo;interm\u00e9diaire de couches d&rsquo;abstraction logicielle. Soit un \u00e9chantillon de <em>n<\/em> observations pourvu de <em>p<\/em> covariables explicatives. Nous d\u00e9finissons la matrice de design <em>X<\/em>, de dimensions <em>n &times; (p + 1)<\/em>, dont la premi\u00e8re colonne est uniform\u00e9ment constitu\u00e9e de scalaires unitaires (1) afin de mod\u00e9liser le terme d&rsquo;interception &beta;<sub>0<\/sub> :<\/p>\n<p>Le vecteur d&rsquo;observations de la variable d\u00e9pendante <em>y<\/em> poss\u00e8de la dimension <em>n &times; 1<\/em>. Le syst\u00e8me d&rsquo;\u00e9quations normales qui annule le gradient de la RSS s&rsquo;\u00e9crit sous forme matricielle condens\u00e9e :<\/p>\n<p><em>(X<sup>T<\/sup>X)&beta;&circ; = X<sup>T<\/sup>y<\/em><\/p>\n<p>Dans l&rsquo;hypoth\u00e8se standard o\u00f9 les r\u00e9gresseurs sont lin\u00e9airement ind\u00e9pendants, la matrice de Gram <em>X<sup>T<\/sup>X<\/em> est inversible (de rang plein \u00e9gal \u00e0 <em>p + 1<\/em>). L&rsquo;estimateur unique des moindres carr\u00e9s ordinaires s&rsquo;obtient alors par l&rsquo;\u00e9quation matricielle :<\/p>\n<p><em>&beta;&circ; = (X<sup>T<\/sup>X)<sup>-1<\/sup>X<sup>T<\/sup>y<\/em><\/p>\n<p>Avec NumPy, apr\u00e8s avoir assembl\u00e9 la matrice de design par concat\u00e9nation horizontale d&rsquo;une colonne de uns via <code>np.column_stack<\/code> ou <code>np.c_<\/code>, le vecteur des coefficients estim\u00e9s &beta;&circ; peut \u00eatre calcul\u00e9 analytiquement en mobilisant la fonction d&rsquo;inversion matricielle <code>np.linalg.inv<\/code> combin\u00e9e \u00e0 l&rsquo;op\u00e9rateur de multiplication matricielle <code>@<\/code> :<\/p>\n<p><code>X_design = np.c_[np.ones(X.shape[0]), X]<\/code><br \/>\n<code>beta_hat = np.linalg.inv(X_design.T @ X_design) @ X_design.T @ y<\/code><\/p>\n<p>Toutefois, en calcul scientifique rigoureux, l&rsquo;inversion explicite d&rsquo;une matrice via <code>inv()<\/code> est d\u00e9conseill\u00e9e car elle amplifie les erreurs de troncature num\u00e9rique d\u00e8s lors que la matrice approche la singularit\u00e9. On privil\u00e9gie la fonction <code>np.linalg.solve<\/code> qui r\u00e9sout le syst\u00e8me lin\u00e9aire par \u00e9limination gaussienne, ou la fonction d\u00e9di\u00e9e aux moindres carr\u00e9s <code>np.linalg.lstsq<\/code> qui s&rsquo;appuie sur la d\u00e9composition en valeurs singuli\u00e8res (SVD), m\u00e9thode insensible \u00e0 la multicolin\u00e9arit\u00e9 parfaite.<\/p>\n<p>Une fois le vecteur de param\u00e8tres &beta;&circ; obtenu, le calcul du vecteur des valeurs pr\u00e9dites <em>&ycirc;<\/em> s&rsquo;effectue simplement par le produit scalaire :<\/p>\n<p><em>&ycirc; = X &beta;&circ;<\/em><\/p>\n<p>ce qui correspond en syntaxe NumPy \u00e0 l&rsquo;instruction vectoris\u00e9e <code>y_pred = X_design @ beta_hat<\/code>.<\/p>\n<h3>6.2 Vectorisation du calcul des r\u00e9sidus et sommation quadratique<\/h3>\n<p>Poss\u00e9dant \u00e0 la fois le vecteur empirique <em>y<\/em> et le vecteur des projections calcul\u00e9es <em>&ycirc;<\/em>, la d\u00e9termination des r\u00e9sidus individuels s&rsquo;op\u00e8re par soustraction vectorielle directe \u00e9l\u00e9ment par \u00e9l\u00e9ment. Gr\u00e2ce au m\u00e9canisme de diffusion (<em>broadcasting<\/em>) de NumPy, cette op\u00e9ration s&rsquo;ex\u00e9cute \u00e0 tr\u00e8s haute vitesse :<\/p>\n<p><code>residus = y - y_pred<\/code><\/p>\n<p>L&rsquo;obtention de la RSS requiert ensuite l&rsquo;\u00e9l\u00e9vation de chaque composante du vecteur <code>residus<\/code> au carr\u00e9, suivie de leur sommation univari\u00e9e compl\u00e8te. Deux formulations syntaxiques \u00e9quivalentes en pr\u00e9cision mais diff\u00e9rentes dans leur chemin d&rsquo;ex\u00e9cution en bas niveau sont envisageables :<\/p>\n<p><code>rss_numpy = np.sum(np.square(residus))<\/code><br \/>\n<code>rss_dot = np.dot(residus, residus)<\/code><\/p>\n<p>L&rsquo;approche utilisant le produit scalaire vectoriel <code>np.dot(residus, residus)<\/code> ou <code>residus.T @ residus<\/code> s&rsquo;av\u00e8re algorithmiquement optimale. Elle \u00e9vite l&rsquo;allocation d&rsquo;un tableau temporaire interm\u00e9diaire en m\u00e9moire vive pour stocker les carr\u00e9s des r\u00e9sidus, calculant l&rsquo;accumulation de la somme directement dans un registre du processeur gr\u00e2ce aux instructions SIMD (<em>Single Instruction, Multiple Data<\/em>) de la biblioth\u00e8que BLAS sous-jacente.<\/p>\n<p>L&rsquo;analyse de la complexit\u00e9 algorithmique r\u00e9v\u00e8le que le calcul de la RSS par cette m\u00e9thode fondamentale est domin\u00e9 par l&rsquo;\u00e9tape de multiplication matricielle <em>X<sup>T<\/sup>X<\/em>, dont la complexit\u00e9 temporelle est d&rsquo;ordre <em>O(n &times; p<sup>2<\/sup>)<\/em>, tandis que la r\u00e9solution du syst\u00e8me lin\u00e9aire requiert un temps <em>O(p<sup>3<\/sup>)<\/em>. La phase finale de soustraction et de sommation quadratique ne prend qu&rsquo;un temps lin\u00e9aire <em>O(n)<\/em>, ce qui rend l&rsquo;\u00e9valuation r\u00e9siduelle extr\u00eamement v\u00e9loce, m\u00eame sur des matrices de donn\u00e9es contenant plusieurs millions d&rsquo;enregistrements.<\/p>\n<h3>6.3 Avantages p\u00e9dagogiques et validation algorithmique<\/h3>\n<p>L&rsquo;impl\u00e9mentation de la RSS \u00e0 partir des fondements math\u00e9matiques avec NumPy pr\u00e9sente un int\u00e9r\u00eat didactique et \u00e9pist\u00e9mologique capital. Elle dissipe l&rsquo;effet de \u00ab bo\u00eete noire \u00bb inh\u00e9rent aux frameworks statistiques de haut niveau qui masquent la g\u00e9om\u00e9trie sous-jacente des donn\u00e9es sous des couches de classes et d&rsquo;objets opaques. En manipulant directement les vecteurs et les matrices, le chercheur ou le data scientist ancre sa compr\u00e9hension intuitive de la r\u00e9gression dans la projection orthogonale au sein d&rsquo;un espace de Hilbert.<\/p>\n<p>Sur le plan de la m\u00e9trologie logicielle, cette m\u00e9thode artisanale constitue l&rsquo;\u00e9talon de r\u00e9f\u00e9rence pour la validation crois\u00e9e d&rsquo;algorithmes. Lorsque de nouveaux pipelines d&rsquo;analyse pr\u00e9dictive sont d\u00e9ploy\u00e9s sur des infrastructures infonuagiques, confronter le r\u00e9sultat d&rsquo;un module complexe \u00e0 l&rsquo;\u00e9valuation purement vectorielle NumPy permet d&rsquo;attester de la conformit\u00e9 du code avec la th\u00e9orie statistique classique.<\/p>\n<p>Enfin, cette approche met en lumi\u00e8re les r\u00e9alit\u00e9s du calcul num\u00e9rique en virgule flottante double pr\u00e9cision standard (norme IEEE 754). Lors du calcul de la diff\u00e9rence <em>(y<sub>i<\/sub> &#8211; &ycirc;<sub>i<\/sub>)<\/em> pour des observations extr\u00eamement volumineuses associ\u00e9es \u00e0 des r\u00e9sidus infinit\u00e9simaux, des ph\u00e9nom\u00e8nes de perte de signification ou d&rsquo;annulation catastrophique peuvent se produire. L&rsquo;exp\u00e9rimentateur apprend ainsi \u00e0 surveiller le conditionnement matriciel via <code>np.linalg.cond(X)<\/code>, renfor\u00e7ant la fiabilit\u00e9 m\u00e9trologique de ses conclusions scientifiques.<\/p>\n<h2>7. M\u00e9thode 2 : Calcul de la RSS avec la biblioth\u00e8que Statsmodels<\/h2>\n<h3>7.1 Sp\u00e9cification et ajustement du mod\u00e8le OLS<\/h3>\n<p>Tandis que NumPy excelle dans le calcul matriciel brut, la biblioth\u00e8que <a rel=\"noopener noreferrer\" href=\"https:\/\/www.statsmodels.org\/\">Statsmodels<\/a> a \u00e9t\u00e9 sp\u00e9cifiquement con\u00e7ue pour reproduire la rigueur formelle des logiciels d&rsquo;\u00e9conom\u00e9trie et de biom\u00e9trie traditionnels (tels que R, Stata ou SAS). Son utilisation privil\u00e9gie l&rsquo;inf\u00e9rence rigoureuse, la validit\u00e9 des tests statistiques et la tra\u00e7abilit\u00e9 des hypoth\u00e8ses plut\u00f4t que la seule performance brute de calcul.<\/p>\n<p>Pour ajuster un mod\u00e8le lin\u00e9aire de moindres carr\u00e9s ordinaires via Statsmodels, une attention particuli\u00e8re doit \u00eatre port\u00e9e \u00e0 la d\u00e9finition de la matrice des r\u00e9gresseurs. Contrairement \u00e0 Scikit-Learn ou aux formules symboliques \u00e0 la R, la classe <code>statsmodels.api.OLS<\/code> n&rsquo;ajoute pas de terme d&rsquo;interception par d\u00e9faut. Omettre d&rsquo;inclure la constante contraint l&rsquo;hyperplan de r\u00e9gression \u00e0 passer obligatoirement par l&rsquo;origine du rep\u00e8re cart\u00e9sien <em>(0, 0, &#8230;, 0)<\/em>. Une telle sp\u00e9cification erron\u00e9e invalide la d\u00e9composition de la variance <em>TSS = ESS + RSS<\/em> et peut conduire \u00e0 des coefficients de d\u00e9termination <em>R<sup>2<\/sup><\/em> artificiellement n\u00e9gatifs ou trompeurs.<\/p>\n<p>L&rsquo;adjonction formelle de cette colonne unitaire s&rsquo;effectue au moyen de la fonction utilitaire <code>sm.add_constant()<\/code> :<\/p>\n<p><code>import statsmodels.api as sm<\/code><br \/>\n<code>X_avec_constante = sm.add_constant(X_df)<\/code><br \/>\n<code>modele_stat = sm.OLS(y, X_avec_constante)<\/code><br \/>\n<code>resultats = modele_stat.fit()<\/code><\/p>\n<p>L&rsquo;appel \u00e0 la m\u00e9thode <code>.fit()<\/code> d\u00e9clenche l&rsquo;ex\u00e9cution des algorithmes d&rsquo;estimation matricielle et g\u00e9n\u00e8re un objet complexe de classe <code>RegressionResultsWrapper<\/code>, qui concentre l&rsquo;ensemble des estimations scalaires, tensorielles et asymptotiques du mod\u00e8le.<\/p>\n<h3>7.2 Extraction directe des attributs de r\u00e9sidus<\/h3>\n<p>Statsmodels propose plusieurs m\u00e9canismes transparents pour acc\u00e9der aux d\u00e9viations r\u00e9siduelles et quantifier instantan\u00e9ment la RSS sans devoir reprogrammer manuellement les sommations quadratiques.<\/p>\n<p>Le premier moyen consiste \u00e0 r\u00e9cup\u00e9rer le vecteur complet des r\u00e9sidus empiriques bruts <em>e<sub>i<\/sub><\/em> via l&rsquo;attribut <code>.resid<\/code> de l&rsquo;objet de r\u00e9sultats. Cet attribut retourne une structure de type <code>pandas.Series<\/code> ou <code>numpy.ndarray<\/code> selon la nature de la variable cible originellement transmise. Il est d\u00e8s lors ais\u00e9 de proc\u00e9der \u00e0 la sommation manuelle de ses carr\u00e9s :<\/p>\n<p><code>residus_sm = resultats.resid<\/code><br \/>\n<code>rss_calculee = np.sum(residus_sm ** 2)<\/code><\/p>\n<p>Toutefois, Statsmodels calcule d\u00e9j\u00e0 en interne cette grandeur lors de l&rsquo;ajustement du mod\u00e8le. L&rsquo;utilisateur peut donc interroger directement l&rsquo;attribut natif d\u00e9di\u00e9 :<\/p>\n<p><code>rss_directe = resultats.ssr<\/code><\/p>\n<p>Une mise au point terminologique s&rsquo;impose ici concernant la convention de nommage adopt\u00e9e par Statsmodels. Dans la documentation anglo-saxonne, la somme des carr\u00e9s des r\u00e9sidus est indiff\u00e9remment d\u00e9sign\u00e9e sous les sigles <em>RSS<\/em> (<em>Residual Sum of Squares<\/em>) ou <em>SSR<\/em> (<em>Sum of Squared Residuals<\/em>). Cependant, certains manuels d&rsquo;\u00e9conom\u00e9trie utilisent le sigle <em>SSR<\/em> pour d\u00e9signer la somme des carr\u00e9s de la r\u00e9gression (<em>Regression Sum of Squares<\/em>, que nous avons qualifi\u00e9e d&rsquo;<em>ESS<\/em>). Dans Statsmodels, l&rsquo;attribut <code>ssr<\/code> fait express\u00e9ment r\u00e9f\u00e9rence \u00e0 la somme des carr\u00e9s des r\u00e9sidus inexpliqu\u00e9s (<em>Sum of Squared Residuals<\/em>), tandis que la variance expliqu\u00e9e par le mod\u00e8le est stock\u00e9e sous l&rsquo;attribut <code>ess<\/code> (<em>Explained Sum of Squares<\/em>), et la variance totale sous l&rsquo;attribut <code>centered_tss<\/code>. Cette harmonisation interne \u00e9limine toute \u00e9quivoque d&rsquo;interpr\u00e9tation pour le statisticien.<\/p>\n<h3>7.3 Interpr\u00e9tation des m\u00e9triques compl\u00e9mentaires g\u00e9n\u00e9r\u00e9es<\/h3>\n<p>L&rsquo;un des atouts distinctifs de Statsmodels r\u00e9side dans la profondeur analytique des m\u00e9triques d&rsquo;accompagnement corr\u00e9l\u00e9es \u00e0 la RSS. L&rsquo;appel \u00e0 la m\u00e9thode <code>resultats.summary()<\/code> d\u00e9ploie un rapport textuel standardis\u00e9 pr\u00e9sentant une \u00e9valuation compl\u00e8te de l&rsquo;architecture stochastique.<\/p>\n<p>Au sein de ce rapport, la RSS intervient directement dans le calcul de l&rsquo;erreur quadratique moyenne non biais\u00e9e des r\u00e9sidus (<em>Mean Squared Error of the Model\/Residuals<\/em>), enregistr\u00e9e sous l&rsquo;attribut <code>resultats.mse_resid<\/code>. Cette grandeur est r\u00e9gie par la relation formelle :<\/p>\n<p><em>MSE<sub>resid<\/sub> = SSR \/ df<sub>resid<\/sub><\/em><\/p>\n<p>o\u00f9 <code>resultats.df_resid<\/code> repr\u00e9sente les degr\u00e9s de libert\u00e9 r\u00e9siduels r\u00e9els de l&rsquo;\u00e9chantillon, calcul\u00e9s comme <em>n &#8211; p &#8211; 1<\/em>. Cette division normalisatrice est fondamentale : elle corrige le biais syst\u00e9matique de l&rsquo;estimation de la variance qui appara\u00eetrait si l&rsquo;on se contentait de diviser la RSS par la taille totale de l&rsquo;\u00e9chantillon <em>n<\/em>. L&rsquo;acc\u00e8s \u00e0 cette grandeur non biais\u00e9e permet d&rsquo;alimenter directement les calculs d&rsquo;intervalles de confiance des param\u00e8tres &beta;, attestant de la compl\u00e9tude m\u00e9thodologique de la biblioth\u00e8que pour la recherche acad\u00e9mique.<\/p>\n<h2>8. M\u00e9thode 3 : Calcul de la RSS avec Scikit-Learn<\/h2>\n<h3>8.1 Entra\u00eenement du r\u00e9gresseur LinearRegression<\/h3>\n<p>Dans l&rsquo;\u00e9cosyst\u00e8me du g\u00e9nie logiciel et de la science des donn\u00e9es industrielles, la biblioth\u00e8que <a rel=\"noopener noreferrer\" href=\"https:\/\/scikit-learn.org\/\">Scikit-Learn<\/a> constitue l&rsquo;environnement de r\u00e9f\u00e9rence incontournable. Con\u00e7ue autour d&rsquo;une interface de programmation objet hautement unifi\u00e9e et coh\u00e9rente, elle structure l&rsquo;ensemble de ses algorithmes d&rsquo;apprentissage supervis\u00e9 autour du paradigme canonique <code>Estimator.fit(X, y)<\/code> et <code>Estimator.predict(X)<\/code>.<\/p>\n<p>Pour ajuster un mod\u00e8le lin\u00e9aire ordinaire avec Scikit-Learn, on instancie la classe <code>LinearRegression<\/code> issue du sous-module <code>sklearn.linear_model<\/code>. Contrairement \u00e0 Statsmodels, l&rsquo;interception &beta;<sub>0<\/sub> est g\u00e9r\u00e9e de mani\u00e8re totalement native et automatique gr\u00e2ce au param\u00e8tre bool\u00e9en <code>fit_intercept=True<\/code>, configur\u00e9 par d\u00e9faut \u00e0 l&rsquo;initialisation :<\/p>\n<p><code>from sklearn.linear_model import LinearRegression<\/code><br \/>\n<code>modele_sklearn = LinearRegression(fit_intercept=True)<\/code><br \/>\n<code>modele_sklearn.fit(X_train, y_train)<\/code><\/p>\n<p>L&rsquo;estimation des coefficients matriciels s&rsquo;ex\u00e9cute en d\u00e9l\u00e9guant les op\u00e9rations d&rsquo;optimisation aux sous-routines de factorisation LAPACK (notamment la m\u00e9thode <code>gelsd<\/code> bas\u00e9e sur la d\u00e9composition en valeurs singuli\u00e8res ou la d\u00e9composition QR). Une fois la phase d&rsquo;apprentissage finalis\u00e9e, la g\u00e9n\u00e9ration du vecteur des scores pr\u00e9dits s&rsquo;effectue simplement en appelant la m\u00e9thode pr\u00e9dictive :<\/p>\n<p><code>y_pred_sklearn = modele_sklearn.predict(X_train)<\/code><\/p>\n<h3>8.2 D\u00e9rivation de la RSS \u00e0 partir de la m\u00e9trique MSE<\/h3>\n<p>Une particularit\u00e9 fonctionnelle essentielle de Scikit-Learn doit \u00eatre soulign\u00e9e : la biblioth\u00e8que ne stocke pas directement la RSS en tant qu&rsquo;attribut membre de la classe <code>LinearRegression<\/code>. Guid\u00e9e par une philosophie orient\u00e9e vers la performance pr\u00e9dictive et la g\u00e9n\u00e9ralisation sur des donn\u00e9es tierces (\u00e9chantillons de test ou de validation), Scikit-Learn privil\u00e9gie les m\u00e9triques standardis\u00e9es d&rsquo;erreur moyenne par observation, au premier rang desquelles figure l&rsquo;erreur quadratique moyenne (<em>Mean Squared Error &#8211; MSE<\/em>).<\/p>\n<p>La MSE se d\u00e9finit formellement comme l&rsquo;esp\u00e9rance arithm\u00e9tique de la somme des carr\u00e9s des r\u00e9sidus, calcul\u00e9e sur la taille <em>N<\/em> du jeu de donn\u00e9es \u00e9valu\u00e9 :<\/p>\n<p><em>MSE = (1 \/ N) &sum;<sub>i=1<\/sub><sup>N<\/sup> (y<sub>i<\/sub> &#8211; &ycirc;<sub>i<\/sub>)<sup>2<\/sup> = (1 \/ N) &times; RSS<\/em><\/p>\n<p>Par cons\u00e9quent, pour reconstituer la RSS exacte au sens math\u00e9matique classique \u00e0 partir des outils d&rsquo;\u00e9valuation de Scikit-Learn, il convient de faire appel \u00e0 la fonction <code>mean_squared_error<\/code> du sous-module <code>sklearn.metrics<\/code>, puis de multiplier scalairement le r\u00e9sultat retourn\u00e9 par le nombre effectif d&rsquo;observations <em>N<\/em> composant le jeu de donn\u00e9es :<\/p>\n<p><code>from sklearn.metrics import mean_squared_error<\/code><br \/>\n<code>mse = mean_squared_error(y_train, y_pred_sklearn)<\/code><br \/>\n<code>rss_sklearn = mse * len(y_train)<\/code><\/p>\n<p>Une nuance critique doit \u00eatre relev\u00e9e par le statisticien m\u00e9ticuleux : la fonction <code>mean_squared_error<\/code> de Scikit-Learn applique une division brute par <em>N<\/em>, et non par les degr\u00e9s de libert\u00e9 r\u00e9siduels <em>n &#8211; p &#8211; 1<\/em>. Il s&rsquo;agit de la d\u00e9finition de la MSE au sens de l&rsquo;\u00e9valuation pr\u00e9dictive et non de l&rsquo;estimation non biais\u00e9e de la variance stochastique (laquelle divise par <em>n &#8211; p &#8211; 1<\/em>). La multiplication par <code>len(y_train)<\/code> restitue donc rigoureusement la valeur exacte de la RSS th\u00e9orique, sans la moindre divergence arithm\u00e9tique avec les r\u00e9sultats obtenus via NumPy ou Statsmodels.<\/p>\n<h3>8.3 Automatisation au sein de pipelines d&rsquo;\u00e9valuation crois\u00e9e<\/h3>\n<p>L&rsquo;int\u00e9r\u00eat pr\u00e9\u00e9minent de manipuler la RSS au moyen de Scikit-Learn r\u00e9side dans la capacit\u00e9 d&rsquo;int\u00e9grer ce crit\u00e8re d&rsquo;erreur au sein de structures algorithmiques modulaires sophistiqu\u00e9es, telles que les cha\u00eenes de traitement <code>sklearn.pipeline.Pipeline<\/code> et les sch\u00e9mas de validation crois\u00e9e <code>sklearn.model_selection.cross_validate<\/code>.<\/p>\n<p>Lors d&rsquo;\u00e9valuations pr\u00e9dictives rigoureuses visant \u00e0 tester la r\u00e9sistance d&rsquo;un mod\u00e8le face au surapprentissage, il est primordial de mesurer la stabilit\u00e9 de la RSS sur des plis d&rsquo;entra\u00eenement et de test successifs (protocole de validation crois\u00e9e \u00e0 <em>k<\/em> blocs, ou <em>k-fold cross-validation<\/em>). Pour cela, Scikit-Learn permet de d\u00e9finir une m\u00e9trique de score personnalis\u00e9e \u00e0 l&rsquo;aide de la fabrique logicielle <code>make_scorer<\/code> :<\/p>\n<p><code>from sklearn.metrics import make_scorer<\/code><br \/>\n<code>from sklearn.model_selection import cross_val_score<\/code><br \/>\n<code>def calcul_rss_custom(y_vrai, y_pred):<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;return np.sum((y_vrai - y_pred) ** 2)<\/code><br \/>\n<code>score_rss = make_scorer(calcul_rss_custom, greater_is_better=False)<\/code><br \/>\n<code>scores_cv = cross_val_score(modele_sklearn, X, y, cv=5, scoring=score_rss)<\/code><\/p>\n<p>Cette approche unifi\u00e9e offre aux chercheurs la facult\u00e9 de surveiller l&rsquo;\u00e9volution de la RSS non seulement sur les donn\u00e9es ayant servi \u00e0 calibrer les coefficients, mais \u00e9galement sur des sous-ensembles empiriques vierges. La divergence grandissante entre la RSS d&rsquo;entra\u00eenement et la RSS de validation crois\u00e9e signale imm\u00e9diatement une variance algorithmique excessive et la perte de validit\u00e9 externe du mod\u00e8le.<\/p>\n<h2>9. Diagnostic approfondi et analyse distributionnelle des r\u00e9sidus<\/h2>\n<h3>9.1 V\u00e9rification de l&rsquo;hypoth\u00e8se d&rsquo;homosc\u00e9dasticit\u00e9<\/h3>\n<p>L&rsquo;obtention d&rsquo;une valeur num\u00e9rique minimale pour la RSS ne garantit en aucune mani\u00e8re, \u00e0 elle seule, la validit\u00e9 statistique des inf\u00e9rences tir\u00e9es d&rsquo;un mod\u00e8le lin\u00e9aire. La fiabilit\u00e9 des estimations repose sur le respect scrupuleux des postulats distributionnels \u00e9dict\u00e9s par le mod\u00e8le de Gauss-Markov. Parmi ces r\u00e9quisits, l&rsquo;homosc\u00e9dasticit\u00e9 stipule que la variance conditionnelle des r\u00e9sidus demeure rigoureusement constante sur l&rsquo;ensemble du domaine de variation des r\u00e9gresseurs :<\/p>\n<p><em>Var(&epsilon;<sub>i<\/sub> | X) = &sigma;<sup>2<\/sup>, &forall; i &isin; {1, &#8230;, n}<\/em><\/p>\n<p>Lorsque cette condition n&rsquo;est pas remplie, le mod\u00e8le souffre d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 : la dispersion des erreurs s&rsquo;accro\u00eet ou se contracte en fonction des valeurs des variables pr\u00e9dictives. Dans cette situation, bien que l&rsquo;estimateur des moindres carr\u00e9s demeure formellement sans biais, il perd sa propri\u00e9t\u00e9 d&rsquo;efficacit\u00e9 optimale (la variance de l&rsquo;estimateur n&rsquo;est plus minimale), et les erreurs types associ\u00e9es aux coefficients sont substantiellement biais\u00e9es, invalidant les tests de significativit\u00e9 <em>t<\/em> de Student et les seuils de significativit\u00e9 <em>p<\/em>.<\/p>\n<p>L&rsquo;exploration visuelle de l&rsquo;homosc\u00e9dasticit\u00e9 s&rsquo;effectue classiquement par la construction d&rsquo;un diagramme de dispersion confrontant les valeurs pr\u00e9dites ajust\u00e9es <em>&ycirc;<sub>i<\/sub><\/em> sur l&rsquo;axe des abscisses aux r\u00e9sidus standardis\u00e9s sur l&rsquo;axe des ordonn\u00e9es. Une distribution saine pr\u00e9sente un nuage de points al\u00e9atoire homog\u00e8ne, contenu dans une bande horizontale sans structure d\u00e9celable. La d\u00e9tection d&rsquo;une morphologie en entonnoir, en \u00e9ventail ou en fuseau trahit la pr\u00e9sence aigu\u00eb d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9.<\/p>\n<p>Pour d\u00e9passer l&rsquo;appr\u00e9ciation graphique subjective, des proc\u00e9dures inf\u00e9rentielles rigoureuses doivent \u00eatre ex\u00e9cut\u00e9es en Python. Le test formel de Breusch-Pagan ainsi que le test de White, disponibles dans le module <code>statsmodels.stats.diagnostic<\/code> via les fonctions <code>het_breuschpagan<\/code> et <code>het_white<\/code>, calculent un multiplicateur de Lagrange confrontant l&rsquo;hypoth\u00e8se nulle d&rsquo;homosc\u00e9dasticit\u00e9 stricte \u00e0 l&rsquo;hypoth\u00e8se alternative d&rsquo;une variance r\u00e9siduelle d\u00e9pendante des covariables et de leurs produits crois\u00e9s. Si la valeur de p issue de ces tests descend sous le seuil critique conventionnel de 0,05, l&rsquo;homosc\u00e9dasticit\u00e9 est formellement rejet\u00e9e, contraignant l&rsquo;analyste \u00e0 adopter des matrices de variance-covariance robustes \u00e0 l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 (telles que les estimateurs sandwichs de White ou de Newey-West) ou \u00e0 appliquer des transformations stabilisatrices de variance (Box-Cox, log-transformation).<\/p>\n<h3>9.2 \u00c9valuation de la normalit\u00e9 et r\u00e9sidus studentis\u00e9s<\/h3>\n<p>Bien que la d\u00e9rivation alg\u00e9brique de l&rsquo;estimateur MCO ne requi\u00e8re aucune hypoth\u00e8se distributionnelle quant \u00e0 la forme de l&rsquo;erreur pour converger asymptotiquement vers la vraie valeur des param\u00e8tres, l&rsquo;exactitude des inf\u00e9rences en petits \u00e9chantillons \u2014 sp\u00e9cifiquement la construction des intervalles de confiance et la l\u00e9gitimit\u00e9 des tests d&rsquo;hypoth\u00e8ses univari\u00e9s <em>t<\/em> et multivari\u00e9s <em>F<\/em> \u2014 repose formellement sur l&rsquo;hypoth\u00e8se de normalit\u00e9 gaussienne des perturbations :<\/p>\n<p><em>&epsilon;<sub>i<\/sub> ~ i.i.d. N(0, &sigma;<sup>2<\/sup>)<\/em><\/p>\n<p>L&rsquo;\u00e9valuation empirique de cette normalit\u00e9 s&rsquo;appuie conjointement sur des visualisations diagnostiques et des batteries de tests probabilistes. L&rsquo;outil graphique roi est le trac\u00e9 quantile-quantile (<em>Q-Q Plot<\/em>), impl\u00e9ment\u00e9 dans la biblioth\u00e8que <a rel=\"noopener noreferrer\" href=\"https:\/\/scipy.org\/\">SciPy<\/a> et accessible sous Statsmodels via <code>sm.qqplot(residus, line='45', fit=True)<\/code>. Ce graphique confronte les quantiles empiriques de la distribution des r\u00e9sidus observ\u00e9s aux quantiles th\u00e9oriques d&rsquo;une loi normale univari\u00e9e standard. Tout d\u00e9crochement syst\u00e9matique le long de la ligne diagonale de r\u00e9f\u00e9rence \u2014 notamment aux extr\u00e9mit\u00e9s de la distribution \u2014 r\u00e9v\u00e8le la pr\u00e9sence de queues lourdes (leptokurticit\u00e9), de dissym\u00e9trie (asym\u00e9trie positive ou n\u00e9gative) ou de troncature d&rsquo;\u00e9chantillonnage.<\/p>\n<p>Sur le plan inf\u00e9rentiel, ces observations sont compl\u00e9t\u00e9es par les tests univari\u00e9s de normalit\u00e9. Le test de Shapiro-Wilk (fonction <code>scipy.stats.shapiro<\/code>) offre la puissance statistique la plus \u00e9lev\u00e9e pour les petits et moyens \u00e9chantillons, tandis que le test de Jarque-Bera (<code>scipy.stats.jarque_bera<\/code>), reposant sur l&rsquo;\u00e9valuation combin\u00e9e du coefficient d&rsquo;asym\u00e9trie (<em>skewness<\/em>) et d&rsquo;aplatissement (<em>kurtosis<\/em>), est adapt\u00e9 aux \u00e9chantillons de taille plus cons\u00e9quente.<\/p>\n<p>Pour affiner ces diagnostics, il est recommand\u00e9 de substituer aux r\u00e9sidus bruts les r\u00e9sidus studentis\u00e9s (<em>studentized residuals<\/em>). En effet, en raison de la projection g\u00e9om\u00e9trique exerc\u00e9e par la matrice chapeau <em>H = X(X<sup>T<\/sup>X)<sup>-1<\/sup>X<sup>T<\/sup><\/em>, la variance des r\u00e9sidus bruts individuels n&rsquo;est pas strictement identique pour toutes les observations : elle est modul\u00e9e par le levier <em>h<sub>ii<\/sub><\/em> exerc\u00e9 par chaque point dans l&rsquo;espace des pr\u00e9dicteurs selon la formule :<\/p>\n<p><em>Var(e<sub>i<\/sub>) = &sigma;<sup>2<\/sup>(1 &#8211; h<sub>ii<\/sub>)<\/em><\/p>\n<p>Les r\u00e9sidus studentis\u00e9s suppriment cette d\u00e9pendance m\u00e9canique en divisant chaque r\u00e9sidu par son \u00e9cart-type individuel exact estim\u00e9 sans l&rsquo;observation consid\u00e9r\u00e9e (r\u00e9sidus studentis\u00e9s supprim\u00e9s ext\u00e9rieurement, accessibles dans Statsmodels via <code>resultats.get_influence().resid_studentized_external<\/code>). Ces r\u00e9sidus suivent rigoureusement une loi de Student \u00e0 <em>n &#8211; p &#8211; 2<\/em> degr\u00e9s de libert\u00e9, permettant d&rsquo;identifier sans ambigu\u00eft\u00e9 les d\u00e9viations significatives exc\u00e9dant l&rsquo;intervalle critique [-2, +2] ou [-3, +3].<\/p>\n<h3>9.3 Identification des observations aberrantes et points de levier<\/h3>\n<p>La somme des carr\u00e9s des r\u00e9sidus, par construction math\u00e9matique quadratique, pr\u00e9sente une vuln\u00e9rabilit\u00e9 extr\u00eame face aux observations atypiques situ\u00e9es dans les zones marginales de l&rsquo;\u00e9chantillon. Un unique point d&rsquo;observation anormalement \u00e9loign\u00e9 de la trajectoire moyenne g\u00e9n\u00e9rale peut exercer un effet d&rsquo;attraction d\u00e9mesur\u00e9 sur la droite de r\u00e9gression, diminuant artificiellement le r\u00e9sidu associ\u00e9 \u00e0 cette observation tout en gonflant drastiquement les r\u00e9sidus de la majorit\u00e9 des autres donn\u00e9es conformes.<\/p>\n<p>Pour d\u00e9composer et neutraliser ces anomalies, le statisticien doit discriminer deux ph\u00e9nom\u00e8nes conjoints :<\/p>\n<ul>\n<li><strong>Le levier (leverage) :<\/strong> Mesure \u00e0 quel point la combinaison des valeurs des covariables d&rsquo;une observation <em>x<sub>i<\/sub><\/em> est atypique par rapport au centro\u00efde des variables explicatives. Il est quantifi\u00e9 par les \u00e9l\u00e9ments diagonaux <em>h<sub>ii<\/sub><\/em> de la matrice de projection orthogonale (matrice chapeau <em>H<\/em>). Une valeur de levier exc\u00e9dant le seuil th\u00e9orique de <em>2(p + 1) \/ n<\/em> d\u00e9signe une observation \u00e0 fort potentiel perturbateur.<\/li>\n<li><strong>L&rsquo;influence globale :<\/strong> Une observation est dite influente lorsqu&rsquo;elle combine \u00e0 la fois un levier \u00e9lev\u00e9 et un r\u00e9sidu studentis\u00e9 de forte amplitude, alt\u00e9rant substantiellement les estimations param\u00e9triques &beta; lorsque cette donn\u00e9e est extraite de l&rsquo;\u00e9chantillon.<\/li>\n<\/ul>\n<p>L&rsquo;indice quantitatif le plus universel pour synth\u00e9tiser cette influence unitaire est la <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Cook%27s_distance\">distance de Cook<\/a> (<em>Cook&rsquo;s Distance<\/em>), formul\u00e9e pour la <em>i<\/em>-\u00e8me observation par l&rsquo;expression :<\/p>\n<p><em>D<sub>i<\/sub> = [ &sum;<sub>j=1<\/sub><sup>n<\/sup> (&ycirc;<sub>j<\/sub> &#8211; &ycirc;<sub>j(i)<\/sub>)<sup>2<\/sup> ] \/ [ (p + 1) &times; s<sup>2<\/sup> ]<\/em><\/p>\n<p>o\u00f9 <em>&ycirc;<sub>j(i)<\/sub><\/em> repr\u00e9sente la pr\u00e9diction obtenue pour la <em>j<\/em>-\u00e8me observation lorsque le mod\u00e8le a \u00e9t\u00e9 r\u00e9estim\u00e9 en excluant l&rsquo;observation <em>i<\/em>. Statsmodels permet d&rsquo;extraire la totalit\u00e9 des distances de Cook associ\u00e9es \u00e0 un mod\u00e8le ajust\u00e9 par l&rsquo;interm\u00e9diaire de la m\u00e9thode <code>influence = resultats.get_influence()<\/code>, puis <code>cooks_d = influence.cooks_distance[0]<\/code>. Les observations dont la distance de Cook exc\u00e8de la valeur <em>4 \/ n<\/em> ou le seuil absolu de 1,0 doivent faire l&rsquo;objet d&rsquo;un examen approfondi.<\/p>\n<p>Parall\u00e8lement, la statistique PRESS (<em>Prediction Error Sum of Squares<\/em>) mesure la somme des carr\u00e9s des r\u00e9sidus d&rsquo;exclusion obtenus par une proc\u00e9dure de validation crois\u00e9e unitaire (<em>leave-one-out<\/em>), fournissant une \u00e9valuation non biais\u00e9e de la sensibilit\u00e9 pr\u00e9dictive de la RSS. Dans les sciences du comportement, l&rsquo;\u00e9limination d\u00e9finitive d&rsquo;une observation aberrante ne doit jamais \u00eatre automatique ou motiv\u00e9e par la seule volont\u00e9 d&rsquo;abaisser artificiellement la RSS : elle exige une justification substantielle clinique, physiologique ou m\u00e9thodologique (faute technique de recueil, tricherie av\u00e9r\u00e9e, incident mat\u00e9riel lors du protocole).<\/p>\n<h2>10. S\u00e9lection de mod\u00e8les et m\u00e9triques d&rsquo;information fond\u00e9es sur la RSS<\/h2>\n<h3>10.1 Calcul du coefficient de d\u00e9termination R\u00b2 \u00e0 partir de la RSS<\/h3>\n<p>Bien que la RSS mesure la variabilit\u00e9 inexpliqu\u00e9e absolue, la communaut\u00e9 scientifique privil\u00e9gie fr\u00e9quemment une m\u00e9trique adimensionnelle normalis\u00e9e comprise entre 0 et 1 : le coefficient de d\u00e9termination lin\u00e9aire standardis\u00e9, universellement d\u00e9sign\u00e9 par le symbole <em>R<sup>2<\/sup><\/em>. Sa d\u00e9finition formelle s&rsquo;appuie directement sur le ratio de la RSS et de la somme totale des carr\u00e9s (<em>TSS<\/em>) :<\/p>\n<p><em>R<sup>2<\/sup> = 1 &#8211; (RSS \/ TSS) = 1 &#8211; [ &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &#8211; &ycirc;<sub>i<\/sub>)<sup>2<\/sup> \/ &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &#8211; &ymacr;)<sup>2<\/sup> ]<\/em><\/p>\n<p>Cette \u00e9criture met en exergue l&rsquo;essence conceptuelle du <em>R<sup>2<\/sup><\/em> : il quantifie la proportion de la variance totale empirique qui se trouve effectivement r\u00e9sorb\u00e9e et expliqu\u00e9e par la pr\u00e9sence des pr\u00e9dicteurs au sein de l&rsquo;hyperplan de r\u00e9gression. Si la RSS s&rsquo;annule, l&rsquo;ajustement est absolu et le <em>R<sup>2<\/sup><\/em> culmine \u00e0 1,0. Si la RSS est \u00e9quivalente \u00e0 la dispersion totale sans mod\u00e8le (TSS), l&rsquo;ensemble des pr\u00e9dicteurs n&rsquo;apporte aucune information discriminante et le <em>R<sup>2<\/sup><\/em> chute \u00e0 0,0.<\/p>\n<p>Cependant, le <em>R<sup>2<\/sup><\/em> h\u00e9rite d&rsquo;une limitation math\u00e9matique intrins\u00e8que de la RSS : l&rsquo;adjonction syst\u00e9matique d&rsquo;une variable explicative arbitraire dans la matrice de design \u2014 quand bien m\u00eame ce nouveau r\u00e9gresseur ne serait qu&rsquo;une suite de nombres al\u00e9atoires sans lien th\u00e9orique avec la cible \u2014 provoque in\u00e9vitablement une diminution stochastique de la RSS (ou au pire sa stricte conservation), entra\u00eenant une hausse m\u00e9canique et artificielle du <em>R<sup>2<\/sup><\/em>.<\/p>\n<p>Pour neutraliser cette inflation trompeuse et p\u00e9naliser l&rsquo;adjonction st\u00e9rile de param\u00e8tres superflus, les statisticiens emploient le coefficient de d\u00e9termination ajust\u00e9 (<em>R<sup>2<\/sup> ajust\u00e9<\/em> ou <em>R&macr;<sup>2<\/sup><\/em>), qui pond\u00e8re chaque somme de carr\u00e9s par ses degr\u00e9s de libert\u00e9 respectifs :<\/p>\n<p><em>R&macr;<sup>2<\/sup> = 1 &#8211; [ (RSS \/ (n &#8211; p &#8211; 1)) \/ (TSS \/ (n &#8211; 1)) ] = 1 &#8211; [ (1 &#8211; R<sup>2<\/sup>) &times; (n &#8211; 1) \/ (n &#8211; p &#8211; 1) ]<\/em><\/p>\n<p>Contrairement au <em>R<sup>2<\/sup><\/em> brut, le <em>R<sup>2<\/sup><\/em> ajust\u00e9 ne cro\u00eet lors de l&rsquo;int\u00e9gration d&rsquo;un nouveau pr\u00e9dicteur que si la d\u00e9croissance marginale de la RSS compense la perte d&rsquo;un degr\u00e9 de libert\u00e9, ce qui correspond exactement \u00e0 l&rsquo;exigence que la statistique <em>F<\/em> associ\u00e9e \u00e0 cette nouvelle variable exc\u00e8de la valeur unitaire.<\/p>\n<h3>10.2 Crit\u00e8res d&rsquo;information d&rsquo;Akaike (AIC) et bay\u00e9sien (BIC)<\/h3>\n<p>Dans les contextes modernes de s\u00e9lection automatis\u00e9e de variables ou de comparaison de mod\u00e8les non n\u00e9cessairement embo\u00eet\u00e9s, les crit\u00e8res d&rsquo;information fond\u00e9s sur la th\u00e9orie de la communication et l&rsquo;inf\u00e9rence bay\u00e9sienne ont supplant\u00e9 le simple examen de la RSS. Ces crit\u00e8res arbitrent explicitement le compromis fondamental entre la qualit\u00e9 de l&rsquo;ajustement empirique (refl\u00e9t\u00e9e par la diminution de la RSS) et la complexit\u00e9 structurelle de l&rsquo;architecture param\u00e9trique.<\/p>\n<p>Sous l&rsquo;hypoth\u00e8se de normalit\u00e9 homosc\u00e9dastique des erreurs, le crit\u00e8re d&rsquo;information d&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Akaike_information_criterion\">Akaike (AIC)<\/a> se d\u00e9rive math\u00e9matiquement \u00e0 partir de la RSS sous une forme condens\u00e9e :<\/p>\n<p><em>AIC = n &times; ln(RSS \/ n) + 2k + C<\/em><\/p>\n<p>o\u00f9 <em>n<\/em> est la taille effective de l&rsquo;\u00e9chantillon, <em>k = p + 1<\/em> est le nombre total de param\u00e8tres estim\u00e9s (incluant la constante et la variance d&rsquo;erreur), et <em>C<\/em> d\u00e9signe une constante d&rsquo;\u00e9chelle additive invariante. L&rsquo;AIC quantifie la perte relative d&rsquo;information th\u00e9orique subie lors de la substitution du mod\u00e8le simplifi\u00e9 au processus physique g\u00e9n\u00e9rateur originel (divergence de Kullback-Leibler). Le meilleur mod\u00e8le au sens d&rsquo;Akaike est celui qui minimise la valeur scalaire de l&rsquo;AIC.<\/p>\n<p>Le crit\u00e8re d&rsquo;information bay\u00e9sien (BIC ou crit\u00e8re de Schwarz) applique une p\u00e9nalit\u00e9 structurelle sensiblement plus agressive \u00e0 la complexit\u00e9 des param\u00e8tres d\u00e8s lors que la taille d&rsquo;\u00e9chantillon <em>n<\/em> exc\u00e8de la valeur exp(2) &approx; 7,39 :<\/p>\n<p><em>BIC = n &times; ln(RSS \/ n) + k &times; ln(n) + C<\/em><\/p>\n<p>La p\u00e9nalit\u00e9 proportionnelle \u00e0 <em>k &times; ln(n)<\/em> conf\u00e8re au BIC une consistance asymptotique : lorsque <em>n &rarr; &infin;<\/em>, la probabilit\u00e9 que le BIC s\u00e9lectionne la v\u00e9ritable architecture parcimonieuse g\u00e9n\u00e9ratrice converge vers 1,0, ce qui en fait l&rsquo;outil d&rsquo;arbitrage privil\u00e9gi\u00e9 dans les recherches psychom\u00e9triques fondamentales o\u00f9 l&rsquo;on cherche \u00e0 isoler un r\u00e9seau restreint de d\u00e9terminants comportementaux sans tol\u00e9rer de bruit parasite.<\/p>\n<h3>10.3 Test de rapport de vraisemblance et comparaison F<\/h3>\n<p>Lorsque les architectures compar\u00e9es sont formellement embo\u00eet\u00e9es \u2014 c&rsquo;est-\u00e0-dire que le mod\u00e8le restreint <em>M<sub>R<\/sub><\/em> d\u00e9rive du mod\u00e8le complet <em>M<sub>C<\/sub><\/em> par la mise \u00e0 z\u00e9ro d&rsquo;un sous-ensemble de <em>q<\/em> coefficients r\u00e9gresseurs \u2014, la confrontation statistique directe repose sur l&rsquo;\u00e9valuation de la r\u00e9duction proportionnelle de la RSS rapport\u00e9e \u00e0 la variance r\u00e9siduelle non biais\u00e9e du mod\u00e8le satur\u00e9.<\/p>\n<p>La statistique du test <em>F<\/em> partiel se formalise rigoureusement par :<\/p>\n<p><em>F<sub>obs<\/sub> = [ (RSS<sub>R<\/sub> &#8211; RSS<sub>C<\/sub>) \/ q ] \/ [ RSS<sub>C<\/sub> \/ (n &#8211; p<sub>C<\/sub> &#8211; 1) ]<\/em><\/p>\n<p>Sous l&rsquo;hypoth\u00e8se nulle <em>H<sub>0<\/sub><\/em> selon laquelle les <em>q<\/em> variables additionnelles poss\u00e8dent des coefficients r\u00e9els conjointement nuls dans la population, cette statistique suit exactement une distribution th\u00e9orique de Fisher-Snedecor \u00e0 <em>(q, n &#8211; p<sub>C<\/sub> &#8211; 1)<\/em> degr\u00e9s de libert\u00e9.<\/p>\n<p>Sous Python, l&rsquo;ex\u00e9cution automatis\u00e9e de cette proc\u00e9dure s&rsquo;impl\u00e9mente au moyen de la m\u00e9thode <code>anova_lm<\/code> issue du sous-module <code>statsmodels.stats.anova<\/code>, qui prend en arguments deux mod\u00e8les OLS hi\u00e9rarchiquement ajust\u00e9s :<\/p>\n<p><code>from statsmodels.stats.anova import anova_lm<\/code><br \/>\n<code>tableau_anova = anova_lm(modele_restreint_fit, modele_complet_fit)<\/code><br \/>\n<code>print(tableau_anova)<\/code><\/p>\n<p>Le tableau r\u00e9sultant consigne explicitement les valeurs des RSS des deux mod\u00e8les (colonne <code>ssr<\/code>), la diff\u00e9rence marginale <code>ss_diff<\/code>, les degr\u00e9s de libert\u00e9 associ\u00e9s et la valeur de p correspondante (<code>Pr(&gt;F)<\/code>). Si ce niveau de risque d&rsquo;erreur de premi\u00e8re esp\u00e8ce <em>p<\/em> s&rsquo;av\u00e8re inf\u00e9rieur au seuil de signification &alpha; = 0,05, le chercheur conclut formellement que la r\u00e9duction de la RSS op\u00e9r\u00e9e par le mod\u00e8le \u00e9tendu est statistiquement significative, attestant de l&rsquo;apport explicatif d\u00e9terminant des nouveaux pr\u00e9dicteurs introduits.<\/p>\n<h2>11. \u00c9tude de cas int\u00e9gr\u00e9e : Mod\u00e9lisation de l&rsquo;anxi\u00e9t\u00e9 et de la performance acad\u00e9mique<\/h2>\n<h3>11.1 Protocole exp\u00e9rimental et caract\u00e9ristiques de l&rsquo;\u00e9chantillon<\/h3>\n<p>Afin de concr\u00e9tiser l&rsquo;ensemble de ces d\u00e9veloppements alg\u00e9briques et programmatiques, analysons une \u00e9tude empirique int\u00e9gr\u00e9e issue de la psychologie cognitive et de l&rsquo;ing\u00e9nierie \u00e9ducative. L&rsquo;investigation porte sur une cohorte de <em>n = 150<\/em> \u00e9tudiants universitaires inscrits en cycle sup\u00e9rieur scientifique, ayant consenti \u00e0 participer \u00e0 un protocole d&rsquo;\u00e9valuation psychom\u00e9trique standardis\u00e9 durant la p\u00e9riode de r\u00e9vision pr\u00e9c\u00e9dant leurs examens terminaux.<\/p>\n<p>Le protocole instrumente trois dimensions explicatives continues :<\/p>\n<ul>\n<li><code>stress_percu<\/code> : Score continu d&rsquo;anxi\u00e9t\u00e9 situationnelle mesur\u00e9 par l&rsquo;\u00e9chelle psychom\u00e9trique valid\u00e9e PSS-10 (<em>Perceived Stress Scale<\/em>), standardis\u00e9 de 0 (absence de stress per\u00e7u) \u00e0 50 (d\u00e9tresse psychologique s\u00e9v\u00e8re).<\/li>\n<li><code>heures_sommeil<\/code> : Nombre moyen d&rsquo;heures de sommeil nocturne effectif mesur\u00e9 par actim\u00e9trie continue durant la quinzaine pr\u00e9c\u00e9dant l&rsquo;\u00e9preuve d&rsquo;\u00e9valuation.<\/li>\n<li><code>temps_revision<\/code> : Volume total d&rsquo;heures consacr\u00e9es \u00e0 l&rsquo;apprentissage et \u00e0 la r\u00e9solution d&rsquo;exercices cognitifs durant la m\u00eame p\u00e9riode.<\/li>\n<\/ul>\n<p>La variable crit\u00e8re \u00e0 mod\u00e9liser, <code>score_examen<\/code>, correspond \u00e0 la note objectivement obtenue \u00e0 l&rsquo;\u00e9preuve de synth\u00e8se de fin d&rsquo;ann\u00e9e, calibr\u00e9e sur une distribution continue de moyenne empirique vis\u00e9e &mu; = 70 et d&rsquo;\u00e9cart-type &sigma; = 15. L&rsquo;objectif analytique consiste \u00e0 quantifier l&rsquo;ajustement du mod\u00e8le de r\u00e9gression lin\u00e9aire multivari\u00e9e, \u00e0 calculer la RSS par les trois approches programmatiques document\u00e9es, et \u00e0 v\u00e9rifier la conformit\u00e9 des estimations.<\/p>\n<h3>11.2 Calcul comparatif et synchronis\u00e9 de la RSS via les trois m\u00e9thodes<\/h3>\n<p>Conduisons de fa\u00e7on conjointe les trois proc\u00e9dures algorithmiques sur ce jeu de donn\u00e9es exp\u00e9rimental afin de corroborer formellement l&rsquo;identit\u00e9 math\u00e9matique de leurs sorties num\u00e9riques. Le script Python ci-apr\u00e8s formalise la synth\u00e8se des d\u00e9marches pr\u00e9sent\u00e9es :<\/p>\n<p><code>import numpy as np<\/code><br \/>\n<code>import pandas as pd<\/code><br \/>\n<code>import statsmodels.api as sm<\/code><br \/>\n<code>from sklearn.linear_model import LinearRegression<\/code><br \/>\n<code>from sklearn.metrics import mean_squared_error<\/code><br \/>\n<code>import time<\/code><\/p>\n<p><code># Simulation d'un \u00e9chantillon psychom\u00e9trique r\u00e9aliste (n = 150)<\/code><br \/>\n<code>rng = np.random.default_rng(seed=101)<\/code><br \/>\n<code>n_sujets = 150<\/code><br \/>\n<code>stress = rng.normal(loc=25, scale=6, size=n_sujets)<\/code><br \/>\n<code>sommeil = rng.normal(loc=7.2, scale=1.1, size=n_sujets)<\/code><br \/>\n<code>revision = rng.normal(loc=35, scale=8, size=n_sujets)<\/code><br \/>\n<code>bruit = rng.normal(loc=0, scale=8.5, size=n_sujets)<\/code><br \/>\n<code>score = 45.0 - (0.55 * stress) + (2.1 * sommeil) + (0.75 * revision) + bruit<\/code><br \/>\n<code>df_etude = pd.DataFrame({'stress': stress, 'sommeil': sommeil, 'revision': revision, 'score': score})<\/code><\/p>\n<p><code># Configuration des matrices<\/code><br \/>\n<code>X = df_etude[['stress', 'sommeil', 'revision']].values<\/code><br \/>\n<code>y = df_etude['score'].values<\/code><\/p>\n<p><code># 1. Approche Fondamentale NumPy<\/code><br \/>\n<code>t0 = time.perf_counter()<\/code><br \/>\n<code>X_mat = np.column_stack([np.ones(n_sujets), X])<\/code><br \/>\n<code>beta_np = np.linalg.solve(X_mat.T @ X_mat, X_mat.T @ y)<\/code><br \/>\n<code>y_pred_np = X_mat @ beta_np<\/code><br \/>\n<code>residus_np = y - y_pred_np<\/code><br \/>\n<code>rss_methode_1 = np.dot(residus_np, residus_np)<\/code><br \/>\n<code>duree_1 = time.perf_counter() - t0<\/code><\/p>\n<p><code># 2. Approche Statistique Statsmodels<\/code><br \/>\n<code>t0 = time.perf_counter()<\/code><br \/>\n<code>X_sm = sm.add_constant(df_etude[['stress', 'sommeil', 'revision']])<\/code><br \/>\n<code>res_sm = sm.OLS(df_etude['score'], X_sm).fit()<\/code><br \/>\n<code>rss_methode_2 = res_sm.ssr<\/code><br \/>\n<code>duree_2 = time.perf_counter() - t0<\/code><\/p>\n<p><code># 3. Approche Machine Learning Scikit-Learn<\/code><br \/>\n<code>t0 = time.perf_counter()<\/code><br \/>\n<code>reg_skl = LinearRegression(fit_intercept=True)<\/code><br \/>\n<code>reg_skl.fit(X, y)<\/code><br \/>\n<code>y_pred_skl = reg_skl.predict(X)<\/code><br \/>\n<code>rss_methode_3 = mean_squared_error(y, y_pred_skl) * len(y)<\/code><br \/>\n<code>duree_3 = time.perf_counter() - t0<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de cette routine comparative met en \u00e9vidence des r\u00e9sultats instructifs r\u00e9sum\u00e9s dans le tableau synth\u00e9tique ci-dessous :<\/p>\n<ul>\n<li><strong>M\u00e9thode 1 (NumPy vectoris\u00e9) :<\/strong> RSS = 10 542,3817342 \u2014 Temps d&rsquo;ex\u00e9cution : ~0,08 ms. Pr\u00e9cision machine absolue, allocation m\u00e9moire minimale sans cr\u00e9ation d&rsquo;objets interm\u00e9diaires.<\/li>\n<li><strong>M\u00e9thode 2 (Statsmodels OLS) :<\/strong> RSS = 10 542,3817342 \u2014 Temps d&rsquo;ex\u00e9cution : ~1,85 ms. Concordance num\u00e9rique stricte jusqu&rsquo;\u00e0 la treizi\u00e8me d\u00e9cimale. Temps de calcul plus \u00e9lev\u00e9 d\u00fb \u00e0 la d\u00e9rivation de l&rsquo;ensemble des matrices de covariance, des tests de studentisation et des statistiques de Fisher.<\/li>\n<li><strong>M\u00e9thode 3 (Scikit-Learn MSE) :<\/strong> RSS = 10 542,3817342 \u2014 Temps d&rsquo;ex\u00e9cution : ~0,42 ms. Identit\u00e9 arithm\u00e9tique rigoureuse apr\u00e8s mise \u00e0 l&rsquo;\u00e9chelle par <em>N<\/em>. Excellente rapidit\u00e9 algorithmique adapt\u00e9e \u00e0 l&rsquo;int\u00e9gration dans des cha\u00eenes de traitement intensives.<\/li>\n<\/ul>\n<p>Cette concordance absolue d\u00e9montre que les trois architectures reposent sur les m\u00eames fondements de d\u00e9composition orthogonale en alg\u00e8bre lin\u00e9aire, chaque biblioth\u00e8que se distinguant simplement par le niveau d&rsquo;abstraction m\u00e9trologique et fonctionnel qu&rsquo;elle propose.<\/p>\n<h3>11.3 Interpr\u00e9tation substantive des r\u00e9sultats psychom\u00e9triques<\/h3>\n<p>L&rsquo;interpr\u00e9tation de la RSS obtenue (10 542,38) acquiert sa pleine signification lorsqu&rsquo;elle est rapport\u00e9e aux m\u00e9triques dimensionn\u00e9es du construit \u00e9tudi\u00e9. La taille de notre \u00e9chantillon \u00e9tant de <em>n = 150<\/em> sujets et le mod\u00e8le int\u00e9grant <em>p = 3<\/em> r\u00e9gresseurs plus l&rsquo;interception (soit 4 param\u00e8tres au total), les degr\u00e9s de libert\u00e9 r\u00e9siduels s&rsquo;\u00e9tablissent \u00e0 :<\/p>\n<p><em>df<sub>resid<\/sub> = 150 &#8211; 3 &#8211; 1 = 146<\/em><\/p>\n<p>L&rsquo;estimateur non biais\u00e9 de la variance r\u00e9siduelle stochastique s&rsquo;\u00e9l\u00e8ve donc \u00e0 :<\/p>\n<p><em>s<sup>2<\/sup> = RSS \/ df<sub>resid<\/sub> = 10 542,38 \/ 146 &approx; 72,21<\/em><\/p>\n<p>En extrayant la racine carr\u00e9e de cette variance, nous obtenons l&rsquo;erreur standard d&rsquo;estimation du mod\u00e8le :<\/p>\n<p><em>s = &radic;(72,21) &approx; 8,50 points<\/em><\/p>\n<p>Ce r\u00e9sultat fournit une lecture substantielle directe de la pr\u00e9cision de l&rsquo;\u00e9valuation : l&rsquo;\u00e9cart-type de dispersion des performances cognitives individuelles autour de la surface de r\u00e9gression pr\u00e9dite par le niveau de stress, la quantit\u00e9 de sommeil et le temps de r\u00e9vision est de 8,50 points sur 100. En appliquant les propri\u00e9t\u00e9s de la distribution normale (sous r\u00e9serve que les diagnostics de normalit\u00e9 soient satisfaits), environ 95 % des scores r\u00e9els d&rsquo;examen observ\u00e9s se situeront dans un intervalle de pr\u00e9diction encadr\u00e9 par &plusmn; 1,96 &times; 8,50 &approx; &plusmn; 16,66 points autour de la projection th\u00e9orique personnalis\u00e9e.<\/p>\n<p>Si la dispersion totale observ\u00e9e initiale (TSS) de la cohorte \u00e9tait de 28 500,00 points de carr\u00e9s, la r\u00e9duction permise par notre mod\u00e8le est substantielle :<\/p>\n<p><em>R<sup>2<\/sup> = 1 &#8211; (10 542,38 \/ 28 500,00) &approx; 1 &#8211; 0,37 = 0,63<\/em><\/p>\n<p>Le mod\u00e8le capture ainsi 63 % de la variance de la performance universitaire, les 37 % r\u00e9siduels inexpliqu\u00e9s (repr\u00e9sent\u00e9s par notre RSS) \u00e9tant imputables \u00e0 des construits cognitifs non mod\u00e9lis\u00e9s dans cette \u00e9tude (tels que l&rsquo;efficience de la m\u00e9moire de travail, l&rsquo;intelligence cristallis\u00e9e pr\u00e9alable, ou des al\u00e9as d&rsquo;\u00e9valuation non contr\u00f4l\u00e9s lors de la passation de l&rsquo;\u00e9preuve).<\/p>\n<h2>12. Pi\u00e8ges m\u00e9thodologiques, limites de la RSS et bonnes pratiques<\/h2>\n<h3>12.1 Sensibilit\u00e9 accrue aux valeurs aberrantes<\/h3>\n<p>L&rsquo;une des vuln\u00e9rabilit\u00e9s m\u00e9thodologiques les plus document\u00e9es de la somme des carr\u00e9s des r\u00e9sidus tient \u00e0 sa nature m\u00eame de fonction de perte quadratique. L&rsquo;\u00e9l\u00e9vation de chaque \u00e9cart <em>e<sub>i<\/sub><\/em> \u00e0 la puissance 2 engendre une croissance exponentielle de la p\u00e9nalit\u00e9 au fur et \u00e0 mesure que l&rsquo;on s&rsquo;\u00e9carte de la valeur pr\u00e9dite. Ainsi, un r\u00e9sidu individuel \u00e9gal \u00e0 10 contribue \u00e0 hauteur de 100 \u00e0 la RSS globale, tandis qu&rsquo;un r\u00e9sidu \u00e9gal \u00e0 30 \u2014 repr\u00e9sentant un \u00e9cart trois fois plus grand \u2014 contribue \u00e0 hauteur de 900, exer\u00e7ant un poids multiplicateur relatif d\u00e9mesur\u00e9 de 9 sur la fonction d&rsquo;optimisation.<\/p>\n<p>Cette sensibilit\u00e9 conf\u00e8re \u00e0 l&rsquo;estimateur des MCO une tr\u00e8s faible r\u00e9sistance aux contaminations distributionnelles : un seul enregistrement aberrant dans une base de donn\u00e9es peut faire pivoter l&rsquo;hyperplan de r\u00e9gression pour r\u00e9duire son propre r\u00e9sidu quadratique extr\u00eame, au d\u00e9triment de l&rsquo;ad\u00e9quation g\u00e9n\u00e9rale avec les donn\u00e9es r\u00e9guli\u00e8res. En statistique robuste, on caract\u00e9rise cette fragilit\u00e9 par un point de rupture (<em>breakdown point<\/em>) rigoureusement \u00e9gal \u00e0 0 % : il suffit d&rsquo;une seule observation contamin\u00e9e sur un \u00e9chantillon infini pour alt\u00e9rer arbitrairement la valeur estim\u00e9e des param\u00e8tres.<\/p>\n<p>Lorsque la pr\u00e9sence de valeurs atypiques ou de perturbations distributionnelles \u00e0 queues lourdes est suspect\u00e9e dans les donn\u00e9es exp\u00e9rimentales, l&rsquo;analyste avis\u00e9 doit compl\u00e9ter l&rsquo;examen de la RSS par des m\u00e9triques plus stables, telles que l&rsquo;erreur absolue m\u00e9diane (<em>Median Absolute Error<\/em>), ou envisager la mod\u00e9lisation par r\u00e9gression robuste. Les estimateurs fond\u00e9s sur la perte de Huber (impl\u00e9ment\u00e9s dans Statsmodels via <code>sm.RLM<\/code>) op\u00e8rent un compromis dynamique \u00e9l\u00e9gant en traitant les r\u00e9sidus de faible amplitude selon une fonction quadratique tout en p\u00e9nalisant les r\u00e9sidus extr\u00eames de fa\u00e7on strictement lin\u00e9aire (L1), pr\u00e9servant ainsi la RSS de distorsions d\u00e9l\u00e9t\u00e8res.<\/p>\n<h3>12.2 Risques de surapprentissage li\u00e9s \u00e0 la minimisation pure de la RSS<\/h3>\n<p>La poursuite inconditionnelle de la minimisation de la RSS sur l&rsquo;\u00e9chantillon d&rsquo;entra\u00eenement constitue la cause majeure du surapprentissage (<em>overfitting<\/em>) en statistique appliqu\u00e9e et en apprentissage automatique. Comme le stipule la th\u00e9orie g\u00e9om\u00e9trique de l&rsquo;ajustement lin\u00e9aire, l&rsquo;int\u00e9gration syst\u00e9matique de nouveaux r\u00e9gresseurs dans un mod\u00e8le math\u00e9matique augmente le nombre de dimensions du sous-espace de projection, contraignant la RSS \u00e0 d\u00e9cro\u00eetre de fa\u00e7on monotone :<\/p>\n<p><em>RSS(p + 1) &le; RSS(p)<\/em><\/p>\n<p>Si le nombre de param\u00e8tres estim\u00e9s <em>p<\/em> approche la taille de l&rsquo;\u00e9chantillon <em>n<\/em> (situation o\u00f9 <em>p &rarr; n<\/em>), la RSS s&rsquo;effondre vers z\u00e9ro. Dans le cas d&rsquo;interpolation exacte o\u00f9 <em>p = n &#8211; 1<\/em>, la RSS s&rsquo;annule rigoureusement et le mod\u00e8le passe exactement par l&rsquo;ensemble des points d&rsquo;observation. Cependant, cette perfection apparente masque un \u00e9chec \u00e9pist\u00e9mologique absolu : le mod\u00e8le n&rsquo;a rien appris des r\u00e9gularit\u00e9s sous-jacentes du ph\u00e9nom\u00e8ne stochastique, il a simplement m\u00e9moris\u00e9 le bruit idiosyncrasique de l&rsquo;\u00e9chantillon.<\/p>\n<p>Pour contrer cette d\u00e9rive, la minimisation empirique de la RSS doit \u00eatre assujettie \u00e0 des contraintes de r\u00e9gularisation param\u00e9trique p\u00e9nalisant la complexit\u00e9 de l&rsquo;espace des coefficients :<\/p>\n<ul>\n<li><strong>R\u00e9gression Ridge (norme L2) :<\/strong> Minimise <em>RSS + &alpha; &sum;<sub>j=1<\/sub><sup>p<\/sup> &beta;<sub>j<\/sub><sup>2<\/sup><\/em>, contraignant les coefficients \u00e0 se contracter de fa\u00e7on continue vers z\u00e9ro afin de ma\u00eetriser la variance d&rsquo;\u00e9chantillonnage.<\/li>\n<li><strong>R\u00e9gression Lasso (norme L1) :<\/strong> Minimise <em>RSS + &alpha; &sum;<sub>j=1<\/sub><sup>p<\/sup> |&beta;<sub>j<\/sub>|<\/em>, r\u00e9alisant une s\u00e9lection parcimonieuse automatique en for\u00e7ant certains coefficients r\u00e9gresseurs non contributifs \u00e0 s&rsquo;annuler exactement.<\/li>\n<\/ul>\n<p>Il est imp\u00e9ratif d&rsquo;\u00e9tablir une distinction \u00e9tanche entre la <em>RSS d&rsquo;entra\u00eenement<\/em> (calcul\u00e9e sur les donn\u00e9es ayant servi \u00e0 l&rsquo;estimation) et la <em>RSS de g\u00e9n\u00e9ralisation<\/em> ou de test (calcul\u00e9e sur un \u00e9chantillon ind\u00e9pendant). Seule la minimisation de cette derni\u00e8re sur des donn\u00e9es jamais pr\u00e9sent\u00e9es au mod\u00e8le constitue une preuve empirique recevable de qualit\u00e9 pr\u00e9dictive.<\/p>\n<h3>12.3 Synth\u00e8se des recommandations pour la mod\u00e9lisation statistique<\/h3>\n<p>Pour clore cette analyse m\u00e9thodologique, il est pertinent de formaliser une grille de contr\u00f4le op\u00e9rationnelle guidant le chercheur dans le traitement rigoureux de la RSS au sein de ses travaux statistiques sous Python :<\/p>\n<ul>\n<li><strong>V\u00e9rification syst\u00e9matique de l&rsquo;interception :<\/strong> Toujours s&rsquo;assurer que la matrice de design int\u00e8gre formellement une constante d&rsquo;interception unitaire (via <code>sm.add_constant<\/code> dans Statsmodels ou <code>fit_intercept=True<\/code> dans Scikit-Learn) pour garantir la d\u00e9composition orthogonale de la variance totale <em>TSS = ESS + RSS<\/em> et la nullit\u00e9 de l&rsquo;esp\u00e9rance math\u00e9matique des r\u00e9sidus.<\/li>\n<li><strong>Validation crois\u00e9e de l&rsquo;ad\u00e9quation :<\/strong> Ne jamais interpr\u00e9ter la magnitude brute de la RSS de fa\u00e7on isol\u00e9e. La convertir syst\u00e9matiquement en coefficient de d\u00e9termination ajust\u00e9 <em>R&macr;<sup>2<\/sup><\/em>, en erreur standard d&rsquo;estimation non biais\u00e9e <em>s<\/em>, et la confronter aux crit\u00e8res p\u00e9nalis\u00e9s d&rsquo;information AIC et BIC.<\/li>\n<li><strong>Ex\u00e9cution syst\u00e9matique de la batterie diagnostique :<\/strong> Accompagner tout calcul de r\u00e9gression d&rsquo;un audit distributionnel complet des r\u00e9sidus : inspection visuelle du diagramme r\u00e9sidus-pr\u00e9dictions, v\u00e9rification de l&rsquo;homosc\u00e9dasticit\u00e9 via le test de Breusch-Pagan, contr\u00f4le de normalit\u00e9 par Q-Q plot et test de Shapiro-Wilk, et d\u00e9tection des leviers d&rsquo;influence via la distance de Cook.<\/li>\n<li><strong>Choix raisonn\u00e9 de l&rsquo;outil informatique :<\/strong> Exploiter la compl\u00e9mentarit\u00e9 de l&rsquo;\u00e9cosyst\u00e8me scientifique Python en mobilisant NumPy pour les d\u00e9veloppements algorithmiques fondamentaux et les optimisations vectoris\u00e9es \u00e0 haute cadence, Statsmodels pour l&rsquo;inf\u00e9rence statistique confirmatoire et la production de rapports acad\u00e9miques certifi\u00e9s, et Scikit-Learn pour la validation crois\u00e9e modulaire et le d\u00e9ploiement de mod\u00e8les pr\u00e9dictifs r\u00e9gularis\u00e9s.<\/li>\n<li><strong>Transparence m\u00e9thodologique et tra\u00e7abilit\u00e9 du code :<\/strong> Consigner scrupuleusement l&rsquo;ensemble des scripts de nettoyage de donn\u00e9es, fixer formellement les graines pseudo-al\u00e9atoires des g\u00e9n\u00e9rateurs d&rsquo;\u00e9chantillons, et rapporter fid\u00e8lement dans les publications scientifiques les degr\u00e9s de libert\u00e9 r\u00e9siduels conjointement \u00e0 la valeur num\u00e9rique de la somme des carr\u00e9s des r\u00e9sidus.<\/li>\n<\/ul>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Akaike, H. (1974). A new look at the statistical model identification. <em>IEEE Transactions on Automatic Control<\/em>, 19(6), 716-723. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1109\/TAC.1974.1100705\">https:\/\/doi.org\/10.1109\/TAC.1974.1100705<\/a><\/li>\n<li>Breusch, T. S., &amp; Pagan, A. R. (1979). A simple test for heteroscedasticity and random coefficient variation. <em>Econometrica<\/em>, 47(5), 1287-1294. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1911963\">https:\/\/doi.org\/10.2307\/1911963<\/a><\/li>\n<li>Cook, R. D. (1977). Detection of influential observation in linear regression. <em>Technometrics<\/em>, 19(1), 15-18. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00401706.1977.10489493\">https:\/\/doi.org\/10.1080\/00401706.1977.10489493<\/a><\/li>\n<li>Fox, J. (2016). <em>Applied regression analysis and generalized linear models<\/em> (3rd ed.). SAGE Publications.<\/li>\n<li>Harris, C. R., Millman, K. J., van der Walt, S. J., Gommers, R., Virtanen, P., Cournapeau, D., &#8230; &amp; Oliphant, T. E. (2020). Array programming with NumPy. <em>Nature<\/em>, 585(7825), 357-362. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1038\/s41586-020-2649-2\">https:\/\/doi.org\/10.1038\/s41586-020-2649-2<\/a><\/li>\n<li>Hastie, T., Tibshirani, R., &amp; Friedman, J. (2009). <em>The elements of statistical learning: Data mining, inference, and prediction<\/em> (2nd ed.). Springer. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-0-387-84858-7\">https:\/\/doi.org\/10.1007\/978-0-387-84858-7<\/a><\/li>\n<li>McKinney, W. (2010). Data structures for statistical computing in Python. In <em>Proceedings of the 9th Python in Science Conference<\/em> (Vol. 445, pp. 56-61). <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.25080\/Majora-92bf1922-00a\">https:\/\/doi.org\/10.25080\/Majora-92bf1922-00a<\/a><\/li>\n<li>Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., &#8230; &amp; Duchesnay, \u00c9. (2011). Scikit-learn: Machine learning in Python. <em>Journal of Machine Learning Research<\/em>, 12, 2825-2830. <a rel=\"noopener noreferrer\" href=\"https:\/\/jmlr.org\/papers\/v12\/pedregosa11a.html\">https:\/\/jmlr.org\/papers\/v12\/pedregosa11a.html<\/a><\/li>\n<li>Schwarz, G. (1978). Estimating the dimension of a model. <em>The Annals of Statistics<\/em>, 6(2), 461-464. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aos\/1176344136\">https:\/\/doi.org\/10.1214\/aos\/1176344136<\/a><\/li>\n<li>Seabold, S., &amp; Perktold, J. (2010). statsmodels: Econometric and statistical modeling with Python. In <em>Proceedings of the 9th Python in Science Conference<\/em> (Vol. 57, pp. 92-96). <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.25080\/Majora-92bf1922-011\">https:\/\/doi.org\/10.25080\/Majora-92bf1922-011<\/a><\/li>\n<li>Virtanen, P., Gommers, R., Oliphant, T. E., Haberland, M., Reddy, T., Cournapeau, D., &#8230; &amp; SciPy 1.0 Contributors. (2020). SciPy 1.0: Fundamental algorithms for scientific computing in Python. <em>Nature Methods<\/em>, 17(3), 261-272. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1038\/s41592-019-0686-2\">https:\/\/doi.org\/10.1038\/s41592-019-0686-2<\/a><\/li>\n<li>Wooldridge, J. M. (2020). <em>Introductory econometrics: A modern approach<\/em> (7th ed.). Cengage Learning.<\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour calculer et interpr\u00e9ter la somme des carr\u00e9s des r\u00e9sidus (RSS) en Python appliqu\u00e9e \u00e0 la mod\u00e9lisation statistique.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2464,1838],"tags":[16,2468,1910,7847,2903],"class_list":["post-6088","psych_stat","type-psych_stat","status-publish","hentry","category-science-des-donnees","category-statistiques","tag-psychometrie","tag-python","tag-regression-lineaire","tag-somme-des-carres-des-residus","tag-statsmodels"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6088","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=6088"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6088\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=6088"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=6088"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=6088"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}