{"id":5907,"date":"2026-09-16T03:29:02","date_gmt":"2026-09-16T03:29:02","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/qu-est-ce-qu-une-distribution-de-probabilite-conjointe\/"},"modified":"2026-09-16T03:29:04","modified_gmt":"2026-09-16T03:29:04","slug":"qu-est-ce-qu-une-distribution-de-probabilite-conjointe","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=qu-est-ce-qu-une-distribution-de-probabilite-conjointe","title":{"rendered":"Qu&rsquo;est-ce qu&rsquo;une distribution de probabilit\u00e9 conjointe ?"},"content":{"rendered":"<p>Dans la gen\u00e8se des sciences contemporaines, la quantification de l&rsquo;incertitude s&rsquo;est heurt\u00e9e \u00e0 une r\u00e9alit\u00e9 fondamentale : aucun ph\u00e9nom\u00e8ne complexe n&rsquo;\u00e9volue de mani\u00e8re strictement isol\u00e9e. Que l&rsquo;on s&rsquo;int\u00e9resse \u00e0 la neurobiologie comportementale, aux fluctuations macro\u00e9conomiques ou \u00e0 la m\u00e9canique quantique, l&rsquo;illusion d&rsquo;une variable al\u00e9atoire solitaire gouvernant un m\u00e9canisme unitaire s&rsquo;effondre invariablement face \u00e0 l&rsquo;enchev\u00eatrement des causalit\u00e9s et des concomitances. La <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Distribution_de_probabilit%C3%A9\">distribution de probabilit\u00e9<\/a> univari\u00e9e, si utile pour d\u00e9crire la dispersion d&rsquo;un attribut singulier, se r\u00e9v\u00e8le incapable de saisir les dynamiques relationnelles, les covariations et les d\u00e9pendances mutuelles qui constituent l&rsquo;essence m\u00eame de l&rsquo;architecture du monde r\u00e9el.<\/p>\n<p>Pour d\u00e9passer cette r\u00e9duction \u00e9pist\u00e9mologique, le calcul des probabilit\u00e9s a d\u00e9velopp\u00e9 le cadre math\u00e9matique de la <strong>distribution de probabilit\u00e9 conjointe<\/strong>. Cet \u00e9difice conceptuel ne se contente pas d&rsquo;aligner en parall\u00e8le plusieurs distributions marginales ind\u00e9pendantes ; il \u00e9rige un espace g\u00e9om\u00e9trique et probabiliste multidimensionnel capable de cartographier la structure intime de d\u00e9pendance liant simultan\u00e9ment deux ou plusieurs variables al\u00e9atoires. En attribuant une mesure de vraisemblance \u00e0 chaque combinaison concevable d&rsquo;\u00e9tats, la loi conjointe devient le socle indispensable \u00e0 toute inf\u00e9rence pr\u00e9dictive, \u00e0 toute mod\u00e9lisation causale et \u00e0 toute d\u00e9composition structurelle des syst\u00e8mes complexes.<\/p>\n<p>Le pr\u00e9sent trait\u00e9 propose une exploration exhaustive de la distribution conjointe, guidant le lecteur depuis ses fondations intuitives dans les tableaux de contingence jusqu&rsquo;aux th\u00e9or\u00e8mes asymptotiques les plus sophistiqu\u00e9s de la mod\u00e9lisation en haute dimension. \u00c0 travers un dialogue constant entre formalisation math\u00e9matique rigoureuse et applications m\u00e9thodologiques issues notamment de la psychom\u00e9trie, des sciences cognitives et de la science des donn\u00e9es, nous exposerons la m\u00e9canique qui permet de condenser, d&rsquo;analyser et d&rsquo;exploiter la totalit\u00e9 de l&rsquo;information contenue dans les al\u00e9as multiples.<\/p>\n<h2>1. Introduction \u00e9pist\u00e9mologique et conceptuelle \u00e0 la distribution conjointe<\/h2>\n<h3>1.1 Origine et n\u00e9cessit\u00e9 de l&rsquo;analyse multivari\u00e9e simultan\u00e9e<\/h3>\n<p>L&rsquo;histoire de la pens\u00e9e statistique est indissociable du passage progressif d&rsquo;une vision r\u00e9ductionniste et unidimensionnelle \u00e0 une appr\u00e9hension globale et multivari\u00e9e des ph\u00e9nom\u00e8nes. Aux pr\u00e9mices de la discipline, l&rsquo;attention des observateurs se focalisait sur des grandeurs isol\u00e9es : la taille moyenne d&rsquo;une population, le taux de natalit\u00e9 annuel ou la vitesse de s\u00e9dimentation d&rsquo;un compos\u00e9 chimique. Cependant, ces distributions univari\u00e9es manifestent rapidement des limites structurelles insurmontables d\u00e8s lors que l&rsquo;on cherche \u00e0 expliquer, et non plus simplement \u00e0 d\u00e9crire. En isolant une variable de son contexte \u00e9cologique, l&rsquo;analyste commet une amputation informationnelle : il ignore d\u00e9lib\u00e9r\u00e9ment les forces d&rsquo;entra\u00eenement r\u00e9ciproques, les r\u00e9troactions et les r\u00e9gulations syst\u00e9miques qui gouvernent l&rsquo;ensemble \u00e9tudi\u00e9.<\/p>\n<p>La transition conceptuelle de l&rsquo;observation isol\u00e9e \u00e0 la covariation de variables multiples s&rsquo;est op\u00e9r\u00e9e au cours du XIXe si\u00e8cle, sous l&rsquo;impulsion de pionniers comme Francis Galton et Karl Pearson. Confront\u00e9s \u00e0 l&rsquo;\u00e9tude des caract\u00e8res biologiques h\u00e9r\u00e9ditaires, ces chercheurs comprirent que la distribution des statures des enfants ne pouvait \u00eatre appr\u00e9hend\u00e9e sans y associer simultan\u00e9ment celle des parents. Consid\u00e9rer s\u00e9par\u00e9ment la distribution univari\u00e9e de la taille des g\u00e9niteurs et celle des descendants dissimulait le ph\u00e9nom\u00e8ne essentiel : le degr\u00e9 de transmission du trait d&rsquo;une g\u00e9n\u00e9ration \u00e0 l&rsquo;autre. Il devint alors imp\u00e9ratif d&rsquo;\u00e9laborer un formalisme math\u00e9matique capable d&rsquo;assigner une probabilit\u00e9 non pas \u00e0 une valeur isol\u00e9e, mais \u00e0 un vecteur d&rsquo;observations conjointes, capturant ainsi la courbure et l&rsquo;intensit\u00e9 des liens probabilistes unissant les composantes du syst\u00e8me.<\/p>\n<p>Dans la mod\u00e9lisation quantitative moderne, la prise en compte simultan\u00e9e des al\u00e9as est devenue un pr\u00e9requis m\u00e9thodologique fondamental. Lorsqu&rsquo;un algorithme \u00e9value le risque de cr\u00e9dit d&rsquo;un individu, qu&rsquo;un mod\u00e8le m\u00e9t\u00e9orologique pr\u00e9voit la trajectoire d&rsquo;une d\u00e9pression ou qu&rsquo;un psychologue cognitif analyse l&rsquo;interaction entre m\u00e9moire de travail et raisonnement abstrait, l&rsquo;objet math\u00e9matique sous-jacent n&rsquo;est jamais une succession de lois univari\u00e9es arbitrairement juxtapos\u00e9es. Il s&rsquo;agit d&rsquo;une fonction de distribution conjointe, qui prend en charge la totalit\u00e9 de la structure d&rsquo;interd\u00e9pendance. Sans cette approche multivari\u00e9e, la statistique inf\u00e9rentielle moderne perdrait sa puissance pr\u00e9dictive et son aptitude \u00e0 isoler les structures sous-jacentes \u00e0 travers des m\u00e9thodes avanc\u00e9es telles que les mod\u00e8les d&rsquo;\u00e9quations structurelles, l&rsquo;analyse factorielle ou les r\u00e9seaux de neurones profonds probabilistes.<\/p>\n<h3>1.2 D\u00e9finition formelle intuitive de l&rsquo;\u00e9v\u00e9nement conjoint<\/h3>\n<p>Sur le plan intuitif, un \u00e9v\u00e9nement probabiliste classique renvoie \u00e0 la r\u00e9alisation d&rsquo;un \u00e9tat particulier au sein d&rsquo;un univers des possibles &Omega;. Lorsque l&rsquo;on s&rsquo;\u00e9l\u00e8ve au rang de l&rsquo;analyse multivari\u00e9e, la notion centrale devient celle d&rsquo;<strong>\u00e9v\u00e9nement conjoint<\/strong>. Consid\u00e9rons deux espaces probabilisables associ\u00e9s respectivement aux variables al\u00e9atoires <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>. Un \u00e9v\u00e9nement conjoint d\u00e9signe la survenue simultan\u00e9e, concomitante et indissociable de la condition <span class=\"arabpsy-math-inline\" data-latex=\"X = x\">$X = x$<\/span> et de la condition <span class=\"arabpsy-math-inline\" data-latex=\"Y = y\">$Y = y$<\/span>. G\u00e9om\u00e9triquement et ensemblistement, cet \u00e9tat correspond pr\u00e9cis\u00e9ment \u00e0 l&rsquo;intersection logique de deux sous-ensembles au sein de la tribu des \u00e9v\u00e9nements : <span class=\"arabpsy-math-inline\" data-latex=\"A cap B\">$A cap B$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"A = {\\omega in \\Omega mid X(\\omega) = x}\">$A = {\\omega in \\Omega mid X(\\omega) = x}$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"B = {\\omega in \\Omega mid Y(\\omega) = y}\">$B = {\\omega in \\Omega mid Y(\\omega) = y}$<\/span>.<\/p>\n<p>Le qualificatif de \u00ab conjoint \u00bb requiert une attention s\u00e9mantique particuli\u00e8re. Il ne traduit pas une simple cons\u00e9cution temporelle ou une causalit\u00e9 directionnelle pr\u00e9\u00e9tablie, mais une co-occurrence dans l&rsquo;espace des phases mesurable. Il convient d&rsquo;\u00e9tablir une distinction \u00e9pist\u00e9mologique stricte entre des \u00e9v\u00e9nements mutuellement exclusifs et des \u00e9v\u00e9nements concomitants. Deux \u00e9v\u00e9nements sont disjoints ou mutuellement exclusifs lorsque leur intersection est vide (<span class=\"arabpsy-math-inline\" data-latex=\"A cap B = emptyset\">$A cap B = emptyset$<\/span>), ce qui rend leur r\u00e9alisation conjointe rigoureusement impossible (<span class=\"arabpsy-math-inline\" data-latex=\"P(A cap B) = 0\">$P(A cap B) = 0$<\/span>). \u00c0 l&rsquo;inverse, l&rsquo;analyse conjointe cible des \u00e9v\u00e9nements compossibles, dont la r\u00e9alisation synchrone r\u00e9v\u00e8le la texture informationnelle liant les attributs observ\u00e9s.<\/p>\n<p>Des illustrations \u00e9l\u00e9mentaires issues de l&rsquo;observation empirique et de la psychom\u00e9trie \u00e9clairent imm\u00e9diatement cette construction. Imaginons un protocole d&rsquo;\u00e9valuation neuropsychologique o\u00f9 un sujet est soumis \u00e0 un test d&rsquo;inhibition cognitive. La performance peut \u00eatre caract\u00e9ris\u00e9e par deux variables : le succ\u00e8s ou l&rsquo;\u00e9chec \u00e0 la t\u00e2che (variable binaire <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span>) et la manifestation d&rsquo;un signe d&rsquo;anxi\u00e9t\u00e9 neurov\u00e9g\u00e9tative mesur\u00e9 par la r\u00e9ponse \u00e9lectrodermale (variable dichotomis\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>). L&rsquo;\u00e9v\u00e9nement isol\u00e9 \u00ab le sujet \u00e9choue \u00bb n&rsquo;offre qu&rsquo;un aper\u00e7u partiel du fonctionnement cognitif ; de m\u00eame, l&rsquo;\u00e9v\u00e9nement isol\u00e9 \u00ab le sujet pr\u00e9sente une activation \u00e9lectrodermale \u00e9lev\u00e9e \u00bb reste ambigu sur le plan clinique. C&rsquo;est l&rsquo;\u00e9v\u00e9nement conjoint \u2014 \u00ab le sujet \u00e9choue ET pr\u00e9sente une r\u00e9ponse \u00e9lectrodermale \u00e9lev\u00e9e \u00bb \u2014 qui constitue l&rsquo;unit\u00e9 d&rsquo;analyse signifiante. La distribution de probabilit\u00e9 conjointe a pour r\u00f4le d&rsquo;affecter un poids stochastique \u00e0 chacune de ces configurations synchrones pour dresser la cartographie exhaustive du comportement \u00e9tudi\u00e9.<\/p>\n<h2>2. Du tableau de contingence \u00e0 la loi de probabilit\u00e9 discr\u00e8te<\/h2>\n<h3>2.1 Structure et lecture d&rsquo;un tableau crois\u00e9 \u00e0 deux entr\u00e9es<\/h3>\n<p>La porte d&rsquo;entr\u00e9e la plus concr\u00e8te vers la formalisation d&rsquo;une loi conjointe discr\u00e8te r\u00e9side dans le tableau de contingence, ou tableau crois\u00e9 \u00e0 deux entr\u00e9es. Imagin\u00e9 pour consigner les observations empiriques de deux variables qualitatives ou discr\u00e8tes, ce dispositif organise les donn\u00e9es selon une architecture matricielle bidimensionnelle. Les lignes du tableau correspondent aux diff\u00e9rentes modalit\u00e9s ou valeurs r\u00e9alisables par la premi\u00e8re variable <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span>, indic\u00e9es par <span class=\"arabpsy-math-inline\" data-latex=\"i in {1, dots, r}\">$i in {1, dots, r}$<\/span>, tandis que les colonnes mat\u00e9rialisent les modalit\u00e9s de la seconde variable <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>, indic\u00e9es par <span class=\"arabpsy-math-inline\" data-latex=\"j in {1, dots, c}\">$j in {1, dots, c}$<\/span>.<\/p>\n<p>Chaque cellule interne \u00e0 l&rsquo;intersection de la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-i\u00e8me ligne et de la <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span>-i\u00e8me colonne h\u00e9berge un effectif conjoint observ\u00e9, conventionnellement not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"n_{ij}\">$n_{ij}$<\/span>. Ce nombre entier repr\u00e9sente la quantit\u00e9 exacte d&rsquo;individus ou d&rsquo;unit\u00e9s exp\u00e9rimentales qui manifestent simultan\u00e9ment la modalit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"x_i\">$x_i$<\/span> et la modalit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"y_j\">$y_j$<\/span>. La rigueur m\u00e9thodologique impose de tracer une fronti\u00e8re conceptuelle infranchissable entre ces effectifs observ\u00e9s <span class=\"arabpsy-math-inline\" data-latex=\"n_{ij}\">$n_{ij}$<\/span>, contingents \u00e0 l&rsquo;\u00e9chantillon pr\u00e9lev\u00e9 et soumis aux al\u00e9as d&rsquo;\u00e9chantillonnage, et les probabilit\u00e9s th\u00e9oriques sous-jacentes <span class=\"arabpsy-math-inline\" data-latex=\"p_{ij} = P(X = x_i, Y = y_j)\">$p_{ij} = P(X = x_i, Y = y_j)$<\/span>, qui gouvernent le m\u00e9canisme g\u00e9n\u00e9rateur de donn\u00e9es au sein de la population cible globale.<\/p>\n<p>Prenons un exemple paradigmatique emprunt\u00e9 \u00e0 la psychologie sociale du consommateur : croisons le style d\u00e9cisionnel d&rsquo;un groupe d&rsquo;individus (<span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span>, prenant les modalit\u00e9s {Rationnel, Intuitif}) avec leur modalit\u00e9 d&rsquo;engagement envers une nouvelle technologie (<span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>, prenant les valeurs {Adoption imm\u00e9diate, H\u00e9sitation, Rejet}). Le tableau de contingence issu d&rsquo;une enqu\u00eate sur un \u00e9chantillon de 500 participants comportera 2 lignes et 3 colonnes, g\u00e9n\u00e9rant 6 cellules d&rsquo;intersection. L&rsquo;effectif <span class=\"arabpsy-math-inline\" data-latex=\"n_{11}\">$n_{11}$<\/span> indiquera combien d&rsquo;individus sont \u00e0 la fois rationnels et adoptants imm\u00e9diats, tandis que <span class=\"arabpsy-math-inline\" data-latex=\"n_{23}\">$n_{23}$<\/span> d\u00e9nombrera les profils intuitifs ayant rejet\u00e9 l&rsquo;innovation. Cette disposition bidimensionnelle condense l&rsquo;int\u00e9gralit\u00e9 du profil d&rsquo;association empirique sans alt\u00e9rer l&rsquo;individualit\u00e9 de chaque combinaison d&rsquo;attributs.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/marginal1-1.png\" alt=\"\" \/><\/figure>\n<h3>2.2 Normalisation et passage aux fr\u00e9quences relatives conjointes<\/h3>\n<p>Pour s&rsquo;extraire de la d\u00e9pendance \u00e0 la taille d&rsquo;un \u00e9chantillon sp\u00e9cifique et amorcer la transition vers la th\u00e9orie des probabilit\u00e9s, les effectifs bruts doivent subir une op\u00e9ration de normalisation arithm\u00e9tique. Cette \u00e9tape consiste \u00e0 diviser chaque effectif conjoint <span class=\"arabpsy-math-inline\" data-latex=\"n_{ij}\">$n_{ij}$<\/span> par l&rsquo;effectif total d&rsquo;observations <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span>, o\u00f9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"N = \\sum_{i=1}^r \\sum_{j=1}^c n_{ij}\">$$N = \\sum_{i=1}^r \\sum_{j=1}^c n_{ij}$$<\/span><\/p>\n<p>Le ratio obtenu d\u00e9finit la fr\u00e9quence relative conjointe, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"f_{ij} = \\frac{n_{ij}}{N}\">$f_{ij} = \\frac{n_{ij}}{N}$<\/span>. Selon la loi forte des grands nombres, sous des conditions d&rsquo;\u00e9chantillonnage al\u00e9atoire et simple, cette fr\u00e9quence empirique converge presque s\u00fbrement vers la probabilit\u00e9 th\u00e9orique conjointe lorsque la taille de l&rsquo;\u00e9chantillon tend vers l&rsquo;infini : <span class=\"arabpsy-math-inline\" data-latex=\"\\lim_{N to \\infty} f_{ij} = P(X = x_i, Y = y_j)\">$\\lim_{N to \\infty} f_{ij} = P(X = x_i, Y = y_j)$<\/span>. La fr\u00e9quence relative constitue ainsi un <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Estimateur_(statistique)\">estimateur non biais\u00e9<\/a> et convergent de la loi de probabilit\u00e9 conjointe du syst\u00e8me bivari\u00e9.<\/p>\n<p>Une propri\u00e9t\u00e9 fondamentale d\u00e9coule m\u00e9caniquement de cette normalisation : la contrainte de sommabilit\u00e9 unitaire. Puisque la somme de toutes les cellules internes reproduit par construction l&rsquo;effectif total <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span>, la double sommation des fr\u00e9quences relatives conjointes sur l&rsquo;ensemble de la grille matricielle est obligatoirement \u00e9gale \u00e0 l&rsquo;unit\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\sum_{i=1}^r \\sum_{j=1}^c f_{ij} = \\sum_{i=1}^r \\sum_{j=1}^c \\frac{n_{ij}}{N} = \\frac{N}{N} = 1\">$$\\sum_{i=1}^r \\sum_{j=1}^c f_{ij} = \\sum_{i=1}^r \\sum_{j=1}^c \\frac{n_{ij}}{N} = \\frac{N}{N} = 1$$<\/span><\/p>\n<p>Sur le plan \u00e9pist\u00e9mologique, cette \u00e9galit\u00e9 traduit l&rsquo;exhaustivit\u00e9 de la partition de l&rsquo;univers des \u00e9tats possibles : l&rsquo;espace produit des modalit\u00e9s de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et de <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> circonscrit l&rsquo;int\u00e9gralit\u00e9 des r\u00e9alisations concevables. Chaque case du tableau normalis\u00e9 n&rsquo;exprime plus une quantit\u00e9 d&rsquo;individus, mais une mesure de masse probabiliste, fournissant l&rsquo;\u00e9valuation empirique de la propension d&rsquo;un sujet al\u00e9atoirement tir\u00e9 \u00e0 occuper cet \u00e9tat bivari\u00e9 pr\u00e9cis.<\/p>\n<h3>2.3 Notation math\u00e9matique standard pour les variables discr\u00e8tes<\/h3>\n<p>L&rsquo;expression probabiliste rigoureuse des structures observ\u00e9es dans les tableaux crois\u00e9s requiert un appareillage formel unifi\u00e9. Lorsque les variables <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> sont des variables al\u00e9atoires discr\u00e8tes prenant respectivement leurs valeurs dans les ensembles d\u00e9nombrables <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{X} = {x_1, x_2, dots}\">$\\mathcal{X} = {x_1, x_2, dots}$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{Y} = {y_1, y_2, dots}\">$\\mathcal{Y} = {y_1, y_2, dots}$<\/span>, la probabilit\u00e9 que <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> prenne la valeur <span class=\"arabpsy-math-inline\" data-latex=\"x\">$x$<\/span> et que <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> prenne la valeur <span class=\"arabpsy-math-inline\" data-latex=\"y\">$y$<\/span> se formule indistinctement sous les notations \u00e9quivalentes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(X = x \\text{ et } Y = y) \\equiv P({X = x} \\cap {Y = y}) \\equiv P(X = x, Y = y)\">$$P(X = x \\text{ et } Y = y) \\equiv P({X = x} \\cap {Y = y}) \\equiv P(X = x, Y = y)$$<\/span><\/p>\n<p>La virgule s\u00e9parant les deux variables dans l&rsquo;\u00e9criture abr\u00e9g\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"P(X = x, Y = y)\">$P(X = x, Y = y)$<\/span> joue le r\u00f4le d&rsquo;un connecteur logique de conjonction conjonctive. Pour manipuler ais\u00e9ment cet objet, on d\u00e9finit la <strong>fonction de masse conjointe<\/strong> (ou loi de probabilit\u00e9 conjointe), traditionnellement d\u00e9sign\u00e9e par <span class=\"arabpsy-math-inline\" data-latex=\"p_{X, Y}(x, y)\">$p_{X, Y}(x, y)$<\/span> ou plus sobrement par <span class=\"arabpsy-math-inline\" data-latex=\"p(x, y)\">$p(x, y)$<\/span> lorsque l&rsquo;assignation des variables ne pr\u00e9sente aucune ambigu\u00eft\u00e9. Cette fonction est une application formelle dont le domaine de d\u00e9part est le produit cart\u00e9sien des supports discrets et dont l&rsquo;espace d&rsquo;arriv\u00e9e est l&rsquo;intervalle r\u00e9el <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p_{X, Y} : \\mathcal{X} \\times \\mathcal{Y} to [0, 1], \\quad (x, y) \\mapsto P(X = x, Y = y)\">$$p_{X, Y} : \\mathcal{X} \\times \\mathcal{Y} to [0, 1], \\quad (x, y) \\mapsto P(X = x, Y = y)$$<\/span><\/p>\n<p>Lorsque les supports <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{X}\">$\\mathcal{X}$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{Y}\">$\\mathcal{Y}$<\/span> sont finis, comprenant respectivement <span class=\"arabpsy-math-inline\" data-latex=\"r\">$r$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"c\">$c$<\/span> \u00e9l\u00e9ments ordonn\u00e9s, la fonction de masse conjointe admet une repr\u00e9sentation matricielle compacte. On note cette matrice stochastique <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{P} in \\mathbb{R}^{r \\times c}\">$\\mathbf{P} in \\mathbb{R}^{r \\times c}$<\/span>, o\u00f9 chaque coefficient <span class=\"arabpsy-math-inline\" data-latex=\"[\\mathbf{P}]_{ij} = p(x_i, y_j)\">$[\\mathbf{P}]_{ij} = p(x_i, y_j)$<\/span> v\u00e9rifie scrupuleusement les contraintes de non-n\u00e9gativit\u00e9 et de normalisation globale. Cette \u00e9criture alg\u00e9brique s&rsquo;av\u00e8re particuli\u00e8rement puissante pour automatiser les calculs d&rsquo;esp\u00e9rances, de moments crois\u00e9s et de projections marginales via des produits matriciels matriciels vectoriels \u00e9l\u00e9mentaires.<\/p>\n<h2>3. Formalisation math\u00e9matique : Le cas des variables discr\u00e8tes<\/h2>\n<h3>3.1 Axiomatique et propri\u00e9t\u00e9s de la fonction de masse conjointe<\/h3>\n<p>L&rsquo;assise math\u00e9matique d&rsquo;une distribution conjointe discr\u00e8te repose sur l&rsquo;axiomatique de Kolmogorov adapt\u00e9e aux espaces produits. Une fonction arbitraire <span class=\"arabpsy-math-inline\" data-latex=\"p(x, y)\">$p(x, y)$<\/span> d\u00e9finie sur un ensemble produit discret <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{S} = \\mathcal{X} \\times \\mathcal{Y}\">$\\mathcal{S} = \\mathcal{X} \\times \\mathcal{Y}$<\/span> constitue une fonction de masse de probabilit\u00e9 conjointe valide si et seulement si elle satisfait imp\u00e9rativement deux axiomes directeurs :<\/p>\n<ol>\n<li><strong>Axiome de non-n\u00e9gativit\u00e9 :<\/strong> Pour tout couple <span class=\"arabpsy-math-inline\" data-latex=\"(x, y) in \\mathcal{X} \\times \\mathcal{Y}\">$(x, y) in \\mathcal{X} \\times \\mathcal{Y}$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"p(x, y) ge 0\">$p(x, y) ge 0$<\/span>. Aucune probabilit\u00e9 ne peut \u00eatre strictement n\u00e9gative, et les couples appartenant au produit cart\u00e9sien mais impossibles dans les faits se voient assigner une valeur nulle (<span class=\"arabpsy-math-inline\" data-latex=\"p(x, y) = 0\">$p(x, y) = 0$<\/span>).<\/li>\n<li><strong>Axiome de normalisation totale :<\/strong> La double somme \u00e9tendue \u00e0 l&rsquo;ensemble du support conjoint doit converger et valoir rigoureusement l&rsquo;unit\u00e9 :<br \/>\n  <span class=\"arabpsy-math-display\" data-latex=\"\\sum_{x in \\mathcal{X}} \\sum_{y in \\mathcal{Y}} p(x, y) = 1\">$$\\sum_{x in \\mathcal{X}} \\sum_{y in \\mathcal{Y}} p(x, y) = 1$$<\/span><\/li>\n<\/ol>\n<p>\u00c0 partir de ces deux contraintes fondamentales, le calcul de la probabilit\u00e9 d&rsquo;un \u00e9v\u00e9nement arbitraire <span class=\"arabpsy-math-inline\" data-latex=\"A \\subseteq \\mathcal{X} \\times \\mathcal{Y}\">$A \\subseteq \\mathcal{X} \\times \\mathcal{Y}$<\/span> s&rsquo;obtient par l&rsquo;application de la propri\u00e9t\u00e9 d&rsquo;additivit\u00e9 d\u00e9nombrable. Si l&rsquo;\u00e9v\u00e9nement d&rsquo;int\u00e9r\u00eat correspond \u00e0 une r\u00e9gion g\u00e9om\u00e9trique sp\u00e9cifique du plan discret (par exemple, la probabilit\u00e9 que la somme des variables exc\u00e8de un seuil critique <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>), la mesure de cet \u00e9v\u00e9nement se calcule en sommant la fonction de masse conjointe sur l&rsquo;ensemble des couples satisfaisant ce crit\u00e8re :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P((X, Y) in A) = \\sum_{(x, y) in A} p(x, y)\">$$P((X, Y) in A) = \\sum_{(x, y) in A} p(x, y)$$<\/span><\/p>\n<p>La distinction math\u00e9matique entre un support discret fini et un support d\u00e9nombrable infini introduit des subtilit\u00e9s d&rsquo;analyse r\u00e9elle. Lorsque les supports <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{X}\">$\\mathcal{X}$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{Y}\">$\\mathcal{Y}$<\/span> sont infinis d\u00e9nombrables (comme l&rsquo;ensemble des entiers naturels <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{N}\">$\\mathbb{N}$<\/span> dans le cadre de lois de Poisson bivari\u00e9es), la double somme doit faire preuve d&rsquo;une convergence absolue. En vertu du <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Th%C3%A9or%C3%A8me_de_Fubini\">th\u00e9or\u00e8me de Fubini<\/a> pour les s\u00e9ries doubles \u00e0 termes positifs, la non-n\u00e9gativit\u00e9 garantit que l&rsquo;ordre de sommation peut \u00eatre interverti sans alt\u00e9rer la valeur de la somme globale, assurant ainsi la robustesse th\u00e9orique des calculs d&rsquo;esp\u00e9rances et de probabilit\u00e9s marginales.<\/p>\n<h3>3.2 Construction d&rsquo;une fonction de r\u00e9partition conjointe discr\u00e8te<\/h3>\n<p>Parall\u00e8lement \u00e0 la fonction de masse, la description d&rsquo;un vecteur al\u00e9atoire discret s&rsquo;effectue au moyen de la <strong>fonction de r\u00e9partition conjointe<\/strong> (ou fonction de distribution cumulative bivari\u00e9e), not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"F_{X, Y}(x, y)\">$F_{X, Y}(x, y)$<\/span>. Pour un vecteur bidimensionnel <span class=\"arabpsy-math-inline\" data-latex=\"(X, Y)\">$(X, Y)$<\/span>, cette fonction est d\u00e9finie sur l&rsquo;espace continu <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^2\">$\\mathbb{R}^2$<\/span> tout entier par l&rsquo;accumulation des probabilit\u00e9s situ\u00e9es dans le quadrant inf\u00e9rieur gauche ind\u00e9fini :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"F_{X, Y}(x, y) = P(X le x, Y le y) = \\sum_{u le x} \\sum_{v le y} p(u, v)\">$$F_{X, Y}(x, y) = P(X le x, Y le y) = \\sum_{u le x} \\sum_{v le y} p(u, v)$$<\/span><\/p>\n<p>Cette fonction de r\u00e9partition cumulative h\u00e9rite de propri\u00e9t\u00e9s structurelles rigoureuses qui refl\u00e8tent la topologie des probabilit\u00e9s sur le plan euclidien. En premier lieu, <span class=\"arabpsy-math-inline\" data-latex=\"F(x, y)\">$F(x, y)$<\/span> est une fonction monotone non d\u00e9croissante par rapport \u00e0 chacune de ses composantes \u00e9valu\u00e9es s\u00e9par\u00e9ment. Fixer une coordonn\u00e9e tout en augmentant l&rsquo;autre ne peut qu&rsquo;\u00e9largir l&rsquo;ensemble d&rsquo;\u00e9v\u00e9nements cumul\u00e9s, induisant une valeur de probabilit\u00e9 sup\u00e9rieure ou \u00e9gale : si <span class=\"arabpsy-math-inline\" data-latex=\"x_1 le x_2\">$x_1 le x_2$<\/span>, alors <span class=\"arabpsy-math-inline\" data-latex=\"F(x_1, y) le F(x_2, y)\">$F(x_1, y) le F(x_2, y)$<\/span> pour tout <span class=\"arabpsy-math-inline\" data-latex=\"y\">$y$<\/span> fix\u00e9.<\/p>\n<p>En second lieu, le comportement asymptotique aux fronti\u00e8res de l&rsquo;espace infini respecte scrupuleusement les contraintes d&rsquo;absorption probabiliste :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\lim_{x to -\\infty} F(x, y) = 0, \\quad \\lim_{y to -\\infty} F(x, y) = 0\">$$\\lim_{x to -\\infty} F(x, y) = 0, \\quad \\lim_{y to -\\infty} F(x, y) = 0$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\lim_{x to \\infty, y to \\infty} F(x, y) = 1\">$$\\lim_{x to \\infty, y to \\infty} F(x, y) = 1$$<\/span><\/p>\n<p>Dans l&rsquo;univers discret, <span class=\"arabpsy-math-inline\" data-latex=\"F(x, y)\">$F(x, y)$<\/span> pr\u00e9sente une allure g\u00e9om\u00e9trique caract\u00e9ristique en \u00ab escalier bidimensionnel \u00bb, form\u00e9e d&rsquo;une mosa\u00efque de plateaux horizontaux discontinus, continus \u00e0 droite le long de chaque dimension. Un int\u00e9r\u00eat majeur de cette fonction r\u00e9side dans la d\u00e9termination de la probabilit\u00e9 d&rsquo;appartenance \u00e0 un pav\u00e9 rectangulaire semi-ouvert de la forme <span class=\"arabpsy-math-inline\" data-latex=\"]a, b] \\times ]c, d]\">$]a, b] \\times ]c, d]$<\/span>. Par le principe d&rsquo;inclusion-exclusion, cette probabilit\u00e9 ne se r\u00e9sume pas \u00e0 une simple diff\u00e9rence bivari\u00e9e, mais implique quatre \u00e9valuations de la fonction de r\u00e9partition cumulative :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(a &lt; X le b, c &lt; Y le d) = F(b, d) - F(a, d) - F(b, c) + F(a, c)\">$$P(a &lt; X le b, c &lt; Y le d) = F(b, d) &#8211; F(a, d) &#8211; F(b, c) + F(a, c)$$<\/span><\/p>\n<p>Le terme correctif <span class=\"arabpsy-math-inline\" data-latex=\"+ F(a, c)\">$+ F(a, c)$<\/span> compense la double soustraction du quadrant <span class=\"arabpsy-math-inline\" data-latex=\"X le a, Y le c\">$X le a, Y le c$<\/span>, garantissant que la mesure r\u00e9sultante demeure strictement non n\u00e9gative pour tout rectangle consid\u00e9r\u00e9.<\/p>\n<h3>3.3 \u00c9tude de cas num\u00e9rique d\u00e9taill\u00e9 d&rsquo;une loi conjointe discr\u00e8te<\/h3>\n<p>Afin de concr\u00e9tiser ces concepts, formalisons un univers d&rsquo;\u00e9chantillonnage psychom\u00e9trique o\u00f9 deux variables binaires d\u00e9crivent les r\u00e9actions d&rsquo;\u00e9tudiants face \u00e0 une \u00e9preuve d&rsquo;\u00e9valuation sous pression chronom\u00e9tr\u00e9e. D\u00e9finissons la variable <span class=\"arabpsy-math-inline\" data-latex=\"X in {0, 1}\">$X in {0, 1}$<\/span> mesurant l&rsquo;\u00e9tat de panique cognitive (<span class=\"arabpsy-math-inline\" data-latex=\"0 = \\text{serein}, 1 = \\text{paniqu\u00e9}\">$0 = \\text{serein}, 1 = \\text{paniqu\u00e9}$<\/span>) et la variable <span class=\"arabpsy-math-inline\" data-latex=\"Y in {0, 1}\">$Y in {0, 1}$<\/span> indiquant la r\u00e9ussite finale \u00e0 l&rsquo;\u00e9valuation (<span class=\"arabpsy-math-inline\" data-latex=\"0 = \\text{\u00e9chec}, 1 = \\text{succ\u00e8s}\">$0 = \\text{\u00e9chec}, 1 = \\text{succ\u00e8s}$<\/span>). Supposons que l&rsquo;administration du test \u00e0 une vaste population calibr\u00e9e \u00e9tablisse la fonction de masse conjointe <span class=\"arabpsy-math-inline\" data-latex=\"p(x, y)\">$p(x, y)$<\/span> suivante :<\/p>\n<ul>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"p(0, 0) = P(X = 0, Y = 0) = 0{,}10\">$p(0, 0) = P(X = 0, Y = 0) = 0{,}10$<\/span> : l&rsquo;\u00e9tudiant est serein mais \u00e9choue par manque de r\u00e9vision.<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"p(0, 1) = P(X = 0, Y = 1) = 0{,}55\">$p(0, 1) = P(X = 0, Y = 1) = 0{,}55$<\/span> : l&rsquo;\u00e9tudiant est serein et valide avec brio l&rsquo;\u00e9preuve.<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"p(1, 0) = P(X = 1, Y = 0) = 0{,}25\">$p(1, 0) = P(X = 1, Y = 0) = 0{,}25$<\/span> : la panique cognitive entra\u00eene la d\u00e9sorganisation et l&rsquo;\u00e9chec.<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"p(1, 1) = P(X = 1, Y = 1) = 0{,}10\">$p(1, 1) = P(X = 1, Y = 1) = 0{,}10$<\/span> : l&rsquo;\u00e9tudiant surmonte sa panique et r\u00e9ussit malgr\u00e9 tout le test.<\/li>\n<\/ul>\n<p>V\u00e9rifions en premier lieu la conformit\u00e9 aux axiomes fondamentaux : chaque probabilit\u00e9 individuelle est strictement positive (<span class=\"arabpsy-math-inline\" data-latex=\"p(x, y) ge 0\">$p(x, y) ge 0$<\/span>), et la somme totale s&rsquo;\u00e9l\u00e8ve \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"0{,}10 + 0{,}55 + 0{,}25 + 0{,}10 = 1{,}00\">$0{,}10 + 0{,}55 + 0{,}25 + 0{,}10 = 1{,}00$<\/span>. Le cadre kolmogorovien est int\u00e9gralement honor\u00e9. D\u00e9terminons \u00e0 pr\u00e9sent la probabilit\u00e9 de l&rsquo;\u00e9v\u00e9nement composite <span class=\"arabpsy-math-inline\" data-latex=\"E\">$E$<\/span> : \u00ab l&rsquo;\u00e9tudiant a au moins \u00e9prouv\u00e9 de la panique ou a \u00e9chou\u00e9 \u00e0 l&rsquo;\u00e9preuve \u00bb. L&rsquo;ensemble des issues favorables s&rsquo;\u00e9crit <span class=\"arabpsy-math-inline\" data-latex=\"E = {(0, 0), (1, 0), (1, 1)}\">$E = {(0, 0), (1, 0), (1, 1)}$<\/span>.<\/p>\n<p>En mobilisant la propri\u00e9t\u00e9 d&rsquo;additivit\u00e9 pour les \u00e9v\u00e9nements disjoints, nous calculons :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(E) = p(0, 0) + p(1, 0) + p(1, 1) = 0{,}10 + 0{,}25 + 0{,}10 = 0{,}45\">$$P(E) = p(0, 0) + p(1, 0) + p(1, 1) = 0{,}10 + 0{,}25 + 0{,}10 = 0{,}45$$<\/span><\/p>\n<p>Alternativement, le recours \u00e0 l&rsquo;\u00e9v\u00e9nement compl\u00e9mentaire simplifie le calcul : <span class=\"arabpsy-math-inline\" data-latex=\"E^c = {(0, 1)}\">$E^c = {(0, 1)}$<\/span>, d&rsquo;o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"P(E) = 1 - p(0, 1) = 1 - 0{,}55 = 0{,}45\">$P(E) = 1 &#8211; p(0, 1) = 1 &#8211; 0{,}55 = 0{,}45$<\/span>. Ce mod\u00e8le bivari\u00e9 simple met imm\u00e9diatement en exergue l&rsquo;interd\u00e9pendance des deux m\u00e9canismes : la proportion d&rsquo;\u00e9chec est de <span class=\"arabpsy-math-inline\" data-latex=\"0{,}10 \/ (0{,}10 + 0{,}55) \\approx 15{,}4%\">$0{,}10 \/ (0{,}10 + 0{,}55) \\approx 15{,}4%$<\/span> chez les \u00e9tudiants sereins, tandis qu&rsquo;elle culmine \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"0{,}25 \/ (0{,}25 + 0{,}10) \\approx 71{,}4%\">$0{,}25 \/ (0{,}25 + 0{,}10) \\approx 71{,}4%$<\/span> chez les \u00e9tudiants paniqu\u00e9s. La loi conjointe r\u00e9v\u00e8le ainsi la vuln\u00e9rabilit\u00e9 diff\u00e9rentielle induite par l&rsquo;\u00e9tat \u00e9motionnel, une r\u00e9alit\u00e9 que la seule observation isol\u00e9e du taux de r\u00e9ussite global (<span class=\"arabpsy-math-inline\" data-latex=\"65%\">$65%$<\/span>) aurait totalement masqu\u00e9e.<\/p>\n<h2>4. La distribution conjointe pour variables continues : Densit\u00e9 conjointe<\/h2>\n<h3>4.1 Concept et d\u00e9finition de la fonction de densit\u00e9 conjointe<\/h3>\n<p>Lorsque nous franchissons la fronti\u00e8re des variables quantitatives continues \u2014 telles que le temps d&rsquo;inf\u00e9rence neuronale, la concentration s\u00e9rique d&rsquo;un neurotransmetteur ou la temp\u00e9rature corporelle \u2014, l&rsquo;appareil conceptuel de la fonction de masse s&rsquo;av\u00e8re caduc. Dans un continuum r\u00e9el, la cardinalit\u00e9 du support d&rsquo;une variable devient infinie non d\u00e9nombrable. En cons\u00e9quence, la probabilit\u00e9 ponctuelle exacte affect\u00e9e \u00e0 n&rsquo;importe quel point pr\u00e9cis <span class=\"arabpsy-math-inline\" data-latex=\"(x, y) in \\mathbb{R}^2\">$(x, y) in \\mathbb{R}^2$<\/span> est rigoureusement nulle :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(X = x, Y = y) = 0, \\quad \\forall (x, y) in \\mathbb{R}^2\">$$P(X = x, Y = y) = 0, \\quad \\forall (x, y) in \\mathbb{R}^2$$<\/span><\/p>\n<p>Ce paradoxe apparent ne signifie nullement que l&rsquo;\u00e9v\u00e9nement est impossible, mais traduit la granularit\u00e9 infinie de l&rsquo;espace continu. Pour quantifier l&rsquo;incertitude dans un tel environnement, la th\u00e9orie de la mesure introduit la <strong>fonction de densit\u00e9 conjointe<\/strong> (ou fonction de densit\u00e9 de probabilit\u00e9 bivari\u00e9e), usuellement not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"f_{X, Y}(x, y)\">$f_{X, Y}(x, y)$<\/span> ou simplement <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y)\">$f(x, y)$<\/span>. Cette fonction se d\u00e9finit rigoureusement comme la d\u00e9riv\u00e9e partielle seconde crois\u00e9e de la fonction de r\u00e9partition conjointe bivari\u00e9e, l\u00e0 o\u00f9 celle-ci est diff\u00e9rentiable :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f(x, y) = \\frac{\\partial^2 F(x, y)}{\\partial x \\partial y}\">$$f(x, y) = \\frac{\\partial^2 F(x, y)}{\\partial x \\partial y}$$<\/span><\/p>\n<p>Le sens physique de <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y)\">$f(x, y)$<\/span> ne correspond donc pas \u00e0 une probabilit\u00e9 brute, mais \u00e0 une <em>densit\u00e9 surfacique de probabilit\u00e9<\/em> sur le plan euclidien. La quantit\u00e9 infinit\u00e9simale <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y) , dx , dy\">$f(x, y) , dx , dy$<\/span> exprime la probabilit\u00e9 que le vecteur stochastique <span class=\"arabpsy-math-inline\" data-latex=\"(X, Y)\">$(X, Y)$<\/span> vienne se loger au sein du rectangle diff\u00e9rentiel \u00e9l\u00e9mentaire <span class=\"arabpsy-math-inline\" data-latex=\"[x, x + dx] \\times [y, y + dy]\">$[x, x + dx] \\times [y, y + dy]$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(x le X le x + dx, y le Y le y + dy) \\approx f(x, y) , dx , dy\">$$P(x le X le x + dx, y le Y le y + dy) \\approx f(x, y) , dx , dy$$<\/span><\/p>\n<p>Pour \u00eatre qualifi\u00e9e de densit\u00e9 conjointe admissible au sens des th\u00e9ories de Riemann et de <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Int%C3%A9grale_de_Lebesgue\">Lebesgue<\/a>, la fonction <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y)\">$f(x, y)$<\/span> doit \u00eatre localement int\u00e9grable, mesurable au sens de Borel, et se plier \u00e0 la condition de non-n\u00e9gativit\u00e9 absolue : <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y) ge 0\">$f(x, y) ge 0$<\/span> pour tout <span class=\"arabpsy-math-inline\" data-latex=\"(x, y) in \\mathbb{R}^2\">$(x, y) in \\mathbb{R}^2$<\/span>. Contrairement aux probabilit\u00e9s ponctuelles discr\u00e8tes, la valeur locale <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y)\">$f(x, y)$<\/span> peut parfaitement exc\u00e9der la valeur 1, \u00e0 la condition expresse que son int\u00e9grale surfacique globale demeure strictement contrainte.<\/p>\n<h3>4.2 Calcul int\u00e9gral des probabilit\u00e9s sur des domaines continus<\/h3>\n<p>L&rsquo;extraction d&rsquo;une mesure de probabilit\u00e9 concr\u00e8te \u00e0 partir d&rsquo;une fonction de densit\u00e9 bivari\u00e9e requiert le passage \u00e0 l&rsquo;int\u00e9gration sur des domaines bidimensionnels. La probabilit\u00e9 que le couple al\u00e9atoire <span class=\"arabpsy-math-inline\" data-latex=\"(X, Y)\">$(X, Y)$<\/span> appartienne \u00e0 une r\u00e9gion g\u00e9om\u00e9trique arbitraire <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{D} \\subset \\mathbb{R}^2\">$\\mathcal{D} \\subset \\mathbb{R}^2$<\/span> est donn\u00e9e par l&rsquo;int\u00e9grale double de la densit\u00e9 conjointe sur ce domaine :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P((X, Y) in \\mathcal{D}) = \\iint_{\\mathcal{D}} f(x, y) , dx , dy\">$$P((X, Y) in \\mathcal{D}) = \\iint_{\\mathcal{D}} f(x, y) , dx , dy$$<\/span><\/p>\n<p>La contrainte de normalisation fondamentale impose que l&rsquo;int\u00e9grale double \u00e9tendue au plan r\u00e9el infini tout entier soit \u00e9gale \u00e0 l&rsquo;unit\u00e9 sans \u00e9quivoque :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\int_{-\\infty}^{\\infty} \\int_{-\\infty}^{\\infty} f(x, y) , dx , dy = 1\">$$\\int_{-\\infty}^{\\infty} \\int_{-\\infty}^{\\infty} f(x, y) , dx , dy = 1$$<\/span><\/p>\n<p>Lorsque le domaine d&rsquo;int\u00e9gration <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{D}\">$\\mathcal{D}$<\/span> n&rsquo;est pas un rectangle euclidien standard \u00e0 bornes ind\u00e9pendantes, mais pr\u00e9sente une g\u00e9om\u00e9trie complexe d\u00e9limit\u00e9e par des fronti\u00e8res fonctionnelles (par exemple, un triangle ou un cercle), la d\u00e9finition minutieuse des bornes d&rsquo;int\u00e9gration devient une op\u00e9ration critique. Si le domaine s&rsquo;\u00e9crit <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{D} = {(x, y) in \\mathbb{R}^2 mid a le x le b, , g_1(x) le y le g_2(x)}\">$\\mathcal{D} = {(x, y) in \\mathbb{R}^2 mid a le x le b, , g_1(x) le y le g_2(x)}$<\/span>, le calcul s&rsquo;op\u00e8re par int\u00e9grales it\u00e9r\u00e9es en respectant l&rsquo;ordre d&rsquo;imbrication des variables :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P((X, Y) in \\mathcal{D}) = \\int_{a}^{b} \\left( \\int_{g_1(x)}^{g_2(x)} f(x, y) , dy \\right) dx\">$$P((X, Y) in \\mathcal{D}) = \\int_{a}^{b} \\left( \\int_{g_1(x)}^{g_2(x)} f(x, y) , dy \\right) dx$$<\/span><\/p>\n<p>Dans de nombreuses configurations analytiques \u2014 notamment lors de l&rsquo;\u00e9valuation de densit\u00e9s normales bivari\u00e9es ou de lois radiales \u2014, la g\u00e9om\u00e9trie cart\u00e9sienne induit des int\u00e9grales insolubles sous leur forme primitive. Le recours \u00e0 un changement de variables bivari\u00e9 s&rsquo;av\u00e8re alors indispensable. Si l&rsquo;on effectue la transformation <span class=\"arabpsy-math-inline\" data-latex=\"u = \\phi_1(x, y)\">$u = \\phi_1(x, y)$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"v = \\phi_2(x, y)\">$v = \\phi_2(x, y)$<\/span>, la distribution conjointe transform\u00e9e fait intervenir la valeur absolue du d\u00e9terminant jacobien de la transformation inverse :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_{U, V}(u, v) = f_{X, Y}(x(u, v), y(u, v)) \\cdot left| det \\begin{\\pmatrix} \\frac{\\partial x}{\\partial u} &amp;a\\mp; \\frac{\\partial x}{\\partial v}  \\frac{\\partial y}{\\partial u} &amp;a\\mp; \\frac{\\partial y}{\\partial v} \\end{\\pmatrix} right|\">$$f_{U, V}(u, v) = f_{X, Y}(x(u, v), y(u, v)) \\cdot left| det \\begin{\\pmatrix} \\frac{\\partial x}{\\partial u} &#038;a\\mp; \\frac{\\partial x}{\\partial v}  \\frac{\\partial y}{\\partial u} &#038;a\\mp; \\frac{\\partial y}{\\partial v} \\end{\\pmatrix} right|$$<\/span><\/p>\n<p>Ce multiplicateur diff\u00e9rentiel compense la distorsion g\u00e9om\u00e9trique locale de l&rsquo;\u00e9l\u00e9ment de surface induite par la transformation, garantissant la stricte conservation de la masse probabiliste.<\/p>\n<h3>4.3 Fonction de r\u00e9partition cumulative bivari\u00e9e continue<\/h3>\n<p>\u00c0 l&rsquo;instar du contexte discret, la <strong>fonction de r\u00e9partition cumulative conjointe<\/strong> dans l&rsquo;espace continu cumule la densit\u00e9 de probabilit\u00e9 depuis l&rsquo;infini n\u00e9gatif jusqu&rsquo;au point d&rsquo;\u00e9valuation <span class=\"arabpsy-math-inline\" data-latex=\"(x, y)\">$(x, y)$<\/span>. Elle s&rsquo;exprime analytiquement par l&rsquo;int\u00e9grale param\u00e9trique bivari\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"F(x, y) = P(X le x, Y le y) = \\int_{-\\infty}^{x} \\int_{-\\infty}^{y} f(u, v) , dv , du\">$$F(x, y) = P(X le x, Y le y) = \\int_{-\\infty}^{x} \\int_{-\\infty}^{y} f(u, v) , dv , du$$<\/span><\/p>\n<p>Sur le plan g\u00e9om\u00e9trique, <span class=\"arabpsy-math-inline\" data-latex=\"F(x, y)\">$F(x, y)$<\/span> repr\u00e9sente le volume sous la surface de densit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"z = f(u, v)\">$z = f(u, v)$<\/span> \u00e9rig\u00e9e au-dessus de la r\u00e9gion quadrantale d\u00e9limit\u00e9e par <span class=\"arabpsy-math-inline\" data-latex=\"{u le x, v le y}\">${u le x, v le y}$<\/span>. Par le th\u00e9or\u00e8me fondamental de l&rsquo;analyse \u00e9tendu au calcul diff\u00e9rentiel multidimensionnel, d\u00e9river contin\u00fbment cette fonction cumulative permet de reconstituer instantan\u00e9ment la densit\u00e9 conjointe locale :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\frac{\\partial^2 F(x, y)}{\\partial x \\partial y} = f(x, y)\">$$\\frac{\\partial^2 F(x, y)}{\\partial x \\partial y} = f(x, y)$$<\/span><\/p>\n<p>La probabilit\u00e9 qu&rsquo;un vecteur stochastique r\u00e9side dans une bo\u00eete euclidienne <span class=\"arabpsy-math-inline\" data-latex=\"[x_1, x_2] \\times [y_1, y_2]\">$[x_1, x_2] \\times [y_1, y_2]$<\/span> s&rsquo;obtient avec une parfaite \u00e9l\u00e9gance alg\u00e9brique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(x_1 le X le x_2, y_1 le Y le y_2) = F(x_2, y_2) - F(x_1, y_2) - F(x_2, y_1) + F(x_1, y_1)\">$$P(x_1 le X le x_2, y_1 le Y le y_2) = F(x_2, y_2) &#8211; F(x_1, y_2) &#8211; F(x_2, y_1) + F(x_1, y_1)$$<\/span><\/p>\n<p>Cette relation met en lumi\u00e8re l&rsquo;\u00e9quivalence op\u00e9rationnelle entre l&rsquo;int\u00e9gration de la densit\u00e9 et la diff\u00e9rentiation de la fonction de r\u00e9partition. Visuellement, une fonction de densit\u00e9 continue se projette dans l&rsquo;espace sous la forme d&rsquo;un paysage tridimensionnel : une surface en d\u00f4me dont l&rsquo;\u00e9l\u00e9vation indique les zones de forte concentration stochastique. L&rsquo;intersection de cette surface par des plans d&rsquo;altitude constante dessine des <em>courbes d&rsquo;isodensit\u00e9<\/em> sur le plan orthogonal <span class=\"arabpsy-math-inline\" data-latex=\"(X, Y)\">$(X, Y)$<\/span>, fournissant une repr\u00e9sentation topologique compacte de la morphologie conjointe et des axes pr\u00e9f\u00e9rentiels de variation stochastique.<\/p>\n<h2>5. Extraction des distributions marginales \u00e0 partir de la loi conjointe<\/h2>\n<h3>5.1 Le processus de marginalisation discr\u00e8te par sommation<\/h3>\n<p>L&rsquo;un des apports majeurs de la th\u00e9orie de la probabilit\u00e9 r\u00e9side dans sa coh\u00e9rence projective : poss\u00e9der la distribution conjointe d&rsquo;un syst\u00e8me multivari\u00e9 permet de retrouver automatiquement les lois r\u00e9gissant chaque variable prise isol\u00e9ment. Cette op\u00e9ration de r\u00e9duction dimensionnelle porte le nom de <strong>marginalisation<\/strong>. Son \u00e9tymologie renvoie directement \u00e0 la pratique des statisticiens qui, dans les tableaux de contingence, notaient la somme des fr\u00e9quences de chaque ligne et de chaque colonne dans les marges blanches du registre papier.<\/p>\n<p>Sur le plan math\u00e9matique, pour un couple de variables al\u00e9atoires discr\u00e8tes <span class=\"arabpsy-math-inline\" data-latex=\"(X, Y)\">$(X, Y)$<\/span> dot\u00e9 d&rsquo;une fonction de masse conjointe <span class=\"arabpsy-math-inline\" data-latex=\"p_{X, Y}(x, y)\">$p_{X, Y}(x, y)$<\/span>, la <strong>loi marginale<\/strong> de la variable <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"p_X(x)\">$p_X(x)$<\/span>, s&rsquo;obtient en sommant les probabilit\u00e9s conjointes sur la totalit\u00e9 des valeurs concevables que peut prendre la variable concomitante <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p_X(x) = P(X = x) = \\sum_{y in \\mathcal{Y}} p_{X, Y}(x, y)\">$$p_X(x) = P(X = x) = \\sum_{y in \\mathcal{Y}} p_{X, Y}(x, y)$$<\/span><\/p>\n<p>De mani\u00e8re sym\u00e9trique, la distribution marginale de la variable <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> s&rsquo;obtient en \u00ab int\u00e9grant \u00bb ou sommant sur l&rsquo;ensemble du support de la variable <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p_Y(y) = P(Y = y) = \\sum_{x in \\mathcal{X}} p_{X, Y}(x, y)\">$$p_Y(y) = P(Y = y) = \\sum_{x in \\mathcal{X}} p_{X, Y}(x, y)$$<\/span><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/marginal3-1.png\" alt=\"Marginal distribution example with two-way table\" \/><figcaption>Marginal distribution example with two-way table<\/figcaption><\/figure>\n<p>L&rsquo;op\u00e9ration de marginalisation dissout litt\u00e9ralement l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 introduite par la variable neutralis\u00e9e. En sommant selon les colonnes pour une ligne <span class=\"arabpsy-math-inline\" data-latex=\"x\">$x$<\/span> donn\u00e9e, on regroupe tous les \u00e9tats o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"X = x\">$X = x$<\/span>, quelle que soit l&rsquo;issue constat\u00e9e pour <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>. Cette agr\u00e9gation pr\u00e9serve fid\u00e8lement la distribution univari\u00e9e de la variable focale tout en \u00e9liminant d\u00e9lib\u00e9r\u00e9ment l&rsquo;information relative \u00e0 sa structure de d\u00e9pendance crois\u00e9e avec le second al\u00e9a.<\/p>\n<h3>5.2 Marginalisation continue par int\u00e9gration partielle<\/h3>\n<p>Dans l&rsquo;univers des variables continues gouvern\u00e9es par une fonction de densit\u00e9 bivari\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"f_{X, Y}(x, y)\">$f_{X, Y}(x, y)$<\/span>, la sommation discr\u00e8te c\u00e8de la place \u00e0 l&rsquo;int\u00e9gration partielle d\u00e9finie sur toute l&rsquo;\u00e9tendue de l&rsquo;espace orthogonal. La <strong>fonction de densit\u00e9 marginale<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"f_X(x)\">$f_X(x)$<\/span> de la variable al\u00e9atoire <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> se calcule en int\u00e9grant la densit\u00e9 conjointe par rapport \u00e0 la variable <span class=\"arabpsy-math-inline\" data-latex=\"y\">$y$<\/span>, le long de son support complet allant de <span class=\"arabpsy-math-inline\" data-latex=\"-\\infty\">$-\\infty$<\/span> \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"+\\infty\">$+\\infty$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_X(x) = \\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dy\">$$f_X(x) = \\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dy$$<\/span><\/p>\n<p>R\u00e9ciproquement, la densit\u00e9 marginale de la variable continue <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> r\u00e9sulte de la projection int\u00e9grale de la surface conjointe le long de l&rsquo;axe des abscisses r\u00e9gies par <span class=\"arabpsy-math-inline\" data-latex=\"x\">$x$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_Y(y) = \\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dx\">$$f_Y(y) = \\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dx$$<\/span><\/p>\n<p>G\u00e9om\u00e9triquement, cette marginalisation s&rsquo;apparente \u00e0 un effondrement dimensionnel : on calcule l&rsquo;aire de la section transversale d\u00e9coup\u00e9e dans le volume de densit\u00e9 conjointe \u00e0 une coordonn\u00e9e fix\u00e9e. En effectuant cette op\u00e9ration pour chaque coordonn\u00e9e, on projette la totalit\u00e9 de la masse tridimensionnelle sur une ligne plane univari\u00e9e. Il est imm\u00e9diat de v\u00e9rifier que les fonctions ainsi obtenues respectent scrupuleusement les exigences r\u00e9gissant les densit\u00e9s univari\u00e9es. L&rsquo;int\u00e9gration de la marginale <span class=\"arabpsy-math-inline\" data-latex=\"f_X(x)\">$f_X(x)$<\/span> sur son domaine r\u00e9tablit l&rsquo;int\u00e9grale double globale du syst\u00e8me :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\int_{-\\infty}^{\\infty} f_X(x) , dx = \\int_{-\\infty}^{\\infty} \\left( \\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dy \\right) dx = 1\">$$\\int_{-\\infty}^{\\infty} f_X(x) , dx = \\int_{-\\infty}^{\\infty} \\left( \\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dy \\right) dx = 1$$<\/span><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/marginal4-1.png\" alt=\"Marginal distribution example\" \/><figcaption>Marginal distribution example<\/figcaption><\/figure>\n<p>Ce proc\u00e9d\u00e9 de d\u00e9rivation analytique permet de retrouver les lois univari\u00e9es classiques (comme une loi b\u00eata, gamma ou gaussienne) \u00e0 partir de constructions conjointes complexes d\u00e9crivant des interactions physiques ou psychophysiques intriqu\u00e9es.<\/p>\n<h3>5.3 Perte d&rsquo;information lors du passage de la loi conjointe aux marginales<\/h3>\n<p>Il est imp\u00e9ratif de souligner un principe m\u00e9thodologique capital : <em>l&rsquo;op\u00e9ration de marginalisation est irr\u00e9versible en l&rsquo;absence d&rsquo;hypoth\u00e8ses structurelles fortes<\/em>. Conna\u00eetre la distribution conjointe <span class=\"arabpsy-math-inline\" data-latex=\"p(x, y)\">$p(x, y)$<\/span> conf\u00e8re une omniscience probabiliste sur le syst\u00e8me bivari\u00e9, permettant d&rsquo;extraire sans ambigu\u00eft\u00e9 les deux distributions marginales <span class=\"arabpsy-math-inline\" data-latex=\"p_X(x)\">$p_X(x)$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"p_Y(y)\">$p_Y(y)$<\/span>. En revanche, la r\u00e9ciproque est rigoureusement fausse : poss\u00e9der les marginales <span class=\"arabpsy-math-inline\" data-latex=\"p_X(x)\">$p_X(x)$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"p_Y(y)\">$p_Y(y)$<\/span> ne permet en aucun cas de reconstituer la loi conjointe dont elles sont issues.<\/p>\n<p>Cette asym\u00e9trie informationnelle d\u00e9coule du fait que la marginalisation condense et dissipe la structure de d\u00e9pendance mutuelle. Pour s&rsquo;en convaincre math\u00e9matiquement, consid\u00e9rons deux variables de Bernoulli prenant les valeurs dans <span class=\"arabpsy-math-inline\" data-latex=\"{0, 1}\">${0, 1}$<\/span>. D\u00e9finissons une premi\u00e8re distribution conjointe o\u00f9 les deux variables sont en d\u00e9pendance d\u00e9terministe parfaite : <span class=\"arabpsy-math-inline\" data-latex=\"P(0, 0) = 0{,}5\">$P(0, 0) = 0{,}5$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"P(1, 1) = 0{,}5\">$P(1, 1) = 0{,}5$<\/span>, et <span class=\"arabpsy-math-inline\" data-latex=\"P(0, 1) = P(1, 0) = 0\">$P(0, 1) = P(1, 0) = 0$<\/span>. Les lois marginales associ\u00e9es attribuent trivialement la probabilit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"0{,}5\">$0{,}5$<\/span> \u00e0 chaque \u00e9tat pour les deux variables (<span class=\"arabpsy-math-inline\" data-latex=\"P(X = 0) = P(X = 1) = 0{,}5\">$P(X = 0) = P(X = 1) = 0{,}5$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"P(Y = 0) = P(Y = 1) = 0{,}5\">$P(Y = 0) = P(Y = 1) = 0{,}5$<\/span>).<\/p>\n<p>Consid\u00e9rons \u00e0 pr\u00e9sent une seconde distribution conjointe r\u00e9gie par l&rsquo;ind\u00e9pendance stochastique totale : <span class=\"arabpsy-math-inline\" data-latex=\"P(0, 0) = 0{,}25\">$P(0, 0) = 0{,}25$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"P(0, 1) = 0{,}25\">$P(0, 1) = 0{,}25$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"P(1, 0) = 0{,}25\">$P(1, 0) = 0{,}25$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"P(1, 1) = 0{,}25\">$P(1, 1) = 0{,}25$<\/span>. Un calcul marginal imm\u00e9diat r\u00e9v\u00e8le que les distributions univari\u00e9es de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et de <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> sont rigoureusement identiques \u00e0 celles du premier sc\u00e9nario, avec <span class=\"arabpsy-math-inline\" data-latex=\"50%\">$50%$<\/span> de chances pour chaque modalit\u00e9. Pourtant, la nature de l&rsquo;association empirique est radicalement divergente : dans un cas, la connaissance de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> pr\u00e9dit avec une certitude absolue l&rsquo;\u00e9tat de <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> ; dans l&rsquo;autre, elle n&rsquo;apporte aucune r\u00e9duction d&rsquo;incertitude.<\/p>\n<p>Les r\u00e9percussions de ce constat en psychom\u00e9trie et en sciences du comportement sont profondes. Deux cohortes cliniques peuvent pr\u00e9senter des distributions parfaitement superposables de sympt\u00f4mes anxieux et de troubles d\u00e9pressifs lorsqu&rsquo;on les \u00e9value \u00e9chelle par \u00e9chelle, tout en abritant des distributions conjointes totalement dissemblables. Dans le premier groupe, l&rsquo;anxi\u00e9t\u00e9 et la d\u00e9pression peuvent coexister syst\u00e9matiquement chez les m\u00eames patients, tandis que dans le second, les pathologies touchent des sous-groupes mutuellement disjoints. Restreindre son analyse aux seules lois marginales induit une c\u00e9cit\u00e9 diagnostique majeure, occultant la dynamique syndromique r\u00e9elle sous-tendue par la distribution conjointe.<\/p>\n<h2>6. Conditionnement probabiliste et liens bay\u00e9siens<\/h2>\n<h3>6.1 D\u00e9finition math\u00e9matique de la distribution conditionnelle<\/h3>\n<p>L&rsquo;observation du monde r\u00e9el s&rsquo;op\u00e8re presque toujours sous contrainte d&rsquo;information partielle : nous apprenons qu&rsquo;un \u00e9v\u00e9nement s&rsquo;est produit, et nous cherchons \u00e0 recalculer la probabilit\u00e9 d&rsquo;une autre variable \u00e0 la lumi\u00e8re de cette acquisition empirique. Ce processus d&rsquo;actualisation est formellement d\u00e9crit par la <strong>distribution conditionnelle<\/strong>. Sur le plan structural, la distribution conditionnelle repr\u00e9sente le quotient normalis\u00e9 de la distribution conjointe par la distribution marginale de la variable conditionnante.<\/p>\n<p>Pour des variables al\u00e9atoires discr\u00e8tes, la probabilit\u00e9 conditionnelle que la variable <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> prenne la valeur <span class=\"arabpsy-math-inline\" data-latex=\"x\">$x$<\/span> sachant que la variable <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> a \u00e9t\u00e9 observ\u00e9e \u00e0 la valeur <span class=\"arabpsy-math-inline\" data-latex=\"y\">$y$<\/span> (sous la condition pr\u00e9alable imp\u00e9rative que <span class=\"arabpsy-math-inline\" data-latex=\"p_Y(y) &gt; 0\">$p_Y(y) &gt; 0$<\/span>) est formul\u00e9e par :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p_{X mid Y}(x mid y) = P(X = x mid Y = y) = \\frac{p_{X, Y}(x, y)}{p_Y(y)} = \\frac{p_{X, Y}(x, y)}{\\sum_{u in \\mathcal{X}} p_{X, Y}(u, y)}\">$$p_{X mid Y}(x mid y) = P(X = x mid Y = y) = \\frac{p_{X, Y}(x, y)}{p_Y(y)} = \\frac{p_{X, Y}(x, y)}{\\sum_{u in \\mathcal{X}} p_{X, Y}(u, y)}$$<\/span><\/p>\n<p>Dans l&rsquo;univers des variables continues, malgr\u00e9 le fait que la probabilit\u00e9 de l&rsquo;\u00e9v\u00e9nement conditionnant soit ponctuellement nulle (<span class=\"arabpsy-math-inline\" data-latex=\"P(Y = y) = 0\">$P(Y = y) = 0$<\/span>), la <strong>fonction de densit\u00e9 conditionnelle<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"f_{X mid Y}(x mid y)\">$f_{X mid Y}(x mid y)$<\/span> est rigoureusement d\u00e9finie comme la limite du rapport de mesures sur des voisinages \u00e9vanescents, aboutissant \u00e0 l&rsquo;expression homotrope :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_{X mid Y}(x mid y) = \\frac{f_{X, Y}(x, y)}{f_Y(y)}, \\quad \\text{pour tout } y \\text{ tel que } f_Y(y) &gt; 0\">$$f_{X mid Y}(x mid y) = \\frac{f_{X, Y}(x, y)}{f_Y(y)}, \\quad \\text{pour tout } y \\text{ tel que } f_Y(y) &gt; 0$$<\/span><\/p>\n<p>Sur le plan morphologique, la distribution conditionnelle correspond \u00e0 l&rsquo;isolation d&rsquo;une \u00ab tranche \u00bb unidimensionnelle d\u00e9coup\u00e9e dans la surface conjointe le long de l&rsquo;hyperplan orthogonal d\u00e9fini par la condition <span class=\"arabpsy-math-inline\" data-latex=\"Y = y\">$Y = y$<\/span>. Le d\u00e9nominateur marginal <span class=\"arabpsy-math-inline\" data-latex=\"f_Y(y)\">$f_Y(y)$<\/span> fait office de facteur d&rsquo;\u00e9chelle r\u00e9gulateur : il dilate ou contracte verticalement cette coupe g\u00e9om\u00e9trique pour contraindre son aire totale \u00e0 valoir rigoureusement 1, respectant ainsi l&rsquo;axiomatique \u00e9l\u00e9mentaire d&rsquo;une mesure de probabilit\u00e9 univari\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\int_{-\\infty}^{\\infty} f_{X mid Y}(x mid y) , dx = \\frac{\\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dx}{f_Y(y)} = \\frac{f_Y(y)}{f_Y(y)} = 1\">$$\\int_{-\\infty}^{\\infty} f_{X mid Y}(x mid y) , dx = \\frac{\\int_{-\\infty}^{\\infty} f_{X, Y}(x, y) , dx}{f_Y(y)} = \\frac{f_Y(y)}{f_Y(y)} = 1$$<\/span><\/p>\n<h3>6.2 La r\u00e8gle du produit et la cha\u00eene des probabilit\u00e9s<\/h3>\n<p>La relation de d\u00e9finition du conditionnement peut \u00eatre r\u00e9arrang\u00e9e pour d\u00e9boucher sur un principe op\u00e9ratoire cardinal : la <strong>r\u00e8gle du produit<\/strong>. Celle-ci stipule que la probabilit\u00e9 conjointe de deux \u00e9v\u00e9nements s&rsquo;exprime comme le produit de la probabilit\u00e9 marginale de l&rsquo;un par la probabilit\u00e9 conditionnelle de l&rsquo;autre sachant le premier :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p_{X, Y}(x, y) = p_Y(y) \\cdot p_{X mid Y}(x mid y) = p_X(x) \\cdot p_{Y mid X}(y mid x)\">$$p_{X, Y}(x, y) = p_Y(y) \\cdot p_{X mid Y}(x mid y) = p_X(x) \\cdot p_{Y mid X}(y mid x)$$<\/span><\/p>\n<p>Cette sym\u00e9trie op\u00e9ratoire met en relief la plasticit\u00e9 de la mod\u00e9lisation probabiliste : une m\u00eame distribution conjointe admet deux factorisations directionnelles distinctes, permettant de choisir l&rsquo;ordre de d\u00e9composition le plus naturel ou le plus ais\u00e9ment estimable selon le protocole exp\u00e9rimental en vigueur.<\/p>\n<p>Cette logique se g\u00e9n\u00e9ralise de mani\u00e8re imm\u00e9diate \u00e0 un syst\u00e8me form\u00e9 d&rsquo;un nombre arbitraire <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> de variables al\u00e9atoires gr\u00e2ce \u00e0 la <strong>r\u00e8gle de factorisation en cha\u00eene<\/strong> (ou formule des probabilit\u00e9s compos\u00e9es) :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p(x_1, x_2, dots, x_n) = p(x_1) \\cdot p(x_2 mid x_1) \\cdot p(x_3 mid x_1, x_2) \\cdots p(x_n mid x_1, dots, x_{n-1}) = \\prod_{k=1}^n p\\left(x_k mid x_1, dots, x_{k-1}\\right)\">$$p(x_1, x_2, dots, x_n) = p(x_1) \\cdot p(x_2 mid x_1) \\cdot p(x_3 mid x_1, x_2) \\cdots p(x_n mid x_1, dots, x_{n-1}) = \\prod_{k=1}^n p\\left(x_k mid x_1, dots, x_{k-1}\\right)$$<\/span><\/p>\n<p>Ce d\u00e9veloppement constitue le fondement architectural de la mod\u00e9lisation graphique probabiliste et des <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/R%C3%A9seau_bay%C3%A9sien\">r\u00e9seaux bay\u00e9siens<\/a>. Dans ces repr\u00e9sentations en graphes acycliques orient\u00e9s (DAG), l&rsquo;immense complexit\u00e9 combinatoire d&rsquo;une loi conjointe de dimension \u00e9lev\u00e9e se d\u00e9compose en un produit de conditionnements locaux simplifi\u00e9s, tirant parti des hypoth\u00e8ses d&rsquo;ind\u00e9pendance conditionnelle inh\u00e9rentes \u00e0 la topologie du r\u00e9seau.<\/p>\n<h3>6.3 Actualisation bay\u00e9sienne et inversion des conditionnements<\/h3>\n<p>L&rsquo;articulation entre distribution conjointe, marginales et lois conditionnelles culmine avec le <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Th%C3%A9or%C3%A8me_de_Bayes\">th\u00e9or\u00e8me de Bayes<\/a>. Ce th\u00e9or\u00e8me orchestre l&rsquo;inversion rigoureuse des conditionnements stochastiques : il formalise la fa\u00e7on dont la connaissance a posteriori sur une variable latente, un \u00e9tat inobserv\u00e9 ou un param\u00e8tre th\u00e9orique <span class=\"arabpsy-math-inline\" data-latex=\"Theta\">$Theta$<\/span> s&rsquo;ajuste apr\u00e8s la prise en compte d&rsquo;un vecteur d&rsquo;observations conjointes <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_{\\Theta mid X}(\\theta mid x) = \\frac{f_{X, \\Theta}(x, \\theta)}{f_X(x)} = \\frac{f_\\Theta(\\theta) \\cdot f_{X mid \\Theta}(x mid \\theta)}{\\int_{\\Omega_\\Theta} f_\\Theta(\\vartheta) \\cdot f_{X mid \\Theta}(x mid \\vartheta) , d\\vartheta}\">$$f_{\\Theta mid X}(\\theta mid x) = \\frac{f_{X, \\Theta}(x, \\theta)}{f_X(x)} = \\frac{f_\\Theta(\\theta) \\cdot f_{X mid \\Theta}(x mid \\theta)}{\\int_{\\Omega_\\Theta} f_\\Theta(\\vartheta) \\cdot f_{X mid \\Theta}(x mid \\vartheta) , d\\vartheta}$$<\/span><\/p>\n<p>Dans ce formalisme continu, <span class=\"arabpsy-math-inline\" data-latex=\"f_\\Theta(\\theta)\">$f_\\Theta(\\theta)$<\/span> repr\u00e9sente la densit\u00e9 de probabilit\u00e9 a priori (prior), traduisant l&rsquo;\u00e9tat de croyance ou de plausibilit\u00e9 \u00e9pist\u00e9mique avant toute mesure. Le terme <span class=\"arabpsy-math-inline\" data-latex=\"f_{X mid \\Theta}(x mid \\theta)\">$f_{X mid \\Theta}(x mid \\theta)$<\/span>, envisag\u00e9 en fonction de <span class=\"arabpsy-math-inline\" data-latex=\"\\theta\">$\\theta$<\/span> pour une r\u00e9alisation empirique <span class=\"arabpsy-math-inline\" data-latex=\"x\">$x$<\/span> observ\u00e9e, correspond \u00e0 la fonction de <em>vraisemblance<\/em> conjointe des donn\u00e9es. Enfin, le d\u00e9nominateur <span class=\"arabpsy-math-inline\" data-latex=\"f_X(x)\">$f_X(x)$<\/span> constitue la densit\u00e9 pr\u00e9dictive marginale des observations \u2014 ou constante de normalisation \u2014, obtenue par int\u00e9gration de la densit\u00e9 conjointe globale <span class=\"arabpsy-math-inline\" data-latex=\"f_{X, \\Theta}(x, \\theta)\">$f_{X, \\Theta}(x, \\theta)$<\/span> sur l&rsquo;espace des param\u00e8tres.<\/p>\n<p>Dans le domaine des sciences cognitives computationnelles, ce cadre bivari\u00e9 r\u00e9git l&rsquo;hypoth\u00e8se du \u00ab cerveau bay\u00e9sien \u00bb. On y mod\u00e9lise la perception comme le calcul continu d&rsquo;une distribution a posteriori sur les causes environnementales distales <span class=\"arabpsy-math-inline\" data-latex=\"Theta\">$Theta$<\/span>, sachant les signaux sensoriels conjoints proximaux <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> capt\u00e9s par les r\u00e9cepteurs biologiques. L&rsquo;exactitude des inf\u00e9rences et des comportements d\u00e9cisionnels d\u00e9coule ainsi de la fid\u00e9lit\u00e9 avec laquelle les circuits neuronaux parviennent \u00e0 manipuler la loi conjointe r\u00e9gissant les interactions entre l&rsquo;organisme et son milieu.<\/p>\n<h2>7. Ind\u00e9pendance stochastique au prisme de la loi conjointe<\/h2>\n<h3>7.1 Crit\u00e8re formel de factorisation de la loi conjointe<\/h3>\n<p>Le concept d&rsquo;<strong>ind\u00e9pendance stochastique<\/strong> est l&rsquo;un des piliers th\u00e9oriques les plus d\u00e9terminants des probabilit\u00e9s multivari\u00e9es. Deux variables al\u00e9atoires <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> sont dites mutuellement ind\u00e9pendantes au sens probabiliste si et seulement si leur fonction de distribution conjointe se factorise exactement en le produit pur de leurs distributions marginales respectives. Pour des variables discr\u00e8tes, cette condition imp\u00e9rative s&rsquo;\u00e9nonce :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p_{X, Y}(x, y) = p_X(x) \\cdot p_Y(y), \\quad \\forall (x, y) in \\mathcal{X} \\times \\mathcal{Y}\">$$p_{X, Y}(x, y) = p_X(x) \\cdot p_Y(y), \\quad \\forall (x, y) in \\mathcal{X} \\times \\mathcal{Y}$$<\/span><\/p>\n<p>Dans le domaine continu, la factorisation s&rsquo;applique de mani\u00e8re isomorphe sur la fonction de densit\u00e9 conjointe presque partout :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_{X, Y}(x, y) = f_X(x) \\cdot f_Y(y), \\quad \\forall (x, y) in \\mathbb{R}^2\">$$f_{X, Y}(x, y) = f_X(x) \\cdot f_Y(y), \\quad \\forall (x, y) in \\mathbb{R}^2$$<\/span><\/p>\n<p>Ce crit\u00e8re multiplicatif entra\u00eene une cons\u00e9quence \u00e9pist\u00e9mique fondamentale : l&rsquo;ind\u00e9pendance stochastique \u00e9quivaut \u00e0 affirmer que l&rsquo;information apport\u00e9e par une variable n&rsquo;alt\u00e8re en rien la structure probabiliste de la seconde. Math\u00e9matiquement, cette stabilit\u00e9 se traduit par la superposition parfaite entre distributions conditionnelles et distributions marginales :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_{X mid Y}(x mid y) = \\frac{f(x, y)}{f_Y(y)} = \\frac{f_X(x) \\cdot f_Y(y)}{f_Y(y)} = f_X(x)\">$$f_{X mid Y}(x mid y) = \\frac{f(x, y)}{f_Y(y)} = \\frac{f_X(x) \\cdot f_Y(y)}{f_Y(y)} = f_X(x)$$<\/span><\/p>\n<p>Au niveau de la fonction de r\u00e9partition cumulative bivari\u00e9e, l&rsquo;ind\u00e9pendance impose sym\u00e9triquement une s\u00e9parabilit\u00e9 multiplicative totale : <span class=\"arabpsy-math-inline\" data-latex=\"F_{X, Y}(x, y) = F_X(x) \\cdot F_Y(y)\">$F_{X, Y}(x, y) = F_X(x) \\cdot F_Y(y)$<\/span>. L&rsquo;existence d&rsquo;une telle factorisation garantit que l&rsquo;esp\u00e9rance du produit crois\u00e9 se r\u00e9duit au produit des esp\u00e9rances individuelles :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbb{E}[XY] = \\iint xy , f(x, y) , dx , dy = \\left(\\int x f_X(x) , dx\\right) \\left(\\int y f_Y(y) , dy\\right) = \\mathbb{E}[X] \\cdot \\mathbb{E}[Y]\">$$\\mathbb{E}[XY] = \\iint xy , f(x, y) , dx , dy = \\left(\\int x f_X(x) , dx\\right) \\left(\\int y f_Y(y) , dy\\right) = \\mathbb{E}[X] \\cdot \\mathbb{E}[Y]$$<\/span><\/p>\n<h3>7.2 Distinction cruciale entre ind\u00e9pendance et d\u00e9corr\u00e9lation<\/h3>\n<p>L&rsquo;une des erreurs conceptuelles les plus tenaces et dommageables dans la pratique de l&rsquo;analyse quantitative r\u00e9side dans l&rsquo;assimilation h\u00e2tive de la d\u00e9corr\u00e9lation \u00e0 l&rsquo;ind\u00e9pendance stochastique. Il existe une asym\u00e9trie logique stricte entre ces deux notions : <em>l&rsquo;ind\u00e9pendance stochastique implique syst\u00e9matiquement une d\u00e9corr\u00e9lation lin\u00e9aire, mais la r\u00e9ciproque est g\u00e9n\u00e9ralement fausse<\/em>.<\/p>\n<p>Pour \u00e9tablir la premi\u00e8re implication, calculons la covariance entre deux variables stochastiquement ind\u00e9pendantes <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathrm{Cov}(X, Y) = \\mathbb{E}[XY] - \\mathbb{E}[X]\\mathbb{E}[Y] = \\mathbb{E}[X]\\mathbb{E}[Y] - \\mathbb{E}[X]\\mathbb{E}[Y] = 0\">$$\\mathrm{Cov}(X, Y) = \\mathbb{E}[XY] &#8211; \\mathbb{E}[X]\\mathbb{E}[Y] = \\mathbb{E}[X]\\mathbb{E}[Y] &#8211; \\mathbb{E}[X]\\mathbb{E}[Y] = 0$$<\/span><\/p>\n<p>La covariance \u00e9tant nulle, le coefficient de corr\u00e9lation lin\u00e9aire s&rsquo;annule obligatoirement (<span class=\"arabpsy-math-inline\" data-latex=\"\\rho_{X, Y} = 0\">$\\rho_{X, Y} = 0$<\/span>).<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;annulation de la corr\u00e9lation ne garantit en aucun cas l&rsquo;ind\u00e9pendance stochastique. Cette dissociation se d\u00e9montre \u00e9l\u00e9gamment au moyen d&rsquo;un contre-exemple analytique classique. Soit <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> une variable al\u00e9atoire continue sym\u00e9trique centr\u00e9e autour de z\u00e9ro, dont la densit\u00e9 v\u00e9rifie <span class=\"arabpsy-math-inline\" data-latex=\"f_X(x) = f_X(-x)\">$f_X(x) = f_X(-x)$<\/span> pour tout <span class=\"arabpsy-math-inline\" data-latex=\"x in [-1, 1]\">$x in [-1, 1]$<\/span>, et d\u00e9finissons la variable <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> par une relation fonctionnelle d\u00e9terministe non lin\u00e9aire : <span class=\"arabpsy-math-inline\" data-latex=\"Y = X^2\">$Y = X^2$<\/span>. Par construction, <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> est totalement d\u00e9pendante de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> : conna\u00eetre <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> permet de pr\u00e9dire la valeur exacte de <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> avec une erreur nulle. \u00c9valuons n\u00e9anmoins leur covariance :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbb{E}[XY] = \\mathbb{E}[X \\cdot X^2] = \\mathbb{E}[X^3] = \\int_{-1}^{1} x^3 f_X(x) , dx = 0\">$$\\mathbb{E}[XY] = \\mathbb{E}[X \\cdot X^2] = \\mathbb{E}[X^3] = \\int_{-1}^{1} x^3 f_X(x) , dx = 0$$<\/span><\/p>\n<p>Puisque la fonction <span class=\"arabpsy-math-inline\" data-latex=\"x^3 f_X(x)\">$x^3 f_X(x)$<\/span> est une fonction impaire int\u00e9gr\u00e9e sur un intervalle sym\u00e9trique, l&rsquo;esp\u00e9rance crois\u00e9e s&rsquo;annule. Par suite, <span class=\"arabpsy-math-inline\" data-latex=\"\\mathrm{Cov}(X, Y) = \\mathbb{E}[XY] - \\mathbb{E}[X]\\mathbb{E}[Y] = 0 - 0 = 0\">$\\mathrm{Cov}(X, Y) = \\mathbb{E}[XY] &#8211; \\mathbb{E}[X]\\mathbb{E}[Y] = 0 &#8211; 0 = 0$<\/span>. Le coefficient de corr\u00e9lation de Pearson entre ces deux variables est rigoureusement nul (<span class=\"arabpsy-math-inline\" data-latex=\"\\rho_{X, Y} = 0\">$\\rho_{X, Y} = 0$<\/span>).<\/p>\n<p>Cet exemple atteste avec \u00e9clat que la corr\u00e9lation lin\u00e9aire ne mesure que l&rsquo;alignement rectiligne des variations. La distribution conjointe <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y)\">$f(x, y)$<\/span> h\u00e9berge quant \u00e0 elle l&rsquo;int\u00e9gralit\u00e9 des d\u00e9pendances d&rsquo;ordres sup\u00e9rieurs, qu&rsquo;elles soient quadratiques, exponentielles ou sinuso\u00efdales. Se fier exclusivement \u00e0 un indicateur lin\u00e9aire pour conclure \u00e0 l&rsquo;absence de lien au sein d&rsquo;une distribution conjointe constitue une impasse analytique majeure.<\/p>\n<h3>7.3 Tests statistiques d&rsquo;ind\u00e9pendance fond\u00e9s sur la distribution conjointe<\/h3>\n<p>La validation empirique de l&rsquo;hypoth\u00e8se d&rsquo;ind\u00e9pendance stochastique n\u00e9cessite la mise en \u0153uvre de tests d&rsquo;hypoth\u00e8ses rigoureux comparant la distribution conjointe observ\u00e9e \u00e0 la distribution th\u00e9orique sous condition d&rsquo;ind\u00e9pendance. Le cadre paradigmatique pour les variables qualitatives ou discr\u00e8tes est le <strong>test du chi-deux (<span class=\"arabpsy-math-inline\" data-latex=\"\\chi^2\">$\\chi^2$<\/span>) d&rsquo;ind\u00e9pendance<\/strong> de Karl Pearson.<\/p>\n<p>Ce test \u00e9value l&rsquo;\u00e9cart quadratique entre les effectifs conjoints observ\u00e9s <span class=\"arabpsy-math-inline\" data-latex=\"n_{ij}\">$n_{ij}$<\/span> au sein du tableau de contingence et les effectifs th\u00e9oriques attendus sous l&rsquo;hypoth\u00e8se nulle <span class=\"arabpsy-math-inline\" data-latex=\"H_0 : p_{ij} = p_{i\\cdot} \\cdot p_{\\cdot j}\">$H_0 : p_{ij} = p_{i\\cdot} \\cdot p_{\\cdot j}$<\/span>. L&rsquo;effectif attendu dans la cellule <span class=\"arabpsy-math-inline\" data-latex=\"(i, j)\">$(i, j)$<\/span> s&rsquo;exprime par :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"e_{ij} = \\frac{n_{i\\cdot} \\cdot n_{\\cdot j}}{N}\">$$e_{ij} = \\frac{n_{i\\cdot} \\cdot n_{\\cdot j}}{N}$$<\/span><\/p>\n<p>La statistique de test s&rsquo;\u00e9crit alors :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\chi^2 = \\sum_{i=1}^r \\sum_{j=1}^c \\frac{(n_{ij} - e_{ij})^2}{e_{ij}}\">$$\\chi^2 = \\sum_{i=1}^r \\sum_{j=1}^c \\frac{(n_{ij} &#8211; e_{ij})^2}{e_{ij}}$$<\/span><\/p>\n<p>Sous <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span> et sous r\u00e9serve du respect des crit\u00e8res asymptotiques de Cochran (effectifs th\u00e9oriques <span class=\"arabpsy-math-inline\" data-latex=\"e_{ij} ge 5\">$e_{ij} ge 5$<\/span> dans la majorit\u00e9 des cellules), cette grandeur converge asymptotiquement vers une loi du <span class=\"arabpsy-math-inline\" data-latex=\"\\chi^2\">$\\chi^2$<\/span> \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"(r - 1)(c - 1)\">$(r &#8211; 1)(c &#8211; 1)$<\/span> degr\u00e9s de libert\u00e9.<\/p>\n<p>Dans le champ de l&rsquo;inf\u00e9rence moderne et de la th\u00e9orie de l&rsquo;information, le test du rapport de vraisemblance (ou statistique <span class=\"arabpsy-math-inline\" data-latex=\"G^2\">$G^2$<\/span>) offre une alternative robuste. Il entretient un lien direct avec la <strong>divergence de Kullback-Leibler<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"D_{\\mathrm{KL}}\">$D_{\\mathrm{KL}}$<\/span>, qui quantifie la distance relative entre la distribution conjointe empirique <span class=\"arabpsy-math-inline\" data-latex=\"P_{X, Y}\">$P_{X, Y}$<\/span> et le mod\u00e8le produit <span class=\"arabpsy-math-inline\" data-latex=\"P_X o\\times P_Y\">$P_X o\\times P_Y$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"I(X; Y) = D_{\\mathrm{KL}}(P_{X, Y} parallel P_X o\\times P_Y) = \\sum_{x in \\mathcal{X}} \\sum_{y in \\mathcal{Y}} p(x, y) \\log \\left( \\frac{p(x, y)}{p_X(x) p_Y(y)} \\right)\">$$I(X; Y) = D_{\\mathrm{KL}}(P_{X, Y} parallel P_X o\\times P_Y) = \\sum_{x in \\mathcal{X}} \\sum_{y in \\mathcal{Y}} p(x, y) \\log \\left( \\frac{p(x, y)}{p_X(x) p_Y(y)} \\right)$$<\/span><\/p>\n<p>Cette mesure, d\u00e9sign\u00e9e sous le nom d&rsquo;<strong>information mutuelle<\/strong>, est strictement positive ou nulle, ne s&rsquo;annulant si et seulement si l&rsquo;ind\u00e9pendance stochastique est rigoureusement v\u00e9rifi\u00e9e sur l&rsquo;ensemble du domaine conjointe. Elle transcende les contraintes de lin\u00e9arit\u00e9 pour capter la totalit\u00e9 des structures d&rsquo;association discr\u00e8tes ou continues.<\/p>\n<h2>8. Moments bivari\u00e9s : Esp\u00e9rance conjointe, covariance et corr\u00e9lation<\/h2>\n<h3>8.1 Esp\u00e9rance math\u00e9matique d&rsquo;une fonction de variables conjointes<\/h3>\n<p>L&rsquo;analyse des distributions conjointes s&rsquo;\u00e9tend naturellement au calcul des caract\u00e9ristiques num\u00e9riques r\u00e9sumant le comportement moyen des syst\u00e8mes bivari\u00e9s. Pour \u00e9valuer l&rsquo;esp\u00e9rance math\u00e9matique d&rsquo;une fonction d\u00e9terministe arbitraire <span class=\"arabpsy-math-inline\" data-latex=\"g(X, Y)\">$g(X, Y)$<\/span> appliqu\u00e9e au couple stochastique, la th\u00e9orie statistique mobilise le <strong>th\u00e9or\u00e8me de l&rsquo;esp\u00e9rance non consciente<\/strong> (connu sous l&rsquo;acronyme anglophone <em>LOTUS<\/em> pour <em>Law of the Unconscious Statistician<\/em>). Ce r\u00e9sultat majeur dispense d&rsquo;\u00e9tablir pr\u00e9alablement la loi univari\u00e9e de la variable transform\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"Z = g(X, Y)\">$Z = g(X, Y)$<\/span>, permettant de calculer l&rsquo;esp\u00e9rance directement \u00e0 partir de la fonction conjointe de d\u00e9part.<\/p>\n<p>Dans le cas discret, cette esp\u00e9rance correspond \u00e0 la sommation pond\u00e9r\u00e9e par la masse de chaque point du support :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbb{E}[g(X, Y)] = \\sum_{x in \\mathcal{X}} \\sum_{y in \\mathcal{Y}} g(x, y) \\cdot p_{X, Y}(x, y)\">$$\\mathbb{E}[g(X, Y)] = \\sum_{x in \\mathcal{X}} \\sum_{y in \\mathcal{Y}} g(x, y) \\cdot p_{X, Y}(x, y)$$<\/span><\/p>\n<p>Dans le cas continu, l&rsquo;expression s&rsquo;\u00e9nonce sous forme d&rsquo;une int\u00e9grale surfacique bivari\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbb{E}[g(X, Y)] = \\int_{-\\infty}^{\\infty} \\int_{-\\infty}^{\\infty} g(x, y) \\cdot f_{X, Y}(x, y) , dx , dy\">$$\\mathbb{E}[g(X, Y)] = \\int_{-\\infty}^{\\infty} \\int_{-\\infty}^{\\infty} g(x, y) \\cdot f_{X, Y}(x, y) , dx , dy$$<\/span><\/p>\n<p>Une propri\u00e9t\u00e9 fondamentale de cet op\u00e9rateur r\u00e9side dans la lin\u00e9arit\u00e9 stricte de l&rsquo;esp\u00e9rance. Quelle que soit la nature de la loi conjointe et sans exiger la moindre hypoth\u00e8se d&rsquo;ind\u00e9pendance, l&rsquo;esp\u00e9rance d&rsquo;une combinaison lin\u00e9aire de variables al\u00e9atoires \u00e9quivaut invariablement \u00e0 la combinaison lin\u00e9aire de leurs esp\u00e9rances individuelles :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbb{E}[aX + bY + c] = a\\mathbb{E}[X] + b\\mathbb{E}[Y] + c, \\quad \\forall (a, b, c) in \\mathbb{R}^3\">$$\\mathbb{E}[aX + bY + c] = a\\mathbb{E}[X] + b\\mathbb{E}[Y] + c, \\quad \\forall (a, b, c) in \\mathbb{R}^3$$<\/span><\/p>\n<p>En th\u00e9orie hilbertienne des espaces probabilis\u00e9s, l&rsquo;esp\u00e9rance conditionnelle <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{E}[X mid Y = y]\">$\\mathbb{E}[X mid Y = y]$<\/span> s&rsquo;interpr\u00e8te comme la projection orthogonale de la variable al\u00e9atoire <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> sur le sous-espace ferm\u00e9 des fonctions de carr\u00e9 int\u00e9grable <span class=\"arabpsy-math-inline\" data-latex=\"L^2(\\Omega)\">$L^2(\\Omega)$<\/span> engendr\u00e9 par la variable <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>. Cette formulation g\u00e9om\u00e9trique \u00e9tablit que la fonction conditionnelle <span class=\"arabpsy-math-inline\" data-latex=\"\\mu(Y) = \\mathbb{E}[X mid Y]\">$\\mu(Y) = \\mathbb{E}[X mid Y]$<\/span> constitue le meilleur pr\u00e9dicteur de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> au sens de la minimisation de l&rsquo;erreur quadratique moyenne.<\/p>\n<h3>8.2 La covariance comme mesure du moment conjoint centr\u00e9<\/h3>\n<p>Pour appr\u00e9hender l&rsquo;amplitude de la covariation lin\u00e9aire conjointe de deux grandeurs, on introduit le deuxi\u00e8me moment conjoint centr\u00e9, universellement d\u00e9sign\u00e9 sous le terme de <strong>covariance<\/strong>. Formellement, la covariance entre les variables al\u00e9atoires <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> se d\u00e9finit comme l&rsquo;esp\u00e9rance du produit de leurs \u00e9carts respectifs \u00e0 leurs moyennes marginales :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathrm{Cov}(X, Y) = \\sigma_{XY} = \\mathbb{E}[(X - \\mathbb{E}[X])(Y - \\mathbb{E}[Y])]\">$$\\mathrm{Cov}(X, Y) = \\sigma_{XY} = \\mathbb{E}[(X &#8211; \\mathbb{E}[X])(Y &#8211; \\mathbb{E}[Y])]$$<\/span><\/p>\n<p>En d\u00e9ployant cette expression par bilin\u00e9arit\u00e9 de l&rsquo;esp\u00e9rance, on obtient la formule calculatoire canonique (identit\u00e9 de Huygens-K\u00f6nig bivari\u00e9e) :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathrm{Cov}(X, Y) = \\mathbb{E}[XY] - \\mathbb{E}[X]\\mathbb{E}[Y]\">$$\\mathrm{Cov}(X, Y) = \\mathbb{E}[XY] &#8211; \\mathbb{E}[X]\\mathbb{E}[Y]$$<\/span><\/p>\n<p>Sur le plan interpr\u00e9tatif, la covariance prend une valeur positive lorsque les r\u00e9alisations sup\u00e9rieures \u00e0 la moyenne de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> ont tendance \u00e0 s&rsquo;associer pr\u00e9f\u00e9rentiellement aux r\u00e9alisations sup\u00e9rieures \u00e0 la moyenne de <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>, et inversement pour les valeurs inf\u00e9rieures. \u00c0 l&rsquo;oppos\u00e9, si les valeurs \u00e9lev\u00e9es de l&rsquo;une co\u00efncident de mani\u00e8re pr\u00e9pond\u00e9rante avec les valeurs basses de l&rsquo;autre, le produit des \u00e9carts est majoritairement n\u00e9gatif, induisant une covariance globale n\u00e9gative.<\/p>\n<p>La covariance ob\u00e9it \u00e0 des propri\u00e9t\u00e9s alg\u00e9briques rigoureuses : elle est sym\u00e9trique (<span class=\"arabpsy-math-inline\" data-latex=\"\\mathrm{Cov}(X, Y) = \\mathrm{Cov}(Y, X)\">$\\mathrm{Cov}(X, Y) = \\mathrm{Cov}(Y, X)$<\/span>), bilin\u00e9aire, et sa valeur appliqu\u00e9e \u00e0 une variable avec elle-m\u00eame restitue imm\u00e9diatement sa variance univari\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"\\mathrm{Cov}(X, X) = \\mathrm{Var}(X) = \\sigma_X^2\">$\\mathrm{Cov}(X, X) = \\mathrm{Var}(X) = \\sigma_X^2$<\/span>). La variance d&rsquo;une somme de variables s&rsquo;exprime ainsi directement par l&rsquo;int\u00e9gration de leur terme de d\u00e9pendance conjointe :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathrm{Var}(aX + bY) = a^2 \\mathrm{Var}(X) + b^2 \\mathrm{Var}(Y) + 2ab , \\mathrm{Cov}(X, Y)\">$$\\mathrm{Var}(aX + bY) = a^2 \\mathrm{Var}(X) + b^2 \\mathrm{Var}(Y) + 2ab , \\mathrm{Cov}(X, Y)$$<\/span><\/p>\n<p>N\u00e9anmoins, la covariance souffre d&rsquo;une limitation m\u00e9thodologique intrins\u00e8que : sa magnitude d\u00e9pend strictement des unit\u00e9s d&rsquo;\u00e9chelle adopt\u00e9es pour mesurer <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>. Si l&rsquo;on multiplie l&rsquo;une des grandeurs par un facteur d&rsquo;\u00e9chelle (par exemple, en convertissant des secondes en millisecondes), la covariance se trouve m\u00e9caniquement multipli\u00e9e par ce m\u00eame facteur, interdisant toute comparaison directe d&rsquo;intensit\u00e9 d&rsquo;association entre des couples de variables distincts.<\/p>\n<h3>8.3 Le coefficient de corr\u00e9lation lin\u00e9aire de Pearson<\/h3>\n<p>Pour affranchir la mesure d&rsquo;association lin\u00e9aire des contingences dimensionnelles et des choix d&rsquo;unit\u00e9s de mesure, Karl Pearson a formalis\u00e9 le <strong>coefficient de corr\u00e9lation lin\u00e9aire<\/strong> bivari\u00e9, conventionnellement not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\rho_{X, Y}\">$\\rho_{X, Y}$<\/span> ou plus sobrement <span class=\"arabpsy-math-inline\" data-latex=\"r\">$r$<\/span>. Ce coefficient r\u00e9sulte de la normalisation rigoureuse de la covariance par le produit des \u00e9carts-types marginaux :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\rho_{X, Y} = \\frac{\\mathrm{Cov}(X, Y)}{\\sigma_X \\sigma_Y} = \\frac{\\mathbb{E}[(X - \\mu_X)(Y - \\mu_Y)]}{\\sqrt{\\mathbb{E}[(X - \\mu_X)^2]} \\sqrt{\\mathbb{E}[(Y - \\mu_Y)^2]}}\">$$\\rho_{X, Y} = \\frac{\\mathrm{Cov}(X, Y)}{\\sigma_X \\sigma_Y} = \\frac{\\mathbb{E}[(X &#8211; \\mu_X)(Y &#8211; \\mu_Y)]}{\\sqrt{\\mathbb{E}[(X &#8211; \\mu_X)^2]} \\sqrt{\\mathbb{E}[(Y &#8211; \\mu_Y)^2]}}$$<\/span><\/p>\n<p>En vertu de l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/In%C3%A9galit%C3%A9_de_Cauchy-Schwarz\">in\u00e9galit\u00e9 de Cauchy-Schwarz<\/a> appliqu\u00e9e aux variables al\u00e9atoires centr\u00e9es consid\u00e9r\u00e9es comme des vecteurs dans l&rsquo;espace pr\u00e9hilbertien <span class=\"arabpsy-math-inline\" data-latex=\"L^2\">$L^2$<\/span>, la valeur de <span class=\"arabpsy-math-inline\" data-latex=\"rho\">$rho$<\/span> est th\u00e9oriquement born\u00e9e dans l&rsquo;intervalle compact ferm\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"[-1, 1]\">$[-1, 1]$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"-1 le \\rho_{X, Y} le 1\">$$-1 le \\rho_{X, Y} le 1$$<\/span><\/p>\n<p>Les bornes extr\u00eames correspondent \u00e0 des situations limites d&rsquo;une clart\u00e9 g\u00e9om\u00e9trique absolue : <span class=\"arabpsy-math-inline\" data-latex=\"\\rho_{X, Y} = 1\">$\\rho_{X, Y} = 1$<\/span> indique une colin\u00e9arit\u00e9 d\u00e9terministe positive exacte (<span class=\"arabpsy-math-inline\" data-latex=\"Y = aX + b\">$Y = aX + b$<\/span> avec <span class=\"arabpsy-math-inline\" data-latex=\"a &gt; 0\">$a &gt; 0$<\/span>), tandis que <span class=\"arabpsy-math-inline\" data-latex=\"\\rho_{X, Y} = -1\">$\\rho_{X, Y} = -1$<\/span> caract\u00e9rise une relation lin\u00e9aire d\u00e9croissante parfaite (<span class=\"arabpsy-math-inline\" data-latex=\"a &lt; 0\">$a &lt; 0$<\/span>). Une valeur <span class=\"arabpsy-math-inline\" data-latex=\"\\rho_{X, Y} = 0\">$\\rho_{X, Y} = 0$<\/span> atteste quant \u00e0 elle d&#039;une d\u00e9corr\u00e9lation lin\u00e9aire compl\u00e8te au sein du couple bivari\u00e9.<\/p>\n<p>Le coefficient de corr\u00e9lation pr\u00e9sente l&rsquo;avantage majeur d&rsquo;\u00eatre invariant face \u00e0 toute transformation affine strictement croissante des variables d&rsquo;origine : si <span class=\"arabpsy-math-inline\" data-latex=\"X&#039; = aX + b\">$X&rsquo; = aX + b$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y&#039; = cY + d\">$Y&rsquo; = cY + d$<\/span> avec <span class=\"arabpsy-math-inline\" data-latex=\"a &gt; 0\">$a &gt; 0$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"c &gt; 0\">$c &gt; 0$<\/span>, alors <span class=\"arabpsy-math-inline\" data-latex=\"\\rho_{X&#039;, Y&#039;} = \\rho_{X, Y}\">$\\rho_{X&rsquo;, Y&rsquo;} = \\rho_{X, Y}$<\/span>. Au sein d&rsquo;une distribution conjointe continue, le coefficient de corr\u00e9lation gouverne directement l&rsquo;ellipticit\u00e9 et l&rsquo;inclinaison des lignes d&rsquo;isodensit\u00e9 projet\u00e9es sur le plan orthogonal : plus <span class=\"arabpsy-math-inline\" data-latex=\"|rho|\">$|rho|$<\/span> converge vers 1, plus les contours concentriques s&rsquo;amincissent jusqu&rsquo;\u00e0 se fondre en un segment lin\u00e9aire unidimensionnel.<\/p>\n<h2>9. Distributions conjointes notables : Mod\u00e8les th\u00e9oriques de r\u00e9f\u00e9rence<\/h2>\n<h3>9.1 La distribution multinomiale : Extension conjointe de la loi binomiale<\/h3>\n<p>Parmi les distributions conjointes discr\u00e8tes r\u00e9gissant les processus combinatoires fondamentaux, la <strong>loi multinomiale<\/strong> repr\u00e9sente le mod\u00e8le de r\u00e9f\u00e9rence par excellence. Elle constitue la g\u00e9n\u00e9ralisation multivari\u00e9e directe de la loi binomiale au cas o\u00f9 chaque \u00e9preuve exp\u00e9rimentale ind\u00e9pendante admet non plus deux issues (succ\u00e8s\/\u00e9chec), mais un ensemble fini de <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> issues mutuellement exclusives et collectivement exhaustives, survenant avec des probabilit\u00e9s respectives <span class=\"arabpsy-math-inline\" data-latex=\"p_1, p_2, dots, p_k\">$p_1, p_2, dots, p_k$<\/span> telles que <span class=\"arabpsy-math-inline\" data-latex=\"\\sum_{j=1}^k p_j = 1\">$\\sum_{j=1}^k p_j = 1$<\/span>.<\/p>\n<p>Si l&rsquo;on r\u00e9it\u00e8re cette exp\u00e9rience al\u00e9atoire <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> fois de fa\u00e7on ind\u00e9pendante et stationnaire, le vecteur al\u00e9atoire <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X} = (X_1, X_2, dots, X_k)^T\">$\\mathbf{X} = (X_1, X_2, dots, X_k)^T$<\/span> consigne le nombre de r\u00e9alisations observ\u00e9es pour chacune des <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> cat\u00e9gories, sous la contrainte stochastique rigide <span class=\"arabpsy-math-inline\" data-latex=\"\\sum_{j=1}^k X_j = n\">$\\sum_{j=1}^k X_j = n$<\/span>. La fonction de masse conjointe de ce vecteur s&rsquo;\u00e9nonce par la formule multinomiale classique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(X_1 = x_1, X_2 = x_2, dots, X_k = x_k) = \\frac{n!}{x_1! x_2! \\cdots x_k!} p_1^{x_1} p_2^{x_2} \\cdots p_k^{x_k}\">$$P(X_1 = x_1, X_2 = x_2, dots, X_k = x_k) = \\frac{n!}{x_1! x_2! \\cdots x_k!} p_1^{x_1} p_2^{x_2} \\cdots p_k^{x_k}$$<\/span><\/p>\n<p>pour tout uplet d&rsquo;entiers naturels v\u00e9rifiant <span class=\"arabpsy-math-inline\" data-latex=\"\\sum_{j=1}^k x_j = n\">$\\sum_{j=1}^k x_j = n$<\/span>, et s&rsquo;annule pour toute autre configuration.<\/p>\n<p>L&rsquo;\u00e9tude des moments conjoints de la loi multinomiale illustre parfaitement les m\u00e9canismes de d\u00e9pendance stochastique induits par des contraintes de conservation globales. Si l&rsquo;esp\u00e9rance marginale de chaque comptage reproduit la moyenne binomiale standard (<span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{E}[X_j] = n p_j\">$\\mathbb{E}[X_j] = n p_j$<\/span>) et sa variance marginale prend la forme <span class=\"arabpsy-math-inline\" data-latex=\"\\mathrm{Var}(X_j) = n p_j (1 - p_j)\">$\\mathrm{Var}(X_j) = n p_j (1 &#8211; p_j)$<\/span>, la covariance crois\u00e9e entre deux composantes distinctes <span class=\"arabpsy-math-inline\" data-latex=\"i ne j\">$i ne j$<\/span> est structurellement n\u00e9gative :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathrm{Cov}(X_i, X_j) = -n p_i p_j &lt; 0\">$$\\mathrm{Cov}(X_i, X_j) = -n p_i p_j &lt; 0$$<\/span><\/p>\n<p>Ce signe n\u00e9gatif traduit le ph\u00e9nom\u00e8ne de comp\u00e9tition inh\u00e9rent \u00e0 un \u00e9chantillon de taille fix\u00e9e : toute r\u00e9alisation surnum\u00e9raire dans la modalit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> s&rsquo;op\u00e8re obligatoirement au d\u00e9triment statistique du nombre d&rsquo;observations disponibles pour les autres cat\u00e9gories. La loi multinomiale constitue le socle indispensable au codage et \u00e0 la mod\u00e9lisation statistique des r\u00e9ponses cat\u00e9gorielles multiples, des tables de contingence polydimensionnelles et des donn\u00e9es d&rsquo;annotation linguistique ou g\u00e9nomique.<\/p>\n<h3>9.2 La distribution normale bivari\u00e9e et multivari\u00e9e<\/h3>\n<p>Au panth\u00e9on des distributions continues multivari\u00e9es, la <strong>loi normale bivari\u00e9e<\/strong> (et son extension multivari\u00e9e g\u00e9n\u00e9rale) occupe une place h\u00e9g\u00e9monique, dict\u00e9e par la puissance r\u00e9gulatrice du th\u00e9or\u00e8me central limite multidimensionnel. Pour un couple continu <span class=\"arabpsy-math-inline\" data-latex=\"(X, Y)\">$(X, Y)$<\/span>, la distribution normale conjointe est enti\u00e8rement d\u00e9termin\u00e9e par un ensemble fini de cinq param\u00e8tres scalaires : les moyennes marginales <span class=\"arabpsy-math-inline\" data-latex=\"\\mu_X\">$\\mu_X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mu_Y\">$\\mu_Y$<\/span>, les \u00e9carts-types marginaux <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_X &gt; 0\">$\\sigma_X &gt; 0$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_Y &gt; 0\">$\\sigma_Y &gt; 0$<\/span>, et le coefficient de corr\u00e9lation lin\u00e9aire <span class=\"arabpsy-math-inline\" data-latex=\"rho in ]-1, 1[\">$rho in ]-1, 1[$<\/span>.<\/p>\n<p>La fonction de densit\u00e9 conjointe analytique se formule par l&rsquo;expression gaussienne bivari\u00e9e canonique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f(x, y) = \\frac{1}{2\\pi \\sigma_X \\sigma_Y \\sqrt{1 - \\rho^2}} \\exp\\left( -\\frac{1}{2(1 - \\rho^2)} \\left[ \\frac{(x - \\mu_X)^2}{\\sigma_X^2} - \\frac{2\\rho(x - \\mu_X)(y - \\mu_Y)}{\\sigma_X \\sigma_Y} + \\frac{(y - \\mu_Y)^2}{\\sigma_Y^2} \\right] \\right)\">$$f(x, y) = \\frac{1}{2\\pi \\sigma_X \\sigma_Y \\sqrt{1 &#8211; \\rho^2}} \\exp\\left( -\\frac{1}{2(1 &#8211; \\rho^2)} \\left[ \\frac{(x &#8211; \\mu_X)^2}{\\sigma_X^2} &#8211; \\frac{2\\rho(x &#8211; \\mu_X)(y &#8211; \\mu_Y)}{\\sigma_X \\sigma_Y} + \\frac{(y &#8211; \\mu_Y)^2}{\\sigma_Y^2} \\right] \\right)$$<\/span><\/p>\n<p>En notation matricielle compacte, pour un vecteur al\u00e9atoire <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X} = (X_1, dots, X_p)^T \\sim \\mathcal{N}_p(boldsymbol{\\mu}, boldsymbol{\\Sigma})\">$\\mathbf{X} = (X_1, dots, X_p)^T \\sim \\mathcal{N}_p(boldsymbol{\\mu}, boldsymbol{\\Sigma})$<\/span>, la densit\u00e9 s&rsquo;\u00e9crit de mani\u00e8re synth\u00e9tique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f(\\mathbf{x}) = \\frac{1}{(2\\pi)^{p\/2} |boldsymbol{\\Sigma}|^{1\/2}} \\exp\\left( -\\frac{1}{2} (\\mathbf{x} - boldsymbol{\\mu})^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} - boldsymbol{\\mu}) \\right)\">$$f(\\mathbf{x}) = \\frac{1}{(2\\pi)^{p\/2} |boldsymbol{\\Sigma}|^{1\/2}} \\exp\\left( -\\frac{1}{2} (\\mathbf{x} &#8211; boldsymbol{\\mu})^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} &#8211; boldsymbol{\\mu}) \\right)$$<\/span><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}\">$boldsymbol{\\mu}$<\/span> est le vecteur des esp\u00e9rances math\u00e9matiques et <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span> d\u00e9signe la matrice de variance-covariance sym\u00e9trique d\u00e9finie positive de dimension <span class=\"arabpsy-math-inline\" data-latex=\"p \\times p\">$p \\times p$<\/span>, dont <span class=\"arabpsy-math-inline\" data-latex=\"|boldsymbol{\\Sigma}|\">$|boldsymbol{\\Sigma}|$<\/span> repr\u00e9sente le d\u00e9terminant.<\/p>\n<p>La loi normale bivari\u00e9e poss\u00e8de des propri\u00e9t\u00e9s analytiques singuli\u00e8res d&rsquo;une \u00e9l\u00e9gance remarquable :<\/p>\n<ul>\n<li>Toute combinaison lin\u00e9aire univari\u00e9e d&rsquo;un vecteur gaussien conjoint est elle-m\u00eame rigoureusement gaussienne.<\/li>\n<li>Les distributions marginales d\u00e9riv\u00e9es d&rsquo;une gaussienne multivari\u00e9e sont purement gaussiennes : <span class=\"arabpsy-math-inline\" data-latex=\"X \\sim \\mathcal{N}(\\mu_X, \\sigma_X^2)\">$X \\sim \\mathcal{N}(\\mu_X, \\sigma_X^2)$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y \\sim \\mathcal{N}(\\mu_Y, \\sigma_Y^2)\">$Y \\sim \\mathcal{N}(\\mu_Y, \\sigma_Y^2)$<\/span>.<\/li>\n<li>Toutes les lois conditionnelles sont rigoureusement gaussiennes univari\u00e9es, et leurs moments s&rsquo;expriment par des fonctions affines simples :<br \/>\n  <span class=\"arabpsy-math-display\" data-latex=\"\\mathbb{E}[Y mid X = x] = \\mu_Y + \\rho \\frac{\\sigma_Y}{\\sigma_X}(x - \\mu_X)\">$$\\mathbb{E}[Y mid X = x] = \\mu_Y + \\rho \\frac{\\sigma_Y}{\\sigma_X}(x &#8211; \\mu_X)$$<\/span><br \/>\n  <span class=\"arabpsy-math-display\" data-latex=\"\\mathrm{Var}(Y mid X = x) = \\sigma_Y^2 (1 - \\rho^2)\">$$\\mathrm{Var}(Y mid X = x) = \\sigma_Y^2 (1 &#8211; \\rho^2)$$<\/span><\/li>\n<li><strong>\u00c9quivalence exceptionnelle :<\/strong> Dans le cadre exclusif d&rsquo;une loi normale conjointe, <em>la d\u00e9corr\u00e9lation (<span class=\"arabpsy-math-inline\" data-latex=\"rho = 0\">$rho = 0$<\/span>) \u00e9quivaut strictement et sans r\u00e9serve \u00e0 l&rsquo;ind\u00e9pendance stochastique mutuelle<\/em>.<\/li>\n<\/ul>\n<p>La topologie de cette distribution se traduit par des ellipso\u00efdes de confiance concentriques : les lieux de niveau d&rsquo;\u00e9gale densit\u00e9 dessinent des ellipses sur le plan r\u00e9el, dont le grand axe s&rsquo;aligne selon la premi\u00e8re composante principale de dispersion du syst\u00e8me.<\/p>\n<h3>9.3 Th\u00e9orie des copules et th\u00e9or\u00e8me de Sklar<\/h3>\n<p>Bien que la distribution normale conjointe offre un confort analytique exceptionnel, elle impose une contrainte structurelle rigide : elle suppose implicitement que les distributions marginales sont toutes deux gaussiennes et que leur interd\u00e9pendance est parfaitement sym\u00e9trique et lin\u00e9aire. Confront\u00e9s \u00e0 des ph\u00e9nom\u00e8nes r\u00e9els pr\u00e9sentant des queues lourdes asym\u00e9triques ou des d\u00e9pendances non lin\u00e9aires complexes, les math\u00e9maticiens ont forg\u00e9 un outil r\u00e9volutionnaire : la <strong>th\u00e9orie des copules<\/strong>.<\/p>\n<p>Une copule est, par d\u00e9finition, une fonction de r\u00e9partition cumulative multivari\u00e9e d\u00e9finie sur l&rsquo;hypercube unit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]^d\">$[0, 1]^d$<\/span>, dont chaque loi marginale univari\u00e9e est uniforme sur l&rsquo;intervalle <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span>. Le socle de cette discipline repose int\u00e9gralement sur le c\u00e9l\u00e8bre <strong>th\u00e9or\u00e8me de Sklar<\/strong> (\u00e9tabli en 1959 par Abe Sklar). Ce th\u00e9or\u00e8me stipule que toute fonction de r\u00e9partition conjointe multivari\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"F(x_1, dots, x_d)\">$F(x_1, dots, x_d)$<\/span> reliant des variables aux marginales quelconques <span class=\"arabpsy-math-inline\" data-latex=\"F_1(x_1), dots, F_d(x_d)\">$F_1(x_1), dots, F_d(x_d)$<\/span> peut \u00eatre repr\u00e9sent\u00e9e sous la forme d&rsquo;une copule <span class=\"arabpsy-math-inline\" data-latex=\"C\">$C$<\/span> op\u00e9rant sur les probabilit\u00e9s cumul\u00e9es :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"F(x_1, x_2, dots, x_d) = C(F_1(x_1), F_2(x_2), dots, F_d(x_d))\">$$F(x_1, x_2, dots, x_d) = C(F_1(x_1), F_2(x_2), dots, F_d(x_d))$$<\/span><\/p>\n<p>Si les distributions marginales <span class=\"arabpsy-math-inline\" data-latex=\"F_i\">$F_i$<\/span> sont toutes rigoureusement continues, la copule <span class=\"arabpsy-math-inline\" data-latex=\"C\">$C$<\/span> associ\u00e9e est unique. Le th\u00e9or\u00e8me de Sklar r\u00e9alise ainsi une s\u00e9paration conceptuelle majeure : il permet de dissocier totalement l&rsquo;\u00e9tude des caract\u00e9ristiques marginales univari\u00e9es de chaque variable (leur forme, leur asym\u00e9trie, leurs moments) de l&rsquo;analyse de leur <em>structure de d\u00e9pendance pure<\/em>, celle-ci \u00e9tant enti\u00e8rement captur\u00e9e par la copule <span class=\"arabpsy-math-inline\" data-latex=\"C\">$C$<\/span>.<\/p>\n<p>Parmi les familles de copules abondamment mobilis\u00e9es en mod\u00e9lisation avanc\u00e9e, citons les copules gaussiennes et de Student (dites copules elliptiques), ainsi que les <strong>copules archim\u00e9diennes<\/strong> (telles que les copules de Clayton, Gumbel ou Frank). Ces derni\u00e8res sont d\u00e9finies au moyen d&rsquo;un g\u00e9n\u00e9rateur continu et poss\u00e8dent la propri\u00e9t\u00e9 cruciale de pouvoir mod\u00e9liser l&rsquo;<em>asym\u00e9trie de d\u00e9pendance dans les queues de distribution<\/em>. Une copule de Clayton, par exemple, permet de formaliser un comportement o\u00f9 deux variables deviennent hautement interd\u00e9pendantes lors d&rsquo;\u00e9pisodes de d\u00e9faillance extr\u00eame (effondrement boursier ou comorbidit\u00e9 clinique s\u00e9v\u00e8re), tout en demeurant quasi-ind\u00e9pendantes dans leur r\u00e9gime de variation mod\u00e9r\u00e9 habituel.<\/p>\n<h2>10. G\u00e9n\u00e9ralisation aux espaces multivari\u00e9s de dimension sup\u00e9rieure<\/h2>\n<h3>10.1 Vecteurs al\u00e9atoires et notation matricielle en dimension p<\/h3>\n<p>L&rsquo;appr\u00e9hension formelle de syst\u00e8mes complexes comportant des dizaines, voire des centaines de descripteurs stochastiques simultan\u00e9s exige une g\u00e9n\u00e9ralisation \u00e9l\u00e9gante de l&rsquo;alg\u00e8bre bivari\u00e9e vers l&rsquo;espace multidimensionnel <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^p\">$\\mathbb{R}^p$<\/span>. On d\u00e9finit alors un <strong>vecteur al\u00e9atoire<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X} = (X_1, X_2, dots, X_p)^T\">$\\mathbf{X} = (X_1, X_2, dots, X_p)^T$<\/span> comme une application mesurable reliant l&rsquo;espace probabilisable universel <span class=\"arabpsy-math-inline\" data-latex=\"(\\Omega, \\mathcal{F}, P)\">$(\\Omega, \\mathcal{F}, P)$<\/span> \u00e0 l&rsquo;espace euclidien muni de sa tribu bor\u00e9lienne <span class=\"arabpsy-math-inline\" data-latex=\"(\\mathbb{R}^p, \\mathcal{B}(\\mathbb{R}^p))\">$(\\mathbb{R}^p, \\mathcal{B}(\\mathbb{R}^p))$<\/span>.<\/p>\n<p>La fonction de r\u00e9partition conjointe multidimensionnelle s&rsquo;exprime par le cumul stochastique sur le c\u00f4ne n\u00e9gatif infini multidimensionnel :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"F(\\mathbf{x}) = F(x_1, x_2, dots, x_p) = P(X_1 le x_1, X_2 le x_2, dots, X_p le x_p)\">$$F(\\mathbf{x}) = F(x_1, x_2, dots, x_p) = P(X_1 le x_1, X_2 le x_2, dots, X_p le x_p)$$<\/span><\/p>\n<p>Pour un vecteur \u00e0 composantes continues, la fonction de densit\u00e9 de probabilit\u00e9 conjointe se formalise par la d\u00e9rivation d&rsquo;ordre <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> crois\u00e9e totale :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f(\\mathbf{x}) = \\frac{\\partial^p F(x_1, dots, x_p)}{\\partial x_1 \\partial x_2 \\cdots \\partial x_p}\">$$f(\\mathbf{x}) = \\frac{\\partial^p F(x_1, dots, x_p)}{\\partial x_1 \\partial x_2 \\cdots \\partial x_p}$$<\/span><\/p>\n<p>L&rsquo;architecture lin\u00e9aire des moments conjoints d&rsquo;ordre 1 et 2 se condense au sein de matrices aux propri\u00e9t\u00e9s spectrales fondamentales. L&rsquo;esp\u00e9rance math\u00e9matique vectorielle s&rsquo;\u00e9crit sous la forme du vecteur moyen :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"boldsymbol{\\mu} = \\mathbb{E}[\\mathbf{X}] = (\\mathbb{E}[X_1], \\mathbb{E}[X_2], dots, \\mathbb{E}[X_p])^T\">$$boldsymbol{\\mu} = \\mathbb{E}[\\mathbf{X}] = (\\mathbb{E}[X_1], \\mathbb{E}[X_2], dots, \\mathbb{E}[X_p])^T$$<\/span><\/p>\n<p>Tandis que la structure des dispersions conjointes s&rsquo;incarne dans la <strong>matrice de variance-covariance<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma} in \\mathbb{R}^{p \\times p}\">$boldsymbol{\\Sigma} in \\mathbb{R}^{p \\times p}$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"boldsymbol{\\Sigma} = \\mathrm{Var}(\\mathbf{X}) = \\mathbb{E}[(\\mathbf{X} - boldsymbol{\\mu})(\\mathbf{X} - boldsymbol{\\mu})^T] = \\begin{\\pmatrix} \\sigma_1^2 &amp;a\\mp; \\sigma_{12} &amp;a\\mp; \\cdots &amp;a\\mp; \\sigma_{1p}  \\sigma_{21} &amp;a\\mp; \\sigma_2^2 &amp;a\\mp; \\cdots &amp;a\\mp; \\sigma_{2p}  \\vdots &amp;a\\mp; \\vdots &amp;a\\mp; \\ddots &amp;a\\mp; \\vdots  \\sigma_{p1} &amp;a\\mp; \\sigma_{p2} &amp;a\\mp; \\cdots &amp;a\\mp; \\sigma_p^2 \\end{\\pmatrix}\">$$boldsymbol{\\Sigma} = \\mathrm{Var}(\\mathbf{X}) = \\mathbb{E}[(\\mathbf{X} &#8211; boldsymbol{\\mu})(\\mathbf{X} &#8211; boldsymbol{\\mu})^T] = \\begin{\\pmatrix} \\sigma_1^2 &#038;a\\mp; \\sigma_{12} &#038;a\\mp; \\cdots &#038;a\\mp; \\sigma_{1p}  \\sigma_{21} &#038;a\\mp; \\sigma_2^2 &#038;a\\mp; \\cdots &#038;a\\mp; \\sigma_{2p}  \\vdots &#038;a\\mp; \\vdots &#038;a\\mp; \\ddots &#038;a\\mp; \\vdots  \\sigma_{p1} &#038;a\\mp; \\sigma_{p2} &#038;a\\mp; \\cdots &#038;a\\mp; \\sigma_p^2 \\end{\\pmatrix}$$<\/span><\/p>\n<p>Cette matrice <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span> est par essence sym\u00e9trique (<span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}^T = boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}^T = boldsymbol{\\Sigma}$<\/span>) et semi-d\u00e9finie positive, ce qui garantit que pour tout vecteur d\u00e9terministe non nul <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{a} in \\mathbb{R}^p\">$\\mathbf{a} in \\mathbb{R}^p$<\/span>, la variance de la combinaison scalaire r\u00e9sultante demeure strictement non n\u00e9gative : <span class=\"arabpsy-math-inline\" data-latex=\"\\mathrm{Var}(\\mathbf{a}^T \\mathbf{X}) = \\mathbf{a}^T boldsymbol{\\Sigma} \\mathbf{a} ge 0\">$\\mathrm{Var}(\\mathbf{a}^T \\mathbf{X}) = \\mathbf{a}^T boldsymbol{\\Sigma} \\mathbf{a} ge 0$<\/span>. En pratique algorithmique et computationnelle, l&rsquo;exploitation de cette propri\u00e9t\u00e9 permet la <strong>d\u00e9composition de Cholesky<\/strong> de la matrice (<span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma} = \\mathbf{L}\\mathbf{L}^T\">$boldsymbol{\\Sigma} = \\mathbf{L}\\mathbf{L}^T$<\/span>), \u00e9tape indispensable pour simuler num\u00e9riquement des vecteurs al\u00e9atoires conjoints corr\u00e9l\u00e9s \u00e0 partir de simples g\u00e9n\u00e9rateurs gaussiens standard ind\u00e9pendants.<\/p>\n<h3>10.2 Marginalisation et conditionnement multidimensionnels<\/h3>\n<p>L&rsquo;un des accomplissements les plus raffin\u00e9s du calcul stochastique multivari\u00e9 r\u00e9side dans le partitionnement analytique des vecteurs al\u00e9atoires. Imaginons un vecteur al\u00e9atoire gaussien conjoint <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X} \\sim \\mathcal{N}_p(boldsymbol{\\mu}, boldsymbol{\\Sigma})\">$\\mathbf{X} \\sim \\mathcal{N}_p(boldsymbol{\\mu}, boldsymbol{\\Sigma})$<\/span> partitionn\u00e9 en deux sous-vecteurs de dimensions respectives <span class=\"arabpsy-math-inline\" data-latex=\"q\">$q$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"p-q\">$p-q$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{X} = \\begin{\\pmatrix} \\mathbf{X}_1  \\mathbf{X}_2 \\end{\\pmatrix}, \\quad boldsymbol{\\mu} = \\begin{\\pmatrix} boldsymbol{\\mu}_1  boldsymbol{\\mu}_2 \\end{\\pmatrix}, \\quad boldsymbol{\\Sigma} = \\begin{\\pmatrix} boldsymbol{\\Sigma}_{11} &amp;a\\mp; boldsymbol{\\Sigma}_{12}  boldsymbol{\\Sigma}_{21} &amp;a\\mp; boldsymbol{\\Sigma}_{22} \\end{\\pmatrix}\">$$\\mathbf{X} = \\begin{\\pmatrix} \\mathbf{X}_1  \\mathbf{X}_2 \\end{\\pmatrix}, \\quad boldsymbol{\\mu} = \\begin{\\pmatrix} boldsymbol{\\mu}_1  boldsymbol{\\mu}_2 \\end{\\pmatrix}, \\quad boldsymbol{\\Sigma} = \\begin{\\pmatrix} boldsymbol{\\Sigma}_{11} &#038;a\\mp; boldsymbol{\\Sigma}_{12}  boldsymbol{\\Sigma}_{21} &#038;a\\mp; boldsymbol{\\Sigma}_{22} \\end{\\pmatrix}$$<\/span><\/p>\n<p>L&rsquo;op\u00e9ration de marginalisation multidimensionnelle pr\u00e9sente une simplicit\u00e9 op\u00e9ratoire remarquable : la loi marginale conjointe du sous-vecteur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X}_1\">$\\mathbf{X}_1$<\/span> s&rsquo;obtient simplement en extrayant ses param\u00e8tres au sein de la partition globale, ignorant purement et simplement le bloc <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X}_2\">$\\mathbf{X}_2$<\/span> sans n\u00e9cessiter la moindre int\u00e9gration analytique fastidieuse :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{X}_1 \\sim \\mathcal{N}_q(boldsymbol{\\mu}_1, boldsymbol{\\Sigma}_{11})\">$$\\mathbf{X}_1 \\sim \\mathcal{N}_q(boldsymbol{\\mu}_1, boldsymbol{\\Sigma}_{11})$$<\/span><\/p>\n<p>Plus remarquable encore, la distribution conditionnelle conjointe multivari\u00e9e du sous-vecteur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X}_1\">$\\mathbf{X}_1$<\/span> sachant la fixation empirique du vecteur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X}_2 = \\mathbf{x}_2\">$\\mathbf{X}_2 = \\mathbf{x}_2$<\/span> demeure rigoureusement une distribution normale multidimensionnelle :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{X}_1 mid \\mathbf{X}_2 = \\mathbf{x}_2 \\sim \\mathcal{N}_q(boldsymbol{\\mu}_{1|2}, boldsymbol{\\Sigma}_{1|2})\">$$\\mathbf{X}_1 mid \\mathbf{X}_2 = \\mathbf{x}_2 \\sim \\mathcal{N}_q(boldsymbol{\\mu}_{1|2}, boldsymbol{\\Sigma}_{1|2})$$<\/span><\/p>\n<p>Les param\u00e8tres de cette loi conditionnelle se d\u00e9duisent analytiquement au moyen du <strong>compl\u00e9ment de Schur<\/strong> de la matrice de covariance :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"boldsymbol{\\mu}_{1|2} = boldsymbol{\\mu}_1 + boldsymbol{\\Sigma}_{12}boldsymbol{\\Sigma}_{22}^{-1}(\\mathbf{x}_2 - boldsymbol{\\mu}_2)\">$$boldsymbol{\\mu}_{1|2} = boldsymbol{\\mu}_1 + boldsymbol{\\Sigma}_{12}boldsymbol{\\Sigma}_{22}^{-1}(\\mathbf{x}_2 &#8211; boldsymbol{\\mu}_2)$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"boldsymbol{\\Sigma}_{1|2} = boldsymbol{\\Sigma}_{11} - boldsymbol{\\Sigma}_{12}boldsymbol{\\Sigma}_{22}^{-1}boldsymbol{\\Sigma}_{21}\">$$boldsymbol{\\Sigma}_{1|2} = boldsymbol{\\Sigma}_{11} &#8211; boldsymbol{\\Sigma}_{12}boldsymbol{\\Sigma}_{22}^{-1}boldsymbol{\\Sigma}_{21}$$<\/span><\/p>\n<p>Ce r\u00e9sultat th\u00e9orique profond sous-tend l&rsquo;int\u00e9gralit\u00e9 de la th\u00e9orie des mod\u00e8les d&rsquo;\u00e9quations structurelles, des s\u00e9ries temporelles multivari\u00e9es (filtres de Kalman) et de l&rsquo;analyse causale. L&rsquo;annulation d&rsquo;un coefficient sp\u00e9cifique dans la matrice de pr\u00e9cision conditionnelle inverse <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Theta} = boldsymbol{\\Sigma}^{-1}\">$boldsymbol{\\Theta} = boldsymbol{\\Sigma}^{-1}$<\/span> \u00e9quivaut formellement \u00e0 une hypoth\u00e8se d&rsquo;<strong>ind\u00e9pendance conditionnelle<\/strong> entre les deux variables concern\u00e9es \u00e9tant donn\u00e9 toutes les autres, mat\u00e9rialisant ainsi l&rsquo;absence d&rsquo;ar\u00eate directe dans un mod\u00e8le graphique gaussien non orient\u00e9 (Markov Random Field).<\/p>\n<h3>10.3 Le fl\u00e9au de la dimensionnalit\u00e9 en mod\u00e9lisation conjointe<\/h3>\n<p>D\u00e8s lors que l&rsquo;on s&rsquo;efforce d&rsquo;\u00e9tendre la mod\u00e9lisation conjointe \u00e0 un nombre \u00e9lev\u00e9 de dimensions, le statisticien se heurte \u00e0 un obstacle th\u00e9orique et pratique redoutable : le <strong>fl\u00e9au de la dimensionnalit\u00e9<\/strong> (expression originellement forg\u00e9e par Richard Bellman). Dans un espace discret, ce ph\u00e9nom\u00e8ne se traduit par une explosion combinatoire incontr\u00f4lable du nombre de cellules \u00e0 estimer. Si un tableau de contingence croise <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> variables binaires, le nombre d&rsquo;\u00e9tats conjoints distincts s&rsquo;\u00e9l\u00e8ve \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"2^p\">$2^p$<\/span>. Pour <span class=\"arabpsy-math-inline\" data-latex=\"p = 30\">$p = 30$<\/span> variables, le nombre de param\u00e8tres de masse libres exc\u00e8de un milliard (<span class=\"arabpsy-math-inline\" data-latex=\"2^{30} \\approx 1{,}07 \\times 10^9\">$2^{30} \\approx 1{,}07 \\times 10^9$<\/span>), rendant illusoire toute tentative d&rsquo;estimation empirique non r\u00e9gularis\u00e9e sur des cohortes humaines r\u00e9elles.<\/p>\n<p>Dans un espace continu, le fl\u00e9au prend une dimension g\u00e9om\u00e9trique : \u00e0 mesure que la dimensionnalit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> cro\u00eet, le volume de l&rsquo;espace euclidien augmente de mani\u00e8re exponentielle, provoquant une dispersion et une rar\u00e9faction drastique des observations empiriques. La quasi-totalit\u00e9 de la masse de probabilit\u00e9 conjointe se retrouve expuls\u00e9e vers les zones p\u00e9riph\u00e9riques et les coquilles externes de l&rsquo;espace, tandis que le centre g\u00e9om\u00e9trique se d\u00e9peuple enti\u00e8rement. En cons\u00e9quence, le nombre d&rsquo;\u00e9chantillons n\u00e9cessaires pour interpoler localement une densit\u00e9 conjointe au moyen d&rsquo;approches non param\u00e9triques cro\u00eet exponentiellement avec la dimension (<span class=\"arabpsy-math-inline\" data-latex=\"O(n^{-4\/(4+p)})\">$O(n^{-4\/(4+p)})$<\/span> pour un estimateur par noyau), d\u00e9truisant l&rsquo;efficacit\u00e9 statistique des estimateurs.<\/p>\n<p>Face \u00e0 cette impasse computationnelle, la recherche statistique s&rsquo;est tourn\u00e9e vers des strat\u00e9gies de parcimonie param\u00e9trique et de r\u00e9duction structurelle de dimension :<\/p>\n<ul>\n<li>L&rsquo;imposition de structures de d\u00e9pendance contraintes (matrices de covariance de rang faible, factorisation de Toeplitz ou mod\u00e8les autor\u00e9gressifs multivari\u00e9s).<\/li>\n<li>L&rsquo;adoption d&rsquo;a priori de r\u00e9gularisation bay\u00e9sienne et de p\u00e9nalisations par norme <span class=\"arabpsy-math-inline\" data-latex=\"ell_1\">$ell_1$<\/span> (Graphical Lasso) for\u00e7ant la matrice de pr\u00e9cision \u00e0 adopter un profil creux (sparse).<\/li>\n<li>Le d\u00e9ploiement de techniques de projection sous-spatiale, telles que l&rsquo;analyse en composantes principales (ACP) probabiliste ou les autoencodeurs variationnels, qui projettent la loi conjointe de haute dimension sur une vari\u00e9t\u00e9 latente de basse dimensionnalit\u00e9 pr\u00e9servant l&rsquo;information topologique essentielle.<\/li>\n<\/ul>\n<h2>11. Applications pratiques en psychom\u00e9trie et m\u00e9thodologie comportementale<\/h2>\n<h3>11.1 Mod\u00e9lisation conjointe des r\u00e9ponses aux items et des traits latents<\/h3>\n<p>L&rsquo;\u00e9dification des instruments psychom\u00e9triques modernes constitue l&rsquo;un des terrains d&rsquo;\u00e9lection les plus achev\u00e9s pour l&rsquo;application des lois conjointes. Dans le cadre th\u00e9orique de la <strong>th\u00e9orie de la r\u00e9ponse \u00e0 l&rsquo;item (TRI)<\/strong>, l&rsquo;ambition centrale est de relier un ensemble de variables manifestes observables (les r\u00e9ponses binaires ou polytomiques <span class=\"arabpsy-math-inline\" data-latex=\"Y_1, dots, Y_m\">$Y_1, dots, Y_m$<\/span> donn\u00e9es \u00e0 un questionnaire de personnalit\u00e9 ou d&rsquo;\u00e9valuation intellectuelle) \u00e0 une ou plusieurs variables continues inobservables appel\u00e9es traits latents, not\u00e9es <span class=\"arabpsy-math-inline\" data-latex=\"Theta\">$Theta$<\/span>.<\/p>\n<p>Le mod\u00e8le math\u00e9matique prend la forme d&rsquo;une distribution de probabilit\u00e9 conjointe param\u00e9tr\u00e9e reliant l&rsquo;ensemble des r\u00e9ponses observables et le continuum psychologique latent du sujet :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f(\\mathbf{y}, \\theta) = f_\\Theta(\\theta) \\cdot P(\\mathbf{Y} = \\mathbf{y} mid \\Theta = \\theta)\">$$f(\\mathbf{y}, \\theta) = f_\\Theta(\\theta) \\cdot P(\\mathbf{Y} = \\mathbf{y} mid \\Theta = \\theta)$$<\/span><\/p>\n<p>Pour rendre calculable ce produit multidimensionnel, la psychom\u00e9trie introduit l&rsquo;hypoth\u00e8se centrale d&rsquo;<strong>ind\u00e9pendance locale<\/strong> (ou principe d&rsquo;ind\u00e9pendance conditionnelle). Cette hypoth\u00e8se postule que pour un niveau de trait latent <span class=\"arabpsy-math-inline\" data-latex=\"\\theta\">$\\theta$<\/span> parfaitement contr\u00f4l\u00e9 et fix\u00e9, les r\u00e9ponses apport\u00e9es aux diff\u00e9rents items du test deviennent stochastiquement ind\u00e9pendantes les unes des autres :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(\\mathbf{Y} = \\mathbf{y} mid \\Theta = \\theta) = \\prod_{j=1}^m P(Y_j = y_j mid \\Theta = \\theta)\">$$P(\\mathbf{Y} = \\mathbf{y} mid \\Theta = \\theta) = \\prod_{j=1}^m P(Y_j = y_j mid \\Theta = \\theta)$$<\/span><\/p>\n<p>Chaque terme individuel <span class=\"arabpsy-math-inline\" data-latex=\"P(Y_j = 1 mid \\Theta = \\theta)\">$P(Y_j = 1 mid \\Theta = \\theta)$<\/span> est mod\u00e9lis\u00e9 par une fonction logistique (courbe caract\u00e9ristique de l&rsquo;item), param\u00e9tr\u00e9e par la difficult\u00e9 de la question et son pouvoir de discrimination. La distribution conjointe globale permet ainsi, par le truchement d&rsquo;algorithmes d&rsquo;estimation de type marginal maximum de vraisemblance (MMLE) utilisant l&rsquo;algorithme EM, d&rsquo;inf\u00e9rer simultan\u00e9ment les caract\u00e9ristiques psychom\u00e9triques des items et les coordonn\u00e9es factorielles des individus sur l&rsquo;espace latent.<\/p>\n<h3>11.2 Distribution conjointe du temps de r\u00e9ponse et de la pr\u00e9cision<\/h3>\n<p>En psychologie cognitive exp\u00e9rimentale, les th\u00e9ories traditionnelles \u00e9valuaient s\u00e9par\u00e9ment le taux d&rsquo;exactitude des r\u00e9ponses (pourcentage de r\u00e9ussite \u00e0 une t\u00e2che motrice ou perceptive) et la latence moyenne de r\u00e9ponse (chronom\u00e9trie mentale univari\u00e9e). Or, ce fractionnement m\u00e9thodologique ignore d\u00e9lib\u00e9r\u00e9ment le ph\u00e9nom\u00e8ne cognitif sous-jacent universel : le <strong>compromis vitesse-pr\u00e9cision<\/strong> (speed-accuracy tradeoff). Un sujet peut artificiellement maximiser son exactitude en adoptant une strat\u00e9gie d&rsquo;accumulation d&rsquo;information lente et pr\u00e9cautionneuse, ou inversement \u00e9courter ses latences au prix d&rsquo;une prolif\u00e9ration d&rsquo;erreurs stochastiques.<\/p>\n<p>Pour d\u00e9passer cette scission r\u00e9ductrice, la mod\u00e9lisation contemporaine s&rsquo;appuie sur des mod\u00e8les de processus stochastiques conjoints continus-discrets, au premier rang desquels figure le <strong>mod\u00e8le de diffusion de drift<\/strong> (Drift Diffusion Model, DDM) formul\u00e9 par Roger Ratcliff. Ce cadre math\u00e9matique d\u00e9crit le processus de d\u00e9lib\u00e9ration d\u00e9cisionnelle binaire comme l&rsquo;\u00e9volution temporelle d&rsquo;un mouvement brownien avec d\u00e9rive (d\u00e9rive de drift <span class=\"arabpsy-math-inline\" data-latex=\"v\">$v$<\/span>, variabilit\u00e9 stochastique <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>), progressant entre deux barri\u00e8res absorbantes distantes de <span class=\"arabpsy-math-inline\" data-latex=\"a\">$a$<\/span> (la fronti\u00e8re sup\u00e9rieure mat\u00e9rialisant la d\u00e9cision correcte, la fronti\u00e8re inf\u00e9rieure codant l&rsquo;erreur).<\/p>\n<p>Le point d&rsquo;arriv\u00e9e analytique du mod\u00e8le est une fonction de densit\u00e9 de probabilit\u00e9 conjointe bivari\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"f(t, r)\">$f(t, r)$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"T\">$T$<\/span> est une variable continue mesurant le temps de r\u00e9action physique et <span class=\"arabpsy-math-inline\" data-latex=\"R in {0, 1}\">$R in {0, 1}$<\/span> est une variable discr\u00e8te indiquant la pr\u00e9cision du choix. Cette densit\u00e9 conjointe calcule avec exactitude le temps de premier passage stochastique \u00e0 travers les fronti\u00e8res d&rsquo;absorption du mod\u00e8le :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f(t, R=1 mid v, a) = \\frac{\\pi}{a^2} e^{-v a - \\frac{v^2 t}{2}} \\sum_{k=1}^\\infty k \\sin(k \\pi) \\exp\\left( -\\frac{k^2 \\pi^2 t}{2 a^2} \\right)\">$$f(t, R=1 mid v, a) = \\frac{\\pi}{a^2} e^{-v a &#8211; \\frac{v^2 t}{2}} \\sum_{k=1}^\\infty k \\sin(k \\pi) \\exp\\left( -\\frac{k^2 \\pi^2 t}{2 a^2} \\right)$$<\/span><\/p>\n<p>Gr\u00e2ce \u00e0 cette loi conjointe, les chercheurs peuvent d\u00e9composer les performances individuelles en composantes neuro-cognitives invariantes : le taux d&rsquo;accumulation de preuve sensorielle (<span class=\"arabpsy-math-inline\" data-latex=\"v\">$v$<\/span>), le degr\u00e9 de conservatisme strat\u00e9gique (<span class=\"arabpsy-math-inline\" data-latex=\"a\">$a$<\/span>) et le temps d&rsquo;encodage perceptif et moteur non d\u00e9cisionnel (<span class=\"arabpsy-math-inline\" data-latex=\"T_{er}\">$T_{er}$<\/span>). L&rsquo;ajustement conjoint sur l&rsquo;ensemble de la surface bidimensionnelle temps-pr\u00e9cision surpasse infiniment les r\u00e9gressions univari\u00e9es conduites sur les seules moyennes arithm\u00e9tiques empiriques.<\/p>\n<h3>11.3 Validit\u00e9 de construit et structures de corr\u00e9lation inter-\u00e9chelles<\/h3>\n<p>L&rsquo;\u00e9valuation des instruments diagnostiques en neuropsychologie clinique et en psychopathologie quantitative (cadre HiTOP \u2014 <em>Hierarchical Taxonomy of Psychopathology<\/em>) fait un usage intensif de l&rsquo;analyse conjointe multidimensionnelle pour asseoir la <strong>validit\u00e9 de construit<\/strong>. Lorsqu&rsquo;un clinicien administre une batterie multi-\u00e9chelles mesurant simultan\u00e9ment le n\u00e9vrosisme, l&rsquo;affectivit\u00e9 n\u00e9gative et la symptomatologie obsessionnelle-compulsive, l&rsquo;enjeu diagnostique ne se limite pas aux scores isol\u00e9s.<\/p>\n<p>Historiquement formalis\u00e9e par Donald Campbell et Donald Fiske \u00e0 travers la matrice multi-traits multi-m\u00e9thodes (MTMM), cette d\u00e9marche examine l&rsquo;architecture de la loi conjointe sous l&rsquo;angle de la validit\u00e9 convergente (les distributions conjointes de construits homologues mesur\u00e9s par des instruments distincts doivent pr\u00e9senter une forte densit\u00e9 sur la diagonale principale) et de la validit\u00e9 discriminante (les distributions conjointes liant des construits th\u00e9oriquement orthogonaux doivent exhiber une circularit\u00e9 quasi parfaite, t\u00e9moignant d&rsquo;une d\u00e9corr\u00e9lation robuste).<\/p>\n<p>De surcro\u00eet, l&rsquo;inspection visuelle et g\u00e9om\u00e9trique de la distribution conjointe permet d&rsquo;isoler des profils nosographiques atypiques situ\u00e9s aux confins des queues de distribution multidimensionnelles. Un patient manifestant simultan\u00e9ment un score extr\u00eame de d\u00e9sinhibition comportementale et un score exceptionnellement bas de r\u00e9activit\u00e9 anxieuse occupera une position excentr\u00e9e singuli\u00e8re au sein du volume probabiliste conjoint. Cette conjonction pr\u00e9cise, invisible si l&rsquo;on examine chaque trait s\u00e9par\u00e9ment, signe une configuration psychopathologique sp\u00e9cifique (comme la psychopathie primaire) dont la prise en charge clinique diff\u00e8re radicalement de celle de sujets occupant d&rsquo;autres quadrants de l&rsquo;espace multivari\u00e9.<\/p>\n<h2>12. Estimation empirique, inf\u00e9rence et \u00e9cueils computationnels<\/h2>\n<h3>12.1 M\u00e9thodes d&rsquo;estimation param\u00e9trique et non param\u00e9trique<\/h3>\n<p>Dans la pratique de la recherche empirique, la v\u00e9ritable loi conjointe sous-jacente d&rsquo;une population demeure inconnue de l&rsquo;observateur. Elle doit \u00eatre inf\u00e9r\u00e9e \u00e0 partir d&rsquo;un \u00e9chantillon fini de r\u00e9alisations vectorielles ind\u00e9pendantes et identiquement distribu\u00e9es (i.i.d.) <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{D} = {\\mathbf{x}_1, \\mathbf{x}_2, dots, \\mathbf{x}_n}\">$\\mathcal{D} = {\\mathbf{x}_1, \\mathbf{x}_2, dots, \\mathbf{x}_n}$<\/span>, o\u00f9 chaque observation est un point de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^p\">$\\mathbb{R}^p$<\/span>. Deux paradigmes th\u00e9oriques s&rsquo;opposent et se compl\u00e8tent pour concr\u00e9tiser cette estimation : l&rsquo;approche param\u00e9trique et l&rsquo;approche non param\u00e9trique.<\/p>\n<p>L&rsquo;<strong>estimation param\u00e9trique<\/strong> postule a priori une famille math\u00e9matique explicite pour la loi conjointe, index\u00e9e par un vecteur de param\u00e8tres finis <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\theta}\">$boldsymbol{\\theta}$<\/span> (par exemple, une loi normale bivari\u00e9e r\u00e9gie par <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\theta} = (\\mu_X, \\mu_Y, \\sigma_X, \\sigma_Y, \\rho)\">$boldsymbol{\\theta} = (\\mu_X, \\mu_Y, \\sigma_X, \\sigma_Y, \\rho)$<\/span>). L&rsquo;estimation s&rsquo;op\u00e8re canoniquement par la m\u00e9thode du <strong>maximum de vraisemblance<\/strong> (MLE), consistant \u00e0 d\u00e9terminer la valeur optimale des param\u00e8tres qui maximise la vraisemblance conjointe globale de l&rsquo;\u00e9chantillon collect\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{boldsymbol{\\theta}}_{\\mathrm{MLE}} = arg\\max_{boldsymbol{\\theta}} \\sum_{i=1}^n \\log f(\\mathbf{x}_i mid boldsymbol{\\theta})\">$$\\hat{boldsymbol{\\theta}}_{\\mathrm{MLE}} = arg\\max_{boldsymbol{\\theta}} \\sum_{i=1}^n \\log f(\\mathbf{x}_i mid boldsymbol{\\theta})$$<\/span><\/p>\n<p>Lorsque des donn\u00e9es sont incompl\u00e8tes ou que des variables latentes inobserv\u00e9es s&rsquo;immiscent dans le syst\u00e8me, le recours \u00e0 l&rsquo;algorithme d&rsquo;esp\u00e9rance-maximisation (EM) devient incontournable pour optimiser it\u00e9rativement l&rsquo;esp\u00e9rance de la log-vraisemblance conjointe conditionnelle.<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;<strong>estimation non param\u00e9trique<\/strong> s&rsquo;affranchit de toute conjecture structurelle rigide quant \u00e0 la g\u00e9om\u00e9trie de la distribution sous-jacente. L&rsquo;approche phare est l&rsquo;<strong>estimation de la densit\u00e9 conjointe par noyau bivari\u00e9<\/strong> (Bivariate Kernel Density Estimation, KDE). Pour un \u00e9chantillon bivari\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"(X_i, Y_i)_{i=1,dots,n}\">$(X_i, Y_i)_{i=1,dots,n}$<\/span>, l&rsquo;estimateur de densit\u00e9 conjointe s&rsquo;\u00e9crit formellement :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{f}_{\\mathbf{H}}(x, y) = \\frac{1}{n |\\mathbf{H}|^{1\/2}} \\sum_{i=1}^n K\\left( \\mathbf{H}^{-1\/2} \\begin{\\pmatrix} x - X_i  y - Y_i \\end{\\pmatrix} \\right)\">$$\\hat{f}_{\\mathbf{H}}(x, y) = \\frac{1}{n |\\mathbf{H}|^{1\/2}} \\sum_{i=1}^n K\\left( \\mathbf{H}^{-1\/2} \\begin{\\pmatrix} x &#8211; X_i  y &#8211; Y_i \\end{\\pmatrix} \\right)$$<\/span><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span> repr\u00e9sente une fonction de noyau bivari\u00e9e multivari\u00e9e sym\u00e9trique (couramment un noyau gaussien standard) et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{H}\">$\\mathbf{H}$<\/span> d\u00e9signe la <strong>matrice de lissage<\/strong> (bandwidth matrix) sym\u00e9trique d\u00e9finie positive de dimension <span class=\"arabpsy-math-inline\" data-latex=\"2 \\times 2\">$2 \\times 2$<\/span>. Le choix de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{H}\">$\\mathbf{H}$<\/span> est d&rsquo;une criticit\u00e9 m\u00e9thodologique absolue : une matrice sous-dimensionn\u00e9e induit un \u00ab sous-lissage \u00bb g\u00e9n\u00e9rant une multitude d&rsquo;asp\u00e9rit\u00e9s et de modes parasites art\u00e9factuels, tandis qu&rsquo;une matrice surdimensionn\u00e9e provoque un \u00ab sur-lissage \u00bb diluant totalement la topologie conjointe et masquant les corr\u00e9lations et sous-populations r\u00e9elles.<\/p>\n<h3>12.2 G\u00e9n\u00e9ration d&rsquo;\u00e9chantillons et simulation de Monte-Carlo<\/h3>\n<p>L&rsquo;exploration des distributions conjointes complexes, notamment lorsqu&rsquo;elles ne disposent pas de formes int\u00e9grales analytiquement solubles, requiert des outils de simulation stochastique computationnelle : les m\u00e9thodes de <strong>Monte-Carlo<\/strong>. Pour simuler informatiquement un vecteur al\u00e9atoire gouvern\u00e9 par une loi conjointe continue, la m\u00e9thode universelle en basse dimension s&rsquo;appuie sur la r\u00e8gle de d\u00e9composition conditionnelle en cha\u00eene combin\u00e9e \u00e0 la <strong>m\u00e9thode de la transformation inverse<\/strong> univari\u00e9e :<\/p>\n<ol>\n<li>On \u00e9chantillonne un premier nombre pseudo-al\u00e9atoire uniforme <span class=\"arabpsy-math-inline\" data-latex=\"u_1 \\sim \\mathcal{U}[0, 1]\">$u_1 \\sim \\mathcal{U}[0, 1]$<\/span> et on d\u00e9termine <span class=\"arabpsy-math-inline\" data-latex=\"x = F_X^{-1}(u_1)\">$x = F_X^{-1}(u_1)$<\/span> via l&rsquo;inverse de la fonction de r\u00e9partition marginale cumulative de <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span>.<\/li>\n<li>On \u00e9chantillonne un second nombre uniforme <span class=\"arabpsy-math-inline\" data-latex=\"u_2 \\sim \\mathcal{U}[0, 1]\">$u_2 \\sim \\mathcal{U}[0, 1]$<\/span> et on r\u00e9sout la valeur <span class=\"arabpsy-math-inline\" data-latex=\"y = F_{Y mid X}^{-1}(u_2 mid x)\">$y = F_{Y mid X}^{-1}(u_2 mid x)$<\/span> au moyen de la fonction de r\u00e9partition conditionnelle calcul\u00e9e au point <span class=\"arabpsy-math-inline\" data-latex=\"x\">$x$<\/span> pr\u00e9alablement fix\u00e9.<\/li>\n<\/ol>\n<p>Dans des espaces de dimension arbitraire ou lorsque la constante de normalisation de la loi conjointe s&rsquo;av\u00e8re intra\u00e7able (comme dans l&rsquo;inf\u00e9rence a posteriori de mod\u00e8les bay\u00e9siens complexes), on d\u00e9ploie les algorithmes de <strong>cha\u00eenes de Markov Monte-Carlo (MCMC)<\/strong>. L&rsquo;<strong>\u00e9chantillonneur de Gibbs<\/strong> est sp\u00e9cifiquement adapt\u00e9 \u00e0 l&rsquo;architecture des lois conjointes : au lieu de tenter de tirer directement le vecteur multidimensionnel <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X} = (X_1, dots, X_p)\">$\\mathbf{X} = (X_1, dots, X_p)$<\/span> en un bloc unique, l&rsquo;algorithme tire cycliquement chaque variable scalaire selon sa distribution conditionnelle univari\u00e9e compl\u00e8te sachant l&rsquo;ensemble de toutes les autres variables actualis\u00e9es :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"X_j^{(t+1)} \\sim f\\left(X_j mid X_1^{(t+1)}, dots, X_{j-1}^{(t+1)}, X_{j+1}^{(t)}, dots, X_p^{(t)}\\right)\">$$X_j^{(t+1)} \\sim f\\left(X_j mid X_1^{(t+1)}, dots, X_{j-1}^{(t+1)}, X_{j+1}^{(t)}, dots, X_p^{(t)}\\right)$$<\/span><\/p>\n<p>Sous r\u00e9serve du respect des conditions de cha\u00eene ergodique et de r\u00e9versibilit\u00e9 microscopique, la distribution stationnaire de cette cha\u00eene de Markov converge exactement vers la loi conjointe cible <span class=\"arabpsy-math-inline\" data-latex=\"f(\\mathbf{x})\">$f(\\mathbf{x})$<\/span>, permettant d&rsquo;approximer num\u00e9riquement toute int\u00e9grale ou moment complexe avec une pr\u00e9cision r\u00e9gie par la loi des grands nombres.<\/p>\n<h3>12.3 \u00c9cueils fr\u00e9quents et erreurs d&rsquo;interpr\u00e9tation m\u00e9thodologique<\/h3>\n<p>L&rsquo;utilisation pratique des distributions conjointes en recherche appliqu\u00e9e s&rsquo;accompagne d&rsquo;une s\u00e9rie de pi\u00e8ges \u00e9pist\u00e9mologiques et computationnels qui peuvent compromettre la validit\u00e9 des conclusions scientifiques :<\/p>\n<p><strong>1. L&rsquo;assimilation erron\u00e9e de la d\u00e9pendance conjointe \u00e0 la causalit\u00e9 :<\/strong> Une distribution conjointe <span class=\"arabpsy-math-inline\" data-latex=\"f(x, y)\">$f(x, y)$<\/span> est, par essence, une construction stochastique purement sym\u00e9trique. L&rsquo;observation d&rsquo;une concentration conjointe \u00e9lev\u00e9e, d&rsquo;une covariance forte ou d&rsquo;un rapport de cotes d\u00e9s\u00e9quilibr\u00e9 n&rsquo;autorise en aucun cas \u00e0 d\u00e9duire qu&rsquo;une variable exerce une action causale sur la seconde. En l&rsquo;absence d&rsquo;assignation randomis\u00e9e ou d&rsquo;intervention au sens du calcul des interventions de Judea Pearl (l&rsquo;op\u00e9rateur <span class=\"arabpsy-math-inline\" data-latex=\"do(X = x)\">$do(X = x)$<\/span>), la covariation mesur\u00e9e au sein d&rsquo;une distribution conjointe peut r\u00e9sulter int\u00e9gralement de l&rsquo;action non mesur\u00e9e d&rsquo;un facteur de confusion latent <span class=\"arabpsy-math-inline\" data-latex=\"Z\">$Z$<\/span> gouvernant simultan\u00e9ment les deux variables observ\u00e9es.<\/p>\n<p><strong>2. Le paradoxe de Simpson :<\/strong> Cet \u00e9cueil spectaculaire met en sc\u00e8ne l&rsquo;inversion pure et simple de l&rsquo;orientation d&rsquo;une d\u00e9pendance conjointe lors de l&rsquo;agr\u00e9gation ou de la marginalisation des donn\u00e9es. Deux variables <span class=\"arabpsy-math-inline\" data-latex=\"X\">$X$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span> peuvent pr\u00e9senter une corr\u00e9lation marginale positive lorsqu&rsquo;on \u00e9value leur distribution conjointe globale, tout en manifestant une corr\u00e9lation conditionnelle rigoureusement n\u00e9gative au sein de <em>chaque<\/em> sous-strate d\u00e9finie par une troisi\u00e8me variable polytomique <span class=\"arabpsy-math-inline\" data-latex=\"Z\">$Z$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathrm{Cov}(X, Y mid Z = z)  0\">$$\\mathrm{Cov}(X, Y mid Z = z)  0$$<\/span><\/p>\n<p>Ce paradoxe rappelle avec force qu&rsquo;une loi conjointe bivari\u00e9e apparente peut \u00eatre une projection trompeuse d\u00e8s lors qu&rsquo;elle marginalise \u00e0 son insu une variable de stratification d\u00e9terminante.<\/p>\n<p><strong>3. La d\u00e9formation stochastique induite par les donn\u00e9es manquantes :<\/strong> Dans les protocoles comportementaux, les observations vectorielles sont fr\u00e9quemment incompl\u00e8tes. Selon la taxonomie fondamentale \u00e9tablie par Donald Rubin, les m\u00e9canismes de perte se divisent en donn\u00e9es manquantes compl\u00e8tement de mani\u00e8re al\u00e9atoire (MCAR), manquantes de mani\u00e8re al\u00e9atoire (MAR) et manquantes de mani\u00e8re non al\u00e9atoire (MNAR). Traiter un profil conjoint lacunaire par simple suppression des listes (listwise deletion) lorsque le m\u00e9canisme sous-jacent est MAR ou MNAR d\u00e9truit la repr\u00e9sentativit\u00e9 du sous-espace d&rsquo;\u00e9chantillonnage et introduit un biais syst\u00e9matique majeur dans l&rsquo;estimation de la matrice de covariance et des param\u00e8tres conjoints.<\/p>\n<p><strong>4. Le biais de troncature et de censure conjointe :<\/strong> Lorsque la collecte des donn\u00e9es impose des seuils de d\u00e9tection physiques ou m\u00e9thodologiques sur une variable (par exemple, un temps de r\u00e9action plafonn\u00e9 \u00e0 1500 millisecondes ou une \u00e9chelle psychom\u00e9trique satur\u00e9e par effet de plafond), la distribution conjointe observ\u00e9e subit une distorsion g\u00e9om\u00e9trique s\u00e9v\u00e8re. Cette coupure artificielle rabote les ellipso\u00efdes de confiance et contractedrastiquement les variances et covariances apparentes (effet d&rsquo;att\u00e9nuation de restriction de variance). L&rsquo;omission de correctifs statistiques appropri\u00e9s (comme les mod\u00e8les de Tobit bivari\u00e9s) aboutit \u00e0 une sous-estimation syst\u00e9matique de l&rsquo;intensit\u00e9 de la d\u00e9pendance conjointe r\u00e9elle liant les attributs examin\u00e9s.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul style=\"list-style-type: none;padding-left: 0\">\n<li style=\"margin-bottom: 1rem\">Ash, R. B. (2000). <em>Probability and measure theory<\/em> (2nd ed.). Academic Press.<\/li>\n<li style=\"margin-bottom: 1rem\">Casella, G., &amp; Berger, R. L. (2002). <em>Statistical inference<\/em> (2nd ed.). Duxbury Thomson Learning.<\/li>\n<li style=\"margin-bottom: 1rem\">Gelman, A., Carlin, J. B., Stern, H. S., Dunson, D. B., Vehtari, A., &amp; Rubin, D. B. (2013). <em>Bayesian data analysis<\/em> (3rd ed.). Chapman and Hall\/CRC. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1201\/b16018\">https:\/\/doi.org\/10.1201\/b16018<\/a><\/li>\n<li style=\"margin-bottom: 1rem\">Joe, H. (2014). <em>Dependence modeling with copulas<\/em>. Chapman and Hall\/CRC. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1201\/b17116\">https:\/\/doi.org\/10.1201\/b17116<\/a><\/li>\n<li style=\"margin-bottom: 1rem\">Kolmogorov, A. N. (1956). <em>Foundations of the theory of probability<\/em> (2nd ed., N. Morrison, Trans.). Chelsea Publishing Co.<\/li>\n<li style=\"margin-bottom: 1rem\">Nelsen, R. B. (2006). <em>An introduction to copulas<\/em> (2nd ed.). Springer. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/0-387-28678-0\">https:\/\/doi.org\/10.1007\/0-387-28678-0<\/a><\/li>\n<li style=\"margin-bottom: 1rem\">Pearl, J. (2009). <em>Causality: Models, reasoning, and inference<\/em> (2nd ed.). Cambridge University Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1017\/CBO9780511803161\">https:\/\/doi.org\/10.1017\/CBO9780511803161<\/a><\/li>\n<li style=\"margin-bottom: 1rem\">Ratcliff, R., &amp; McKoon, G. (2008). The diffusion decision model: Theory and data for two-choice decision tasks. <em>Neural Computation<\/em>, 20(4), 873\u2013922. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1162\/neco.2008.12-06-420\">https:\/\/doi.org\/10.1162\/neco.2008.12-06-420<\/a><\/li>\n<li style=\"margin-bottom: 1rem\">Sklar, A. (1959). Fonctions de r\u00e9partition \u00e0 n dimensions et leurs marges. <em>Publications de l&rsquo;Institut de Statistique de l&rsquo;Universit\u00e9 de Paris<\/em>, 8, 229\u2013231.<\/li>\n<li style=\"margin-bottom: 1rem\">Wasserman, L. (2004). <em>All of statistics: A concise course in statistical inference<\/em>. Springer. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-0-387-21736-9\">https:\/\/doi.org\/10.1007\/978-0-387-21736-9<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>D\u00e9couvrez la distribution de probabilit\u00e9 conjointe : d\u00e9finition rigoureuse, formalisme math\u00e9matique, propri\u00e9t\u00e9s fondamentales et applications m\u00e9thodologiques.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2],"tags":[3791,7722,4452,7720,7721],"class_list":["post-5907","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-psychometrie","tag-analyse-multivariee","tag-densite-conjointe","tag-modelisation-stochastique","tag-probabilite-conjointe","tag-statistiques-bivariees"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/5907","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=5907"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/5907\/revisions"}],"predecessor-version":[{"id":5911,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/5907\/revisions\/5911"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=5907"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=5907"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=5907"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}