{"id":1286,"date":"2026-10-14T18:00:38","date_gmt":"2026-10-14T18:00:38","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-calculer-le-cp-de-mallows-dans-r\/"},"modified":"2026-10-14T18:00:38","modified_gmt":"2026-10-14T18:00:38","slug":"comment-calculer-le-cp-de-mallows-dans-r","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-calculer-le-cp-de-mallows-dans-r","title":{"rendered":"Comment calculer le Cp de Mallows dans R"},"content":{"rendered":"<p>La mod\u00e9lisation statistique contemporaine se trouve confront\u00e9e \u00e0 un d\u00e9fi \u00e9pist\u00e9mologique permanent : concilier l&rsquo;ajustement empirique aux donn\u00e9es observ\u00e9es avec la parcimonie structurelle indispensable \u00e0 la g\u00e9n\u00e9ralisation pr\u00e9dictive. Dans le cadre de la r\u00e9gression lin\u00e9aire multiple, l&rsquo;ajout syst\u00e9matique de covariables tend in\u00e9vitablement \u00e0 r\u00e9duire ou maintenir constante la somme des carr\u00e9s des r\u00e9sidus, procurant l&rsquo;illusion fallacieuse d&rsquo;une am\u00e9lioration continue de la qualit\u00e9 du mod\u00e8le. Cette inflation artificielle des performances sur l&rsquo;\u00e9chantillon d&rsquo;apprentissage masque un ph\u00e9nom\u00e8ne d\u00e9l\u00e9t\u00e8re bien connu des biostatisticiens et des psychom\u00e9triciens : le sur-apprentissage ou sur-ajustement (<em>overfitting<\/em>), caract\u00e9ris\u00e9 par une explosion de la variance des estimateurs et une d\u00e9gradation dramatique du pouvoir pr\u00e9dictif sur des donn\u00e9es ind\u00e9pendantes.<\/p>\n<p>Pour arbitrer rigoureusement ce compromis fondamental entre biais et variance, la th\u00e9orie statistique a vu \u00e9merger au d\u00e9but des ann\u00e9es 1970 plusieurs crit\u00e8res d&rsquo;information et de p\u00e9nalisation. Parmi ces instruments de d\u00e9cision algorithmique, le crit\u00e8re <strong>Cp de Mallows<\/strong> occupe une place matricielle. Con\u00e7u sp\u00e9cifiquement pour comparer des sous-mod\u00e8les d\u00e9riv\u00e9s d&rsquo;un mod\u00e8le complet pr\u00e9sum\u00e9 non biais\u00e9, le Cp \u00e9value l&rsquo;erreur quadratique totale moyenne de pr\u00e9diction normalis\u00e9e en p\u00e9nalisant explicitement la dimensionnalit\u00e9 de l&rsquo;espace des pr\u00e9dicteurs. Contrairement aux coefficients de d\u00e9termination classiques qui r\u00e9compensent aveugl\u00e9ment la complexit\u00e9 param\u00e9trique, le crit\u00e8re de Mallows offre une boussole analytique permettant d&rsquo;isoler l&rsquo;architecture optimale garantissant une estimation robuste et non biais\u00e9e des relations sous-jacentes.<\/p>\n<p>L&rsquo;impl\u00e9mentation de cette m\u00e9trique au sein de l&rsquo;environnement de calcul statistique <a rel=\"noopener noreferrer\" href=\"https:\/\/www.r-project.org\/\">R<\/a> constitue aujourd&rsquo;hui une comp\u00e9tence fondamentale pour tout chercheur ou analyste de donn\u00e9es op\u00e9rant en \u00e9conom\u00e9trie, en \u00e9pid\u00e9miologie ou en psychom\u00e9trie. Ce guide exhaustif explore les fondements math\u00e9matiques, les subtilit\u00e9s interpr\u00e9tatives et les protocoles logiciels n\u00e9cessaires pour calculer, visualiser et contextualiser le crit\u00e8re Cp de Mallows \u00e0 l&rsquo;aide de routines R \u00e9prouv\u00e9es, allant des biblioth\u00e8ques sp\u00e9cialis\u00e9es comme <em>olsrr<\/em> et <em>leaps<\/em> jusqu&rsquo;\u00e0 l&rsquo;\u00e9criture d&rsquo;algorithmes vectoriels autonomes.<\/p>\n<h2>1. Introduction th\u00e9orique au crit\u00e8re Cp de Mallows en mod\u00e9lisation statistique<\/h2>\n<h3>1.1 Origine historique et place dans l&rsquo;\u00e9conom\u00e9trie et la psychom\u00e9trie<\/h3>\n<p>Propos\u00e9 formellement par le statisticien britannique Colin Lingwood Mallows dans son article s\u00e9minal publi\u00e9 en 1973 dans la revue <em>Technometrics<\/em> (\u00ab Some Comments on Cp \u00bb), le crit\u00e8re Cp est n\u00e9 des r\u00e9flexions entourant la s\u00e9lection de variables dans les contextes de r\u00e9gression lin\u00e9aire gaussienne. \u00c0 une \u00e9poque o\u00f9 les capacit\u00e9s de calcul informatique commen\u00e7aient \u00e0 permettre l&rsquo;estimation s\u00e9quentielle de multiples \u00e9quations de r\u00e9gression, la communaut\u00e9 scientifique manquait d&rsquo;une m\u00e9trique standardis\u00e9e pour sanctionner l&rsquo;introduction superflue de r\u00e9gresseurs non informatifs. Mallows a d\u00e9velopp\u00e9 cette approche non pas comme un test d&rsquo;hypoth\u00e8se formel reposant sur des valeurs critiques rigides, mais comme un dispositif graphique et diagnostique destin\u00e9 \u00e0 guider le praticien dans l&rsquo;exploration raisonn\u00e9e des sous-espaces vectoriels de pr\u00e9diction.<\/p>\n<p>Au fil des d\u00e9cennies, cette m\u00e9thodologie s&rsquo;est profond\u00e9ment enracin\u00e9e dans l&rsquo;\u00e9conom\u00e9trie appliqu\u00e9e et la psychom\u00e9trie quantitative. Dans les sciences du comportement, l&rsquo;investigation des construits latents implique fr\u00e9quemment l&rsquo;administration de batteries de questionnaires g\u00e9n\u00e9rant des dizaines d&rsquo;indicateurs psychologiques potentiellement redondants. Face \u00e0 des matrices de donn\u00e9es o\u00f9 le nombre de pr\u00e9dicteurs cro\u00eet rapidement par rapport au nombre de participants, le risque d&rsquo;incorporer du bruit stochastique dans les \u00e9quations structurelles devient critique. Le crit\u00e8re de Mallows a offert aux chercheurs un cadre formel pour \u00e9laguer les inventaires de personnalit\u00e9, les \u00e9chelles cliniques et les indices cognitifs sans sacrifier la pr\u00e9cision globale de mesure.<\/p>\n<p>Dans le paysage contemporain de la science des donn\u00e9es, marqu\u00e9 par l&rsquo;h\u00e9g\u00e9monie apparente des algorithmes d&rsquo;apprentissage automatique de type bo\u00eete noire, la pertinence du Cp de Mallows demeure intacte. Alors que les macro-donn\u00e9es psychologiques issues du pistage num\u00e9rique et des enqu\u00eates en ligne massives multiplient les covariables disponibles, la n\u00e9cessit\u00e9 de pr\u00e9server des mod\u00e8les explicatifs hautement interpr\u00e9tables r\u00e9habilite les proc\u00e9dures de s\u00e9lection parcimonieuse. Le Cp de Mallows s&rsquo;affirme ainsi comme un pont conceptuel fondamental entre la rigueur de l&rsquo;inf\u00e9rence statistique classique et les imp\u00e9ratifs modernes de validation pr\u00e9dictive.<\/p>\n<h3>1.2 D\u00e9finition conceptuelle et utilit\u00e9 en s\u00e9lection de mod\u00e8les<\/h3>\n<p>Sur le plan conceptuel, le crit\u00e8re Cp de Mallows mesure l&rsquo;erreur quadratique totale moyenne standardis\u00e9e des pr\u00e9dictions g\u00e9n\u00e9r\u00e9es par un mod\u00e8le restreint contenant un sous-ensemble de variables, comparativement \u00e0 un mod\u00e8le satur\u00e9 ou complet englobant l&rsquo;int\u00e9gralit\u00e9 des pr\u00e9dicteurs disponibles. L&rsquo;intuition fondamentale qui sous-tend ce crit\u00e8re repose sur le constat suivant : l&rsquo;omission d&rsquo;un pr\u00e9dicteur v\u00e9ritablement associ\u00e9 \u00e0 la variable r\u00e9ponse introduit un biais syst\u00e9matique dans l&rsquo;estimation des coefficients restants, tandis que l&rsquo;inclusion d&rsquo;un pr\u00e9dicteur superflu augmente inutilement la variance d&rsquo;\u00e9chantillonnage de toutes les pr\u00e9dictions. Le Cp op\u00e8re ainsi un arbitrage objectif et transparent entre ces deux forces antagonistes.<\/p>\n<p>Dans une d\u00e9marche de mod\u00e9lisation par r\u00e9gression, l&rsquo;analyste se trouve confront\u00e9 \u00e0 un ensemble de mod\u00e8les candidats g\u00e9n\u00e9r\u00e9s soit par des consid\u00e9rations th\u00e9oriques a priori, soit par des algorithmes combinatoires exhaustifs. L&rsquo;utilit\u00e9 du Cp r\u00e9side dans sa capacit\u00e9 \u00e0 projeter l&rsquo;ensemble de ces sp\u00e9cifications comp\u00e9titives sur une \u00e9chelle unidimensionnelle standardis\u00e9e, directement comparable au nombre de param\u00e8tres estim\u00e9s. Cette m\u00e9trique ne requiert pas que les mod\u00e8les compar\u00e9s soient strictement embo\u00eet\u00e9s au sens des tests de rapport de vraisemblance, pourvu qu&rsquo;ils d\u00e9rivent tous du m\u00eame espace de variables d\u00e9fini par le mod\u00e8le de r\u00e9f\u00e9rence complet.<\/p>\n<p>Il convient n\u00e9anmoins de souligner la distinction \u00e9pist\u00e9mologique cruciale entre la pr\u00e9diction empirique et l&rsquo;inf\u00e9rence causale lors de l&rsquo;usage du Cp de Mallows. Ce crit\u00e8re est explicitement calibr\u00e9 pour minimiser l&rsquo;erreur quadratique moyenne de pr\u00e9diction (<em>Mean Squared Error of Prediction<\/em>). Par cons\u00e9quent, un mod\u00e8le s\u00e9lectionn\u00e9 par le Cp optimise la pr\u00e9cision avec laquelle de futures observations issues de la m\u00eame population pourront \u00eatre estim\u00e9es. Cela n&rsquo;implique nullement que les coefficients de r\u00e9gression retenus poss\u00e8dent une interpr\u00e9tabilit\u00e9 causale non biais\u00e9e, particuli\u00e8rement si des variables de confusion majeures ont \u00e9t\u00e9 omises d\u00e8s la formulation du mod\u00e8le complet de r\u00e9f\u00e9rence.<\/p>\n<h3>1.3 La r\u00e8gle d\u00e9cisionnelle fondamentale du Cp de Mallows<\/h3>\n<p>La mise en \u0153uvre pratique du crit\u00e8re de Mallows repose sur une r\u00e8gle de d\u00e9cision \u00e9l\u00e9gante et imm\u00e9diatement visualisable. Soit un mod\u00e8le candidat comportant <em>p<\/em> param\u00e8tres (incluant usuellement la constante d&rsquo;ajustement ou <em>intercept<\/em>). L&rsquo;esp\u00e9rance math\u00e9matique de la statistique Cp d\u00e9montre que pour un mod\u00e8le correctement sp\u00e9cifi\u00e9, exempt de biais d&rsquo;omission de variables pertinentes, la valeur attendue de Cp est approximativement \u00e9gale \u00e0 <em>p<\/em>. Cette \u00e9galit\u00e9 fondamentale, formalis\u00e9e par l&rsquo;expression math\u00e9matique E[Cp | biais nul] = p, constitue la pierre angulaire de l&rsquo;interpr\u00e9tation diagnostique.<\/p>\n<p>D\u00e8s lors, l&rsquo;analyste \u00e9value les mod\u00e8les candidats en confrontant la valeur num\u00e9rique de leur Cp \u00e0 leur cardinalit\u00e9 param\u00e9trique <em>p<\/em>. Deux trajectoires de d\u00e9viation apparaissent fr\u00e9quemment dans les applications concr\u00e8tes :<\/p>\n<ul>\n<li><strong>Cas o\u00f9 Cp est nettement sup\u00e9rieur \u00e0 p (Cp &gt;&gt; p) :<\/strong> Cette situation traduit l&rsquo;existence d&rsquo;un biais substantiel d\u00fb \u00e0 la sous-sp\u00e9cification du mod\u00e8le. Le sous-ensemble de variables s\u00e9lectionn\u00e9 \u00e9choue \u00e0 capturer une part essentielle de la variance structur\u00e9e du ph\u00e9nom\u00e8ne \u00e9tudi\u00e9, conduisant \u00e0 des pr\u00e9dictions syst\u00e9matiquement d\u00e9cal\u00e9es par rapport \u00e0 la r\u00e9alit\u00e9.<\/li>\n<li><strong>Cas o\u00f9 Cp est proche ou inf\u00e9rieur \u00e0 p (Cp &le; p) :<\/strong> Cette configuration indique que le biais r\u00e9siduel est n\u00e9gligeable ou statistiquement indiscernable du bruit d&rsquo;\u00e9chantillonnage. Lorsque le Cp tombe l\u00e9g\u00e8rement en dessous de <em>p<\/em>, ce ph\u00e9nom\u00e8ne r\u00e9sulte le plus souvent de fluctuations stochastiques d&rsquo;\u00e9chantillonnage favorables dans l&rsquo;estimation de la variance de l&rsquo;erreur.<\/li>\n<\/ul>\n<p>La r\u00e8gle heuristique universellement reconnue en mod\u00e9lisation consiste donc \u00e0 identifier, parmi l&rsquo;ensemble des sp\u00e9cifications candidates, les mod\u00e8les pour lesquels la valeur de Cp est \u00e0 la fois faible dans l&rsquo;absolu et extr\u00eamement proche de <em>p<\/em>. En pr\u00e9sence de plusieurs mod\u00e8les satisfaisant \u00e0 cette proximit\u00e9 (Cp &asymp; p), le principe de parcimonie scientifique enjoint de retenir le mod\u00e8le pr\u00e9sentant la valeur de <em>p<\/em> la plus faible, r\u00e9duisant ainsi la complexit\u00e9 cognitive et op\u00e9rationnelle du syst\u00e8me explicatif.<\/p>\n<h2>2. Fondements math\u00e9matiques et d\u00e9rivation du Cp de Mallows<\/h2>\n<h3>2.1 Formulation alg\u00e9brique standard du crit\u00e8re<\/h3>\n<p>Pour appr\u00e9hender la m\u00e9canique interne du crit\u00e8re Cp de Mallows, il est n\u00e9cessaire de formaliser le cadre de la r\u00e9gression lin\u00e9aire sous forme matricielle. Consid\u00e9rons un \u00e9chantillon compos\u00e9 de <em>n<\/em> observations ind\u00e9pendantes. Le mod\u00e8le complet, suppos\u00e9 contenir l&rsquo;ensemble des <em>K<\/em> r\u00e9gresseurs potentiels plus la constante (soit <em>P = K + 1<\/em> param\u00e8tres), s&rsquo;\u00e9crit sous la forme vectorielle standard Y = X&beta; + &epsilon;, o\u00f9 Y est le vecteur de dimension (n &times; 1) des r\u00e9ponses observ\u00e9es, X est la matrice de dessein de dimension (n &times; P), &beta; est le vecteur des param\u00e8tres de r\u00e9gression et &epsilon; est le vecteur des perturbations al\u00e9atoires suivant une loi normale multivari\u00e9e d&rsquo;esp\u00e9rance nulle et de matrice de covariance &sigma;<sup>2<\/sup>I<sub>n<\/sub>.<\/p>\n<p>Soit d\u00e9sormais un sous-mod\u00e8le restreint n&rsquo;incorporant qu&rsquo;un sous-ensemble de pr\u00e9dicteurs correspondant \u00e0 <em>p<\/em> param\u00e8tres (avec <em>p &lt; P<\/em>). La matrice de projection associ\u00e9e \u00e0 ce sous-mod\u00e8le projette le vecteur Y sur le sous-espace vectoriel engendr\u00e9 par les colonnes s\u00e9lectionn\u00e9es de X, produisant le vecteur des valeurs ajust\u00e9es Y_chapeau<sub>p<\/sub>. La somme des carr\u00e9s des r\u00e9sidus associ\u00e9e \u00e0 ce mod\u00e8le r\u00e9duit, usuellement d\u00e9sign\u00e9e par l&rsquo;acronyme RSS<sub>p<\/sub> (pour <em>Residual Sum of Squares<\/em>), s&rsquo;exprime formellement par :<\/p>\n<p>RSS<sub>p<\/sub> = &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &#8211; y_chapeau<sub>i,p<\/sub>)<sup>2<\/sup> = Y<sup>T<\/sup> (I<sub>n<\/sub> &#8211; H<sub>p<\/sub>) Y<\/p>\n<p>o\u00f9 H<sub>p<\/sub> = X<sub>p<\/sub>(X<sub>p<\/sub><sup>T<\/sup>X<sub>p<\/sub>)<sup>-1<\/sup>X<sub>p<\/sub><sup>T<\/sup> repr\u00e9sente la matrice de projection orthogonale (ou matrice chapeau) de rang <em>p<\/em>. Le crit\u00e8re Cp de Mallows normalise cette quantit\u00e9 r\u00e9siduelle en exploitant un estimateur sans biais de la variance de l&rsquo;erreur issu exclusivement du mod\u00e8le complet, not\u00e9 s<sup>2<\/sup> (ou &sigma;_chapeau<sup>2<\/sup>). L&rsquo;\u00e9quation alg\u00e9brique d\u00e9finitive du crit\u00e8re s&rsquo;\u00e9nonce comme suit :<\/p>\n<p>Cp = (RSS<sub>p<\/sub> \/ s<sup>2<\/sup>) &#8211; n + 2p<\/p>\n<p>Dans cette expression fondamentale, le terme (RSS<sub>p<\/sub> \/ s<sup>2<\/sup>) refl\u00e8te l&rsquo;ad\u00e9quation brute du mod\u00e8le r\u00e9duit aux donn\u00e9es, pond\u00e9r\u00e9e par l&rsquo;incertitude intrins\u00e8que du syst\u00e8me. Le terme constant (- n) op\u00e8re un recentrage dimensionnel, tandis que la composante (+ 2p) agit comme un m\u00e9canisme de p\u00e9nalisation lin\u00e9aire proportionnel au nombre de coefficients allou\u00e9s dans le sous-mod\u00e8le.<\/p>\n<h3>2.2 Hypoth\u00e8ses sous-jacentes du mod\u00e8le complet<\/h3>\n<p>La validit\u00e9 math\u00e9matique et l&rsquo;efficacit\u00e9 op\u00e9rationnelle du crit\u00e8re Cp reposent sur une condition pr\u00e9alable imp\u00e9rative : le mod\u00e8le complet \u00e0 partir duquel s<sup>2<\/sup> est d\u00e9riv\u00e9 doit \u00eatre rigoureusement exempt de biais de sp\u00e9cification. En termes statistiques, cela signifie que la matrice de r\u00e9gresseurs du mod\u00e8le satur\u00e9 doit inclure toutes les variables explicatives structurales pertinentes pour d\u00e9crire l&rsquo;esp\u00e9rance de Y, garantissant ainsi que l&rsquo;esp\u00e9rance math\u00e9matique de la somme des carr\u00e9s r\u00e9siduelle du mod\u00e8le complet divis\u00e9e par ses degr\u00e9s de libert\u00e9 (n &#8211; P) constitue un estimateur sans biais et convergent de la variance th\u00e9orique &sigma;<sup>2<\/sup> :<\/p>\n<p>E[s<sup>2<\/sup>] = E[RSS<sub>P<\/sub> \/ (n &#8211; P)] = &sigma;<sup>2<\/sup><\/p>\n<p>Si cette hypoth\u00e8se de base est viol\u00e9e \u2014 par exemple si le chercheur a omis des variables d\u00e9terminantes d\u00e8s la phase de formulation th\u00e9orique globale \u2014, l&rsquo;estimateur s<sup>2<\/sup> absorbera une fraction de la variance syst\u00e9matique non expliqu\u00e9e, conduisant \u00e0 une sur\u00e9valuation grossi\u00e8re de la variance de l&rsquo;erreur r\u00e9siduelle. Dans une telle conjoncture, le d\u00e9nominateur de la fraction (RSS<sub>p<\/sub> \/ s<sup>2<\/sup>) devient artificiellement \u00e9lev\u00e9, ce qui fausse l&rsquo;ensemble des calculs de Cp, comprime les valeurs obtenues vers le bas et entra\u00eene une s\u00e9lection erron\u00e9e de sous-mod\u00e8les sous-dimensionn\u00e9s.<\/p>\n<p>Par ailleurs, les hypoth\u00e8ses canoniques de la r\u00e9gression lin\u00e9aire doivent \u00eatre satisfaites : homosc\u00e9dasticit\u00e9 parfaite (variance constante des r\u00e9sidus \u00e0 travers toutes les strates de pr\u00e9diction), ind\u00e9pendance stochastique des perturbations r\u00e9siduelles et normalit\u00e9 gaussienne des erreurs. De surcro\u00eet, la pr\u00e9sence de valeurs aberrantes extr\u00eames (<em>outliers<\/em>) ou de points \u00e0 fort effet de levier (<em>high leverage points<\/em>) dans la matrice X peut alt\u00e9rer drastiquement l&rsquo;estimation de s<sup>2<\/sup> et de RSS<sub>p<\/sub>, d\u00e9stabilisant compl\u00e8tement la hi\u00e9rarchie relative des valeurs de Cp calcul\u00e9es sur les diff\u00e9rents sous-ensembles.<\/p>\n<h3>2.3 Composantes de l&rsquo;esp\u00e9rance math\u00e9matique de Cp<\/h3>\n<p>L&rsquo;\u00e9l\u00e9gance th\u00e9orique de la m\u00e9trique de Mallows se r\u00e9v\u00e8le pleinement \u00e0 travers la d\u00e9composition analytique de son esp\u00e9rance math\u00e9matique. Soit l&rsquo;erreur quadratique totale moyenne de pr\u00e9diction normalis\u00e9e, not\u00e9e &Gamma;<sub>p<\/sub> (Gamma de p), d\u00e9finie par la sommation des esp\u00e9rances des \u00e9carts au carr\u00e9 entre la v\u00e9ritable fonction d&rsquo;esp\u00e9rance conditionnelle &mu; = E[Y|X] et les pr\u00e9dictions vectorielles Y_chapeau<sub>p<\/sub> issues du sous-mod\u00e8le :<\/p>\n<p>&Gamma;<sub>p<\/sub> = (1 \/ &sigma;<sup>2<\/sup>) &sum;<sub>i=1<\/sub><sup>n<\/sup> E[(y_chapeau<sub>i,p<\/sub> &#8211; &mu;<sub>i<\/sub>)<sup>2<\/sup>]<\/p>\n<p>En d\u00e9composant cet \u00e9cart en termes de biais quadratique d&rsquo;omission et de variance d&rsquo;estimation par les propri\u00e9t\u00e9s des projecteurs matriciels, on d\u00e9montre rigoureusement l&rsquo;\u00e9galit\u00e9 suivante :<\/p>\n<p>&Gamma;<sub>p<\/sub> = (Biais<sub>p<\/sub> \/ &sigma;<sup>2<\/sup>) + p<\/p>\n<p>o\u00f9 Biais<sub>p<\/sub> correspond \u00e0 la quantit\u00e9 &sum;<sub>i=1<\/sub><sup>n<\/sup> (E[y_chapeau<sub>i,p<\/sub>] &#8211; &mu;<sub>i<\/sub>)<sup>2<\/sup>. Lorsque le mod\u00e8le r\u00e9duit ne souffre d&rsquo;aucun biais de sp\u00e9cification \u2014 c&rsquo;est-\u00e0-dire que tous les r\u00e9gresseurs omis dans le passage de <em>P<\/em> \u00e0 <em>p<\/em> poss\u00e9daient des coefficients th\u00e9oriques rigoureusement \u00e9gaux \u00e0 z\u00e9ro \u2014, la composante Biais<sub>p<\/sub> s&rsquo;annule identiquement. Dans ces conditions de sp\u00e9cification ad\u00e9quate, &Gamma;<sub>p<\/sub> se r\u00e9duit strictement \u00e0 sa composante de variance, soit exactement la valeur <em>p<\/em>.<\/p>\n<p>La statistique Cp a \u00e9t\u00e9 pr\u00e9cis\u00e9ment configur\u00e9e par Colin Mallows pour constituer un estimateur non biais\u00e9 de cette quantit\u00e9 th\u00e9orique &Gamma;<sub>p<\/sub>. En prenant l&rsquo;esp\u00e9rance math\u00e9matique conditionnelle de Cp sous l&rsquo;hypoth\u00e8se que s<sup>2<\/sup> est un estimateur coh\u00e9rent de &sigma;<sup>2<\/sup> ind\u00e9pendant des r\u00e9sidus du sous-mod\u00e8le, on obtient :<\/p>\n<p>E[Cp] = E[(RSS<sub>p<\/sub> \/ &sigma;<sup>2<\/sup>) &#8211; n + 2p] = ((n &#8211; p + Biais<sub>p<\/sub> \/ &sigma;<sup>2<\/sup>)) &#8211; n + 2p = p + (Biais<sub>p<\/sub> \/ &sigma;<sup>2<\/sup>)<\/p>\n<p>Cette d\u00e9rivation formelle d\u00e9montre sans \u00e9quivoque que si le sous-mod\u00e8le souffre d&rsquo;un biais d&rsquo;omission, E[Cp] exc\u00e9dera syst\u00e9matiquement <em>p<\/em> d&rsquo;une quantit\u00e9 strictement proportionnelle \u00e0 l&rsquo;ampleur du biais quadratique standardis\u00e9. Inversement, si le mod\u00e8le n&rsquo;introduit aucun biais d&rsquo;omission, l&rsquo;esp\u00e9rance s&rsquo;\u00e9tablit exactement \u00e0 <em>p<\/em>. Toute d\u00e9viation empirique observable o\u00f9 Cp est inf\u00e9rieur \u00e0 <em>p<\/em> s&rsquo;explique par la nature al\u00e9atoire de l&rsquo;\u00e9chantillonnage, la variable al\u00e9atoire RSS<sub>p<\/sub> suivant une distribution du chi-deux non centr\u00e9e.<\/p>\n<h2>3. Le compromis biais-variance et l&rsquo;interpr\u00e9tation des param\u00e8tres<\/h2>\n<h3>3.1 Dilemme classique entre complexit\u00e9 et g\u00e9n\u00e9ralisation<\/h3>\n<p>Le d\u00e9veloppement de mod\u00e8les statistiques explicatifs ou pr\u00e9dictifs en psychologie exp\u00e9rimentale et dans les sciences sociales s&rsquo;inscrit au c\u0153ur du dilemme classique entre biais et variance. Un mod\u00e8le sur-param\u00e9tr\u00e9 pr\u00e9sente une tr\u00e8s faible erreur d&rsquo;ajustement sur l&rsquo;\u00e9chantillon initial parce qu&rsquo;il poss\u00e8de une flexibilit\u00e9 excessive qui lui permet d&rsquo;\u00e9pouser non seulement les tendances structurelles profondes des variables, mais \u00e9galement les fluctuations al\u00e9atoires idiosyncrasiques propres au recueil de donn\u00e9es. Cette perm\u00e9abilit\u00e9 au bruit statistique g\u00e9n\u00e8re une variance d&rsquo;\u00e9chantillonnage \u00e9lev\u00e9e : si l&rsquo;on r\u00e9p\u00e9tait l&rsquo;\u00e9tude sur une autre cohorte de sujets, les coefficients estim\u00e9s subiraient d&rsquo;importantes oscillations.<\/p>\n<p>\u00c0 l&rsquo;oppos\u00e9 du spectre, un mod\u00e8le \u0447\u0440\u0435\u0437\u043c\u0435\u0440ement parcimonieux compresse l&rsquo;espace explicatif en \u00e9liminant des covariables pertinentes. Si cette simplification drastique a pour effet de stabiliser les coefficients restants \u2014 conf\u00e9rant ainsi une faible variance d&rsquo;\u00e9chantillonnage aux estimateurs \u2014, elle injecte un biais syst\u00e9matique non nul dans toutes les inf\u00e9rences subs\u00e9quentes. Les pr\u00e9dictions \u00e9mises par un tel mod\u00e8le seront invariablement d\u00e9cal\u00e9es par rapport aux v\u00e9ritables processus g\u00e9n\u00e9rateurs de donn\u00e9es.<\/p>\n<p>Le crit\u00e8re Cp de Mallows op\u00e9rationnalise pr\u00e9cis\u00e9ment ce compromis en formalisant le point d&rsquo;inflexion o\u00f9 le co\u00fbt marginal li\u00e9 \u00e0 l&rsquo;augmentation de la variance (p\u00e9nalis\u00e9 par le terme + 2p) compense exactement le gain marginal obtenu en termes de r\u00e9duction de la somme des carr\u00e9s des r\u00e9sidus (RSS<sub>p<\/sub>). Il prot\u00e8ge l&rsquo;investigateur contre les mirages de l&rsquo;optimisation intra-\u00e9chantillon en s\u00e9lectionnant les architectures de r\u00e9gression qui maximisent les capacit\u00e9s de r\u00e9plication externe et de g\u00e9n\u00e9ralisation hors-\u00e9chantillon.<\/p>\n<h3>3.2 Convention de notation : p versus p + 1<\/h3>\n<p>Une source r\u00e9currente de confusion et d&rsquo;erreurs m\u00e9thodologiques dans la litt\u00e9rature empirique et les manuels d&rsquo;\u00e9conom\u00e9trie r\u00e9side dans la d\u00e9nomination formelle du nombre de param\u00e8tres, sp\u00e9cifiquement en ce qui concerne le statut de la constante d&rsquo;ajustement (le terme d&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine ou <em>intercept<\/em>). Selon les disciplines et les traditions logicielles, le symbole <em>p<\/em> peut d\u00e9signer soit le nombre total de pr\u00e9dicteurs stochastiques (hors constante), soit le nombre absolu de param\u00e8tres libres de l&rsquo;\u00e9quation de r\u00e9gression.<\/p>\n<p>Dans la formulation th\u00e9orique originale de Mallows ainsi que dans la majorit\u00e9 des progiciels statistiques sous R, <strong>p repr\u00e9sente l&rsquo;ensemble des param\u00e8tres estim\u00e9s<\/strong>, ce qui inclut syst\u00e9matiquement la constante si celle-ci figure dans l&rsquo;\u00e9quation. Ainsi, pour une r\u00e9gression lin\u00e9aire standard mod\u00e9lisant une variable r\u00e9ponse continue \u00e0 partir de 4 variables ind\u00e9pendantes et d&rsquo;une constante, le nombre de param\u00e8tres s&rsquo;\u00e9l\u00e8ve rigoureusement \u00e0 :<\/p>\n<p>p = 4 r\u00e9gresseurs + 1 intercept = 5 param\u00e8tres<\/p>\n<p>Cette pr\u00e9cision est fondamentale pour l&rsquo;application des r\u00e8gles de d\u00e9cision. Si un chercheur confronte une valeur de Cp calcul\u00e9e \u00e0 4 au lieu de 5 en omettant de comptabiliser l&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine, il conclura \u00e0 tort \u00e0 l&rsquo;existence d&rsquo;une inad\u00e9quation ou d&rsquo;un biais syst\u00e9matique. Les packages R modernes, tels que <em>olsrr<\/em> et <em>leaps<\/em>, adoptent heureusement cette convention int\u00e9gr\u00e9e : la ligne d&rsquo;\u00e9quivalence sur laquelle le Cp th\u00e9orique id\u00e9al s&rsquo;aligne correspond toujours au nombre total de colonnes lin\u00e9airement ind\u00e9pendantes engag\u00e9es dans la sous-matrice de projection X<sub>p<\/sub>.<\/p>\n<h3>3.3 Typologie des sc\u00e9narios de s\u00e9lection<\/h3>\n<p>L&rsquo;analyse des valeurs de Cp issues d&rsquo;une proc\u00e9dure d&rsquo;\u00e9valuation de mod\u00e8les candidats permet de cat\u00e9goriser les sous-ensembles de variables en plusieurs sc\u00e9narios distincts, dictant la conduite diagnostique du statisticien :<\/p>\n<ul>\n<li><strong>Le sc\u00e9nario d&rsquo;ad\u00e9quation optimale (Cp &asymp; p avec p minimal) :<\/strong> Le mod\u00e8le r\u00e9duit pr\u00e9sente une statistique Cp extr\u00eamement proche de son cardinal de param\u00e8tres (par exemple un Cp de 4.12 pour p = 4). Cela signifie que le biais d&rsquo;omission est virtuellement nul tout en ayant \u00e9pur\u00e9 l&rsquo;\u00e9quation de r\u00e9gression des dimensions parasitaires. Ce mod\u00e8le constitue la cible privil\u00e9gi\u00e9e du chercheur.<\/li>\n<li><strong>Le sc\u00e9nario de sous-param\u00e9trisation s\u00e9v\u00e8re (Cp &gt;&gt; p) :<\/strong> Le mod\u00e8le test\u00e9 affiche une valeur Cp disproportionnellement \u00e9lev\u00e9e par rapport \u00e0 sa dimension (par exemple un Cp de 38.6 pour p = 3). Un tel foss\u00e9 num\u00e9rique atteste qu&rsquo;une ou plusieurs variables d\u00e9terminantes ont \u00e9t\u00e9 exclues de la sp\u00e9cification. Les estimations sont grev\u00e9es d&rsquo;un biais structurel substantiel qui rend le mod\u00e8le inapte \u00e0 la pr\u00e9diction rigoureuse.<\/li>\n<li><strong>Le sc\u00e9nario de sur-param\u00e9trisation redondante (Cp &asymp; p avec p \u00e9lev\u00e9) :<\/strong> Lorsque l&rsquo;on examine des mod\u00e8les dont la taille approche celle du mod\u00e8le complet <em>P<\/em>, les valeurs de Cp convergent n\u00e9cessairement vers <em>P<\/em> (puisque pour le mod\u00e8le complet, Cp est par d\u00e9finition \u00e9gal \u00e0 P). Cependant, conserver ces variables surnum\u00e9raires n&rsquo;apporte aucune r\u00e9duction du biais et d\u00e9grade la variance des estimations.<\/li>\n<li><strong>Le sc\u00e9nario sous-unitaire ou d&rsquo;\u00e9chantillonnage favorable (Cp &lt; p) :<\/strong> Dans certains cas empiriques, l&rsquo;analyste observe un Cp inf\u00e9rieur \u00e0 la dimension <em>p<\/em> (par exemple un Cp de 2.8 pour p = 4). Loin d&rsquo;indiquer une anomalie, ce r\u00e9sultat traduit le fait que le sous-mod\u00e8le b\u00e9n\u00e9ficie d&rsquo;une variance d&rsquo;\u00e9chantillonnage r\u00e9siduelle l\u00e9g\u00e8rement inf\u00e9rieure \u00e0 la pr\u00e9diction asymptotique moyenne sous distribution gaussienne.<\/li>\n<\/ul>\n<p>Face \u00e0 plusieurs sp\u00e9cifications pr\u00e9sentant simultan\u00e9ment un \u00e9cart |Cp &#8211; p| minimal, le chercheur devra arbitrer non seulement sur la base du crit\u00e8re num\u00e9rique, mais \u00e9galement \u00e0 la lumi\u00e8re de la coh\u00e9rence th\u00e9orique, de la faisabilit\u00e9 logistique de collecte des variables et de la parcimonie conceptuelle.<\/p>\n<h2>4. Environnement informatique R et pr\u00e9paration des donn\u00e9es psychom\u00e9triques<\/h2>\n<h3>4.1 Installation et configuration des biblioth\u00e8ques n\u00e9cessaires<\/h3>\n<p>La mise en \u0153uvre computationnelle du crit\u00e8re Cp de Mallows sous l&rsquo;environnement <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/\">CRAN (Comprehensive R Archive Network)<\/a> b\u00e9n\u00e9ficie d&rsquo;un riche \u00e9cosyst\u00e8me logiciel. Afin de garantir un environnement de travail reproductible et fluide, plusieurs biblioth\u00e8ques sp\u00e9cialis\u00e9es doivent \u00eatre d\u00e9ploy\u00e9es. Le package <em>olsrr<\/em> (Ordinary Least Squares Regression Reports) s&rsquo;impose comme la r\u00e9f\u00e9rence contemporaine pour l&rsquo;analyse p\u00e9dagogique et visuelle des r\u00e9gressions lin\u00e9aires ordinaires. Parall\u00e8lement, le package <em>leaps<\/em>, con\u00e7u originellement par Thomas Lumley, offre des routines algorithmiques compil\u00e9es en Fortran permettant l&rsquo;exploration exhaustive ultra-rapide des espaces combinatoires de pr\u00e9dicteurs.<\/p>\n<p>Pour orchestrer la manipulation des structures de donn\u00e9es, la production de sorties graphiques de haute qualit\u00e9 \u00e9ditoriale et l&rsquo;\u00e9valuation psychom\u00e9trique pr\u00e9alable des \u00e9chelles de mesure, nous mobiliserons conjointement l&rsquo;\u00e9cosyst\u00e8me <em>tidyverse<\/em> (notamment <em>dplyr<\/em> et <em>ggplot2<\/em>), le package <em>car<\/em> (Companion to Applied Regression) pour le calcul des diagnostics de colin\u00e9arit\u00e9, ainsi que le package <em>psych<\/em> d\u00e9di\u00e9 \u00e0 l&rsquo;analyse descriptive des construits psychologiques. L&rsquo;instruction d&rsquo;installation et de chargement coordonn\u00e9 de ces extensions s&rsquo;articule via le script R suivant :<\/p>\n<p><code>install.packages(c(\"olsrr\", \"leaps\", \"car\", \"tidyverse\", \"psych\"))<\/code><\/p>\n<p><code>library(olsrr)<\/code><br \/>\n<code>library(leaps)<\/code><br \/>\n<code>library(car)<\/code><br \/>\n<code>library(tidyverse)<\/code><br \/>\n<code>library(psych)<\/code><\/p>\n<p>Il est recommand\u00e9 de s&rsquo;assurer que la version active de R soit \u00e9gale ou sup\u00e9rieure \u00e0 la version 4.1.0, garantissant l&rsquo;acc\u00e8s aux optimisations natives de l&rsquo;op\u00e9rateur de tuyau (<em>pipe<\/em>) et une gestion m\u00e9moire optimis\u00e9e lors de la manipulation de matrices de r\u00e9gression de taille interm\u00e9diaire \u00e0 grande.<\/p>\n<h3>4.2 Structure et pr\u00e9traitement du jeu de donn\u00e9es<\/h3>\n<p>En psychom\u00e9trie et en sciences sociales appliqu\u00e9es, les donn\u00e9es empiriques proviennent typiquement d&rsquo;\u00e9chelles d&rsquo;auto-\u00e9valuation de type Likert agr\u00e9g\u00e9es en scores composites continus, combin\u00e9es \u00e0 des mesures biom\u00e9triques, comportementales ou sociod\u00e9mographiques. Consid\u00e9rons, \u00e0 titre d&rsquo;illustration op\u00e9rationnelle tout au long de ce guide, une cohorte de recherche en psychologie de la sant\u00e9 \u00e9valuant le niveau d&rsquo;\u00e9puisement professionnel (<em>burnout<\/em>) aupr\u00e8s d&rsquo;un \u00e9chantillon d&rsquo;individus suivis en milieu organisationnel.<\/p>\n<p>La variable d&rsquo;int\u00e9r\u00eat continue \u00e0 pr\u00e9dire est le score global d&rsquo;\u00e9puisement \u00e9motionnel (mesur\u00e9 par l&rsquo;inventaire MBI). Les covariables candidates \u00e0 l&rsquo;inclusion dans le syst\u00e8me pr\u00e9dictif comprennent : le stress per\u00e7u (PSS-10), le n\u00e9vrosisme (facette du Big Five), la charge horaire hebdomadaire moyenne, l&rsquo;indice de soutien social per\u00e7u au travail, l&rsquo;auto-efficacit\u00e9 g\u00e9n\u00e9rale, la qualit\u00e9 subjective du sommeil (indice PSQI), et l&rsquo;anciennet\u00e9 professionnelle en ann\u00e9es. Avant toute entreprise de mod\u00e9lisation math\u00e9matique, une inspection rigoureuse de l&rsquo;int\u00e9grit\u00e9 des vecteurs de donn\u00e9es est imp\u00e9rative.<\/p>\n<p>Le pr\u00e9traitement doit adresser de mani\u00e8re syst\u00e9matique les valeurs manquantes (<em>missing values<\/em>). Dans le contexte sp\u00e9cifique de la s\u00e9lection de variables via le crit\u00e8re Cp de Mallows, une r\u00e8gle m\u00e9thodologique fondamentale s&rsquo;impose : <strong>tous les sous-mod\u00e8les candidats doivent imp\u00e9rativement \u00eatre estim\u00e9s sur un \u00e9chantillon d&rsquo;observations rigoureusement identique<\/strong>. Si des lignes sont \u00e9cart\u00e9es dynamiquement en raison de valeurs manquantes dispers\u00e9es \u00e0 travers les diff\u00e9rents r\u00e9gresseurs, la variance de r\u00e9f\u00e9rence s<sup>2<\/sup> et les sommes des carr\u00e9s des r\u00e9sidus RSS<sub>p<\/sub> deviendront math\u00e9matiquement incommensurables, invalidant totalement la proc\u00e9dure de comparaison. Une imputation prudente ou l&rsquo;exclusion cas-par-cas globale (<em>listwise deletion<\/em>) pr\u00e9alable est donc requise.<\/p>\n<h3>4.3 D\u00e9finition formelle de l&rsquo;\u00e9chantillon d&rsquo;analyse<\/h3>\n<p>Pour assurer la solidit\u00e9 inf\u00e9rentielle de l&rsquo;exercice, nous g\u00e9n\u00e9rons synth\u00e9tiquement un jeu de donn\u00e9es psychom\u00e9triques r\u00e9aliste comportant 250 participants (n = 250), dimensionn\u00e9 de fa\u00e7on ad\u00e9quate pour respecter un ratio sujets\/variables d&rsquo;au moins 25 observations par pr\u00e9dicteur, minimisant ainsi les risques de sur-ajustement accidentel d\u00e8s le stade de la conception. Le script R ci-dessous configure la matrice de donn\u00e9es avec une structure de covariance empirique ma\u00eetris\u00e9e :<\/p>\n<p><code>set.seed(1973) # Hommage a l'annee de publication de Colin Mallows<\/code><br \/>\n<code>n &lt;- 250<\/code><br \/>\n<code>stress &lt;- rnorm(n, mean = 25, sd = 6)<\/code><br \/>\n<code>nevrosisme &lt;- rnorm(n, mean = 50, sd = 10) + 0.4 * stress<\/code><br \/>\n<code>charge_travail &lt;- rnorm(n, mean = 40, sd = 8)<\/code><br \/>\n<code>soutien_social &lt;- rnorm(n, mean = 60, sd = 12) - 0.3 * stress<\/code><br \/>\n<code>auto_efficacite &lt;- rnorm(n, mean = 30, sd = 5) - 0.2 * nevrosisme<\/code><br \/>\n<code>sommeil_degrade &lt;- rnorm(n, mean = 8, sd = 3) + 0.3 * stress<\/code><br \/>\n<code>anciennete &lt;- runif(n, min = 1, max = 25)<\/code><\/p>\n<p>Dans ce mod\u00e8le de simulation g\u00e9n\u00e9rative, concevons la variable r\u00e9ponse <em>epuisement<\/em> de telle sorte que seules certaines variables exercent un effet structural direct r\u00e9el, les autres n&rsquo;introduisant que des corr\u00e9lations collat\u00e9rales ou du bruit stochastique :<\/p>\n<p><code>epuisement &lt;- 15 + 0.8 * stress + 0.5 * nevrosisme - 0.4 * soutien_social + 1.2 * sommeil_degrade + rnorm(n, mean = 0, sd = 8)<\/code><br \/>\n<code>df_psycho &lt;- data.frame(epuisement, stress, nevrosisme, charge_travail, soutien_social, auto_efficacite, sommeil_degrade, anciennete)<\/code><\/p>\n<p>Cette configuration matricielle servira de substrat exp\u00e9rimental transparent : nous savons a priori que la charge de travail, l&rsquo;auto-efficacit\u00e9 et l&rsquo;anciennet\u00e9 sont structurellement redondantes dans l&rsquo;\u00e9quation de g\u00e9n\u00e9ration une fois contr\u00f4l\u00e9s le stress, le n\u00e9vrosisme, le soutien social et le sommeil. L&rsquo;objectif analytique consistera \u00e0 v\u00e9rifier si le crit\u00e8re Cp de Mallows sous R isole avec pr\u00e9cision ce sous-ensemble th\u00e9orique exact.<\/p>\n<h2>5. Ajustement du mod\u00e8le complet et des sous-mod\u00e8les via la fonction lm()<\/h2>\n<h3>5.1 Sp\u00e9cification rigoureuse du mod\u00e8le satur\u00e9<\/h3>\n<p>L&rsquo;estimation sous R d\u00e9bute obligatoirement par la sp\u00e9cification formelle du mod\u00e8le lin\u00e9aire complet au moyen de la primitive fondamentale <code>lm()<\/code>. Ce mod\u00e8le complet, parfois qualifi\u00e9 de mod\u00e8le satur\u00e9 ou de r\u00e9f\u00e9rence, doit incorporer simultan\u00e9ment l&rsquo;ensemble des r\u00e9gresseurs identifi\u00e9s comme scientifiquement plausibles. La syntaxe compacte de R permet de d\u00e9clarer la r\u00e9gression de la variable r\u00e9ponse sur l&rsquo;int\u00e9gralit\u00e9 des colonnes r\u00e9siduelles de la matrice de donn\u00e9es gr\u00e2ce \u00e0 la notation par point :<\/p>\n<p><code>modele_complet &lt;- lm(epuisement ~ ., data = df_psycho)<\/code><\/p>\n<p>L&rsquo;ajustement de cet objet de classe <code>lm<\/code> constitue la cl\u00e9 de vo\u00fbte de toute l&rsquo;\u00e9valuation subs\u00e9quente. C&rsquo;est \u00e0 partir de cet objet que l&rsquo;estimateur sans biais de la variance de l&rsquo;erreur r\u00e9siduelle sera extrait. Dans l&rsquo;environnement R, cette variance correspond au carr\u00e9 de la statistique d\u00e9sign\u00e9e sous le terme de <em>residual standard error<\/em> pr\u00e9sente dans la sortie textuelle standard. Alg\u00e9briquement, cette valeur est accessible directement via la structure interne de la liste R :<\/p>\n<p><code>sigma2_complet &lt;- summary(modele_complet)$sigma^2<\/code><\/p>\n<p>Il est capital d&rsquo;inspecter visuellement et num\u00e9riquement ce mod\u00e8le initial. Si le mod\u00e8le complet souffre de singularit\u00e9s math\u00e9matiques \u2014 par exemple si le rang de la matrice X est inf\u00e9rieur au nombre de colonnes en raison d&rsquo;une colin\u00e9arit\u00e9 parfaite \u2014, la fonction <code>lm()<\/code> produira des coefficients non d\u00e9finis (marqu\u00e9s par la valeur <code>NA<\/code>). Une telle situation invaliderait imm\u00e9diatement l&rsquo;estimation de &sigma;<sup>2<\/sup> et doit \u00eatre r\u00e9solue avant de proc\u00e9der au calcul du crit\u00e8re de Mallows.<\/p>\n<h3>5.2 Construction manuelle de mod\u00e8les concurrents embo\u00eet\u00e9s et non embo\u00eet\u00e9s<\/h3>\n<p>Pour illustrer la m\u00e9canique compar\u00e9e, proc\u00e9dons \u00e0 l&rsquo;ajustement explicite de plusieurs mod\u00e8les concurrents repr\u00e9sentant des hypoth\u00e8ses th\u00e9oriques alternatives quant aux d\u00e9terminants de l&rsquo;\u00e9puisement professionnel. Cette approche manuelle refl\u00e8te fid\u00e8lement la pratique de recherche confirmatoire o\u00f9 le statisticien teste un ensemble restreint de sp\u00e9cifications dict\u00e9es par la litt\u00e9rature scientifique :<\/p>\n<ul>\n<li><strong>Mod\u00e8le 1 (Mod\u00e8le purement biologique et contextuel) :<\/strong> \u00c9puisement pr\u00e9dit par la charge de travail, l&rsquo;anciennet\u00e9 et le sommeil d\u00e9grad\u00e9.<br \/>\n    <code>m1 &lt;- lm(epuisement ~ charge_travail + anciennete + sommeil_degrade, data = df_psycho)<\/code><\/li>\n<li><strong>Mod\u00e8le 2 (Mod\u00e8le cognitivo-\u00e9motionnel) :<\/strong> \u00c9puisement pr\u00e9dit par le stress per\u00e7u, le n\u00e9vrosisme et l&rsquo;auto-efficacit\u00e9.<br \/>\n    <code>m2 &lt;- lm(epuisement ~ stress + nevrosisme + auto_efficacite, data = df_psycho)<\/code><\/li>\n<li><strong>Mod\u00e8le 3 (Mod\u00e8le structurel congruent) :<\/strong> \u00c9puisement pr\u00e9dit par le stress per\u00e7u, le n\u00e9vrosisme, le soutien social et le sommeil d\u00e9grad\u00e9.<br \/>\n    <code>m3 &lt;- lm(epuisement ~ stress + nevrosisme + soutien_social + sommeil_degrade, data = df_psycho)<\/code><\/li>\n<li><strong>Mod\u00e8le 4 (Mod\u00e8le quasi-complet avec bruit) :<\/strong> \u00c9puisement pr\u00e9dit par toutes les variables \u00e0 l&rsquo;exclusion unique de l&rsquo;anciennet\u00e9.<br \/>\n    <code>m4 &lt;- lm(epuisement ~ stress + nevrosisme + charge_travail + soutien_social + auto_efficacite + sommeil_degrade, data = df_psycho)<\/code><\/li>\n<\/ul>\n<p>Chacun de ces objets encapsule une structure matricielle sp\u00e9cifique, caract\u00e9ris\u00e9e par une somme des carr\u00e9s r\u00e9siduelle distincte (accessible par la commande <code>deviance(m)<\/code>) et un nombre de degr\u00e9s de libert\u00e9 r\u00e9siduels propre. Le tableau d&rsquo;analyse compar\u00e9e consistera \u00e0 \u00e9valuer comment chacune de ces configurations th\u00e9oriques se positionne par rapport au mod\u00e8le complet de r\u00e9f\u00e9rence.<\/p>\n<h3>5.3 Contr\u00f4le des sorties interm\u00e9diaires<\/h3>\n<p>Avant d&rsquo;engager le calcul formel du Cp, l&rsquo;investigateur doit soumettre chaque mod\u00e8le candidat \u00e0 un examen rigoureux via la fonction <code>summary()<\/code>. Cet examen interm\u00e9diaire ne vise pas \u00e0 s\u00e9lectionner les variables sur la base exclusive de leur significativit\u00e9 marginale (valeur de <em>p<\/em> individuelle du test de Student), car cette pratique courante souffre notoirement du pi\u00e8ge de l&rsquo;omission conjointe : deux pr\u00e9dicteurs individuellement non significatifs en raison d&rsquo;une corr\u00e9lation mutuelle mod\u00e9r\u00e9e peuvent s&rsquo;av\u00e9rer conjointement indispensables pour r\u00e9duire la variance d&rsquo;erreur.<\/p>\n<p>L&rsquo;inspection des sorties interm\u00e9diaires permet de relever les coefficients de d\u00e9termination classiques R<sup>2<\/sup> et ajust\u00e9s R<sup>2<\/sup><sub>adj<\/sub>. On observe syst\u00e9matiquement que le R<sup>2<\/sup> ordinaire progresse de fa\u00e7on monotone \u00e0 mesure que l&rsquo;on \u00e9largit la sp\u00e9cification param\u00e9trique : le Mod\u00e8le 4 affichera invariablement un R<sup>2<\/sup> sup\u00e9rieur au Mod\u00e8le 3, qui lui-m\u00eame surpassera les mod\u00e8les 1 et 2. Cependant, cette progression brute est incapable d&rsquo;indiquer si le gain explicatif compense la perte de parcimonie.<\/p>\n<p>Le contr\u00f4le dimensionnel s&rsquo;op\u00e8re en extrayant le nombre effectif de param\u00e8tres pour chaque mod\u00e8le. Sous R, la dimension <em>p<\/em> d&rsquo;un mod\u00e8le estim\u00e9 par <code>lm()<\/code> s&rsquo;obtient avec exactitude via la longueur du vecteur des coefficients : <code>length(coef(m))<\/code>. Il convient de s&rsquo;assurer scrupuleusement que ce d\u00e9compte inclut bien le terme <code>(Intercept)<\/code>, conform\u00e9ment aux fondements math\u00e9matiques du crit\u00e8re explicit\u00e9s \u00e0 la section 3.2.<\/p>\n<h2>6. Calcul du Cp de Mallows avec le package olsrr et la fonction ols_mallows_cp()<\/h2>\n<h3>6.1 Mise en \u0153uvre technique de ols_mallows_cp()<\/h3>\n<p>Le package <em>olsrr<\/em> a standardis\u00e9 l&rsquo;extraction du crit\u00e8re de Mallows au travers d&rsquo;une fonction d\u00e9di\u00e9e, au nom explicite et \u00e0 la signature syntaxique intuitive : <code>ols_mallows_cp()<\/code>. Cette fonction mat\u00e9rialise l&rsquo;interaction computationnelle directe entre un mod\u00e8le candidat restreint et le mod\u00e8le satur\u00e9 de r\u00e9f\u00e9rence. Sa signature logicielle requiert deux arguments positionnels obligatoires : l&rsquo;objet <code>lm<\/code> repr\u00e9sentant le sous-mod\u00e8le analys\u00e9 et l&rsquo;objet <code>lm<\/code> incarnant le mod\u00e8le de r\u00e9f\u00e9rence complet.<\/p>\n<p>L&rsquo;appel canonique s&rsquo;ex\u00e9cute selon le sch\u00e9ma syntaxique suivant :<\/p>\n<p><code>ols_mallows_cp(model, full_model)<\/code><\/p>\n<p>Lors de l&rsquo;\u00e9valuation de cette expression, R ex\u00e9cute en arri\u00e8re-plan une s\u00e9quence d&rsquo;op\u00e9rations vectorielles hautement optimis\u00e9es : il extrait la d\u00e9viance r\u00e9siduelle du sous-mod\u00e8le, calcule la variance r\u00e9siduelle non biais\u00e9e du mod\u00e8le satur\u00e9, quantifie le nombre d&rsquo;observations effectives (<em>n<\/em>) ainsi que la dimension param\u00e9trique (<em>p<\/em>) de l&rsquo;\u00e9quation r\u00e9duite, puis applique rigoureusement l&rsquo;algorithme :<\/p>\n<p>Cp = (deviance(model) \/ (deviance(full_model) \/ df.residual(full_model))) &#8211; nobs(model) + 2 * length(coef(model))<\/p>\n<p>La sortie renvoy\u00e9e par la fonction se pr\u00e9sente sous la forme d&rsquo;une valeur scalaire num\u00e9rique continue \u00e0 double pr\u00e9cision. Cette simplicit\u00e9 d&rsquo;utilisation autorise une int\u00e9gration ais\u00e9e au sein de boucles de traitement, de fonctions personnalis\u00e9es de la famille <em>apply<\/em> ou de flux de calcul construits avec le package <em>purrr<\/em>.<\/p>\n<h3>6.2 Exemple pas \u00e0 pas bas\u00e9 sur un jeu de donn\u00e9es de r\u00e9f\u00e9rence<\/h3>\n<p>Appliquons imm\u00e9diatement cette fonction op\u00e9rationnelle \u00e0 notre ensemble de mod\u00e8les psychom\u00e9triques concurrents initialis\u00e9s \u00e0 la section 5.2, en utilisant l&rsquo;objet <code>modele_complet<\/code> comme \u00e9talon de r\u00e9f\u00e9rence pour la variance de l&rsquo;erreur :<\/p>\n<p><code>cp_m1 &lt;- ols_mallows_cp(m1, modele_complet)<\/code><br \/>\n<code>cp_m2 &lt;- ols_mallows_cp(m2, modele_complet)<\/code><br \/>\n<code>cp_m3 &lt;- ols_mallows_cp(m3, modele_complet)<\/code><br \/>\n<code>cp_m4 &lt;- ols_mallows_cp(m4, modele_complet)<\/code><\/p>\n<p>Pour structurer l&rsquo;analyse comparative et faciliter l&rsquo;interpr\u00e9tation substantielle des r\u00e9sultats, synth\u00e9tisons ces sorties num\u00e9riques au sein d&rsquo;un tableau r\u00e9capitulatif structur\u00e9 (<em>tibble<\/em> ou <em>data.frame<\/em>) int\u00e9grant la dimension <em>p<\/em> correspondante ainsi que l&rsquo;\u00e9cart absolu |Cp &#8211; p| :<\/p>\n<p><code>resultats_cp &lt;- data.frame(<br \/>\n  Modele = c(\"M1: Bio-contexte\", \"M2: Cognitivo-emotif\", \"M3: Structurel congruent\", \"M4: Quasi-complet\"),<br \/>\n  Parametres_p = c(length(coef(m1)), length(coef(m2)), length(coef(m3)), length(coef(m4))),<br \/>\n  Mallows_Cp = c(cp_m1, cp_m2, cp_m3, cp_m4)<br \/>\n)<\/code><br \/>\n<code>resultats_cp<span class=\"arabpsy-math-inline\" data-latex=\"Ecart_Absolu &lt;- \\abs(resultats_cp\">$Ecart_Absolu &lt;- \\abs(resultats_cp$<\/span>Mallows_Cp - resultats_cp$Parametres_p)<\/code><br \/>\n<code>print(resultats_cp)<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de cette s\u00e9quence met en relief la puissance discriminative du crit\u00e8re. Le Mod\u00e8le 1 (p = 4) produit une valeur de Cp massivement \u00e9lev\u00e9e (d\u00e9passant typiquement le seuil de 80), traduisant un biais d&rsquo;omission s\u00e9v\u00e8re cons\u00e9cutif \u00e0 l&rsquo;exclusion du stress et du n\u00e9vrosisme. Le Mod\u00e8le 2 (p = 4), bien que conceptuellement s\u00e9duisant, exhibe \u00e9galement un Cp nettement sup\u00e9rieur \u00e0 4 en raison de l&rsquo;absence du soutien social et du sommeil. Le Mod\u00e8le 4 (p = 7) affiche quant \u00e0 lui un Cp tr\u00e8s proche de sa dimension (approchant 7.0), mais au prix d&rsquo;une sur-param\u00e9trisation inutile. Enfin, le Mod\u00e8le 3 (p = 5) se singularise par une valeur de Cp gravitant rigoureusement autour de 4.8 \u00e0 5.1, minimisant simultan\u00e9ment le crit\u00e8re global et l&rsquo;\u00e9cart |Cp &#8211; p|, confirmant empiriquement la structure g\u00e9n\u00e9ratrice sous-jacente.<\/p>\n<h3>6.3 Erreurs classiques de manipulation et r\u00e9solution<\/h3>\n<p>Lors de l&rsquo;exploitation de la fonction <code>ols_mallows_cp()<\/code> en contexte appliqu\u00e9, plusieurs erreurs de manipulation fr\u00e9quentes sont susceptibles d&rsquo;interrompre l&rsquo;ex\u00e9cution du script ou de corrompre silencieusement la validit\u00e9 des inf\u00e9rences :<\/p>\n<ul>\n<li><strong>L&rsquo;inversion involontaire des arguments fonctionnels :<\/strong> Une m\u00e9prise classique consiste \u00e0 sp\u00e9cifier le mod\u00e8le complet en premier argument et le sous-mod\u00e8le en second argument (<code>ols_mallows_cp(modele_complet, m3)<\/code>). Dans ce cas, R calcule un Cp math\u00e9matiquement absurde o\u00f9 le mod\u00e8le complet est p\u00e9nalis\u00e9 par rapport \u00e0 la variance d&rsquo;un mod\u00e8le sous-sp\u00e9cifi\u00e9. L&rsquo;analyste doit m\u00e9moriser la s\u00e9quence imp\u00e9rative : premier argument = mod\u00e8le r\u00e9duit \u00e0 tester ; second argument = mod\u00e8le complet satur\u00e9.<\/li>\n<li><strong>La divergence du nombre d&rsquo;observations (taille d&rsquo;\u00e9chantillon in\u00e9gale) :<\/strong> Si le jeu de donn\u00e9es initial contient des valeurs manquantes non trait\u00e9es et que le mod\u00e8le complet int\u00e8gre une variable pr\u00e9sentant des donn\u00e9es manquantes que le sous-mod\u00e8le ne contient pas, la fonction <code>lm()<\/code> appliquera des suppressions d&rsquo;observations asym\u00e9triques. R renverra alors un message d&rsquo;erreur explicite signalant une non-concordance de dimensionnalit\u00e9 vectorielle : <em>\u00ab models were not fitted to the same size of dataset \u00bb<\/em>. La solution consiste \u00e0 cr\u00e9er un sous-ensemble pr\u00e9alablement filtr\u00e9 par <code>na.omit()<\/code> sur l&rsquo;ensemble des variables du mod\u00e8le satur\u00e9.<\/li>\n<li><strong>L&rsquo;omission de termes d&rsquo;interaction dans le mod\u00e8le de r\u00e9f\u00e9rence :<\/strong> Si le chercheur \u00e9value un sous-mod\u00e8le incorporant une interaction multiplicative (ex. <code>stress * soutien_social<\/code>), mais que cette interaction a \u00e9t\u00e9 omise de la d\u00e9claration du mod\u00e8le complet, l&rsquo;hypoth\u00e8se de non-biais du mod\u00e8le satur\u00e9 s&rsquo;effondre. Le mod\u00e8le de r\u00e9f\u00e9rence doit obligatoirement inclure l&rsquo;espace vectoriel le plus vaste envisag\u00e9 dans l&rsquo;\u00e9tude.<\/li>\n<\/ul>\n<h2>7. Approches alternatives sous R : calcul manuel et packages compl\u00e9mentaires<\/h2>\n<h3>7.1 Impl\u00e9mentation d&rsquo;une fonction R personnalis\u00e9e pour le calcul de Cp<\/h3>\n<p>Bien que les biblioth\u00e8ques sp\u00e9cialis\u00e9es offrent un confort op\u00e9ratoire appr\u00e9ciable, la ma\u00eetrise conceptuelle approfondie d&rsquo;un crit\u00e8re statistique s&rsquo;enrichit consid\u00e9rablement de sa programmation sous forme d&rsquo;algorithme matriciel ou vectoriel autonome (<em>from scratch<\/em>). L&rsquo;impl\u00e9mentation d&rsquo;une fonction personnalis\u00e9e sous R permet de s&rsquo;affranchir des d\u00e9pendances logicielles externes, d&rsquo;acc\u00e9l\u00e9rer les temps de calcul lors de simulations Monte Carlo intensives et d&rsquo;int\u00e9grer des routines de contr\u00f4le d&rsquo;int\u00e9grit\u00e9 sur mesure.<\/p>\n<p>Consid\u00e9rons l&rsquo;\u00e9criture d&rsquo;une fonction robuste baptis\u00e9e <code>custom_mallows_cp()<\/code>, con\u00e7ue pour extraire les scalaires analytiques indispensables directement depuis les environnements d&rsquo;\u00e9valuation des objets <code>lm<\/code> :<\/p>\n<p><code>custom_mallows_cp &lt;- function(sub_model, full_model) {<br \/>\n  # Verification formelle de la classe des objets en entree<br \/>\n  if (!inherits(sub_model, \"lm\") || !inherits(full_model, \"lm\")) {<br \/>\n    stop(\"Les deux arguments doivent obligatoirement etre des objets de classe 'lm'.\")<br \/>\n  }<\/p>\n<p>  # Controle de l'identite rigoureuse du nombre d'observations<br \/>\n  n_sub &lt;- stats::nobs(sub_model)<br \/>\n  n_full &lt;- stats::nobs(full_model)<br \/>\n  if (n_sub != n_full) {<br \/>\n    stop(paste(\"Incompatibilite de taille d'echantillon : modele reduit n =\", n_sub,<br \/>\n               \"| modele complet n =\", n_full))<br \/>\n  }<\/p>\n<p>  # Extraction de la somme des carres des residus (RSS) et des degres de liberte<br \/>\n  rss_p &lt;- stats::deviance(sub_model)<br \/>\n  rss_full &lt;- stats::deviance(full_model)<br \/>\n  df_full &lt;- stats::df.residual(full_model)<\/p>\n<p>  # Estimation non biaisee de la variance de l'erreur residuelle du modele complet<br \/>\n  sigma2_est &lt;- rss_full \/ df_full<\/p>\n<p>  # Nombre total de parametres du sous-modele (incluant l'intercept)<br \/>\n  p &lt;- length(stats::coef(sub_model))<\/p>\n<p>  # Application formelle de la formule algebrique de Mallows<br \/>\n  cp_val &lt;- (rss_p \/ sigma2_est) - n_sub + (2 * p)<\/p>\n<p>  return(cp_val)<br \/>\n}<\/code><\/p>\n<p>Cette impl\u00e9mentation native garantit une transparence totale sur le flux de calcul. La v\u00e9rification crois\u00e9e entre <code>custom_mallows_cp(m3, modele_complet)<\/code> et <code>ols_mallows_cp(m3, modele_complet)<\/code> produit des r\u00e9sultats rigoureusement identiques au niveau de pr\u00e9cision de la virgule flottante, attestant de la concordance th\u00e9orique absolue du processus.<\/p>\n<h3>7.2 Utilisation de la biblioth\u00e8que leaps pour la r\u00e9gression exhaustive<\/h3>\n<p>Lorsque le nombre de variables pr\u00e9dictives potentielles exc\u00e8de quatre ou cinq r\u00e9gresseurs, l&rsquo;ajustement manuel individuel de toutes les combinaisons possibles devient rapidement irr\u00e9alisable sur le plan combinatoire. Pour un ensemble de <em>K<\/em> pr\u00e9dicteurs, le nombre total d&rsquo;\u00e9quations candidates envisageables s&rsquo;\u00e9l\u00e8ve \u00e0 2<sup>K<\/sup> (soit 128 mod\u00e8les pour K = 7, mais 1 048 576 pour K = 20). Pour relever ce d\u00e9fi computationnel, la biblioth\u00e8que <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=leaps\">leaps<\/a> impl\u00e9mente l&rsquo;algorithme hautement efficace de Furnival et Wilson (1974), reposant sur une exploration arborescente par s\u00e9paration et \u00e9valuation progressive (<em>branch-and-bound<\/em>).<\/p>\n<p>La fonction ma\u00eetresse <code>regsubsets()<\/code> permet d&rsquo;ex\u00e9cuter cette recherche combinatoire int\u00e9grale en une seule ligne de code, en isolant instantan\u00e9ment la meilleure combinaison de pr\u00e9dicteurs pour chaque taille de sous-mod\u00e8le donn\u00e9e :<\/p>\n<p><code>selection_exhaustive &lt;- regsubsets(epuisement ~ ., data = df_psycho, nvmax = 7, method = \"exhaustive\")<\/code><br \/>\n<code>resume_leaps &lt;- summary(selection_exhaustive)<\/code><\/p>\n<p>L&rsquo;objet de synth\u00e8se produit par <code>summary(selection_exhaustive)<\/code> rec\u00e8le une matrice logique d\u00e9taillant les pr\u00e9dicteurs engag\u00e9s, ainsi qu&rsquo;un vecteur num\u00e9rique contenant les valeurs de Cp associ\u00e9es \u00e0 la meilleure combinaison identifi\u00e9e pour chaque strate de dimensionnalit\u00e9 :<\/p>\n<p><code>resume_leaps$cp<\/code><\/p>\n<p>Pour convertir ces sorties vectorielles brutes en un tableau analytique contemporain et manipulable, nous pouvons structurer les r\u00e9sultats \u00e0 l&rsquo;aide des fonctions de manipulation de donn\u00e9es modernes :<\/p>\n<p><code>synthese_leaps &lt;- tibble(<br \/>\n  Nombre_Predicteurs = 1:length(resume_leaps$cp),<br \/>\n  Parametres_p = (1:length(resume_leaps$cp)) + 1, # Prise en compte de la constante<br \/>\n  Cp_Mallows = resume_leaps$cp,<br \/>\n  R2 = resume_leaps$rsq,<br \/>\n  R2_Ajuste = resume_leaps$adjr2,<br \/>\n  BIC = resume_leaps$bic<br \/>\n)<\/code><br \/>\n<code>print(synthese_leaps)<\/code><\/p>\n<p>Cette approche syst\u00e9matique met en \u00e9vidence la trajectoire d&rsquo;\u00e9volution du Cp \u00e0 travers l&rsquo;espace des dimensions, identifiant instantan\u00e9ment \u00e0 quel niveau de cardinalit\u00e9 le crit\u00e8re plonge pour la premi\u00e8re fois \u00e0 proximit\u00e9 imm\u00e9diate de la ligne d&rsquo;\u00e9quipartition <em>p<\/em>.<\/p>\n<h3>7.3 S\u00e9lection multimorphique avec MuMIn et car<\/h3>\n<p>L&rsquo;\u00e9cosyst\u00e8me logiciel de R dispose d&rsquo;autres extensions sophistiqu\u00e9es d\u00e9di\u00e9es \u00e0 l&rsquo;inf\u00e9rence multimorphique et \u00e0 la s\u00e9lection d&rsquo;\u00e9quations structurales. Parmi celles-ci, le package <em>MuMIn<\/em> (Multi-Model Inference) fournit un cadre de travail avanc\u00e9 pour le classement, le partitionnement et la combinaison bay\u00e9sienne ou fr\u00e9quentiste de mod\u00e8les de r\u00e9gression. Bien que <em>MuMIn<\/em> soit principalement orient\u00e9 vers le calcul de l&rsquo;AICc (crit\u00e8re d&rsquo;Akaike corrig\u00e9 pour petits \u00e9chantillons), il permet d&rsquo;int\u00e9grer des crit\u00e8res de p\u00e9nalisation personnalis\u00e9s au sein de ses tables d&rsquo;ordonnancement global.<\/p>\n<p>La fonction <code>dredge()<\/code> de <em>MuMIn<\/em> permet de g\u00e9n\u00e9rer automatiquement la totalit\u00e9 des sous-mod\u00e8les d\u00e9riv\u00e9s d&rsquo;un mod\u00e8le global pr\u00e9alablement ajust\u00e9. Elle requiert toutefois que l&rsquo;argument d&rsquo;ajustement global soit configur\u00e9 de fa\u00e7on \u00e0 pr\u00e9server l&rsquo;acc\u00e8s int\u00e8gre aux donn\u00e9es d&rsquo;origine (via l&rsquo;instruction <code>options(na.action = \"na.fail\")<\/code>) :<\/p>\n<p><code>library(MuMIn)<\/code><br \/>\n<code>options(na.action = \"na.fail\")<\/code><br \/>\n<code>modele_global &lt;- lm(epuisement ~ ., data = df_psycho)<\/code><br \/>\n<code>table_selection &lt;- dredge(modele_global)<\/code><br \/>\n<code>head(table_selection, 5)<\/code><\/p>\n<p>Parall\u00e8lement, la biblioth\u00e8que <em>car<\/em> apporte un support m\u00e9thodologique d\u00e9terminant en permettant de tester si l&rsquo;omission conjointe d&rsquo;un ensemble de pr\u00e9dicteurs \u00e9cart\u00e9s par l&rsquo;arbitrage du Cp induit une d\u00e9gradation de variance statistiquement significative. L&#8217;emploi conjoint de la fonction <code>linearHypothesis()<\/code> ou de l&rsquo;analyse de variance comparative <code>anova(modele_reduit, modele_complet)<\/code> permet ainsi d&rsquo;\u00e9tablir un pont diagnostique formel entre la r\u00e8gle d\u00e9cisionnelle heuristique de Colin Mallows et les tests de rapport de variance de Fisher-Snedecor traditionnels.<\/p>\n<h2>8. Proc\u00e9dure de s\u00e9lection automatis\u00e9e par sous-ensembles et comparaison graphique<\/h2>\n<h3>8.1 Ex\u00e9cution de la s\u00e9lection globale avec ols_step_best_subset()<\/h3>\n<p>La biblioth\u00e8que <em>olsrr<\/em> propose une interface d&rsquo;\u00e9valuation automatis\u00e9e de niveau sup\u00e9rieur gr\u00e2ce \u00e0 la fonction <code>ols_step_best_subset()<\/code>. Contrairement aux approches s\u00e9quentielles pas \u00e0 pas (<em>stepwise regression<\/em>) qui souffrent de limitations analytiques s\u00e9v\u00e8res en pr\u00e9sence de colin\u00e9arit\u00e9, cette routine examine rigoureusement l&rsquo;ensemble des 2<sup>K<\/sup> configurations possibles et synth\u00e9tise les m\u00e9triques d&rsquo;\u00e9valuation d\u00e9cisionnelles au sein d&rsquo;une structure de donn\u00e9es unifi\u00e9e.<\/p>\n<p>L&rsquo;ex\u00e9cution de la proc\u00e9dure s&rsquo;op\u00e8re directement \u00e0 partir du mod\u00e8le complet initial :<\/p>\n<p><code>selection_sous_ensembles &lt;- ols_step_best_subset(modele_complet)<\/code><br \/>\n<code>print(selection_sous_ensembles)<\/code><\/p>\n<p>La console affiche un rapport tabulaire d&rsquo;une richesse exceptionnelle. Pour chaque niveau de dimension param\u00e9trique (de 1 \u00e0 K pr\u00e9dicteurs), la routine consigne la combinaison gagnante, son coefficient de d\u00e9termination R<sup>2<\/sup>, son R<sup>2<\/sup> ajust\u00e9, le crit\u00e8re d&rsquo;information d&rsquo;Akaike (AIC), le crit\u00e8re de Schwarz-Bayes (BIC), le crit\u00e8re de Sawa, et de fa\u00e7on centrale, le <strong>Cp de Mallows<\/strong>. L&rsquo;analyste peut ainsi retracer visuellement le cheminement de la p\u00e9nalisation : \u00e0 mesure que le mod\u00e8le int\u00e8gre les variables structurellement d\u00e9terminantes (le stress, puis le n\u00e9vrosisme, puis le soutien social, puis le sommeil), le Cp s&rsquo;effondre de valeurs stratosph\u00e9riques (sup\u00e9rieures \u00e0 200 pour le mod\u00e8le \u00e0 une variable) jusqu&rsquo;\u00e0 atteindre un point d&rsquo;\u00e9quilibre stabilis\u00e9.<\/p>\n<h3>8.2 Visualisation graphique : le diagramme de Mallows<\/h3>\n<p>L&rsquo;apport m\u00e9thodologique fondamental de Colin Mallows en 1973 ne r\u00e9sidait pas uniquement dans la d\u00e9rivation alg\u00e9brique de sa formule, mais dans la proposition d&rsquo;un dispositif d&rsquo;exploration visuelle d\u00e9sormais c\u00e9l\u00e8bre : le <strong>diagramme de Cp<\/strong> (ou <em>Cp plot<\/em>). Ce graphique cartographie sur l&rsquo;axe des ordonn\u00e9es la valeur num\u00e9rique de la statistique Cp calcul\u00e9e pour chaque mod\u00e8le candidat, et sur l&rsquo;axe des abscisses le nombre de param\u00e8tres totaux <em>p<\/em> engag\u00e9s dans l&rsquo;\u00e9quation.<\/p>\n<p>La fonction g\u00e9n\u00e9rique int\u00e9gr\u00e9e \u00e0 <em>olsrr<\/em> g\u00e9n\u00e8re automatiquement ce type de rendu :<\/p>\n<p><code>plot(selection_sous_ensembles)<\/code><\/p>\n<p>Cependant, afin de r\u00e9pondre aux exigences de clart\u00e9 et de personnalisation des publications scientifiques \u00e0 fort facteur d&rsquo;impact, il est hautement pr\u00e9f\u00e9rable de concevoir ce trac\u00e9 directement sous l&rsquo;environnement <em>ggplot2<\/em>. Ce contr\u00f4le programmatique permet d&rsquo;incorporer explicitement la ligne th\u00e9orique d&rsquo;\u00e9quipartition y = x symbolisant l&rsquo;absence de biais :<\/p>\n<p><code>df_graph &lt;- as.data.frame(selection_sous_ensembles$metrics)<\/code><br \/>\n<code>df_graph<span class=\"arabpsy-math-inline\" data-latex=\"p &lt;- df_graph\">$p &lt;- df_graph$<\/span>n + 1 # n representant ici le nombre de predicteurs hors intercept<\/code><\/p>\n<p><code>ggplot(df_graph, aes(x = p, y = mallows_cp)) +<br \/>\n  geom_point(color = \"#1f77b4\", size = 3.5, alpha = 0.8) +<br \/>\n  geom_abline(intercept = 0, slope = 1, linetype = \"dashed\", color = \"#d62728\", linewidth = 0.9) +<br \/>\n  geom_text(aes(label = paste0(\"p=\", p, \"nCp=\", round(mallows_cp, 1))),<br \/>\n            vjust = -0.7, size = 3, check_overlap = TRUE) +<br \/>\n  scale_x_continuous(breaks = 2:(ncol(df_psycho))) +<br \/>\n  labs(<br \/>\n    title = \"Diagramme de selection de Mallows (Cp Plot)\",<br \/>\n    subtitle = \"La droite en pointille rouge figure la reference d'absence theorique de biais (y = p)\",<br \/>\n    x = \"Nombre total de parametres estimes p (constante incluse)\",<br \/>\n    y = \"Statistique Cp de Mallows\"<br \/>\n  ) +<br \/>\n  theme_minimal(base_size = 12) +<br \/>\n  theme(<br \/>\n    panel.grid.minor = element_blank(),<br \/>\n    plot.title = element_text(face = \"bold\", hjust = 0.5),<br \/>\n    plot.subtitle = element_text(hjust = 0.5)<br \/>\n  )<\/code><\/p>\n<h3>8.3 Interpr\u00e9tation s\u00e9miologique du graphique de s\u00e9lection<\/h3>\n<p>L&rsquo;interpr\u00e9tation s\u00e9miologique de cette cartographie visuelle repose sur la position relative des points empiriques par rapport \u00e0 la droite diagonale de r\u00e9f\u00e9rence y = x :<\/p>\n<ul>\n<li><strong>Zone sup\u00e9rieure gauche (Cp &gt;&gt; p) :<\/strong> Les mod\u00e8les positionn\u00e9s haut au-dessus de la ligne diagonale souffrent d&rsquo;un biais d&rsquo;omission s\u00e9v\u00e8re. Bien que parcimonieux en termes de param\u00e8tres, leur erreur quadratique totale de pr\u00e9diction est inacceptable. Ces mod\u00e8les doivent \u00eatre cat\u00e9goriquement rejet\u00e9s.<\/li>\n<li><strong>Trajectoire d&rsquo;inflexion (Le point de coude) :<\/strong> En examinant le trac\u00e9 de gauche \u00e0 droite, l&rsquo;analyste observe g\u00e9n\u00e9ralement une chute abrupte de la courbe du Cp jusqu&rsquo;\u00e0 atteindre pour la premi\u00e8re fois la proximit\u00e9 imm\u00e9diate de la droite y = x. Ce point d&rsquo;inflexion signale la saturation de l&rsquo;information explicative majeure.<\/li>\n<li><strong>Zone de concordance optimale (Cp &asymp; p) :<\/strong> Le premier mod\u00e8le situ\u00e9 sur la ligne diagonale ou imm\u00e9diatement en dessous de celle-ci avec le plus petit nombre de param\u00e8tres identifie la sp\u00e9cification id\u00e9ale. Dans notre simulation, le point correspondant \u00e0 p = 5 (stress, nevrosisme, soutien_social, sommeil_degrade + intercept) s&rsquo;ancre pr\u00e9cis\u00e9ment sur la droite.<\/li>\n<li><strong>Plateau asymptotique \u00e0 droite (Cp &asymp; p pour p croissant) :<\/strong> \u00c0 mesure que <em>p<\/em> progresse au-del\u00e0 du point optimal vers le mod\u00e8le satur\u00e9 <em>P<\/em>, les points continuent de coller \u00e0 la diagonale. Toutefois, cette trajectoire n&rsquo;apporte aucune r\u00e9duction substantielle du biais : elle ne fait qu&rsquo;ajouter des param\u00e8tres redondants dont le co\u00fbt en variance compense st\u00e9rilement l&rsquo;infime gain r\u00e9siduel.<\/li>\n<\/ul>\n<h2>9. Analyse comparative : Cp de Mallows versus AIC, BIC et R-carr\u00e9 ajust\u00e9<\/h2>\n<h3>9.1 Propri\u00e9t\u00e9s asymptotiques face \u00e0 l&rsquo;Akaike Information Criterion (AIC)<\/h3>\n<p>Le crit\u00e8re Cp de Mallows entretient des affinit\u00e9s formelles profondes avec le c\u00e9l\u00e8bre crit\u00e8re d&rsquo;information d&rsquo;Akaike (<a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1109\/TAC.1974.1100705\">AIC<\/a>). Sous l&rsquo;hypoth\u00e8se de r\u00e9sidus gaussiens distribu\u00e9s selon une loi normale d&rsquo;esp\u00e9rance nulle et de variance constante &sigma;<sup>2<\/sup>, il est math\u00e9matiquement d\u00e9montrable que le crit\u00e8re Cp de Mallows et l&rsquo;AIC sont asymptotiquement \u00e9quivalents. En effet, l&rsquo;expression de la log-vraisemblance d&rsquo;un mod\u00e8le de r\u00e9gression lin\u00e9aire conduit \u00e0 la formulation suivante de l&rsquo;AIC :<\/p>\n<p>AIC = n &times; ln(RSS<sub>p<\/sub> \/ n) + 2p + Cte<\/p>\n<p>En effectuant un d\u00e9veloppement limit\u00e9 au premier ordre de la fonction logarithmique au voisinage de la variance de r\u00e9f\u00e9rence s<sup>2<\/sup>, la p\u00e9nalit\u00e9 de l&rsquo;AIC converge vers la fonction de p\u00e9nalisation lin\u00e9aire 2p caract\u00e9ristique du Cp de Mallows. Dans les deux cas, le co\u00fbt marginal allou\u00e9 \u00e0 l&rsquo;introduction d&rsquo;un param\u00e8tre suppl\u00e9mentaire correspond exactement \u00e0 deux unit\u00e9s de variance r\u00e9siduelle normalis\u00e9e.<\/p>\n<p>N\u00e9anmoins, des nuances op\u00e9rationnelles notables subsistent \u00e0 petite \u00e9chelle d&rsquo;\u00e9chantillonnage. Le Cp de Mallows utilise explicitement l&rsquo;estimateur sans biais de la variance r\u00e9siduelle d\u00e9riv\u00e9 du <em>mod\u00e8le complet de r\u00e9f\u00e9rence<\/em> au d\u00e9nominateur, tandis que l&rsquo;AIC classique \u00e9value la log-vraisemblance interne de chaque sous-mod\u00e8le de mani\u00e8re autonome (en estimant sa variance par le maximum de vraisemblance RSS<sub>p<\/sub> \/ n, qui constitue un estimateur biais\u00e9 pour les faibles valeurs de n). Par cons\u00e9quent, sur des \u00e9chantillons restreints, le Cp de Mallows offre fr\u00e9quemment une \u00e9valuation plus stable et moins sensible aux fluctuations locales d&rsquo;ajustement que l&rsquo;AIC non corrig\u00e9.<\/p>\n<h3>9.2 Cp de Mallows face au Bayesian Information Criterion (BIC)<\/h3>\n<p>La confrontation entre le Cp de Mallows et le crit\u00e8re d&rsquo;information bay\u00e9sien (<a rel=\"noopener noreferrer\" href=\"https:\/\/projecteuclid.org\/journals\/annals-of-statistics\/volume-6\/issue-2\/Estimating-the-Dimension-of-a-Model\/10.1214\/aos\/1176344136.full\">BIC<\/a> ou crit\u00e8re de Schwarz) illustre deux philosophies statistiques distinctes de la s\u00e9lection de mod\u00e8les. L&rsquo;\u00e9quation canonique du BIC s&rsquo;\u00e9nonce selon :<\/p>\n<p>BIC = n &times; ln(RSS<sub>p<\/sub> \/ n) + ln(n) &times; p<\/p>\n<p>La diff\u00e9rence fondamentale r\u00e9side dans la s\u00e9v\u00e9rit\u00e9 de la p\u00e9nalisation param\u00e9trique. Alors que le Cp de Mallows et l&rsquo;AIC appliquent une p\u00e9nalit\u00e9 fixe \u00e9gale \u00e0 2 unit\u00e9s par param\u00e8tre quel que soit le volume d&rsquo;observations, le BIC applique un coefficient multiplicateur logarithmique ln(n). D\u00e8s lors que la taille d&rsquo;\u00e9chantillon <em>n<\/em> exc\u00e8de e<sup>2<\/sup> &asymp; 7.39 observations (ce qui est syst\u00e9matiquement le cas en pratique), ln(n) devient strictement sup\u00e9rieur \u00e0 2. Le BIC p\u00e9nalise donc la complexit\u00e9 structurelle avec une rigueur exponentiellement croissante par rapport au Cp.<\/p>\n<p>Sur le plan des propri\u00e9t\u00e9s fondamentales, le BIC est <strong>asymptotiquement convergent<\/strong> (<em>consistent<\/em>) : si le \u00ab vrai \u00bb mod\u00e8le g\u00e9n\u00e9rateur de donn\u00e9es appartient formellement \u00e0 la famille des sous-mod\u00e8les examin\u00e9s, la probabilit\u00e9 que le BIC s\u00e9lectionne ce mod\u00e8le exact tend vers 1 lorsque <em>n<\/em> tend vers l&rsquo;infini. En revanche, le Cp de Mallows est <strong>efficace pour la pr\u00e9diction<\/strong> (<em>asymptotically efficient<\/em>) : il minimise l&rsquo;erreur quadratique moyenne de pr\u00e9diction sur de nouvelles observations, quitte \u00e0 conserver une l\u00e9g\u00e8re propension \u00e0 sur-s\u00e9lectionner des pr\u00e9dicteurs marginaux sur les tr\u00e8s grands \u00e9chantillons. En psychologie et en \u00e9pid\u00e9miologie, o\u00f9 les ph\u00e9nom\u00e8nes comportementaux sont multifactoriels et o\u00f9 la notion de \u00ab vrai \u00bb mod\u00e8le parcimonieux est une abstraction th\u00e9orique, l&rsquo;orientation pr\u00e9dictive du Cp s&rsquo;av\u00e8re souvent hautement pragmatique.<\/p>\n<h3>9.3 Limites critiques du R-carr\u00e9 ajust\u00e9<\/h3>\n<p>Le coefficient de d\u00e9termination ajust\u00e9 (R<sup>2<\/sup><sub>adj<\/sub> ou <em>adjusted R-squared<\/em>), propos\u00e9 initialement par Henri Theil, constitue traditionnellement la m\u00e9trique la plus famili\u00e8re des \u00e9tudiants et des praticiens de la r\u00e9gression lin\u00e9aire. Con\u00e7u pour pallier l&rsquo;accroissement m\u00e9canique du R<sup>2<\/sup> standard, il pond\u00e8re le rapport des sommes des carr\u00e9s par leurs degr\u00e9s de libert\u00e9 respectifs :<\/p>\n<p>R<sup>2<\/sup><sub>adj<\/sub> = 1 &#8211; [ (RSS<sub>p<\/sub> \/ (n &#8211; p)) \/ (TSS \/ (n &#8211; 1)) ]<\/p>\n<p>Malgr\u00e9 sa popularit\u00e9 p\u00e9dagogique, le R<sup>2<\/sup> ajust\u00e9 pr\u00e9sente des limites critiques structurelles qui le rendent tr\u00e8s inf\u00e9rieur au Cp de Mallows pour la s\u00e9lection rigoureuse de mod\u00e8les :<\/p>\n<ul>\n<li><strong>P\u00e9nalisation marginale insuffisante :<\/strong> On d\u00e9montre analytiquement que le R<sup>2<\/sup> ajust\u00e9 augmente d\u00e8s lors qu&rsquo;un pr\u00e9dicteur additionnel pr\u00e9sente une statistique de Student en valeur absolue strictement sup\u00e9rieure \u00e0 1 (|t| &gt; 1, soit une valeur de <em>p<\/em> approximativement inf\u00e9rieure \u00e0 0.32). Cette condition d&rsquo;inclusion est particuli\u00e8rement laxiste et conduit fr\u00e9quemment le chercheur \u00e0 retenir des pr\u00e9dicteurs purement parasitaires qui d\u00e9t\u00e9riorent la variance de pr\u00e9diction.<\/li>\n<li><strong>Absence de quantification explicite du biais :<\/strong> Le R<sup>2<\/sup> ajust\u00e9 normalise la variance intra-mod\u00e8le sans jamais la confronter \u00e0 un \u00e9talon de r\u00e9f\u00e9rence ext\u00e9rieur impartial. Il est incapable d&rsquo;indiquer si le sous-mod\u00e8le retenu g\u00e9n\u00e8re un biais syst\u00e9matique d&rsquo;omission.<\/li>\n<li><strong>Divergences fr\u00e9quentes en pr\u00e9sence de bruit :<\/strong> Dans des contextes de colin\u00e9arit\u00e9 mod\u00e9r\u00e9e, le R<sup>2<\/sup> ajust\u00e9 tend \u00e0 s\u00e9lectionner le mod\u00e8le contenant le plus grand nombre de variables, l\u00e0 o\u00f9 le diagramme de Mallows met en exergue la stabilisation du Cp bien en amont de cette saturation artificielle.<\/li>\n<\/ul>\n<h2>10. Cas pratique approfondi en recherche psychologique et comportementale<\/h2>\n<h3>10.1 Pr\u00e9sentation du cas d&rsquo;\u00e9tude et hypoth\u00e8ses scientifiques<\/h3>\n<p>Afin d&rsquo;ancrer ces concepts dans la r\u00e9alit\u00e9 d&rsquo;un laboratoire de recherche quantitative, d\u00e9veloppons une \u00e9tude de cas d\u00e9taill\u00e9e portant sur l&rsquo;\u00e9valuation psychom\u00e9trique de la d\u00e9tresse psychologique en contexte post-traumatique ou de stress chronique. Une \u00e9quipe de chercheurs souhaite concevoir un outil de d\u00e9pistage rapide (<em>screening<\/em>) permettant de pr\u00e9dire le score d&rsquo;anxi\u00e9t\u00e9 clinique g\u00e9n\u00e9ralis\u00e9e (mesur\u00e9 par l&rsquo;\u00e9chelle continue GAD-7) aupr\u00e8s d&rsquo;une cohorte de patients admis en consultation de soutien.<\/p>\n<p>Le protocole d&rsquo;admission initial administre un protocole exhaustif documentant dix pr\u00e9dicteurs psychom\u00e9triques et physiologiques potentiels :<\/p>\n<ul>\n<li><em>N\u00e9vrosisme<\/em> (inventaire NEO-PI-R)<\/li>\n<li><em>R\u00e9activit\u00e9 \u00e9motionnelle<\/em> (score composite d&rsquo;intensit\u00e9 affective)<\/li>\n<li><em>Soutien familial per\u00e7u<\/em> (\u00e9chelle MSPSS)<\/li>\n<li><em>Soutien amical per\u00e7u<\/em> (\u00e9chelle MSPSS)<\/li>\n<li><em>Strat\u00e9gies de coping actif<\/em> (\u00e9chelle Brief-COPE)<\/li>\n<li><em>Strat\u00e9gies de coping par \u00e9vitement<\/em> (\u00e9chelle Brief-COPE)<\/li>\n<li><em>Fr\u00e9quence cardiaque au repos<\/em> (biomarqueur continu)<\/li>\n<li><em>Dur\u00e9e moyenne de sommeil paradoxal<\/em> (mesure polysomnographique)<\/li>\n<li><em>Niveau de revenu annuel d\u00e9clar\u00e9<\/em> (variable socio-\u00e9conomique)<\/li>\n<li><em>\u00c2ge chronologique<\/em> du participant<\/li>\n<\/ul>\n<p>L&rsquo;hypoth\u00e8se scientifique fondamentale postule que la d\u00e9tresse anxieuse clinique peut \u00eatre mod\u00e9lis\u00e9e avec une pr\u00e9cision pr\u00e9dictive optimale au moyen d&rsquo;un sous-ensemble restreint combinant le n\u00e9vrosisme, le coping par \u00e9vitement et la perturbation du sommeil, sans qu&rsquo;il soit n\u00e9cessaire de mobiliser les mesures biom\u00e9triques co\u00fbteuses ou les indicateurs socio-\u00e9conomiques p\u00e9riph\u00e9riques. L&rsquo;objectif m\u00e9thodologique consiste \u00e0 appliquer la proc\u00e9dure de Mallows pour arbitrer rigoureusement cette hypoth\u00e8se.<\/p>\n<h3>10.2 Pipeline de traitement R complet pas \u00e0 pas<\/h3>\n<p>D\u00e9ployons l&rsquo;int\u00e9gralit\u00e9 du pipeline computationnel sous R, depuis la simulation contr\u00f4l\u00e9e du jeu de donn\u00e9es cliniques (n = 300) jusqu&rsquo;\u00e0 la s\u00e9lection automatis\u00e9e guid\u00e9e par le Cp :<\/p>\n<p><code># 1. Generation de la matrice clinique simulee<\/code><br \/>\n<code>set.seed(2024)<\/code><br \/>\n<code>n_clinique &lt;- 300<\/code><br \/>\n<code>nevrosisme &lt;- rnorm(n_clinique, mean = 55, sd = 12)<\/code><br \/>\n<code>reactivite &lt;- rnorm(n_clinique, mean = 50, sd = 10) + 0.35 * nevrosisme<\/code><br \/>\n<code>soutien_famille &lt;- rnorm(n_clinique, mean = 28, sd = 6)<\/code><br \/>\n<code>soutien_amis &lt;- rnorm(n_clinique, mean = 26, sd = 5) + 0.25 * soutien_famille<\/code><br \/>\n<code>coping_actif &lt;- rnorm(n_clinique, mean = 32, sd = 7)<\/code><br \/>\n<code>coping_evitement &lt;- rnorm(n_clinique, mean = 24, sd = 6) + 0.3 * nevrosisme<\/code><br \/>\n<code>freq_cardiaque &lt;- rnorm(n_clinique, mean = 72, sd = 9)<\/code><br \/>\n<code>sommeil_rem &lt;- rnorm(n_clinique, mean = 90, sd = 20) - 0.4 * nevrosisme<\/code><br \/>\n<code>revenu &lt;- runif(n_clinique, min = 15, max = 80)<\/code><br \/>\n<code>age &lt;- rnorm(n_clinique, mean = 42, sd = 11)<\/code><\/p>\n<p><code># 2. Generation de la variable reponse : Anxiete GAD-7<\/code><br \/>\n<code>anxiete_gad7 &lt;- 10 + 0.45 * nevrosisme + 0.65 * coping_evitement - 0.15 * sommeil_rem + rnorm(n_clinique, mean = 0, sd = 6)<\/code><\/p>\n<p><code># 3. Assemblage du dataframe d'etude<\/code><br \/>\n<code>df_clinique &lt;- data.frame(<br \/>\n  anxiete_gad7, nevrosisme, reactivite, soutien_famille, soutien_amis,<br \/>\n  coping_actif, coping_evitement, freq_cardiaque, sommeil_rem, revenu, age<br \/>\n)<\/code><\/p>\n<p><code># 4. Ajustement du modele complet sature<\/code><br \/>\n<code>mod_complet_clinique &lt;- lm(anxiete_gad7 ~ ., data = df_clinique)<\/code><\/p>\n<p><code># 5. Selection exhaustive des sous-ensembles optimaux<\/code><br \/>\n<code>selection_clinique &lt;- ols_step_best_subset(mod_complet_clinique)<\/code><\/p>\n<p><code># 6. Extraction et structuration des metriques comparatives<\/code><br \/>\n<code>metrics_tbl &lt;- as_tibble(selection_clinique$metrics) %&gt;%<br \/>\n  mutate(p_total = n + 1,<br \/>\n         ecart_mallows = abs(mallows_cp - p_total)) %&gt;%<br \/>\n  select(n, p_total, predictors, rsquare, adjr, mallows_cp, ecart_mallows, aic, bic)<\/code><\/p>\n<p><code># 7. Isolement de la solution minimisant le compromis de Mallows<\/code><br \/>\n<code>modele_selectionne &lt;- metrics_tbl %&gt;%<br \/>\n  filter(mallows_cp &lt;= p_total + 1.5) %&gt;%<br \/>\n  arrange(p_total) %&gt;%<br \/>\n  slice(1)<\/code><\/p>\n<p><code>print(modele_selectionne)<\/code><\/p>\n<h3>10.3 Synth\u00e8se et discussion des r\u00e9sultats psychom\u00e9triques<\/h3>\n<p>L&rsquo;inspection analytique de la table g\u00e9n\u00e9r\u00e9e met en \u00e9vidence des constats psychom\u00e9triques d&rsquo;une grande limpidit\u00e9. Pour n = 1 r\u00e9gresseur (le n\u00e9vrosisme seul, p = 2), la statistique Cp s&rsquo;\u00e9tablit \u00e0 une valeur massivement prohibitive (Cp &gt; 120), d\u00e9montrant sans \u00e9quivoque que la prise en compte exclusive du trait de personnalit\u00e9 n&rsquo;\u00e9puise nullement la variance clinique de l&rsquo;anxi\u00e9t\u00e9. Pour n = 2 r\u00e9gresseurs (n\u00e9vrosisme et coping par \u00e9vitement, p = 3), le Cp r\u00e9gresse spectaculairement autour de 22.4, mais demeure substantiellement au-dessus du seuil de non-biais.<\/p>\n<p>La bascule d\u00e9cisive s&rsquo;op\u00e8re rigoureusement \u00e0 <strong>n = 3 pr\u00e9dicteurs (p = 4 param\u00e8tres)<\/strong>, correspondant \u00e0 l&rsquo;int\u00e9gration conjointe de <em>nevrosisme<\/em>, <em>coping_evitement<\/em> et <em>sommeil_rem<\/em>. \u00c0 cet \u00e9chelon pr\u00e9cis, le crit\u00e8re Cp de Mallows chute \u00e0 <strong>3.78<\/strong>. Cette valeur est \u00e0 la fois remarquablement basse et strictement inf\u00e9rieure \u00e0 la dimension param\u00e9trique nominale (p = 4), attestant d&rsquo;une absence compl\u00e8te de biais de sp\u00e9cification coupl\u00e9e \u00e0 une variance minimale.<\/p>\n<p>L&rsquo;examen des strates ult\u00e9rieures (de n = 4 \u00e0 n = 10 pr\u00e9dicteurs) confirme la pertinence de la r\u00e8gle de d\u00e9cision. L&rsquo;adjonction de la fr\u00e9quence cardiaque, du revenu ou des facettes de soutien social ne provoque qu&rsquo;une r\u00e9duction marginale et stochastique de la somme des carr\u00e9s r\u00e9siduelle : le Cp oscille st\u00e9rilement entre 4.8 et 11.0 en collant \u00e0 la droite d&rsquo;\u00e9quipartition sans jamais am\u00e9liorer le pouvoir explicatif net. Gr\u00e2ce au crit\u00e8re de Mallows, le protocole clinique peut valablement \u00e9laguer 7 des 10 questionnaires initiaux, all\u00e9geant drastiquement le fardeau cognitif impos\u00e9 aux patients sans aucune d\u00e9gradation de la validit\u00e9 de d\u00e9pistage.<\/p>\n<h2>11. Diagnostics r\u00e9siduels et validation post-s\u00e9lection du mod\u00e8le retenu<\/h2>\n<h3>11.1 Validation formelle des postulats de r\u00e9gression gaussienne<\/h3>\n<p>L&rsquo;isolement d&rsquo;un mod\u00e8le parcimonieux par le crit\u00e8re Cp de Mallows ne dispense aucunement le statisticien de l&rsquo;obligation m\u00e9thodologique de soumettre la sp\u00e9cification retenue \u00e0 une batterie compl\u00e8te de diagnostics r\u00e9siduels. L&rsquo;ensemble de la d\u00e9rivation math\u00e9matique du Cp reposant sur l&rsquo;hypoth\u00e8se de perturbations r\u00e9siduelles sph\u00e9riques et gaussiennes, toute entorse s\u00e9v\u00e8re \u00e0 ces principes fragiliserait les conclusions tir\u00e9es.<\/p>\n<p>Ajustons de fa\u00e7on autonome le mod\u00e8le clinique optimis\u00e9 isol\u00e9 \u00e0 la section pr\u00e9c\u00e9dente :<\/p>\n<p><code>modele_optimal &lt;- lm(anxiete_gad7 ~ nevrosisme + coping_evitement + sommeil_rem, data = df_clinique)<\/code><\/p>\n<p>La validation formelle implique l&rsquo;ex\u00e9cution des tests diagnostiques canoniques :<\/p>\n<ul>\n<li><strong>Normalit\u00e9 des r\u00e9sidus :<\/strong> \u00c9valu\u00e9e conjointement par l&rsquo;inspection du diagramme quantile-quantile (QQ-plot) et par le test formel de Shapiro-Wilk :<br \/>\n    <code>shapiro.test(residuals(modele_optimal))<\/code><br \/>\n    Une valeur de <em>p<\/em> sup\u00e9rieure au seuil usuel de 0.05 confirme que la distribution des erreurs ne s&rsquo;\u00e9carte pas substantiellement d&rsquo;une loi normale univari\u00e9e.<\/li>\n<li><strong>Homosc\u00e9dasticit\u00e9 des perturbations :<\/strong> Contr\u00f4l\u00e9e par le test d&rsquo;homog\u00e9n\u00e9it\u00e9 des variances de Breusch-Pagan impl\u00e9ment\u00e9 dans le package <em>car<\/em> :<br \/>\n    <code>car::ncvTest(modele_optimal)<\/code><br \/>\n    L&rsquo;absence de significativit\u00e9 statistique valide l&rsquo;hypoth\u00e8se de variance constante de l&rsquo;erreur &sigma;<sup>2<\/sup> \u00e0 travers tout le spectre pr\u00e9dictif.<\/li>\n<li><strong>Lin\u00e9arit\u00e9 structurelle :<\/strong> V\u00e9rifi\u00e9e par la confrontation graphique des r\u00e9sidus ordinaires standardis\u00e9s aux valeurs ajust\u00e9es (<em>fitted values<\/em>), garantissant l&rsquo;absence de courbures polynomiales non mod\u00e9lis\u00e9es.<\/li>\n<\/ul>\n<h3>11.2 \u00c9valuation de la colin\u00e9arit\u00e9 r\u00e9siduelle (VIF)<\/h3>\n<p>L&rsquo;un des b\u00e9n\u00e9fices collat\u00e9raux majeurs de l&rsquo;\u00e9lagage parcimonieux orchestr\u00e9 par le crit\u00e8re de Mallows r\u00e9side dans l&rsquo;assainissement de la matrice de covariance des pr\u00e9dicteurs. Dans le mod\u00e8le complet initial, la pr\u00e9sence simultan\u00e9e de multiples indicateurs psychologiques corr\u00e9l\u00e9s engendre fr\u00e9quemment une instabilit\u00e9 num\u00e9rique d\u00e9l\u00e9t\u00e8re d\u00e9sign\u00e9e sous le nom de <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/web\/packages\/car\/index.html\">multicolin\u00e9arit\u00e9<\/a>.<\/p>\n<p>Cette instabilit\u00e9 se quantifie au moyen du Facteur d&rsquo;Inflation de la Variance (<em>Variance Inflation Factor<\/em> ou VIF). Pour un pr\u00e9dicteur <em>j<\/em> donn\u00e9, le VIF mesure \u00e0 quel point la variance de l&rsquo;estimateur &beta;_chapeau<sub>j<\/sub> est gonfl\u00e9e par sa corr\u00e9lation lin\u00e9aire avec l&rsquo;ensemble des autres r\u00e9gresseurs :<\/p>\n<p>VIF<sub>j<\/sub> = 1 \/ (1 &#8211; R<sub>j<\/sub><sup>2<\/sup>)<\/p>\n<p>Comparons les indices VIF calcul\u00e9s avant et apr\u00e8s s\u00e9lection \u00e0 l&rsquo;aide de la primitive <code>vif()<\/code> :<\/p>\n<p><code>vif_complet &lt;- car::vif(mod_complet_clinique)<\/code><br \/>\n<code>vif_optimal &lt;- car::vif(modele_optimal)<\/code><\/p>\n<p><code>print(\"VIF du modele complet :\")<\/code><br \/>\n<code>print(vif_complet)<\/code><br \/>\n<code>print(\"VIF du modele optimal apres arbitrage par Mallows :\")<\/code><br \/>\n<code>print(vif_optimal)<\/code><\/p>\n<p>Alors que le mod\u00e8le complet pr\u00e9sentait des facteurs VIF \u00e9lev\u00e9s sur les \u00e9chelles de coping et de r\u00e9activit\u00e9 en raison de redondances conceptuelles manifestes, le mod\u00e8le restreint affiche des indices VIF strictement inf\u00e9rieurs au seuil critique de 2.0. Cela confirme que le crit\u00e8re Cp a naturellement \u00e9limin\u00e9 les pr\u00e9dicteurs redondants, restituant une base vectorielle quasi-orthogonale garantissant des estimations de coefficients stables et ais\u00e9ment interpr\u00e9tables.<\/p>\n<h3>11.3 Validation crois\u00e9e et stabilit\u00e9 de l&rsquo;\u00e9chantillonnage<\/h3>\n<p>Bien que le Cp de Mallows constitue un estimateur asymptotique sans biais de l&rsquo;erreur quadratique de pr\u00e9diction, son comportement effectif hors-\u00e9chantillon gagne \u00e0 \u00eatre corrobor\u00e9 par une proc\u00e9dure de <strong>validation crois\u00e9e \u00e0 K blocs<\/strong> (<em>K-fold cross-validation<\/em>). Cette technique empirique consiste \u00e0 partitionner l&rsquo;\u00e9chantillon global en <em>K<\/em> sous-groupes ind\u00e9pendants (usuellement K = 10), \u00e0 entra\u00eener it\u00e9rativement le mod\u00e8le sur 9 blocs et \u00e0 mesurer son erreur quadratique moyenne de pr\u00e9diction (RMSEP) sur le dixi\u00e8me bloc non vu.<\/p>\n<p>Moyennant l&rsquo;utilisation des structures vectorielles de base de R, la validation crois\u00e9e s&rsquo;impl\u00e9mente en quelques instructions limpides :<\/p>\n<p><code>set.seed(42)<\/code><br \/>\n<code>k &lt;- 10<\/code><br \/>\n<code>folds &lt;- sample(rep(1:k, length.out = nrow(df_clinique)))<\/code><br \/>\n<code>cv_errors &lt;- numeric(k)<\/code><\/p>\n<p><code>for (i in 1:k) {<br \/>\n  train_data &lt;- df_clinique[folds != i, ]<br \/>\n  test_data  &lt;- df_clinique[folds == i, ]<\/p>\n<p>  fit_cv &lt;- lm(anxiete_gad7 ~ nevrosisme + coping_evitement + sommeil_rem, data = train_data)<br \/>\n  preds &lt;- predict(fit_cv, newdata = test_data)<br \/>\n  cv_errors[i] &lt;- mean((test_data$anxiete_gad7 - preds)^2)<br \/>\n}<\/code><\/p>\n<p><code>rmsep_cv &lt;- sqrt(mean(cv_errors))<\/code><br \/>\n<code>cat(\"Erreur quadratique moyenne de prediction hors-echantillon (RMSEP) :\", round(rmsep_cv, 3))<\/code><\/p>\n<p>Parall\u00e8lement, une d\u00e9marche de r\u00e9\u00e9chantillonnage non param\u00e9trique (<em>bootstrap<\/em>) peut \u00eatre engag\u00e9e pour \u00e9prouver la stabilit\u00e9 structurelle du choix de variables op\u00e9r\u00e9 par le Cp. En g\u00e9n\u00e9rant 1000 r\u00e9plicats de l&rsquo;\u00e9chantillon initial et en r\u00e9-ex\u00e9cutant l&rsquo;algorithme <code>ols_step_best_subset()<\/code> sur chaque tirage avec remise, le statisticien peut quantifier la fr\u00e9quence relative \u00e0 laquelle la combinaison ternaire {n\u00e9vrosisme, coping par \u00e9vitement, sommeil} est formellement isol\u00e9e. Une fr\u00e9quence de r\u00e9plication \u00e9lev\u00e9e (par exemple &gt; 85%) confirme la robustesse \u00e9cologique du mod\u00e8le face aux fluctuations al\u00e9atoires du recrutement clinique.<\/p>\n<h2>12. Bonnes pratiques, limites m\u00e9thodologiques et consid\u00e9rations \u00e9ditoriales en recherche quantitative<\/h2>\n<h3>12.1 Le biais induit par la s\u00e9lection de variables (post-selection inference)<\/h3>\n<p>L&rsquo;utilisation d&rsquo;algorithmes de s\u00e9lection de mod\u00e8les \u2014 qu&rsquo;ils reposent sur le crit\u00e8re de Mallows, sur l&rsquo;AIC ou sur des proc\u00e9dures pas \u00e0 pas \u2014 soul\u00e8ve une pr\u00e9occupation \u00e9pist\u00e9mologique majeure souvent pass\u00e9e sous silence dans les manuels appliqu\u00e9s : le probl\u00e8me de l&rsquo;<strong>inf\u00e9rence post-s\u00e9lection<\/strong> (<em>post-selection inference<\/em> ou biais d&rsquo;optimisation par les donn\u00e9es). Lorsque la m\u00eame matrice empirique est exploit\u00e9e simultan\u00e9ment pour identifier le sous-ensemble optimal de covariables et pour tester la significativit\u00e9 statistique des coefficients retenus, les postulats fondamentaux des tests d&rsquo;hypoth\u00e8se classiques sont viol\u00e9s.<\/p>\n<p>En s\u00e9lectionnant les r\u00e9gresseurs qui maximisent l&rsquo;ajustement empirique tout en minimisant RSS<sub>p<\/sub>, l&rsquo;algorithme capitalise in\u00e9vitablement sur le hasard d&rsquo;\u00e9chantillonnage. En cons\u00e9quence directe :<\/p>\n<ul>\n<li>Les valeurs de <em>p<\/em> issues du tableau <code>summary()<\/code> du mod\u00e8le final sont artificiellement trop optimistes (inflation du taux d&rsquo;erreur de type I, assimilable \u00e0 une forme involontaire de <em>p-hacking<\/em>).<\/li>\n<li>Les amplitudes des coefficients de r\u00e9gression standardis\u00e9s sont syst\u00e9matiquement surestim\u00e9es en valeur absolue (ph\u00e9nom\u00e8ne connu sous le nom d&rsquo;illusion du gagnant ou <em>winner&rsquo;s curse<\/em>).<\/li>\n<li>Les intervalles de confiance construits selon la loi de Student sont trop \u00e9troits, produisant un faux sentiment de pr\u00e9cision m\u00e9trique.<\/li>\n<\/ul>\n<p>Conform\u00e9ment aux recommandations \u00e9ditoriales de l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\/\">American Psychological Association (APA)<\/a> et de l&rsquo;<em>Open Science Framework<\/em>, les chercheurs doivent expliciter sans ambigu\u00eft\u00e9 dans leurs manuscrits si la sp\u00e9cification finale r\u00e9sulte d&rsquo;une d\u00e9marche purement exploratoire assist\u00e9e par le Cp de Mallows ou d&rsquo;un sch\u00e9ma d&rsquo;analyse confirmatoire pr\u00e9-enregistr\u00e9. Id\u00e9alement, une partition initiale de l&rsquo;\u00e9chantillon (<em>data splitting<\/em>) r\u00e9servant une fraction des donn\u00e9es (ex. 30%) \u00e0 la confirmation ind\u00e9pendante des inf\u00e9rences post-s\u00e9lection permet de restaurer la validit\u00e9 nominale des tests d&rsquo;hypoth\u00e8se.<\/p>\n<h3>12.2 Limites structurelles du crit\u00e8re de Mallows<\/h3>\n<p>En d\u00e9pit de son \u00e9l\u00e9gance matricielle et de son universalit\u00e9 dans le cadre de la r\u00e9gression lin\u00e9aire standard, le crit\u00e8re Cp de Mallows pr\u00e9sente des fronti\u00e8res structurelles d&rsquo;applicabilit\u00e9 que le mod\u00e9lisateur doit imp\u00e9rativement respecter :<\/p>\n<ul>\n<li><strong>Inapplicabilit\u00e9 aux Mod\u00e8les Lin\u00e9aires G\u00e9n\u00e9ralis\u00e9s (GLM) non gaussiens :<\/strong> La d\u00e9rivation alg\u00e9brique de l&rsquo;esp\u00e9rance math\u00e9matique E[Cp] = p repose explicitement sur l&rsquo;hypoth\u00e8se d&rsquo;additivit\u00e9 des r\u00e9sidus gaussiens homosc\u00e9dastiques. Dans le contexte de la r\u00e9gression logistique binaire, de la r\u00e9gression de Poisson pour donn\u00e9es de comptage ou des mod\u00e8les de survie \u00e0 risques proportionnels, le crit\u00e8re de Mallows n&rsquo;est pas transposable directement sous sa forme classique. Le chercheur devra lui substituer l&rsquo;AIC, le crit\u00e8re de d\u00e9viance d&rsquo;information (DIC) ou la d\u00e9viance r\u00e9siduelle g\u00e9n\u00e9ralis\u00e9e.<\/li>\n<li><strong>Vuln\u00e9rabilit\u00e9 face aux non-lin\u00e9arit\u00e9s complexes :<\/strong> Le crit\u00e8re de Mallows \u00e9value des sous-espaces vectoriels lin\u00e9aires. Si la relation v\u00e9ritable entre un r\u00e9gresseur et la variable r\u00e9ponse est hautement non lin\u00e9aire (ex. fonctions sinuso\u00efdales, seuils critiques, effets de saturation non lin\u00e9aires) et que celle-ci n&rsquo;a pas \u00e9t\u00e9 pr\u00e9alablement lin\u00e9aris\u00e9e via des splines ou des transformations polynomiales explicites dans le mod\u00e8le complet, le Cp \u00e9chouera \u00e0 quantifier ad\u00e9quatement le biais de sp\u00e9cification.<\/li>\n<li><strong>Sensibilit\u00e9 extr\u00eame \u00e0 la sp\u00e9cification du mod\u00e8le complet initial :<\/strong> Le crit\u00e8re Cp n&rsquo;\u00e9value jamais la validit\u00e9 intrins\u00e8que absolue d&rsquo;un mod\u00e8le dans l&rsquo;absolu ; il n&rsquo;\u00e9value qu&rsquo;une <em>distance relative<\/em> par rapport \u00e0 l&rsquo;\u00e9talon constitu\u00e9 par le mod\u00e8le complet. Si ce mod\u00e8le complet est lui-m\u00eame d\u00e9ficient, mal mesur\u00e9 ou amput\u00e9 de r\u00e9gresseurs structurels indispensables, l&rsquo;ensemble de l&rsquo;\u00e9difice comparatif est corrompu d\u00e8s sa base.<\/li>\n<li><strong>Concurrence des m\u00e9thodes de r\u00e9gularisation moderne :<\/strong> Face \u00e0 des matrices de tr\u00e8s haute dimension o\u00f9 le nombre de pr\u00e9dicteurs <em>P<\/em> approche ou d\u00e9passe la taille d&rsquo;\u00e9chantillon <em>n<\/em> (probl\u00e9matique P &ge; n courante en neuro-imagerie ou en g\u00e9nomique comportementale), la r\u00e9gression compl\u00e8te ne peut plus \u00eatre estim\u00e9e par moindres carr\u00e9s ordinaires en raison de la singularit\u00e9 de X<sup>T<\/sup>X. Le crit\u00e8re de Mallows devient math\u00e9matiquement inop\u00e9rant. Dans de telles configurations, les techniques de r\u00e9gularisation p\u00e9nalis\u00e9e continue \u2014 telles que la r\u00e9gression <em>Lasso<\/em> (L1), la r\u00e9gression <em>Ridge<\/em> (L2) ou l&rsquo;<em>Elastic Net<\/em> \u2014 s&rsquo;imposent comme les alternatives contemporaines incontournables.<\/li>\n<\/ul>\n<h3>12.3 Recommandations de r\u00e9daction pour publications acad\u00e9miques<\/h3>\n<p>La pr\u00e9sentation des r\u00e9sultats issus d&rsquo;une s\u00e9lection de mod\u00e8les par le crit\u00e8re Cp de Mallows au sein d&rsquo;une revue internationale \u00e0 comit\u00e9 de lecture requiert une transparence m\u00e9thodologique totale. L&rsquo;espace \u00e9ditorial ne doit pas se limiter \u00e0 mentionner laconiquement le mod\u00e8le s\u00e9lectionn\u00e9, mais doit documenter le processus d&rsquo;arbitrage ayant guid\u00e9 l&rsquo;investigateur.<\/p>\n<p>Les standards de communication scientifique imposent de respecter les consignes suivantes :<\/p>\n<ol>\n<li><strong>Sp\u00e9cification exhaustive du mod\u00e8le de r\u00e9f\u00e9rence :<\/strong> D\u00e9clarer formellement l&rsquo;ensemble des variables incluses dans le mod\u00e8le complet initial, leur \u00e9chelle de mesure, le nombre d&rsquo;observations effectives (apr\u00e8s traitement contr\u00f4l\u00e9 des donn\u00e9es manquantes), ainsi que la variance r\u00e9siduelle de r\u00e9f\u00e9rence s<sup>2<\/sup> (avec ses degr\u00e9s de libert\u00e9 r\u00e9siduels n &#8211; P).<\/li>\n<li><strong>Pr\u00e9sentation d&rsquo;un tableau synth\u00e9tique comparatif :<\/strong> Structurer les r\u00e9sultats sous forme d&rsquo;une table normalis\u00e9e alignant, pour les meilleurs sous-mod\u00e8les comp\u00e9titifs de chaque strate dimensionnelle : la liste des pr\u00e9dicteurs engag\u00e9s, le nombre total de param\u00e8tres estim\u00e9s <em>p<\/em> (incluant explicitement l&rsquo;intercept), le R<sup>2<\/sup>, le R<sup>2<\/sup> ajust\u00e9, la statistique Cp de Mallows, et l&rsquo;\u00e9cart absolu |Cp &#8211; p|.<\/li>\n<li><strong>Justification formelle du choix final :<\/strong> \u00c9noncer explicitement la r\u00e8gle heuristique mobilis\u00e9e pour d\u00e9partager les sp\u00e9cifications concurrentes (ex. s\u00e9lection du mod\u00e8le pr\u00e9sentant le Cp le plus bas sous la condition de proximit\u00e9 Cp &le; p). Argumenter substantiellement la plausibilit\u00e9 th\u00e9orique des variables \u00e9cart\u00e9es par l&rsquo;arbitrage math\u00e9matique.<\/li>\n<li><strong>Mise \u00e0 disposition des scripts et des donn\u00e9es :<\/strong> Conform\u00e9ment aux principes de la science ouverte (<em>Open Science<\/em>), d\u00e9poser les scripts de calcul R rigoureusement comment\u00e9s ainsi que les jeux de donn\u00e9es anonymis\u00e9s sur des d\u00e9p\u00f4ts institutionnels publics p\u00e9rennes tels que l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/osf.io\/\">Open Science Framework (OSF)<\/a>, Zenodo ou GitHub, garantissant la reproductibilit\u00e9 computationnelle imm\u00e9diate des analyses par les pairs.<\/li>\n<\/ol>\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\u2013723. <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>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7th ed.). American Psychological Association. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/li>\n<li>Burnham, K. P., &amp; Anderson, D. R. (2002). <em>Model selection and multimodel inference: A practical information-theoretic approach<\/em> (2nd ed.). Springer-Verlag. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/b97636\">https:\/\/doi.org\/10.1007\/b97636<\/a><\/li>\n<li>Fox, J., &amp; Weisberg, S. (2019). <em>An R companion to applied regression<\/em> (3rd ed.). SAGE Publications. <a rel=\"noopener noreferrer\" href=\"https:\/\/socialsciences.mcmaster.ca\/jfox\/Books\/Companion\/\">https:\/\/socialsciences.mcmaster.ca\/jfox\/Books\/Companion\/<\/a><\/li>\n<li>Furnival, G. M., &amp; Wilson, R. W. (1974). Regressions by leaps and bounds. <em>Technometrics<\/em>, 16(4), 499\u2013511. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00401706.1974.10489231\">https:\/\/doi.org\/10.1080\/00401706.1974.10489231<\/a><\/li>\n<li>Hebbali, A. (2020). <em>olsrr: Tools for building OLS regression models<\/em> (R package version 0.5.3). Comprehensive R Archive Network. <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=olsrr\">https:\/\/cran.r-project.org\/package=olsrr<\/a><\/li>\n<li>James, G., Witten, D., Hastie, T., &amp; Tibshirani, R. (2021). <em>An introduction to statistical learning: With applications in R<\/em> (2nd ed.). Springer. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-1-0716-1418-1\">https:\/\/doi.org\/10.1007\/978-1-0716-1418-1<\/a><\/li>\n<li>Lumley, T. (2020). <em>leaps: Regression subset selection<\/em> (R package version 3.1). Comprehensive R Archive Network. <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=leaps\">https:\/\/cran.r-project.org\/package=leaps<\/a><\/li>\n<li>Mallows, C. L. (1973). Some comments on Cp. <em>Technometrics<\/em>, 15(4), 661\u2013675. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00401706.1973.10489103\">https:\/\/doi.org\/10.1080\/00401706.1973.10489103<\/a><\/li>\n<li>R Core Team. (2023). <em>R: A language and environment for statistical computing<\/em>. R Foundation for Statistical Computing. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.R-project.org\/\">https:\/\/www.R-project.org\/<\/a><\/li>\n<li>Schwarz, G. (1978). Estimating the dimension of a model. <em>The Annals of Statistics<\/em>, 6(2), 461\u2013464. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aos\/1176344136\">https:\/\/doi.org\/10.1214\/aos\/1176344136<\/a><\/li>\n<li>Wickham, H., Averick, M., Bryan, J., Chang, W., McGowan, L. D., Fran\u00e7ois, R., Grolemund, G., Hayes, A., Henry, L., Hester, J., Kuhn, M., Pedersen, T. L., Miller, E., Bache, S. M., M\u00fcller, K., Ooms, J., Robinson, D., Seidel, D. P., Spinu, V., \u2026 Yutani, H. (2019). Welcome to the Tidyverse. <em>Journal of Open Source Software<\/em>, 4(43), Article 1686. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.21105\/joss.01686\">https:\/\/doi.org\/10.21105\/joss.01686<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet pour comprendre, calculer et interpr\u00e9ter le Cp de Mallows dans R pour la s\u00e9lection rigoureuse de mod\u00e8les de r\u00e9gression lin\u00e9aire.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2,1838],"tags":[2870,2881,1910,2871,2882],"class_list":["post-1286","psych_stat","type-psych_stat","status-publish","hentry","category-psychometrie","category-statistiques","tag-cp-de-mallows","tag-package-olsrr","tag-regression-lineaire","tag-selection-de-modele","tag-statistiques-r"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1286","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=1286"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1286\/revisions"}],"predecessor-version":[{"id":16215,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/1286\/revisions\/16215"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1286"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1286"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1286"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}