{"id":7504,"date":"2026-09-16T20:41:22","date_gmt":"2026-09-16T20:41:22","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/regression-binomiale-negative-vs-poisson-choix-modele\/"},"modified":"2026-09-16T20:41:25","modified_gmt":"2026-09-16T20:41:25","slug":"regression-binomiale-negative-vs-poisson-choix-modele","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=regression-binomiale-negative-vs-poisson-choix-modele","title":{"rendered":"R\u00e9gression binomiale n\u00e9gative vs Poisson : comment choisir un mod\u00e8le de r\u00e9gression"},"content":{"rendered":"<p>Dans le domaine des sciences comportementales, de la psychologie clinique, de la psychiatrie quantitative et de la recherche biom\u00e9dicale, les chercheurs sont constamment confront\u00e9s \u00e0 des variables qui mesurent le nombre de fois o\u00f9 un \u00e9v\u00e9nement sp\u00e9cifique se manifeste au cours d&rsquo;une p\u00e9riode temporelle ou dans un espace d\u00e9limit\u00e9. Qu&rsquo;il s&rsquo;agisse de comptabiliser les \u00e9pisodes d&rsquo;automutilation chez des adolescents souffrant de troubles de la personnalit\u00e9 limite, les r\u00e9cidives d&rsquo;abus de substances chez des patients en cours de sevrage, les crises de panique hebdomadaires ou le nombre d&rsquo;erreurs d&rsquo;inattention dans une t\u00e2che de neuropsychologie informatis\u00e9e, ces mesures partagent une structure math\u00e9matique commune : elles se composent exclusivement d&rsquo;entiers naturels positifs ou nuls. Malgr\u00e9 cette caract\u00e9ristique fondamentale, une proportion alarmante d&rsquo;\u00e9tudes publi\u00e9es continue de recourir \u00e0 des mod\u00e8les de r\u00e9gression lin\u00e9aire par les moindres carr\u00e9s ordinaires (MCO), au prix d&rsquo;artifices m\u00e9thodologiques tels que la transformation logarithmique arbitraire des scores.<\/p>\n<p>L&rsquo;application inappropri\u00e9e de mod\u00e8les lin\u00e9aires continus et gaussiens \u00e0 des distributions discr\u00e8tes et asym\u00e9triques engendre des biais m\u00e9thodologiques majeurs. Ces failles analytiques se traduisent par des estimations d&rsquo;intervalles de confiance invalides, des pr\u00e9dictions absurdes de valeurs n\u00e9gatives pour des ph\u00e9nom\u00e8nes pourtant non n\u00e9gatifs, et une distorsion syst\u00e9matique des erreurs-types. Afin de pr\u00e9server l&rsquo;int\u00e9grit\u00e9 de la mod\u00e9lisation statistique, l&rsquo;adoption des <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Generalized_linear_model\">mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s (GLM)<\/a> s&rsquo;impose comme une n\u00e9cessit\u00e9 m\u00e9thodologique incontournable. Au c\u0153ur de ce paradigme, deux cadres de mod\u00e9lisation param\u00e9trique se disputent la primaut\u00e9 : la r\u00e9gression de Poisson et la r\u00e9gression binomiale n\u00e9gative. Bien que la premi\u00e8re soit souvent per\u00e7ue comme le point de d\u00e9part canonique pour les donn\u00e9es de comptage, la seconde constitue fr\u00e9quemment l&rsquo;ajustement empirique le plus r\u00e9aliste face \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des comportements humains.<\/p>\n<p>L&rsquo;arbitrage formel entre ces deux approches ne rel\u00e8ve pas d&rsquo;une simple pr\u00e9f\u00e9rence computationnelle, mais engage la validit\u00e9 m\u00eame des conclusions scientifiques tir\u00e9es d&rsquo;une \u00e9tude. Choisir \u00e0 tort le mod\u00e8le de Poisson en pr\u00e9sence de surdispersion expose l&rsquo;analyste \u00e0 un effondrement dramatique de ses erreurs-types, multipliant artificiellement les d\u00e9clarations de significativit\u00e9 statistique et alimentant directement la crise de reproductibilit\u00e9 qui frappe les sciences empiriques. \u00c0 l&rsquo;inverse, l&#8217;emploi syst\u00e9matique et non justifi\u00e9 de mod\u00e8les plus complexes peut engendrer une instabilit\u00e9 param\u00e9trique inutile lorsque la parcimonie est de mise. Ce guide exhaustif a pour vocation d&rsquo;analyser en profondeur les fondements math\u00e9matiques, les hypoth\u00e8ses stochastiques, les m\u00e9thodes diagnostiques et les strat\u00e9gies de s\u00e9lection r\u00e9gissant la confrontation entre r\u00e9gression de Poisson et r\u00e9gression binomiale n\u00e9gative, offrant ainsi aux chercheurs un cadre rigoureux pour \u00e9tayer leurs choix analytiques.<\/p>\n<h2>1. 1. Introduction aux donn\u00e9es de comptage en sciences comportementales et psychologiques<\/h2>\n<h3>1.1 D\u00e9finition et propri\u00e9t\u00e9s math\u00e9matiques des donn\u00e9es de comptage<\/h3>\n<p>Les donn\u00e9es de comptage correspondent \u00e0 des variables al\u00e9atoires quantitatives discr\u00e8tes dont le support math\u00e9matique est restreint \u00e0 l&rsquo;ensemble des entiers naturels, not\u00e9 &#x2115; = {0, 1, 2, 3, &#8230;}. Contrairement aux variables continues, une variable de comptage ne peut adopter de valeurs fractionnaires ou d\u00e9cimales : elle d\u00e9nombre des unit\u00e9s d&rsquo;occurrence distinctes et indivisibles d&rsquo;un ph\u00e9nom\u00e8ne d\u00e9fini. De plus, ces donn\u00e9es poss\u00e8dent une fronti\u00e8re structurelle stricte en z\u00e9ro, ce qui interdit rigoureusement l&rsquo;existence de valeurs n\u00e9gatives. Dans la quasi-totalit\u00e9 des contextes comportementaux et psychologiques, les distributions empiriques de ces variables pr\u00e9sentent une forte asym\u00e9trie positive (dextrogire), caract\u00e9ris\u00e9e par une masse probabiliste imposante concentr\u00e9e sur les petites valeurs enti\u00e8res \u2014 et particuli\u00e8rement sur le z\u00e9ro \u2014 ainsi qu&rsquo;une longue tra\u00eene s&rsquo;\u00e9tendant vers des valeurs \u00e9lev\u00e9es plus rares.<\/p>\n<p>Cette structure distributionnelle viole intrins\u00e8quement les deux postulats fondamentaux sur lesquels repose la r\u00e9gression lin\u00e9aire classique par les moindres carr\u00e9s ordinaires (OLS) : la normalit\u00e9 des r\u00e9sidus et l&rsquo;homosc\u00e9dasticit\u00e9, c&rsquo;est-\u00e0-dire la constance de la variance de l&rsquo;erreur conditionnellement aux pr\u00e9dicteurs. Dans les processus de comptage r\u00e9els, la variance empirique est intrins\u00e8quement d\u00e9pendante de la moyenne conditionnelle : plus le nombre moyen d&rsquo;\u00e9v\u00e9nements pr\u00e9dits est grand, plus l&rsquo;incertitude et l&rsquo;\u00e9talement des observations augmentent. Ignorer cette liaison moyenne-variance conduit in\u00e9vitablement l&rsquo;estimateur OLS \u00e0 produire des erreurs-types h\u00e9t\u00e9rosc\u00e9dastiques et non fiables, ruinant toute tentative d&rsquo;inf\u00e9rence valide bas\u00e9e sur les statistiques de Student ou de Fisher.<\/p>\n<p>Sur le plan pr\u00e9dictif, l&rsquo;application aveugle d&rsquo;un mod\u00e8le lin\u00e9aire standard \u00e0 une variable de comptage aboutit fr\u00e9quemment \u00e0 des aberrations manifestes, telles que la g\u00e9n\u00e9ration de valeurs pr\u00e9dites n\u00e9gatives pour des sous-groupes d&rsquo;individus pr\u00e9sentant de faibles scores sur les covariables explicatives. En recherche clinique, il est impensable d&rsquo;estimer qu&rsquo;un patient ayant suivi une th\u00e9rapie comportementale r\u00e9alisera, en moyenne, \u00ab -1,4 crise d&rsquo;angoisse \u00bb par mois. Par ailleurs, des variables critiques comme le nombre de tentatives de suicide, la fr\u00e9quence des r\u00e9cidives chez les usagers d&rsquo;opiac\u00e9s, ou le d\u00e9nombrement des lapsus attentionnels lors d&rsquo;un protocole cognitif sous privation de sommeil, exhibent une nature discr\u00e8te ind\u00e9niable qui interdit l&rsquo;assimilation \u00e0 un continuum gaussien sans introduire de biais substantiels dans l&rsquo;estimation de la pente des effets.<\/p>\n<h3>1.2 Le r\u00f4le fondamental des mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s (GLM)<\/h3>\n<p>Introduits dans les travaux fondateurs de <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2344614\">Nelder et Wedderburn (1972)<\/a>, les mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s ont profond\u00e9ment r\u00e9volutionn\u00e9 l&rsquo;analyse des variables d\u00e9pendantes non gaussiennes en proposant une architecture param\u00e9trique unifi\u00e9e. Un GLM repose formellement sur une structure tripartite : une composante al\u00e9atoire, un pr\u00e9dicteur lin\u00e9aire et une fonction de lien. La composante al\u00e9atoire sp\u00e9cifie que la variable d\u00e9pendante <em>Y<\/em>, conditionnellement aux pr\u00e9dicteurs, suit une distribution appartenant \u00e0 la famille exponentielle de dispersions (englobant les lois gaussienne, de Poisson, binomiale, gamma, et sous certaines formulations, binomiale n\u00e9gative). Le pr\u00e9dicteur lin\u00e9aire, not\u00e9 &eta;<sub>i<\/sub>, agr\u00e8ge lin\u00e9airement l&rsquo;ensemble des covariables associ\u00e9es \u00e0 leurs coefficients respectifs selon la formulation vectorielle standard :<\/p>\n<p>&eta;<sub>i<\/sub> = &beta;<sub>0<\/sub> + &beta;<sub>1<\/sub>X<sub>1i<\/sub> + &beta;<sub>2<\/sub>X<sub>2i<\/sub> + &#8230; + &beta;<sub>k<\/sub>X<sub>ki<\/sub><\/p>\n<p>La composante pivotale qui articule ce pr\u00e9dicteur lin\u00e9aire avec l&rsquo;esp\u00e9rance conditionnelle &mu;<sub>i<\/sub> = E(Y<sub>i<\/sub>|X<sub>i<\/sub>) de la variable de comptage est la fonction de lien, not\u00e9e <em>g(&middot;)<\/em>, telle que &eta;<sub>i<\/sub> = <em>g(&mu;<sub>i<\/sub>)<\/em>. Pour les mod\u00e8les de comptage, la fonction de lien canonique est la fonction logarithmique naturelle, <em>g(&mu;)<\/em> = ln(&mu;). D\u00e8s lors, le mod\u00e8le s&rsquo;exprime sous la forme ln(&mu;<sub>i<\/sub>) = X<sub>i<\/sub>&beta;, ce qui \u00e9quivaut, par inversion de la fonction de lien, \u00e0 une relation multiplicative non lin\u00e9aire pour la moyenne conditionnelle :<\/p>\n<p>&mu;<sub>i<\/sub> = exp(X<sub>i<\/sub>&beta;) = exp(&beta;<sub>0<\/sub>) &times; exp(&beta;<sub>1<\/sub>X<sub>1i<\/sub>) &times; &#8230; &times; exp(&beta;<sub>k<\/sub>X<sub>ki<\/sub>)<\/p>\n<p>Cette formulation garantit rigoureusement que l&rsquo;esp\u00e9rance math\u00e9matique &mu;<sub>i<\/sub> demeure strictement positive pour toute combinaison de covariables, \u00e9liminant ainsi toute possibilit\u00e9 de pr\u00e9diction n\u00e9gative absurde. De surcro\u00eet, le recours \u00e0 un GLM log-lin\u00e9aire surpasse cat\u00e9goriquement la pratique artisanale et d\u00e9pass\u00e9e consistant \u00e0 appliquer une transformation empirique directe de type log(Y + 1) sur la variable brute avant d&rsquo;ajuster un mod\u00e8le gaussien. Les travaux m\u00e9thodologiques de <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2041-210X.2010.00021.x\">O&rsquo;Hara et Kotze (2010)<\/a> ont d\u00e9montr\u00e9 de mani\u00e8re indiscutable que la transformation log(Y + 1) induit des biais d&rsquo;estimation massifs, distord les relations fonctionnelles sous-jacentes, et engendre une perte s\u00e9v\u00e8re de puissance statistique, en particulier lorsque la moyenne conditionnelle est faible ou que la proportion de z\u00e9ros est pr\u00e9pond\u00e9rante. Le GLM pr\u00e9serve quant \u00e0 lui l&rsquo;int\u00e9grit\u00e9 de l&rsquo;espace d&rsquo;\u00e9chantillonnage discret tout en mod\u00e9lisant ad\u00e9quatement la structure h\u00e9t\u00e9rosc\u00e9dastique des r\u00e9sidus.<\/p>\n<h3>1.3 Probl\u00e9matique centrale du choix de mod\u00e9lisation<\/h3>\n<p>L&rsquo;acc\u00e8s aux mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s ouvre une interrogation d\u00e9cisive pour le statisticien et le chercheur en sciences humaines : quel mod\u00e8le param\u00e9trique sp\u00e9cifique convient-il de privil\u00e9gier pour d\u00e9crire la composante al\u00e9atoire ? L&rsquo;arbitrage s&rsquo;articule principalement entre la r\u00e9gression de Poisson, mod\u00e8le de r\u00e9f\u00e9rence par excellence pour les d\u00e9nombrements, et la r\u00e9gression binomiale n\u00e9gative, sa g\u00e9n\u00e9ralisation flexible con\u00e7ue pour rel\u00e2cher les contraintes stochastiques strictes. La r\u00e9gression de Poisson se distingue par son \u00e9l\u00e9gante parcimonie : elle postule que la variance conditionnelle est rigoureusement identique \u00e0 la moyenne conditionnelle, un ph\u00e9nom\u00e8ne math\u00e9matique d\u00e9sign\u00e9 sous le terme d&rsquo;\u00e9quidispersion.<\/p>\n<p>Cependant, dans la r\u00e9alit\u00e9 complexe des donn\u00e9es observationnelles et exp\u00e9rimentales en psychologie, cette \u00e9quivalence th\u00e9orique se v\u00e9rifie presque jamais. L&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 biologique, cognitive et sociologique non captur\u00e9e par les variables explicatives g\u00e9n\u00e8re une dispersion empirique bien plus large que celle tol\u00e9r\u00e9e par le cadre de Poisson, un ph\u00e9nom\u00e8ne omnipr\u00e9sent qualifi\u00e9 de surdispersion. Lorsque l&rsquo;analyste persiste \u00e0 employer un mod\u00e8le de Poisson face \u00e0 des donn\u00e9es surdispers\u00e9es, les estimateurs des coefficients de r\u00e9gression demeurent convergents, mais leurs erreurs-types associ\u00e9es sont s\u00e9v\u00e8rement sous-estim\u00e9es. Cette contraction artificielle de l&rsquo;incertitude entra\u00eene une explosion des erreurs de premi\u00e8re esp\u00e8ce (taux de faux positifs), conduisant les chercheurs \u00e0 conclure abusivement \u00e0 l&rsquo;efficacit\u00e9 d&rsquo;une th\u00e9rapie ou \u00e0 la significativit\u00e9 de corr\u00e9lations psychopathologiques.<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;utilisation automatique et non contextualis\u00e9e de la r\u00e9gression binomiale n\u00e9gative soul\u00e8ve des questions de robustesse et de sur-param\u00e9trisation dans les \u00e9chantillons de taille restreinte. Estimer conjointement les pentes de r\u00e9gression et un param\u00e8tre de surdispersion requiert des algorithmes d&rsquo;optimisation num\u00e9rique it\u00e9ratifs qui peuvent peiner \u00e0 converger si la taille d&rsquo;\u00e9chantillon est faible ou si les donn\u00e9es ne contiennent pas une variance exc\u00e9dentaire suffisante. Le choix raisonn\u00e9 entre ces deux approches ne peut \u00eatre r\u00e9duit \u00e0 une simple habitude logicielle : il conditionne la validit\u00e9 inf\u00e9rentielle, la r\u00e9plicabilit\u00e9 des r\u00e9sultats empiriques et la port\u00e9e clinique des inf\u00e9rences formul\u00e9es.<\/p>\n<h2>2. 2. Fondements th\u00e9oriques du mod\u00e8le de r\u00e9gression de Poisson<\/h2>\n<h3>2.1 Distribution de Poisson et formulation stochastique<\/h3>\n<p>La distribution univari\u00e9e de Poisson constitue le socle probabiliste de la mod\u00e9lisation des \u00e9v\u00e9nements discrets se d\u00e9roulant dans un continuum spatio-temporel fixe. Soit <em>Y<\/em> une variable al\u00e9atoire d\u00e9signant le nombre d&rsquo;occurrences d&rsquo;un \u00e9v\u00e9nement au cours d&rsquo;une fen\u00eatre d&rsquo;observation standardis\u00e9e. On dit que <em>Y<\/em> suit une loi de Poisson de param\u00e8tre &lambda; &gt; 0, not\u00e9e <em>Y ~ P(&lambda;)<\/em>, si sa fonction de masse de probabilit\u00e9 discr\u00e8te s&rsquo;exprime, pour toute valeur enti\u00e8re <em>y<\/em> &isin; &#x2115;, par l&rsquo;\u00e9quation suivante :<\/p>\n<p>P(Y = y | &lambda;) = [e<sup>-&lambda;<\/sup> &times; &lambda;<sup>y<\/sup>] \/ y!<\/p>\n<p>o\u00f9 <em>e<\/em> repr\u00e9sente la constante n\u00e9p\u00e9rienne et <em>y!<\/em> d\u00e9signe la factorielle de l&rsquo;entier <em>y<\/em>. L&rsquo;\u00e9mergence d&rsquo;un processus de Poisson pur repose sur trois postulats fondamentaux : l&rsquo;homog\u00e9n\u00e9it\u00e9 temporelle (le taux d&rsquo;occurrence instantan\u00e9 &lambda; est constant durant toute la dur\u00e9e d&rsquo;observation), la non-simultan\u00e9it\u00e9 des \u00e9v\u00e9nements (la probabilit\u00e9 que deux \u00e9v\u00e9nements distincts se produisent dans un intervalle de temps infinit\u00e9simal est n\u00e9gligeable), et l&rsquo;ind\u00e9pendance stochastique (le nombre d&rsquo;\u00e9v\u00e9nements observ\u00e9s dans un intervalle donn\u00e9 est totalement ind\u00e9pendant du nombre d&rsquo;\u00e9v\u00e9nements survenus dans tout intervalle disjoint ant\u00e9rieur).<\/p>\n<p>Dans le cadre de la r\u00e9gression de Poisson, l&rsquo;intensit\u00e9 moyenne d&rsquo;occurrence n&rsquo;est plus consid\u00e9r\u00e9e comme constante pour l&rsquo;ensemble des individus, mais varie en fonction d&rsquo;un vecteur de pr\u00e9dicteurs <em>X<sub>i<\/sub><\/em> selon la sp\u00e9cification log-lin\u00e9aire &lambda;<sub>i<\/sub> = exp(X<sub>i<\/sub>&beta;). Pour un \u00e9chantillon de <em>n<\/em> observations ind\u00e9pendantes, la fonction de vraisemblance conjointe <em>L(&beta;)<\/em> s&rsquo;obtient par le produit des probabilit\u00e9s individuelles :<\/p>\n<p>L(&beta;) = &prod;<sub>i=1<\/sub><sup>n<\/sup> [exp(-exp(X<sub>i<\/sub>&beta;)) &times; (exp(X<sub>i<\/sub>&beta;))<sup>y<sub>i<\/sub><\/sup>] \/ y<sub>i<\/sub>!<\/p>\n<p>En prenant le logarithme naturel de cette fonction, nous obtenons la log-vraisemblance standard :<\/p>\n<p>ln L(&beta;) = &sum;<sub>i=1<\/sub><sup>n<\/sup> [y<sub>i<\/sub> (X<sub>i<\/sub>&beta;) &#8211; exp(X<sub>i<\/sub>&beta;) &#8211; ln(y<sub>i<\/sub>!)]<\/p>\n<p>L&rsquo;estimation des coefficients vectoriels &beta; s&rsquo;effectue par la m\u00e9thode du <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Maximum_likelihood_estimation\">maximum de vraisemblance (MLE)<\/a>. En d\u00e9rivant la log-vraisemblance par rapport au vecteur &beta;, on obtient le vecteur gradient, ou \u00e9quations de score :<\/p>\n<p>&part;(ln L) \/ &part;&beta; = &sum;<sub>i=1<\/sub><sup>n<\/sup> [y<sub>i<\/sub> &#8211; exp(X<sub>i<\/sub>&beta;)] X<sub>i<\/sub> = 0<\/p>\n<p>Ce syst\u00e8me d&rsquo;\u00e9quations transcendantes ne poss\u00e8de pas de solution analytique ferm\u00e9e dans le cas g\u00e9n\u00e9ral. Sa r\u00e9solution requiert l&#8217;emploi d&rsquo;algorithmes d&rsquo;optimisation num\u00e9rique tels que l&rsquo;algorithme des moindres carr\u00e9s it\u00e9rativement repond\u00e9r\u00e9s (IRLS) ou la m\u00e9thode de Newton-Raphson. Sous les conditions standards de r\u00e9gularit\u00e9 asymptotique, l&rsquo;estimateur r\u00e9sultant &beta;&#770; est convergent, asymptotiquement normal et asymptotiquement efficace, avec une matrice de variance-covariance asymptotique donn\u00e9e par l&rsquo;inverse de la matrice d&rsquo;information de Fisher observ\u00e9e.<\/p>\n<h3>2.2 Interpr\u00e9tation des coefficients et ratios de taux d&rsquo;incidence (IRR)<\/h3>\n<p>L&rsquo;interpr\u00e9tation des coefficients estim\u00e9s &beta;&#770; dans un mod\u00e8le log-lin\u00e9aire diff\u00e8re radicalement de celle des mod\u00e8les lin\u00e9aires gaussiens en raison de la nature multiplicative impos\u00e9e par la fonction de lien. Soit un mod\u00e8le univari\u00e9 comportant une seule covariable continue <em>X<\/em> : ln(E(Y|X)) = &beta;<sub>0<\/sub> + &beta;<sub>1<\/sub>X. Si nous augmentons la valeur de la variable explicative d&rsquo;une unit\u00e9, la nouvelle esp\u00e9rance conditionnelle devient E(Y|X + 1) = exp(&beta;<sub>0<\/sub> + &beta;<sub>1<\/sub>(X + 1)). En examinant le rapport de ces deux esp\u00e9rances math\u00e9matiques, nous constatons :<\/p>\n<p>E(Y|X + 1) \/ E(Y|X) = exp(&beta;<sub>0<\/sub> + &beta;<sub>1<\/sub>X + &beta;<sub>1<\/sub>) \/ exp(&beta;<sub>0<\/sub> + &beta;<sub>1<\/sub>X) = exp(&beta;<sub>1<\/sub>)<\/p>\n<p>Cette quantit\u00e9 exponenti\u00e9e, not\u00e9e exp(&beta;), est universellement appel\u00e9e le <strong>Ratio de Taux d&rsquo;Incidence<\/strong> (<em>Incidence Rate Ratio<\/em>, ou IRR). L&rsquo;IRR exprime la modification proportionnelle et multiplicative du taux attendu d&rsquo;occurrence de l&rsquo;\u00e9v\u00e9nement induite par une augmentation d&rsquo;une unit\u00e9 de la variable pr\u00e9dictive correspondante, toutes les autres covariables du mod\u00e8le \u00e9tant maintenues rigoureusement constantes.<\/p>\n<p>Dans la pratique clinique, la conversion des coefficients logarithmiques bruts en IRR simplifie consid\u00e9rablement la restitution des conclusions empiriques :<\/p>\n<ul>\n<li>Si exp(&beta;<sub>j<\/sub>) = 1 (ou &beta;<sub>j<\/sub> = 0), la covariable <em>X<sub>j<\/sub><\/em> n&rsquo;exerce aucune influence sur le taux moyen de survenue des \u00e9v\u00e9nements.<\/li>\n<li>Si exp(&beta;<sub>j<\/sub>) &gt; 1 (ou &beta;<sub>j<\/sub> &gt; 0), chaque accroissement d&rsquo;une unit\u00e9 de <em>X<sub>j<\/sub><\/em> multiplie le taux moyen d&rsquo;incidents par un facteur exp(&beta;<sub>j<\/sub>). Par exemple, un IRR de 1,35 se traduit cliniquement par une augmentation de 35 % du nombre attendu d&rsquo;\u00e9v\u00e9nements : [exp(&beta;) &#8211; 1] &times; 100 %.<\/li>\n<li>Si exp(&beta;<sub>j<\/sub>) &lt; 1 (ou &beta;<sub>j<\/sub> &lt; 0), la covariable est associ\u00e9e \u00e0 une r\u00e9duction du taux moyen d&rsquo;incidents. Un IRR de 0,72 indique ainsi une diminution relative de 28 % du taux attendu de l&rsquo;\u00e9v\u00e9nement clinique consid\u00e9r\u00e9.<\/li>\n<\/ul>\n<p>L&rsquo;estimation des intervalles de confiance \u00e0 95 % entourant ces ratios de taux d&rsquo;incidence s&rsquo;effectue usuellement en appliquant la transformation exponentielle directement aux bornes de l&rsquo;intervalle de confiance de Wald calcul\u00e9 sur le coefficient additif brut : [exp(&beta;&#770;<sub>j<\/sub> &#8211; 1.96 &times; SE(&beta;&#770;<sub>j<\/sub>)) ; exp(&beta;&#770;<sub>j<\/sub> + 1.96 &times; SE(&beta;&#770;<sub>j<\/sub>))]. Alternativement, pour les param\u00e8tres transform\u00e9s, l&rsquo;usage de la <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Delta_method\">m\u00e9thode delta<\/a> permet d&rsquo;approximer la variance asymptotique de exp(&beta;&#770;), garantissant une pr\u00e9cision inf\u00e9rentielle ad\u00e9quate m\u00eame sous de l\u00e9g\u00e8res asym\u00e9tries de l&rsquo;\u00e9chantillonnage fini.<\/p>\n<h3>2.3 Inclusion des termes d&rsquo;exposition statistique (Offsets)<\/h3>\n<p>L&rsquo;une des contraintes analytiques fr\u00e9quentes en psychologie exp\u00e9rimentale et \u00e9pid\u00e9miologique r\u00e9side dans le fait que les unit\u00e9s observationnelles ne sont pas toujours examin\u00e9es durant des intervalles de temps identiques, ou ne sont pas expos\u00e9es \u00e0 un volume \u00e9gal d&rsquo;opportunit\u00e9s de g\u00e9n\u00e9rer l&rsquo;\u00e9v\u00e9nement \u00e9tudi\u00e9. Par exemple, mesurer le nombre d&rsquo;actes d&rsquo;agressivit\u00e9 chez des patients hospitalis\u00e9s en psychiatrie n&rsquo;a de sens que si l&rsquo;on prend en consid\u00e9ration le nombre exact de jours d&rsquo;hospitalisation de chaque patient. Un patient hospitalis\u00e9 180 jours pr\u00e9sentant 10 actes agressifs pr\u00e9sente une intensit\u00e9 comportementale nettement inf\u00e9rieure \u00e0 celle d&rsquo;un patient cumulant 8 actes au cours d&rsquo;un s\u00e9jour de seulement 14 jours.<\/p>\n<p>Pour int\u00e9grer rigoureusement cette variabilit\u00e9 de l&rsquo;exposition, le mod\u00e8le doit mod\u00e9liser un taux par unit\u00e9 de temps ou d&rsquo;effort, soit &mu;<sub>i<\/sub> \/ t<sub>i<\/sub>, o\u00f9 <em>t<sub>i<\/sub><\/em> d\u00e9signe la mesure d&rsquo;exposition quantitative propre \u00e0 l&rsquo;individu <em>i<\/em>. En formalisant cette relation par le lien logarithmique, nous obtenons :<\/p>\n<p>ln(&mu;<sub>i<\/sub> \/ t<sub>i<\/sub>) = X<sub>i<\/sub>&beta; &amp;DoubleRightArrow; ln(&mu;<sub>i<\/sub>) &#8211; ln(t<sub>i<\/sub>) = X<sub>i<\/sub>&beta; &amp;DoubleRightArrow; ln(&mu;<sub>i<\/sub>) = ln(t<sub>i<\/sub>) + X<sub>i<\/sub>&beta;<\/p>\n<p>Le terme ln(t<sub>i<\/sub>) est d\u00e9sign\u00e9 sous le terme technique d&rsquo;<strong>offset<\/strong> statistique. Dans l&rsquo;architecture du GLM, l&rsquo;offset correspond \u00e0 une covariable quantitative dont le coefficient de r\u00e9gression est structurellement fix\u00e9 \u00e0 la valeur 1, et n&rsquo;est donc pas estim\u00e9 par l&rsquo;algorithme d&rsquo;optimisation. L&rsquo;omission d&rsquo;un offset requis alt\u00e8re substantiellement la signification th\u00e9orique des estimateurs : elle confond la variation r\u00e9elle du taux comportemental intrins\u00e8que avec la simple variabilit\u00e9 des dur\u00e9es de mesure ou du nombre d&rsquo;items administr\u00e9s aux participants, produisant des artefacts statistiques majeurs.<\/p>\n<h2>3. 3. L&rsquo;hypoth\u00e8se fondamentale d&rsquo;\u00e9quidispersion et ses limites pratiques<\/h2>\n<h3>3.1 Le principe formel de l&rsquo;\u00e9quidispersion<\/h3>\n<p>La distribution de Poisson poss\u00e8de une singularit\u00e9 math\u00e9matique remarquable qui constitue \u00e0 la fois sa force conceptuelle et sa plus grande vuln\u00e9rabilit\u00e9 empirique : elle est enti\u00e8rement caract\u00e9ris\u00e9e par un param\u00e8tre unique, &lambda;. De cette architecture \u00e0 param\u00e8tre unique d\u00e9coule directement la propri\u00e9t\u00e9 fondamentale d&rsquo;<strong>\u00e9quidispersion<\/strong>, stipulant l&rsquo;\u00e9galit\u00e9 absolue entre l&rsquo;esp\u00e9rance math\u00e9matique du processus et sa variance :<\/p>\n<p>E(Y) = Var(Y) = &lambda;<\/p>\n<p>Transpos\u00e9e au contexte de la r\u00e9gression de Poisson, l&rsquo;\u00e9quidispersion impose que pour n&rsquo;importe quelle configuration donn\u00e9e de covariables explicatives <em>X<\/em>, la variance conditionnelle des r\u00e9ponses individuelles soit strictement \u00e9gale \u00e0 la moyenne conditionnelle pr\u00e9dite par le mod\u00e8le :<\/p>\n<p>Var(Y<sub>i<\/sub> | X<sub>i<\/sub>) = E(Y<sub>i<\/sub> | X<sub>i<\/sub>) = &mu;<sub>i<\/sub> = exp(X<sub>i<\/sub>&beta;)<\/p>\n<p>Cette hypoth\u00e8se implique une rigidit\u00e9 structurelle extr\u00eame. Elle postule que le m\u00e9canisme stochastique r\u00e9gissant la fr\u00e9quence moyenne des \u00e9v\u00e9nements gouverne \u00e9galement, de fa\u00e7on d\u00e9terministe et invariable, la dispersion de ces \u00e9v\u00e9nements autour de cette moyenne. Si cette propri\u00e9t\u00e9 offre une \u00e9l\u00e9gance math\u00e9matique ind\u00e9niable et simplifie le calcul de la vraisemblance, elle s&rsquo;accorde tr\u00e8s mal avec la dynamique des comportements humains. En sciences sociales et comportementales, il est exceptionnel que la variabilit\u00e9 r\u00e9siduelle observ\u00e9e se conforme \u00e0 cette contrainte th\u00e9orique \u00e9troite.<\/p>\n<h3>3.2 Facteurs induisant la rupture de l&rsquo;\u00e9quidispersion<\/h3>\n<p>L&rsquo;effondrement syst\u00e9matique de l&rsquo;\u00e9quidispersion dans les plans de recherche exp\u00e9rimentaux et observationnels d\u00e9coule de plusieurs m\u00e9canismes sous-jacents bien document\u00e9s en m\u00e9thodologie quantitative. Le premier r\u00e9side dans l&rsquo;existence d&rsquo;une <strong>h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 individuelle non observ\u00e9e<\/strong>. Dans toute \u00e9tude psychologique, des facteurs latents majeurs \u2014 variations g\u00e9n\u00e9tiques, susceptibilit\u00e9s neurobiologiques, traits de personnalit\u00e9 ou contextes familiaux non mesur\u00e9s \u2014 introduisent des disparit\u00e9s substantielles de propension de base entre les individus. D\u00e8s lors que cette variabilit\u00e9 latente n&rsquo;est pas explicitement mod\u00e9lis\u00e9e par des pr\u00e9dicteurs dans la matrice <em>X<\/em>, elle se diffuse dans le terme d&rsquo;erreur r\u00e9siduel, augmentant la variabilit\u00e9 bien au-del\u00e0 de la valeur pr\u00e9dite par le mod\u00e8le conditionnel de Poisson.<\/p>\n<p>Le second m\u00e9canisme rel\u00e8ve de la <strong>violation du postulat d&rsquo;ind\u00e9pendance des \u00e9v\u00e9nements<\/strong>, fr\u00e9quemment d\u00e9sign\u00e9e sous l&rsquo;appellation de contagion temporelle ou de renforcement comportemental. Le mod\u00e8le de Poisson exige que la survenue d&rsquo;un \u00e9v\u00e9nement n&rsquo;alt\u00e8re en rien la probabilit\u00e9 d&rsquo;occurrence des \u00e9v\u00e9nements futurs. Or, dans les comportements addictifs, les conduites suicidaires ou les compulsions associ\u00e9es au trouble obsessionnel-compulsif, la r\u00e9alisation d&rsquo;un premier acte modifie profond\u00e9ment la dynamique psychique du sujet : elle brise des barri\u00e8res inhibitrices ou active des circuits de renforcement neurochimique qui rendent les \u00e9v\u00e9nements ult\u00e9rieurs bien plus probables (contagion positive). \u00c0 l&rsquo;inverse, l&rsquo;apparition d&rsquo;un \u00e9v\u00e9nement peut induire une p\u00e9riode r\u00e9fractaire au cours de laquelle toute nouvelle occurrence est impossible (contagion n\u00e9gative, induisant alors une sous-dispersion).<\/p>\n<p>Enfin, le regroupement des observations au sein de clusters \u00e9cologiques \u2014 tels que des patients trait\u00e9s au sein de diff\u00e9rentes cliniques, des \u00e9l\u00e8ves nich\u00e9s dans des classes scolaires ou des familles soumises \u00e0 un m\u00eame environnement \u2014 introduit des corr\u00e9lations intra-groupes qui gonflent artificiellement la variance empirique globale. Associ\u00e9e aux in\u00e9vitables erreurs de mesure inh\u00e9rentes aux inventaires auto-rapport\u00e9s (qui augmentent le bruit r\u00e9siduel), cette structure complexe d&rsquo;agr\u00e9gation d&rsquo;\u00e9v\u00e9nements rompt d\u00e9finitivement l&rsquo;\u00e9quidispersion th\u00e9orique.<\/p>\n<h3>3.3 Cons\u00e9quences statistiques de l&rsquo;ignorance du postulat<\/h3>\n<p>Lorsque la condition d&rsquo;\u00e9quidispersion est viol\u00e9e par la pr\u00e9sence d&rsquo;une surdispersion non trait\u00e9e, le recours na\u00eff au mod\u00e8le de r\u00e9gression de Poisson standard produit des cons\u00e9quences statistiques d\u00e9sastreuses sur la validit\u00e9 de l&rsquo;inf\u00e9rence. En vertu des propri\u00e9t\u00e9s th\u00e9oriques de l&rsquo;estimation quasi-vraisemblance d\u00e9montr\u00e9es par <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1982.10477803\">Gouri\u00e9roux, Monfort et Trognon (1984)<\/a>, les estimateurs ponctuels des coefficients de r\u00e9gression &beta;&#770; conservent leur propri\u00e9t\u00e9 de convergence et demeurent asymptotiquement sans biais, \u00e0 condition expresse que la fonction de lien logarithmique soit correctement sp\u00e9cifi\u00e9e : E(Y|X) = exp(X&beta;).<\/p>\n<p>Cependant, l&rsquo;ensemble de l&rsquo;inf\u00e9rence statistique s&rsquo;effondre au niveau de l&rsquo;estimation de l&rsquo;incertitude. La matrice de variance-covariance asymptotique standard, calcul\u00e9e sous l&rsquo;hypoth\u00e8se stricte Var(Y|X) = &mu;, sous-estime drastiquement les v\u00e9ritables erreurs-types des estimateurs :<\/p>\n<p>SE<sub>Poisson<\/sub>(&beta;&#770;<sub>j<\/sub>) &amp;Lt; SE<sub>R\u00e9elle<\/sub>(&beta;&#770;<sub>j<\/sub>)<\/p>\n<p>Cette contraction artificielle des erreurs-types a pour cons\u00e9quence directe une surestimation massive des statistiques de test de Wald :<\/p>\n<p>z = &beta;&#770;<sub>j<\/sub> \/ SE(&beta;&#770;<sub>j<\/sub>)<\/p>\n<p>Les p-valeurs calcul\u00e9es deviennent fallacieusement minuscules, provoquant une inflation incontr\u00f4l\u00e9e du taux d&rsquo;erreur de premi\u00e8re esp\u00e8ce (&alpha;). L&rsquo;analyste conclut alors \u00e0 l&rsquo;existence d&rsquo;effets th\u00e9rapeutiques hautement significatifs ou valide de pr\u00e9tendus corr\u00e9lats \u00e9tiologiques qui ne constituent en r\u00e9alit\u00e9 que des artefacts statistiques dus \u00e0 la surdispersion non contr\u00f4l\u00e9e. De surcro\u00eet, les intervalles de confiance construits sous ces erreurs-types sous-\u00e9valu\u00e9es affichent un taux de recouvrement empirique bien inf\u00e9rieur au seuil nominal de 95 %, sapant directement la reproductibilit\u00e9 scientifique des travaux publi\u00e9s.<\/p>\n<h2>4. La surdispersion statistique : m\u00e9canismes, formes et implications<\/h2>\n<h3>4.1 Ph\u00e9nom\u00e9nologie de la surdispersion (Overdispersion)<\/h3>\n<p>La surdispersion (<em>overdispersion<\/em> en anglais) d\u00e9signe l&rsquo;\u00e9tat statistique dans lequel la variance observ\u00e9e au sein des donn\u00e9es d\u00e9passe de fa\u00e7on statistiquement et substantiellement significative la variance postul\u00e9e par le mod\u00e8le th\u00e9orique sous-jacent. Dans le contexte d&rsquo;une distribution de Poisson, cela se traduit par la constatation empirique suivante :<\/p>\n<p>Var(Y | X) &gt; E(Y | X)<\/p>\n<p>Il est fondamental d&rsquo;\u00e9tablir une distinction analytique rigoureuse entre la surdispersion marginale et la surdispersion conditionnelle. La surdispersion marginale se rapporte aux donn\u00e9es brutes non ajust\u00e9es : lorsque l&rsquo;on calcule la moyenne et la variance univari\u00e9es de la variable d\u00e9pendante <em>Y<\/em>, il est fr\u00e9quent d&rsquo;observer une variance tr\u00e8s largement sup\u00e9rieure \u00e0 la moyenne. Toutefois, cette surdispersion brute peut simplement refl\u00e9ter le fait que des covariables puissantes n&rsquo;ont pas encore \u00e9t\u00e9 prises en compte. Une fois ces pr\u00e9dicteurs introduits dans le pr\u00e9dicteur lin\u00e9aire, la variance <em>conditionnelle<\/em> peut, en th\u00e9orie, s&rsquo;aligner sur la moyenne conditionnelle &mu;<sub>i<\/sub>. On parle alors de surdispersion apparente, caus\u00e9e par l&rsquo;omission de variables explicatives ou par une sp\u00e9cification erron\u00e9e de la forme fonctionnelle (comme l&rsquo;omission de termes polynomiaux ou d&rsquo;interactions non lin\u00e9aires).<\/p>\n<p>En revanche, la surdispersion r\u00e9elle (ou intrins\u00e8que) persiste m\u00eame apr\u00e8s avoir sp\u00e9cifi\u00e9 un mod\u00e8le complet et ad\u00e9quatement param\u00e9tr\u00e9. Elle traduit la pr\u00e9sence d&rsquo;une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 stochastique inh\u00e9rente au processus g\u00e9n\u00e9rateur de donn\u00e9es lui-m\u00eame. \u00c0 l&rsquo;oppos\u00e9, la sous-dispersion (<em>underdispersion<\/em>), situation o\u00f9 Var(Y|X) &lt; E(Y|X), demeure un ph\u00e9nom\u00e8ne excessivement rare dans les \u00e9tudes comportementales humaines. Elle ne s&rsquo;observe typiquement que dans des protocoles de laboratoire ultra-contr\u00f4l\u00e9s impliquant des processus d&rsquo;exclusion mutuelle ou des temps de latence physiologiquement contraignants limitant artificiellement la variabilit\u00e9 des d\u00e9nombrements.<\/p>\n<h3>4.2 Mod\u00e9lisation de l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 par un m\u00e9lange continu<\/h3>\n<p>La m\u00e9thode la plus rigoureuse et \u00e9l\u00e9gante pour formaliser la surdispersion r\u00e9elle consiste \u00e0 conceptualiser le processus de comptage individuel comme un mod\u00e8le de m\u00e9lange continu hi\u00e9rarchique, plus sp\u00e9cifiquement un m\u00e9lange continu de Poisson. Supposons que, conditionnellement \u00e0 la moyenne observ\u00e9e &mu;<sub>i<\/sub> et \u00e0 un facteur d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 latente inobserv\u00e9 &nu;<sub>i<\/sub> propre \u00e0 chaque individu, la variable d\u00e9pendante <em>Y<sub>i<\/sub><\/em> suive rigoureusement une loi de Poisson :<\/p>\n<p>Y<sub>i<\/sub> | &mu;<sub>i<\/sub>, &nu;<sub>i<\/sub> ~ P(&mu;<sub>i<\/sub> &times; &nu;<sub>i<\/sub>)<\/p>\n<p>o\u00f9 &nu;<sub>i<\/sub> &gt; 0 repr\u00e9sente un effet al\u00e9atoire multiplicatif non observ\u00e9, r\u00e9sumant l&rsquo;ensemble des pr\u00e9dispositions g\u00e9n\u00e9tiques, cognitives ou environnementales sp\u00e9cifiques \u00e0 l&rsquo;individu <em>i<\/em> qui n&rsquo;ont pu \u00eatre captur\u00e9es par les covariables mesur\u00e9es. Pour pr\u00e9server l&rsquo;identifiabilit\u00e9 du mod\u00e8le sans distordre la moyenne globale, nous imposons que ce facteur latent poss\u00e8de une esp\u00e9rance math\u00e9matique \u00e9gale \u00e0 l&rsquo;unit\u00e9 : E(&nu;<sub>i<\/sub>) = 1. Sa variance, en revanche, est un param\u00e8tre libre inconnu not\u00e9 &alpha;, tel que Var(&nu;<sub>i<\/sub>) = &alpha;, refl\u00e9tant l&rsquo;ampleur de l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 interindividuelle inobserv\u00e9e au sein de la population.<\/p>\n<p>Sur le plan math\u00e9matique, le choix de la distribution probabiliste pour ce terme latent &nu;<sub>i<\/sub> est guid\u00e9 par les principes de la th\u00e9orie bay\u00e9sienne et des distributions conjugu\u00e9es. La distribution continue la plus ad\u00e9quate pour mod\u00e9liser une variable al\u00e9atoire strictement positive de moyenne unitaire est la loi Gamma, param\u00e9tr\u00e9e par une forme <em>k<\/em> = 1\/&alpha; et un taux &theta; = 1\/&alpha; :<\/p>\n<p>&nu;<sub>i<\/sub> ~ Gamma(k = 1\/&alpha;, &theta; = 1\/&alpha;)<\/p>\n<p>En int\u00e9grant le param\u00e8tre latent &nu;<sub>i<\/sub> sur l&rsquo;ensemble de son domaine de d\u00e9finition stochastique [0, +&infin;[, la distribution marginale de <em>Y<sub>i<\/sub><\/em> (inconditionnelle par rapport \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 inobserv\u00e9e) s&rsquo;obtient analytiquement par convolution :<\/p>\n<p>P(Y<sub>i<\/sub> = y<sub>i<\/sub> | &mu;<sub>i<\/sub>, &alpha;) = &int;<sub>0<\/sub><sup>&infin;<\/sup> P(Y<sub>i<\/sub> = y<sub>i<\/sub> | &mu;<sub>i<\/sub>, &nu;<sub>i<\/sub>) &times; f(&nu;<sub>i<\/sub>; &alpha;) d&nu;<sub>i<\/sub><\/p>\n<p>Cette int\u00e9grale de m\u00e9lange Poisson-Gamma poss\u00e8de une solution ferm\u00e9e exacte : elle g\u00e9n\u00e8re pr\u00e9cis\u00e9ment la distribution de probabilit\u00e9 de la <strong>r\u00e9gression binomiale n\u00e9gative<\/strong>. Ce cadre math\u00e9matique conf\u00e8re une l\u00e9gitimit\u00e9 th\u00e9orique solide au mod\u00e8le binomial n\u00e9gatif, en l&rsquo;interpr\u00e9tant directement comme un mod\u00e8le de Poisson enrichi pour int\u00e9grer formellement la variabilit\u00e9 r\u00e9siduelle interindividuelle.<\/p>\n<h3>4.3 Impacts sur la puissance statistique et la pr\u00e9diction<\/h3>\n<p>Les cons\u00e9quences d&rsquo;une surdispersion non mod\u00e9lis\u00e9e ne se limitent pas \u00e0 l&rsquo;inflation des erreurs de type I : elles d\u00e9gradent profond\u00e9ment l&rsquo;exactitude des projections pr\u00e9dictives individuelles. Lorsqu&rsquo;un mod\u00e8le de Poisson est erron\u00e9ment ajust\u00e9 \u00e0 des donn\u00e9es surdispers\u00e9es, sa variance sous-estim\u00e9e produit des intervalles de pr\u00e9diction beaucoup trop \u00e9troits. Les cas cliniques extr\u00eames (par exemple, les patients \u00e0 tr\u00e8s haut risque de r\u00e9cidive ou \u00e0 tr\u00e8s forte charge symptomatique) sont alors consid\u00e9r\u00e9s comme des anomalies probabilistes impossibles par le mod\u00e8le, alors qu&rsquo;ils rel\u00e8vent simplement de la distribution r\u00e9elle \u00e0 queues lourdes.<\/p>\n<p>De plus, cette mauvaise sp\u00e9cification affaiblit la puissance statistique r\u00e9elle des tests d&rsquo;hypoth\u00e8ses conduits dans le cadre d&rsquo;essais cliniques randomis\u00e9s. En sous-estimant la variabilit\u00e9 r\u00e9siduelle authentique du ph\u00e9nom\u00e8ne sous \u00e9tude, les calculs de taille d&rsquo;\u00e9chantillon <em>a priori<\/em> fond\u00e9s sur l&rsquo;hypoth\u00e8se de Poisson sous-dimensionnent syst\u00e9matiquement les cohortes requises pour mettre en \u00e9vidence un effet th\u00e9rapeutique r\u00e9el avec une puissance ad\u00e9quate de 80 % ou 90 %. La prise en compte explicite de la surdispersion via un mod\u00e8le binomial n\u00e9gatif s&rsquo;av\u00e8re donc indispensable d\u00e8s la conception m\u00e9thodologique du protocole de recherche.<\/p>\n<h2>5. 5. Le mod\u00e8le de r\u00e9gression binomiale n\u00e9gative : cadre math\u00e9matique et param\u00e9trage<\/h2>\n<h3>5.1 Formulation probabiliste et fonctions de variance<\/h3>\n<p>Issue de la d\u00e9rivation du m\u00e9lange continu Poisson-Gamma, la r\u00e9gression binomiale n\u00e9gative r\u00e9introduit la flexibilit\u00e9 n\u00e9cessaire dans l&rsquo;ajustement des d\u00e9nombrements. Soit <em>Y<\/em> une variable de comptage dont l&rsquo;esp\u00e9rance conditionnelle est not\u00e9e &mu; = E(Y|X) = exp(X&beta;). La fonction de masse de probabilit\u00e9 univari\u00e9e de la distribution binomiale n\u00e9gative, exprim\u00e9e sous sa param\u00e9trisation moyenne-dispersion, prend la forme suivante :<\/p>\n<p>P(Y = y | &mu;, &alpha;) = [ &Gamma;(y + &alpha;<sup>-1<\/sup>) \/ ( y! &times; &Gamma;(&alpha;<sup>-1<\/sup>) ) ] &times; [ &alpha;<sup>-1<\/sup> \/ (&mu; + &alpha;<sup>-1<\/sup>) ]<sup>&alpha;<sup>-1<\/sup><\/sup> &times; [ &mu; \/ (&mu; + &alpha;<sup>-1<\/sup>) ]<sup>y<\/sup><\/p>\n<p>o\u00f9 &Gamma;(&middot;) d\u00e9signe la fonction gamma d&rsquo;Euler, qui g\u00e9n\u00e9ralise la factorielle aux nombres r\u00e9els telle que &Gamma;(n) = (n &#8211; 1)!. L&rsquo;\u00e9l\u00e9ment structural d\u00e9terminant de ce mod\u00e8le r\u00e9side dans sa <strong>fonction de variance quadratique<\/strong>, caract\u00e9ristique du mod\u00e8le canonique qualifi\u00e9 dans la litt\u00e9rature \u00e9conom\u00e9trique et statistique de mod\u00e8le <strong>NB2<\/strong> (<a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1017\/CBO9781139016650\">Cameron et Trivedi, 2013<\/a>) :<\/p>\n<p>Var(Y | X) = &mu; + &alpha; &times; &mu;<sup>2<\/sup><\/p>\n<p>Dans cette \u00e9quation, la variance conditionnelle est d\u00e9compos\u00e9e en deux termes distincts : le premier terme (&mu;) correspond \u00e0 la composante de variance de Poisson pure (variabilit\u00e9 intrins\u00e8que du d\u00e9nombrement), tandis que le second terme (&alpha;&mu;<sup>2<\/sup>) repr\u00e9sente la variabilit\u00e9 exc\u00e9dentaire g\u00e9n\u00e9r\u00e9e par l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 inobserv\u00e9e. Il existe \u00e9galement une formulation alternative, d\u00e9sign\u00e9e sous le nom de mod\u00e8le <strong>NB1<\/strong>, dans laquelle la fonction de variance est lin\u00e9aire : Var(Y|X) = &mu; + &alpha;&mu; = (1 + &alpha;)&mu;. Toutefois, dans les sciences du comportement, le mod\u00e8le quadratique NB2 s&rsquo;impose de fa\u00e7on pr\u00e9dominante en raison de son ad\u00e9quation remarquable avec l&rsquo;amplification non lin\u00e9aire de la variance observ\u00e9e aux valeurs moyennes \u00e9lev\u00e9es.<\/p>\n<p>Il est capital d&rsquo;observer le comportement analytique limite de la fonction de variance lorsque le param\u00e8tre de dispersion &alpha; tend vers z\u00e9ro :<\/p>\n<p>lim<sub>&alpha; &rarr; 0<\/sub> [Var(Y | X)] = lim<sub>&alpha; &rarr; 0<\/sub> (&mu; + &alpha;&mu;<sup>2<\/sup>) = &mu;<\/p>\n<p>De fa\u00e7on rigoureusement identique, la fonction de masse de probabilit\u00e9 binomiale n\u00e9gative converge ponctuellement vers la distribution de Poisson lorsque &alpha; &rarr; 0. Le mod\u00e8le de r\u00e9gression de Poisson appara\u00eet donc formellement comme un cas particulier nich\u00e9 (<em>nested<\/em>) au sein de la famille plus g\u00e9n\u00e9rale de la r\u00e9gression binomiale n\u00e9gative.<\/p>\n<h3>5.2 Le param\u00e8tre de dispersion alpha<\/h3>\n<p>Le param\u00e8tre &alpha; (souvent document\u00e9 comme le coefficient de surdispersion) quantifie l&rsquo;intensit\u00e9 de la variabilit\u00e9 non mod\u00e9lis\u00e9e. Une confusion terminologique et logicielle r\u00e9currente m\u00e9rite d&rsquo;\u00eatre soulign\u00e9e pour \u00e9viter toute erreur d&rsquo;interpr\u00e9tation lors de la manipulation des packages statistiques contemporains. Alors que la litt\u00e9rature th\u00e9orique privil\u00e9gie le param\u00e8tre de surdispersion directe &alpha; (o\u00f9 &alpha; = 0 d\u00e9signe l&rsquo;absence totale de surdispersion), des environnements majeurs tels que la fonction <code>glm.nb()<\/code> du package <code>MASS<\/code> sous le logiciel R param\u00e8trent le mod\u00e8le selon le terme &theta; = 1\/&alpha;, d\u00e9sign\u00e9 comme le param\u00e8tre d&rsquo;\u00e9chelle ou de taille (<em>shape\/size parameter<\/em>).<\/p>\n<p>Dans ce cadre param\u00e9trique alternatif, la relation limite s&rsquo;inverse : une absence de surdispersion (mod\u00e8le de Poisson pur) correspond \u00e0 &theta; &rarr; +&infin;. Inversement, une valeur empirique de &theta; faible (&theta; &lt; 1 ou &theta; proche de z\u00e9ro) signale une surdispersion particuli\u00e8rement s\u00e9v\u00e8re. L&rsquo;estimation de ce param\u00e8tre conjointement aux coefficients vectoriels &beta; s&rsquo;effectue usuellement par un algorithme d&rsquo;estimation par maximum de vraisemblance altern\u00e9 ou par profilage de vraisemblance : le syst\u00e8me optimise alternativement les coefficients de r\u00e9gression conditionnellement \u00e0 &alpha;, puis r\u00e9estime &alpha; conditionnellement aux coefficients jusqu&rsquo;\u00e0 convergence num\u00e9rique stable de la hessienne.<\/p>\n<p>Dans les cohortes cliniques de psychiatrie quantitative, des valeurs estim\u00e9es de &alpha; oscillant entre 0,20 et 1,50 sont extr\u00eamement courantes. Elles t\u00e9moignent de l&rsquo;impossibilit\u00e9 pratique de capter l&rsquo;int\u00e9gralit\u00e9 de la variabilit\u00e9 humaine \u00e0 l&rsquo;aide de simples pr\u00e9dicteurs lin\u00e9aires auto-rapport\u00e9s, l\u00e9gitimant l&rsquo;estimation empirique de ce param\u00e8tre d&rsquo;\u00e9chelle additionnel.<\/p>\n<h3>5.3 Interpr\u00e9tation des estimateurs en r\u00e9gression binomiale n\u00e9gative<\/h3>\n<p>L&rsquo;un des avantages conceptuels majeurs du mod\u00e8le binomial n\u00e9gatif r\u00e9side dans la pr\u00e9servation int\u00e9grale de la fonction de lien logarithmique. D\u00e8s lors que l&rsquo;on mod\u00e9lise ln(&mu;) = X&beta;, l&rsquo;interpr\u00e9tation s\u00e9mantique et clinique des coefficients de r\u00e9gression &beta;&#770; est <strong>rigoureusement identique<\/strong> \u00e0 celle expos\u00e9e pr\u00e9c\u00e9demment pour la r\u00e9gression de Poisson. L&rsquo;exponentiation des coefficients, exp(&beta;&#770;), produit exactement les Ratios de Taux d&rsquo;Incidence (IRR), autorisant une comparabilit\u00e9 directe des amplitudes d&rsquo;effet entre \u00e9tudes.<\/p>\n<p>La divergence analytique fondamentale s&rsquo;exprime exclusivement dans l&rsquo;\u00e9valuation de la pr\u00e9cision de ces estimations : les erreurs-types estim\u00e9es sous le mod\u00e8le binomial n\u00e9gatif int\u00e8grent la composante de variance exc\u00e9dentaire &alpha;&mu;<sup>2<\/sup>. En cons\u00e9quence, les erreurs-types sont plus larges et r\u00e9alistes que celles issues d&rsquo;un mod\u00e8le de Poisson mal ajust\u00e9. Les intervalles de confiance \u00e0 95 % s&rsquo;\u00e9largissent pour refl\u00e9ter la v\u00e9ritable incertitude exp\u00e9rimentale, et les p-valeurs augmentent jusqu&rsquo;\u00e0 des seuils scientifiquement d\u00e9fendables. N\u00e9anmoins, il importe de rappeler que la r\u00e9gression binomiale n\u00e9gative ne constitue pas une panac\u00e9e universelle : si la surdispersion est exclusivement engendr\u00e9e par une abondance massive de scores nuls r\u00e9sultant d&rsquo;un processus binaire distinct, le mod\u00e8le NB2 peut demeurer sous-optimal face aux mod\u00e8les d\u00e9di\u00e9s \u00e0 l&rsquo;inflation de z\u00e9ros.<\/p>\n<h2>6. 6. Comparaison empirique et m\u00e9thodologique : Poisson versus Binomiale n\u00e9gative<\/h2>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/negativeBinom1.png\" alt=\"Negative binomial vs. Poisson regression\" \/><figcaption>Negative binomial vs. Poisson regression<\/figcaption><\/figure>\n<h3>6.1 Tableau comparatif des propri\u00e9t\u00e9s structurelles<\/h3>\n<p>Afin de synth\u00e9tiser les divergences structurelles majeures distinguant les deux mod\u00e8les fondamentaux de donn\u00e9es de comptage, le tableau suivant r\u00e9capitule leurs propri\u00e9t\u00e9s math\u00e9matiques, computationnelles et inf\u00e9rentielles fondamentales :<\/p>\n<ul>\n<li><strong>Hypoth\u00e8se distributionnelle sous-jacente :<\/strong> Le mod\u00e8le de Poisson repose sur une loi univari\u00e9e \u00e0 un seul param\u00e8tre <em>P(&lambda;)<\/em>. Le mod\u00e8le binomial n\u00e9gatif (NB2) postule une distribution de m\u00e9lange continu bivari\u00e9e Poisson-Gamma.<\/li>\n<li><strong>Formulation de la fonction de variance :<\/strong> Poisson contraint la variance \u00e0 l&rsquo;\u00e9galit\u00e9 formelle avec la moyenne, Var(Y|X) = &mu;. La binomiale n\u00e9gative incorpore une relation quadratique extensible, Var(Y|X) = &mu; + &alpha;&mu;<sup>2<\/sup>.<\/li>\n<li><strong>Param\u00e8tres estim\u00e9s :<\/strong> Poisson n&rsquo;estime que le vecteur des coefficients de r\u00e9gression &beta; (dimension <em>p<\/em>). La binomiale n\u00e9gative estime le vecteur &beta; ainsi que le param\u00e8tre de dispersion d&rsquo;\u00e9chelle &alpha; (dimension <em>p + 1<\/em>).<\/li>\n<li><strong>Comportement asymptotique :<\/strong> Lorsque la taille d&rsquo;\u00e9chantillon <em>N<\/em> cro\u00eet et que &alpha; &rarr; 0, le mod\u00e8le binomial n\u00e9gatif converge asymptotiquement vers l&rsquo;estimateur de Poisson. En petits \u00e9chantillons homog\u00e8nes sans surdispersion, l&rsquo;ajustement de la binomiale n\u00e9gative peut induire une instabilit\u00e9 num\u00e9rique mineure li\u00e9e \u00e0 l&rsquo;estimation de &alpha; sur la fronti\u00e8re param\u00e9trique.<\/li>\n<li><strong>Sensibilit\u00e9 aux valeurs extr\u00eames :<\/strong> Le mod\u00e8le de Poisson est excessivement sensible aux valeurs isol\u00e9es de comptage \u00e9lev\u00e9, lesquelles tirent anormalement la moyenne conditionnelle et g\u00e9n\u00e8rent des r\u00e9sidus disproportionn\u00e9s. Le mod\u00e8le binomial n\u00e9gatif absorbe ces observations \u00e9loign\u00e9es avec une robustesse substantiellement accrue gr\u00e2ce \u00e0 sa fonction de variance tol\u00e9rante.<\/li>\n<\/ul>\n<h3>6.2 Biais d&rsquo;estimation et robustesse comparative<\/h3>\n<p>Sur le plan des propri\u00e9t\u00e9s statistiques formelles des estimateurs, un paradoxe m\u00e9thodologique apparent \u00e9merge des comparaisons th\u00e9oriques. D&rsquo;un c\u00f4t\u00e9, le mod\u00e8le de Poisson jouit d&rsquo;une robustesse structurelle reconnue par le th\u00e9or\u00e8me de Gouri\u00e9roux et coll\u00e8gues : si la relation moyenne est correctement mod\u00e9lis\u00e9e sous la forme exp(X&beta;), l&rsquo;estimateur &beta;&#770;<sub>Poisson<\/sub> demeure rigoureusement sans biais et convergent, m\u00eame si la distribution r\u00e9elle s&rsquo;av\u00e8re binomiale n\u00e9gative, gamma ou log-normale. Son talon d&rsquo;Achille r\u00e9side exclusivement dans sa matrice de variance empirique.<\/p>\n<p>D&rsquo;un autre c\u00f4t\u00e9, bien que le mod\u00e8le binomial n\u00e9gatif r\u00e9solve magistralement le probl\u00e8me de l&rsquo;estimation de la variance r\u00e9siduelle en ajustant formellement le terme &alpha;, il repose th\u00e9oriquement sur l&rsquo;exactitude de la distribution de m\u00e9lange Gamma. Si le v\u00e9ritable processus g\u00e9n\u00e9rateur d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 latente d\u00e9vie drastiquement d&rsquo;une loi Gamma (en suivant par exemple une distribution bimodale ou log-normale discontinue), les estimations conjointes de &beta; et &alpha; peuvent theoretically \u00eatre affect\u00e9es dans des conditions extr\u00eames d&rsquo;\u00e9chantillonnage fini. Cependant, dans la quasi-totalit\u00e9 des simulations mont\u00e9es en psychom\u00e9trie, le co\u00fbt statistique li\u00e9 \u00e0 l&rsquo;estimation du param\u00e8tre additionnel &alpha; dans le mod\u00e8le NB2 est infinit\u00e9simal en regard du gain spectaculaire r\u00e9alis\u00e9 en mati\u00e8re de validit\u00e9 de l&rsquo;inf\u00e9rence par rapport au mod\u00e8le de Poisson non corrig\u00e9.<\/p>\n<h3>6.3 Sc\u00e9narios typiques orientant vers l&rsquo;un ou l&rsquo;autre mod\u00e8le<\/h3>\n<p>Dans la pratique de la recherche acad\u00e9mique, le choix entre ces deux paradigmes s&rsquo;aligne \u00e9troitement sur le design d&rsquo;investigation et la nature des construits psychologiques investigu\u00e9s. Le mod\u00e8le de <strong>Poisson standard<\/strong> conserve toute sa pertinence et sa sup\u00e9riorit\u00e9 parcimonieuse dans les contextes hautement contr\u00f4l\u00e9s :<\/p>\n<ul>\n<li>Protocoles exp\u00e9rimentaux rigoureux en laboratoire o\u00f9 chaque sujet subit un protocole de stimulation parfaitement standardis\u00e9 et temporalis\u00e9 (ex. t\u00e2ches de temps de r\u00e9action avec d\u00e9compte d&rsquo;omissions rares).<\/li>\n<li>\u00c9tudes d&rsquo;\u00e9pid\u00e9miologie clinique sur des cohortes d&rsquo;une extr\u00eame homog\u00e9n\u00e9it\u00e9 portant sur des \u00e9v\u00e9nements rares survenant au sein de fen\u00eatres temporelles courtes.<\/li>\n<li>Situations o\u00f9 la variance empirique calcul\u00e9e au sein des cellules factorielles du plan exp\u00e9rimental se r\u00e9v\u00e8le sensiblement \u00e9quivalente \u00e0 leurs moyennes arithm\u00e9tiques respectives.<\/li>\n<\/ul>\n<p>\u00c0 l&rsquo;oppos\u00e9, la <strong>r\u00e9gression binomiale n\u00e9gative<\/strong> s&rsquo;av\u00e8re imp\u00e9rative et m\u00e9thodologiquement obligatoire dans les contextes suivants :<\/p>\n<ul>\n<li>\u00c9tudes en psychopathologie d\u00e9veloppementale mesurant des \u00e9pisodes comportementaux au sein de la population g\u00e9n\u00e9rale (ex. nombre de crises de col\u00e8re, comportements d\u00e9linquants, fugues), caract\u00e9ris\u00e9s par une immense variabilit\u00e9 interindividuelle et une pr\u00e9dominance de faibles scores accompagn\u00e9e d&rsquo;une fraction d&rsquo;individus cumulant des centaines d&rsquo;occurrences.<\/li>\n<li>Enqu\u00eates longitudinales ou transversales fond\u00e9es sur des questionnaires auto-administr\u00e9s \u00e9valuant la fr\u00e9quence de consommation de substances licites ou illicites (alcool, nicotine, cannabis).<\/li>\n<li>Donn\u00e9es d&rsquo;observation en milieu naturel (\u00e9coles, services d&rsquo;urgence, prisons) o\u00f9 des ph\u00e9nom\u00e8nes d&rsquo;agr\u00e9gation institutionnelle et d&rsquo;influence mutuelle gonflent irr\u00e9m\u00e9diablement la dispersion r\u00e9siduelle.<\/li>\n<\/ul>\n<h2>7. 7. Techniques diagnostiques 1 : Analyse des r\u00e9sidus et graphiques de r\u00e9sidus standardis\u00e9s<\/h2>\n<h3>7.1 Typologie des r\u00e9sidus pour mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s<\/h3>\n<p>Dans le cadre des mod\u00e8les de r\u00e9gression lin\u00e9aire classique par OLS, l&rsquo;analyse des r\u00e9sidus est conceptuellement directe : le r\u00e9sidu brut e<sub>i<\/sub> = y<sub>i<\/sub> &#8211; y&#770;<sub>i<\/sub> poss\u00e8de, sous l&rsquo;hypoth\u00e8se de validit\u00e9 du mod\u00e8le, une variance constante &sigma;<sup>2<\/sup>. Dans le p\u00e9rim\u00e8tre des GLM, cette d\u00e9finition brute est totalement obsol\u00e8te, car la variance th\u00e9orique d\u00e9pend intrins\u00e8quement de la valeur pr\u00e9dite &mu;<sub>i<\/sub>. Pour \u00e9valuer la qualit\u00e9 d&rsquo;ajustement d&rsquo;un mod\u00e8le de comptage, les statisticiens emploient deux grandes cat\u00e9gories de r\u00e9sidus transform\u00e9s : les r\u00e9sidus de Pearson et les r\u00e9sidus de d\u00e9viance.<\/p>\n<p>Le <strong>r\u00e9sidu de Pearson<\/strong> ajuste la d\u00e9viation brute en la divisant par l&rsquo;\u00e9cart-type conditionnel pr\u00e9dit par la fonction de variance <em>V(&mu;)<\/em> du mod\u00e8le examin\u00e9. Pour le mod\u00e8le de Poisson, o\u00f9 V(&mu;<sub>i<\/sub>) = &mu;<sub>i<\/sub>, et pour le mod\u00e8le binomial n\u00e9gatif (NB2), o\u00f9 V(&mu;<sub>i<\/sub>) = &mu;<sub>i<\/sub> + &alpha;&mu;<sub>i<\/sub><sup>2<\/sup>, le r\u00e9sidu de Pearson individuel s&rsquo;\u00e9crit respectivement :<\/p>\n<p>r<sub>P, i<\/sub><sup>Poisson<\/sup> = (y<sub>i<\/sub> &#8211; &mu;&#770;<sub>i<\/sub>) \/ &radic;(&mu;&#770;<sub>i<\/sub>) &emsp; et &emsp; r<sub>P, i<\/sub><sup>NB<\/sup> = (y<sub>i<\/sub> &#8211; &mu;&#770;<sub>i<\/sub>) \/ &radic;(&mu;&#770;<sub>i<\/sub> + &alpha;&#770; &mu;&#770;<sub>i<\/sub><sup>2<\/sup>)<\/p>\n<p>Le <strong>r\u00e9sidu de d\u00e9viance<\/strong>, not\u00e9 <em>d<sub>i<\/sub><\/em>, quantifie quant \u00e0 lui la contribution alg\u00e9brique individuelle de chaque observation \u00e0 la statistique de d\u00e9viance r\u00e9siduelle globale du mod\u00e8le (la d\u00e9viance \u00e9tant le double de la diff\u00e9rence de log-vraisemblance entre le mod\u00e8le satur\u00e9 parfait et le mod\u00e8le ajust\u00e9). Il s&rsquo;exprime sous la forme <em>d<sub>i<\/sub> = sign(y<sub>i<\/sub> &#8211; &mu;&#770;<sub>i<\/sub>) &times; &radic;(d<sub>i<\/sub><sup>2<\/sup>)<\/em>, offrant une distribution souvent plus proche de la normalit\u00e9 asymptotique que les r\u00e9sidus de Pearson.<\/p>\n<p>Pour affiner ces diagnostics, il est n\u00e9cessaire de calculer les <strong>r\u00e9sidus standardis\u00e9s<\/strong> (souvent appel\u00e9s r\u00e9sidus studentis\u00e9s en OLS). Cette op\u00e9ration pond\u00e8re le r\u00e9sidu par la racine de <em>(1 &#8211; h<sub>ii<\/sub>)<\/em>, o\u00f9 <em>h<sub>ii<\/sub><\/em> d\u00e9signe le <em>leverage<\/em> (l&rsquo;\u00e9l\u00e9ment diagonal correspondant de la matrice chapeau <em>H<\/em>) refl\u00e9tant l&rsquo;effet de levier exerc\u00e9 par l&rsquo;individu dans l&rsquo;espace multidimensionnel des covariables :<\/p>\n<p>r<sub>P, i<\/sub><sup>std<\/sup> = r<sub>P, i<\/sub> \/ &radic;(1 &#8211; h<sub>ii<\/sub>)<\/p>\n<h3>7.2 Construction et interpr\u00e9tation du graphique des r\u00e9sidus standardis\u00e9s<\/h3>\n<p>L&rsquo;inspection visuelle du graphique mettant en relation les <strong>r\u00e9sidus de Pearson standardis\u00e9s<\/strong> (en ordonn\u00e9e) avec les <strong>valeurs pr\u00e9dites<\/strong> &mu;&#770; (ou leur transformation logarithmique, en abscisse) constitue la proc\u00e9dure diagnostique empirique la plus r\u00e9v\u00e9latrice de la surdispersion. Si le mod\u00e8le de r\u00e9gression de Poisson sous-jacent est parfaitement sp\u00e9cifi\u00e9 et que l&rsquo;hypoth\u00e8se d&rsquo;\u00e9quidispersion est rigoureusement satisfaite, les r\u00e9sidus standardis\u00e9s doivent se disperser al\u00e9atoirement autour de la ligne horizontale de r\u00e9f\u00e9rence z\u00e9ro, sans structure apparente, et leur variance empirique doit \u00eatre approximativement \u00e9gale \u00e0 1.<\/p>\n<p>Deux r\u00e8gles heuristiques formelles r\u00e9gissent la lecture critique de ce trac\u00e9 diagnostique :<\/p>\n<ul>\n<li><strong>La r\u00e8gle de l&rsquo;intervalle canonique [-2, +2] :<\/strong> Dans un mod\u00e8le ad\u00e9quatement \u00e9quidispens\u00e9, approximativement 95 % des r\u00e9sidus standardis\u00e9s doivent imp\u00e9rativement s&rsquo;inscrire \u00e0 l&rsquo;int\u00e9rieur de la bande d\u00e9limit\u00e9e par les seuils ordonn\u00e9s [-2, +2], conform\u00e9ment aux propri\u00e9t\u00e9s d&rsquo;une loi normale centr\u00e9e r\u00e9duite de r\u00e9f\u00e9rence.<\/li>\n<li><strong>Le seuil critique d&rsquo;alerte [-3, +3] :<\/strong> La pr\u00e9sence r\u00e9currente de multiples observations exc\u00e9dant la plage normative [-3, +3] constitue un signal diagnostique indiscutable de surdispersion s\u00e9v\u00e8re. Si les r\u00e9sidus s&rsquo;\u00e9talent sur des amplitudes atteignant des valeurs de +5, +8 ou +15, l&rsquo;\u00e9quidispersion de Poisson est cat\u00e9goriquement invalid\u00e9e.<\/li>\n<\/ul>\n<p>De plus, l&rsquo;analyste doit traquer l&rsquo;apparition de profils morphologiques caract\u00e9ristiques tels que les <em>motifs en entonnoir<\/em> (\u00e9vasement progressif de l&rsquo;\u00e9talement vertical des r\u00e9sidus \u00e0 mesure que les valeurs pr\u00e9dites augmentent). Ce motif en entonnoir trahit l&rsquo;incapacit\u00e9 de la fonction de variance lin\u00e9aire de Poisson \u00e0 contenir l&rsquo;amplification d&rsquo;amplitude des r\u00e9sidus, signalant qu&rsquo;une fonction de variance quadratique de type NB2 est structurellement requise.<\/p>\n<h3>7.3 Diagnostics graphiques compl\u00e9mentaires<\/h3>\n<p>Au-del\u00e0 du traditionnel trac\u00e9 des r\u00e9sidus de Pearson, la bo\u00eete \u00e0 outils diagnostique moderne met \u00e0 disposition des repr\u00e9sentations graphiques plus sophistiqu\u00e9es. Les <strong>r\u00e9sidus quantiles al\u00e9atoires<\/strong>, d\u00e9mocratis\u00e9s par le package de r\u00e9f\u00e9rence <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=DHARMa\">DHARMa<\/a> dans l&rsquo;environnement R, r\u00e9solvent le probl\u00e8me de la discr\u00e9tisation des r\u00e9sidus dans les mod\u00e8les GLM. En transformant les r\u00e9sidus discrets en quantiles continus uniform\u00e9ment distribu\u00e9s via des simulations probabilistes du mod\u00e8le sous hypoth\u00e8se nulle, ils permettent la g\u00e9n\u00e9ration de diagrammes quantile-quantile (Q-Q plots) d&rsquo;une clart\u00e9 remarquable. Tout \u00e9cart syst\u00e9matique de la trajectoire des points par rapport \u00e0 la bissectrice th\u00e9orique fournit un test visuel infaillible de surdispersion ou d&rsquo;inflation de z\u00e9ros.<\/p>\n<p>Une autre technique visuelle remarquable est le <strong>rootogramme suspendu<\/strong> (<em>hanging rootogram<\/em>), conceptualis\u00e9 \u00e0 l&rsquo;origine par John Tukey et adapt\u00e9 aux GLM de comptage par <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00031305.2016.1173590\">Kleiber et Zeileis (2016)<\/a>. En tra\u00e7ant la racine carr\u00e9e des fr\u00e9quences observ\u00e9es et attendues de telle sorte que les barres pendent litt\u00e9ralement depuis la courbe des valeurs attendues jusqu&rsquo;\u00e0 l&rsquo;axe horizontal, le rootogramme met imm\u00e9diatement en relief les d\u00e9nombrements pour lesquels le mod\u00e8le s&rsquo;av\u00e8re d\u00e9ficitaire ou exc\u00e9dentaire. Il visualise instantan\u00e9ment si le mod\u00e8le sous-estime syst\u00e9matiquement la fr\u00e9quence des z\u00e9ros ou \u00e9choue \u00e0 mod\u00e9liser la tra\u00eene lourde des observations massives.<\/p>\n<h2>8. 8. Techniques diagnostiques 2 : Tests statistiques formels de surdispersion et comparaison de mod\u00e8les<\/h2>\n<h3>8.1 Le test du rapport de vraisemblance (LRT) appliqu\u00e9 \u00e0 la dispersion<\/h3>\n<p>Si l&rsquo;inspection visuelle des r\u00e9sidus fournit des indices pr\u00e9cieux, la rigueur m\u00e9thodologique impose la conduite de tests d&rsquo;hypoth\u00e8ses formels. Le test le plus largement d\u00e9ploy\u00e9 pour discriminer le mod\u00e8le de Poisson du mod\u00e8le binomial n\u00e9gatif est le <strong>Test du Rapport de Vraisemblance<\/strong> (<em>Likelihood Ratio Test<\/em>, ou LRT). Ce test \u00e9value formellement l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9quidispersion de Poisson contre l&rsquo;hypoth\u00e8se alternative de surdispersion binomiale n\u00e9gative :<\/p>\n<p>H<sub>0<\/sub> : &alpha; = 0 &emsp; versus &emsp; H<sub>1<\/sub> : &alpha; &gt; 0<\/p>\n<p>La statistique de test standard repose sur la diff\u00e9rence des log-vraisemblances maximales obtenues respectivement pour le mod\u00e8le binomial n\u00e9gatif restreint (ln L<sub>NB<\/sub>) et le mod\u00e8le de Poisson embo\u00eet\u00e9 (ln L<sub>Poisson<\/sub>) :<\/p>\n<p>LRT = 2 &times; [ ln L<sub>NB<\/sub>(&beta;&#770;<sub>NB<\/sub>, &alpha;&#770;) &#8211; ln L<sub>Poisson<\/sub>(&beta;&#770;<sub>Poisson<\/sub>) ]<\/p>\n<p>Cependant, ce test soul\u00e8ve une difficult\u00e9 statistique majeure trop souvent ignor\u00e9e : l&rsquo;hypoth\u00e8se nulle &alpha; = 0 se situe exactement sur la <strong>fronti\u00e8re<\/strong> de l&rsquo;espace param\u00e9trique admissible, puisque le param\u00e8tre de dispersion &alpha; ne peut math\u00e9matiquement pas adopter de valeurs strictement n\u00e9gatives dans la formulation standard NB2. En vertu des d\u00e9monstrations de Chernoff (1954) et Self et Liang (1987), la distribution asymptotique sous l&rsquo;hypoth\u00e8se nulle d&rsquo;un test sur la fronti\u00e8re n&rsquo;est pas une loi classique du &chi;<sup>2<\/sup> \u00e0 1 degr\u00e9 de libert\u00e9, mais un m\u00e9lange probabiliste 50:50 d&rsquo;une masse ponctuelle en z\u00e9ro (Dirac &delta;<sub>0<\/sub>) et d&rsquo;un &chi;<sup>2<\/sup> \u00e0 1 degr\u00e9 de libert\u00e9, not\u00e9 formellement :<\/p>\n<p>LRT ~ 0.5 &times; &chi;<sub>0<\/sub><sup>2<\/sup> + 0.5 &times; &chi;<sub>1<\/sub><sup>2<\/sup><\/p>\n<p>Pour corriger cette contrainte de fronti\u00e8re et ne pas employer un test excessivement conservateur, la v\u00e9ritable p-valeur du LRT corrig\u00e9 pour tester l&rsquo;adjonction du param\u00e8tre de dispersion s&rsquo;obtient simplement en divisant par deux la p-valeur d\u00e9riv\u00e9e de la distribution usuelle du &chi;<sup>2<\/sup>(1) :<\/p>\n<p>p-valeur<sub>corrig\u00e9e<\/sub> = 0.5 &times; P(&chi;<sub>1<\/sub><sup>2<\/sup> &gt; LRT<sub>observ\u00e9<\/sub>)<\/p>\n<p>Si cette p-valeur corrig\u00e9e est inf\u00e9rieure au seuil critique nominal (typiquement 0,05), l&rsquo;hypoth\u00e8se nulle &alpha; = 0 est rejet\u00e9e avec confiance, attestant formellement de la sup\u00e9riorit\u00e9 statistique du mod\u00e8le binomial n\u00e9gatif.<\/p>\n<h3>8.2 Tests de score et statistiques de Wald pour la dispersion<\/h3>\n<p>L&rsquo;ex\u00e9cution du test de rapport de vraisemblance exige d&rsquo;ajuster cons\u00e9cutivement les deux mod\u00e8les, ce qui peut se r\u00e9v\u00e9ler computationnellement co\u00fbteux face \u00e0 de gigantesques bases de donn\u00e9es. Pour pallier cette contrainte, Cameron et Trivedi ont d\u00e9riv\u00e9 un <strong>test de score<\/strong> \u00e9l\u00e9gant qui ne n\u00e9cessite que l&rsquo;ajustement pr\u00e9alable du mod\u00e8le de r\u00e9gression de Poisson standard. Le test de Cameron-Trivedi \u00e9value la relation entre les r\u00e9sidus de Poisson et les valeurs pr\u00e9dites en ajustant une r\u00e9gression lin\u00e9aire auxiliaire sans constante :<\/p>\n<p>[ (y<sub>i<\/sub> &#8211; &mu;&#770;<sub>i<\/sub>)<sup>2<\/sup> &#8211; y<sub>i<\/sub> ] \/ &mu;&#770;<sub>i<\/sub> = &alpha; &times; g(&mu;&#770;<sub>i<\/sub>) + &epsilon;<sub>i<\/sub><\/p>\n<p>Pour tester sp\u00e9cifiquement l&rsquo;architecture de variance quadratique NB2, la fonction choisie est <em>g(&mu;&#770;<sub>i<\/sub>) = &mu;&#770;<sub>i<\/sub><\/em>. La significativit\u00e9 de la pente estim\u00e9e &alpha;&#770; au moyen d&rsquo;un simple test <em>t<\/em> de Student unilat\u00e9ral permet de conclure instantan\u00e9ment \u00e0 la pr\u00e9sence de surdispersion, sans avoir \u00e0 faire converger l&rsquo;algorithme d&rsquo;optimisation de la binomiale n\u00e9gative.<\/p>\n<p>Par ailleurs, les chercheurs inspectent couramment l&rsquo;indice heuristique de dispersion bas\u00e9 sur la <strong>d\u00e9viance r\u00e9siduelle<\/strong> ou le <strong>&chi;<sup>2<\/sup> de Pearson<\/strong> divis\u00e9 par les degr\u00e9s de libert\u00e9 r\u00e9siduels (<em>df<\/em>) :<\/p>\n<p>&Phi;&#770; = &chi;<sub>Pearson<\/sub><sup>2<\/sup> \/ (n &#8211; k) = [ &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &#8211; &mu;&#770;<sub>i<\/sub>)<sup>2<\/sup> \/ &mu;&#770;<sub>i<\/sub> ] \/ (n &#8211; k)<\/p>\n<p>Sous l&rsquo;hypoth\u00e8se d&rsquo;\u00e9quidispersion stricte, ce ratio &Phi;&#770; doit \u00eatre sensiblement \u00e9gal \u00e0 1. Une valeur de &Phi;&#770; exc\u00e9dant substantiellement 1,20 ou 1,50 fournit une indication diagnostique empirique forte d&rsquo;une violation de l&rsquo;hypoth\u00e8se de Poisson. Toutefois, il convient d&rsquo;\u00eatre vigilant : ce ratio est sensible aux erreurs de sp\u00e9cification de la fonction de lien et \u00e0 la taille d&rsquo;\u00e9chantillon, et ne doit jamais se substituer aux tests formels de score ou de rapport de vraisemblance.<\/p>\n<h3>8.3 Crit\u00e8res d&rsquo;information et validation crois\u00e9e<\/h3>\n<p>Lorsque l&rsquo;on compare des architectures de mod\u00e9lisation concurrentes qui ne sont pas n\u00e9cessairement parfaitement embo\u00eet\u00e9es, les crit\u00e8res d&rsquo;information p\u00e9nalis\u00e9s constituent des rep\u00e8res d\u00e9cisionnels indispensables. Le crit\u00e8re d&rsquo;information d&rsquo;Akaike (<a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1109\/TAC.1974.1100705\">AIC<\/a>) et le crit\u00e8re d&rsquo;information bay\u00e9sien (<a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aos\/1176344136\">BIC<\/a>) op\u00e8rent un arbitrage quantitatif rigoureux entre la qualit\u00e9 de l&rsquo;ajustement aux donn\u00e9es (captur\u00e9e par la log-vraisemblance maximale ln L) et la complexit\u00e9 du mod\u00e8le (p\u00e9nalis\u00e9e par le nombre de param\u00e8tres estim\u00e9s <em>K<\/em>) :<\/p>\n<p>AIC = -2 &times; ln L + 2K &emsp; et &emsp; BIC = -2 &times; ln L + K &times; ln(n)<\/p>\n<p>Dans la comp\u00e9tition entre Poisson (K param\u00e8tres) et Binomiale N\u00e9gative (K + 1 param\u00e8tres), la binomiale n\u00e9gative n&rsquo;est p\u00e9nalis\u00e9e que par l&rsquo;adjonction d&rsquo;une seule unit\u00e9 param\u00e9trique (&Delta;K = 1). D\u00e8s lors que la surdispersion est pr\u00e9sente, le gain en log-vraisemblance procur\u00e9 par l&rsquo;int\u00e9gration de la fonction de variance quadratique d\u00e9passe presque syst\u00e9matiquement la faible p\u00e9nalit\u00e9 de complexit\u00e9 de 2 points d&rsquo;AIC. Selon les directives classiques de Burnham et Anderson (2002), un diff\u00e9rentiel &Delta;AIC = AIC<sub>Poisson<\/sub> &#8211; AIC<sub>NB<\/sub> sup\u00e9rieur \u00e0 10 d\u00e9montre un support empirique \u00e9crasant en faveur du mod\u00e8le binomial n\u00e9gatif, rel\u00e9guant le mod\u00e8le de Poisson au rang de sp\u00e9cification erron\u00e9e.<\/p>\n<p>Pour les chercheurs privil\u00e9giant l&rsquo;exactitude pr\u00e9dictive hors-\u00e9chantillon plut\u00f4t que la seule th\u00e9orie de l&rsquo;information, la <strong>validation crois\u00e9e par <em>k<\/em>-blocs<\/strong> (<em>k-fold cross-validation<\/em>) apporte une validation empirique suppl\u00e9mentaire. En scindant les donn\u00e9es en sous-ensembles d&rsquo;entra\u00eenement et de validation, elle permet de comparer la racine de l&rsquo;erreur quadratique moyenne de pr\u00e9diction (RMSEP) ou le score de probabilit\u00e9 logarithmique (<em>log-score<\/em> probabiliste). Cette d\u00e9marche confirme syst\u00e9matiquement que face \u00e0 des donn\u00e9es surdispers\u00e9es, le mod\u00e8le binomial n\u00e9gatif g\u00e9n\u00e8re des pr\u00e9dictions de densit\u00e9 de probabilit\u00e9 bien mieux \u00e9talonn\u00e9es sur les donn\u00e9es de test non observ\u00e9es.<\/p>\n<h2>9. 9. Gestion des cas particuliers : Inflation de z\u00e9ros, troncature et alternatives semi-param\u00e9triques<\/h2>\n<h3>9.1 Le ph\u00e9nom\u00e8ne d&rsquo;exc\u00e8s de z\u00e9ros (Zero-Inflation)<\/h3>\n<p>L&rsquo;un des d\u00e9fis analytiques les plus r\u00e9pandus en psychopathologie quantitative r\u00e9side dans la prolif\u00e9ration massive de scores nuls au sein de la variable d\u00e9pendante. Par exemple, lors de l&rsquo;\u00e9valuation du nombre d&rsquo;arrestations polici\u00e8res au cours des six derniers mois dans un \u00e9chantillon de la population g\u00e9n\u00e9rale, une \u00e9crasante majorit\u00e9 des participants pr\u00e9sentera un score rigoureusement \u00e9gal \u00e0 z\u00e9ro. Si une proportion importante de z\u00e9ros peut naturellement \u00eatre absorb\u00e9e par le mod\u00e8le binomial n\u00e9gatif gr\u00e2ce \u00e0 sa grande flexibilit\u00e9 d&rsquo;\u00e9chelle, il survient des configurations o\u00f9 la masse en z\u00e9ro transcende ce que n&rsquo;importe quelle distribution univari\u00e9e standard peut mod\u00e9liser.<\/p>\n<p>Il est fondamental de distinguer deux sources \u00e9tiologiques de z\u00e9ros :<\/p>\n<ul>\n<li><strong>Les z\u00e9ros structurels :<\/strong> Ils correspondent \u00e0 des individus qui ne sont conceptuellement pas \u00e0 risque de produire l&rsquo;\u00e9v\u00e9nement \u00e9tudi\u00e9 (par exemple, des participants totalement abstinents de toute consommation d&rsquo;alcool dans une \u00e9tude sur les \u00e9pisodes de fr\u00e9n\u00e9sie \u00e9thylique).<\/li>\n<li><strong>Les z\u00e9ros d&rsquo;\u00e9chantillonnage (ou accidentels) :<\/strong> Ils concernent des individus qui appartiennent bel et bien \u00e0 la population susceptible de manifester le comportement, mais qui se trouve n&rsquo;avoir commis aucun acte durant la fen\u00eatre temporelle restreinte de l&rsquo;\u00e9tude par simple fluctuation stochastique.<\/li>\n<\/ul>\n<p>Pour mod\u00e9liser formellement cette double dynamique g\u00e9n\u00e9ratrice, Lambert (1992) a introduit les mod\u00e8les gonfl\u00e9s en z\u00e9ros : le mod\u00e8le <strong>ZIP<\/strong> (<em>Zero-Inflated Poisson<\/em>) et le mod\u00e8le <strong>ZINB<\/strong> (<em>Zero-Inflated Negative Binomial<\/em>). Ces mod\u00e8les superposent deux processus : une composante logistique binaire qui pr\u00e9dit l&rsquo;appartenance \u00e0 la classe des z\u00e9ros structurels, et une composante de comptage (Poisson ou Binomiale N\u00e9gative) r\u00e9gissant la fr\u00e9quence de survenue chez les individus \u00e9ligibles. Alternativement, les <strong>mod\u00e8les barri\u00e8re<\/strong> (<em>Hurdle Models<\/em>), th\u00e9oris\u00e9s par Mullahy (1986), d\u00e9couplent totalement le franchissement du premier \u00e9v\u00e9nement (mod\u00e8le binaire 0 vs &gt; 0) de la dynamique de r\u00e9currence conditionnelle aux valeurs strictement positives.<\/p>\n<p>La discrimination formelle entre un mod\u00e8le standard (Poisson ou NB2) et son homologue gonfl\u00e9 en z\u00e9ros (ZIP ou ZINB) ne peut s&rsquo;op\u00e9rer par un simple LRT classique, car les mod\u00e8les ne sont pas rigoureusement nich\u00e9s. On recourt alors traditionnellement au <strong>test non embo\u00eet\u00e9 de Vuong (1989)<\/strong>, compl\u00e9t\u00e9 par les ajustements m\u00e9thodologiques contemporains de Desmarais et Harden (2013). Une statistique de Vuong significativement positive valide la sup\u00e9riorit\u00e9 de l&rsquo;architecture avec inflation de z\u00e9ros.<\/p>\n<h3>9.2 Donn\u00e9es tronqu\u00e9es et censur\u00e9es<\/h3>\n<p>Un autre d\u00e9fi d&rsquo;\u00e9chantillonnage concerne les donn\u00e9es de d\u00e9nombrement tronqu\u00e9es \u00e0 z\u00e9ro (<em>Zero-Truncated Counts<\/em>). Ce sc\u00e9nario \u00e9merge d\u00e8s lors que le processus de recrutement d&rsquo;une \u00e9tude clinique conditionne l&rsquo;inclusion des patients \u00e0 la condition formelle d&rsquo;avoir manifest\u00e9 le comportement au moins une fois. Par exemple, une \u00e9tude examinant la dur\u00e9e de s\u00e9jour hospitalier (en nombre de jours entiers) chez des patients psychiatriques admis en unit\u00e9 ferm\u00e9e, ou le d\u00e9nombrement des consultations m\u00e9dicales chez des patients recrut\u00e9s en salle d&rsquo;attente, exclut structurellement la possibilit\u00e9 d&rsquo;observer la valeur z\u00e9ro.<\/p>\n<p>Ajuster un mod\u00e8le de Poisson ou binomial n\u00e9gatif standard \u00e0 de telles donn\u00e9es biaise lourdement l&rsquo;estimation de l&rsquo;esp\u00e9rance conditionnelle, car l&rsquo;espace d&rsquo;\u00e9chantillonnage r\u00e9el d\u00e9bute \u00e0 <em>Y<\/em> = 1. Il est alors requis d&#8217;employer des formulations distributionnelles tronqu\u00e9es : la <strong>Poisson tronqu\u00e9e \u00e0 z\u00e9ro (ZTP)<\/strong> ou la <strong>Binomiale n\u00e9gative tronqu\u00e9e \u00e0 z\u00e9ro (ZTNB)<\/strong>. La fonction de masse probabiliste est alors re-normalis\u00e9e en divisant la probabilit\u00e9 standard par la probabilit\u00e9 compl\u00e9mentaire de ne pas observer z\u00e9ro :<\/p>\n<p>P(Y = y | Y &gt; 0, &mu;) = P(Y = y | &mu;) \/ [1 &#8211; P(Y = 0 | &mu;)]<\/p>\n<p>De fa\u00e7on sym\u00e9trique, lorsque la fen\u00eatre temporelle d&rsquo;observation est interrompue arbitrairement avant que certains participants n&rsquo;aient termin\u00e9 leur s\u00e9quence comportementale (censure \u00e0 droite des comptages), le recours \u00e0 des extensions de mod\u00e8les de survie discrets ou \u00e0 des r\u00e9gressions de Poisson censur\u00e9es devient imp\u00e9ratif pour \u00e9viter des estimations artificiellement tir\u00e9es vers le bas.<\/p>\n<h3>9.3 Approches alternatives : Quasi-Poisson et erreurs-types robustes<\/h3>\n<p>Face \u00e0 la surdispersion, l&rsquo;analyste n&rsquo;est pas strictement restreint au mod\u00e8le binomial n\u00e9gatif param\u00e9trique standard. Deux alternatives m\u00e9thodologiques majeures m\u00e9ritent d&rsquo;\u00eatre consid\u00e9r\u00e9es dans l&rsquo;arsenal du statisticien : le mod\u00e8le Quasi-Poisson et l&rsquo;estimation de variance sandwich robuste.<\/p>\n<p>L&rsquo;approche du <strong>Quasi-Poisson<\/strong>, formul\u00e9e par Wedderburn (1974) dans le cadre de la th\u00e9orie des quasi-vraisemblances, refuse de postuler une distribution probabiliste enti\u00e8rement sp\u00e9cifi\u00e9e pour la composante al\u00e9atoire. Elle retient uniquement la fonction de lien logarithmique ln(&mu;) = X&beta;, mais mod\u00e9lise la variance empirique comme une fonction lin\u00e9aire de la moyenne proportionnelle \u00e0 un param\u00e8tre d&rsquo;\u00e9chelle libre &Phi; :<\/p>\n<p>Var(Y | X) = &Phi; &times; &mu;<\/p>\n<p>Dans ce cadre, les coefficients estim\u00e9s &beta;&#770; sont strictement identiques \u00e0 ceux obtenus par la r\u00e9gression de Poisson standard, mais leurs erreurs-types associ\u00e9es sont syst\u00e9matiquement corrig\u00e9es en \u00e9tant multipli\u00e9es par le facteur multiplicatif &radic;&Phi;&#770;. Le quasi-Poisson offre une flexibilit\u00e9 int\u00e9ressante lorsqu&rsquo;on h\u00e9site sur la nature exacte de la distribution de probabilit\u00e9, mais il poss\u00e8de une limite conceptuelle majeure : n&rsquo;\u00e9tant pas fond\u00e9 sur une v\u00e9ritable fonction de vraisemblance compl\u00e8te, il interdit le calcul de l&rsquo;AIC et du BIC, compliquant grandement les proc\u00e9dures de s\u00e9lection formelle de mod\u00e8les.<\/p>\n<p>La seconde alternative repose sur le conservatisme du mod\u00e8le de Poisson coupl\u00e9 \u00e0 une <strong>estimation de variance sandwich robuste de Huber-White<\/strong> (parfois d\u00e9sign\u00e9e sous le nom d&rsquo;estimateur sandwich ou de quasi-maximum de vraisemblance, QMLE). Comme l&rsquo;ont formalis\u00e9 Wooldridge (1999) et Cameron et Trivedi (2013), l&rsquo;estimateur sandwich calcule la variance asymptotique selon la formule matricielle :<\/p>\n<p>V<sub>sandwich<\/sub> = (X&rsquo; M X)<sup>-1<\/sup> [ &sum;<sub>i=1<\/sub><sup>n<\/sup> (y<sub>i<\/sub> &#8211; &mu;&#770;<sub>i<\/sub>)<sup>2<\/sup> X<sub>i<\/sub> X<sub>i<\/sub>&lsquo; ] (X&rsquo; M X)<sup>-1<\/sup><\/p>\n<p>o\u00f9 <em>M<\/em> est la matrice diagonale des variances pr\u00e9dites. Cette correction procure des erreurs-types asymptotiquement valides m\u00eame en pr\u00e9sence d&rsquo;une surdispersion arbitraire ou de certaines formes d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 non mod\u00e9lis\u00e9es. N\u00e9anmoins, pour des \u00e9chantillons cliniques de taille mod\u00e9r\u00e9e (N &lt; 200), la mod\u00e9lisation param\u00e9trique explicite via la r\u00e9gression binomiale n\u00e9gative demeure syst\u00e9matiquement plus efficace et offre une puissance statistique sup\u00e9rieure \u00e0 l&rsquo;approche sandwich de Poisson.<\/p>\n<h2>10. 10. Guide d\u00e9cisionnel s\u00e9quentiel pour la s\u00e9lection du mod\u00e8le<\/h2>\n<h3>10.1 \u00c9tape 1 : Analyse exploratoire approfondie des donn\u00e9es brutes<\/h3>\n<p>Toute d\u00e9marche de mod\u00e9lisation de d\u00e9nombrements doit imp\u00e9rativement d\u00e9buter par une analyse exploratoire minutieuse de la distribution univari\u00e9e de la variable d\u00e9pendante <em>Y<\/em>, ainsi que de ses relations bivari\u00e9es avec les pr\u00e9dicteurs cl\u00e9s. La premi\u00e8re op\u00e9ration quantitative fondamentale consiste \u00e0 calculer l&rsquo;<strong>indice de dispersion brut non ajust\u00e9<\/strong>, d\u00e9fini comme le ratio de la variance \u00e9chantillonnale sur la moyenne empirique :<\/p>\n<p>ID = s<sub>Y<\/sub><sup>2<\/sup> \/ Y&#773;<\/p>\n<p>Si la variable d\u00e9pendait rigoureusement d&rsquo;un processus de Poisson univari\u00e9, cet indice oscillerait \u00e9troitement autour de la valeur 1,0. Dans les cohortes de psychopathologie ou de neurosciences comportementales, il n&rsquo;est pas rare de constater un ratio brut exc\u00e9dant 5, 10 ou 20, fournissant un premier signal pr\u00e9liminaire massif de surdispersion, bien que cette dispersion puisse \u00eatre en partie r\u00e9sorb\u00e9e par l&rsquo;adjonction future de covariables pertinentes.<\/p>\n<p>Parall\u00e8lement, l&rsquo;analyste doit dresser l&rsquo;histogramme des comptages bruts pour scruter la forme distributionnelle globale. Trois \u00e9l\u00e9ments visuels critiques doivent \u00eatre inventori\u00e9s :<\/p>\n<ul>\n<li>La proportion exacte d&rsquo;observations nulles (P(Y = 0)). Si cette fraction exc\u00e8de 50 % ou 60 % de l&rsquo;\u00e9chantillon global, l&rsquo;\u00e9ventualit\u00e9 d&rsquo;un mod\u00e8le avec inflation de z\u00e9ros devra \u00eatre gard\u00e9e en m\u00e9moire pour les \u00e9tapes ult\u00e9rieures.<\/li>\n<li>La pr\u00e9sence de valeurs extr\u00eames isol\u00e9es dans la queue sup\u00e9rieure droite de la distribution. Ces points doivent \u00eatre audit\u00e9s afin de s&rsquo;assurer qu&rsquo;il ne s&rsquo;agit pas d&rsquo;erreurs d&rsquo;encodage de donn\u00e9es ou d&rsquo;artefacts exp\u00e9rimentaux n\u00e9cessitant une exclusion ou une analyse de sensibilit\u00e9.<\/li>\n<li>L&rsquo;identification formelle de toute variable d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 d&rsquo;exposition (temps pass\u00e9 en consultation, dur\u00e9e du suivi, nombre de questions pr\u00e9sent\u00e9es), imposant la construction imm\u00e9diate du terme d&rsquo;offset logarithmique ln(t<sub>i<\/sub>).<\/li>\n<\/ul>\n<h3>10.2 \u00c9tape 2 : Ajustement comparatif initial et diagnostic des r\u00e9sidus<\/h3>\n<p>La seconde phase consiste \u00e0 ajuster en parall\u00e8le le mod\u00e8le de r\u00e9gression de Poisson canonique de r\u00e9f\u00e9rence et le mod\u00e8le de r\u00e9gression binomiale n\u00e9gative (NB2) int\u00e9grant l&rsquo;ensemble th\u00e9orique des covariables <em>X<\/em> s\u00e9lectionn\u00e9es pour r\u00e9pondre \u00e0 la question scientifique de recherche. D\u00e8s la fin de l&rsquo;estimation de Poisson, le premier r\u00e9flexe m\u00e9thodologique ne doit pas \u00eatre d&rsquo;examiner la significativit\u00e9 des p-valeurs des coefficients, mais d&rsquo;extraire la matrice des r\u00e9sidus de Pearson standardis\u00e9s.<\/p>\n<p>L&rsquo;exp\u00e9rimentateur construit imm\u00e9diatement le nuage de points mettant en relation ces r\u00e9sidus de Pearson standardis\u00e9s avec les valeurs pr\u00e9dites par le mod\u00e8le de Poisson. Si le graphique r\u00e9v\u00e8le une diss\u00e9mination homog\u00e8ne et sym\u00e9trique \u00e0 l&rsquo;int\u00e9rieur de la bande normative [-2, +2], sans forme g\u00e9om\u00e9trique d&rsquo;entonnoir et sans valeur s&rsquo;\u00e9cartant au-del\u00e0 du corridor critique [-3, +3], l&rsquo;hypoth\u00e8se de travail de Poisson peut \u00eatre provisoirement retenue.<\/p>\n<p>En revanche, si l&rsquo;on constate qu&rsquo;une proportion substantielle de r\u00e9sidus (par exemple 8 %, 15 % ou plus) d\u00e9borde hors de l&rsquo;intervalle [-2, +2], ou que des r\u00e9sidus extr\u00eames atteignent des valeurs de +4 ou +8, la pr\u00e9sence d&rsquo;une surdispersion conditionnelle est graphiquement confirm\u00e9e. \u00c0 ce stade, la comparaison visuelle avec le trac\u00e9 des r\u00e9sidus issu du mod\u00e8le binomial n\u00e9gatif NB2 d\u00e9montre g\u00e9n\u00e9ralement une normalisation spectaculaire de la dispersion r\u00e9siduelle, attestant que la fonction de variance quadratique absorbe l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 exc\u00e9dentaire.<\/p>\n<h3>10.3 \u00c9tape 3 : Tests formels, arbitrage et s\u00e9lection finale<\/h3>\n<p>La phase conclusive valide la s\u00e9lection du mod\u00e8le optimal par le croisement convergent des tests inf\u00e9rentiels formels et des crit\u00e8res de p\u00e9nalisation th\u00e9orique de l&rsquo;information :<\/p>\n<ul>\n<li><strong>Ex\u00e9cution du test LRT sur la fronti\u00e8re :<\/strong> L&rsquo;analyste r\u00e9alise le test du rapport de vraisemblance entre le mod\u00e8le de Poisson et le mod\u00e8le binomial n\u00e9gatif. En appliquant la correction de Self-Liang (p-valeur divis\u00e9e par 2), si la statistique de test se r\u00e9v\u00e8le hautement significative (p &lt; 0,001), l&rsquo;hypoth\u00e8se nulle &alpha; = 0 est d\u00e9finitivement abandonn\u00e9e au profit du mod\u00e8le binomial n\u00e9gatif NB2.<\/li>\n<li><strong>Confrontation des m\u00e9triques AIC et BIC :<\/strong> L&rsquo;examen des crit\u00e8res d&rsquo;information doit corroborer ce choix. Une r\u00e9duction de l&rsquo;AIC sup\u00e9rieure \u00e0 10 points sous le mod\u00e8le binomial n\u00e9gatif ent\u00e9rine formellement sa sup\u00e9riorit\u00e9 parcimonieuse.<\/li>\n<li><strong>Arbitrage final sur l&rsquo;exc\u00e8s de z\u00e9ros :<\/strong> Si la proportion de z\u00e9ros au sein des donn\u00e9es \u00e9tait particuli\u00e8rement massive lors de la phase exploratoire, l&rsquo;exp\u00e9rimentateur compare le mod\u00e8le NB2 retenu \u00e0 sa formulation avec inflation de z\u00e9ros (ZINB) au moyen du test non embo\u00eet\u00e9 de Vuong. Si le test de Vuong n&rsquo;atteint pas le seuil de significativit\u00e9 statistique et que l&rsquo;AIC ne s&rsquo;am\u00e9liore pas de fa\u00e7on notable, le mod\u00e8le binomial n\u00e9gatif standard NB2 est proclam\u00e9 mod\u00e8le d\u00e9finitif.<\/li>\n<\/ul>\n<p>Cet arbre d\u00e9cisionnel unifi\u00e9 garantit une d\u00e9marche rigoureuse, reproductible et m\u00e9thodologiquement \u00e9tay\u00e9e, immunisant le chercheur contre les biais inf\u00e9rentiels cons\u00e9cutifs \u00e0 un choix arbitraire de mod\u00e9lisation.<\/p>\n<h2>11. 11. Impl\u00e9mentation pratique et interpr\u00e9tation logicielle (R et Python)<\/h2>\n<h3>11.1 Mise en \u0153uvre compl\u00e8te sous l&rsquo;environnement R<\/h3>\n<p>L&rsquo;\u00e9cosyst\u00e8me statistique du logiciel libre R offre l&rsquo;environnement computationnel le plus mature et le plus complet pour l&rsquo;ajustement et le diagnostic comparatif des mod\u00e8les de donn\u00e9es de comptage. Le mod\u00e8le de r\u00e9gression de Poisson s&rsquo;ajuste nativement via la fonction standard <code>glm()<\/code> en sp\u00e9cifiant l&rsquo;argument distributionnel <code>family = poisson(link = \"log\")<\/code>. L&rsquo;ajustement du mod\u00e8le binomial n\u00e9gatif quadratique (NB2) avec estimation conjointe par maximum de vraisemblance du param\u00e8tre de taille &theta; requiert quant \u00e0 lui l&#8217;emploi de la fonction optimis\u00e9e <code>glm.nb()<\/code> fournie par l&rsquo;incontournable package <code>MASS<\/code>.<\/p>\n<p>Le flux analytique de r\u00e9f\u00e9rence se d\u00e9ploie selon la syntaxe suivante :<\/p>\n<p>Pour le mod\u00e8le de base : <code>modele_poi &lt;- glm(crises ~ therapie + score_anxiete + offset(log(semaines_suivi)), data = donnees_cliniques, family = poisson(link = \"log\"))<\/code><\/p>\n<p>Pour le mod\u00e8le binomial n\u00e9gatif : <code>library(MASS)<br \/>modele_nb &lt;- glm.nb(crises ~ therapie + score_anxiete + offset(log(semaines_suivi)), data = donnees_cliniques)<\/code><\/p>\n<p>La comparaison formelle des ajustements par le test de rapport de vraisemblance corrig\u00e9 sur la fronti\u00e8re est ex\u00e9cut\u00e9e via la fonction <code>lrtest()<\/code> du package <code>lmtest<\/code>, ou par l&rsquo;appel g\u00e9n\u00e9rique <code>anova(modele_poi, modele_nb)<\/code>. Afin d&rsquo;extraire les diagnostics r\u00e9siduels bas\u00e9s sur les quantiles simul\u00e9s d&rsquo;une fiabilit\u00e9 absolue, l&rsquo;analyste sollicitera le package <code>DHARMa<\/code> :<\/p>\n<p><code>library(DHARMa)<br \/>simulation_res &lt;- simulateResiduals(fittedModel = modele_nb, n = 1000)<br \/>plot(simulation_res)<br \/>testDispersion(simulation_res)<\/code><\/p>\n<p>La fonction <code>testDispersion()<\/code> calcule une simulation non param\u00e9trique de l&rsquo;indice de dispersion r\u00e9siduel et d\u00e9livre une p-valeur empirique exacte attestant de l&rsquo;\u00e9radication compl\u00e8te de la surdispersion sous la param\u00e9trisation binomiale n\u00e9gative choisie.<\/p>\n<h3>11.2 Mise en \u0153uvre compl\u00e8te sous Python<\/h3>\n<p>Dans l&rsquo;\u00e9cosyst\u00e8me Python orient\u00e9 vers la science des donn\u00e9es, la biblioth\u00e8que de r\u00e9f\u00e9rence pour la mod\u00e9lisation statistique inf\u00e9rentielle est <a rel=\"noopener noreferrer\" href=\"https:\/\/www.statsmodels.org\/\">statsmodels<\/a>. Contrairement \u00e0 Scikit-learn, qui priorise la minimisation de fonctions de perte pour la pr\u00e9diction sans d\u00e9livrance d&rsquo;erreurs-types ou de p-valeurs, Statsmodels offre des tables statistiques exhaustives conformes aux standards acad\u00e9miques.<\/p>\n<p>L&rsquo;impl\u00e9mentation repose sur le module <code>statsmodels.formula.api<\/code>, qui autorise la d\u00e9claration de formules de r\u00e9gression inspir\u00e9es de la syntaxe de Wilkinson-Rogers :<\/p>\n<p><code>import statsmodels.formula.api as smf<br \/>import statsmodels.api as sm<\/p>\n<p># Ajustement du mod\u00e8le de r\u00e9gression de Poisson<br \/>mod_poisson = smf.poisson('crises ~ C(therapie) + score_anxiete', data=df, offset=np.log(df['semaines_suivi'])).fit()<br \/>print(mod_poisson.summary())<\/code><\/p>\n<p>Pour ajuster le mod\u00e8le binomial n\u00e9gatif NB2 sous Python, Statsmodels propose deux approches distinctes. La premi\u00e8re recourt \u00e0 la classe g\u00e9n\u00e9rique <code>smf.negativebinomial(formula, data)<\/code>, en sp\u00e9cifiant l&rsquo;argument <code>loglike_method='nb2'<\/code>. La seconde utilise le module GLM g\u00e9n\u00e9ralis\u00e9 :<\/p>\n<p><code># R\u00e9gression binomiale n\u00e9gative via Statsmodels GLM<br \/>mod_nb2 = smf.glm('crises ~ C(therapie) + score_anxiete', data=df, offset=np.log(df['semaines_suivi']), family=sm.families.NegativeBinomial(alpha=sm.families.links.nbinom2())).fit()<br \/>print(mod_nb2.summary())<\/code><\/p>\n<p>Pour \u00e9valuer formellement le test de surdispersion de Cameron et Trivedi sous Python, un script automatis\u00e9 extrait les valeurs pr\u00e9dites &mu;&#770; et calcule la r\u00e9gression lin\u00e9aire des r\u00e9sidus transform\u00e9s :<\/p>\n<p><code>y_hat = mod_poisson.predict()<br \/>y_val = df['crises']<br \/>res_aux = ((y_val - y_hat)**2 - y_val) \/ y_hat<br \/>reg_aux = sm.OLS(res_aux, y_hat).fit()<br \/>print('P-valeur du test auxiliaire de surdispersion :', reg_aux.pvalues[0])<\/code><\/p>\n<p>L&rsquo;obtention d&rsquo;une p-valeur unilat\u00e9rale hautement significative confirme l&rsquo;imp\u00e9ratif computationnel de basculer vers le mod\u00e8le <code>NegativeBinomial<\/code>.<\/p>\n<h3>11.3 Lecture critique des sorties console et tableaux statistiques<\/h3>\n<p>L&rsquo;analyse d&rsquo;une sortie console produite par R ou Statsmodels requiert une grille de lecture m\u00e9thodique. Prenons l&rsquo;exemple d&rsquo;un tableau d&rsquo;estimation issu d&rsquo;un mod\u00e8le binomial n\u00e9gatif mod\u00e9lisant le nombre d&rsquo;incidents comportementaux :<\/p>\n<ul>\n<li><strong>D\u00e9viance r\u00e9siduelle et Degr\u00e9s de libert\u00e9 :<\/strong> Une lecture rapide compare la <em>Residual Deviance<\/em> au nombre de <em>Degrees of Freedom<\/em>. Si le ratio s&rsquo;\u00e9tablit \u00e0 proximit\u00e9 de 1,00 \u00e0 1,15, cela signale un \u00e9talonnage harmonieux de la variance conditionnelle.<\/li>\n<li><strong>Log-Likelihood :<\/strong> La log-vraisemblance maximale atteinte par le mod\u00e8le binomial n\u00e9gatif doit afficher une valeur alg\u00e9brique substantiellement moins n\u00e9gative que celle du mod\u00e8le de Poisson (par exemple, -412,5 contre -580,2), traduisant un gain d&rsquo;ajustement consid\u00e9rable.<\/li>\n<li><strong>Le param\u00e8tre d&rsquo;\u00e9chelle (Theta ou Alpha) :<\/strong> Dans la sortie de <code>glm.nb()<\/code> sous R, la derni\u00e8re ligne affiche : <code>Theta: 2.145 (Std. Err.: 0.284)<\/code>. La valeur de surdispersion \u00e9quivalente est &alpha;&#770; = 1 \/ 2,145 = 0,466. L&rsquo;erreur-type associ\u00e9e d\u00e9montre que ce param\u00e8tre est born\u00e9 tr\u00e8s loin de l&rsquo;infini, attestant d&rsquo;une surdispersion structurelle robuste.<\/li>\n<li><strong>Lecture des pentes &beta; et conversion en IRR :<\/strong> Si le coefficient brut associ\u00e9 \u00e0 la covariable <code>therapie_TCC<\/code> affiche <code>Estimate = -0.446, Std. Error = 0.112, z = -3.98, p = 6.8e-05<\/code>, l&rsquo;analyste proc\u00e8de imm\u00e9diatement \u00e0 l&rsquo;exponentiation : IRR = exp(-0,446) = 0,640. L&rsquo;intervalle de confiance \u00e0 95 % s&rsquo;\u00e9tablit \u00e0 [exp(-0,446 &#8211; 1,96 &times; 0,112) ; exp(-0,446 + 1,96 &times; 0,112)] = [0,514 ; 0,797]. La conclusion clinique s&rsquo;\u00e9nonce ainsi : \u00ab Toutes choses \u00e9gales par ailleurs, les patients ayant b\u00e9n\u00e9fici\u00e9 de la th\u00e9rapie cognitivo-comportementale manifestent une r\u00e9duction relative statistiquement significative de 36 % du taux attendu d&rsquo;incidents par rapport au groupe t\u00e9moin (IRR = 0,64, IC 95 % [0,51, 0,80], p &lt; 0,001) \u00bb.<\/li>\n<\/ul>\n<h2>12. 12. Bonnes pratiques de reporting acad\u00e9mique et consid\u00e9rations avanc\u00e9es<\/h2>\n<h3>12.1 Normes de publication scientifique (APA 7e \u00e9dition et \u00e9quivalents)<\/h3>\n<p>La publication de mod\u00e8les de r\u00e9gression pour donn\u00e9es de comptage dans les revues internationales de premier rang (telles que celles affili\u00e9es \u00e0 l&rsquo;American Psychological Association, APA) impose un niveau \u00e9lev\u00e9 de transparence m\u00e9thodologique. Trop d&rsquo;articles omettent encore de justifier le choix de leur fonction de variance, se contentant d&rsquo;\u00e9voquer l&#8217;emploi vague d&rsquo;un \u00ab mod\u00e8le de Poisson \u00bb. Les recommandations actuelles (normes APA 7e \u00e9dition et lignes directrices STROBE\/CONSORT) exigent une rigueur formalis\u00e9e d\u00e8s la section M\u00e9thode :<\/p>\n<ul>\n<li><strong>Justification th\u00e9orique et empirique pr\u00e9alable :<\/strong> L&rsquo;auteur doit expliciter pourquoi la r\u00e9gression lin\u00e9aire standard OLS a \u00e9t\u00e9 \u00e9cart\u00e9e (nature discr\u00e8te, borne en z\u00e9ro, asym\u00e9trie marqu\u00e9e) et d\u00e9crire la proc\u00e9dure statistique pr\u00e9cise employ\u00e9e pour arbitrer entre la distribution de Poisson et la distribution binomiale n\u00e9gative.<\/li>\n<li><strong>Restitution des m\u00e9triques diagnostiques de dispersion :<\/strong> La section des R\u00e9sultats doit rapporter textuellement le test de surdispersion r\u00e9alis\u00e9, la valeur estim\u00e9e du param\u00e8tre d&rsquo;\u00e9chelle &alpha; (ou &theta;) avec son intervalle de confiance \u00e0 95 %, ainsi que la comparaison des crit\u00e8res d&rsquo;information (&Delta;AIC). Une formulation canonique recommand\u00e9e est la suivante : \u00ab L&rsquo;analyse pr\u00e9liminaire a mis en \u00e9vidence une surdispersion conditionnelle substantielle (&chi;<sup>2<\/sup>\/df = 3,42 ; test du rapport de vraisemblance corrig\u00e9 sur la fronti\u00e8re : &chi;<sup>2<\/sup><sub>0:1<\/sub> = 148,6, p &lt; 0,001), justifiant le rejet du mod\u00e8le de Poisson au profit d&rsquo;une r\u00e9gression binomiale n\u00e9gative quadratique (NB2 ; param\u00e8tre d&rsquo;\u00e9chelle &theta; = 2,15, IC 95 % [1,64, 2,79] ; &Delta;AIC = 146,6) \u00bb.<\/li>\n<li><strong>Pr\u00e9sentation tabulaire normalis\u00e9e :<\/strong> Les tableaux de r\u00e9gression doivent comporter de fa\u00e7on syst\u00e9matique : les coefficients de r\u00e9gression bruts non standardis\u00e9s (<em>B<\/em>), les erreurs-types associ\u00e9es (<em>SE<\/em>), les Ratios de Taux d&rsquo;Incidence (<em>IRR<\/em>), leurs intervalles de confiance \u00e0 95 % exacts, et les p-valeurs pr\u00e9cises (sans recourir \u00e0 la simple notation binaire p &lt; 0,05). La taille de l&rsquo;\u00e9chantillon <em>N<\/em>, la log-vraisemblance maximale et la d\u00e9viance r\u00e9siduelle doivent imp\u00e9rativement figurer en note de bas de tableau.<\/li>\n<\/ul>\n<h3>12.2 Extension aux mod\u00e8les mixtes et hi\u00e9rarchiques (Multilevel Count Models)<\/h3>\n<p>Dans la recherche contemporaine en psychologie cognitive, exp\u00e9rimentale ou sociale, les protocoles empiriques impliquent presque constamment des mesures r\u00e9p\u00e9t\u00e9es dans le temps (\u00e9tudes longitudinales, m\u00e9thodes d&rsquo;\u00e9chantillonnage des exp\u00e9riences de vie ou <em>Ecological Momentary Assessment<\/em>) ou des structures d&rsquo;\u00e9chantillonnage en grappes (\u00e9l\u00e8ves imbriqu\u00e9s dans des classes, patients nich\u00e9s au sein de th\u00e9rapeutes). Dans de telles configurations, les observations individuelles violent l&rsquo;hypoth\u00e8se d&rsquo;ind\u00e9pendance stochastique r\u00e9siduelle.<\/p>\n<p>Le cadre th\u00e9orique d\u00e9velopp\u00e9 dans cet article s&rsquo;\u00e9tend alors naturellement aux <strong>Mod\u00e8les Lin\u00e9aires G\u00e9n\u00e9ralis\u00e9s Mixtes (GLMM)<\/strong>. La question de l&rsquo;arbitrage entre Poisson et Binomiale N\u00e9gative se complexifie, car il convient d\u00e9sormais de distinguer conceptuellement deux sources distinctes de variabilit\u00e9 exc\u00e9dentaire :<\/p>\n<ul>\n<li>La surdispersion induite par la structure de regroupement (<em>cluster-level variance<\/em>), qui peut \u00eatre \u00e9l\u00e9gamment captur\u00e9e par l&rsquo;introduction d&rsquo;un effet al\u00e9atoire gaussien sur l&rsquo;intercept au niveau du groupe : ln(&mu;<sub>ij<\/sub>) = X<sub>ij<\/sub>&beta; + u<sub>0j<\/sub>, avec u<sub>0j<\/sub> ~ N(0, &sigma;<sub>u<\/sub><sup>2<\/sup>).<\/li>\n<li>La surdispersion observationnelle pure (<em>observation-level overdispersion<\/em>), persistant au niveau individuel le plus fin m\u00eame apr\u00e8s conditionnement sur les effets al\u00e9atoires de grappe.<\/li>\n<\/ul>\n<p>Sous l&rsquo;environnement R, le package sp\u00e9cialis\u00e9 <code>glmmTMB<\/code> (<a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.32614\/RJ-2017-066\">Brooks et al., 2017<\/a>) s&rsquo;est impos\u00e9 comme la r\u00e9f\u00e9rence absolue pour estimer ces mod\u00e8les complexes, surpassant largement les limites computationnelles de la fonction <code>glmer()<\/code> du package <code>lme4<\/code>. <code>glmmTMB<\/code> autorise l&rsquo;ajustement sans faille de mod\u00e8les hi\u00e9rarchiques longitudinaux de Poisson, de quasi-Poisson, de binomiale n\u00e9gative (NB1 et NB2), ainsi que leurs extensions avec inflation de z\u00e9ros, tout en mod\u00e9lisant des structures de corr\u00e9lation temporelle r\u00e9siduelle de type autor\u00e9gressif d&rsquo;ordre 1 (AR1).<\/p>\n<h3>12.3 Synth\u00e8se \u00e9pist\u00e9mologique pour la recherche clinique<\/h3>\n<p>Le choix d&rsquo;une architecture de mod\u00e9lisation statistique ne constitue en aucun cas un d\u00e9tail technique subalterne ou un simple choix de commodit\u00e9 calculatoire. Il engage directement la responsabilit\u00e9 \u00e9thique et soci\u00e9tale du chercheur. Dans le champ de la sant\u00e9 mentale et des neurosciences comportementales, valider \u00e0 tort l&rsquo;efficacit\u00e9 clinique d&rsquo;une mol\u00e9cule pharmacologique ou d&rsquo;une intervention psychoth\u00e9rapeutique sur la foi de faux positifs engendr\u00e9s par un mod\u00e8le de Poisson surdispers\u00e9 peut orienter inutilement des ann\u00e9es de recherche translationnelle et exposer des patients r\u00e9els \u00e0 des th\u00e9rapeutiques inefficaces.<\/p>\n<p>Il est donc temps d&rsquo;abandonner d\u00e9finitivement les pratiques m\u00e9thodologiques obsol\u00e8tes que constituent la r\u00e9gression lin\u00e9aire ordinaire sur comptages bruts et les transformations artificielles par le logarithme de Y + 1. La ma\u00eetrise rigoureuse des mod\u00e8les lin\u00e9aires g\u00e9n\u00e9ralis\u00e9s, articul\u00e9e autour de la distinction fondamentale entre \u00e9quidispersion de Poisson et flexibilit\u00e9 quadratique binomiale n\u00e9gative, constitue un pilier m\u00e9thodologique indispensable \u00e0 tout chercheur soucieux de rigueur inf\u00e9rentielle, de robustesse empirique et de reproductibilit\u00e9 scientifique.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Brooks, M. E., Kristensen, K., van Benthem, K. J., Magnusson, A., Berg, C. W., Nielsen, A., Skaug, H. J., M\u00e4chler, M., &amp; Bolker, B. M. (2017). glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. <em>The R Journal<\/em>, 9(2), 378\u2013400. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.32614\/RJ-2017-066\">https:\/\/doi.org\/10.32614\/RJ-2017-066<\/a><\/li>\n<li>Burnham, K. P., &amp; Anderson, D. R. (2002). <em>Model selection and multimodel inference: A practical information-theoretic approach<\/em> (2e \u00e9d.). Springer-Verlag. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/b97636\">https:\/\/doi.org\/10.1007\/b97636<\/a><\/li>\n<li>Cameron, A. C., &amp; Trivedi, P. K. (2013). <em>Regression analysis of count data<\/em> (2e \u00e9d.). Cambridge University Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1017\/CBO9781139016650\">https:\/\/doi.org\/10.1017\/CBO9781139016650<\/a><\/li>\n<li>Chernoff, H. (1954). On the distribution of the likelihood ratio test of one-sided hypotheses associated with a normally distributed variable. <em>The Annals of Mathematical Statistics<\/em>, 25(3), 573\u2013578. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aoms\/1177728723\">https:\/\/doi.org\/10.1214\/aoms\/1177728723<\/a><\/li>\n<li>Desmarais, B. A., &amp; Harden, J. J. (2013). Testing for zero inflation in count models: Bias correction for the Vuong test. <em>The Stata Journal<\/em>, 13(4), 810\u2013835. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1177\/1536867X1301300408\">https:\/\/doi.org\/10.1177\/1536867X1301300408<\/a><\/li>\n<li>Gouri\u00e9roux, C., Monfort, A., &amp; Trognon, A. (1984). Pseudo maximum likelihood methods: Applications to Poisson models. <em>Econometrica<\/em>, 52(3), 701\u2013720. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1913472\">https:\/\/doi.org\/10.2307\/1913472<\/a><\/li>\n<li>Kleiber, C., &amp; Zeileis, A. (2016). Visualizing count data regressions using rootograms. <em>The American Statistician<\/em>, 70(3), 296\u2013303. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/00031305.2016.1173590\">https:\/\/doi.org\/10.1080\/00031305.2016.1173590<\/a><\/li>\n<li>Lambert, D. (1992). Zero-inflated Poisson regression, with an application to defects in manufacturing. <em>Technometrics<\/em>, 34(1), 1\u201314. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1269547\">https:\/\/doi.org\/10.2307\/1269547<\/a><\/li>\n<li>Mullahy, J. (1986). Specification and testing of some modified count data models. <em>Journal of Econometrics<\/em>, 33(3), 341\u2013365. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/0304-4076(86)90002-3\">https:\/\/doi.org\/10.1016\/0304-4076(86)90002-3<\/a><\/li>\n<li>Nelder, J. A., &amp; Wedderburn, R. W. M. (1972). Generalized linear models. <em>Journal of the Royal Statistical Society: Series A (General)<\/em>, 135(3), 370\u2013384. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2344614\">https:\/\/doi.org\/10.2307\/2344614<\/a><\/li>\n<li>O&rsquo;Hara, R. B., &amp; Kotze, D. J. (2010). Do not log-transform count data. <em>Methods in Ecology and Evolution<\/em>, 1(2), 118\u2013122. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2041-210X.2010.00021.x\">https:\/\/doi.org\/10.1111\/j.2041-210X.2010.00021.x<\/a><\/li>\n<li>Self, S. G., &amp; Liang, K.-Y. (1987). Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. <em>Journal of the American Statistical Association<\/em>, 82(398), 605\u2013610. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1987.10478472\">https:\/\/doi.org\/10.1080\/01621459.1987.10478472<\/a><\/li>\n<li>Vuong, Q. H. (1989). Likelihood ratio tests for model selection and non-nested hypotheses. <em>Econometrica<\/em>, 57(2), 307\u2013333. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1912557\">https:\/\/doi.org\/10.2307\/1912557<\/a><\/li>\n<li>Wedderburn, R. W. M. (1974). Quasi-likelihood functions, generalized linear models, and the Gauss-Newton method. <em>Biometrika<\/em>, 61(3), 439\u2013447. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/61.3.439\">https:\/\/doi.org\/10.1093\/biomet\/61.3.439<\/a><\/li>\n<li>Wooldridge, J. M. (1999). Distribution-free estimation of some nonlinear panel data models. <em>Journal of Econometrics<\/em>, 90(1), 77\u201397. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/S0304-4076(98)00033-5\">https:\/\/doi.org\/10.1016\/S0304-4076(98)00033-5<\/a><\/li>\n<li>Zeileis, A., Kleiber, C., &amp; Jackman, S. (2008). Regression models for count data in R. <em>Journal of Statistical Software<\/em>, 27(8), 1\u201325. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.18637\/jss.v027.i08\">https:\/\/doi.org\/10.18637\/jss.v027.i08<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide m\u00e9thodologique complet pour choisir entre la r\u00e9gression de Poisson et binomiale n\u00e9gative face aux donn\u00e9es de comptage en recherche psychologique.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2169],"tags":[5481,7422,8883,4990,5354],"class_list":["post-7504","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-psychologie-quantitative","tag-donnees-de-comptage","tag-modeles-lineaires-generalises","tag-regression-binomiale-negative","tag-regression-de-poisson","tag-surdispersion"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7504","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=7504"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7504\/revisions"}],"predecessor-version":[{"id":7506,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7504\/revisions\/7506"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=7504"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=7504"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=7504"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}