{"id":9892,"date":"2026-10-14T18:02:35","date_gmt":"2026-10-14T18:02:35","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/dissimilarite-de-bray-curtis-definition-et-exemples\/"},"modified":"2026-10-14T18:02:35","modified_gmt":"2026-10-14T18:02:35","slug":"dissimilarite-de-bray-curtis-definition-et-exemples","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=dissimilarite-de-bray-curtis-definition-et-exemples","title":{"rendered":"Dissimilarit\u00e9 de Bray-Curtis : D\u00e9finition et Exemples"},"content":{"rendered":"<p>L&rsquo;analyse multivari\u00e9e contemporaine repose sur une pr\u00e9misse fondamentale : la capacit\u00e9 de quantifier avec exactitude et pertinence th\u00e9orique la divergence ou la ressemblance entre des entit\u00e9s complexes d\u00e9crites par de multiples variables. Qu&rsquo;il s&rsquo;agisse de comparer la structure taxonomique de deux \u00e9cosyst\u00e8mes forestiers, de mesurer l&rsquo;\u00e9volution spatiotemporelle d&rsquo;un microbiote intestinal sous traitement th\u00e9rapeutique, ou de caract\u00e9riser la divergence clinique entre les r\u00e9pertoires symptomatiques de deux patients en psychiatrie computationnelle, le choix de la m\u00e9trique d&rsquo;association constitue le pivot m\u00e9thodologique d\u00e9terminant l&rsquo;ensemble des inf\u00e9rences statistiques ult\u00e9rieures. Parmi la vaste pl\u00e9thore d&rsquo;indices \u00e9labor\u00e9s au cours du XXe si\u00e8cle, la dissimilarit\u00e9 de Bray-Curtis occupe une place pr\u00e9pond\u00e9rante, s&rsquo;\u00e9rigeant en v\u00e9ritable standard de facto dans les disciplines manipulant des matrices d&rsquo;abondances clairsem\u00e9es et h\u00e9t\u00e9rog\u00e8nes.<\/p>\n<p>Con\u00e7ue initialement pour r\u00e9pondre aux impasses m\u00e9thodologiques des distances g\u00e9om\u00e9triques classiques face aux inventaires floristiques, cette mesure a rapidement transcend\u00e9 son cadre disciplinaire d&rsquo;origine. Les \u00e9cologues J. Roger Bray et John Thomas Curtis ont introduit une formulation alg\u00e9brique \u00e0 la fois intuitive et d&rsquo;une redoutable efficacit\u00e9, capable d&rsquo;int\u00e9grer l&rsquo;intensit\u00e9 quantitative des descripteurs tout en neutralisant le biais structurel induit par les absences conjointes. Alors que la distance euclidienne tend \u00e0 distordre la r\u00e9alit\u00e9 \u00e9cologique en traitant l&rsquo;absence simultan\u00e9e d&rsquo;un taxon comme une preuve de ressemblance, la m\u00e9trique de Bray-Curtis concentre son pouvoir discriminant sur l&rsquo;information biologique positivement partag\u00e9e et sur les diff\u00e9rentiels r\u00e9els d&rsquo;abondance.<\/p>\n<p>La robustesse conceptuelle de cette mesure explique son adoption g\u00e9n\u00e9ralis\u00e9e bien au-del\u00e0 de la phyto\u00e9cologie, trouvant des applications d\u00e9cisives en bio-informatique, en m\u00e9tag\u00e9nomique environnementale et, plus r\u00e9cemment, dans la psychom\u00e9trie des r\u00e9seaux et l&rsquo;analyse dimensionnelle du comportement. Comprendre la dissimilarit\u00e9 de Bray-Curtis implique d&rsquo;explorer ses fondements historiques, de diss\u00e9quer son appareil math\u00e9matique, d&rsquo;analyser ses propri\u00e9t\u00e9s m\u00e9triques particuli\u00e8res \u2014 notamment son statut formel de semi-m\u00e9trique \u2014 et de ma\u00eetriser son impl\u00e9mentation algorithmique au sein des architectures modernes d&rsquo;ordination multivari\u00e9e. Cette analyse exhaustive propose un parcours th\u00e9orique, m\u00e9thodologique et appliqu\u00e9 \u00e0 travers l&rsquo;une des m\u00e9triques les plus influentes de la statistique multivari\u00e9e appliqu\u00e9e.<\/p>\n<h2>1. Origines historiques et fondements th\u00e9oriques de la mesure de Bray-Curtis<\/h2>\n<h3>1.1 Le contexte d&rsquo;\u00e9mergence en \u00e9cologie quantitative<\/h3>\n<p>L&rsquo;\u00e9cologie v\u00e9g\u00e9tale des ann\u00e9es 1950 traversait une profonde r\u00e9volution conceptuelle et m\u00e9thodologique. Sous l&rsquo;impulsion de chercheurs visionnaires \u00e0 l&rsquo;Universit\u00e9 du Wisconsin \u00e0 Madison, la vision holistique et organismale des communaut\u00e9s v\u00e9g\u00e9tales, longtemps d\u00e9fendue par Frederic Clements, commen\u00e7ait \u00e0 c\u00e9der le pas face au concept individualiste de l&rsquo;association v\u00e9g\u00e9tale formalis\u00e9 par Henry Gleason. Dans ce paradigme \u00e9mergent, les esp\u00e8ces v\u00e9g\u00e9tales ne sont plus per\u00e7ues comme des entit\u00e9s indissociables formant des supra-organismes discrets, mais comme des populations ind\u00e9pendantes se distribuant de mani\u00e8re continue et individualis\u00e9e le long de gradients \u00e9cologiques et environnementaux complexes. C&rsquo;est dans ce bouillonnement intellectuel que <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1942268\">J. Roger Bray et John Thomas Curtis<\/a> publi\u00e8rent en 1957 leur \u00e9tude fondatrice consacr\u00e9e aux communaut\u00e9s foresti\u00e8res des hautes terres du sud du Wisconsin, intitul\u00e9e <em>An Ordination of the Upland Forest Communities of Southern Wisconsin<\/em>.<\/p>\n<p>\u00c0 cette \u00e9poque, les \u00e9cologues quantitatifs \u00e9taient confront\u00e9s \u00e0 une barri\u00e8re m\u00e9thodologique majeure : l&rsquo;absence d&rsquo;outils analytiques adapt\u00e9s pour quantifier de fa\u00e7on r\u00e9aliste la ressemblance entre des stations d&rsquo;\u00e9chantillonnage \u00e9cologiquement diversifi\u00e9es. Les coefficients d&rsquo;association disponibles \u00e9taient pour la plupart restreints \u00e0 des formulations qualitatives binaires d&rsquo;occurrence, telles que l&rsquo;indice de Jaccard ou le coefficient de S\u00f8rensen. Bien que ces mesures permettent d&rsquo;\u00e9valuer le simple partage de pr\u00e9sence ou d&rsquo;absence d&rsquo;esp\u00e8ces, elles s&rsquo;av\u00e9raient totalement impuissantes \u00e0 capturer les variations subtiles mais d\u00e9terminantes de l&rsquo;abondance relative, de la dominance d\u00e9mographique et de la biomasse v\u00e9g\u00e9tale. Deux peuplements forestiers partageant les m\u00eames essences arbor\u00e9es mais dans des proportions radicalement inverses \u2014 l&rsquo;un domin\u00e9 par le ch\u00eane blanc avec de rares \u00e9rables, l&rsquo;autre \u00e9cras\u00e9 par la canop\u00e9e de l&rsquo;\u00e9rable avec de sporadiques ch\u00eanes r\u00e9siduels \u2014 \u00e9taient jug\u00e9s math\u00e9matiquement indiscernables par les m\u00e9triques qualitatives de l&rsquo;\u00e9poque.<\/p>\n<p>Conscients de cette lacune critique, Bray et Curtis comprirent qu&rsquo;une v\u00e9ritable analyse de gradient, appel\u00e9e ordination polaire dans leurs travaux fondateurs, n\u00e9cessitait une fonction de calcul capable d&rsquo;int\u00e9grer harmonieusement les comptages r\u00e9els, les densit\u00e9s relatives ou les fr\u00e9quences de recouvrement des esp\u00e8ces. En adaptant les principes de comparaison proportionnelle aux abondances absolues ou standardis\u00e9es, les deux chercheurs \u00e9labor\u00e8rent une mesure qui allait devenir la pierre angulaire de l&rsquo;ordination floristique, ouvrant la voie \u00e0 la transition d\u00e9finitive de l&rsquo;\u00e9cologie descriptive vers une \u00e9cologie num\u00e9rique rigoureuse, math\u00e9matis\u00e9e et empiriquement v\u00e9rifiable.<\/p>\n<h3>1.2 Philosophie de la mesure des communaut\u00e9s et des profils<\/h3>\n<p>La philosophie sous-jacente \u00e0 la dissimilarit\u00e9 de Bray-Curtis repose sur un postulat \u00e9cologique fondamental : deux communaut\u00e9s vivantes se ressemblent d&rsquo;autant plus qu&rsquo;elles partagent des volumes concrets d&rsquo;individus appartenant aux m\u00eames taxons, et non pas simplement des identit\u00e9s taxonomiques abstraites. Dans la nature, la structure d&rsquo;un assemblage \u00e9cologique est d\u00e9finie par le flux d&rsquo;\u00e9nergie et la mati\u00e8re que chaque population mobilise au sein du biotope. Ainsi, la v\u00e9ritable ressemblance entre deux \u00e9cosyst\u00e8mes d\u00e9coule du partage mat\u00e9riel de leurs abondances r\u00e9elles. Cette conception s&rsquo;oppose radicalement aux approches purement g\u00e9om\u00e9triques qui consid\u00e8rent les descripteurs comme de simples dimensions orthogonales au sein d&rsquo;un espace vectoriel abstrait.<\/p>\n<p>Sur le plan \u00e9pist\u00e9mologique, Bray et Curtis ont introduit une distinction cruciale entre la dissemblance d&rsquo;\u00e9chantillonnage \u00e9cologique et la distance spatiale euclidienne. Dans un rep\u00e8re physique tridimensionnel, la distance entre deux points est invariante par rapport aux propri\u00e9t\u00e9s intrins\u00e8ques de l&rsquo;espace. En revanche, dans l&rsquo;espace des profils \u00e9cologiques, la distance per\u00e7ue est conditionn\u00e9e par la biologie des organismes : un site o\u00f9 aucune esp\u00e8ce n&rsquo;est observ\u00e9e ne peut \u00eatre consid\u00e9r\u00e9 comme infiniment proche d&rsquo;un autre site vide sur la base d&rsquo;une proximit\u00e9 g\u00e9om\u00e9trique. La mesure de Bray-Curtis pose pour principe que la quantit\u00e9 d&rsquo;information partag\u00e9e r\u00e9side dans le minimum d&rsquo;abondance d&rsquo;un descripteur commun aux deux unit\u00e9s d&rsquo;\u00e9chantillonnage. Si une esp\u00e8ce est pr\u00e9sente avec une abondance de 50 individus sur le site A et de 20 individus sur le site B, le niveau d&rsquo;\u00e9change, d&rsquo;interaction partag\u00e9e ou de compatibilit\u00e9 \u00e9cologique pour cette esp\u00e8ce est strictement born\u00e9 par les 20 unit\u00e9s communes.<\/p>\n<p>Cette approche relationnelle a exerc\u00e9 une influence colossale sur le d\u00e9veloppement de l&rsquo;\u00e9cologie num\u00e9rique moderne, th\u00e9oris\u00e9e ult\u00e9rieurement par des figures tut\u00e9laires telles que <a rel=\"noopener noreferrer\" href=\"https:\/\/www.elsevier.com\/books\/numerical-ecology\/legendre\/978-0-444-53868-0\">Pierre Legendre et Louis Legendre<\/a>. Elle a permis de conceptualiser les assemblages biologiques non pas comme des vecteurs inertes plong\u00e9s dans une g\u00e9om\u00e9trie rigide, mais comme des profils dynamiques dont la proximit\u00e9 mutuelle doit refl\u00e9ter fid\u00e8lement le chevauchement fonctionnel et d\u00e9mographique r\u00e9el au sein des \u00e9cosyst\u00e8mes observ\u00e9s.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/bray_curtis.png\" alt=\"\" \/><\/figure>\n<h3>1.3 Du paradigme \u00e9cologique aux sciences comportementales<\/h3>\n<p>Bien que forg\u00e9e au c\u0153ur des for\u00eats feuillues nord-am\u00e9ricaines, la dissimilarit\u00e9 de Bray-Curtis poss\u00e8de des propri\u00e9t\u00e9s math\u00e9matiques intrins\u00e8ques qui la rendent universellement applicable \u00e0 tout syst\u00e8me complexe caract\u00e9ris\u00e9 par des r\u00e9pertoires h\u00e9t\u00e9rog\u00e8nes et des profils d&rsquo;intensit\u00e9. C&rsquo;est ainsi que la psychom\u00e9trie quantitative et les sciences comportementales se sont progressivement appropri\u00e9 cette m\u00e9trique au cours des derni\u00e8res d\u00e9cennies. En psychologie clinique et en psychopathologie quantitative, la description d&rsquo;un individu s&rsquo;op\u00e8re traditionnellement via l&rsquo;\u00e9valuation d&rsquo;une multiplicit\u00e9 de sympt\u00f4mes, dont la s\u00e9v\u00e9rit\u00e9 est scor\u00e9e sur des \u00e9chelles d&rsquo;intensit\u00e9 continues ou ordinales (par exemple les \u00e9chelles de d\u00e9pression, d&rsquo;anxi\u00e9t\u00e9 ou de d\u00e9sorganisation cognitive).<\/p>\n<p>Dans ce contexte, la transposition conceptuelle est directe et f\u00e9conde : un individu ou un patient peut \u00eatre mod\u00e9lis\u00e9 math\u00e9matiquement de la m\u00eame mani\u00e8re qu&rsquo;un site d&rsquo;\u00e9chantillonnage \u00e9cologique, tandis que les sympt\u00f4mes cliniques correspondent aux diff\u00e9rentes esp\u00e8ces v\u00e9g\u00e9tales ou animales composant le profil global. Deux patients consultant pour des troubles de l&rsquo;humeur ne pr\u00e9sentent presque jamais des constellations s\u00e9miologiques binaires ; ils exhibent des variations d&rsquo;amplitude consid\u00e9rables sur des dimensions telles que l&rsquo;insomnie, l&rsquo;anh\u00e9donie, l&rsquo;agitation motrice ou les id\u00e9ations suicidaires. L&rsquo;adoption de la m\u00e9trique de Bray-Curtis permet de formaliser la distance inter-individuelle en se fondant sur l&rsquo;intensit\u00e9 exacte des traits cliniques manifest\u00e9s, tout en r\u00e9solvant un probl\u00e8me psychom\u00e9trique majeur : le traitement des sympt\u00f4mes absents.<\/p>\n<p>La validit\u00e9 conceptuelle de la comparaison de r\u00e9pertoires comportementaux au moyen de Bray-Curtis s&rsquo;av\u00e8re particuli\u00e8rement pertinente dans l&rsquo;analyse des patrons d&rsquo;activation comportementale et l&rsquo;\u00e9thologie humaine. Lorsque des chercheurs enregistrent les dur\u00e9es quotidiennes allou\u00e9es par diff\u00e9rents individus \u00e0 une typologie \u00e9tendue d&rsquo;activit\u00e9s (travail, loisirs passifs, interactions sociales, pratiques sportives, sommeil), la similarit\u00e9 de mode de vie ne d\u00e9pend pas seulement de la nature des activit\u00e9s pratiqu\u00e9es, mais du volume temporel minimal partag\u00e9 au sein de chaque cat\u00e9gorie comportementale. La m\u00e9trique de Bray-Curtis capture cette r\u00e9alit\u00e9 empirique avec une pr\u00e9cision que les coefficients de corr\u00e9lation lin\u00e9aire ou les distances euclidiennes ne peuvent atteindre, consolidant ainsi son statut d&rsquo;outil transdisciplinaire incontournable.<\/p>\n<h2>2. D\u00e9finition math\u00e9matique et formulation formelle de la dissimilarit\u00e9<\/h2>\n<h3>2.1 \u00c9quation canonique et d\u00e9composition des termes<\/h3>\n<p>La formulation alg\u00e9brique canonique de la dissimilarit\u00e9 de Bray-Curtis entre deux unit\u00e9s d&rsquo;\u00e9chantillonnage, d\u00e9sign\u00e9es conventionnellement par les indices <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span>, s&rsquo;articule autour du rapport entre le volume d&rsquo;abondances minimales partag\u00e9es et la masse globale des observations recens\u00e9es sur les deux sites. Soit un ensemble de <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> descripteurs (esp\u00e8ces, g\u00e8nes, items psychologiques ou comportements) mesur\u00e9s sur deux unit\u00e9s distinctes. Le vecteur d&rsquo;observation associ\u00e9 \u00e0 l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> est not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_i = (x_{i1}, x_{i2}, dots, x_{ip})\">$\\mathbf{x}_i = (x_{i1}, x_{i2}, dots, x_{ip})$<\/span> et celui de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span> est not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_j = (x_{j1}, x_{j2}, dots, x_{jp})\">$\\mathbf{x}_j = (x_{j1}, x_{j2}, dots, x_{jp})$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"x_{ik}\">$x_{ik}$<\/span> repr\u00e9sente l&rsquo;abondance, la fr\u00e9quence ou l&rsquo;intensit\u00e9 du descripteur <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> au sein de l&rsquo;unit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>.<\/p>\n<p>L&rsquo;expression math\u00e9matique standardis\u00e9e de la dissimilarit\u00e9 de Bray-Curtis, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij}\">$BC_{ij}$<\/span>, s&rsquo;\u00e9nonce selon la relation canonique suivante :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij} = 1 - \\frac{2 C_{ij}}{S_i + S_j}\">$$BC_{ij} = 1 &#8211; \\frac{2 C_{ij}}{S_i + S_j}$$<\/span><\/p>\n<p>Dans cette \u00e9quation, le terme <span class=\"arabpsy-math-inline\" data-latex=\"C_{ij}\">$C_{ij}$<\/span> repr\u00e9sente la covariance minimale ou la somme des minima d&rsquo;abondances partag\u00e9es entre les deux profils pour l&rsquo;ensemble des <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> descripteurs recens\u00e9s. Ce terme fondamental est formellement d\u00e9fini par la sommation :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"C_{ij} = \\sum_{k=1}^{p} \\min(x_{ik}, x_{jk})\">$$C_{ij} = \\sum_{k=1}^{p} \\min(x_{ik}, x_{jk})$$<\/span><\/p>\n<p>Cette sommation isole rigoureusement la quantit\u00e9 d&rsquo;abondance effective que les deux \u00e9chantillons poss\u00e8dent en commun pour chaque variable <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>. Parall\u00e8lement, les termes <span class=\"arabpsy-math-inline\" data-latex=\"S_i\">$S_i$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"S_j\">$S_j$<\/span> repr\u00e9sentent respectivement les sommes marginales totales des abondances observ\u00e9es au sein de chacune des deux unit\u00e9s compar\u00e9es :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_i = \\sum_{k=1}^{p} x_{ik} \\quad \\text{et} \\quad S_j = \\sum_{k=1}^{p} x_{jk}\">$$S_i = \\sum_{k=1}^{p} x_{ik} \\quad \\text{et} \\quad S_j = \\sum_{k=1}^{p} x_{jk}$$<\/span><\/p>\n<p>La somme combin\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"S_i + S_j\">$S_i + S_j$<\/span> exprime donc l&rsquo;effectif cumul\u00e9 total de l&rsquo;ensemble des sp\u00e9cimens, intensit\u00e9s ou occurrences d\u00e9nombr\u00e9s sur les deux profils. Le ratio fractionnaire <span class=\"arabpsy-math-inline\" data-latex=\"(2 C_{ij}) \/ (S_i + S_j)\">$(2 C_{ij}) \/ (S_i + S_j)$<\/span> constitue quant \u00e0 lui l&rsquo;indice de similarit\u00e9 de Bray-Curtis (ou coefficient de similarit\u00e9 de S\u00f8rensen quantitatif). En soustrayant ce ratio de l&rsquo;unit\u00e9 scalaire, on obtient la mesure compl\u00e9mentaire directe : la dissimilarit\u00e9 de Bray-Curtis, dont les valeurs sont rigoureusement confin\u00e9es dans l&rsquo;intervalle ferm\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span>.<\/p>\n<h3>2.2 Formulations alg\u00e9briques alternatives et \u00e9quivalentes<\/h3>\n<p>Bien que la formulation canonique expos\u00e9e ci-dessus mette l&rsquo;accent sur la notion de partage d&rsquo;abondance via l&rsquo;op\u00e9rateur du minimum, il existe une formulation alg\u00e9brique alternative strictement \u00e9quivalente qui explicite la m\u00e9trique sous la forme d&rsquo;un rapport normalis\u00e9 de diff\u00e9rences absolues. Cette formulation d\u00e9montre que la dissimilarit\u00e9 de Bray-Curtis est l&rsquo;expression directe de la somme des \u00e9carts absolus composante par composante, divis\u00e9e par la somme totale des abondances des deux vecteurs :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij} = \\frac{\\sum_{k=1}^{p} |x_{ik} - x_{jk}|}{\\sum_{k=1}^{p} (x_{ik} + x_{jk})} = \\frac{\\sum_{k=1}^{p} |x_{ik} - x_{jk}|}{S_i + S_j}\">$$BC_{ij} = \\frac{\\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}|}{\\sum_{k=1}^{p} (x_{ik} + x_{jk})} = \\frac{\\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}|}{S_i + S_j}$$<\/span><\/p>\n<p>L&rsquo;\u00e9quivalence math\u00e9matique entre ces deux expressions se d\u00e9montre avec aisance en exploitant l&rsquo;identit\u00e9 fondamentale liant la valeur absolue, la somme et le minimum de deux r\u00e9els positifs. Pour toute paire de scalaires r\u00e9els <span class=\"arabpsy-math-inline\" data-latex=\"a ge 0\">$a ge 0$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"b ge 0\">$b ge 0$<\/span>, la relation alg\u00e9brique universelle suivante est v\u00e9rifi\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"|a - b| = (a + b) - 2 \\min(a, b)\">$$|a &#8211; b| = (a + b) &#8211; 2 \\min(a, b)$$<\/span><\/p>\n<p>En substituant cette propri\u00e9t\u00e9 au sein du num\u00e9rateur de l&rsquo;expression des diff\u00e9rences absolues, on retrouve imm\u00e9diatement la formulation canonique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\sum_{k=1}^{p} |x_{ik} - x_{jk}| = \\sum_{k=1}^{p} (x_{ik} + x_{jk}) - 2 \\sum_{k=1}^{p} \\min(x_{ik}, x_{jk}) = (S_i + S_j) - 2 C_{ij}\">$$\\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}| = \\sum_{k=1}^{p} (x_{ik} + x_{jk}) &#8211; 2 \\sum_{k=1}^{p} \\min(x_{ik}, x_{jk}) = (S_i + S_j) &#8211; 2 C_{ij}$$<\/span><\/p>\n<p>En divisant chaque terme par <span class=\"arabpsy-math-inline\" data-latex=\"(S_i + S_j)\">$(S_i + S_j)$<\/span>, l&rsquo;expression se scinde en <span class=\"arabpsy-math-inline\" data-latex=\"1 - (2 C_{ij}) \/ (S_i + S_j)\">$1 &#8211; (2 C_{ij}) \/ (S_i + S_j)$<\/span>, confirmant l&rsquo;identit\u00e9 math\u00e9matique parfaite des deux approches.<\/p>\n<p>Cette \u00e9criture alg\u00e9brique d\u00e9voile la filiation \u00e9troite entre la mesure de Bray-Curtis et d&rsquo;autres coefficients d&rsquo;association historiques. En statistique num\u00e9rique, cette m\u00e9trique est formellement identique au coefficient de dissemblance de Czekanowski quantitatif ainsi qu&rsquo;\u00e0 l&rsquo;indice de Steinhaus appliqu\u00e9 aux vecteurs d&rsquo;abondance. De plus, lorsque les donn\u00e9es d&rsquo;entr\u00e9e sont pr\u00e9alablement r\u00e9duites \u00e0 des variables binaires d&rsquo;occurrence (o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"x_{ik} in {0, 1}\">$x_{ik} in {0, 1}$<\/span>), la formule se simplifie instantan\u00e9ment pour donner le coefficient classique de S\u00f8rensen pour donn\u00e9es qualitatives, soulignant la continuit\u00e9 \u00e9l\u00e9gante entre les m\u00e9triques de pr\u00e9sence-absence et leurs g\u00e9n\u00e9ralisations quantitatives continues.<\/p>\n<h3>2.3 Hypoth\u00e8ses math\u00e9matiques sous-jacentes<\/h3>\n<p>L&rsquo;application rigoureuse du coefficient de dissimilarit\u00e9 de Bray-Curtis repose sur un ensemble de contraintes et de postulats math\u00e9matiques stricts qu&rsquo;il est indispensable de v\u00e9rifier lors du traitement de donn\u00e9es exp\u00e9rimentales. La condition premi\u00e8re et incontournable est l&rsquo;exigence de positivit\u00e9 stricte ou de nullit\u00e9 des donn\u00e9es d&rsquo;entr\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\forall i, \\forall k, \\quad x_{ik} ge 0\">$$\\forall i, \\forall k, \\quad x_{ik} ge 0$$<\/span><\/p>\n<p>L&rsquo;algorithme de Bray-Curtis est structurellement incompatible avec les valeurs n\u00e9gatives. Si une variable pr\u00e9sente des scores n\u00e9gatifs (comme cela peut survenir apr\u00e8s un centrage statistique autour de la moyenne, ou lors de l&rsquo;utilisation d&rsquo;\u00e9chelles psychom\u00e9triques bipolaires s&rsquo;\u00e9tendant de -3 \u00e0 +3), la signification g\u00e9om\u00e9trique et \u00e9cologique de l&rsquo;op\u00e9rateur <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{ik}, x_{jk})\">$\\min(x_{ik}, x_{jk})$<\/span> s&rsquo;effondre totalement, conduisant \u00e0 des dissimilarit\u00e9s qui peuvent sortir de l&rsquo;intervalle canonique <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span> et perdre toute coh\u00e9rence interpr\u00e9tative.<\/p>\n<p>Une seconde caract\u00e9ristique math\u00e9matique cruciale r\u00e9side dans la sensibilit\u00e9 structurelle de la m\u00e9trique \u00e0 l&rsquo;\u00e9chelle absolue des abondances brutes. Contrairement \u00e0 des coefficients normalis\u00e9s par descripteur (tels que la distance de Gower ou le coefficient de corr\u00e9lation de Pearson), le calcul de Bray-Curtis accorde un poids math\u00e9matique directement proportionnel \u00e0 l&rsquo;amplitude num\u00e9rique brute de chaque descripteur. Par cons\u00e9quent, une esp\u00e8ce ou un item clinique affichant des comptages se chiffrant en centaines ou en milliers d&rsquo;unit\u00e9s exercera une influence \u00e9crasante sur la valeur finale de <span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij}\">$BC_{ij}$<\/span>, rel\u00e9guant les variables \u00e0 faibles effectifs \u00e0 un r\u00f4le marginal, voire n\u00e9gligeable, dans l&rsquo;\u00e9valuation de la ressemblance globale.<\/p>\n<p>Enfin, le comportement asymptotique de la m\u00e9trique face aux descripteurs nuls partag\u00e9s m\u00e9rite une attention particuli\u00e8re. Si pour une composante <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> donn\u00e9e, les deux \u00e9chantillons pr\u00e9sentent une valeur nulle (<span class=\"arabpsy-math-inline\" data-latex=\"x_{ik} = 0\">$x_{ik} = 0$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"x_{jk} = 0\">$x_{jk} = 0$<\/span>), cette composante contribue pour exactement 0 au num\u00e9rateur (<span class=\"arabpsy-math-inline\" data-latex=\"|0 - 0| = 0\">$|0 &#8211; 0| = 0$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"min(0, 0) = 0\">$min(0, 0) = 0$<\/span>) ainsi qu&rsquo;au d\u00e9nominateur (<span class=\"arabpsy-math-inline\" data-latex=\"0 + 0 = 0\">$0 + 0 = 0$<\/span>). Les doubles z\u00e9ros n&rsquo;alt\u00e8rent donc en rien les sommes marginales ni la covariance minimale. Il en r\u00e9sulte que le calcul de Bray-Curtis est strictement invariant face \u00e0 l&rsquo;adjonction arbitraire de descripteurs totalement absents des deux unit\u00e9s compar\u00e9es, une propri\u00e9t\u00e9 fondatrice qui le distingue radicalement des distances spatiales conventionnelles.<\/p>\n<h2>3. Propri\u00e9t\u00e9s m\u00e9triques et caract\u00e9ristiques de l&rsquo;espace de dissimilarit\u00e9<\/h2>\n<h3>3.1 Non-respect de l&rsquo;in\u00e9galit\u00e9 triangulaire et statut de semi-m\u00e9trique<\/h3>\n<p>En topologie g\u00e9n\u00e9rale et en g\u00e9om\u00e9trie m\u00e9trique, pour qu&rsquo;une fonction de dissimilarit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"d(x, y)\">$d(x, y)$<\/span> soit formellement qualifi\u00e9e de v\u00e9ritable m\u00e9trique ou de distance au sens strict, elle doit satisfaire simultan\u00e9ment \u00e0 quatre axiomes fondamentaux pour tous points <span class=\"arabpsy-math-inline\" data-latex=\"x, y, z\">$x, y, z$<\/span> appartenant \u00e0 l&rsquo;ensemble consid\u00e9r\u00e9 :<\/p>\n<ul>\n<li><strong>La non-n\u00e9gativit\u00e9 :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"d(x, y) ge 0\">$d(x, y) ge 0$<\/span><\/li>\n<li><strong>L&rsquo;identit\u00e9 des indiscernables :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"d(x, y) = 0 iff x = y\">$d(x, y) = 0 iff x = y$<\/span><\/li>\n<li><strong>La sym\u00e9trie :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"d(x, y) = d(y, x)\">$d(x, y) = d(y, x)$<\/span><\/li>\n<li><strong>L&rsquo;in\u00e9galit\u00e9 triangulaire :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"d(x, z) le d(x, y) + d(y, z)\">$d(x, z) le d(x, y) + d(y, z)$<\/span><\/li>\n<\/ul>\n<p>Si la dissimilarit\u00e9 de Bray-Curtis satisfait pleinement aux trois premiers crit\u00e8res (sous r\u00e9serve que l&rsquo;on traite des profils contenant au moins une observation non nulle), elle \u00e9choue \u00e0 respecter syst\u00e9matiquement le quatri\u00e8me axiome fondamental : l&rsquo;in\u00e9galit\u00e9 triangulaire. Pour cette raison, la mesure de Bray-Curtis n&rsquo;est pas une m\u00e9trique pure, mais poss\u00e8de le statut math\u00e9matique de <strong>semi-m\u00e9trique<\/strong> (ou pseudo-m\u00e9trique de ressemblance).<\/p>\n<p>L&rsquo;\u00e9chec de l&rsquo;in\u00e9galit\u00e9 triangulaire peut \u00eatre ais\u00e9ment d\u00e9montr\u00e9 par un contre-exemple num\u00e9rique simple. Consid\u00e9rons trois profils univari\u00e9s ou multivari\u00e9s r\u00e9duits \u00e0 trois esp\u00e8ces distinctes, repr\u00e9sent\u00e9s par les vecteurs d&rsquo;abondance suivants :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{x} = (10, 0, 0), \\quad \\mathbf{y} = (10, 10, 0), \\quad \\mathbf{z} = (0, 10, 10)\">$$\\mathbf{x} = (10, 0, 0), \\quad \\mathbf{y} = (10, 10, 0), \\quad \\mathbf{z} = (0, 10, 10)$$<\/span><\/p>\n<p>Calculons successivement les dissimilarit\u00e9s de Bray-Curtis pour chaque paire :<\/p>\n<p>Pour la paire <span class=\"arabpsy-math-inline\" data-latex=\"(\\mathbf{x}, \\mathbf{y})\">$(\\mathbf{x}, \\mathbf{y})$<\/span> : les sommes valent <span class=\"arabpsy-math-inline\" data-latex=\"S_x = 10\">$S_x = 10$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"S_y = 20\">$S_y = 20$<\/span>. Le minimum commun pour l&rsquo;esp\u00e8ce 1 est <span class=\"arabpsy-math-inline\" data-latex=\"min(10, 10) = 10\">$min(10, 10) = 10$<\/span>. D&rsquo;o\u00f9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC(\\mathbf{x}, \\mathbf{y}) = 1 - \\frac{2 \\times 10}{10 + 20} = 1 - \\frac{20}{30} = \\frac{1}{3} \\approx 0{,}333\">$$BC(\\mathbf{x}, \\mathbf{y}) = 1 &#8211; \\frac{2 \\times 10}{10 + 20} = 1 &#8211; \\frac{20}{30} = \\frac{1}{3} \\approx 0{,}333$$<\/span><\/p>\n<p>Pour la paire <span class=\"arabpsy-math-inline\" data-latex=\"(\\mathbf{y}, \\mathbf{z})\">$(\\mathbf{y}, \\mathbf{z})$<\/span> : les sommes valent <span class=\"arabpsy-math-inline\" data-latex=\"S_y = 20\">$S_y = 20$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"S_z = 20\">$S_z = 20$<\/span>. Le minimum commun pour l&rsquo;esp\u00e8ce 2 est <span class=\"arabpsy-math-inline\" data-latex=\"min(10, 10) = 10\">$min(10, 10) = 10$<\/span>. D&rsquo;o\u00f9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC(\\mathbf{y}, \\mathbf{z}) = 1 - \\frac{2 \\times 10}{20 + 20} = 1 - \\frac{20}{40} = \\frac{1}{2} = 0{,}500\">$$BC(\\mathbf{y}, \\mathbf{z}) = 1 &#8211; \\frac{2 \\times 10}{20 + 20} = 1 &#8211; \\frac{20}{40} = \\frac{1}{2} = 0{,}500$$<\/span><\/p>\n<p>Pour la paire <span class=\"arabpsy-math-inline\" data-latex=\"(\\mathbf{x}, \\mathbf{z})\">$(\\mathbf{x}, \\mathbf{z})$<\/span> : les vecteurs n&rsquo;ont rigoureusement aucune esp\u00e8ce en commun. Le num\u00e9rateur de similarit\u00e9 est nul (<span class=\"arabpsy-math-inline\" data-latex=\"C_{xz} = 0\">$C_{xz} = 0$<\/span>). D&rsquo;o\u00f9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC(\\mathbf{x}, \\mathbf{z}) = 1 - \\frac{0}{10 + 20} = 1{,}000\">$$BC(\\mathbf{x}, \\mathbf{z}) = 1 &#8211; \\frac{0}{10 + 20} = 1{,}000$$<\/span><\/p>\n<p>Examinons \u00e0 pr\u00e9sent l&rsquo;in\u00e9galit\u00e9 triangulaire pour le trajet direct entre <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}\">$\\mathbf{x}$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{z}\">$\\mathbf{z}$<\/span> en passant par l&rsquo;interm\u00e9diaire <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{y}\">$\\mathbf{y}$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC(\\mathbf{x}, \\mathbf{y}) + BC(\\mathbf{y}, \\mathbf{z}) = \\frac{1}{3} + \\frac{1}{2} = \\frac{5}{6} \\approx 0{,}833\">$$BC(\\mathbf{x}, \\mathbf{y}) + BC(\\mathbf{y}, \\mathbf{z}) = \\frac{1}{3} + \\frac{1}{2} = \\frac{5}{6} \\approx 0{,}833$$<\/span><\/p>\n<p>Or, nous constatons que :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC(\\mathbf{x}, \\mathbf{z}) = 1{,}000 &gt; 0{,}833 = BC(\\mathbf{x}, \\mathbf{y}) + BC(\\mathbf{y}, \\mathbf{z})\">$$BC(\\mathbf{x}, \\mathbf{z}) = 1{,}000 &gt; 0{,}833 = BC(\\mathbf{x}, \\mathbf{y}) + BC(\\mathbf{y}, \\mathbf{z})$$<\/span><\/p>\n<p>L&rsquo;in\u00e9galit\u00e9 triangulaire est ici formellement viol\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"1{,}000 &gt; 0{,}833\">$1{,}000 &gt; 0{,}833$<\/span>). Les cons\u00e9quences g\u00e9om\u00e9triques de cette violation sont majeures : l&rsquo;espace relationnel engendr\u00e9 par une matrice de dissimilarit\u00e9 de Bray-Curtis est non euclidien. Il est impossible de projeter directement et sans distorsion l&rsquo;ensemble des points d&rsquo;une telle matrice dans un espace euclidien rigide \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> dimensions sans g\u00e9n\u00e9rer des valeurs propres n\u00e9gatives lors d&rsquo;une analyse en coordonn\u00e9es principales (PCoA). C&rsquo;est pr\u00e9cis\u00e9ment cette non-lin\u00e9arit\u00e9 fondamentale qui impose le recours pr\u00e9f\u00e9rentiel \u00e0 des m\u00e9thodes d&rsquo;ordination non m\u00e9triques, au premier rang desquelles figure le positionnement multidimensionnel non m\u00e9trique (NMDS).<\/p>\n<h3>3.2 Propri\u00e9t\u00e9 d&rsquo;invariance et ind\u00e9pendance des doubles z\u00e9ros<\/h3>\n<p>L&rsquo;une des propri\u00e9t\u00e9s les plus admirables et \u00e9cologiquement signifiantes de la dissimilarit\u00e9 de Bray-Curtis est son invariance absolue face \u00e0 l&rsquo;adjonction de doubles z\u00e9ros, souvent qualifi\u00e9e de propri\u00e9t\u00e9 de neutralit\u00e9 des absences conjointes. Dans l&rsquo;analyse quantitative de communaut\u00e9s naturelles ou de profils cliniques, l&rsquo;absence simultan\u00e9e d&rsquo;une esp\u00e8ce ou d&rsquo;un sympt\u00f4me chez deux individus distincts constitue une information hautement ambigu\u00eb. Deux for\u00eats temp\u00e9r\u00e9es \u00e9chantillonn\u00e9es en hiver ne pr\u00e9sentent ni palmiers tropicaux, ni cactus d\u00e9sertiques, ni mousses alpines ; cela justifie-t-il pour autant de consid\u00e9rer que ces deux for\u00eats sont hautement similaires parce qu&rsquo;elles partagent l&rsquo;absence de milliers d&rsquo;esp\u00e8ces tropicales ?<\/p>\n<p>En \u00e9cologie num\u00e9rique, r\u00e9pondre affirmativement \u00e0 cette question constitue une erreur m\u00e9thodologique gravissime connue sous le nom de <strong>paradoxe des doubles z\u00e9ros<\/strong>. Les m\u00e9triques g\u00e9om\u00e9triques conventionnelles, telles que la c\u00e9l\u00e8bre distance euclidienne ou la distance de Manhattan non normalis\u00e9e, sont lourdement sujettes \u00e0 ce biais :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"d_{\\text{Euclidienne}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\sqrt{\\sum_{k=1}^{p} (x_{ik} - x_{jk})^2}\">$$d_{\\text{Euclidienne}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\sqrt{\\sum_{k=1}^{p} (x_{ik} &#8211; x_{jk})^2}$$<\/span><\/p>\n<p>Si l&rsquo;on enrichit la matrice d&rsquo;\u00e9tude de dix esp\u00e8ces suppl\u00e9mentaires totalement absentes des deux sites (<span class=\"arabpsy-math-inline\" data-latex=\"x_{ik} = 0\">$x_{ik} = 0$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"x_{jk} = 0\">$x_{jk} = 0$<\/span>), chaque terme <span class=\"arabpsy-math-inline\" data-latex=\"(0 - 0)^2 = 0\">$(0 &#8211; 0)^2 = 0$<\/span> s&rsquo;additionne sous le radical sans alt\u00e9rer la distance absolue. Toutefois, si l&rsquo;on calcule des indices de ressemblance proportionnels ou des distances statistiques normalis\u00e9es, la pr\u00e9sence massive de z\u00e9ros fait gonfler artificiellement la proximit\u00e9 globale per\u00e7ue entre les sites. La distance euclidienne tend \u00e0 regrouper fallacieusement deux sites d\u00e9sertiques ou deux stations polaires biologiquement tr\u00e8s pauvres simplement parce qu&rsquo;ils partagent une immense majorit\u00e9 de z\u00e9ros taxonomiques.<\/p>\n<p>La dissimilarit\u00e9 de Bray-Curtis neutralise math\u00e9matiquement ce probl\u00e8me par construction. Reprenons sa formulation alternative :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij} = \\frac{\\sum_{k=1}^{p} |x_{ik} - x_{jk}|}{\\sum_{k=1}^{p} (x_{ik} + x_{jk})}\">$$BC_{ij} = \\frac{\\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}|}{\\sum_{k=1}^{p} (x_{ik} + x_{jk})}$$<\/span><\/p>\n<p>Si pour un descripteur arbitraire <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, <span class=\"arabpsy-math-inline\" data-latex=\"x_{ik} = 0\">$x_{ik} = 0$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"x_{jk} = 0\">$x_{jk} = 0$<\/span>, nous obtenons <span class=\"arabpsy-math-inline\" data-latex=\"|0 - 0| = 0\">$|0 &#8211; 0| = 0$<\/span> au num\u00e9rateur et <span class=\"arabpsy-math-inline\" data-latex=\"0 + 0 = 0\">$0 + 0 = 0$<\/span> au d\u00e9nominateur. Ce descripteur n&rsquo;apporte strictement aucune contribution \u00e0 la valeur de la fraction. Que la matrice initiale contienne 10 esp\u00e8ces ou que l&rsquo;on y adjoigne 10 000 esp\u00e8ces hypoth\u00e9tiques jamais observ\u00e9es sur le terrain d&rsquo;\u00e9tude, la valeur num\u00e9rique de <span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij}\">$BC_{ij}$<\/span> demeure rigoureusement identique au milli\u00e8me pr\u00e8s. Bray-Curtis focalise donc l&rsquo;analyse sur les pr\u00e9sences confirm\u00e9es et les intensit\u00e9s partag\u00e9es, \u00e9cartant avec une rigueur absolue le bruit d&rsquo;illusion d&rsquo;association induit par les absences ubiquitaires.<\/p>\n<h3>3.3 Sym\u00e9trie et bornes de l&rsquo;espace d&rsquo;association<\/h3>\n<p>Sur le plan structural, l&rsquo;op\u00e9rateur de dissimilarit\u00e9 de Bray-Curtis est caract\u00e9ris\u00e9 par une sym\u00e9trie math\u00e9matique parfaite. Pour toute paire de vecteurs d&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_i\">$\\mathbf{x}_i$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_j\">$\\mathbf{x}_j$<\/span>, la relation d&rsquo;invariance par permutation est strictement satisfaite :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij} = BC_{ji}\">$$BC_{ij} = BC_{ji}$$<\/span><\/p>\n<p>La preuve formelle d\u00e9coule de la commutativit\u00e9 inh\u00e9rente aux op\u00e9rateurs sous-jacents d&rsquo;addition et de minimisation dans l&rsquo;ensemble des nombres r\u00e9els :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\min(x_{ik}, x_{jk}) = \\min(x_{jk}, x_{ik}) implies C_{ij} = C_{ji}\">$$\\min(x_{ik}, x_{jk}) = \\min(x_{jk}, x_{ik}) implies C_{ij} = C_{ji}$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"x_{ik} + x_{jk} = x_{jk} + x_{ik} implies S_i + S_j = S_j + S_i\">$$x_{ik} + x_{jk} = x_{jk} + x_{ik} implies S_i + S_j = S_j + S_i$$<\/span><\/p>\n<p>Il s&rsquo;ensuit que la matrice d&rsquo;association de Bray-Curtis g\u00e9n\u00e9r\u00e9e pour un jeu de donn\u00e9es contenant <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> unit\u00e9s d&rsquo;\u00e9chantillonnage est une matrice carr\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{D}\">$\\mathbf{D}$<\/span> d&rsquo;ordre <span class=\"arabpsy-math-inline\" data-latex=\"n \\times n\">$n \\times n$<\/span>, dont la diagonale principale est exclusivement compos\u00e9e de z\u00e9ros (<span class=\"arabpsy-math-inline\" data-latex=\"BC_{ii} = 0\">$BC_{ii} = 0$<\/span> pour tout <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>), et qui est parfaitement sym\u00e9trique par rapport \u00e0 cette diagonale (<span class=\"arabpsy-math-inline\" data-latex=\"D_{ij} = D_{ji}\">$D_{ij} = D_{ji}$<\/span>). Cette propri\u00e9t\u00e9 de sym\u00e9trie garantit la coh\u00e9rence des structures topologiques d\u00e9riv\u00e9es lors des proc\u00e9dures de classification hi\u00e9rarchique ou d&rsquo;ordination globale.<\/p>\n<p>Concernant le domaine de d\u00e9finition et les bornes scalaires, la dissimilarit\u00e9 de Bray-Curtis est une fonction strictement born\u00e9e dans l&rsquo;intervalle r\u00e9el <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span>, \u00e0 la condition absolue que les vecteurs d&rsquo;entr\u00e9e contiennent au moins une valeur strictement positive (afin d&rsquo;\u00e9viter une ind\u00e9termination math\u00e9matique du type <span class=\"arabpsy-math-inline\" data-latex=\"0\/0\">$0\/0$<\/span> au d\u00e9nominateur). Lorsque deux profils sont quantitativement identiques dans leur structure et leurs amplitudes, la dissimilarit\u00e9 atteint sa borne inf\u00e9rieure absolue : <span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij} = 0\">$BC_{ij} = 0$<\/span>. \u00c0 l&rsquo;autre extr\u00e9mit\u00e9, lorsque deux profils ne partagent absolument aucun descripteur commun (<span class=\"arabpsy-math-inline\" data-latex=\"C_{ij} = 0\">$C_{ij} = 0$<\/span>), le num\u00e9rateur de similarit\u00e9 s&rsquo;annule compl\u00e8tement et la dissimilarit\u00e9 atteint sa borne sup\u00e9rieure unitaire : <span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij} = 1\">$BC_{ij} = 1$<\/span>.<\/p>\n<p>Le comportement analytique de la m\u00e9trique lors de variations unilat\u00e9rales extr\u00eames s&rsquo;av\u00e8re particuli\u00e8rement stable. Consid\u00e9rons un syst\u00e8me o\u00f9 un \u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> conserve une composition fixe tandis que l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span> voit l&rsquo;abondance d&rsquo;un taxon exclusif tendre vers l&rsquo;infini (<span class=\"arabpsy-math-inline\" data-latex=\"x_{jm} to +\\infty\">$x_{jm} to +\\infty$<\/span>). Dans ce cas de figure, le d\u00e9nominateur <span class=\"arabpsy-math-inline\" data-latex=\"S_i + S_j\">$S_i + S_j$<\/span> cro\u00eet de mani\u00e8re illimit\u00e9e tandis que le num\u00e9rateur des minima partag\u00e9s <span class=\"arabpsy-math-inline\" data-latex=\"C_{ij}\">$C_{ij}$<\/span> demeure plafonn\u00e9 par les abondances limit\u00e9es du vecteur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_i\">$\\mathbf{x}_i$<\/span>. Le ratio de similarit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"(2 C_{ij}) \/ (S_i + S_j)\">$(2 C_{ij}) \/ (S_i + S_j)$<\/span> tend alors asymptotiquement vers z\u00e9ro, et la dissimilarit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij}\">$BC_{ij}$<\/span> converge de fa\u00e7on fluide et monotone vers sa borne unitaire $1$, sans jamais diverger vers l&rsquo;infini comme le ferait une distance euclidienne non normalis\u00e9e.<\/p>\n<h2>4. Guide de calcul pas \u00e0 pas : Exemple empirique complet<\/h2>\n<h3>4.1 Pr\u00e9sentation du jeu de donn\u00e9es exp\u00e9rimentales<\/h3>\n<p>Afin d&rsquo;ancrer les concepts th\u00e9oriques dans une r\u00e9alit\u00e9 calculatoire tangible, d\u00e9roulons l&rsquo;int\u00e9gralit\u00e9 des op\u00e9rations alg\u00e9briques sur un jeu de donn\u00e9es empiriques concret. Imaginons une \u00e9tude comparative portant sur deux profils d&rsquo;observation distincts. Selon le prisme disciplinaire choisi, ces profils peuvent repr\u00e9senter :<\/p>\n<ul>\n<li>En <strong>\u00e9cologie v\u00e9g\u00e9tale<\/strong> : deux parcelles foresti\u00e8res distinctes (Site 1 et Site 2) au sein desquelles on a recens\u00e9 le nombre d&rsquo;individus matures pour cinq essences ligneuses majeures (Arbres A, B, C, D et E).<\/li>\n<li>En <strong>psychopathologie quantitative<\/strong> : deux profils cliniques de patients (Sujet A et Sujet B) \u00e9valu\u00e9s sur une \u00e9chelle de s\u00e9v\u00e9rit\u00e9 symptomatique gradu\u00e9e de 0 \u00e0 50 pour cinq dimensions psychopathologiques distinctes (Sympt\u00f4mes A, B, C, D et E).<\/li>\n<\/ul>\n<p>Structurons notre matrice de donn\u00e9es brutes contenant les vecteurs d&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_1\">$\\mathbf{x}_1$<\/span> (repr\u00e9sentant le Site 1 ou Sujet A) et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}_2\">$\\mathbf{x}_2$<\/span> (repr\u00e9sentant le Site 2 ou Sujet B) pour les <span class=\"arabpsy-math-inline\" data-latex=\"p = 5\">$p = 5$<\/span> variables d&rsquo;int\u00e9r\u00eat :<\/p>\n<ul>\n<li><strong>Variable A :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"x_{1A} = 10 \\quad ; \\quad x_{2A} = 20\">$x_{1A} = 10 \\quad ; \\quad x_{2A} = 20$<\/span><\/li>\n<li><strong>Variable B :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"x_{1B} = 25 \\quad ; \\quad x_{2B} = 10\">$x_{1B} = 25 \\quad ; \\quad x_{2B} = 10$<\/span><\/li>\n<li><strong>Variable C :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"x_{1C} = 0 \\quad ; \\quad x_{2C} = 5\">$x_{1C} = 0 \\quad ; \\quad x_{2C} = 5$<\/span><\/li>\n<li><strong>Variable D :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"x_{1D} = 15 \\quad ; \\quad x_{2D} = 15\">$x_{1D} = 15 \\quad ; \\quad x_{2D} = 15$<\/span><\/li>\n<li><strong>Variable E :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"x_{1E} = 50 \\quad ; \\quad x_{2E} = 30\">$x_{1E} = 50 \\quad ; \\quad x_{2E} = 30$<\/span><\/li>\n<\/ul>\n<p>Ces donn\u00e9es pr\u00e9sentent un \u00e9ventail repr\u00e9sentatif de situations concr\u00e8tes : des variables o\u00f9 le profil 1 domine (Variable B et E), des variables o\u00f9 le profil 2 domine (Variable A), une variable marqu\u00e9e par une stricte \u00e9galit\u00e9 d&rsquo;abondance (Variable D), et une variable illustrant une absence unilat\u00e9rale (Variable C, totalement absente du profil 1 mais observ\u00e9e dans le profil 2).<\/p>\n<h3>4.2 Calcul m\u00e9thodique des sommes d&rsquo;abondance marginales<\/h3>\n<p>La premi\u00e8re phase calculatoire de la proc\u00e9dure de Bray-Curtis consiste \u00e0 \u00e9tablir les sommes marginales vectorielles pour chacun des deux profils d&rsquo;observation. Cette op\u00e9ration permet de quantifier la masse d&rsquo;\u00e9chantillonnage globale ou l&rsquo;intensit\u00e9 totale manifest\u00e9e par chaque unit\u00e9 statistique ind\u00e9pendamment de l&rsquo;autre.<\/p>\n<p>Calculons en premier lieu la somme d&rsquo;abondance marginale pour le profil 1, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"S_1\">$S_1$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_1 = \\sum_{k=A}^{E} x_{1k} = x_{1A} + x_{1B} + x_{1C} + x_{1D} + x_{1E}\">$$S_1 = \\sum_{k=A}^{E} x_{1k} = x_{1A} + x_{1B} + x_{1C} + x_{1D} + x_{1E}$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_1 = 10 + 25 + 0 + 15 + 50 = 100\">$$S_1 = 10 + 25 + 0 + 15 + 50 = 100$$<\/span><\/p>\n<p>Le profil 1 totalise donc exactement 100 unit\u00e9s cumul\u00e9es (qu&rsquo;il s&rsquo;agisse de tiges d&rsquo;arbres inventori\u00e9es ou de points cumul\u00e9s de charge symptomatique globale).<\/p>\n<p>Calculons ensuite la somme d&rsquo;abondance marginale pour le profil 2, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"S_2\">$S_2$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_2 = \\sum_{k=A}^{E} x_{2k} = x_{2A} + x_{2B} + x_{2C} + x_{2D} + x_{2E}\">$$S_2 = \\sum_{k=A}^{E} x_{2k} = x_{2A} + x_{2B} + x_{2C} + x_{2D} + x_{2E}$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_2 = 20 + 10 + 5 + 15 + 30 = 80\">$$S_2 = 20 + 10 + 5 + 15 + 30 = 80$$<\/span><\/p>\n<p>Le profil 2 totalise ainsi 80 unit\u00e9s cumul\u00e9es. L&rsquo;effectif global du second profil est donc inf\u00e9rieur de 20 % \u00e0 celui du premier, une situation classique de variation de l&rsquo;effort d&rsquo;\u00e9chantillonnage ou de la r\u00e9activit\u00e9 globale du sujet.<\/p>\n<p>La sommation marginale combin\u00e9e, qui constituera le d\u00e9nominateur de notre fraction canonique, se r\u00e9sout imm\u00e9diatement par l&rsquo;addition des deux scalaires obtenus :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_1 + S_2 = 100 + 80 = 180\">$$S_1 + S_2 = 100 + 80 = 180$$<\/span><\/p>\n<p>La masse totale des observations d\u00e9ploy\u00e9es sur l&rsquo;ensemble du syst\u00e8me comparatif s&rsquo;\u00e9tablit pr\u00e9cis\u00e9ment \u00e0 180 unit\u00e9s.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/bray_curtis2.png\" alt=\"Bray-Curtis Dissimilarity\" \/><figcaption>Bray-Curtis Dissimilarity<\/figcaption><\/figure>\n<h3>4.3 Identification et sommation des valeurs minimales partag\u00e9es<\/h3>\n<p>La deuxi\u00e8me \u00e9tape analytique constitue le c\u0153ur m\u00e9canique de la m\u00e9trique : elle vise \u00e0 extraire pour chaque variable <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> la composante commune aux deux profils, mat\u00e9rialis\u00e9e par l&rsquo;infimum ou la valeur minimale partag\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{1k}, x_{2k})\">$\\min(x_{1k}, x_{2k})$<\/span>. Proc\u00e9dons \u00e0 cette extraction m\u00e9thodique item par item :<\/p>\n<ul>\n<li><strong>Pour la Variable A :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{1A}, x_{2A}) = \\min(10, 20) = 10\">$\\min(x_{1A}, x_{2A}) = \\min(10, 20) = 10$<\/span>. Les deux profils partagent 10 unit\u00e9s de cette variable ; les 10 unit\u00e9s exc\u00e9dentaires observ\u00e9es sur le profil 2 constituent une divergence unilat\u00e9rale.<\/li>\n<li><strong>Pour la Variable B :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{1B}, x_{2B}) = \\min(25, 10) = 10\">$\\min(x_{1B}, x_{2B}) = \\min(25, 10) = 10$<\/span>. La valeur commune est born\u00e9e \u00e0 10 unit\u00e9s par la limitation du profil 2.<\/li>\n<li><strong>Pour la Variable C :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{1C}, x_{2C}) = \\min(0, 5) = 0\">$\\min(x_{1C}, x_{2C}) = \\min(0, 5) = 0$<\/span>. Le profil 1 pr\u00e9sentant un effectif nul, il ne peut exister aucune quantit\u00e9 commune partag\u00e9e pour cette dimension. La contribution \u00e0 la similarit\u00e9 est rigoureusement nulle.<\/li>\n<li><strong>Pour la Variable D :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{1D}, x_{2D}) = \\min(15, 15) = 15\">$\\min(x_{1D}, x_{2D}) = \\min(15, 15) = 15$<\/span>. La co\u00efncidence num\u00e9rique \u00e9tant parfaite, l&rsquo;int\u00e9gralit\u00e9 de l&rsquo;abondance de cette variable (15 unit\u00e9s) entre dans le domaine du partage effectif.<\/li>\n<li><strong>Pour la Variable E :<\/strong> <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{1E}, x_{2E}) = \\min(50, 30) = 30\">$\\min(x_{1E}, x_{2E}) = \\min(50, 30) = 30$<\/span>. Le volume commun maximal est d\u00e9fini par l&rsquo;effectif du profil 2, soit 30 unit\u00e9s partag\u00e9es.<\/li>\n<\/ul>\n<p>Nous pouvons \u00e0 pr\u00e9sent sommer l&rsquo;ensemble de ces contributions partielles pour obtenir la valeur consolid\u00e9e du param\u00e8tre fondamental <span class=\"arabpsy-math-inline\" data-latex=\"C_{12}\">$C_{12}$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"C_{12} = \\sum_{k=A}^{E} \\min(x_{1k}, x_{2k}) = 10 + 10 + 0 + 15 + 30 = 65\">$$C_{12} = \\sum_{k=A}^{E} \\min(x_{1k}, x_{2k}) = 10 + 10 + 0 + 15 + 30 = 65$$<\/span><\/p>\n<p>Le volume cumul\u00e9 des abondances positivement partag\u00e9es entre les profils 1 et 2 s&rsquo;\u00e9l\u00e8ve \u00e0 65 unit\u00e9s. Dans l&rsquo;architecture de la formule canonique, ce terme doit \u00eatre doubl\u00e9 afin de rendre compte de la pr\u00e9sence de ces unit\u00e9s dans chacun des deux vecteurs d&rsquo;observation :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"2 \\times C_{12} = 2 \\times 65 = 130\">$$2 \\times C_{12} = 2 \\times 65 = 130$$<\/span><\/p>\n<h3>4.4 R\u00e9solution finale et d\u00e9duction du coefficient de dissemblance<\/h3>\n<p>Disposant de toutes les composantes scalaires n\u00e9cessaires, nous pouvons proc\u00e9der \u00e0 la r\u00e9solution arithm\u00e9tique de la similarit\u00e9, puis \u00e0 la d\u00e9duction imm\u00e9diate de la dissimilarit\u00e9 de Bray-Curtis.<\/p>\n<p>Calculons en premier lieu l&rsquo;indice de similarit\u00e9 de S\u00f8rensen quantitatif :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Similarit\u00e9}_{12} = \\frac{2 C_{12}}{S_1 + S_2} = \\frac{130}{180}\">$$\\text{Similarit\u00e9}_{12} = \\frac{2 C_{12}}{S_1 + S_2} = \\frac{130}{180}$$<\/span><\/p>\n<p>En simplifiant cette fraction par leur plus grand commun diviseur (10), nous obtenons :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Similarit\u00e9}_{12} = \\frac{13}{18} \\approx 0{,}7222\">$$\\text{Similarit\u00e9}_{12} = \\frac{13}{18} \\approx 0{,}7222$$<\/span><\/p>\n<p>Les deux profils partagent ainsi environ 72,22 % de leur structure quantitative globale d&rsquo;abondance.<\/p>\n<p>La dissimilarit\u00e9 de Bray-Curtis se d\u00e9duit directement en soustrayant ce taux de l&rsquo;unit\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{12} = 1 - \\text{Similarit\u00e9}_{12} = 1 - \\frac{130}{180} = \\frac{180 - 130}{180} = \\frac{50}{180} = \\frac{5}{18} \\approx 0{,}2778\">$$BC_{12} = 1 &#8211; \\text{Similarit\u00e9}_{12} = 1 &#8211; \\frac{130}{180} = \\frac{180 &#8211; 130}{180} = \\frac{50}{180} = \\frac{5}{18} \\approx 0{,}2778$$<\/span><\/p>\n<p>Pour attester de la parfaite rigueur de notre calcul, appliquons la seconde formulation alg\u00e9brique bas\u00e9e sur la somme des diff\u00e9rences absolues :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\sum_{k=A}^{E} |x_{1k} - x_{2k}| = |10 - 20| + |25 - 10| + |0 - 5| + |15 - 15| + |50 - 30|\">$$\\sum_{k=A}^{E} |x_{1k} &#8211; x_{2k}| = |10 &#8211; 20| + |25 &#8211; 10| + |0 &#8211; 5| + |15 &#8211; 15| + |50 &#8211; 30|$$<\/span><\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\sum_{k=A}^{E} |x_{1k} - x_{2k}| = |-10| + |15| + |-5| + |0| + |20| = 10 + 15 + 5 + 0 + 20 = 50\">$$\\sum_{k=A}^{E} |x_{1k} &#8211; x_{2k}| = |-10| + |15| + |-5| + |0| + |20| = 10 + 15 + 5 + 0 + 20 = 50$$<\/span><\/p>\n<p>Divisons cette somme d&rsquo;\u00e9carts absolus par la masse d&rsquo;abondance totale <span class=\"arabpsy-math-inline\" data-latex=\"(S_1 + S_2 = 180)\">$(S_1 + S_2 = 180)$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{12} = \\frac{50}{180} = \\frac{5}{18} \\approx 0{,}2778\">$$BC_{12} = \\frac{50}{180} = \\frac{5}{18} \\approx 0{,}2778$$<\/span><\/p>\n<p>La concordance algorithmique est totale et parfaite. La dissimilarit\u00e9 de Bray-Curtis entre les deux unit\u00e9s s&rsquo;\u00e9tablit rigoureusement \u00e0 <strong>0,2778<\/strong> (soit une divergence structurelle de 27,78 % de l&rsquo;ensemble du profil d&rsquo;abondance partag\u00e9).<\/p>\n<h2>5. Interpr\u00e9tation analytique du continuum de valeurs de 0 \u00e0 1<\/h2>\n<h3>5.1 Signification th\u00e9orique et clinique d&rsquo;une dissimilarit\u00e9 nulle (BC = 0)<\/h3>\n<p>L&rsquo;obtention d&rsquo;une valeur de dissimilarit\u00e9 de Bray-Curtis strictement \u00e9gale \u00e0 z\u00e9ro (<span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij} = 0\">$BC_{ij} = 0$<\/span>) correspond \u00e0 une similarit\u00e9 parfaite (<span class=\"arabpsy-math-inline\" data-latex=\"1{,}00\">$1{,}00$<\/span> ou 100 %) entre les deux unit\u00e9s d&rsquo;\u00e9chantillonnage compar\u00e9es. Sur le plan formel, cette condition n&rsquo;est satisfaite que si et seulement si les deux vecteurs d&rsquo;observation pr\u00e9sentent une \u00e9galit\u00e9 num\u00e9rique absolue composante par composante :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij} = 0 iff \\forall k in {1, dots, p}, \\quad x_{ik} = x_{jk}\">$$BC_{ij} = 0 iff \\forall k in {1, dots, p}, \\quad x_{ik} = x_{jk}$$<\/span><\/p>\n<p>Dans un contexte phyto\u00e9cologique, un coefficient <span class=\"arabpsy-math-inline\" data-latex=\"BC = 0\">$BC = 0$<\/span> indique que les deux parcelles recens\u00e9es abritent rigoureusement les m\u00eames esp\u00e8ces avec des nombres d&rsquo;individus ou des biomasses strictement identiques pour chacune d&rsquo;entre elles. En psychom\u00e9trie clinique, cela signifie que deux patients expriment un tableau psychopathologique parfaitement superposable, manifestant la m\u00eame intensit\u00e9 exacte sur chacun des items de l&rsquo;\u00e9chelle d&rsquo;\u00e9valuation dimensionnelle.<\/p>\n<p>Il est primordial de souligner une nuance conceptuelle critique : la dissimilarit\u00e9 de Bray-Curtis n&rsquo;est pas invariante par changement d&rsquo;\u00e9chelle scalaire unilat\u00e9rale. Cela signifie que si le profil <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span> est un multiple scalaire strict du profil <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> (par exemple <span class=\"arabpsy-math-inline\" data-latex=\"x_{jk} = 2 \\times x_{ik}\">$x_{jk} = 2 \\times x_{ik}$<\/span> pour tout <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, traduisant une colin\u00e9arit\u00e9 vectorielle parfaite avec un angle nul entre les vecteurs), la dissimilarit\u00e9 de Bray-Curtis <strong>ne sera pas \u00e9gale \u00e0 z\u00e9ro<\/strong>. Calculons la valeur th\u00e9orique pour un tel cas : si <span class=\"arabpsy-math-inline\" data-latex=\"x_{jk} = 2 x_{ik}\">$x_{jk} = 2 x_{ik}$<\/span>, alors <span class=\"arabpsy-math-inline\" data-latex=\"S_j = 2 S_i\">$S_j = 2 S_i$<\/span>, et <span class=\"arabpsy-math-inline\" data-latex=\"\\min(x_{ik}, x_{jk}) = x_{ik}\">$\\min(x_{ik}, x_{jk}) = x_{ik}$<\/span>, donc <span class=\"arabpsy-math-inline\" data-latex=\"C_{ij} = S_i\">$C_{ij} = S_i$<\/span>. Le d\u00e9nominateur vaut <span class=\"arabpsy-math-inline\" data-latex=\"S_i + 2 S_i = 3 S_i\">$S_i + 2 S_i = 3 S_i$<\/span>. La similarit\u00e9 s&rsquo;\u00e9l\u00e8ve \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"(2 S_i) \/ (3 S_i) = 2\/3 \\approx 0{,}667\">$(2 S_i) \/ (3 S_i) = 2\/3 \\approx 0{,}667$<\/span>, et la dissimilarit\u00e9 de Bray-Curtis s&rsquo;\u00e9tablit \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"1 - 2\/3 = 1\/3 \\approx 0{,}333\">$1 &#8211; 2\/3 = 1\/3 \\approx 0{,}333$<\/span>. Cette sensibilit\u00e9 aux amplitudes absolues distingue fondamentalement Bray-Curtis des mesures angulaires ou de corr\u00e9lation (comme le cosinus de Salton ou le coefficient de Pearson), qui consid\u00e8rent deux profils homoth\u00e9tiques comme strictement identiques.<\/p>\n<h3>5.2 Compr\u00e9hension de la dissimilarit\u00e9 maximale unitaire (BC = 1)<\/h3>\n<p>\u00c0 l&rsquo;autre extr\u00e9mit\u00e9 du continuum d&rsquo;association, une valeur unitaire (<span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij} = 1\">$BC_{ij} = 1$<\/span>) signale une divergence totale et absolue entre les deux entit\u00e9s observ\u00e9es. Alg\u00e9briquement, cette valeur limite est atteinte d\u00e8s lors que la covariance minimale partag\u00e9e est strictement nulle :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij} = 1 iff C_{ij} = 0 iff \\forall k in {1, dots, p}, \\quad \\min(x_{ik}, x_{jk}) = 0\">$$BC_{ij} = 1 iff C_{ij} = 0 iff \\forall k in {1, dots, p}, \\quad \\min(x_{ik}, x_{jk}) = 0$$<\/span><\/p>\n<p>Cette condition signifie que pour chaque descripteur <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> \u00e9valu\u00e9, l&rsquo;un au moins des deux \u00e9chantillons pr\u00e9sente une abondance nulle. En d&rsquo;autres termes, les deux unit\u00e9s compar\u00e9es ne partagent strictement aucune composante positive en commun ; leurs r\u00e9pertoires d&rsquo;observation sont totalement disjoints :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Support}(\\mathbf{x}_i) \\cap \\text{Support}(\\mathbf{x}_j) = \\emptyset\">$$\\text{Support}(\\mathbf{x}_i) \\cap \\text{Support}(\\mathbf{x}_j) = \\emptyset$$<\/span><\/p>\n<p>En \u00e9cologie des communaut\u00e9s, un indice de <span class=\"arabpsy-math-inline\" data-latex=\"1{,}00\">$1{,}00$<\/span> caract\u00e9rise deux biotopes appartenant \u00e0 des biomes totalement irr\u00e9conciliables, par exemple un \u00e9cosyst\u00e8me benthique abyssal et une toundra arctique, o\u00f9 aucun taxon n&rsquo;est capable de franchir la barri\u00e8re physiologique de l&rsquo;autre habitat. Sur le plan de la nosologie psychiatrique, un <span class=\"arabpsy-math-inline\" data-latex=\"BC = 1\">$BC = 1$<\/span> mat\u00e9rialise une divergence symptomatique absolue : l&rsquo;un des patients pr\u00e9sente un syndrome d&rsquo;agitation maniaque pur avec logorrh\u00e9e et fuite des id\u00e9es, tandis que le second souffre d&rsquo;un coma catatonique akin\u00e9tique sans aucune manifestation clinique commune. Ils \u00e9voluent dans des sph\u00e8res pathologiques distinctes ne partageant aucune dimension commune.<\/p>\n<h3>5.3 Graduation et lecture critique des valeurs interm\u00e9diaires<\/h3>\n<p>Dans la pratique empirique des sciences naturelles et sociales, les valeurs extr\u00eames ($0$ ou $1$) ne repr\u00e9sentent qu&rsquo;une fraction marginale des observations ; la quasi-totalit\u00e9 des confrontations num\u00e9riques s&rsquo;inscrit au sein du continuum ouvert <span class=\"arabpsy-math-inline\" data-latex=\"]0, 1[\">$]0, 1[$<\/span>. L&rsquo;interpr\u00e9tation de ces valeurs interm\u00e9diaires requiert une solide compr\u00e9hension de la distribution sous-jacente des donn\u00e9es et des seuils conventionnels \u00e9tablis par la litt\u00e9rature sp\u00e9cialis\u00e9e.<\/p>\n<p>Bien que toute grille d&rsquo;\u00e9valuation comporte une part d&rsquo;arbitraire contextuel, les praticiens de l&rsquo;\u00e9cologie num\u00e9rique s&rsquo;accordent g\u00e9n\u00e9ralement sur la stratification heuristique suivante :<\/p>\n<ul>\n<li><strong><span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij} &lt; 0{,}30\">$BC_{ij} &lt; 0{,}30$<\/span> (Forte similarit\u00e9 \/ Faible dissimilarit\u00e9) :<\/strong> Les profils compar\u00e9s pr\u00e9sentent un recouvrement taxonomique ou structural pr\u00e9pond\u00e9rant. Les esp\u00e8ces dominantes sont partag\u00e9es dans des proportions comparables, et les \u00e9carts constat\u00e9s r\u00e9sultent principalement de fluctuations d\u00e9mographiques mineures, de variabilit\u00e9 intra-stationnelle ou de l\u00e9g\u00e8res impr\u00e9cisions de mesure.<\/li>\n<li><strong><span class=\"arabpsy-math-inline\" data-latex=\"0{,}30 le BC_{ij} le 0{,}70\">$0{,}30 le BC_{ij} le 0{,}70$<\/span> (Zone de divergence mod\u00e9r\u00e9e \u00e0 marqu\u00e9e) :<\/strong> Cette plage m\u00e9diane refl\u00e8te des transitions \u00e9cologiques significatives (\u00e9cotones) ou des profils symptomatiques partiellement concordants. Les entit\u00e9s partagent un socle substantiel de variables communes, mais exhibent des inversions de dominance notables ou des remplacements progressifs de descripteurs interm\u00e9diaires.<\/li>\n<li><strong><span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij} &gt; 0{,}70\">$BC_{ij} &gt; 0{,}70$<\/span> (Forte dissimilarit\u00e9 \/ Renouvellement majeur) :<\/strong> Les communaut\u00e9s biologiques op\u00e8rent un renouvellement quasi-complet (diversit\u00e9 b\u00eata tr\u00e8s \u00e9lev\u00e9e). Seules quelques esp\u00e8ces ubiquistes, pionni\u00e8res ou opportunistes maintiennent un lien de continuit\u00e9 entre les sites. En clinique, cela caract\u00e9rise des sous-types nosologiques clairement distincts partageant tout au plus quelques sympt\u00f4mes non sp\u00e9cifiques de d\u00e9tresse psychologique g\u00e9n\u00e9rale.<\/li>\n<\/ul>\n<p>Il importe toutefois de manipuler ces seuils avec pr\u00e9caution. La valeur absolue de <span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij}\">$BC_{ij}$<\/span> d\u00e9pend fortement de la richesse totale en esp\u00e8ces du jeu de donn\u00e9es, du degr\u00e9 d&rsquo;agr\u00e9gation taxonomique (niveau de l&rsquo;esp\u00e8ce, du genre ou de la famille) et du protocole d&rsquo;\u00e9chantillonnage employ\u00e9. Dans les syst\u00e8mes hautement diversifi\u00e9s comme les for\u00eats n\u00e9otropicales ou les profils m\u00e9tag\u00e9nomiques du microbiote bact\u00e9rien comprenant plusieurs milliers d&rsquo;unit\u00e9s taxonomiques op\u00e9rationnelles (OTU), la valeur moyenne de Bray-Curtis entre sites naturels avoisine couramment des niveaux \u00e9lev\u00e9s (<span class=\"arabpsy-math-inline\" data-latex=\"0{,}75\">$0{,}75$<\/span> \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"0{,}90\">$0{,}90$<\/span>) en raison de la raret\u00e9 inh\u00e9rente de l&rsquo;immense majorit\u00e9 des micro-organismes recens\u00e9s.<\/p>\n<h2>6. Comparaison m\u00e9thodologique : Bray-Curtis versus autres m\u00e9triques<\/h2>\n<h3>6.1 Confrontation avec la distance euclidienne et ses variantes<\/h3>\n<p>Pour saisir toute la pertinence m\u00e9thodologique de la dissimilarit\u00e9 de Bray-Curtis, il est \u00e9clairant de la confronter aux propri\u00e9t\u00e9s de la distance euclidienne canonique. Dans l&rsquo;histoire de la morphom\u00e9trie et de la statistique descriptive classique, la distance euclidienne s&rsquo;est impos\u00e9e naturellement par analogie directe avec la g\u00e9om\u00e9trie de l&rsquo;espace physique tridimensionnel :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"d_E(\\mathbf{x}_i, \\mathbf{x}_j) = \\sqrt{\\sum_{k=1}^{p} (x_{ik} - x_{jk})^2}\">$$d_E(\\mathbf{x}_i, \\mathbf{x}_j) = \\sqrt{\\sum_{k=1}^{p} (x_{ik} &#8211; x_{jk})^2}$$<\/span><\/p>\n<p>Toutefois, lorsqu&rsquo;elle est appliqu\u00e9e \u00e0 des matrices d&rsquo;abondances biologiques ou de r\u00e9pertoires comportementaux hautement h\u00e9t\u00e9rog\u00e8nes, la distance euclidienne engendre ce que Legendre et Gallagher ont formellement qualifi\u00e9 de <strong>paradoxe de l&rsquo;espace euclidien<\/strong>. Ce paradoxe d\u00e9coule de deux tares conceptuelles majeures :<\/p>\n<p>En premier lieu, comme examin\u00e9 pr\u00e9c\u00e9demment, la distance euclidienne viole le principe de neutralit\u00e9 des absences conjointes. Pour la distance euclidienne, le fait que deux sites ne poss\u00e8dent aucun individu de l&rsquo;esp\u00e8ce <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> (<span class=\"arabpsy-math-inline\" data-latex=\"0 - 0 = 0\">$0 &#8211; 0 = 0$<\/span>) contribue \u00e0 les rapprocher g\u00e9om\u00e9triquement avec le m\u00eame poids que s&rsquo;ils poss\u00e9daient tous deux 100 individus de cette esp\u00e8ce (<span class=\"arabpsy-math-inline\" data-latex=\"100 - 100 = 0\">$100 &#8211; 100 = 0$<\/span>). Cela conduit \u00e0 une absurdit\u00e9 \u00e9cologique : des stations totalement arides et st\u00e9riles sont regroup\u00e9es comme hautement apparent\u00e9es uniquement par la gr\u00e2ce de leur vide biologique partag\u00e9.<\/p>\n<p>En second lieu, l&rsquo;\u00e9l\u00e9vation au carr\u00e9 des \u00e9carts diff\u00e9rentiels <span class=\"arabpsy-math-inline\" data-latex=\"(x_{ik} - x_{jk})^2\">$(x_{ik} &#8211; x_{jk})^2$<\/span> conf\u00e8re une pr\u00e9pond\u00e9rance disproportionn\u00e9e aux descripteurs hypersurcharg\u00e9s, tout en \u00e9crasant litt\u00e9ralement le signal apport\u00e9 par les descripteurs mod\u00e9r\u00e9s ou rares. Une diff\u00e9rence de 50 individus sur une esp\u00e8ce arbor\u00e9e hyper-dominante p\u00e8sera <span class=\"arabpsy-math-inline\" data-latex=\"(50)^2 = 2500\">$(50)^2 = 2500$<\/span> unit\u00e9s dans la somme sous le radical, tandis qu&rsquo;une divergence de pr\u00e9sence\/absence sur 5 esp\u00e8ces rares g\u00e9n\u00e9rera un impact d\u00e9risoire de <span class=\"arabpsy-math-inline\" data-latex=\"5 \\times (1 - 0)^2 = 5\">$5 \\times (1 &#8211; 0)^2 = 5$<\/span> unit\u00e9s. La distance euclidienne tend donc \u00e0 masquer l&rsquo;organisation globale de la communaut\u00e9 au profit de fluctuations d\u00e9mographiques al\u00e9atoires survenant sur les quelques taxons ubiquistes. En divisant la somme des \u00e9carts absolus par la masse totale des abondances, Bray-Curtis assure une mise \u00e0 l&rsquo;\u00e9chelle proportionnelle qui pr\u00e9serve la texture \u00e9cologique et psychologique fine des jeux de donn\u00e9es r\u00e9els.<\/p>\n<h3>6.2 Distinction avec la dissimilarit\u00e9 de Jaccard<\/h3>\n<p>La dissimilarit\u00e9 de Jaccard est sans conteste l&rsquo;indice de ressemblance qualitatif le plus c\u00e9l\u00e8bre et le plus ancien de la litt\u00e9rature biologique. Con\u00e7u par le botaniste suisse Paul Jaccard en 1901, cet indice s&rsquo;applique originellement \u00e0 des donn\u00e9es binaires de pr\u00e9sence-absence cod\u00e9es sous forme d&rsquo;une table de contingence bivari\u00e9e comprenant :<\/p>\n<ul>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"a\">$a$<\/span> : le nombre d&rsquo;esp\u00e8ces pr\u00e9sentes simultan\u00e9ment sur les deux sites (<span class=\"arabpsy-math-inline\" data-latex=\"1, 1\">$1, 1$<\/span>)<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"b\">$b$<\/span> : le nombre d&rsquo;esp\u00e8ces pr\u00e9sentes uniquement sur le premier site (<span class=\"arabpsy-math-inline\" data-latex=\"1, 0\">$1, 0$<\/span>)<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"c\">$c$<\/span> : le nombre d&rsquo;esp\u00e8ces pr\u00e9sentes uniquement sur le second site (<span class=\"arabpsy-math-inline\" data-latex=\"0, 1\">$0, 1$<\/span>)<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"d\">$d$<\/span> : le nombre d&rsquo;esp\u00e8ces absentes des deux sites (<span class=\"arabpsy-math-inline\" data-latex=\"0, 0\">$0, 0$<\/span>)<\/li>\n<\/ul>\n<p>L&rsquo;indice qualitatif de Jaccard s&rsquo;\u00e9nonce selon le ratio bien connu :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"J_{\\text{binaire}} = 1 - \\frac{a}{a + b + c} = \\frac{b + c}{a + b + c}\">$$J_{\\text{binaire}} = 1 &#8211; \\frac{a}{a + b + c} = \\frac{b + c}{a + b + c}$$<\/span><\/p>\n<p>Bien que Jaccard respecte l&rsquo;invariance aux doubles absences (le terme <span class=\"arabpsy-math-inline\" data-latex=\"d\">$d$<\/span> est totalement exclu du calcul), il pr\u00e9sente une limite majeure : il est strictement aveugle aux variations d&rsquo;abondance quantitative. Pour rem\u00e9dier \u00e0 cette restriction, l&rsquo;\u00e9cologue Milan Ru\u017ei\u010dka a formul\u00e9 en 1958 l&rsquo;extension quantitative continue de l&rsquo;indice de Jaccard, connue sous le nom de distance de Ru\u017ei\u010dka :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"R_{ij} = 1 - \\frac{\\sum_{k=1}^{p} \\min(x_{ik}, x_{jk})}{\\sum_{k=1}^{p} \\max(x_{ik}, x_{jk})} = \\frac{\\sum_{k=1}^{p} |x_{ik} - x_{jk}|}{\\sum_{k=1}^{p} \\max(x_{ik}, x_{jk})}\">$$R_{ij} = 1 &#8211; \\frac{\\sum_{k=1}^{p} \\min(x_{ik}, x_{jk})}{\\sum_{k=1}^{p} \\max(x_{ik}, x_{jk})} = \\frac{\\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}|}{\\sum_{k=1}^{p} \\max(x_{ik}, x_{jk})}$$<\/span><\/p>\n<p>La confrontation entre Bray-Curtis et Ru\u017ei\u010dka d\u00e9voile une parent\u00e9 math\u00e9matique tr\u00e8s intime. Rappelons que dans Bray-Curtis, le d\u00e9nominateur est la somme simple des abondances <span class=\"arabpsy-math-inline\" data-latex=\"S_i + S_j\">$S_i + S_j$<\/span>. Or, pour toute paire de nombres positifs, nous savons que :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"x_{ik} + x_{jk} = \\min(x_{ik}, x_{jk}) + \\max(x_{ik}, x_{jk})\">$$x_{ik} + x_{jk} = \\min(x_{ik}, x_{jk}) + \\max(x_{ik}, x_{jk})$$<\/span><\/p>\n<p>Il existe par cons\u00e9quent une relation analytique exacte permettant de convertir directement l&rsquo;indice de similarit\u00e9 de Bray-Curtis (not\u00e9 ici <span class=\"arabpsy-math-inline\" data-latex=\"S_{BC}\">$S_{BC}$<\/span>) en similarit\u00e9 de Ru\u017ei\u010dka (not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"S_R\">$S_R$<\/span>) :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_R = \\frac{S_{BC}}{2 - S_{BC}} iff S_{BC} = \\frac{2 S_R}{1 + S_R}\">$$S_R = \\frac{S_{BC}}{2 &#8211; S_{BC}} iff S_{BC} = \\frac{2 S_R}{1 + S_R}$$<\/span><\/p>\n<p>Sur le plan topologique, la distance de Ru\u017ei\u010dka poss\u00e8de un avantage th\u00e9orique formel sur Bray-Curtis : elle satisfait pleinement l&rsquo;in\u00e9galit\u00e9 triangulaire et constitue donc une <strong>v\u00e9ritable m\u00e9trique<\/strong> dans l&rsquo;espace multidimensionnel. Toutefois, dans les analyses d&rsquo;ordination empirique (notamment la NMDS), les r\u00e9sultats fournis par Bray-Curtis et Ru\u017ei\u010dka sont virtuellement indiscernables en raison de la relation de monotonicit\u00e9 stricte qui lie ces deux fonctions sur l&rsquo;intervalle <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span>. Bray-Curtis demeure n\u00e9anmoins historiquement plus pl\u00e9biscit\u00e9 en raison de la simplicit\u00e9 intuitive de son d\u00e9nominateur marginal.<\/p>\n<h3>6.3 Comparaison avec la distance de Manhattan et l&rsquo;indice de Morisita-Horn<\/h3>\n<p>Une autre perspective comparative essentielle consiste \u00e0 analyser les relations entre Bray-Curtis, la m\u00e9trique L1 de Manhattan et les coefficients bas\u00e9s sur les probabilit\u00e9s de rencontre inter-individuelle comme l&rsquo;indice de Morisita-Horn.<\/p>\n<p>La distance de Manhattan (ou distance g\u00e9od\u00e9sique en damier) s&rsquo;exprime comme la somme brute des \u00e9carts absolus non pond\u00e9r\u00e9s :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"d_M(\\mathbf{x}_i, \\mathbf{x}_j) = \\sum_{k=1}^{p} |x_{ik} - x_{jk}|\">$$d_M(\\mathbf{x}_i, \\mathbf{x}_j) = \\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}|$$<\/span><\/p>\n<p>En observant cette \u00e9quation, on constate imm\u00e9diatement que le num\u00e9rateur de la formulation alternative de Bray-Curtis est rigoureusement identique \u00e0 la distance de Manhattan. Ainsi, la dissimilarit\u00e9 de Bray-Curtis peut \u00eatre th\u00e9oriquement interpr\u00e9t\u00e9e comme une <strong>distance de Manhattan normalis\u00e9e \u00e0 l&rsquo;\u00e9chelle locale<\/strong> par la biomasse ou l&rsquo;effectif d&rsquo;\u00e9chantillonnage cumul\u00e9 de la paire consid\u00e9r\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"S_i + S_j\">$S_i + S_j$<\/span>). La distance de Manhattan classique n&rsquo;est pas born\u00e9e sup\u00e9rieurement ; elle diverge proportionnellement \u00e0 la richesse brute du milieu, rendant impossible la comparaison d&rsquo;assemblages floristiques pr\u00e9lev\u00e9s \u00e0 des intensit\u00e9s divergentes.<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;indice de <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/BF02518851\">Morisita (1959)<\/a>, adapt\u00e9 ult\u00e9rieurement par Horn (1966) pour des densit\u00e9s num\u00e9riques continues, aborde la mesure de ressemblance sous un angle probabiliste inspir\u00e9 de l&rsquo;indice de diversit\u00e9 de Simpson :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"C_H = \\frac{2 \\sum_{k=1}^{p} x_{ik} x_{jk}}{\\left( \\frac{\\sum x_{ik}^2}{S_i^2} + \\frac{\\sum x_{jk}^2}{S_j^2} \\right) S_i S_j}\">$$C_H = \\frac{2 \\sum_{k=1}^{p} x_{ik} x_{jk}}{\\left( \\frac{\\sum x_{ik}^2}{S_i^2} + \\frac{\\sum x_{jk}^2}{S_j^2} \\right) S_i S_j}$$<\/span><\/p>\n<p>L&rsquo;indice de Morisita-Horn quantifie la probabilit\u00e9 que deux individus tir\u00e9s au sort au hasard, l&rsquo;un dans l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> et l&rsquo;autre dans l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span>, appartiennent \u00e0 la m\u00eame esp\u00e8ce, standardis\u00e9e par les probabilit\u00e9s d&rsquo;identit\u00e9 intra-site. Cet indice poss\u00e8de une r\u00e9sistance math\u00e9matique exceptionnelle aux variations d&rsquo;effort d&rsquo;\u00e9chantillonnage et aux fluctuations de taille globale des collectes, surpassant sur ce point pr\u00e9cis la m\u00e9trique de Bray-Curtis. En contrepartie, Morisita-Horn est hypersensible \u00e0 la composition de l&rsquo;unique esp\u00e8ce la plus dominante du lot, son terme quadratique (<span class=\"arabpsy-math-inline\" data-latex=\"x_{ik}^2\">$x_{ik}^2$<\/span>) amplifiant l&rsquo;influence des taxons surrepr\u00e9sent\u00e9s bien plus violemment que l&rsquo;op\u00e9rateur lin\u00e9aire du minimum pr\u00e9sent dans Bray-Curtis.<\/p>\n<h2>7. Sensibilit\u00e9 aux transformations math\u00e9matiques et aux biais d&rsquo;\u00e9chantillonnage<\/h2>\n<h3>7.1 Effet de dominance des variables hypersurcharg\u00e9es<\/h3>\n<p>En d\u00e9pit de sa remarquable robustesse face aux absences conjointes, le calcul non corrig\u00e9 de la dissimilarit\u00e9 de Bray-Curtis pr\u00e9sente une vuln\u00e9rabilit\u00e9 inh\u00e9rente li\u00e9e \u00e0 sa formulation lin\u00e9aire : sa propension m\u00e9canique \u00e0 laisser les descripteurs \u00e0 tr\u00e8s fortes abondances dominer unilat\u00e9ralement le r\u00e9sultat num\u00e9rique. Dans les syst\u00e8mes naturels, la distribution des esp\u00e8ces au sein d&rsquo;une communaut\u00e9 suit presque universellement une loi log-normale ou une distribution g\u00e9om\u00e9trique en s\u00e9rie de Motomura, caract\u00e9ris\u00e9e par une infime minorit\u00e9 d&rsquo;esp\u00e8ces hyper-abondantes c\u00f4toyant une longue tra\u00eene d&rsquo;esp\u00e8ces rares ou accidentelles.<\/p>\n<p>Imaginons un inventaire entomologique o\u00f9 une esp\u00e8ce de moustique opportuniste est captur\u00e9e \u00e0 raison de 5 000 individus sur le site 1 et 2 000 individus sur le site 2, tandis que vingt esp\u00e8ces rares de col\u00e9opt\u00e8res sp\u00e9cialistes pr\u00e9sentent des abondances variant entre 0 et 3 individus par station. L&rsquo;\u00e9cart sur l&rsquo;esp\u00e8ce commune g\u00e9n\u00e9rera \u00e0 lui seul une diff\u00e9rence absolue de <span class=\"arabpsy-math-inline\" data-latex=\"|5000 - 2000| = 3000\">$|5000 &#8211; 2000| = 3000$<\/span> unit\u00e9s, \u00e9crasant totalement la contribution combin\u00e9e des vingt col\u00e9opt\u00e8res rares dont la divergence cumul\u00e9e n&rsquo;exc\u00e9dera pas 30 unit\u00e9s. La m\u00e9trique de Bray-Curtis brute sera quasiment \u00e0 100 % le reflet de la dynamique d\u00e9mographique du moustique, ignorant l&rsquo;information bio-indicatrice pr\u00e9cieuse v\u00e9hicul\u00e9e par le cort\u00e8ge des esp\u00e8ces rares.<\/p>\n<p>Ce ph\u00e9nom\u00e8ne d&rsquo;\u00e9crasement se retrouve \u00e0 l&rsquo;identique dans l&rsquo;\u00e9valuation psychom\u00e9trique ou comportementale. Si un questionnaire int\u00e8gre des variables comportementales \u00e0 fr\u00e9quence quotidienne \u00e9lev\u00e9e (par exemple le nombre de minutes pass\u00e9es sur un smartphone, s&rsquo;\u00e9levant \u00e0 300 ou 400 unit\u00e9s) aux c\u00f4t\u00e9s de cotations d&rsquo;actes d&rsquo;auto-agression ou d&rsquo;attaques de panique (variant de 1 \u00e0 3 \u00e9pisodes par semaine), le score de Bray-Curtis calcul\u00e9 sur les donn\u00e9es brutes sera un artefact presque exclusif du temps d&rsquo;\u00e9cran, occultant la divergence critique touchant les variables psychiatriques s\u00e9v\u00e8res. Cette distorsion impose la mise en \u0153uvre syst\u00e9matique de transformations math\u00e9matiques pr\u00e9alables.<\/p>\n<h3>7.2 Transformations pr\u00e9alables recommand\u00e9es : racine carr\u00e9e et log<\/h3>\n<p>Pour restaurer l&rsquo;\u00e9quit\u00e9 informationnelle entre les composantes \u00e0 forte charge num\u00e9rique et les descripteurs de faible amplitude, la pratique biostatistique moderne recommande d&rsquo;appliquer une transformation monotone concave sur la matrice brute des abondances avant le calcul de la dissimilarit\u00e9 de Bray-Curtis. Trois approches dominent la m\u00e9thodologie quantitative :<\/p>\n<p>La transformation par <strong>racine carr\u00e9e<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"y_{ik} = \\sqrt{x_{ik}}\">$y_{ik} = \\sqrt{x_{ik}}$<\/span>) constitue la mod\u00e9ration la plus douce. Elle r\u00e9duit l&rsquo;\u00e9cart relatif entre ordres de grandeur sans alt\u00e9rer la hi\u00e9rarchie fondamentale des descripteurs. Pour des asym\u00e9tries de distribution plus prononc\u00e9es, la transformation par <strong>racine quatri\u00e8me<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"y_{ik} = x_{ik}^{1\/4} = \\sqrt{\\sqrt{x_{ik}}}\">$y_{ik} = x_{ik}^{1\/4} = \\sqrt{\\sqrt{x_{ik}}}$<\/span>) est couramment d\u00e9sign\u00e9e comme le standard d&rsquo;excellence dans les analyses d&rsquo;\u00e9cologie benthique marine selon les protocoles \u00e9tablis par Clarke et Warwick.<\/p>\n<p>La seconde alternative majeure est la <strong>transformation logarithmique<\/strong>, formalis\u00e9e sous la forme :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"y_{ik} = \\log(x_{ik} + 1) \\quad \\text{ou} \\quad y_{ik} = \\log_{10}(c \\cdot x_{ik} + 1)\">$$y_{ik} = \\log(x_{ik} + 1) \\quad \\text{ou} \\quad y_{ik} = \\log_{10}(c \\cdot x_{ik} + 1)$$<\/span><\/p>\n<p>L&rsquo;adjonction de la constante pseudo-compte <span class=\"arabpsy-math-inline\" data-latex=\"+1\">$+1$<\/span> est strictement indispensable pour g\u00e9rer harmonieusement les valeurs nulles du tableau d&rsquo;origine, garantissant que <span class=\"arabpsy-math-inline\" data-latex=\"log(0 + 1) = log(1) = 0\">$log(0 + 1) = log(1) = 0$<\/span>, et pr\u00e9servant ainsi la nature des absences fondamentales. L&rsquo;op\u00e9rateur logarithmique comprime drastiquement les \u00e9chelles num\u00e9riques \u00e9tendues, permettant aux taxons d&rsquo;abondance interm\u00e9diaire et faible d&rsquo;exercer un poids statistique proportionnel \u00e0 leur valeur diagnostique.<\/p>\n<p>Enfin, la <strong>transformation de double standardisation du Wisconsin<\/strong> (Wisconsin double standardization) repr\u00e9sente un protocole d&rsquo;harmonisation particuli\u00e8rement \u00e9l\u00e9gant. Elle consiste \u00e0 diviser dans un premier temps l&rsquo;abondance de chaque descripteur par le maximum observ\u00e9 pour ce descripteur \u00e0 travers l&rsquo;ensemble des \u00e9chantillons (alignant ainsi toutes les variables sur une \u00e9chelle continue de 0 \u00e0 1), puis \u00e0 diviser dans un second temps ces valeurs par la somme totale du profil de l&rsquo;\u00e9chantillon. L&rsquo;application de Bray-Curtis sur une matrice pr\u00e9trait\u00e9e par le Wisconsin neutralise \u00e0 la fois les disparit\u00e9s d&rsquo;\u00e9chelle entre variables et les distorsions li\u00e9es au volume global d&rsquo;\u00e9chantillonnage, fournissant l&rsquo;une des ordinations les plus stables et reproductibles de la statistique multivari\u00e9e.<\/p>\n<h3>7.3 Normalisation et standardisation des profils<\/h3>\n<p>Une question m\u00e9thodologique r\u00e9currente concerne l&rsquo;opportunit\u00e9 de transformer la matrice des comptages bruts en une matrice d&rsquo;abondances relatives ou de pourcentages intra-\u00e9chantillon avant de calculer la dissimilarit\u00e9 de Bray-Curtis. Cette standardisation marginale par ligne s&rsquo;exprime math\u00e9matiquement par :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"p_{ik} = \\frac{x_{ik}}{\\sum_{m=1}^{p} x_{im}} = \\frac{x_{ik}}{S_i}\">$$p_{ik} = \\frac{x_{ik}}{\\sum_{m=1}^{p} x_{im}} = \\frac{x_{ik}}{S_i}$$<\/span><\/p>\n<p>Lorsque cette op\u00e9ration est appliqu\u00e9e, la somme totale de chaque vecteur d&rsquo;observation devient rigoureusement \u00e9gale \u00e0 l&rsquo;unit\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\forall i, \\quad S_i^* = \\sum_{k=1}^{p} p_{ik} = 1{,}00\">$$\\forall i, \\quad S_i^* = \\sum_{k=1}^{p} p_{ik} = 1{,}00$$<\/span><\/p>\n<p>Injectons cette propri\u00e9t\u00e9 d&rsquo;\u00e9quilibrage dans la formule canonique de Bray-Curtis :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij}^* = 1 - \\frac{2 \\sum_{k=1}^{p} \\min(p_{ik}, p_{jk})}{1{,}00 + 1{,}00} = 1 - \\sum_{k=1}^{p} \\min(p_{ik}, p_{jk})\">$$BC_{ij}^* = 1 &#8211; \\frac{2 \\sum_{k=1}^{p} \\min(p_{ik}, p_{jk})}{1{,}00 + 1{,}00} = 1 &#8211; \\sum_{k=1}^{p} \\min(p_{ik}, p_{jk})$$<\/span><\/p>\n<p>Ou, de mani\u00e8re \u00e9quivalente via la formulation des diff\u00e9rences absolues :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"BC_{ij}^* = \\frac{1}{2} \\sum_{k=1}^{p} |p_{ik} - p_{jk}|\">$$BC_{ij}^* = \\frac{1}{2} \\sum_{k=1}^{p} |p_{ik} &#8211; p_{jk}|$$<\/span><\/p>\n<p>L&rsquo;\u00e9quation obtenue d\u00e9montre que l&rsquo;application de Bray-Curtis sur des profils de pourcentages standardis\u00e9s se simplifie en une <strong>demi-distance de Manhattan<\/strong> (\u00e9galement connue sous le vocable d&rsquo;indice de dissemblance de Whittaker). Les cons\u00e9quences conceptuelles de ce choix sont substantielles : la standardisation par les totaux \u00e9limine d\u00e9finitivement toute influence li\u00e9e \u00e0 la taille brute absolue des \u00e9chantillons. Si deux stations abritent des esp\u00e8ces dans des proportions strictement identiques mais que l&rsquo;une a fait l&rsquo;objet d&rsquo;un inventaire dix fois plus exhaustif que l&rsquo;autre (<span class=\"arabpsy-math-inline\" data-latex=\"N = 1000\">$N = 1000$<\/span> contre <span class=\"arabpsy-math-inline\" data-latex=\"N = 100\">$N = 100$<\/span>), leur dissimilarit\u00e9 de Bray-Curtis brute sera de <span class=\"arabpsy-math-inline\" data-latex=\"0{,}818\">$0{,}818$<\/span>, tandis que leur dissimilarit\u00e9 sur profils normalis\u00e9s s&rsquo;effondrera \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"0{,}000\">$0{,}000$<\/span>. Le chercheur doit donc arbitrer en pleine conscience : souhaite-t-il \u00e9valuer la divergence de la taille absolue des populations ou isoler la pure structure compositionnelle des proportions relatives ?<\/p>\n<h2>8. Impl\u00e9mentation algorithmique en langages R et Python<\/h2>\n<h3>8.1 Traitement sous R avec le package vegan<\/h3>\n<p>L&rsquo;environnement statistique open-source R s&rsquo;impose comme l&rsquo;\u00e9cosyst\u00e8me de r\u00e9f\u00e9rence historique pour l&rsquo;analyse \u00e9cologique et multivari\u00e9e, en grande partie gr\u00e2ce au d\u00e9veloppement continu du package sp\u00e9cialis\u00e9 <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=vegan\">vegan<\/a>, maintenu par Jari Oksanen et ses collaborateurs. Au sein de cette biblioth\u00e8que, la fonction ma\u00eetresse d\u00e9di\u00e9e au calcul des matrices d&rsquo;association est <code>vegdist<\/code>.<\/p>\n<p>L&rsquo;ex\u00e9cution du calcul de Bray-Curtis s&rsquo;op\u00e8re de fa\u00e7on hautement optimis\u00e9e. La fonction accepte directement en entr\u00e9e une matrice de donn\u00e9es ou un data frame o\u00f9 les lignes mat\u00e9rialisent les unit\u00e9s d&rsquo;\u00e9chantillonnage et les colonnes repr\u00e9sentent les diff\u00e9rents descripteurs mesur\u00e9s. L&rsquo;appel standardis\u00e9 requiert la sp\u00e9cification de l&rsquo;argument de m\u00e9thode sous la forme <code>method = \"bray\"<\/code>. Par d\u00e9faut, la routine prend en charge les abondances brutes, mais il est d&rsquo;usage courant d&rsquo;encapsuler l&rsquo;appel au sein de transformations pr\u00e9alables en utilisant la fonction compl\u00e9mentaire <code>decostand<\/code>.<\/p>\n<p>Pour des architectures de donn\u00e9es massives \u2014 telles que les tables d&rsquo;inventaire g\u00e9n\u00e9r\u00e9es par le s\u00e9quen\u00e7age \u00e0 haut d\u00e9bit d&rsquo;amplicons m\u00e9tag\u00e9nomiques (16S rRNA ou ITS) comprenant plusieurs milliers d&rsquo;\u00e9chantillons et des dizaines de milliers d&rsquo;OTU \u2014 la fonction <code>vegdist<\/code> d\u00e9montre une v\u00e9locit\u00e9 computationnelle remarquable en s&rsquo;appuyant sur des sous-routines internes compil\u00e9es en langage C, \u00e9vitant ainsi les goulots d&rsquo;\u00e9tranglement inh\u00e9rents aux boucles interpr\u00e9t\u00e9es en R pur.<\/p>\n<h3>8.2 Calcul en Python via Scipy et Scikit-Bio<\/h3>\n<p>Dans l&rsquo;\u00e9cosyst\u00e8me Python, d\u00e9di\u00e9 aux pipelines de science des donn\u00e9es, d&rsquo;apprentissage automatique et de bio-informatique computationnelle, la dissimilarit\u00e9 de Bray-Curtis b\u00e9n\u00e9ficie d&rsquo;impl\u00e9mentations standardis\u00e9es au sein des biblioth\u00e8ques fondamentales <a rel=\"noopener noreferrer\" href=\"https:\/\/scipy.org\/\">SciPy<\/a> et <a rel=\"noopener noreferrer\" href=\"https:\/\/scikit.bio\/\">Scikit-Bio<\/a>.<\/p>\n<p>Au niveau le plus \u00e9l\u00e9mentaire, le module spatial de SciPy offre la fonction <code>scipy.spatial.distance.braycurtis<\/code>. Cette fonction prend en arguments deux tableaux unidimensionnels NumPy de m\u00eame dimension et renvoie la valeur scalaire de la dissimilarit\u00e9 de Bray-Curtis. L&rsquo;impl\u00e9mentation algorithmique sous-jacente de SciPy exploite formellement la version de la somme des \u00e9carts absolus normalis\u00e9e par la somme totale des \u00e9l\u00e9ments :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{braycurtis}(u, v) = \\frac{\\sum |u_k - v_k|}{\\sum |u_k + v_k|}\">$$\\text{braycurtis}(u, v) = \\frac{\\sum |u_k &#8211; v_k|}{\\sum |u_k + v_k|}$$<\/span><\/p>\n<p>Il convient de noter que si les vecteurs d&rsquo;entr\u00e9e contiennent des valeurs n\u00e9gatives, SciPy n&rsquo;\u00e9met pas d&rsquo;erreur bloquante mais applique la formule litt\u00e9rale ci-dessus, ce qui confirme l&rsquo;imp\u00e9rieuse n\u00e9cessit\u00e9 pour le praticien de veiller \u00e0 la non-n\u00e9gativit\u00e9 stricte de ses donn\u00e9es avant tout appel de routine.<\/p>\n<p>Pour le calcul de matrices compl\u00e8tes d&rsquo;association \u00e0 grande \u00e9chelle au sein de structures de donn\u00e9es multidimensionnelles, la biblioth\u00e8que <code>scikit-bio<\/code> s&rsquo;av\u00e8re nettement plus adapt\u00e9e. Sa fonction <code>skbio.diversity.beta.beta_diversity<\/code> permet d&rsquo;ex\u00e9cuter le calcul parall\u00e9lis\u00e9 sur des matrices creuses (objets <code>scipy.sparse<\/code>), optimisant radicalement l&#8217;empreinte m\u00e9moire vive lors du traitement de r\u00e9pertoires hautement clairsem\u00e9s contenant jusqu&rsquo;\u00e0 95 % de z\u00e9ros structurels, typiques des profils de s\u00e9quen\u00e7age m\u00e9tag\u00e9nomique moderne.<\/p>\n<h3>8.3 Cr\u00e9ation d&rsquo;une fonction personnalis\u00e9e en script<\/h3>\n<p>La compr\u00e9hension intime d&rsquo;un algorithme passe id\u00e9alement par sa reconstruction int\u00e9grale sans le recours \u00e0 des biblioth\u00e8ques externes opaques. La mise en \u0153uvre d&rsquo;une routine vectoris\u00e9e calculant la matrice sym\u00e9trique compl\u00e8te de Bray-Curtis illustre la transition directe de la th\u00e9orie math\u00e9matique vers le code op\u00e9rationnel.<\/p>\n<p>Dans une perspective algorithmique vectoris\u00e9e, le calcul optimis\u00e9 pour une paire de vecteurs <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{u}\">$\\mathbf{u}$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{v}\">$\\mathbf{v}$<\/span> se condense en trois \u00e9tapes cl\u00e9s :<\/p>\n<ul>\n<li>\u00c9valuation de la condition d&rsquo;arr\u00eat : v\u00e9rification de la positivit\u00e9 des donn\u00e9es et confirmation que la somme combin\u00e9e des deux vecteurs est strictement sup\u00e9rieure \u00e0 z\u00e9ro (afin de pr\u00e9munir l&rsquo;interpr\u00e9teur contre une division par z\u00e9ro non d\u00e9finie).<\/li>\n<li>Calcul simultan\u00e9 du num\u00e9rateur d&rsquo;\u00e9cart et du d\u00e9nominateur d&rsquo;\u00e9chelle au moyen d&rsquo;op\u00e9rations vectorielles hautement parall\u00e9lisables sur microprocesseur (op\u00e9rations SIMD) : soustraction \u00e9l\u00e9ment par \u00e9l\u00e9ment, calcul des valeurs absolues, et sommation par r\u00e9duction rapide.<\/li>\n<li>Division scalaire finale du num\u00e9rateur par le d\u00e9nominateur pour d\u00e9duire la valeur exacte de la semi-m\u00e9trique.<\/li>\n<\/ul>\n<p>L&rsquo;\u00e9limination formelle des boucles it\u00e9ratives au profit d&rsquo;op\u00e9rations d&rsquo;alg\u00e8bre matricielle (broadcasting en Python ou vectorisation matricielle sous R) permet d&rsquo;acc\u00e9l\u00e9rer le calcul des matrices de dissimilarit\u00e9 d&rsquo;un facteur de 50 \u00e0 100 fois par rapport \u00e0 une double boucle imbriqu\u00e9e conventionnelle, garantissant des performances adapt\u00e9es aux exigences actuelles du traitement des donn\u00e9es massives.<\/p>\n<h2>9. Applications en \u00e9cologie quantitative et dynamique des communaut\u00e9s<\/h2>\n<h3>9.1 Suivi des successions \u00e9cologiques temporelles<\/h3>\n<p>L&rsquo;une des contributions les plus spectaculaires de la m\u00e9trique de Bray-Curtis \u00e0 l&rsquo;\u00e9cologie fondamentale r\u00e9side dans l&rsquo;analyse quantitative des successions \u00e9cologiques temporelles et des processus de restauration environnementale. Lorsqu&rsquo;un \u00e9cosyst\u00e8me subit une perturbation majeure (incendie de for\u00eat, crue d\u00e9vastatrice, abandon d&rsquo;une parcelle agricole ou pollution industrielle des sols), la communaut\u00e9 biologique r\u00e9sidente s&rsquo;engage dans une trajectoire de recolonisation dynamique marqu\u00e9e par le renouvellement progressif des esp\u00e8ces au fil des ans.<\/p>\n<p>En \u00e9chantillonnant r\u00e9guli\u00e8rement un m\u00eame biotope sur une chronos\u00e9quence pluriannuelle, les \u00e9cologues calculent la dissimilarit\u00e9 de Bray-Curtis entre la composition observ\u00e9e \u00e0 l&rsquo;instant initial (<span class=\"arabpsy-math-inline\" data-latex=\"t_0\">$t_0$<\/span>) et les \u00e9tats successifs franchis par la communaut\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"t_1, t_2, dots, t_n\">$t_1, t_2, dots, t_n$<\/span>). Cette d\u00e9marche permet de tracer des <strong>vecteurs de trajectoire \u00e9cologique<\/strong> dans l&rsquo;espace multidimensionnel. L&rsquo;allure de la courbe de dissimilarit\u00e9 renseigne directement sur la vitesse et la dynamique du retour vers l&rsquo;\u00e9tat de r\u00e9f\u00e9rence : une d\u00e9croissance progressive et monotone de la dissimilarit\u00e9 t\u00e9moigne d&rsquo;une r\u00e9silience \u00e9cologique r\u00e9ussie, tandis qu&rsquo;une stagnation du coefficient \u00e0 des valeurs tr\u00e8s \u00e9lev\u00e9es (<span class=\"arabpsy-math-inline\" data-latex=\"BC &gt; 0{,}80\">$BC &gt; 0{,}80$<\/span>) signale l&rsquo;enfermement de l&rsquo;\u00e9cosyst\u00e8me dans un \u00e9tat alternatif stable d\u00e9grad\u00e9.<\/p>\n<p>De surcro\u00eet, la m\u00e9trique de Bray-Curtis offre une sensibilit\u00e9 remarquable pour la d\u00e9tection pr\u00e9coce des points de basculement \u00e9cologique (tipping points). Avant qu&rsquo;un \u00e9cosyst\u00e8me ne s&rsquo;effondre brutalement sous l&rsquo;effet du changement climatique ou de la surfertilisation anthropique, les abondances relatives des populations p\u00e9riph\u00e9riques se mettent \u00e0 fluctuer de mani\u00e8re asynchrone, induisant une instabilit\u00e9 de la dissimilarit\u00e9 de Bray-Curtis d&rsquo;une ann\u00e9e sur l&rsquo;autre, qui op\u00e8re ainsi comme un bio-indicateur d&rsquo;alerte pr\u00e9coce hautement diagnostique.<\/p>\n<h3>9.2 Analyse de gradients biog\u00e9ographiques et environnementaux<\/h3>\n<p>L&rsquo;analyse des patrons spatiaux de la biodiversit\u00e9 repose au niveau fondamental sur le concept de <strong>diversit\u00e9 b\u00eata<\/strong>, d\u00e9fini par Robert Whittaker comme le taux de renouvellement de la composition en esp\u00e8ces d&rsquo;un habitat \u00e0 un autre le long d&rsquo;un gradient g\u00e9ographique ou environnemental continu. La dissimilarit\u00e9 de Bray-Curtis s&rsquo;affirme comme le descripteur op\u00e9rationnel d&rsquo;excellence pour quantifier cette diversit\u00e9 b\u00eata dans l&rsquo;espace physique.<\/p>\n<p>Lors d&rsquo;\u00e9tudes men\u00e9es le long de gradients altitudinaux en montagne ou de transects d&rsquo;aridit\u00e9 en milieu sah\u00e9lien, la matrice de dissimilarit\u00e9 de Bray-Curtis entre paires de stations est confront\u00e9e \u00e0 la matrice des distances physiques (en kilom\u00e8tres) et \u00e0 la matrice des \u00e9carts environnementaux (diff\u00e9rences de temp\u00e9rature moyenne, d&rsquo;humidit\u00e9 \u00e9daphique ou de concentration en mati\u00e8re organique). Par le truchement de tests statistiques d&rsquo;association matricielle non param\u00e9triques, tels que le <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.1469-1809.1967.tb01267.x\">test de Mantel<\/a>, les \u00e9cologues quantifient l&rsquo;impact respectif de l&rsquo;isolement par la distance (dispersion limit\u00e9e) et du filtrage environnemental sur la structuration des communaut\u00e9s.<\/p>\n<p>En outre, l&rsquo;int\u00e9gration de Bray-Curtis dans des proc\u00e9dures avanc\u00e9es de partitionnement de la variance multivari\u00e9e permet de d\u00e9composer pr\u00e9cis\u00e9ment la variabilit\u00e9 spatiale totale observ\u00e9e : quelle proportion de la divergence compositionnelle d\u00e9coule directement des facteurs climatiques abiotiques, quelle part revient \u00e0 la configuration spatiale autocorr\u00e9l\u00e9e du paysage, et quelle fraction r\u00e9siduelle demeure attribuable aux processus stochastiques d&rsquo;\u00e9chantillonnage ? Cette puissance explicative conf\u00e8re \u00e0 la m\u00e9trique un r\u00f4le moteur dans la conception contemporaine des r\u00e9seaux d&rsquo;aires naturelles prot\u00e9g\u00e9es.<\/p>\n<h3>9.3 \u00c9tudes de m\u00e9tag\u00e9nomique et du microbiote<\/h3>\n<p>L&rsquo;av\u00e8nement des technologies de s\u00e9quen\u00e7age de nouvelle g\u00e9n\u00e9ration (NGS) a propuls\u00e9 la dissimilarit\u00e9 de Bray-Curtis au premier rang de la recherche biom\u00e9dicale via l&rsquo;analyse du microbiome humain et environnemental. Dans une exp\u00e9rience typique de m\u00e9tag\u00e9nomique cibl\u00e9e (par exemple le s\u00e9quen\u00e7age de la r\u00e9gion hypervariable V4 du g\u00e8ne de l&rsquo;ARN ribosomique 16S), les lectures brutes sont regroup\u00e9es en unit\u00e9s taxonomiques op\u00e9rationnelles (OTU) ou en variants de s\u00e9quences d&rsquo;amplicons (ASV). Les jeux de donn\u00e9es qui en r\u00e9sultent se mat\u00e9rialisent sous la forme de tables d&rsquo;abondance massives contenant des millions de comptages r\u00e9partis sur des centaines d&rsquo;\u00e9chantillons cliniques.<\/p>\n<p>Dans ce contexte, la matrice de Bray-Curtis sert \u00e0 mesurer la distance taxonomique globale s\u00e9parant les profils microbiens de cohortes de sujets sains et de patients atteints de pathologies chroniques telles que les maladies inflammatoires chroniques de l&rsquo;intestin (MICI), le diab\u00e8te de type 2 ou l&rsquo;ob\u00e9sit\u00e9 morbide. Le calcul de la dissimilarit\u00e9 de Bray-Curtis permet de formaliser math\u00e9matiquement le concept clinique de <strong>dysbiose intestinale<\/strong> : les profils microbiens des patients malades s&rsquo;\u00e9loignent significativement du nuage de points resserr\u00e9 form\u00e9 par les t\u00e9moins sains, traduisant une rupture profonde de l&rsquo;hom\u00e9ostasie bact\u00e9rienne symbiotique.<\/p>\n<p>De plus, l&rsquo;utilisation conjointe de la dissimilarit\u00e9 de Bray-Curtis avec des variables physiologiques ou s\u00e9riques (taux de cytokines pro-inflammatoires, perm\u00e9abilit\u00e9 intestinale, indices de r\u00e9sistance \u00e0 l&rsquo;insuline) permet d&rsquo;\u00e9tablir des corr\u00e9lations robustes reliant des modifications structurelles globales de la flore microbienne \u00e0 l&rsquo;\u00e9volution des marqueurs syst\u00e9miques de l&rsquo;h\u00f4te, ouvrant la voie \u00e0 des approches novatrices de m\u00e9decine translationnelle personnalis\u00e9e.<\/p>\n<h2>10. Transposition et applications en psychom\u00e9trie et sciences du comportement<\/h2>\n<h3>10.1 Mesure de la dissimilarit\u00e9 inter-individuelle des profils symptomatiques<\/h3>\n<p>La nosologie psychiatrique contemporaine fait face \u00e0 une crise \u00e9pist\u00e9mologique majeure quant \u00e0 la validit\u00e9 des cat\u00e9gories diagnostiques discr\u00e8tes formalis\u00e9es dans les manuels de r\u00e9f\u00e9rence tels que le DSM-5. L&rsquo;extr\u00eame h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 clinique constat\u00e9e au sein d&rsquo;une m\u00eame \u00e9tiquette diagnostique \u2014 deux patients diagnostiqu\u00e9s avec un \u00e9pisode d\u00e9pressif caract\u00e9ris\u00e9 pouvant ne partager qu&rsquo;un seul sympt\u00f4me marginal en commun \u2014 incite les chercheurs \u00e0 privil\u00e9gier une mod\u00e9lisation purement quantitative et dimensionnelle de la souffrance psychique.<\/p>\n<p>Dans cette perspective dimensionnelle novatrice, un individu est th\u00e9oriquement assimil\u00e9 \u00e0 un \u00e9cosyst\u00e8me biologique complexe, et chaque sympt\u00f4me \u00e9valu\u00e9 correspond \u00e0 une entit\u00e9 pond\u00e9r\u00e9e dont la s\u00e9v\u00e9rit\u00e9 refl\u00e8te l&rsquo;abondance. En appliquant la dissimilarit\u00e9 de Bray-Curtis \u00e0 des batteries cliniques valid\u00e9es comprenant de multiples items continus \u2014 telles que l&rsquo;\u00e9chelle des syndromes positif et n\u00e9gatif de la schizophr\u00e9nie (PANSS) ou l&rsquo;inventaire de d\u00e9pression de Beck (BDI) \u2014 les psychopathologues construisent des espaces g\u00e9om\u00e9triques d&rsquo;association inter-individuelle d&rsquo;une fid\u00e9lit\u00e9 sans pr\u00e9c\u00e9dent.<\/p>\n<p>Contrairement aux approches par score total cumul\u00e9 (qui masquent des configurations cliniques radicalement divergentes sous un m\u00eame chiffre r\u00e9capitulatif arbitraire), la matrice de Bray-Curtis capture l&rsquo;architecture fine des profils de s\u00e9v\u00e9rit\u00e9. L&rsquo;analyse typologique issue de cette m\u00e9trique permet d&rsquo;isoler des <strong>sous-groupes cliniques homog\u00e8nes<\/strong> (ou endoph\u00e9notypes comportementaux), facilitant ainsi le rep\u00e9rage de trajectoires \u00e9volutives diff\u00e9rentielles et l&rsquo;identification de profils de vuln\u00e9rabilit\u00e9 sp\u00e9cifiques face aux interventions pharmacologiques ou psychoth\u00e9rapeutiques.<\/p>\n<h3>10.2 Mod\u00e9lisation en r\u00e9seaux psychologiques et \u00e9cologie comportementale<\/h3>\n<p>La th\u00e9orie des r\u00e9seaux psychologiques, formalis\u00e9e notamment par Denny Borsboom, postule que les troubles mentaux ne r\u00e9sultent pas d&rsquo;une entit\u00e9 causale latente sous-jacente (comme un virus ou une l\u00e9sion c\u00e9r\u00e9brale unique), mais \u00e9mergent des interactions dynamiques directes et auto-entretenues s&rsquo;exer\u00e7ant entre les sympt\u00f4mes eux-m\u00eames (par exemple : l&rsquo;insomnie engendre la fatigue, qui pr\u00e9cipite la perte de concentration, laquelle induit l&rsquo;anxi\u00e9t\u00e9 de performance, renfor\u00e7ant en retour l&rsquo;insomnie).<\/p>\n<p>L&rsquo;int\u00e9gration de la m\u00e9thode d&rsquo;\u00e9chantillonnage de l&rsquo;exp\u00e9rience (Experience Sampling Method ou ESM) \u2014 consistant \u00e0 administrer aux sujets des micro-questionnaires sur smartphone plusieurs fois par jour au cours de leur vie quotidienne \u2014 g\u00e9n\u00e8re des s\u00e9ries temporelles multivari\u00e9es \u00e0 haute r\u00e9solution temporelle. La dissimilarit\u00e9 de Bray-Curtis s&rsquo;affirme ici comme un op\u00e9rateur privil\u00e9gi\u00e9 pour \u00e9valuer la distance entre les patrons dynamiques d&rsquo;activation comportementale et affective de diff\u00e9rents sujets ou d&rsquo;un m\u00eame sujet \u00e0 travers diff\u00e9rentes phases de vie.<\/p>\n<p>En calculant le coefficient de Bray-Curtis entre les r\u00e9pertoires d&rsquo;affects et d&rsquo;activit\u00e9s quotidiens, les chercheurs parviennent \u00e0 quantifier le degr\u00e9 d&rsquo;inflexibilit\u00e9 psychologique ou de rigidit\u00e9 comportementale des sujets. Deux individus affichant les m\u00eames scores d&rsquo;anxi\u00e9t\u00e9 moyenne peuvent pr\u00e9senter des dynamiques de r\u00e9activit\u00e9 environnementale totalement disjointes : l&rsquo;un conservant une dispersion harmonieuse de ses temps d&rsquo;activit\u00e9s, l&rsquo;autre voyant son r\u00e9pertoire comportemental se r\u00e9tracter massivement autour de routines d&rsquo;\u00e9vitement, une divergence que la m\u00e9trique de Bray-Curtis isole avec une remarquable acuit\u00e9.<\/p>\n<h3>10.3 Analyse de r\u00e9ponses aux tests de personnalit\u00e9 et questionnaires d&rsquo;attitudes<\/h3>\n<p>L&rsquo;\u00e9valuation psychom\u00e9trique de la personnalit\u00e9 normale et pathologique repose usuellement sur des inventaires auto-rapport\u00e9s \u00e0 choix multiples explorant un large \u00e9ventail de facettes psychologiques, \u00e0 l&rsquo;instar du mod\u00e8le en cinq facteurs (Big Five Inventory \/ NEO-PI-R) ou de l&rsquo;inventaire multiphasique de personnalit\u00e9 du Minnesota (MMPI). Dans ce cadre analytique, la comparaison des profils individuels au moyen de m\u00e9triques de distance traditionnelles se heurte fr\u00e9quemment \u00e0 des biais statistiques majeurs.<\/p>\n<p>Le recours \u00e0 la dissimilarit\u00e9 de Bray-Curtis \u00e9limine de fa\u00e7on d\u00e9cisive l&rsquo;art\u00e9fact d&rsquo;absence commune de traits pathologiques. Lors de l&rsquo;administration du MMPI \u00e0 une population g\u00e9n\u00e9rale non clinique, la vaste majorit\u00e9 des items mesurant des traits s\u00e9v\u00e8res d&rsquo;ali\u00e9nation mentale, de parano\u00efa d\u00e9lirante ou de d\u00e9viance psychopathique sont cot\u00e9s \u00e0 z\u00e9ro par l&rsquo;\u00e9crasante majorit\u00e9 des r\u00e9pondants. Une distance euclidienne conventionnelle tendrait \u00e0 d\u00e9cr\u00e9ter que deux individus ont des personnalit\u00e9s quasi identiques simplement parce qu&rsquo;aucun des deux ne manifeste de d\u00e9lires hallucinatoires. Bray-Curtis concentre au contraire l&rsquo;analyse de ressemblance sur les dimensions de personnalit\u00e9 positivement exprim\u00e9es (extraversion, ouverture, amabilit\u00e9, conscience).<\/p>\n<p>De surcro\u00eet, la m\u00e9trique se r\u00e9v\u00e8le particuli\u00e8rement pr\u00e9cieuse dans les protocoles longitudinaux d&rsquo;\u00e9valuation de l&rsquo;efficacit\u00e9 psychoth\u00e9rapeutique. En mesurant la dissimilarit\u00e9 de Bray-Curtis entre le profil d&rsquo;attitudes d&rsquo;un patient \u00e0 l&rsquo;admission en th\u00e9rapie et ses \u00e9valuations successives au fil des mois, les cliniciens disposent d&rsquo;un indice scalaire objectif documentant la restructuration globale de sa personnalit\u00e9 et l&rsquo;att\u00e9nuation graduelle de ses m\u00e9canismes de d\u00e9fense dysfonctionnels.<\/p>\n<h2>11. Visualisation multivari\u00e9e : Int\u00e9gration dans l&rsquo;ordonnancement NMDS<\/h2>\n<h3>11.1 Le positionnement multidimensionnel non m\u00e9trique (NMDS)<\/h3>\n<p>Comme d\u00e9montr\u00e9 formellement dans la section consacr\u00e9e aux propri\u00e9t\u00e9s m\u00e9triques, la dissimilarit\u00e9 de Bray-Curtis constitue une semi-m\u00e9trique qui viole occasionnellement l&rsquo;in\u00e9galit\u00e9 triangulaire, rendant son espace topologique incompatible avec les contraintes d&rsquo;une projection lin\u00e9aire directe en g\u00e9om\u00e9trie euclidienne pure. Tenter d&rsquo;analyser une matrice de Bray-Curtis au moyen d&rsquo;une analyse en composantes principales (ACP) conventionnelle rel\u00e8ve de l&rsquo;h\u00e9r\u00e9sie math\u00e9matique. La technique d&rsquo;ordination reine, universellement reconnue comme la compagne m\u00e9thodologique naturelle de Bray-Curtis, est le <strong>positionnement multidimensionnel non m\u00e9trique<\/strong> (Non-Metric Multidimensional Scaling ou NMDS), conceptualis\u00e9 par Shepard et formalis\u00e9 par <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/BF02289565\">Kruskal (1964)<\/a>.<\/p>\n<p>Le principe algorithmique de la NMDS se distingue radicalement des m\u00e9thodes d&rsquo;ordination factorielle classiques : au lieu de chercher \u00e0 pr\u00e9server fid\u00e8lement les grandeurs absolues des distances g\u00e9om\u00e9triques, la NMDS s&rsquo;assigne pour objectif exclusif de pr\u00e9server le <strong>rang relatif<\/strong> des dissimilarit\u00e9s au sein d&rsquo;un espace g\u00e9om\u00e9trique de dimension r\u00e9duite arbitrairement fix\u00e9 \u00e0 l&rsquo;avance (g\u00e9n\u00e9ralement <span class=\"arabpsy-math-inline\" data-latex=\"k = 2\">$k = 2$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"k = 3\">$k = 3$<\/span> pour des raisons \u00e9videntes de visualisation graphique).<\/p>\n<p>L&rsquo;ajustement optimal de la configuration spatiale des points dans le rep\u00e8re bidimensionnel est guid\u00e9 par une fonction de perte it\u00e9rative non lin\u00e9aire, appel\u00e9e le <strong>stress de Kruskal<\/strong> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Stress} = \\sqrt{\\frac{\\sum_{i &lt; j} (d_{ij} - \\hat{d}_{ij})^2}{\\sum_{i &lt; j} d_{ij}^2}}\">$$\\text{Stress} = \\sqrt{\\frac{\\sum_{i &lt; j} (d_{ij} &#8211; \\hat{d}_{ij})^2}{\\sum_{i &lt; j} d_{ij}^2}}$$<\/span><\/p>\n<p>Dans cette formulation, <span class=\"arabpsy-math-inline\" data-latex=\"d_{ij}\">$d_{ij}$<\/span> repr\u00e9sente la distance euclidienne mesur\u00e9e entre les points <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span> sur le plan graphique d&rsquo;ordination, tandis que <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{d}_{ij}\">$\\hat{d}_{ij}$<\/span> repr\u00e9sente la distance pr\u00e9dite d\u00e9duite d&rsquo;une r\u00e9gression monotone (isotonique) ajust\u00e9e sur les rangs des dissimilarit\u00e9s de Bray-Curtis r\u00e9elles d&rsquo;origine (<span class=\"arabpsy-math-inline\" data-latex=\"BC_{ij}\">$BC_{ij}$<\/span>). L&rsquo;\u00e9valuation de la qualit\u00e9 de la repr\u00e9sentation graphique s&rsquo;effectue via l&rsquo;examen du <strong>diagramme de Shepard<\/strong> (qui confronte visuellement les distances d&rsquo;ordination aux dissimilarit\u00e9s observ\u00e9es) et l&rsquo;interpr\u00e9tation de la valeur finale du stress :<\/p>\n<ul>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"\\text{Stress} &lt; 0{,}05\">$\\text{Stress} &lt; 0{,}05$<\/span> : Repr\u00e9sentation d&#039;une pr\u00e9cision exceptionnelle, exempte de tout risque de mauvaise interpr\u00e9tation graphique.<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"0{,}05 le \\text{Stress} &lt; 0{,}10\">$0{,}05 le \\text{Stress} &lt; 0{,}10$<\/span> : Excellent ajustement, configuration spatiale pleinement repr\u00e9sentative des structures \u00e9cologiques ou cliniques sous-jacentes.<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"0{,}10 le \\text{Stress} &lt; 0{,}20\">$0{,}10 le \\text{Stress} &lt; 0{,}20$<\/span> : Repr\u00e9sentation globale exploitable et informative, bien que les d\u00e9tails des proximit\u00e9s locales entre points voisins doivent \u00eatre interpr\u00e9t\u00e9s avec circonspection.<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"\\text{Stress} ge 0{,}20\">$\\text{Stress} ge 0{,}20$<\/span> : Ajustement m\u00e9diocre \u00e0 trompeur ; il devient imp\u00e9ratif d&rsquo;augmenter le nombre de dimensions de projection (<span class=\"arabpsy-math-inline\" data-latex=\"k = 3\">$k = 3$<\/span>) pour \u00e9viter des distorsions visuelles majeures.<\/li>\n<\/ul>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/bray_curtis3.png\" alt=\"\" \/><\/figure>\n<h3>11.2 Classification ascendante hi\u00e9rarchique (CAH) bas\u00e9e sur Bray-Curtis<\/h3>\n<p>Outre l&rsquo;ordination par NMDS, les matrices de dissimilarit\u00e9 de Bray-Curtis constituent le mat\u00e9riau d&rsquo;entr\u00e9e privil\u00e9gi\u00e9 pour les algorithmes de <strong>classification ascendante hi\u00e9rarchique<\/strong> (CAH), visant \u00e0 structurer les unit\u00e9s d&rsquo;\u00e9chantillonnage en clusters discrets embo\u00eet\u00e9s sous forme d&rsquo;arbres phylog\u00e9n\u00e9tiques ou \u00e9cologiques appel\u00e9s dendrogrammes.<\/p>\n<p>Dans ce cadre classificatoire, le choix du crit\u00e8re d&rsquo;agr\u00e9gation liant les groupes de profils rev\u00eat une importance cruciale :<\/p>\n<ul>\n<li>Le crit\u00e8re de la <strong>liaison moyenne<\/strong> (UPGMA &#8211; Unweighted Pair Group Method with Arithmetic Mean) est incontestablement le plus respectueux de la topologie de Bray-Curtis. \u00c0 chaque \u00e9tape fusionnelle, la dissimilarit\u00e9 entre deux clusters <span class=\"arabpsy-math-inline\" data-latex=\"A\">$A$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> est d\u00e9finie comme la moyenne arithm\u00e9tique exacte des coefficients de Bray-Curtis s\u00e9parant l&rsquo;ensemble des paires d&rsquo;\u00e9chantillons intra-groupes. Ce crit\u00e8re maximise le coefficient de corr\u00e9lation ph\u00e9n\u00e9tique entre la matrice de distance originale et la matrice d\u00e9riv\u00e9e du dendrogramme.<\/li>\n<li>La m\u00e9thode de <strong>Ward<\/strong>, bien qu&rsquo;extr\u00eamement populaire pour sa tendance \u00e0 former des groupes compacts et d&rsquo;effectifs \u00e9quilibr\u00e9s, repose originellement sur le calcul de la somme des carr\u00e9s des \u00e9carts \u00e0 la moyenne dans un espace strictement euclidien. L&rsquo;application directe de l&rsquo;algorithme classique de Ward \u00e0 une matrice de Bray-Curtis non corrig\u00e9e est th\u00e9oriquement contestable en raison de la nature semi-m\u00e9trique du coefficient, et expose au risque d&rsquo;inversions dendrogrammiques. Pour employer l\u00e9gitimement un crit\u00e8re de type Ward sur des dissimilarit\u00e9s de Bray-Curtis, il convient d&rsquo;appliquer au pr\u00e9alable une transformation par racine carr\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"\\sqrt{BC_{ij}}\">$\\sqrt{BC_{ij}}$<\/span>) ou d&rsquo;ex\u00e9cuter une correction par coordonn\u00e9es principales de Cailliez ou de Lingoes garantissant l&rsquo;euclidianit\u00e9 du tableau d&rsquo;entr\u00e9e.<\/li>\n<\/ul>\n<p>Pour d\u00e9terminer de mani\u00e8re objective le nombre optimal de clusters \u00e0 partitionner le long du dendrogramme, les praticiens couplent usuellement la CAH \u00e0 l&rsquo;analyse des largeurs de silhouette de Rousseeuw ou \u00e0 l&rsquo;indice de Calinski-Harabasz adapt\u00e9 aux dissimilarit\u00e9s arbitraires.<\/p>\n<h3>11.3 Interpr\u00e9tation graphique et inf\u00e9rence statistique multivari\u00e9e<\/h3>\n<p>L&rsquo;exploration visuelle d&rsquo;un plan d&rsquo;ordination NMDS issu d&rsquo;une matrice de Bray-Curtis gagne immens\u00e9ment en puissance lorsqu&rsquo;elle est combin\u00e9e \u00e0 des techniques d&rsquo;analyse explicative et d&rsquo;inf\u00e9rence statistique formelle. Dans le cadre de l&rsquo;analyse multivari\u00e9e moderne, la simple inspection qualitative d&rsquo;un nuage de points est syst\u00e9matiquement corrobor\u00e9e par des tests d&rsquo;hypoth\u00e8ses rigoureux.<\/p>\n<p>La premi\u00e8re technique d&rsquo;enrichissement graphique consiste \u00e0 projeter des vecteurs environnementaux ou cliniques sur le plan d&rsquo;ordination (proc\u00e9dure <code>envfit<\/code> sous R). Les gradients continus (pH du sol, temp\u00e9rature, score d&rsquo;\u00e2ge ou taux sanguin de cortisol) sont projet\u00e9s sous forme de fl\u00e8ches orient\u00e9es dont la direction indique la pente d&rsquo;accroissement maximal de la variable et dont la longueur est strictement proportionnelle \u00e0 la force de corr\u00e9lation (<span class=\"arabpsy-math-inline\" data-latex=\"r^2\">$r^2$<\/span>) \u00e9tablie entre le gradient et la configuration d&rsquo;ordination globale.<\/p>\n<p>Sur le plan de l&rsquo;inf\u00e9rence confirmatoire, la m\u00e9thode de r\u00e9f\u00e9rence absolue pour tester des hypoth\u00e8ses de diff\u00e9renciation de groupes structur\u00e9s a priori est l&rsquo;<strong>analyse de variance multivari\u00e9e par permutations<\/strong> (PERMANOVA), formalis\u00e9e par <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.1442-9993.2001.01070.pp.x\">Marti J. Anderson (2001)<\/a> et impl\u00e9ment\u00e9e via la fonction <code>adonis2<\/code>. Contrairement \u00e0 la MANOVA param\u00e9trique classique qui exige la stricte multinormalit\u00e9 des r\u00e9sidus et l&rsquo;\u00e9galit\u00e9 des matrices de variance-covariance, la PERMANOVA op\u00e8re directement sur la matrice de dissimilarit\u00e9 de Bray-Curtis :<\/p>\n<p>Elle calcule une statistique pseudo-<span class=\"arabpsy-math-inline\" data-latex=\"F\">$F$<\/span> quantifiant le rapport de dispersion entre les diff\u00e9rents groupes exp\u00e9rimentaux par rapport \u00e0 la dispersion intra-groupe, puis \u00e9value la significativit\u00e9 statistique (<span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>-value) par le biais de milliers de permutations al\u00e9atoires des \u00e9tiquettes des \u00e9chantillons. Ce test non param\u00e9trique permet de statuer sans \u00e9quivoque : l&rsquo;appartenance \u00e0 un habitat distinct ou \u00e0 une sous-cat\u00e9gorie pathologique engendre-t-elle une s\u00e9paration statistiquement significative des structures de composition biologique ou clinique ?<\/p>\n<p>Enfin, pour garantir la validit\u00e9 d&rsquo;une PERMANOVA, il est indispensable de v\u00e9rifier que la dispersion multivari\u00e9e intra-groupe (l&rsquo;\u00e9tendue du nuage de points au sein de chaque modalit\u00e9) est relativement homog\u00e8ne, une \u00e9tape r\u00e9alis\u00e9e par le test non param\u00e9trique de dispersion multivari\u00e9e <strong>PERMDISP<\/strong> (analogue multivari\u00e9 direct du test d&rsquo;homosc\u00e9dasticit\u00e9 de Levene). L&rsquo;articulation coordonn\u00e9e de Bray-Curtis, de la NMDS, de la PERMANOVA et de PERMDISP forme aujourd&rsquo;hui le socle m\u00e9thodologique standard de l&rsquo;analyse multivari\u00e9e non param\u00e9trique dans les sciences de la vie et du comportement.<\/p>\n<h2>12. Limites m\u00e9thodologiques, pi\u00e8ges fr\u00e9quents et recommandations pratiques<\/h2>\n<h3>12.1 Le pi\u00e8ge de la d\u00e9pendance \u00e0 l&rsquo;\u00e9chelle des unit\u00e9s de mesure<\/h3>\n<p>En d\u00e9pit de son \u00e9l\u00e9gance math\u00e9matique et de son universalit\u00e9 op\u00e9rationnelle, la dissimilarit\u00e9 de Bray-Curtis rec\u00e8le plusieurs chausse-trappes m\u00e9thodologiques qu&rsquo;un analyste rigoureux doit imp\u00e9rativement ma\u00eetriser afin d&rsquo;\u00e9viter des erreurs d&rsquo;interpr\u00e9tation critiques. Le premier pi\u00e8ge r\u00e9side dans l&rsquo;obligation formelle d&rsquo;<strong>homog\u00e9n\u00e9it\u00e9 dimensionnelle absolue<\/strong> entre l&rsquo;ensemble des descripteurs composant la matrice d&rsquo;observation.<\/p>\n<p>La formulation alg\u00e9brique de Bray-Curtis additionne les valeurs brutes de diff\u00e9rentes colonnes au sein d&rsquo;une m\u00eame somme marginale (<span class=\"arabpsy-math-inline\" data-latex=\"S_i = \\sum x_{ik}\">$S_i = \\sum x_{ik}$<\/span>). Cette op\u00e9ration de sommation n&rsquo;a de sens physique et statistique que si et seulement si toutes les variables <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> sont mesur\u00e9es rigoureusement dans la m\u00eame unit\u00e9 scalaire et partagent le m\u00eame ordre de grandeur fondamental. Si une matrice agr\u00e8ge par inadvertance des comptages d&rsquo;individus entiers (variant de 0 \u00e0 500 individus), des mesures de biomasse exprim\u00e9es en milligrammes (variant de 0 \u00e0 10 000 mg) et des pourcentages de recouvrement foliaire (born\u00e9s de 0 \u00e0 100 %), la m\u00e9trique de Bray-Curtis perd toute coh\u00e9rence math\u00e9matique : elle additionne litt\u00e9ralement des grandeurs incommensurables.<\/p>\n<p>Face \u00e0 des jeux de donn\u00e9es composites regroupant des variables environnementales ou psychologiques h\u00e9t\u00e9rog\u00e8nes (par exemple le revenu annuel en euros, le niveau d&rsquo;anxi\u00e9t\u00e9 sur \u00e9chelle de Likert de 1 \u00e0 7 et la dur\u00e9e de sommeil en minutes), l&rsquo;utilisation directe de Bray-Curtis doit \u00eatre formellement proscrite. Dans ce cas de figure pr\u00e9cis, il convient de s&rsquo;orienter vers des coefficients sp\u00e9cifiquement con\u00e7us pour l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 dimensionnelle, au premier rang desquels figure le <strong>coefficient de dissimilarit\u00e9 de Gower<\/strong>, qui standardise de mani\u00e8re autonome chaque descripteur par son amplitude d&rsquo;\u00e9chantillonnage avant d&rsquo;en calculer la moyenne pond\u00e9r\u00e9e.<\/p>\n<h3>12.2 Probl\u00e8mes de saturation pour les sites ou sujets totalement distincts<\/h3>\n<p>Une seconde limite structurelle inh\u00e9rente \u00e0 la dissimilarit\u00e9 de Bray-Curtis est le ph\u00e9nom\u00e8ne d&rsquo;<strong>asymptote de saturation \u00e0 l&rsquo;unit\u00e9<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"BC = 1\">$BC = 1$<\/span>). D\u00e8s lors que deux unit\u00e9s d&rsquo;\u00e9chantillonnage ne partagent strictement aucune variable commune en raison d&rsquo;un remplacement complet des esp\u00e8ces ou des sympt\u00f4mes, la valeur de dissimilarit\u00e9 atteint son plafond absolu de <span class=\"arabpsy-math-inline\" data-latex=\"1{,}00\">$1{,}00$<\/span>. Or, sur le plan \u00e9cologique ou clinique, tous les couples d&rsquo;\u00e9chantillons disjoints ne sont pas n\u00e9cessairement s\u00e9par\u00e9s par la m\u00eame distance biologique ou fonctionnelle sous-jacente.<\/p>\n<p>Consid\u00e9rons trois sites dispos\u00e9s le long d&rsquo;un vaste gradient thermique continental : le site A h\u00e9berge une communaut\u00e9 de toundra alpine, le site B abrite une steppe temp\u00e9r\u00e9e sans aucune esp\u00e8ce commune avec A (<span class=\"arabpsy-math-inline\" data-latex=\"BC_{AB} = 1{,}00\">$BC_{AB} = 1{,}00$<\/span>), et le site C est une for\u00eat pluviale \u00e9quatoriale n&rsquo;ayant aucune esp\u00e8ce commune ni avec B (<span class=\"arabpsy-math-inline\" data-latex=\"BC_{BC} = 1{,}00\">$BC_{BC} = 1{,}00$<\/span>), ni avec A (<span class=\"arabpsy-math-inline\" data-latex=\"BC_{AC} = 1{,}00\">$BC_{AC} = 1{,}00$<\/span>). Dans l&rsquo;espace de Bray-Curtis brut, le site A est jug\u00e9 exactement \u00e0 la m\u00eame distance du site B que du site C (<span class=\"arabpsy-math-inline\" data-latex=\"1{,}00\">$1{,}00$<\/span>). L&rsquo;indice est incapable de retranscrire le fait que le site B est \u00e9cologiquement bien plus proche du site A que ne l&rsquo;est l&rsquo;\u00e9cosyst\u00e8me \u00e9quatorial C.<\/p>\n<p>Ce ph\u00e9nom\u00e8ne d&rsquo;\u00e9crasement ou d&rsquo;effet de troncature conduit \u00e0 des distorsions s\u00e9v\u00e8res lors des ordinations spatiales, les sites totalement disjoints \u00e9tant artificiellement repli\u00e9s les uns sur les autres sur le pourtour p\u00e9riph\u00e9rique des repr\u00e9sentations NMDS (artefact dit en fer \u00e0 cheval ou effet d&rsquo;arche). Pour pallier cette saturation lors de l&rsquo;\u00e9tude de gradients d&rsquo;une ampleur biog\u00e9ographique majeure, il est recommand\u00e9 de recourir \u00e0 l&rsquo;algorithme des <strong>dissimilarit\u00e9s \u00e9tendues<\/strong> (Extended Dissimilarities de De&rsquo;ath) ou \u00e0 la distance g\u00e9od\u00e9sique de plus court chemin calcul\u00e9e sur un graphe de voisinage \u00e9cologique (algorithme Isomap), qui r\u00e9tablissent la m\u00e9trique sous-jacente en cheminant \u00e0 travers une cha\u00eene d&rsquo;\u00e9chantillons interm\u00e9diaires connect\u00e9s.<\/p>\n<h3>12.3 Bonnes pratiques pour une utilisation rigoureuse en recherche<\/h3>\n<p>Pour conclure cette analyse approfondie, il est utile de synth\u00e9tiser les recommandations m\u00e9thodologiques incontournables sous la forme d&rsquo;un protocole d\u00e9cisionnel destin\u00e9 \u00e0 s\u00e9curiser l&#8217;emploi de la dissimilarit\u00e9 de Bray-Curtis dans les protocoles de recherche contemporains :<\/p>\n<ul>\n<li><strong>V\u00e9rification de la nature des variables :<\/strong> S&rsquo;assurer scrupuleusement que les donn\u00e9es d&rsquo;entr\u00e9e sont toutes non n\u00e9gatives (<span class=\"arabpsy-math-inline\" data-latex=\"x_{ik} ge 0\">$x_{ik} ge 0$<\/span>) et mesur\u00e9es rigoureusement dans la m\u00eame unit\u00e9 physique ou psychom\u00e9trique d&rsquo;intensit\u00e9.<\/li>\n<li><strong>Diagnostic de la distribution des abondances :<\/strong> R\u00e9aliser une analyse exploratoire de l&rsquo;asym\u00e9trie de distribution des variables. En pr\u00e9sence de descripteurs hyper-dominants exer\u00e7ant un effet d&rsquo;\u00e9crasement statistique, appliquer imp\u00e9rativement une transformation pr\u00e9alable par racine carr\u00e9e, racine quatri\u00e8me ou logarithme (<span class=\"arabpsy-math-inline\" data-latex=\"log(x+1)\">$log(x+1)$<\/span>).<\/li>\n<li><strong>Arbitrage sur l&rsquo;effort d&rsquo;\u00e9chantillonnage :<\/strong> D\u00e9terminer en amont si la taille totale absolue des collectes constitue une information \u00e9cologique signifiante ou un biais d&rsquo;\u00e9chantillonnage \u00e0 neutraliser. Dans le second cas, convertir les profils bruts en proportions relatives standardis\u00e9es (<span class=\"arabpsy-math-inline\" data-latex=\"\\sum x_{ik} = 1\">$\\sum x_{ik} = 1$<\/span>) avant le calcul de Bray-Curtis.<\/li>\n<li><strong>Choix r\u00e9fl\u00e9chi du mode d&rsquo;ordination :<\/strong> Bannir formellement les projections lin\u00e9aires euclidiennes directes (ACP) au profit de m\u00e9thodes d&rsquo;ordination non m\u00e9triques robustes bas\u00e9es sur les rangs (NMDS) ou d&rsquo;analyses en coordonn\u00e9es principales (PCoA) int\u00e9grant des corrections d&rsquo;euclidianit\u00e9 pour valeurs propres n\u00e9gatives.<\/li>\n<li><strong>Transparence computationnelle absolue :<\/strong> Documenter explicitement dans la section m\u00e9thodologique des publications le logiciel utilis\u00e9, le package employ\u00e9 (ex. vegan sous R, scipy sous Python), la fonction exacte appel\u00e9e ainsi que la cha\u00eene compl\u00e8te des transformations appliqu\u00e9es \u00e0 la matrice de donn\u00e9es brutes.<\/li>\n<\/ul>\n<p>En respectant scrupuleusement ces canons m\u00e9thodologiques, la dissimilarit\u00e9 de Bray-Curtis d\u00e9ploie toute sa puissance analytique, offrant aux chercheurs un instrument statistique d&rsquo;une sensibilit\u00e9 et d&rsquo;une pertinence th\u00e9orique incomparables pour d\u00e9m\u00ealer la complexit\u00e9 multivari\u00e9e du vivant et du comportement humain.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. <em>Austral Ecology<\/em>, 26(1), 32\u201346. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.1442-9993.2001.01070.pp.x\">https:\/\/doi.org\/10.1111\/j.1442-9993.2001.01070.pp.x<\/a><\/li>\n<li>Borsboom, D. (2017). A network approach to psychopathology. <em>World Psychiatry<\/em>, 16(1), 5\u201313. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/wps.20375\">https:\/\/doi.org\/10.1002\/wps.20375<\/a><\/li>\n<li>Bray, J. R., &amp; Curtis, J. T. (1957). An ordination of the upland forest communities of southern Wisconsin. <em>Ecological Monographs<\/em>, 27(4), 325\u2013349. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1942268\">https:\/\/doi.org\/10.2307\/1942268<\/a><\/li>\n<li>Clarke, K. R., &amp; Warwick, R. M. (2001). <em>Change in marine communities: An approach to statistical analysis and interpretation<\/em> (2nd ed.). PRIMER-E, Plymouth.<\/li>\n<li>De&rsquo;ath, G. (1999). Extended dissimilarity: A method of robust estimation of ecological distances from high beta diversity data. <em>Plant Ecology<\/em>, 144(2), 191\u2013199. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1023\/A:1009848526547\">https:\/\/doi.org\/10.1023\/A:1009848526547<\/a><\/li>\n<li>Gleason, H. A. (1926). The individualistic concept of the plant association. <em>Bulletin of the Torrey Botanical Club<\/em>, 53(1), 7\u201326. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2479933\">https:\/\/doi.org\/10.2307\/2479933<\/a><\/li>\n<li>Horn, H. S. (1966). Measurement of \u00ab\u00a0overlap\u00a0\u00bb in comparative ecological studies. <em>The American Naturalist<\/em>, 100(914), 419\u2013424. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1086\/282436\">https:\/\/doi.org\/10.1086\/282436<\/a><\/li>\n<li>Jaccard, P. (1901). \u00c9tude comparative de la distribution florale dans une portion des Alpes et du Jura. <em>Bulletin de la Soci\u00e9t\u00e9 Vaudoise des Sciences Naturelles<\/em>, 37, 547\u2013579.<\/li>\n<li>Kruskal, J. B. (1964). Nonmetric multidimensional scaling: A numerical method. <em>Psychometrika<\/em>, 29(2), 115\u2013129. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/BF02289565\">https:\/\/doi.org\/10.1007\/BF02289565<\/a><\/li>\n<li>Legendre, P., &amp; Gallagher, E. D. (2001). Ecologically meaningful transformations for ordination of species data. <em>Oecologia<\/em>, 129(2), 271\u2013280. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/s004420100716\">https:\/\/doi.org\/10.1007\/s004420100716<\/a><\/li>\n<li>Legendre, P., &amp; Legendre, L. (2012). <em>Numerical ecology<\/em> (3rd English ed.). Elsevier Science BV. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.elsevier.com\/books\/numerical-ecology\/legendre\/978-0-444-53868-0\">https:\/\/www.elsevier.com\/books\/numerical-ecology\/legendre\/978-0-444-53868-0<\/a><\/li>\n<li>Mantel, N. (1967). The detection of disease clustering and a generalized regression approach. <em>Cancer Research<\/em>, 27(2), 209\u2013220. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.1469-1809.1967.tb01267.x\">https:\/\/doi.org\/10.1111\/j.1469-1809.1967.tb01267.x<\/a><\/li>\n<li>Morisita, M. (1959). Measuring of interspecific association and similarity between communities. <em>Memoirs of the Faculty of Science, Kyushu University, Series E (Biology)<\/em>, 3(1), 65\u201380. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/BF02518851\">https:\/\/doi.org\/10.1007\/BF02518851<\/a><\/li>\n<li>Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., Minchin, P. R., O&rsquo;Hara, R. B., Simpson, G. L., Solymos, P., Stevens, M. H. H., Szoecs, E., &amp; Wagner, H. (2020). <em>vegan: Community ecology package<\/em> (R package version 2.5-7). <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=vegan\">https:\/\/cran.r-project.org\/package=vegan<\/a><\/li>\n<li>Ru\u017ei\u010dka, M. (1958). Anwendung mathematisch-statistischer Methoden in der Geobotanik (synthetische Bearbeitung von Aufnahmen). <em>Biol\u00f3gia (Bratislava)<\/em>, 13, 647\u2013661.<\/li>\n<li>S\u00f8rensen, T. (1948). A method of establishing groups of equal amplitude in plant sociology based on similarity of species content and its application to analyses of the vegetation on Danish commons. <em>Kongelige Danske Videnskabernes Selskab, Biologiske Skrifter<\/em>, 5(4), 1\u201334.<\/li>\n<li>Whittaker, R. H. (1960). Vegetation of the Siskiyou Mountains, Oregon and California. <em>Ecological Monographs<\/em>, 30(3), 279\u2013338. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/1943563\">https:\/\/doi.org\/10.2307\/1943563<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet sur la dissimilarit\u00e9 de Bray-Curtis : formulation math\u00e9matique, calculs pas \u00e0 pas, applications \u00e9cologiques et psychom\u00e9triques.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[1829,2620],"tags":[3791,11312,11314,11313,16],"class_list":["post-9892","psych_stat","type-psych_stat","status-publish","hentry","category-methodologie-statistique","category-psychometrie-et-analyse-de-donnees","tag-analyse-multivariee","tag-dissimilarite-de-bray-curtis","tag-ecologie-quantitative","tag-nmds","tag-psychometrie"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/9892","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=9892"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/9892\/revisions"}],"predecessor-version":[{"id":9896,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/9892\/revisions\/9896"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=9892"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=9892"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=9892"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}