{"id":5549,"date":"2026-09-11T22:45:40","date_gmt":"2026-09-11T22:45:40","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/classification-hierarchique-r-exemple-etape-par-etape\/"},"modified":"2026-09-11T22:45:43","modified_gmt":"2026-09-11T22:45:43","slug":"classification-hierarchique-r-exemple-etape-par-etape","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=classification-hierarchique-r-exemple-etape-par-etape","title":{"rendered":"Classification hi\u00e9rarchique en R : exemple \u00e9tape par \u00e9tape"},"content":{"rendered":"<p>L&rsquo;analyse des donn\u00e9es contemporaine, confront\u00e9e \u00e0 l&rsquo;explosion de la dimensionnalit\u00e9 et \u00e0 la complexit\u00e9 croissante des architectures d&rsquo;\u00e9chantillonnage, exige des m\u00e9thodologies d&rsquo;investigation capables d&rsquo;extraire du sens sans imposer de carcan th\u00e9orique restrictif. Dans les sciences comportementales, la psychom\u00e9trie, la sociologie quantitative et la bio-informatique, les chercheurs se heurtent fr\u00e9quemment \u00e0 des corpus empiriques au sein desquels aucune variable d&rsquo;\u00e9tiquetage n&rsquo;est pr\u00e9alablement d\u00e9finie. Face \u00e0 ces espaces multivari\u00e9s, la taxonomie computationnelle s&rsquo;\u00e9rige en discipline cardinale pour structurer l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 sous-jacente des observations et r\u00e9v\u00e9ler des groupements naturels jusqu&rsquo;alors occult\u00e9s par la masse informationnelle brute.<\/p>\n<p>Au c\u0153ur de cet arsenal exploratoire, la classification ascendante hi\u00e9rarchique (CAH) occupe une position \u00e9pist\u00e9mologique et pratique pr\u00e9pond\u00e9rante. Contrairement aux approches de partitionnement direct, la d\u00e9marche hi\u00e9rarchique ne force pas les observations dans un nombre de cat\u00e9gories arbitrairement fig\u00e9 d\u00e8s l&rsquo;amorce de l&rsquo;analyse. Elle d\u00e9ploie au contraire une perspective morphologique continue, organisant les entit\u00e9s \u00e9tudi\u00e9es selon une arborescence gradu\u00e9e qui capture les proximit\u00e9s proximales et les divergences globales. Cette structure g\u00e9n\u00e9alogique offre une finesse d&rsquo;analyse in\u00e9gal\u00e9e pour appr\u00e9hender les organisations multi-\u00e9chelles, des micropopulations homog\u00e8nes aux macro-typologies g\u00e9n\u00e9rales.<\/p>\n<p>L&rsquo;environnement statistique R s&rsquo;impose historiquement comme l&rsquo;\u00e9cosyst\u00e8me de pr\u00e9dilection pour conduire ces investigations typologiques de haute pr\u00e9cision. Dot\u00e9 d&rsquo;une impl\u00e9mentation matricielle native optimis\u00e9e et enrichi par une communaut\u00e9 prolifique d&rsquo;algorithmiciennes et de statisticiens, R propose une suite compl\u00e8te d&rsquo;outils d\u00e9di\u00e9s \u00e0 l&rsquo;ing\u00e9nierie des distances, \u00e0 l&rsquo;agr\u00e9gation hi\u00e9rarchique, \u00e0 la mod\u00e9lisation s\u00e9miologique des arbres et \u00e0 la validation inf\u00e9rentielle des topologies obtenues. Cet article propose une d\u00e9construction exhaustive, th\u00e9orique et op\u00e9rationnelle, de la classification hi\u00e9rarchique sous R, guidant le chercheur depuis les axiomes fondamentaux de la dissimilarit\u00e9 jusqu&rsquo;\u00e0 la validation robuste et l&rsquo;interpr\u00e9tation substantielle des profils comportementaux.<\/p>\n<h2>1. Fondements conceptuels de la classification non supervis\u00e9e et pertinence acad\u00e9mique<\/h2>\n<h3>1.1 D\u00e9finition et principes de l&rsquo;apprentissage non supervis\u00e9<\/h3>\n<p>L&rsquo;apprentissage non supervis\u00e9 se distingue fondamentalement du paradigme supervis\u00e9 par l&rsquo;absence d\u00e9lib\u00e9r\u00e9e d&rsquo;une variable cible, d&rsquo;un label externe ou d&rsquo;un crit\u00e8re de v\u00e9rit\u00e9 terrain pr\u00e9alablement assign\u00e9 aux unit\u00e9s statistiques. Dans le cadre supervis\u00e9, l&rsquo;objectif algorithmique consiste \u00e0 optimiser une fonction de perte param\u00e9trique ou non param\u00e9trique reliant un vecteur de covariables \u00e0 une r\u00e9ponse empirique connue, qu&rsquo;il s&rsquo;agisse d&rsquo;une r\u00e9gression continue ou d&rsquo;une discrimination cat\u00e9gorielle. En revanche, la classification non supervis\u00e9e \u00e9volue au sein d&rsquo;un espace multidimensionnel o\u00f9 seules les relations endog\u00e8nes entre observations guident le processus de d\u00e9couverte statistique. Il s&rsquo;agit d&rsquo;une d\u00e9marche purement inductive dans laquelle l&rsquo;algorithme doit synth\u00e9tiser la topologie intrins\u00e8que du nuage de points sans guidage t\u00e9l\u00e9ologique.<\/p>\n<p>Cette approche vise la d\u00e9couverte de structures latentes au sein d&rsquo;ensembles de donn\u00e9es multidimensionnels. L&rsquo;hypoth\u00e8se directrice postule que les entit\u00e9s observ\u00e9es ne se distribuent pas de mani\u00e8re al\u00e9atoire ou uniforme dans l&rsquo;espace des variables, mais qu&rsquo;elles se condensent autour de sous-vari\u00e9t\u00e9s g\u00e9om\u00e9triques ou de bassins de densit\u00e9 sp\u00e9cifiques. Ces concentrations refl\u00e8tent des r\u00e9gularit\u00e9s sous-jacentes, des m\u00e9canismes g\u00e9n\u00e9rateurs communs ou des contraintes \u00e9cologiques et psychologiques partag\u00e9es. L&rsquo;enjeu m\u00e9thodologique consiste d\u00e8s lors \u00e0 projeter la lumi\u00e8re sur ces configurations dissimul\u00e9es par l&rsquo;enchev\u00eatrement des covariations, sans introduire de biais de confirmation li\u00e9 aux attentes du mod\u00e9lisateur.<\/p>\n<p>D&rsquo;un point de vue \u00e9pist\u00e9mologique, la distinction entre analyse pr\u00e9dictive et analyse exploratoire de donn\u00e9es, telle que th\u00e9oris\u00e9e initialement par <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/www.jstor.org\/stable\/2684361\">John Tukey<\/a>, prend ici tout son sens. Tandis que la mod\u00e9lisation pr\u00e9dictive se focalise sur l&rsquo;extrapolation et la performance d&rsquo;assignation sur de nouveaux \u00e9chantillons, l&rsquo;analyse exploratoire vise l&rsquo;intelligibilit\u00e9 des ph\u00e9nom\u00e8nes et la g\u00e9n\u00e9ration d&rsquo;hypoth\u00e8ses novatrices. La classification non supervis\u00e9e ne pr\u00e9tend pas d\u00e9cr\u00e9ter une v\u00e9rit\u00e9 ontologique absolue, mais fournir une cartographie heuristique simplifi\u00e9e de la r\u00e9alit\u00e9 empirique. Elle permet au chercheur d&rsquo;appr\u00e9hender la complexit\u00e9 d&rsquo;un corpus de donn\u00e9es, d&rsquo;identifier des r\u00e9gularit\u00e9s morphologiques et d&rsquo;\u00e9laborer des cadres typologiques rigoureux susceptibles d&rsquo;\u00eatre soumis ult\u00e9rieurement \u00e0 des tests confirmatoires.<\/p>\n<h3>1.2 Divergences m\u00e9thodologiques entre classification k-means et classification hi\u00e9rarchique<\/h3>\n<p>La distinction op\u00e9rationnelle entre le partitionnement non hi\u00e9rarchique, exemplifi\u00e9 par l&rsquo;algorithme des <em>k-means<\/em> (ou <em>k-moyennes<\/em>), et la classification hi\u00e9rarchique r\u00e9side dans la g\u00e9om\u00e9trie computationnelle sous-jacente et les contraintes structurelles impos\u00e9es \u00e0 l&rsquo;espace des solutions. Dans le partitionnement <em>k-means<\/em>, l&rsquo;analyste est contraint de postuler <em>a priori<\/em> une valeur exacte pour le nombre de classes, d\u00e9sign\u00e9 par le param\u00e8tre <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>. Cette exigence repr\u00e9sente une limite substantielle en phase exploratoire, car la v\u00e9ritable granularit\u00e9 taxinomique du ph\u00e9nom\u00e8ne \u00e9tudi\u00e9 est g\u00e9n\u00e9ralement ignor\u00e9e. Une sous-estimation de <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> comprime l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 empirique et masque des sous-groupes cruciaux, tandis qu&rsquo;une sur-estimation fracture artificiellement des cat\u00e9gories continues et coh\u00e9rentes.<\/p>\n<p>Au-del\u00e0 de cette contrainte initiale, l&rsquo;algorithme des <em>k-means<\/em> repose sur un processus d&rsquo;optimisation it\u00e9ratif hautement sensible aux initialisations al\u00e9atoires des centro\u00efdes. Cette instabilit\u00e9 structurelle enferme fr\u00e9quemment la solution num\u00e9rique dans des minima locaux de la fonction de variance intra-classe, n\u00e9cessitant de multiples r\u00e9plications al\u00e9atoires pour stabiliser une solution optimale. De surcro\u00eet, le partitionnement par centro\u00efdes impose implicitement une hypoth\u00e8se d&rsquo;hyper-sph\u00e9ricit\u00e9 et d&rsquo;homog\u00e9n\u00e9it\u00e9 de volume sur les clusters g\u00e9n\u00e9r\u00e9s, rendant l&rsquo;algorithme particuli\u00e8rement vuln\u00e9rable face \u00e0 des structures de donn\u00e9es asym\u00e9triques, allong\u00e9es ou imbriqu\u00e9es de mani\u00e8re non convexe.<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;approche hi\u00e9rarchique s&rsquo;affranchit totalement de la fixation pr\u00e9alable d&rsquo;un nombre <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> d&rsquo;agr\u00e9gats. Elle d\u00e9ploie une trajectoire d&rsquo;assemblage continue et d\u00e9terministe, mat\u00e9rialis\u00e9e sous la forme d&rsquo;un arbre g\u00e9n\u00e9alogique appel\u00e9 dendrogramme. Chaque palier d&rsquo;agr\u00e9gation ou de scission refl\u00e8te une intensit\u00e9 math\u00e9matique de proximit\u00e9, permettant au praticien d&rsquo;examiner simultan\u00e9ment de multiples niveaux de granularit\u00e9 sans devoir recalculer int\u00e9gralement les partitions. Cette absence de contrainte initiale, coupl\u00e9e \u00e0 un d\u00e9terminisme absolu des calculs d\u00e8s lors que la matrice de dissimilarit\u00e9 est fix\u00e9e, conf\u00e8re \u00e0 la classification hi\u00e9rarchique une transparence conceptuelle et une flexibilit\u00e9 analytique sup\u00e9rieures pour la caract\u00e9risation fine des donn\u00e9es comportementales.<\/p>\n<h3>1.3 Applications psychologiques et comportementales des m\u00e9thodes de clustering<\/h3>\n<p>Dans les disciplines psychologiques et comportementales, la variabilit\u00e9 interindividuelle constitue \u00e0 la fois un d\u00e9fi analytique majeur et le c\u0153ur m\u00eame de l&rsquo;investigation empirique. La classification hi\u00e9rarchique y trouve une application privil\u00e9gi\u00e9e pour l&rsquo;identification de profils de personnalit\u00e9 et de typologies comportementales. Plut\u00f4t que de postuler l&rsquo;existence de types de personnalit\u00e9s rigides et mutuellement exclusifs sur une base purement intuitive, les chercheurs utilisent les m\u00e9thodes de clustering pour agr\u00e9ger des individus sur la base de patterns multidimensionnels complexes. Par exemple, au sein du mod\u00e8le des <em>Big Five<\/em>, la classification permet de v\u00e9rifier l&rsquo;\u00e9mergence empirique de constellations de traits sp\u00e9cifiques, comme les profils r\u00e9silients, surcontr\u00f4l\u00e9s ou sous-contr\u00f4l\u00e9s, validant ou infirmant les th\u00e9ories morphologiques de la personnalit\u00e9 humaine.<\/p>\n<p>L&rsquo;exploitation de la classification hi\u00e9rarchique s&rsquo;\u00e9tend avec une efficacit\u00e9 remarquable \u00e0 la segmentation des r\u00e9ponses apport\u00e9es \u00e0 des \u00e9chelles psychom\u00e9triques et des inventaires cliniques standardis\u00e9s. Les questionnaires d&rsquo;\u00e9valuation de la symptomatologie d\u00e9pressive, des manifestations anxieuses ou du syndrome d&rsquo;\u00e9puisement professionnel (burnout) g\u00e9n\u00e8rent des matrices de covariance denses o\u00f9 les individus manifestent des profils de d\u00e9tresse h\u00e9t\u00e9roclites. Le clustering permet d&rsquo;isoler des sous-populations cliniques pr\u00e9sentant des configurations symptomatiques distinctes, telles que des formes de d\u00e9pression \u00e0 pr\u00e9dominance somatique versus cognitive, autorisant ainsi une cat\u00e9gorisation diagnostique plus fine que celle propos\u00e9e par les scores globaux univari\u00e9s.<\/p>\n<p>Enfin, la recherche translationnelle et les sciences de l&rsquo;\u00e9ducation mobilisent l&rsquo;approche arborescente pour explorer l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des trajectoires th\u00e9rapeutiques, d\u00e9veloppementales ou d&rsquo;apprentissage. En int\u00e9grant des mesures longitudinales r\u00e9p\u00e9t\u00e9es au sein du calcul de dissimilarit\u00e9, la classification hi\u00e9rarchique permet de cartographier la divergence des \u00e9volutions individuelles au fil du temps. Les cliniciens peuvent d\u00e8s lors rep\u00e9rer des trajectoires d&rsquo;\u00e9chec pr\u00e9coce, des profils de r\u00e9mission spontan\u00e9e ou des patterns de r\u00e9sistance aux psychoth\u00e9rapies conventionnelles, fournissant ainsi une base empirique pour l&rsquo;adaptation personnalis\u00e9e des protocoles d&rsquo;intervention.<\/p>\n<h2>2. Taxonomie algorithmique : approches ascendantes et descendantes<\/h2>\n<h3>2.1 Classification ascendante hi\u00e9rarchique (CAH) : algorithme agglom\u00e9ratif<\/h3>\n<p>La classification ascendante hi\u00e9rarchique, d\u00e9sign\u00e9e sous l&rsquo;appellation g\u00e9n\u00e9rique d&rsquo;algorithme agglom\u00e9ratif, suit une trajectoire computationnelle partant de la base de l&rsquo;arborescence pour s&rsquo;\u00e9lever progressivement vers son sommet. Au stade initial de la proc\u00e9dure, not\u00e9 \u00e9tape z\u00e9ro, l&rsquo;algorithme consid\u00e8re que chaque observation individuelle constitue un cluster unitaire \u00e0 part enti\u00e8re, d\u00e9sign\u00e9 dans la litt\u00e9rature sous le terme de singleton. Pour un ensemble de donn\u00e9es comprenant <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> entit\u00e9s d&rsquo;observation, la configuration de d\u00e9part comporte donc rigoureusement <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> clusters disjoints, chacun caract\u00e9ris\u00e9 par une inertie interne nulle et une distance r\u00e9flexive \u00e9gale \u00e0 z\u00e9ro.<\/p>\n<p>Le processus se d\u00e9roule ensuite de mani\u00e8re it\u00e9rative selon un sch\u00e9ma d&rsquo;agr\u00e9gation glouton (<em>greedy algorithm<\/em>). \u00c0 chaque \u00e9tape de la boucle computationnelle, l&rsquo;algorithme scanne l&rsquo;ensemble de la matrice de dissimilarit\u00e9 pour identifier les deux groupes pr\u00e9sentant la plus faible distance globale, selon le crit\u00e8re de liaison pr\u00e9alablement s\u00e9lectionn\u00e9. Ces deux entit\u00e9s sont imm\u00e9diatement fusionn\u00e9es au sein d&rsquo;un nouveau cluster composite, r\u00e9duisant le nombre total de groupes de <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"k - 1\">$k &#8211; 1$<\/span>. La matrice de distance est alors actualis\u00e9e pour int\u00e9grer cette nouvelle entit\u00e9 et recalculer ses proximit\u00e9s avec l&rsquo;ensemble des clusters restants. Ce cycle de d\u00e9tection et de coalescence se r\u00e9p\u00e8te exactement <span class=\"arabpsy-math-inline\" data-latex=\"n - 1\">$n &#8211; 1$<\/span> fois, jusqu&rsquo;\u00e0 ce que la totalit\u00e9 des observations soit agglom\u00e9r\u00e9e au sein d&rsquo;un macro-cluster universel unique renfermant la population globale.<\/p>\n<p>Cette \u00e9l\u00e9gance conceptuelle s&rsquo;accompagne d&rsquo;un co\u00fbt computationnel substantiel qui gouverne la praticabilit\u00e9 de l&rsquo;approche. Sur le plan de la complexit\u00e9 temporelle, l&rsquo;impl\u00e9mentation na\u00efve de l&rsquo;algorithme agglom\u00e9ratif pr\u00e9sente une complexit\u00e9 asymptotique de l&rsquo;ordre de <span class=\"arabpsy-math-inline\" data-latex=\"O(n^3)\">$O(n^3)$<\/span>, bien que des formulations algorithmiques optimis\u00e9es, \u00e0 l&rsquo;instar de l&rsquo;algorithme de <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/comjnl\/20.4.364\">Murtagh et Legendre<\/a>, parviennent \u00e0 abaisser cette charge \u00e0 une complexit\u00e9 quadratique en <span class=\"arabpsy-math-inline\" data-latex=\"O(n^2)\">$O(n^2)$<\/span>. Sur le plan spatial, la n\u00e9cessit\u00e9 de stocker en m\u00e9moire vive la matrice triangulaire inf\u00e9rieure des distances impose une allocation m\u00e9moire en <span class=\"arabpsy-math-inline\" data-latex=\"O(n^2)\">$O(n^2)$<\/span>. Pour des cohortes d\u00e9passant plusieurs dizaines de milliers d&rsquo;individus, ces exigences mat\u00e9rielles entra\u00eenent une saturation rapide de la RAM, \u00e9rigeant la taille d&rsquo;\u00e9chantillon en contrainte technique d\u00e9terminante.<\/p>\n<h3>2.2 Classification descendante hi\u00e9rarchique : algorithme divisif<\/h3>\n<p>\u00c0 l&rsquo;oppos\u00e9 diam\u00e9tral de la logique agglom\u00e9rative, la classification descendante hi\u00e9rarchique, incarn\u00e9e historiquement par l&rsquo;algorithme DIANA (<em>Divisive Analysis<\/em>) d\u00e9velopp\u00e9 par Leonard Kaufman et Peter Rousseeuw, amorce son parcours taxinomique au sommet de la hi\u00e9rarchie. Dans cette configuration initiale, l&rsquo;int\u00e9gralit\u00e9 des <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> observations est condens\u00e9e dans un macro-cluster universel unique. L&rsquo;objectif algorithmique consiste d\u00e8s lors \u00e0 fractionner r\u00e9cursivement cet ensemble h\u00e9t\u00e9rog\u00e8ne en sous-groupes successifs, de mani\u00e8re \u00e0 maximiser la variance ou la distance inter-classes \u00e0 chaque scission spatiale.<\/p>\n<p>\u00c0 chaque it\u00e9ration de la proc\u00e9dure divisive, l&rsquo;algorithme identifie le cluster pr\u00e9sentant le diam\u00e8tre le plus \u00e9lev\u00e9, c&rsquo;est-\u00e0-dire l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 interne maximale. Au sein de ce groupe, l&rsquo;observation d\u00e9montrant la plus grande dissimilarit\u00e9 moyenne par rapport \u00e0 l&rsquo;ensemble de ses pairs est extraite pour initier un groupe dissident (<em>splinter group<\/em>). Les autres entit\u00e9s de l&rsquo;agr\u00e9gat sont ensuite r\u00e9attribu\u00e9es s\u00e9quentiellement \u00e0 ce nouveau noyau dissident si leur proximit\u00e9 moyenne avec celui-ci surpasse leur proximit\u00e9 avec le groupe d&rsquo;origine. Cette m\u00e9canique de bifurcation binaire se poursuit jusqu&rsquo;\u00e0 l&rsquo;isolement complet de chaque individu sous forme de singleton terminal.<\/p>\n<p>En d\u00e9pit de son \u00e9l\u00e9gance math\u00e9matique et de son aptitude th\u00e9orique \u00e0 capturer avec une grande acuit\u00e9 les macro-structures globales des donn\u00e9es d\u00e8s les premiers d\u00e9coupages, l&rsquo;approche divisive demeure marginale dans la recherche appliqu\u00e9e en raison de son co\u00fbt de calcul exponentiel. Le nombre de partitions binaires possibles pour scinder un ensemble de <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> \u00e9l\u00e9ments en deux sous-ensembles non vides est donn\u00e9 par la combinatoire <span class=\"arabpsy-math-inline\" data-latex=\"2^{n-1} - 1\">$2^{n-1} &#8211; 1$<\/span>. L&rsquo;exploration exhaustive de cet espace d\u00e9cisionnel devenant instantan\u00e9ment prohibitive d\u00e8s lors que <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> exc\u00e8de quelques dizaines d&rsquo;unit\u00e9s, les impl\u00e9mentations pratiques doivent recourir \u00e0 des heuristiques d&rsquo;optimisation locale strictes. Cette limitation structurelle explique la pr\u00e9pond\u00e9rance quasi h\u00e9g\u00e9monique de la classification ascendante dans les progiciels statistiques contemporains.<\/p>\n<h3>2.3 Repr\u00e9sentation conceptuelle du dendrogramme<\/h3>\n<p>Le dendrogramme, du grec <em>dendron<\/em> signifiant arbre et <em>gramma<\/em> d\u00e9signant le trac\u00e9, constitue la transcription graphique bidimensionnelle int\u00e9grale du processus de classification hi\u00e9rarchique. Sa topologie arborescente met en sc\u00e8ne l&rsquo;histoire s\u00e9quentielle des regroupements ou des divisions. Les feuilles terminales de l&rsquo;arbre correspondent aux unit\u00e9s statistiques individuelles ou aux variables analys\u00e9es. \u00c0 mesure que l&rsquo;on progresse le long de la structure, des branches horizontales et verticales s&rsquo;entrecroisent pour mat\u00e9rialiser les points de jonction o\u00f9 deux entit\u00e9s distinctes s&rsquo;unissent pour former un groupe morphologique unifi\u00e9.<\/p>\n<p>L&rsquo;axe des ordonn\u00e9es (ou des abscisses selon l&rsquo;orientation adopt\u00e9e) poss\u00e8de une s\u00e9mantique fondamentale : il quantifie la distance ultram\u00e9trique \u00e0 laquelle s&rsquo;op\u00e8re chaque \u00e9v\u00e9nement de fusion. Contrairement \u00e0 une simple distance g\u00e9om\u00e9trique, la distance ultram\u00e9trique satisfait une in\u00e9galit\u00e9 triangulaire renforc\u00e9e, formul\u00e9e par la relation <span class=\"arabpsy-math-inline\" data-latex=\"d(x, z) le max[d(x, y), d(y, z)]\">$d(x, z) le max[d(x, y), d(y, z)]$<\/span>. Cette hauteur d&rsquo;interconnexion indique l&rsquo;ampleur de la dissimilarit\u00e9 s\u00e9parant les sous-ensembles au moment pr\u00e9cis de leur agr\u00e9gation. Plus un n\u0153ud d&rsquo;articulation se situe \u00e0 une altitude \u00e9lev\u00e9e dans le dendrogramme, plus la divergence morphologique entre les groupes fusionn\u00e9s est prononc\u00e9e, signalant une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 substantielle.<\/p>\n<p>Une propri\u00e9t\u00e9 math\u00e9matique cardinale de cette repr\u00e9sentation r\u00e9side dans la monotonie des hauteurs de fusion successives. Dans des conditions standard de liaison, l&rsquo;altitude des n\u0153uds cro\u00eet de mani\u00e8re strictement monotone au fil des it\u00e9rations agglom\u00e9ratives, garantissant l&rsquo;absence d&rsquo;inversions visuelles perturbatrices (<em>crossovers<\/em>). Le dendrogramme r\u00e9ussit la prouesse graphique de compresser une structure relationnelle d\u00e9finie au sein d&rsquo;un hyper-espace multidimensionnel complexe pour la projeter fid\u00e8lement dans un plan interpr\u00e9table, pr\u00e9servant la hi\u00e9rarchisation globale des voisinages sans d\u00e9naturer les proximit\u00e9s relatives originelles.<\/p>\n<h2>3. Formalisation math\u00e9matique des m\u00e9triques de distance et matrices de dissimilarit\u00e9<\/h2>\n<h3>3.1 Distances g\u00e9om\u00e9triques pour donn\u00e9es continues<\/h3>\n<p>Le point de d\u00e9part analytique de toute classification hi\u00e9rarchique r\u00e9side dans la conversion d&rsquo;une matrice d&rsquo;observations de dimensions <span class=\"arabpsy-math-inline\" data-latex=\"n \\times p\">$n \\times p$<\/span>, comportant <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> individus d\u00e9crits par <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> covariables quantitatives, en une matrice de dissimilarit\u00e9 carr\u00e9e et sym\u00e9trique de dimensions <span class=\"arabpsy-math-inline\" data-latex=\"n \\times n\">$n \\times n$<\/span>. Dans le domaine des variables strictement continues, les distances g\u00e9om\u00e9triques issues de l&rsquo;espace m\u00e9trique r\u00e9el repr\u00e9sentent le socle classique de calcul. La m\u00e9trique la plus universellement employ\u00e9e est la distance euclidienne canonique. Pour deux vecteurs d&rsquo;observations individuelles not\u00e9s <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 <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>, la distance euclidienne se formalise ainsi :<\/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>Cette formulation mesure la longueur du segment rectiligne reliant deux points dans un espace euclidien <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>-dimensionnel. Dans de nombreux algorithmes d&rsquo;optimisation, la distance euclidienne au carr\u00e9 est pr\u00e9f\u00e9r\u00e9e afin d&rsquo;all\u00e9ger les calculs en s&rsquo;affranchissant de l&rsquo;op\u00e9rateur de racine carr\u00e9e et de conf\u00e9rer une p\u00e9nalit\u00e9 quadratique croissante aux \u00e9carts dimensionnels marqu\u00e9s. Toutefois, la m\u00e9trique euclidienne manifeste une sensibilit\u00e9 aigu\u00eb aux valeurs aberrantes (<em>outliers<\/em>), car la sommation des diff\u00e9rences \u00e9lev\u00e9es au carr\u00e9 amplifie de fa\u00e7on disproportionn\u00e9e l&rsquo;impact des variations extr\u00eames situ\u00e9es sur une seule dimension.<\/p>\n<p>Pour att\u00e9nuer cette vuln\u00e9rabilit\u00e9 structurelle, la distance de Manhattan, \u00e9galement connue sous le terme de norme <span class=\"arabpsy-math-inline\" data-latex=\"L_1\">$L_1$<\/span> ou distance en blocs de ville (<em>city-block distance<\/em>), propose une formulation g\u00e9om\u00e9trique plus robuste. Elle correspond \u00e0 la somme des valeurs absolues des \u00e9carts unidimensionnels :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"d_{\\text{Man\\hat\\tan}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\sum_{k=1}^{p} |x_{ik} - x_{jk}|\">$$d_{\\text{Man\\hat\\tan}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}|$$<\/span><\/p>\n<p>En renon\u00e7ant \u00e0 l&rsquo;\u00e9l\u00e9vation quadratique, la norme <span class=\"arabpsy-math-inline\" data-latex=\"L_1\">$L_1$<\/span> assure une invariance de pond\u00e9ration face aux \u00e9carts marginaux, ce qui en fait un instrument de choix lorsque les distributions empiriques pr\u00e9sentent des queues \u00e9paisses. Ces deux mesures constituent des cas particuliers d&rsquo;une formalisation param\u00e9trique plus globale : la distance de Minkowski d&rsquo;ordre <span class=\"arabpsy-math-inline\" data-latex=\"m\">$m$<\/span>, d\u00e9finie par :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"d_{\\text{Minkowski}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\left( \\sum_{k=1}^{p} |x_{ik} - x_{jk}|^m \\right)^{\\frac{1}{m}}\">$$d_{\\text{Minkowski}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\left( \\sum_{k=1}^{p} |x_{ik} &#8211; x_{jk}|^m \\right)^{\\frac{1}{m}}$$<\/span><\/p>\n<p>Cette g\u00e9n\u00e9ralisation permet \u00e0 l&rsquo;analyste de moduler la s\u00e9v\u00e9rit\u00e9 accord\u00e9e aux d\u00e9viations dimensionnelles en manipulant le param\u00e8tre continu <span class=\"arabpsy-math-inline\" data-latex=\"m\">$m$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"m=1\">$m=1$<\/span> restitue la distance de Manhattan, <span class=\"arabpsy-math-inline\" data-latex=\"m=2\">$m=2$<\/span> \u00e9quivaut \u00e0 la distance euclidienne, et la limite quand <span class=\"arabpsy-math-inline\" data-latex=\"m to \\infty\">$m to \\infty$<\/span> converge vers la distance de Chebyshev, domin\u00e9e par le diff\u00e9rentiel unidimensionnel maximal absolu.<\/p>\n<h3>3.2 Distances corr\u00e9lationnelles et angulaires pour variables psychom\u00e9triques<\/h3>\n<p>Dans de nombreuses investigations psychologiques et psychom\u00e9triques, la proximit\u00e9 d&rsquo;int\u00e9r\u00eat entre deux sujets ne r\u00e9side pas dans l&rsquo;\u00e9galit\u00e9 arithm\u00e9tique absolue de leurs r\u00e9ponses, mais dans l&rsquo;homologie structurale de leurs profils comportementaux. Deux individus \u00e9valuant une s\u00e9rie d&rsquo;items sur une \u00e9chelle ordinale peuvent manifester un d\u00e9calage d&rsquo;ancrage syst\u00e9matique li\u00e9 \u00e0 des styles de r\u00e9ponse contrast\u00e9s (l&rsquo;un acquies\u00e7ant globalement de fa\u00e7on g\u00e9n\u00e9reuse tandis que l&rsquo;autre mod\u00e8re ses pointages), tout en partageant une hi\u00e9rarchisation rigoureusement identique de leurs pr\u00e9f\u00e9rences. Dans un tel contexte, les distances g\u00e9om\u00e9triques classiques s&rsquo;av\u00e8rent inadapt\u00e9es, car elles interpr\u00e8tent ce d\u00e9calage d&rsquo;\u00e9chelle additif comme une dissemblance substantielle.<\/p>\n<p>La distance bas\u00e9e sur le coefficient de corr\u00e9lation lin\u00e9aire de Pearson pallie cette difficult\u00e9 en mesurant l&rsquo;affinit\u00e9 morphologique des patterns de r\u00e9ponse. Elle s&rsquo;obtient g\u00e9n\u00e9ralement en retranchant le coefficient de corr\u00e9lation lin\u00e9aire bivari\u00e9 de la valeur unitaire, sous la forme <span class=\"arabpsy-math-inline\" data-latex=\"d_{\\text{Pearson}}(\\mathbf{x}_i, \\mathbf{x}_j) = 1 - r_{ij}\">$d_{\\text{Pearson}}(\\mathbf{x}_i, \\mathbf{x}_j) = 1 &#8211; r_{ij}$<\/span>, ou par sa transformation m\u00e9trique normalis\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"d_{\\text{Pearson_metrique}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\sqrt{2(1 - r_{ij})}\">$d_{\\text{Pearson_metrique}}(\\mathbf{x}_i, \\mathbf{x}_j) = \\sqrt{2(1 &#8211; r_{ij})}$<\/span>. Cette mesure s&rsquo;annule lorsque deux profils sont strictement proportionnels et covariants, neutralisant ainsi les variations d&rsquo;ordonn\u00e9e \u00e0 l&rsquo;origine et les amplitudes globales de r\u00e9ponse pour se concentrer exclusivement sur la dynamique des variations relatives \u00e0 travers les items.<\/p>\n<p>Lorsque les donn\u00e9es empiriques d\u00e9coulent d&rsquo;\u00e9chelles de Likert ordinales ou violent les pr\u00e9somptions de normalit\u00e9 bivari\u00e9e, la distance bas\u00e9e sur le coefficient de corr\u00e9lation de rangs de Spearman constitue l&rsquo;alternative non param\u00e9trique id\u00e9ale. Elle s&rsquo;\u00e9value en convertissant pr\u00e9alablement les valeurs dimensionnelles de chaque individu en rangs internes avant d&rsquo;appliquer la formulation corr\u00e9lationnelle. D&rsquo;autre part, la distance cosinus \u00e9value l&rsquo;angle vectoriel s\u00e9parant deux observations au sein de l&rsquo;espace multidimensionnel :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\cos(\\theta) = \\frac{\\mathbf{x}_i \\cdot \\mathbf{x}_j}{|\\mathbf{x}_i|_2 |\\mathbf{x}_j|_2}\">$$\\cos(\\theta) = \\frac{\\mathbf{x}_i \\cdot \\mathbf{x}_j}{|\\mathbf{x}_i|_2 |\\mathbf{x}_j|_2}$$<\/span><\/p>\n<p>La distance cosinus, d\u00e9finie par <span class=\"arabpsy-math-inline\" data-latex=\"1 - \\cos(\\theta)\">$1 &#8211; \\cos(\\theta)$<\/span>, fait abstraction totale de la norme ou de la magnitude vectorielle pour focaliser la quantification de similarit\u00e9 sur l&rsquo;orientation directionnelle pure. Cette m\u00e9trique trouve un \u00e9cho fondamental dans les analyses textuelles, l&rsquo;exploration de corpus s\u00e9mantiques psychologiques et le traitement de donn\u00e9es comportementales de haute dimension o\u00f9 la sparsit\u00e9 matricielle est \u00e9lev\u00e9e.<\/p>\n<h3>3.3 Dissimilarit\u00e9s mixtes et distance de Gower<\/h3>\n<p>Les protocoles de recherche en sciences humaines et sociales collectent rarement des jeux de donn\u00e9es d&rsquo;une puret\u00e9 m\u00e9trique unilat\u00e9rale. Les bases de donn\u00e9es agr\u00e8gent couramment des covariables quantitatives continues (\u00e2ge, indices physiologiques), des variables cat\u00e9gorielles nominales polytomiques (appartenance culturelle, statut marital), des variables dichotomiques binaires sym\u00e9triques ou asym\u00e9triques (pr\u00e9sence ou absence de traits nosographiques) et des mesures d&rsquo;ordre ordinal (degr\u00e9s de s\u00e9v\u00e9rit\u00e9 clinique). L&rsquo;application directe d&rsquo;une m\u00e9trique euclidienne sur de tels ensembles hybrides constitue une aberration m\u00e9thodologique d\u00e9truisant la validit\u00e9 de l&rsquo;analyse.<\/p>\n<p>Pour surmonter cette contrainte complexe, <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2528823\">John C. Gower<\/a> a formul\u00e9 en 1971 un coefficient de dissimilarit\u00e9 g\u00e9n\u00e9rale permettant de combiner simultan\u00e9ment ces diverses \u00e9chelles au sein d&rsquo;une m\u00e9trique unifi\u00e9e. Le coefficient de dissimilarit\u00e9 de Gower entre l&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> et l&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span>, not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"S_{ij}\">$S_{ij}$<\/span>, est d\u00e9fini par la moyenne pond\u00e9r\u00e9e des dissimilarit\u00e9s partielles calcul\u00e9es sur l&rsquo;ensemble des <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> dimensions variables :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"S_{ij} = \\frac{\\sum_{k=1}^{p} w_{ijk} s_{ijk}}{\\sum_{k=1}^{p} w_{ijk}}\">$$S_{ij} = \\frac{\\sum_{k=1}^{p} w_{ijk} s_{ijk}}{\\sum_{k=1}^{p} w_{ijk}}$$<\/span><\/p>\n<p>Dans cette formalisation math\u00e9matique d&rsquo;une rare \u00e9l\u00e9gance, le terme <span class=\"arabpsy-math-inline\" data-latex=\"s_{ijk}\">$s_{ijk}$<\/span> repr\u00e9sente la contribution de dissimilarit\u00e9 de la <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>-i\u00e8me variable. Pour une variable quantitative ou ordinale, <span class=\"arabpsy-math-inline\" data-latex=\"s_{ijk} = |x_{ik} - x_{jk}| \/ R_k\">$s_{ijk} = |x_{ik} &#8211; x_{jk}| \/ R_k$<\/span>, o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"R_k\">$R_k$<\/span> d\u00e9signe l&rsquo;\u00e9tendue (l&rsquo;amplitude maximale observ\u00e9e) de la variable <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> au sein de la cohorte \u00e9chantillonnale globale. Pour une variable cat\u00e9gorielle nominale, <span class=\"arabpsy-math-inline\" data-latex=\"s_{ijk}\">$s_{ijk}$<\/span> adopte une logique bivalente : la valeur est fix\u00e9e \u00e0 $0$ si les deux sujets partagent la m\u00eame modalit\u00e9 et \u00e0 $1$ en cas de discordance cat\u00e9gorielle. Le param\u00e8tre de pond\u00e9ration <span class=\"arabpsy-math-inline\" data-latex=\"w_{ijk}\">$w_{ijk}$<\/span> prend g\u00e9n\u00e9ralement la valeur $1$, mais il bascule automatiquement \u00e0 $0$ si l&rsquo;une des observations pr\u00e9sente une valeur manquante pour la variable consid\u00e9r\u00e9e, ou dans le cas d&rsquo;une variable binaire asym\u00e9trique o\u00f9 l&rsquo;absence mutuelle d&rsquo;un trait rare ne doit pas \u00eatre interpr\u00e9t\u00e9e comme un indice de proximit\u00e9 structurelle.<\/p>\n<p>L&rsquo;impl\u00e9mentation op\u00e9rationnelle de cette m\u00e9trique sous l&rsquo;environnement R est encapsul\u00e9e de mani\u00e8re optimale dans la fonction <em>daisy<\/em> du package <em>cluster<\/em>. Cette fonction identifie automatiquement le type statistique de chaque colonne au sein d&rsquo;un dataframe ou autorise l&rsquo;utilisateur \u00e0 sp\u00e9cifier explicitement la nature m\u00e9trique de chaque covariable via l&rsquo;argument <em>type<\/em>. Elle g\u00e9n\u00e8re en sortie un objet de classe <em>dissimilarity<\/em>, pr\u00eat \u00e0 \u00eatre inject\u00e9 de fa\u00e7on transparente dans les algorithmes de classification hi\u00e9rarchique.<\/p>\n<h2>4. M\u00e9canismes de liaison (Linkage Methods) et dynamique d&rsquo;agr\u00e9gation<\/h2>\n<h3>4.1 Liaison simple (Single Linkage) et liaison compl\u00e8te (Complete Linkage)<\/h3>\n<p>Lors de chaque it\u00e9ration agglom\u00e9rative d&rsquo;une classification hi\u00e9rarchique, l&rsquo;algorithme doit \u00e9valuer la distance non plus entre deux points isol\u00e9s, mais entre deux grappes composites renfermant potentiellement de multiples observations. Le m\u00e9canisme r\u00e9gissant cette \u00e9valuation est formalis\u00e9 par la m\u00e9thode de liaison (<em>linkage method<\/em>). La m\u00e9thode de la liaison simple, ou r\u00e8gle du plus proche voisin, stipule que la dissimilarit\u00e9 s\u00e9parant deux clusters <span class=\"arabpsy-math-inline\" data-latex=\"A\">$A$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> \u00e9quivaut \u00e0 la distance minimale absolue existant entre n&rsquo;importe quel \u00e9l\u00e9ment de <span class=\"arabpsy-math-inline\" data-latex=\"A\">$A$<\/span> et n&rsquo;importe quel \u00e9l\u00e9ment de <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"D_{\\text{single}}(A, B) = \\min { d(x, y) : x in A, , y in B }\">$$D_{\\text{single}}(A, B) = \\min { d(x, y) : x in A, , y in B }$$<\/span><\/p>\n<p>Cette modalit\u00e9 de calcul conf\u00e8re \u00e0 l&rsquo;algorithme la capacit\u00e9 th\u00e9orique de d\u00e9tecter des clusters pr\u00e9sentant des g\u00e9om\u00e9tries non lin\u00e9aires complexes, telles que des structures concentriques ou sinueuses. Cependant, la liaison simple souffre d&rsquo;un \u00e9cueil structurel notoire d\u00e9sign\u00e9 sous l&rsquo;appellation d&rsquo;effet de cha\u00eene (<em>chaining effect<\/em>). Il suffit qu&rsquo;une s\u00e9quence continue de points \u00e9pars fasse office de pont spatial entre deux concentrations de densit\u00e9 denses pour que la liaison simple les fusionne pr\u00e9matur\u00e9ment, g\u00e9n\u00e9rant des dendrogrammes \u00e9tir\u00e9s, asym\u00e9triques et pratiquement inexploitables pour l&rsquo;interpr\u00e9tation substantielle.<\/p>\n<p>\u00c0 l&rsquo;autre extr\u00e9mit\u00e9 du spectre se situe la m\u00e9thode de la liaison compl\u00e8te, ou r\u00e8gle du voisin le plus \u00e9loign\u00e9. Ici, la dissimilarit\u00e9 inter-clusters est d\u00e9finie par la distance maximale observable entre une entit\u00e9 issue du cluster <span class=\"arabpsy-math-inline\" data-latex=\"A\">$A$<\/span> et une entit\u00e9 issue du cluster <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"D_{\\text{complete}}(A, B) = \\max { d(x, y) : x in A, , y in B }\">$$D_{\\text{complete}}(A, B) = \\max { d(x, y) : x in A, , y in B }$$<\/span><\/p>\n<p>Cette approche impose une contrainte de compacit\u00e9 spatiale particuli\u00e8rement stricte. Deux groupes ne peuvent \u00eatre fusionn\u00e9s que si le diam\u00e8tre de l&rsquo;agr\u00e9gat r\u00e9sultant n&rsquo;exc\u00e8de pas un certain seuil de proximit\u00e9. La liaison compl\u00e8te neutralise radicalement l&rsquo;effet de cha\u00eene et tend \u00e0 produire des grappes compactes, bien sph\u00e9riques et de dimensions g\u00e9om\u00e9triques \u00e9quivalentes. Toutefois, sa sensibilit\u00e9 aux valeurs aberrantes est extr\u00eame : une seule observation atypique situ\u00e9e en p\u00e9riph\u00e9rie d&rsquo;un groupe peut enfler consid\u00e9rablement la distance globale et emp\u00eacher une fusion l\u00e9gitime.<\/p>\n<h3>4.2 Liaison moyenne (Average Linkage \/ UPGMA) et liaison m\u00e9diane<\/h3>\n<p>Pour concilier la permissivit\u00e9 de la liaison simple et le conservatisme de la liaison compl\u00e8te, les m\u00e9thodologistes recourent fr\u00e9quemment \u00e0 la liaison moyenne, connue historiquement sous l&rsquo;acronyme UPGMA (<em>Unweighted Pair Group Method with Arithmetic Mean<\/em>). Cette approche \u00e9value 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> comme la moyenne arithm\u00e9tique de l&rsquo;ensemble des distances reliant toutes les paires possibles d&rsquo;observations inter-groupes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"D_{\\text{average}}(A, B) = \\frac{1}{|A| \\cdot |B|} \\sum_{x in A} \\sum_{y in B} d(x, y)\">$$D_{\\text{average}}(A, B) = \\frac{1}{|A| \\cdot |B|} \\sum_{x in A} \\sum_{y in B} d(x, y)$$<\/span><\/p>\n<p>En int\u00e9grant la globalit\u00e9 des interactions par paires, la liaison moyenne amortit l&rsquo;impact d\u00e9stabilisateur des observations singuli\u00e8res et produit des arbres \u00e9quilibr\u00e9s qui refl\u00e8tent fid\u00e8lement la densit\u00e9 relative des nuages de points. Cette m\u00e9thode est r\u00e9put\u00e9e pour sa robustesse empirique et son excellente pr\u00e9servation des distances m\u00e9triques originelles au sein de la matrice coph\u00e9n\u00e9tique.<\/p>\n<p>Il convient de distinguer l&rsquo;UPGMA de sa variante pond\u00e9r\u00e9e, le WPGMA (<em>Weighted Pair Group Method with Arithmetic Mean<\/em>), \u00e9galement d\u00e9sign\u00e9 sous le terme de m\u00e9thode de McQuitty. Dans la formule UPGMA, chaque observation individuelle contribue de mani\u00e8re \u00e9quiprobable au calcul final, ce qui conf\u00e8re un poids sup\u00e9rieur aux clusters comportant un effectif \u00e9lev\u00e9. Dans le WPGMA, lors de la fusion d&rsquo;un cluster <span class=\"arabpsy-math-inline\" data-latex=\"A\">$A$<\/span> avec un groupe r\u00e9sultant de l&rsquo;union pr\u00e9alable de <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"C\">$C$<\/span>, la distance est calcul\u00e9e simplement par la demi-somme des distances <span class=\"arabpsy-math-inline\" data-latex=\"D(A, B)\">$D(A, B)$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"D(A, C)\">$D(A, C)$<\/span>, ind\u00e9pendamment du nombre d&rsquo;individus h\u00e9berg\u00e9s dans chaque entit\u00e9. Cette pond\u00e9ration uniforme des sous-branches est particuli\u00e8rement utile lorsque l&rsquo;\u00e9chantillonnage pr\u00e9sente des d\u00e9s\u00e9quilibres artificiels de repr\u00e9sentativit\u00e9.<\/p>\n<p>Enfin, la liaison m\u00e9diane et la m\u00e9thode du centro\u00efde (UPGMC) projettent les clusters au travers de leurs points focaux g\u00e9om\u00e9triques. La m\u00e9thode du centro\u00efde mesure la distance s\u00e9parant les centres de gravit\u00e9 respectifs des deux groupes. Bien qu&rsquo;intuitive, cette formulation rec\u00e8le un risque math\u00e9matique majeur : elle peut enfreindre la monotonie de la hauteur de fusion, provoquant des inversions de branches sur le dendrogramme (un cluster parent \u00e9tant agr\u00e9g\u00e9 \u00e0 une altitude inf\u00e9rieure \u00e0 celle de ses propres sous-composants), ce qui complique s\u00e9v\u00e8rement sa lecture s\u00e9miologique.<\/p>\n<h3>4.3 M\u00e9thode de Ward : minimisation de l&rsquo;inertie intra-classe<\/h3>\n<p>Introduite par <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1963.10500845\">Joe H. Ward Jr. en 1963<\/a>, la m\u00e9thode de minimisation de la variance constitue sans conteste l&rsquo;algorithme d&rsquo;agr\u00e9gation le plus pl\u00e9biscit\u00e9 dans les sciences sociales et comportementales. Contrairement aux approches g\u00e9om\u00e9triques fond\u00e9es sur des extr\u00eames ou des moyennes locales, Ward ancre sa logique d\u00e9cisionnelle dans la th\u00e9orie de l&rsquo;analyse de la variance (ANOVA) et le th\u00e9or\u00e8me de d\u00e9composition de l&rsquo;inertie de Huygens. Selon ce principe, l&rsquo;inertie totale <span class=\"arabpsy-math-inline\" data-latex=\"I_{\\text{tot}}\">$I_{\\text{tot}}$<\/span> d&rsquo;un nuage d&rsquo;observations se d\u00e9compose rigoureusement en la somme de l&rsquo;inertie intra-classe (<span class=\"arabpsy-math-inline\" data-latex=\"I_{\\text{intra}}\">$I_{\\text{intra}}$<\/span> ou <em>Within-cluster Sum of Squares<\/em>, WSS) et de l&rsquo;inertie inter-classes (<span class=\"arabpsy-math-inline\" data-latex=\"I_{\\text{inter}}\">$I_{\\text{inter}}$<\/span> ou <em>Between-cluster Sum of Squares<\/em>, BSS) :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"I_{\\text{tot}} = I_{\\text{intra}} + I_{\\text{inter}}\">$$I_{\\text{tot}} = I_{\\text{intra}} + I_{\\text{inter}}$$<\/span><\/p>\n<p>\u00c0 chaque \u00e9tape d&rsquo;agr\u00e9gation ascendante, la fusion de deux groupes entra\u00eene n\u00e9cessairement une d\u00e9gradation de l&rsquo;homog\u00e9n\u00e9it\u00e9 interne globale, se traduisant par un accroissement direct de l&rsquo;inertie intra-classe. Le crit\u00e8re de Ward stipule que l&rsquo;algorithme doit fusionner, \u00e0 chaque it\u00e9ration, la paire de clusters dont la r\u00e9union minimise l&rsquo;augmentation de cette inertie intra-classe, ce qui \u00e9quivaut rigoureusement \u00e0 maximiser la variance inter-classes r\u00e9sultante. Le surcro\u00eet de variance induit par la fusion des groupes <span class=\"arabpsy-math-inline\" data-latex=\"A\">$A$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span>, poss\u00e9dant respectivement des effectifs <span class=\"arabpsy-math-inline\" data-latex=\"n_A\">$n_A$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"n_B\">$n_B$<\/span> et des centres de gravit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{g}_A\">$\\mathbf{g}_A$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{g}_B\">$\\mathbf{g}_B$<\/span>, s&rsquo;exprime selon la formule m\u00e9trique :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\Delta I(A, B) = \\frac{n_A n_B}{n_A + n_B} |\\mathbf{g}_A - \\mathbf{g}_B|^2\">$$\\Delta I(A, B) = \\frac{n_A n_B}{n_A + n_B} |\\mathbf{g}_A &#8211; \\mathbf{g}_B|^2$$<\/span><\/p>\n<p>Cette formalisation conf\u00e8re \u00e0 la m\u00e9thode de Ward une propension naturelle \u00e0 g\u00e9n\u00e9rer des classes remarquablement compactes, denses, homog\u00e8nes et de dimensions comparables. Dans l&rsquo;\u00e9cosyst\u00e8me R, une subtilit\u00e9 algorithmique cardinale doit \u00eatre ma\u00eetris\u00e9e lors de l&rsquo;appel de la fonction <em>hclust<\/em> : la distinction entre les arguments <em>method = \u00ab\u00a0ward.D\u00a0\u00bb<\/em> et <em>method = \u00ab\u00a0ward.D2\u00a0\u00bb<\/em>. L&rsquo;algorithme math\u00e9matique original formul\u00e9 par Ward exige que les distances inject\u00e9es soient des distances euclidiennes \u00e9lev\u00e9es au carr\u00e9. Lorsque l&rsquo;analyste utilise <em>ward.D2<\/em>, le logiciel prend en entr\u00e9e une matrice de distances euclidiennes standards <span class=\"arabpsy-math-inline\" data-latex=\"d\">$d$<\/span> et applique automatiquement l&rsquo;\u00e9l\u00e9vation au carr\u00e9 au cours de la formule de r\u00e9currence de Lance-Williams. \u00c0 l&rsquo;inverse, l&rsquo;option historique <em>ward.D<\/em> applique la formule directement sur les distances non carr\u00e9es, g\u00e9n\u00e9rant un \u00e9cart par rapport \u00e0 la v\u00e9ritable minimisation de l&rsquo;inertie de Ward. Dans la quasi-totalit\u00e9 des contextes de recherche contemporains, l&rsquo;usage exclusif de <em>ward.D2<\/em> est donc formellement prescrit.<\/p>\n<h2>5. Configuration de l&rsquo;environnement de travail et s\u00e9lection du jeu de donn\u00e9es en R<\/h2>\n<h3>5.1 S\u00e9lection et chargement des biblioth\u00e8ques R sp\u00e9cialis\u00e9es<\/h3>\n<p>Pour ex\u00e9cuter un pipeline complet de classification hi\u00e9rarchique alliant rigueur algorithmique, puissance graphique et profondeur diagnostique, le chargement pr\u00e9alable d&rsquo;une suite modulaire de packages sp\u00e9cialis\u00e9s est imp\u00e9ratif. Bien que le moteur de base de R contienne des fonctions fonctionnelles comme <em>hclust<\/em> et <em>dist<\/em>, l&rsquo;ing\u00e9nierie moderne s&rsquo;appuie sur des extensions con\u00e7ues pour enrichir l&rsquo;analytique et standardiser les flux de donn\u00e9es. Le package fondamental <em>cluster<\/em> fournit l&rsquo;infrastructure computationnelle avanc\u00e9e pour le calcul de dissimilarit\u00e9s mixtes via <em>daisy<\/em> et les m\u00e9thodes d&rsquo;agr\u00e9gation alternatives comme <em>agnes<\/em>. Parall\u00e8lement, le package <em>factoextra<\/em> propose une couche d&rsquo;abstraction moderne d\u00e9di\u00e9e \u00e0 la visualisation \u00e9l\u00e9gante des mod\u00e8les factoriels et des partitions de clustering.<\/p>\n<p>La biblioth\u00e8que <em>dendextend<\/em> repr\u00e9sente l&rsquo;outil le plus abouti pour la chirurgie et la manipulation fine des objets dendrogrammiques sous R. Elle permet d&rsquo;acc\u00e9der \u00e0 chaque n\u0153ud, branche et \u00e9tiquette de l&rsquo;arbre pour en modifier les propri\u00e9t\u00e9s s\u00e9miologiques de mani\u00e8re programmatique. En compl\u00e9ment, la collection de packages r\u00e9unie sous la banni\u00e8re du <em>tidyverse<\/em>, incluant <em>dplyr<\/em>, <em>tidyr<\/em> et <em>ggplot2<\/em>, garantit une manipulation fluide des structures tabulaires, simplifiant l&rsquo;int\u00e9gration ult\u00e9rieure des scores de clustering au sein du dataframe analytique original.<\/p>\n<p>Enfin, la question de la reproductibilit\u00e9 scientifique impose une gestion scrupuleuse des d\u00e9pendances et de l&rsquo;al\u00e9a logiciel. Bien que la classification ascendante hi\u00e9rarchique soit d\u00e9terministe dans son calcul une fois la matrice de distance \u00e9tablie, de multiples phases adjacentes de l&rsquo;analyse \u2013 telles que le r\u00e9\u00e9chantillonnage statistique (<em>bootstrap<\/em>), les algorithmes d&rsquo;imputation stochastique ou l&rsquo;initialisation de projections factorielles \u2013 d\u00e9pendent directement du g\u00e9n\u00e9rateur de nombres pseudo-al\u00e9atoires de R. L&rsquo;initialisation syst\u00e9matique de la graine computationnelle via l&rsquo;instruction <em>set.seed()<\/em> doit figurer au d\u00e9but absolu de tout script de recherche, scellant l&rsquo;invariance exacte des sorties num\u00e9riques et graphiques \u00e0 travers les r\u00e9plications.<\/p>\n<h3>5.2 Pr\u00e9sentation et inspection descriptive du jeu de donn\u00e9es<\/h3>\n<p>Pour ancrer cette d\u00e9monstration dans un r\u00e9alisme empirique tangible, nous s\u00e9lectionnons un jeu de donn\u00e9es multidimensionnel issu de la psychom\u00e9trie appliqu\u00e9e et des sciences comportementales. Consid\u00e9rons un \u00e9chantillon de <span class=\"arabpsy-math-inline\" data-latex=\"n = 300\">$n = 300$<\/span> individus ayant compl\u00e9t\u00e9 un inventaire d&rsquo;\u00e9valuation de la symptomatologie de l&rsquo;\u00e9puisement professionnel et de la sant\u00e9 psychologique au travail. Ce protocole \u00e9value six dimensions quantitatives fondamentales : l&rsquo;\u00e9puisement \u00e9motionnel (EE), la d\u00e9personnalisation ou cynisme professionnel (DP), la perte d&rsquo;accomplissement personnel (AP), le niveau d&rsquo;anxi\u00e9t\u00e9 situationnelle (ANX), les manifestations somatiques du stress (SOM) et l&rsquo;engagement cognitif au travail (ENG). L&rsquo;enjeu est d&rsquo;identifier si la population se structure spontan\u00e9ment en typologies cliniquement intelligibles, sans orienter les cat\u00e9gories par des coupures arbitraires.<\/p>\n<p>L&rsquo;amorce de la proc\u00e9dure sous R requiert une auscultation minutieuse de la morphologie du jeu de donn\u00e9es \u00e0 l&rsquo;aide des primitives d&rsquo;inspection globale. L&rsquo;application de la fonction <em>str()<\/em> renseigne sur l&rsquo;architecture interne de l&rsquo;objet, confirmant la nature matricielle ou tabulaire des donn\u00e9es et l&rsquo;ad\u00e9quation des typages de colonnes (valeurs num\u00e9riques continues). La fonction <em>glimpse()<\/em>, issue de l&rsquo;\u00e9cosyst\u00e8me <em>dplyr<\/em>, offre une visualisation condens\u00e9e des premiers vecteurs individuels, permettant de v\u00e9rifier l&rsquo;alignement des variables.<\/p>\n<p>Simultan\u00e9ment, l&rsquo;ex\u00e9cution de la commande <em>summary()<\/em> g\u00e9n\u00e8re les statistiques descriptives param\u00e9triques et non param\u00e9triques \u00e9l\u00e9mentaires : minimums, premiers quartiles, m\u00e9dianes, moyennes arithm\u00e9tiques, troisi\u00e8mes quartiles et valeurs maximales pour chaque \u00e9chelle d&rsquo;\u00e9valuation. Cette premi\u00e8re lecture permet de rep\u00e9rer d&rsquo;\u00e9ventuelles anomalies de saisie, telles que des scores exc\u00e9dant les bornes th\u00e9oriques de l&rsquo;inventaire psychom\u00e9trique, et d&rsquo;observer les disparit\u00e9s d&rsquo;amplitude globale s\u00e9parant les diff\u00e9rentes dimensions d&rsquo;investigation.<\/p>\n<h3>5.3 Exploration visuelle pr\u00e9liminaire des corr\u00e9lations<\/h3>\n<p>Avant d&rsquo;engager le calcul de la dissimilarit\u00e9 interindividuelle, il est m\u00e9thodologiquement indispensable de sonder la structure d&rsquo;intercorr\u00e9lation des variables psychom\u00e9triques. En effet, dans le calcul de la plupart des m\u00e9triques g\u00e9om\u00e9triques comme la distance euclidienne, l&rsquo;existence d&rsquo;une forte colin\u00e9arit\u00e9 entre deux ou plusieurs variables introduit une redondance informationnelle majeure. Deux \u00e9chelles hautement corr\u00e9l\u00e9es (par exemple, <span class=\"arabpsy-math-inline\" data-latex=\"r &gt; 0.85\">$r &gt; 0.85$<\/span> entre l&rsquo;\u00e9puisement \u00e9motionnel et l&rsquo;anxi\u00e9t\u00e9) capturent un construit psychologique sous-jacent quasi identique. En int\u00e9grant ces deux colonnes sans ajustement, le chercheur accorde implicitement un poids math\u00e9matique double \u00e0 cette dimension symptomatique, biaisant la topologie finale du dendrogramme au d\u00e9triment des variables ind\u00e9pendantes.<\/p>\n<p>Pour cartographier ces d\u00e9pendances multivari\u00e9es sous R, la matrice de corr\u00e9lation de Pearson ou de Spearman est g\u00e9n\u00e9r\u00e9e \u00e0 l&rsquo;aide de la fonction <em>cor()<\/em>, puis projet\u00e9e visuellement via le package <em>corrplot<\/em>. La commande <em>corrplot(cor(data), method = \u00ab\u00a0circle\u00a0\u00bb, type = \u00ab\u00a0upper\u00a0\u00bb)<\/em> d\u00e9ploie une s\u00e9miologie graphique o\u00f9 la couleur et la dimension des cercles quantifient l&rsquo;amplitude et la direction des liens de covariance. Une telle matrice permet d&rsquo;identifier d&rsquo;un coup d&rsquo;\u0153il les blocs dimensionnels synergiques et d&rsquo;anticiper la structure factorielle qui guidera la formation des clusters.<\/p>\n<p>En pr\u00e9sence d&rsquo;une matrice de dispersion bivari\u00e9e, trac\u00e9e \u00e0 l&rsquo;aide de la fonction <em>pairs()<\/em> ou de biblioth\u00e8ques avanc\u00e9es telles que <em>GGally::ggpairs()<\/em>, l&rsquo;analyste peut \u00e9galement contr\u00f4ler la lin\u00e9arit\u00e9 des relations, l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 des r\u00e9sidus et l&rsquo;existence potentielle de sous-populations multimodales d\u00e9j\u00e0 visibles sur des projections bidimensionnelles. Cette inspection confirme si certaines variables dominantes risquent d&rsquo;\u00e9craser la mesure d&rsquo;\u00e9loignement global.<\/p>\n<h2>6. Pr\u00e9traitement rigoureux : normalisation, centrage et traitement des donn\u00e9es manquantes<\/h2>\n<h3>6.1 Impact critique de l&rsquo;\u00e9chelle de mesure sur les m\u00e9triques de distance<\/h3>\n<p>L&rsquo;\u00e9cueil m\u00e9thodologique le plus d\u00e9vastateur dans l&rsquo;application de la classification hi\u00e9rarchique r\u00e9side dans l&rsquo;asym\u00e9trie non contr\u00f4l\u00e9e des variances et des \u00e9chelles d&rsquo;expression des variables. Consid\u00e9rons une analyse comportementale combinant une variable psychophysiologique, telle que la conductance cutan\u00e9e ou le rythme cardiaque mesur\u00e9 sur une amplitude s&rsquo;\u00e9tendant de $60$ \u00e0 $180$ battements par minute, et une \u00e9chelle psychom\u00e9trique standardis\u00e9e born\u00e9e entre $1$ et $5$. Si la distance euclidienne est calcul\u00e9e directement sur ces donn\u00e9es brutes, l&rsquo;\u00e9cart maximal potentiel sur la dimension cardiaque est de $120$ unit\u00e9s, tandis qu&rsquo;il n&rsquo;est que de $4$ unit\u00e9s sur l&rsquo;\u00e9chelle psychom\u00e9trique.<\/p>\n<p>Dans cette configuration, les diff\u00e9rences arithm\u00e9tiques observ\u00e9es sur le rythme cardiaque vont totalement \u00e9craser les \u00e9carts observ\u00e9s sur l&rsquo;\u00e9chelle psychom\u00e9trique lors de l&rsquo;\u00e9l\u00e9vation des diff\u00e9rences au carr\u00e9. La distance globale deviendra un reflet quasi exclusif de la variable physiologique, r\u00e9duisant l&rsquo;\u00e9chelle psychom\u00e9trique \u00e0 un bruit statistique marginal. L&rsquo;algorithme de clustering n&rsquo;op\u00e9rera d\u00e8s lors pas une synth\u00e8se multidimensionnelle \u00e9quilibr\u00e9e, mais une simple segmentation univari\u00e9e d\u00e9guis\u00e9e sous une forme multivari\u00e9e.<\/p>\n<p>Ce ph\u00e9nom\u00e8ne d&rsquo;\u00e9crasement m\u00e9trique impose une r\u00e8gle m\u00e9thodologique intangible : toute classification hi\u00e9rarchique reposant sur des distances g\u00e9om\u00e9triques doit imp\u00e9rativement \u00eatre pr\u00e9c\u00e9d\u00e9e d&rsquo;une \u00e9tape de standardisation rigoureuse. Cette transformation math\u00e9matique a pour objectif d&rsquo;homog\u00e9n\u00e9iser l&rsquo;influence relative de chaque variable dans l&rsquo;espace multidimensionnel, garantissant que chaque trait psychologique dispose d&rsquo;un pouvoir discriminant \u00e9quitable lors de la construction des topologies dendrogrammiques.<\/p>\n<h3>6.2 Techniques de centrage et r\u00e9duction des variables sous R<\/h3>\n<p>L&rsquo;environnement R propose une panoplie de techniques de transformation pour uniformiser les distributions empiriques. L&rsquo;approche standard la plus employ\u00e9e dans la litt\u00e9rature scientifique est la standardisation par le score <span class=\"arabpsy-math-inline\" data-latex=\"z\">$z$<\/span> (<em>z-score scaling<\/em>), ex\u00e9cut\u00e9e nativement par la fonction <em>scale()<\/em>. Cette proc\u00e9dure soustrait la moyenne arithm\u00e9tique de chaque variable et divise la diff\u00e9rence par l&rsquo;\u00e9cart-type empirique correspondant :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"z_{ik} = \\frac{x_{ik} - \\bar{x}_k}{s_k}\">$$z_{ik} = \\frac{x_{ik} &#8211; \\bar{x}_k}{s_k}$$<\/span><\/p>\n<p>Apr\u00e8s application de la commande <em>data_scaled &lt;- scale(data)<\/em>, chaque covariable se caract\u00e9rise par une moyenne rigoureusement \u00e9gale \u00e0 z\u00e9ro et une variance (ainsi qu&rsquo;un \u00e9cart-type) unitaire. La distance euclidienne calcul\u00e9e sur cette matrice standardis\u00e9e refl\u00e8te d\u00e8s lors les d\u00e9viations des individus par rapport au barycentre de l&rsquo;\u00e9chantillon, exprim\u00e9es en unit\u00e9s d&rsquo;\u00e9cart-type, r\u00e9tablissant une \u00e9quit\u00e9 m\u00e9trique totale entre variables.<\/p>\n<p>Cependant, la standardisation par le score <span class=\"arabpsy-math-inline\" data-latex=\"z\">$z$<\/span> s&rsquo;appuie sur la moyenne et l&rsquo;\u00e9cart-type, deux indicateurs param\u00e9triques notoirement vuln\u00e9rables \u00e0 la pr\u00e9sence de valeurs extr\u00eames. Si un \u00e9chantillon contient des valeurs aberrantes massives, l&rsquo;\u00e9cart-type sera gonfl\u00e9 artificiellement, comprimant la variance du reste de la cohorte. Pour immuniser le pr\u00e9traitement contre cette distorsion, l&rsquo;analyste peut recourir \u00e0 la standardisation robuste bas\u00e9e sur la m\u00e9diane et l&rsquo;\u00e9cart absolu m\u00e9dian (MAD, <em>Median Absolute Deviation<\/em>) :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"x_{ik}^{\\text{robuste}} = \\frac{x_{ik} - \\text{Med}_k}{\\text{MAD}_k}\">$$x_{ik}^{\\text{robuste}} = \\frac{x_{ik} &#8211; \\text{Med}_k}{\\text{MAD}_k}$$<\/span><\/p>\n<p>Sous R, cette transformation s&rsquo;impl\u00e9mente ais\u00e9ment via des fonctions personnalis\u00e9es ou par le package <em>robustbase<\/em>. Une autre strat\u00e9gie alternative, la normalisation min-max, contraint l&rsquo;ensemble des observations dans un intervalle ferm\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span> par la formule <span class=\"arabpsy-math-inline\" data-latex=\"(x_{ik} - \\min_k) \/ (\\max_k - \\min_k)\">$(x_{ik} &#8211; \\min_k) \/ (\\max_k &#8211; \\min_k)$<\/span>. Cette derni\u00e8re est particuli\u00e8rement indiqu\u00e9e lorsque la distribution des donn\u00e9es est born\u00e9e par construction et que l&rsquo;on d\u00e9sire pr\u00e9server rigoureusement l&rsquo;homoth\u00e9tie originelle des supports num\u00e9riques sans assumer une distribution gaussienne sous-jacente.<\/p>\n<h3>6.3 Gestion des valeurs manquantes et des valeurs aberrantes<\/h3>\n<p>Les fonctions algorithmiques de classification hi\u00e9rarchique sous R, telles que <em>dist()<\/em> et <em>hclust()<\/em>, ne tol\u00e8rent structurellement aucune cellule manquante (<em>NA<\/em>). D\u00e8s lors qu&rsquo;une seule composante d&rsquo;un vecteur d&rsquo;observation fait d\u00e9faut, le calcul de la distance euclidienne entre ce sujet et l&rsquo;ensemble de ses pairs devient ind\u00e9termin\u00e9. L&rsquo;analyste se trouve confront\u00e9 \u00e0 l&rsquo;arbitrage critique entre la suppression des cas non complets (<em>listwise deletion<\/em>) et l&rsquo;imputation statistique.<\/p>\n<p>L&rsquo;\u00e9limination pure et simple des sujets pr\u00e9sentant des valeurs manquantes, bien que triviale via <em>na.omit()<\/em>, introduit des biais d&rsquo;\u00e9chantillonnage s\u00e9v\u00e8res si le m\u00e9canisme g\u00e9n\u00e9rateur des donn\u00e9es manquantes n&rsquo;est pas strictement ind\u00e9pendant du ph\u00e9nom\u00e8ne \u00e9tudi\u00e9 (m\u00e9canismes MAR, <em>Missing at Random<\/em>, ou MNAR, <em>Missing Not at Random<\/em>). De plus, elle ampute la puissance statistique de l&rsquo;\u00e9tude. Il est d\u00e8s lors largement pr\u00e9f\u00e9rable de mobiliser des protocoles d&rsquo;imputation multivari\u00e9e sophistiqu\u00e9s. Dans les cadres multidimensionnels, l&rsquo;imputation par les <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> plus proches voisins (<em>k-NN<\/em>) ou l&rsquo;imputation factorielle par le package sp\u00e9cialis\u00e9 <em>missMDA<\/em> (via la fonction <em>imputePCA()<\/em> de Josse et Husson) permettent d&rsquo;estimer les valeurs absentes en exploitant la structure de corr\u00e9lation globale des variables, sans d\u00e9former la matrice de variance-covariance empirique.<\/p>\n<p>Parall\u00e8lement, la d\u00e9tection des valeurs aberrantes multivari\u00e9es (<em>multivariate outliers<\/em>) constitue une \u00e9tape non n\u00e9gociable. Un individu peut pr\u00e9senter des valeurs univari\u00e9es acceptables sur chaque item pris isol\u00e9ment, mais afficher une combinaison hautement improbable au regard de la structure g\u00e9n\u00e9rale de la population. L&rsquo;identification de ces profils s&rsquo;op\u00e8re par le calcul de la distance de Mahalanobis pour chaque observation :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"D_M(\\mathbf{x}_i) = \\sqrt{(\\mathbf{x}_i - boldsymbol{\\mu})^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x}_i - boldsymbol{\\mu})}\">$$D_M(\\mathbf{x}_i) = \\sqrt{(\\mathbf{x}_i &#8211; boldsymbol{\\mu})^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x}_i &#8211; boldsymbol{\\mu})}$$<\/span><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}\">$boldsymbol{\\mu}$<\/span> d\u00e9signe le vecteur des moyennes et <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span> la matrice de covariance de l&rsquo;\u00e9chantillon. Les observations dont la distance de Mahalanobis exc\u00e8de la valeur critique issue de la distribution du Chi-carr\u00e9 au seuil nominal (souvent <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0.001\">$\\alpha = 0.001$<\/span>) doivent \u00eatre soumises \u00e0 un examen approfondi. Si ces entit\u00e9s r\u00e9sultent d&rsquo;erreurs de saisie ou d&rsquo;\u00e9chantillonnage, leur retrait s&rsquo;impose; s&rsquo;il s&rsquo;agit de singularit\u00e9s r\u00e9elles, leur influence sur le dendrogramme doit \u00eatre surveill\u00e9e pour \u00e9viter qu&rsquo;elles ne forment des branches isol\u00e9es artificielles qui d\u00e9graderaient l&rsquo;interpr\u00e9tation des clusters principaux.<\/p>\n<h2>7. Impl\u00e9mentation pas \u00e0 pas de la classification ascendante hi\u00e9rarchique avec hclust<\/h2>\n<h3>7.1 Calcul de la matrice de distance empirique<\/h3>\n<p>Une fois les variables nettoy\u00e9es et correctement standardis\u00e9es au sein d&rsquo;un dataframe ou d&rsquo;une matrice num\u00e9rique, la premi\u00e8re phase op\u00e9rationnelle du pipeline de classification consiste \u00e0 calculer la matrice de dissimilarit\u00e9 pairwise. Dans R, l&rsquo;outil canonique d\u00e9di\u00e9 \u00e0 cette t\u00e2che est la fonction <em>dist()<\/em>. Cette primitive logicielle, hautement optimis\u00e9e en langage C sous-jacent, prend en argument le jeu de donn\u00e9es transform\u00e9 et calcule les \u00e9carts entre toutes les paires d&rsquo;observations.<\/p>\n<p>L&rsquo;argument <em>method<\/em> sp\u00e9cifie la g\u00e9om\u00e9trie m\u00e9trique retenue pour le traitement. Les options disponibles incluent <em>\u00ab\u00a0euclidean\u00a0\u00bb<\/em>, <em>\u00ab\u00a0maximum\u00a0\u00bb<\/em>, <em>\u00ab\u00a0manhattan\u00a0\u00bb<\/em>, <em>\u00ab\u00a0canberra\u00a0\u00bb<\/em>, <em>\u00ab\u00a0binary\u00a0\u00bb<\/em> ou <em>\u00ab\u00a0minkowski\u00a0\u00bb<\/em>. Pour notre \u00e9tude psychom\u00e9trique standardis\u00e9e, l&rsquo;instruction s&rsquo;ex\u00e9cute selon la syntaxe suivante :<\/p>\n<p><code>d_matrix &lt;- dist(data_scaled, method = \"euclidean\")<\/code><\/p>\n<p>Le r\u00e9sultat produit par cet appel n&rsquo;est pas un tableau rectangulaire conventionnel, mais un objet sp\u00e9cifique de classe <em>dist<\/em>. Afin d&rsquo;\u00e9conomiser l&rsquo;espace en m\u00e9moire vive, R ne stocke pas les <span class=\"arabpsy-math-inline\" data-latex=\"n^2\">$n^2$<\/span> entr\u00e9es de la matrice compl\u00e8te, mais uniquement la partie triangulaire inf\u00e9rieure stricte, excluant la diagonale principale dont les valeurs sont identiquement nulles par axiome d&rsquo;auto-similarit\u00e9. Pour examiner visuellement les premiers coefficients d&rsquo;\u00e9loignement, l&rsquo;analyste peut convertir temporairement cet objet en matrice classique via la commande <em>as.matrix(d_matrix)[1:5, 1:5]<\/em>.<\/p>\n<p>Il importe de souligner les contraintes d&rsquo;\u00e9chelle li\u00e9es \u00e0 ce calcul. La dimension de l&rsquo;objet <em>dist<\/em> cro\u00eet selon la combinatoire <span class=\"arabpsy-math-inline\" data-latex=\"n(n - 1) \/ 2\">$n(n &#8211; 1) \/ 2$<\/span>. Si pour <span class=\"arabpsy-math-inline\" data-latex=\"n = 300\">$n = 300$<\/span> individus, la matrice ne requiert que $44,850$ valeurs \u00e0 stocker (une charge d\u00e9risoire pour les syst\u00e8mes informatiques actuels), un \u00e9chantillon de <span class=\"arabpsy-math-inline\" data-latex=\"n = 50,000\">$n = 50,000$<\/span> sujets exige l&rsquo;allocation d&rsquo;environ <span class=\"arabpsy-math-inline\" data-latex=\"1.25 \\times 10^9\">$1.25 \\times 10^9$<\/span> \u00e9l\u00e9ments en double pr\u00e9cision, repr\u00e9sentant pr\u00e8s de $10$ gigaoctets de m\u00e9moire vive rien que pour la structure de distance. Cette r\u00e9alit\u00e9 computationnelle rappelle l&rsquo;importance de dimensionner judicieusement ses protocoles exp\u00e9rimentaux sous R.<\/p>\n<h3>7.2 Ex\u00e9cution de l&rsquo;algorithme agglom\u00e9ratif hclust<\/h3>\n<p>La matrice de distance empirique \u00e9tant instanci\u00e9e, le d\u00e9ploiement de l&rsquo;agglom\u00e9ration hi\u00e9rarchique s&rsquo;op\u00e8re par l&rsquo;appel de la fonction <em>hclust()<\/em>. Cette fonction prend en entr\u00e9e prioritaire l&rsquo;objet de distance pr\u00e9alablement calcul\u00e9 et applique de fa\u00e7on r\u00e9cursive le crit\u00e8re de fusion d\u00e9sign\u00e9 par le param\u00e8tre <em>method<\/em>. Dans le respect des recommandations relatives \u00e0 la minimisation de l&rsquo;inertie intra-classe, nous formulons l&rsquo;instruction :<\/p>\n<p><code>hc_ward &lt;- hclust(d_matrix, method = \"ward.D2\")<\/code><\/p>\n<p>L&rsquo;objet retourn\u00e9, de classe <em>hclust<\/em>, encapsule la trajectoire taxinomique int\u00e9grale sous une architecture structurelle extr\u00eamement compacte. Pour comprendre l&rsquo;anatomie de cette sortie algorithmique, il convient d&rsquo;en d\u00e9composer les trois composantes vectorielles et matricielles fondamentales accessibles via l&rsquo;op\u00e9rateur dollar :<\/p>\n<ul>\n<li><strong>Le composant <code>hc_ward$merge<\/code> :<\/strong> Il s&rsquo;agit d&rsquo;une matrice d&rsquo;entiers de dimensions <span class=\"arabpsy-math-inline\" data-latex=\"(n - 1) \\times 2\">$(n &#8211; 1) \\times 2$<\/span>. Chaque ligne <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> d\u00e9crit l&rsquo;op\u00e9ration de fusion ex\u00e9cut\u00e9e \u00e0 l&rsquo;\u00e9tape <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> de l&rsquo;algorithme. Un indice n\u00e9gatif <span class=\"arabpsy-math-inline\" data-latex=\"-i\">$-i$<\/span> indique que l&rsquo;entit\u00e9 agr\u00e9g\u00e9e est l&rsquo;observation singleton <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> issue du jeu de donn\u00e9es original. Un entier positif <span class=\"arabpsy-math-inline\" data-latex=\"+j\">$+j$<\/span> signale que l&rsquo;\u00e9l\u00e9ment agr\u00e9g\u00e9 est le cluster composite d\u00e9j\u00e0 constitu\u00e9 \u00e0 l&rsquo;\u00e9tape <span class=\"arabpsy-math-inline\" data-latex=\"j\">$j$<\/span> de la boucle agglom\u00e9rative. Cette organisation retrace l&rsquo;arbre g\u00e9n\u00e9alogique exact de l&rsquo;ensemble du corpus.<\/li>\n<li><strong>Le composant <code>hc_ward$height<\/code> :<\/strong> Ce vecteur num\u00e9rique de longueur <span class=\"arabpsy-math-inline\" data-latex=\"n - 1\">$n &#8211; 1$<\/span> consigne la distance ultram\u00e9trique exacte \u00e0 laquelle s&rsquo;est op\u00e9r\u00e9e chaque fusion successive. Dans le cas de l&rsquo;option <em>ward.D2<\/em>, ces hauteurs correspondent rigoureusement \u00e0 l&rsquo;accroissement de la d\u00e9composition d&rsquo;inertie induit par chaque union. L&rsquo;inspection de ce vecteur permet d&rsquo;observer les sauts qualitatifs de dissemblance au fil des \u00e9tapes.<\/li>\n<li><strong>Le composant <code>hc_ward$order<\/code> :<\/strong> Ce vecteur contient la permutation optimale des <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> identifiants individuels, calcul\u00e9e pour disposer les feuilles le long de l&rsquo;axe horizontal du trac\u00e9 sans croisement de branches, assurant ainsi une lisibilit\u00e9 topologique fluide lors du rendu graphique.<\/li>\n<\/ul>\n<h3>7.3 Impl\u00e9mentation alternative via la fonction agnes du package cluster<\/h3>\n<p>Bien que <em>hclust()<\/em> demeure la commande historique de r\u00e9f\u00e9rence, la fonction <em>agnes()<\/em> (<em>Agglomerative Nesting<\/em>), int\u00e9gr\u00e9e dans le package <em>cluster<\/em>, propose une approche fonctionnelle plus riche et polyvalente. Une singularit\u00e9 remarquable d&rsquo;<em>agnes()<\/em> est sa capacit\u00e9 \u00e0 absorber indiff\u00e9remment une matrice de distance pr\u00e9calcul\u00e9e ou directement une matrice brute d&rsquo;observations, prenant alors en charge l&rsquo;\u00e9valuation de la dissimilarit\u00e9 de mani\u00e8re transparente, y compris par le coefficient mixte de Gower si sp\u00e9cifi\u00e9.<\/p>\n<p>L&rsquo;int\u00e9r\u00eat statistique majeur d&rsquo;<em>agnes()<\/em> r\u00e9side dans le calcul automatique du coefficient d&rsquo;agglom\u00e9ration (<em>Agglomerative Coefficient<\/em>, AC). Cet indicateur scalaire, born\u00e9 strictement entre $0$ et $1$, quantifie la robustesse globale de la structure de classification mise au jour par l&rsquo;arbre. Pour chaque observation <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>, l&rsquo;algorithme \u00e9value le rapport entre sa distance de premi\u00e8re fusion et sa hauteur de rattachement finale au macro-cluster global. Plus les observations sont agr\u00e9g\u00e9es pr\u00e9cocement \u00e0 des niveaux bas par rapport \u00e0 l&rsquo;altitude terminale de l&rsquo;arbre, plus le coefficient d&rsquo;agglom\u00e9ration converge vers $1$, attestant d&rsquo;une structure de donn\u00e9es hautement polaris\u00e9e en groupes naturels distincts. Inversement, une valeur d&rsquo;AC proche de $0$ signale un nuage de points isotrope ou uniforme d\u00e9pourvu d&rsquo;organisations typologiques manifestes.<\/p>\n<p>Cette propri\u00e9t\u00e9 m\u00e9trique permet de conduire une comparaison objective et quantitative de l&rsquo;efficacit\u00e9 d&rsquo;agr\u00e9gation entre plusieurs m\u00e9thodes de liaison sur un m\u00eame jeu de donn\u00e9es. En it\u00e9rant une boucle computationnelle sous R \u00e9valuant le coefficient AC \u00e0 travers les m\u00e9thodes <em>\u00ab\u00a0ward\u00a0\u00bb<\/em>, <em>\u00ab\u00a0complete\u00a0\u00bb<\/em>, <em>\u00ab\u00a0average\u00a0\u00bb<\/em> et <em>\u00ab\u00a0single\u00a0\u00bb<\/em>, le chercheur dispose d&rsquo;une m\u00e9trique formelle pour s\u00e9lectionner le m\u00e9canisme d&rsquo;assemblage le plus repr\u00e9sentatif de la g\u00e9om\u00e9trie de son \u00e9chantillon :<\/p>\n<p><code>methods &lt;- c(\"average\", \"single\", \"complete\", \"ward\")<\/code><br \/>\n<code>sapply(methods, function(m) agnes(data_scaled, method = m)$ac)<\/code><\/p>\n<p>Dans la grande majorit\u00e9 des topologies empiriques psychom\u00e9triques, la m\u00e9thode de Ward d\u00e9gage le coefficient d&rsquo;agglom\u00e9ration le plus \u00e9lev\u00e9, confirmant sa capacit\u00e9 sup\u00e9rieure \u00e0 structurer des partitions contrast\u00e9es.<\/p>\n<h2>8. Visualisation, personnalisation et analyse s\u00e9miologique du dendrogramme<\/h2>\n<h3>8.1 G\u00e9n\u00e9ration de visualisations dendrogrammiques de base et avanc\u00e9es<\/h3>\n<p>La transcription visuelle de l&rsquo;objet <em>hclust<\/em> constitue la premi\u00e8re confrontation directe de l&rsquo;analyste avec l&rsquo;architecture classificatoire mise au jour. Dans le moteur de base de R, la simple application de la m\u00e9thode g\u00e9n\u00e9rique <em>plot(hc_ward)<\/em> engendre un dendrogramme fonctionnel. Bien que cette primitive suffise pour un diagnostic exploratoire sommaire, sa s\u00e9miologie graphique demeure aust\u00e8re : les \u00e9tiquettes se chevauchent fr\u00e9quemment sur de grands \u00e9chantillons, l&rsquo;absence de couleur obscurcit la s\u00e9paration des sous-arbres, et les possibilit\u00e9s de mise en page sont contraintes.<\/p>\n<p>Pour d\u00e9passer ces limites esth\u00e9tiques et analytiques, le package <em>factoextra<\/em> propose la fonction <em>fviz_dend()<\/em>, qui s&rsquo;appuie sur le moteur graphique de <em>ggplot2<\/em>. Cette fonction transforme l&rsquo;arbre hi\u00e9rarchique en une visualisation haute d\u00e9finition, totalement personnalisable. Elle permet d&rsquo;assigner des palettes colorim\u00e9triques distinctes aux clusters identifi\u00e9s, d&rsquo;ajuster l&rsquo;\u00e9paisseur des segments, de moduler la typographie des labels et de projeter des rectangles englobants autour des partitions choisies pour focaliser le regard sur la structure de groupe.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/hclust1.png\" alt=\"Dendrogram in R\" \/><figcaption>Dendrogram in R<\/figcaption><\/figure>\n<p>L&rsquo;ajustement g\u00e9om\u00e9trique de l&rsquo;orientation spatiale de l&rsquo;arbre participe \u00e9galement \u00e0 l&rsquo;intelligibilit\u00e9 de la figure. En param\u00e9trant <em>horiz = TRUE<\/em> au sein de <em>fviz_dend()<\/em>, l&rsquo;arbre bascule d&rsquo;une disposition verticale classique \u00e0 une disposition horizontale. Cette r\u00e9orientation est particuli\u00e8rement profitable lorsque les identifiants d&rsquo;individus ou les intitul\u00e9s de variables sont textuellement longs, permettant une lecture de gauche \u00e0 droite sans rotation incommodante de la t\u00eate ou des labels. L&rsquo;esth\u00e9tique des visualisations dendrogrammiques contemporaines permet ainsi d&rsquo;allier rigueur math\u00e9matique et efficience de transmission scientifique lors de la phase de publication.<\/p>\n<h3>8.2 Manipulation avanc\u00e9e avec la biblioth\u00e8que dendextend<\/h3>\n<p>Lorsque la structure de donn\u00e9es exige une granularit\u00e9 de personnalisation chirurgicale, le package <em>dendextend<\/em>, d\u00e9velopp\u00e9 par Tal Galili, s&rsquo;affirme comme le standard algorithmique de pointe. <em>dendextend<\/em> op\u00e8re une mutation de l&rsquo;objet natif <em>hclust<\/em> pour le convertir en un objet de classe <em>dendrogram<\/em>, ouvrant la voie \u00e0 une syntaxe fonctionnelle fluide utilisant l&rsquo;op\u00e9rateur de composition (<em>pipe<\/em>).<\/p>\n<p>Gr\u00e2ce \u00e0 la commande polyvalente <em>set()<\/em>, l&rsquo;analyste acc\u00e8de \u00e0 un contr\u00f4le complet sur chaque composant microscopique de l&rsquo;arbre :<\/p>\n<p><code>dend &lt;- as.dendrogram(hc_ward) %&gt;%<\/code><br \/>\n<code>&nbsp;&nbsp;set(\"branches_k_color\", k = 3, value = c(\"#E41A1C\", \"#377EB8\", \"#4DAF4A\")) %&gt;%<\/code><br \/>\n<code>&nbsp;&nbsp;set(\"branches_lwd\", 1.2) %&gt;%<\/code><br \/>\n<code>&nbsp;&nbsp;set(\"labels_cex\", 0.6) %&gt;%<\/code><br \/>\n<code>&nbsp;&nbsp;set(\"labels_colors\", k = 3)<\/code><\/p>\n<p>Cette flexibilit\u00e9 autorise une coloration diff\u00e9renci\u00e9e des branches selon le d\u00e9coupage en <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> classes, une graduation de l&rsquo;\u00e9paisseur des traits pour souligner des chemins critiques dans la g\u00e9n\u00e9alogie des fusions, et une gestion cibl\u00e9e des n\u0153uds d&rsquo;articulation. Pour les cohortes volumineuses regroupant plusieurs centaines d&rsquo;individus, l&rsquo;affichage lin\u00e9aire traditionnel devient satur\u00e9 et difficilement lisible. <em>dendextend<\/em> r\u00e9sout \u00e9l\u00e9gamment cette contrainte de compacit\u00e9 en autorisant le trac\u00e9 de dendrogrammes circulaires ou phylog\u00e9n\u00e9tiques via la fonction <em>circlize<\/em>. Les branches se d\u00e9ploient alors de fa\u00e7on radiale \u00e0 partir d&rsquo;un centre commun vers une circonf\u00e9rence ext\u00e9rieure, maximisant l&rsquo;utilisation de la surface graphique bidimensionnelle et restituant la structure relationnelle avec une densit\u00e9 esth\u00e9tique et informative optimale.<\/p>\n<h3>8.3 Comparaison graphique de deux dendrogrammes (Tanglegram)<\/h3>\n<p>Dans la recherche typologique, l&rsquo;analyste se trouve r\u00e9guli\u00e8rement confront\u00e9 \u00e0 la n\u00e9cessit\u00e9 de confronter deux solutions de classification concurrentes appliqu\u00e9es sur la m\u00eame cohorte. Cette confrontation peut viser \u00e0 mesurer la convergence morphologique obtenue par deux m\u00e9triques de distance divergentes (par exemple, la distance euclidienne face \u00e0 la distance de Manhattan), \u00e0 sonder la sensibilit\u00e9 de l&rsquo;\u00e9chantillon face \u00e0 deux m\u00e9thodes de liaison distinctes (Ward versus liaison moyenne), ou \u00e0 comparer les structures issues de deux vagues temporelles successives dans un protocole longitudinal.<\/p>\n<p>La biblioth\u00e8que <em>dendextend<\/em> formalise cette d\u00e9marche comparative par la g\u00e9n\u00e9ration d&rsquo;un tanglegramme (<em>tanglegram<\/em>). Cette repr\u00e9sentation visuelle dispose les deux dendrogrammes analys\u00e9s en miroir face \u00e0 face, reliant chaque observation commune situ\u00e9e \u00e0 gauche \u00e0 son homologue situ\u00e9 \u00e0 droite par des lignes de connexion transversales :<\/p>\n<p><code>dend1 &lt;- as.dendrogram(hclust(d_matrix, method = \"ward.D2\"))<\/code><br \/>\n<code>dend2 &lt;- as.dendrogram(hclust(d_matrix, method = \"complete\"))<\/code><br \/>\n<code>tanglegram(dend1, dend2, sort = TRUE, common_subtrees_color_lines = TRUE)<\/code><\/p>\n<p>L&rsquo;\u00e9valuation s\u00e9miologique repose sur l&rsquo;inspection de ces segments d&rsquo;interconnexion. Lorsque les deux arbres partagent une topologie rigoureusement congruente, les lignes de liaison s&rsquo;alignent de mani\u00e8re strictement horizontale et parall\u00e8le. \u00c0 l&rsquo;inverse, des divergences structurales induisent un enchev\u00eatrement chaotique des lignes transversales. Ce degr\u00e9 d&rsquo;enchev\u00eatrement est formalis\u00e9 math\u00e9matiquement par le coefficient d&rsquo;enchev\u00eatrement (<em>entanglement<\/em>), une valeur scalaire oscillant entre $0$ (absence totale de croisement, congruence topologique absolue) et $1$ (enchev\u00eatrement maximal, incompatibilit\u00e9 des arbres).<\/p>\n<p>Sur le plan m\u00e9trique quantitatif, cette concordance graphique est corrobor\u00e9e par le calcul du coefficient de corr\u00e9lation coph\u00e9n\u00e9tique via la fonction <em>cophenetic()<\/em>. La matrice coph\u00e9n\u00e9tique consigne les distances de fusion lues directement sur la hauteur de l&rsquo;arbre pour chaque paire d&rsquo;observations. En appliquant la commande <em>cor(d_matrix, cophenetic(hc_ward))<\/em>, le chercheur \u00e9value le coefficient de fid\u00e9lit\u00e9 avec lequel l&rsquo;arborescence dendrogrammique simplifi\u00e9e restitue les v\u00e9ritables dissimilarit\u00e9s g\u00e9om\u00e9triques multidimensionnelles de l&rsquo;espace de d\u00e9part, une corr\u00e9lation sup\u00e9rieure \u00e0 $0.75$ t\u00e9moignant g\u00e9n\u00e9ralement d&rsquo;une excellente repr\u00e9sentativit\u00e9 de la projection arborescente.<\/p>\n<h2>9. \u00c9lagage de l&rsquo;arbre et crit\u00e8res objectifs de d\u00e9termination du nombre optimal de clusters<\/h2>\n<h3>9.1 \u00c9lagage pratique du dendrogramme avec cutree<\/h3>\n<p>Le dendrogramme dans sa totalit\u00e9 offre un panorama continu allant de l&rsquo;individualit\u00e9 absolue (<span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span> singletons) \u00e0 l&rsquo;universalit\u00e9 collective (un macro-groupe unique). Pour transformer cette trajectoire relationnelle en une typologie op\u00e9rationnelle exploitable par les praticiens, il est n\u00e9cessaire de trancher transversalement l&rsquo;arbre \u00e0 un niveau d\u00e9termin\u00e9. Cette op\u00e9ration de partitionnement binaire, appel\u00e9e \u00e9lagage ou coupure de l&rsquo;arbre, est ex\u00e9cut\u00e9e sous R par la fonction <em>cutree()<\/em>.<\/p>\n<p>Deux modalit\u00e9s d&rsquo;\u00e9lagage mutuellement exclusives s&rsquo;offrent \u00e0 l&rsquo;analyste lors de l&rsquo;appel de la fonction :<\/p>\n<ul>\n<li><strong>La coupure par sp\u00e9cification du nombre de classes (argument <code>k<\/code>) :<\/strong> Le chercheur ordonne \u00e0 l&rsquo;algorithme de sectionner l&rsquo;arbre de mani\u00e8re \u00e0 isoler exactement <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> sous-arbres disjoints. La commande <em>clusters_k3 &lt;- cutree(hc_ward, k = 3)<\/em> balaye les fusions du sommet vers la base et mat\u00e9rialise une partition contenant rigoureusement trois groupes d&rsquo;observations.<\/li>\n<li><strong>La coupure par seuil de hauteur ultram\u00e9trique (argument <code>h<\/code>) :<\/strong> L&rsquo;analyste fixe une altitude de dissemblance critique sur l&rsquo;axe des hauteurs. La commande <em>clusters_h &lt;- cutree(hc_ward, h = 15)<\/em> tire une ligne s\u00e9cante horizontale \u00e0 la coordonn\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"h = 15\">$h = 15$<\/span>. Tous les sous-arbres s&rsquo;articulant en de\u00e7\u00e0 de ce plafond forment des clusters ind\u00e9pendants. Cette seconde option est particuli\u00e8rement utile lorsqu&rsquo;une distance th\u00e9orique maximale d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 acceptable est d\u00e9finie <em>a priori<\/em> dans le protocole exp\u00e9rimental.<\/li>\n<\/ul>\n<p>Le vecteur d&rsquo;entiers retourn\u00e9 par <em>cutree()<\/em>, dont la longueur co\u00efncide exactement avec la taille de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"n\">$n$<\/span>, assigne \u00e0 chaque sujet l&rsquo;identifiant num\u00e9rique de son groupe d&rsquo;appartenance (<span class=\"arabpsy-math-inline\" data-latex=\"1, 2, dots, k\">$1, 2, dots, k$<\/span>). L&rsquo;\u00e9tape finale de cette phase consiste \u00e0 fusionner de mani\u00e8re \u00e9tanche ce vecteur de labels avec le dataframe initial, formalis\u00e9e typiquement par la commande <em>data_final &lt;- data %&gt;% mutate(Cluster = as.factor(clusters_k3))<\/em>. Cette assignation permet d&rsquo;ancrer les regroupements abstraits dans le profil descriptif d&rsquo;origine pour engager les analyses d&rsquo;impact et de caract\u00e9risation statistique.<\/p>\n<h3>9.2 M\u00e9thodes heuristiques et statistiques pour identifier k optimal<\/h3>\n<p>La d\u00e9termination de la hauteur de coupe optimale, ou le choix de la valeur canonique du nombre <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> de classes, repr\u00e9sente le pivot critique conditionnant la validit\u00e9 des inf\u00e9rences subs\u00e9quentes. Laisser cette d\u00e9cision au seul arbitrage subjectif de l&rsquo;observateur face au dendrogramme expose la recherche au risque d&rsquo;arbitraire et de confirmation d&rsquo;attentes th\u00e9oriques pr\u00e9alables. Pour rationaliser ce processus, une s\u00e9rie de crit\u00e8res heuristiques et de m\u00e9triques formelles doit \u00eatre mobilis\u00e9e.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/hclust2.png\" alt=\"Optimal clusters for hierarchical clustering\" \/><figcaption>Optimal clusters for hierarchical clustering<\/figcaption><\/figure>\n<p>La premi\u00e8re approche, d&rsquo;inspiration heuristique, repose sur l&rsquo;analyse de l&rsquo;\u00e9boulis des hauteurs d&rsquo;agr\u00e9gation (<em>Elbow Method<\/em> ou m\u00e9thode du coude). En observant la courbe retra\u00e7ant la hauteur de fusion en fonction du nombre d&rsquo;\u00e9tapes (ou la d\u00e9croissance de l&rsquo;inertie intra-classe selon les <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> successifs), le chercheur identifie une rupture de pente abrupte. Un grand saut vertical entre deux paliers d&rsquo;agr\u00e9gation cons\u00e9cutifs indique que l&rsquo;algorithme tente d&rsquo;unir deux macro-groupes profond\u00e9ment h\u00e9t\u00e9rog\u00e8nes. La coupure optimale se positionne pr\u00e9cis\u00e9ment juste avant cette d\u00e9gradation majeure de l&rsquo;homog\u00e9n\u00e9it\u00e9 interne.<\/p>\n<p>Une seconde approche statistique, particuli\u00e8rement robuste, s&rsquo;appuie sur l&rsquo;optimisation du coefficient moyen de silhouette (Rousseeuw, 1987). Pour chaque observation individuelle <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>, la largeur de silhouette <span class=\"arabpsy-math-inline\" data-latex=\"s(i)\">$s(i)$<\/span> \u00e9value la compacit\u00e9 de son groupe d&rsquo;assignation en regard de son \u00e9loignement par rapport au cluster limitrophe le plus proche :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"s(i) = \\frac{b(i) - a(i)}{\\max(a(i), b(i))}\">$$s(i) = \\frac{b(i) &#8211; a(i)}{\\max(a(i), b(i))}$$<\/span><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"a(i)\">$a(i)$<\/span> mesure la distance moyenne de l&rsquo;entit\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> \u00e0 l&rsquo;ensemble des autres \u00e9l\u00e9ments de son propre cluster, et <span class=\"arabpsy-math-inline\" data-latex=\"b(i)\">$b(i)$<\/span> repr\u00e9sente la distance moyenne minimale s\u00e9parant <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span> des \u00e9l\u00e9ments d&rsquo;un autre cluster rival. La valeur de <span class=\"arabpsy-math-inline\" data-latex=\"s(i)\">$s(i)$<\/span> oscille entre <span class=\"arabpsy-math-inline\" data-latex=\"-1\">$-1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"+1\">$+1$<\/span>. Un score approchant <span class=\"arabpsy-math-inline\" data-latex=\"+1\">$+1$<\/span> indique que le sujet est remarquablement positionn\u00e9 au sein d&rsquo;un groupe dense et bien s\u00e9gr\u00e9gu\u00e9; un score nul signale que l&rsquo;individu est localis\u00e9 sur la fronti\u00e8re de d\u00e9cision spatiale; un score n\u00e9gatif atteste d&rsquo;une assignation erron\u00e9e. La fonction <em>fviz_nbclust(data_scaled, hcut, method = \u00ab\u00a0silhouette\u00a0\u00bb)<\/em> it\u00e8re ce calcul pour des valeurs croissantes de <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> et trace la moyenne globale de la cohorte, d\u00e9signant le sommet de la courbe comme le nombre <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> maximisant la qualit\u00e9 g\u00e9om\u00e9trique de la partition.<\/p>\n<p>Enfin, la statistique de l&rsquo;\u00e9cart (<em>Gap Statistic<\/em>), introduite par Tibshirani, Walther et Hastie, confronte la d\u00e9croissance de l&rsquo;inertie intra-classe empirique \u00e0 celle issue d&rsquo;une distribution de r\u00e9f\u00e9rence nulle simul\u00e9e par tirages uniformes de Monte-Carlo au sein d&rsquo;un hyper-rectangle englobant :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\text{Gap}(k) = E_n^*{\\log(W_k)} - \\log(W_k)\">$$\\text{Gap}(k) = E_n^*{\\log(W_k)} &#8211; \\log(W_k)$$<\/span><\/p>\n<p>Accessible sous R par la commande <em>clusGap()<\/em>, la statistique identifie la valeur minimale de <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> pour laquelle le gain de structuration obtenu d\u00e9passe significativement l&rsquo;esp\u00e9rance math\u00e9matique g\u00e9n\u00e9r\u00e9e par le hasard uniforme, garantissant ainsi que les regroupements isol\u00e9s ne sont pas le sous-produit d&rsquo;une illusion par\u00e9idolique sur du bruit pur.<\/p>\n<h3>9.3 \u00c9valuation exhaustive par consensus multi-indices<\/h3>\n<p>Dans de nombreuses configurations empiriques complexes, les crit\u00e8res d&rsquo;\u00e9valuation univari\u00e9s comme la silhouette, le coude ou la statistique du gap formulent des recommandations discordantes. L&rsquo;un sugg\u00e8re une coupure parcimonieuse \u00e0 deux macro-classes, tandis qu&rsquo;un autre privil\u00e9gie une granularit\u00e9 fine \u00e0 cinq micro-profils. Cette instabilit\u00e9 diagnostique plonge le chercheur dans une incertitude m\u00e9thodologique inconfortable quant au partitionnement le plus rigoureux.<\/p>\n<p>Pour d\u00e9passer ces contradictions locales, le package sp\u00e9cialis\u00e9 <em>NbClust<\/em>, mis au point par Malika Charrad et ses collaborateurs, impl\u00e9mente une proc\u00e9dure d&rsquo;arbitrage par consensus computationnel exhaustif. La fonction ma\u00eetresse <em>NbClust()<\/em> \u00e9value simultan\u00e9ment jusqu&rsquo;\u00e0 30 indices de validation interne historiques issus de la litt\u00e9rature statistique internationale, parmi lesquels l&rsquo;indice de Calinski-Harabasz, l&rsquo;indice de Davies-Bouldin, l&rsquo;indice de Dunn, le crit\u00e8re de Frey et van Groenewoud, la statistique de Hartigan, l&rsquo;indice de McClain-Rao, ainsi que les crit\u00e8res d&rsquo;information pseudo-<span class=\"arabpsy-math-inline\" data-latex=\"F\">$F$<\/span> et pseudo-<span class=\"arabpsy-math-inline\" data-latex=\"T^2\">$T^2$<\/span> :<\/p>\n<p><code>res_nbclust &lt;- NbClust(data = data_scaled, distance = \"euclidean\", min.nc = 2, max.nc = 8, method = \"ward.D2\", index = \"all\")<\/code><\/p>\n<p>Le package soumet l&rsquo;\u00e9chantillon \u00e0 l&rsquo;int\u00e9gralit\u00e9 de ces m\u00e9triques pour chaque niveau d&rsquo;\u00e9lagage potentiel entre <em>min.nc<\/em> et <em>max.nc<\/em>, puis consigne le vote individuel de chaque indicateur math\u00e9matique. En appliquant la r\u00e8gle de la majorit\u00e9 (<em>majority rule<\/em>), <em>NbClust<\/em> synth\u00e9tise les r\u00e9sultats dans un histogramme r\u00e9capitulatif identifiant le nombre exact de classes d\u00e9sign\u00e9 par le plus large consensus m\u00e9thodologique.<\/p>\n<p>Cette validation coll\u00e9giale fournit une l\u00e9gitimit\u00e9 empirique indiscutable \u00e0 la d\u00e9cision de coupure finale. Elle permet d&rsquo;\u00e9quilibrer judicieusement le compromis classique entre la parcimonie th\u00e9orique (pr\u00e9f\u00e9rer un nombre r\u00e9duit de groupes pour en faciliter la communication clinique ou organisationnelle) et la granularit\u00e9 descriptive (conserver suffisamment de nuances pour ne pas dissoudre des sous-profils symptomatiques critiques).<\/p>\n<h2>10. Caract\u00e9risation multidimensionnelle et interpr\u00e9tation substantielle des profils obtenus<\/h2>\n<h3>10.1 Profilage statistique et calcul des indicateurs de tendance centrale par classe<\/h3>\n<p>Une partition algorithmique n&rsquo;acquiert une valeur scientifique que lorsque ses groupes constituants sont traduits en propositions th\u00e9oriques tangibles. Une fois les labels de clusters assign\u00e9s aux observations via <em>cutree()<\/em>, l&rsquo;analyste engage le profilage descriptif en calculant les indicateurs statistiques conditionnels pour chaque variable d&rsquo;int\u00e9r\u00eat. Sous l&rsquo;\u00e9cosyst\u00e8me <em>dplyr<\/em>, ce calcul s&rsquo;ex\u00e9cute avec fluidit\u00e9 par le pipeline :<\/p>\n<p><code>profils_stats &lt;- data_final %&gt;%<\/code><br \/>\n<code>&nbsp;&nbsp;group_by(Cluster) %&gt;%<\/code><br \/>\n<code>&nbsp;&nbsp;summarise(across(where(is.numeric), list(mean = mean, sd = sd, median = median)))<\/code><\/p>\n<p>Ce tableau d&rsquo;agr\u00e9gation dresse l&rsquo;anatomie param\u00e9trique de chaque groupe. Pour d\u00e9terminer quelles dimensions d\u00e9finissent de mani\u00e8re statistiquement significative la sp\u00e9cificit\u00e9 de chaque cluster, le recours aux valeurs-tests (<em>v-test<\/em>), encapsul\u00e9es dans la fonction <em>catdes()<\/em> du package <em>FactoMineR<\/em>, est particuli\u00e8rement f\u00e9cond. Le <em>v-test<\/em> compare la moyenne d&rsquo;une variable continue au sein d&rsquo;un cluster sp\u00e9cifique <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{x}_k\">$\\bar{x}_k$<\/span> \u00e0 la moyenne g\u00e9n\u00e9rale de l&rsquo;ensemble de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{X}\">$\\bar{X}$<\/span>. Sous l&rsquo;hypoth\u00e8se nulle d&rsquo;homog\u00e9n\u00e9it\u00e9 globale, la statistique suit asymptotiquement une loi normale centr\u00e9e r\u00e9duite :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"v = \\frac{\\bar{x}_k - \\bar{X}}{\\sqrt{\\frac{s^2}{n_k} \\left( \\frac{N - n_k}{N - 1} \\right)}}\">$$v = \\frac{\\bar{x}_k &#8211; \\bar{X}}{\\sqrt{\\frac{s^2}{n_k} \\left( \\frac{N &#8211; n_k}{N &#8211; 1} \\right)}}$$<\/span><\/p>\n<p>Une valeur-test sup\u00e9rieure en valeur absolue \u00e0 $1.96$ ou $2.58$ signale que la variable consid\u00e9r\u00e9e est respectivement sur-repr\u00e9sent\u00e9e ou sous-repr\u00e9sent\u00e9e au sein du profil \u00e0 un seuil d&rsquo;erreur de <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0.05\">$\\alpha = 0.05$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0.01\">$\\alpha = 0.01$<\/span>. Ce diagnostic ordonne automatiquement les variables par puissance de caract\u00e9risation, guidant l&rsquo;analyste vers les traits cardinaux qui diff\u00e9rencient substantiellement chaque entit\u00e9 sociologique ou clinique.<\/p>\n<p>Lorsque la base d&rsquo;analyse h\u00e9berge \u00e9galement des covariables cat\u00e9gorielles d&rsquo;arri\u00e8re-plan (sexe, cat\u00e9gorie socioprofessionnelle, comorbidit\u00e9s d\u00e9clar\u00e9es), la caract\u00e9risation s&rsquo;enrichit de tableaux crois\u00e9s de contingence croisant ces variables avec le label de clustering. L&rsquo;ex\u00e9cution de tests d&rsquo;ind\u00e9pendance du Chi-carr\u00e9 de Pearson via <em>chisq.test()<\/em> et l&rsquo;extraction des r\u00e9sidus standardis\u00e9s ajust\u00e9s r\u00e9v\u00e8lent si certaines cat\u00e9gories d\u00e9mographiques manifestent une attraction ou une r\u00e9pulsion \u00e9lective envers certains profils comportementaux.<\/p>\n<h3>10.2 Projection visuelle des clusters dans l&rsquo;espace multidimensionnel r\u00e9duit<\/h3>\n<p>L&rsquo;inspection tabulaire des moyennes conditionnelles, bien qu&rsquo;indispensable, \u00e9choue \u00e0 transmettre intuitivement la mani\u00e8re dont les clusters coexistent et se s\u00e9parent dans l&rsquo;espace multidimensionnel. La conjonction synergique de la classification hi\u00e9rarchique et de l&rsquo;Analyse en Composantes Principales (ACP) offre une solution g\u00e9om\u00e9trique remarquable pour r\u00e9soudre cette contrainte visuelle. Cette approche hybride est conceptualis\u00e9e sous l&rsquo;intitul\u00e9 de Classification Hi\u00e9rarchique sur Composantes Principales (HCPC), promue par l&rsquo;\u00e9cole d&rsquo;analyse des donn\u00e9es fran\u00e7aise.<\/p>\n<p>La d\u00e9marche consiste \u00e0 r\u00e9aliser une ACP sur les variables standardis\u00e9es pr\u00e9alablement \u00e0 la visualisation, puis \u00e0 projeter les assignations de clusters directement sur le premier plan factoriel constitu\u00e9 par les axes principaux $1$ et $2$. \u00c0 l&rsquo;aide de la fonction <em>fviz_cluster()<\/em> de <em>factoextra<\/em>, le mod\u00e9lisateur g\u00e9n\u00e8re un graphique de dispersion enrichi :<\/p>\n<p><code>fviz_cluster(list(data = data_scaled, cluster = clusters_k3), ellipse.type = \"convex\", palette = \"jco\", ggtheme = theme_minimal())<\/code><\/p>\n<p>Sur ce trac\u00e9 factoriel, chaque point repr\u00e9sente un individu, color\u00e9 selon son appartenance typologique. Les groupes sont ceints par des ellipses de confiance param\u00e9triques (par exemple des ellipses de concentration \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"95%\">$95%$<\/span>) ou par leurs enveloppes convexes respectives. L&rsquo;examen des axes factoriels fournit une cl\u00e9 de lecture substantielle imm\u00e9diate : en confrontant les corr\u00e9lations des variables initiales avec les facteurs latents (le cercle des corr\u00e9lations), l&rsquo;analyste d\u00e9crypte la nature des forces sous-jacentes qui provoquent l&rsquo;\u00e9loignement des clusters sur le plan cart\u00e9sien.<\/p>\n<p>Cette visualisation permet en outre d&rsquo;\u00e9valuer la s\u00e9parabilit\u00e9 effective des classes : des ellipses totalement disjointes attestent de typologies nettes, \u00e9tanches et exemptes de recouvrement, tandis que des ellipses se chevauchant partiellement illustrent des zones de transition morphologique o\u00f9 s&rsquo;agr\u00e8gent des profils comportementaux interm\u00e9diaires ou atypiques.<\/p>\n<h3>10.3 Traduction th\u00e9orique des profils dans le domaine d&rsquo;application<\/h3>\n<p>La derni\u00e8re \u00e9tape du cycle analytique r\u00e9side dans la conversion des constats statistiques en constructs th\u00e9oriques identifi\u00e9s, nomm\u00e9s et contextualis\u00e9s. Dans le cadre de notre exemple portant sur l&rsquo;inventaire de sant\u00e9 au travail, supposons que l&rsquo;\u00e9lagage consensuel \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"k = 3\">$k = 3$<\/span> r\u00e9v\u00e8le les caract\u00e9ristiques de profilage suivantes :<\/p>\n<ul>\n<li><strong>Le Groupe 1 (45 % de la cohorte) :<\/strong> D\u00e9montre des scores d&rsquo;\u00e9puisement \u00e9motionnel (EE) et d&rsquo;anxi\u00e9t\u00e9 (ANX) significativement inf\u00e9rieurs \u00e0 la moyenne (<span class=\"arabpsy-math-inline\" data-latex=\"v\\text{-test}  +6\">$v\\text{-test}  +6$<\/span>). Ce profil est baptis\u00e9 en termes substantiels : <em>\u00ab Profil R\u00e9silient et Engag\u00e9 \u00bb<\/em>.<\/li>\n<li><strong>Le Groupe 2 (35 % de la cohorte) :<\/strong> Manifeste un \u00e9puisement \u00e9motionnel mod\u00e9r\u00e9, mais se distingue par des valeurs exceptionnellement \u00e9lev\u00e9es sur l&rsquo;\u00e9chelle de cynisme professionnel et de d\u00e9personnalisation (DP, <span class=\"arabpsy-math-inline\" data-latex=\"v\\text{-test} &gt; +7\">$v\\text{-test} &gt; +7$<\/span>), avec un engagement effondr\u00e9. Ce cluster incarne la typologie : <em>\u00ab Profil de D\u00e9sengagement D\u00e9fensif \u00bb<\/em>.<\/li>\n<li><strong>Le Groupe 3 (20 % de la cohorte) :<\/strong> Cumule des valeurs critiques explosives sur l&rsquo;\u00e9puisement \u00e9motionnel, les sympt\u00f4mes somatiques du stress (SOM) et l&rsquo;anxi\u00e9t\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"v\\text{-test} &gt; +8\">$v\\text{-test} &gt; +8$<\/span>), accompagn\u00e9es d&rsquo;un effondrement total de l&rsquo;accomplissement personnel. Cette classe repr\u00e9sente indubitablement : <em>\u00ab Profil de D\u00e9tresse S\u00e9v\u00e8re et Burnout Av\u00e9r\u00e9 \u00bb<\/em>.<\/li>\n<\/ul>\n<p>Cette d\u00e9nomination conceptuelle ancre l&rsquo;analyse de donn\u00e9es dans la litt\u00e9rature clinique et comportementale \u00e9tablie. Le chercheur peut d\u00e8s lors confronter cette tripartition aux mod\u00e8les th\u00e9oriques existants (par exemple, le mod\u00e8le Exigences-Ressources au travail de Demerouti et Bakker), \u00e9valuant si les regroupements empiriques corroborent les m\u00e9canismes \u00e9tiologiques postul\u00e9s.<\/p>\n<p>De surcro\u00eet, ces typologies d\u00e9ploient des implications directes pour l&rsquo;\u00e9laboration d&rsquo;interventions cibl\u00e9es. Au lieu d&rsquo;administrer des protocoles pr\u00e9ventifs uniformes \u00e0 l&rsquo;ensemble d&rsquo;une organisation, les services de sant\u00e9 au travail peuvent moduler leurs r\u00e9ponses : les individus du profil de d\u00e9sengagement b\u00e9n\u00e9ficieront d&rsquo;un r\u00e9am\u00e9nagement manag\u00e9rial et d&rsquo;une red\u00e9finition du sens des t\u00e2ches, tandis que les membres du profil de d\u00e9tresse s\u00e9v\u00e8re feront l&rsquo;objet d&rsquo;une prise en charge m\u00e9dico-psychologique individualis\u00e9e d&rsquo;urgence. La classification hi\u00e9rarchique transcende ainsi le statut d&rsquo;exercice statistique pour s&rsquo;\u00e9riger en instrument d\u00e9cisionnel \u00e0 haute valeur ajout\u00e9e op\u00e9rationnelle.<\/p>\n<h2>11. Validation statistique de la stabilit\u00e9 et r\u00e9plicabilit\u00e9 des clusters<\/h2>\n<h3>11.1 \u00c9valuation de la stabilit\u00e9 par r\u00e9\u00e9chantillonnage bootstrap<\/h3>\n<p>Une vuln\u00e9rabilit\u00e9 structurelle inh\u00e9rente \u00e0 la classification ascendante hi\u00e9rarchique tient \u00e0 son d\u00e9terminisme glouton : d\u00e8s lors qu&rsquo;une fusion est scell\u00e9e \u00e0 une altitude basse, elle est irr\u00e9versible et contraint d\u00e9finitivement l&rsquo;ensemble des ramifications sup\u00e9rieures. Si un sous-groupe se forme en raison d&rsquo;une fluctuation d&rsquo;\u00e9chantillonnage locale ou d&rsquo;un artefact marginal dans les donn\u00e9es, cette distorsion est propag\u00e9e jusqu&rsquo;au sommet du dendrogramme. Il est d\u00e8s lors capital d&rsquo;\u00e9valuer la stabilit\u00e9 inf\u00e9rentielle de chaque n\u0153ud arborescent par des protocoles de r\u00e9\u00e9chantillonnage stochastique (<em>bootstrapping<\/em>).<\/p>\n<p>La biblioth\u00e8que <em>pvclust<\/em>, d\u00e9velopp\u00e9e par Ryota Suzuki et Hidetoshi Shimodaira, s&rsquo;impose comme la r\u00e9f\u00e9rence m\u00e9thodologique absolue sous R pour mener cette validation. Par l&rsquo;appel de la fonction <em>pvclust()<\/em>, le logiciel g\u00e9n\u00e8re un nombre massif d&rsquo;\u00e9chantillons bootstrap (typiquement <span class=\"arabpsy-math-inline\" data-latex=\"B = 1000\">$B = 1000$<\/span> r\u00e9plications) par tirage avec remise au sein du jeu d&rsquo;observations multidimensionnel, recalculant \u00e0 chaque \u00e9tape l&rsquo;arborescence compl\u00e8te :<\/p>\n<p><code>fit_pv &lt;- pvclust(t(data_scaled), method.hclust = \"ward.D2\", method.dist = \"euclidean\", nboot = 1000)<\/code><\/p>\n<p>L&rsquo;algorithme de <em>pvclust<\/em> calcule pour chaque n\u0153ud du dendrogramme deux indicateurs de significativit\u00e9 statistique aux propri\u00e9t\u00e9s fondamentales :<\/p>\n<ul>\n<li><strong>La valeur p au bootstrap standard (BP, <em>Bootstrap Probability<\/em>) :<\/strong> Elle correspond \u00e0 la fr\u00e9quence relative simple avec laquelle le sous-cluster sp\u00e9cifique r\u00e9appara\u00eet \u00e0 l&rsquo;identique au fil des <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> r\u00e9plications bootstrap. Bien qu&rsquo;intuitive, cette mesure pr\u00e9sente un biais syst\u00e9matique n\u00e9gatif, sous-estimant la v\u00e9ritable significativit\u00e9 du cluster lorsque la taille d&rsquo;\u00e9chantillon ou la variabilit\u00e9 des donn\u00e9es est mod\u00e9r\u00e9e.<\/li>\n<li><strong>La valeur p approximativement non biais\u00e9e (AU, <em>Approximately Unbiased p-value<\/em>) :<\/strong> Calcul\u00e9e au moyen d&rsquo;un r\u00e9\u00e9chantillonnage bootstrap multi-\u00e9chelle faisant varier de mani\u00e8re contr\u00f4l\u00e9e la taille des pseudo-\u00e9chantillons tir\u00e9s, la valeur p AU corrige le biais inh\u00e9rent au bootstrap classique par une m\u00e9thode d&rsquo;approximation g\u00e9om\u00e9trique asymptotique.<\/li>\n<\/ul>\n<p>Lors de la visualisation via <em>plot(fit_pv)<\/em>, les valeurs p sont imprim\u00e9es au-dessus de chaque articulation arborescente. L&rsquo;analyste se focalise sur les n\u0153uds arborant une valeur p AU sup\u00e9rieure ou \u00e9gale \u00e0 $0.95$ (correspondant au seuil conventionnel <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0.05\">$\\alpha = 0.05$<\/span>). Les sous-arbres satisfaisant cette contrainte, pouvant \u00eatre encadr\u00e9s automatiquement par l&rsquo;instruction <em>pvrect(fit_pv, alpha = 0.95)<\/em>, sont consid\u00e9r\u00e9s comme des configurations empiriquement robustes et hautement reproductibles, tandis que les embranchements \u00e0 faible valeur AU doivent \u00eatre interpr\u00e9t\u00e9s avec une extr\u00eame circonspection car susceptibles de r\u00e9sulter du seul bruit d&rsquo;\u00e9chantillonnage.<\/p>\n<h3>11.2 Indices de coh\u00e9rence interne et compacit\u00e9 des partitions<\/h3>\n<p>Au-del\u00e0 de la stabilit\u00e9 stochastique globale des n\u0153uds, la validation d&rsquo;une partition typologique requiert l&rsquo;\u00e9valuation de sa qualit\u00e9 morphologique interne. Une partition est statistiquement solide si elle conjugue une compacit\u00e9 maximale au sein de chaque groupe (faible dispersion interne) et une s\u00e9paration maximale entre les groupes (isolement spatial net).<\/p>\n<p>L&rsquo;inspection d\u00e9taill\u00e9e des largeurs de silhouette individuelles, obtenue par l&rsquo;extraction des donn\u00e9es via <em>cluster::silhouette(clusters_k3, d_matrix)<\/em> et visualis\u00e9e par <em>factoextra::fviz_silhouette()<\/em>, met en lumi\u00e8re la coh\u00e9rence individuelle des sujets. Le trac\u00e9 range les individus par cluster et affiche horizontalement leurs coefficients <span class=\"arabpsy-math-inline\" data-latex=\"s(i)\">$s(i)$<\/span>. Cette visualisation identifie imm\u00e9diatement les unit\u00e9s statistiques marginales arborant des barres n\u00e9gatives : ces individus correspondent \u00e0 des profils aberrants ou mal class\u00e9s, positionn\u00e9s plus pr\u00e8s d&rsquo;un cluster voisin que du centre de leur propre groupe. Un cluster exhibant une proportion notable de coefficients n\u00e9gatifs doit \u00eatre consid\u00e9r\u00e9 comme fragile ou artificiel.<\/p>\n<p>D&rsquo;autres m\u00e9triques math\u00e9matiques formalisent ce ratio compacit\u00e9\/s\u00e9paration :<\/p>\n<ul>\n<li><strong>L&rsquo;indice de Dunn :<\/strong> D\u00e9fini par le rapport entre la distance inter-clusters minimale (la s\u00e9paration la plus \u00e9troite entre deux groupes) et le diam\u00e8tre intra-cluster maximal (la plus grande dispersion interne observ\u00e9e au sein d&rsquo;un cluster). Un indice de Dunn \u00e9lev\u00e9 signale une partition exemplaire, car m\u00eame les deux classes les plus proches demeurent plus \u00e9loign\u00e9es l&rsquo;une de l&rsquo;autre que l&rsquo;extension maximale de n&rsquo;importe quel groupe.<\/li>\n<li><strong>L&rsquo;indice de Davies-Bouldin :<\/strong> Il \u00e9value pour chaque paire de classes le ratio de leur dispersion interne cumul\u00e9e sur la distance s\u00e9parant leurs centro\u00efdes, retenant pour chaque groupe la pire configuration d&rsquo;adjacence avant d&rsquo;en calculer la moyenne g\u00e9n\u00e9rale. Un indice de Davies-Bouldin faible indique que les classes sont globalement bien s\u00e9par\u00e9es par rapport \u00e0 leur taille respective.<\/li>\n<\/ul>\n<p>Ces indicateurs, disponibles dans les packages <em>clValid<\/em> ou <em>clusterCrit<\/em>, permettent de soumettre la partition \u00e0 des tests de sensibilit\u00e9 en injectant des perturbations al\u00e9atoires mineures (bruit gaussien l\u00e9ger) dans la matrice de donn\u00e9es originales afin de surveiller la r\u00e9silience des coefficients de qualit\u00e9 morphologique.<\/p>\n<h3>11.3 Reproductibilit\u00e9 sur sous-\u00e9chantillons ind\u00e9pendants<\/h3>\n<p>Le test ultime de la validit\u00e9 d&rsquo;une typologie r\u00e9side dans sa capacit\u00e9 de g\u00e9n\u00e9ralisation en dehors de l&rsquo;\u00e9chantillon exact qui a servi \u00e0 son estimation algorithmique. Si des profils comportementaux repr\u00e9sentent des entit\u00e9s \u00e9cologiques r\u00e9elles et non de simples idiosyncrasies d&rsquo;un jeu de donn\u00e9es particulier, ils doivent \u00e9merger avec une structure congruente lorsque l&rsquo;analyse est r\u00e9pliqu\u00e9e sur une population ind\u00e9pendante issue de la m\u00eame population cible.<\/p>\n<p>Pour operationaliser ce contr\u00f4le sous R, un protocole de validation crois\u00e9e par partitionnement al\u00e9atoire de l&rsquo;\u00e9chantillon (<em>split-sample validation<\/em>) est instaur\u00e9. L&rsquo;\u00e9chantillon initial de <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span> observations est scind\u00e9 al\u00e9atoirement en deux sous-ensembles \u00e9quivalents : un \u00e9chantillon de d\u00e9rivation (ou d&rsquo;apprentissage, <span class=\"arabpsy-math-inline\" data-latex=\"N_1 = N\/2\">$N_1 = N\/2$<\/span>) et un \u00e9chantillon de validation (ou de r\u00e9plication, <span class=\"arabpsy-math-inline\" data-latex=\"N_2 = N\/2\">$N_2 = N\/2$<\/span>). Le pipeline complet de classification hi\u00e9rarchique (standardisation, calcul de distance, <em>hclust<\/em> avec Ward, \u00e9lagage \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> classes) est ex\u00e9cut\u00e9 de mani\u00e8re totalement autonome sur l&rsquo;\u00e9chantillon de d\u00e9rivation <span class=\"arabpsy-math-inline\" data-latex=\"N_1\">$N_1$<\/span>.<\/p>\n<p>Pour \u00e9valuer la r\u00e9plicabilit\u00e9 sur l&rsquo;\u00e9chantillon de validation <span class=\"arabpsy-math-inline\" data-latex=\"N_2\">$N_2$<\/span>, deux strat\u00e9gies sont couramment combin\u00e9es :<\/p>\n<ul>\n<li>On ex\u00e9cute une classification hi\u00e9rarchique ind\u00e9pendante sur <span class=\"arabpsy-math-inline\" data-latex=\"N_2\">$N_2$<\/span> pour obtenir une partition empirique directe.<\/li>\n<li>Parall\u00e8lement, on projette les individus de <span class=\"arabpsy-math-inline\" data-latex=\"N_2\">$N_2$<\/span> au sein des clusters d\u00e9riv\u00e9s de <span class=\"arabpsy-math-inline\" data-latex=\"N_1\">$N_1$<\/span> en assignant chaque nouveau sujet au centre de gravit\u00e9 le plus proche issu du premier mod\u00e8le (classification par plus proche centro\u00efde).<\/li>\n<\/ul>\n<p>Le degr\u00e9 de concordance entre la partition empirique directe de <span class=\"arabpsy-math-inline\" data-latex=\"N_2\">$N_2$<\/span> et l&rsquo;assignation issue du mod\u00e8le de <span class=\"arabpsy-math-inline\" data-latex=\"N_1\">$N_1$<\/span> est ensuite rigoureusement quantifi\u00e9 \u00e0 l&rsquo;aide de l&rsquo;Indice de Rand Ajust\u00e9 (ARI, <em>Adjusted Rand Index<\/em>), accessible via le package <em>mclust<\/em> ou <em>fossil<\/em> :<\/p>\n<p><code>ari_score &lt;- adjustedRandIndex(partition_directe, partition_projetee)<\/code><\/p>\n<p>L&rsquo;indice de Rand ajust\u00e9 mesure la proportion de paires d&rsquo;individus class\u00e9es de mani\u00e8re concordante (soit group\u00e9es ensemble, soit s\u00e9par\u00e9es dans les deux partitions), corrig\u00e9e de l&rsquo;effet du hasard. Un score ARI sup\u00e9rieur \u00e0 $0.70$ atteste d&rsquo;une excellente reproductibilit\u00e9 structurelle, conf\u00e9rant une assise scientifique irr\u00e9futable \u00e0 la typologie et validant sa pertinence pour une extrapolation \u00e9largie aux populations cibles.<\/p>\n<h2>12. Extensions m\u00e9thodologiques : grands jeux de donn\u00e9es et approches hybrides<\/h2>\n<h3>12.1 Classification hi\u00e9rarchique sur grandes bases de donn\u00e9es : l&rsquo;algorithme CLARA<\/h3>\n<p>L&rsquo;expansion vertigineuse des volumes de donn\u00e9es dans les registres m\u00e9dico-administratifs, les banques de donn\u00e9es g\u00e9nomiques ou les plateformes num\u00e9riques comportementales place les algorithmes hi\u00e9rarchiques conventionnels face au \u00ab mur de la m\u00e9moire \u00bb. D\u00e8s lors qu&rsquo;une base atteint <span class=\"arabpsy-math-inline\" data-latex=\"N = 100,000\">$N = 100,000$<\/span> observations, le calcul et le stockage de la matrice de dissimilarit\u00e9 pairwise exigent pr\u00e8s de <span class=\"arabpsy-math-inline\" data-latex=\"5 \\times 10^9\">$5 \\times 10^9$<\/span> coefficients \u00e0 m\u00e9moriser en double pr\u00e9cision, une charge qui exc\u00e8de instantan\u00e9ment les capacit\u00e9s des stations de travail standards et provoque l&rsquo;\u00e9chec computationnel de la fonction <em>dist()<\/em>.<\/p>\n<p>Pour briser ce goulot d&rsquo;\u00e9tranglement tout en pr\u00e9servant l&rsquo;int\u00e9grit\u00e9 de la philosophie de regroupement, Leonard Kaufman et Peter Rousseeuw ont con\u00e7u l&rsquo;algorithme CLARA (<em>Clustering Large Applications<\/em>). Impl\u00e9ment\u00e9 nativement dans le package <em>cluster<\/em> via la fonction <em>clara()<\/em>, cet algorithme contourne le probl\u00e8me de complexit\u00e9 spatiale en s&rsquo;appuyant sur un sch\u00e9ma d&rsquo;\u00e9chantillonnage probabiliste it\u00e9ratif :<\/p>\n<p><code>clara_res &lt;- clara(big_data_scaled, k = 4, metric = \"euclidean\", samples = 50, pamLike = TRUE)<\/code><\/p>\n<p>Plut\u00f4t que d&rsquo;\u00e9valuer la matrice globale sur l&rsquo;int\u00e9gralit\u00e9 du corpus, CLARA extrait de mani\u00e8re r\u00e9p\u00e9t\u00e9e de petits sous-\u00e9chantillons al\u00e9atoires repr\u00e9sentatifs (typiquement de taille <span class=\"arabpsy-math-inline\" data-latex=\"n_{\\text{sub}} = 40 + 2k\">$n_{\\text{sub}} = 40 + 2k$<\/span>). Sur chaque sous-\u00e9chantillon, l&rsquo;algorithme d\u00e9ploie un partitionnement robuste autour des m\u00e9do\u00efdes (PAM, <em>Partitioning Around Medoids<\/em>), une approche moins sensible aux outliers que les centro\u00efdes classiques. Les m\u00e9do\u00efdes optimaux extraits de ces it\u00e9rations sont ensuite utilis\u00e9s pour assigner l&rsquo;ensemble de la m\u00e9ga-cohorte globale selon la r\u00e8gle de la distance minimale.<\/p>\n<p>Le compromis r\u00e9side dans une l\u00e9g\u00e8re perte de finesse topologique au niveau des micro-branches locales en \u00e9change d&rsquo;une r\u00e9duction drastique de la complexit\u00e9 temporelle et spatiale, ramen\u00e9e \u00e0 une \u00e9chelle lin\u00e9aire <span class=\"arabpsy-math-inline\" data-latex=\"O(k \\cdot n_{\\text{sub}}^2 + k(N - k))\">$O(k \\cdot n_{\\text{sub}}^2 + k(N &#8211; k))$<\/span>, autorisant la segmentation fluide de cohortes comprenant plusieurs centaines de milliers d&rsquo;unit\u00e9s statistiques.<\/p>\n<h3>12.2 Approches hybrides : combinaison k-means et classification hi\u00e9rarchique<\/h3>\n<p>Pour concilier la rapidit\u00e9 computationnelle du partitionnement non hi\u00e9rarchique et la profondeur d&rsquo;interpr\u00e9tation visuelle de l&rsquo;arborescence dendrogrammique, les m\u00e9thodologistes recourent fr\u00e9quemment \u00e0 des architectures de classification hybrides. L&rsquo;une des strat\u00e9gies les plus f\u00e9condes consiste \u00e0 utiliser les <em>k-means<\/em> ou le partitionnement par m\u00e9do\u00efdes comme une phase de pr\u00e9traitement et de compression g\u00e9om\u00e9trique pr\u00e9liminaire.<\/p>\n<p>Dans ce protocole en deux temps, l&rsquo;analyste applique d&rsquo;abord un partitionnement <em>k-means<\/em> param\u00e9tr\u00e9 avec un nombre \u00e9lev\u00e9 de classes (par exemple <span class=\"arabpsy-math-inline\" data-latex=\"k = 100\">$k = 100$<\/span> pour un jeu de donn\u00e9es de <span class=\"arabpsy-math-inline\" data-latex=\"N = 20,000\">$N = 20,000$<\/span> individus). Cette premi\u00e8re \u00e9tape agr\u00e8ge rapidement les individus au sein de micro-clusters denses et compacts, r\u00e9duisant la taille du probl\u00e8me de $20,000$ observations \u00e9l\u00e9mentaires \u00e0 $100$ centro\u00efdes virtuels sans d\u00e9naturer la distribution sous-jacente du nuage de points.<\/p>\n<p>Dans un second temps, la classification ascendante hi\u00e9rarchique avec la m\u00e9thode de Ward est appliqu\u00e9e non plus sur les donn\u00e9es brutes individuelles, mais exclusivement sur la matrice des $100$ centro\u00efdes pond\u00e9r\u00e9s par les effectifs respectifs de leurs classes d&rsquo;origine. Cette d\u00e9marche est encapsul\u00e9e de mani\u00e8re totalement automatis\u00e9e au sein de la fonction <em>HCPC()<\/em> (<em>Hierarchical Clustering on Principal Components<\/em>) du package <em>FactoMineR<\/em> via l&rsquo;argument <em>kk<\/em> :<\/p>\n<p><code>res_hcpc &lt;- HCPC(res_pca, kk = 100, method = \"ward\")<\/code><\/p>\n<p>Cette synergie computationnelle \u00e9limine le bruit statistique microscopique, contourne les limitations de m\u00e9moire vive, et g\u00e9n\u00e8re un dendrogramme parfaitement tra\u00e7able et manipulable par le chercheur. Elle d\u00e9montre comment l&rsquo;ing\u00e9nierie algorithmique contemporaine sous R parvient \u00e0 transcender les clivages m\u00e9thodologiques traditionnels pour offrir des solutions scalables et robustes.<\/p>\n<h3>12.3 Bonnes pratiques, pi\u00e8ges r\u00e9currents et directives pour le chercheur<\/h3>\n<p>La conduite d&rsquo;une classification hi\u00e9rarchique sous R, si elle est d&rsquo;une grande accessibilit\u00e9 logicielle, n&rsquo;est pas exempte d&rsquo;\u00e9cueils m\u00e9thodologiques majeurs que tout chercheur rigoureux doit s&rsquo;efforcer de neutraliser. Le premier pi\u00e8ge r\u00e9side dans la r\u00e9ification abusive des regroupements : un algorithme de classification produit structurellement des clusters, y compris lorsqu&rsquo;il est aliment\u00e9 par des matrices de nombres purement al\u00e9atoires issus d&rsquo;un bruit blanc uniforme. Le devoir d\u00e9ontologique de l&rsquo;analyste consiste \u00e0 prouver que la partition finale capture une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 substantielle r\u00e9elle en s&rsquo;appuyant imp\u00e9rativement sur des crit\u00e8res de validation stochastiques (tels que <em>pvclust<\/em>, la Gap Statistic ou l&rsquo;indice de Rand ajust\u00e9) avant de promulguer l&rsquo;existence de typologies empiriques.<\/p>\n<p>Le second \u00e9cueil m\u00e9thodologique concerne l&rsquo;omission de la standardisation des variables ou l&rsquo;utilisation inappropri\u00e9e d&rsquo;une m\u00e9trique de distance en contradiction avec la nature d&rsquo;\u00e9chelle des donn\u00e9es. L&rsquo;application m\u00e9canique de la distance euclidienne sur des \u00e9chelles qualitatives nominales sans recours \u00e0 la m\u00e9trique de Gower, ou sur des variables continues non r\u00e9duites pr\u00e9sentant des variances asym\u00e9triques, invalide int\u00e9gralement les conclusions scientifiques des travaux.<\/p>\n<p>De m\u00eame, la tra\u00e7abilit\u00e9 et la reproductibilit\u00e9 computationnelle exigent la consignation minutieuse des environnements d&rsquo;ex\u00e9cution. Les variations incr\u00e9mentales entre versions de packages (comme les \u00e9volutions historiques de la m\u00e9thode <em>ward.D<\/em> vers <em>ward.D2<\/em> dans R) peuvent induire des divergences de regroupement sur des cas limites. L&rsquo;analyste doit documenter explicitement la version de R exploit\u00e9e, les versions des biblioth\u00e8ques mobilis\u00e9es (via la commande <em>sessionInfo()<\/em>), fixer scrupuleusement la graine al\u00e9atoire (<em>set.seed<\/em>) et mettre \u00e0 disposition des scripts d&rsquo;analyse modulaires et autonomes.<\/p>\n<p>Pour encadrer la valorisation scientifique des r\u00e9sultats typologiques, la grille d&rsquo;\u00e9valuation suivante formalise les exigences minimales d&rsquo;une publication acad\u00e9mique de haut niveau :<\/p>\n<ul>\n<li><strong>Justification formelle du pr\u00e9traitement :<\/strong> D\u00e9monstration explicite de l&rsquo;impact des \u00e9chelles de mesure, strat\u00e9gie d&rsquo;imputation document\u00e9e pour les donn\u00e9es manquantes et justification du choix de standardisation (z-score versus robuste).<\/li>\n<li><strong>Explicitation du choix m\u00e9trique :<\/strong> Argumentation th\u00e9orique reliant la distance choisie (euclidienne, Manhattan, Pearson, Gower) \u00e0 la morphologie attendue des construits \u00e9tudi\u00e9s.<\/li>\n<li><strong>D\u00e9tail transparent du m\u00e9canisme de liaison :<\/strong> Mention explicite de l&rsquo;algorithme d&rsquo;agr\u00e9gation retenu (notamment <em>ward.D2<\/em>) et contr\u00f4le des coefficients d&rsquo;agglom\u00e9ration (AC).<\/li>\n<li><strong>D\u00e9termination rationnelle du nombre <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> :<\/strong> Pr\u00e9sentation document\u00e9e d&rsquo;au moins deux m\u00e9triques de consensus objectif (silhouette, gap, indice de Davies-Bouldin, ou consensus <em>NbClust<\/em>) pour justifier le seuil d&rsquo;\u00e9lagage.<\/li>\n<li><strong>Preuves de stabilit\u00e9 et de validit\u00e9 :<\/strong> Rapport syst\u00e9matique des valeurs de significativit\u00e9 par r\u00e9\u00e9chantillonnage bootstrap (valeurs p AU de <em>pvclust<\/em>) et validation de la r\u00e9plicabilit\u00e9 sur sous-\u00e9chantillon (ARI).<\/li>\n<li><strong>Caract\u00e9risation multidimensionnelle compl\u00e8te :<\/strong> Fourniture des tableaux de moyennes, des \u00e9carts-types, des valeurs-tests et projections visuelles sur plans factoriels permettant au lecteur d&rsquo;\u00e9valuer la s\u00e9parabilit\u00e9 substantielle des groupes.<\/li>\n<\/ul>\n<p>En respectant rigoureusement ces pr\u00e9ceptes op\u00e9ratoires, le chercheur s&rsquo;assure que la classification ascendante hi\u00e9rarchique sous R ne se r\u00e9duit pas \u00e0 une heuristique descriptive subjective, mais d\u00e9ploie sa v\u00e9ritable puissance : celle d&rsquo;un instrument d&rsquo;investigation scientifique rigoureux, capable d&rsquo;isoler des typologies comportementales authentiques et de renouveler l&rsquo;\u00e9difice th\u00e9orique des disciplines explor\u00e9es.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Charrad, M., Ghazzali, N., Boiteau, V., &amp; Niknafs, A. (2014). NbClust: An R package for determining the relevant number of clusters in a data set. <em>Journal of Statistical Software<\/em>, 61(6), 1\u201336. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.18637\/jss.v061.i06\">https:\/\/doi.org\/10.18637\/jss.v061.i06<\/a><\/li>\n<li>Everitt, B. S., Landau, S., Leese, M., &amp; Stahl, D. (2011). <em>Cluster analysis<\/em> (5th ed.). John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/9780470977811\">https:\/\/doi.org\/10.1002\/9780470977811<\/a><\/li>\n<li>Galili, T. (2015). dendextend: an R package for visualizing, adjusting and comparing trees of hierarchical clustering. <em>Bioinformatics<\/em>, 31(22), 3718\u20133720. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/bioinformatics\/btv428\">https:\/\/doi.org\/10.1093\/bioinformatics\/btv428<\/a><\/li>\n<li>Gower, J. C. (1971). A general coefficient of similarity and some of its properties. <em>Biometrics<\/em>, 27(4), 857\u2013871. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2528823\">https:\/\/doi.org\/10.2307\/2528823<\/a><\/li>\n<li>Husson, F., L\u00ea, S., &amp; Pag\u00e8s, J. (2017). <em>Exploratory multivariate analysis by example using R<\/em> (2nd ed.). Chapman &amp; Hall\/CRC. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1201\/b21874\">https:\/\/doi.org\/10.1201\/b21874<\/a><\/li>\n<li>Kassambara, A. (2017). <em>Practical guide to cluster analysis in R: Unsupervised machine learning<\/em> (Vol. 1). STHDA.<\/li>\n<li>Kaufman, L., &amp; Rousseeuw, P. J. (1990). <em>Finding groups in data: An introduction to cluster analysis<\/em>. John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/9780470316801\">https:\/\/doi.org\/10.1002\/9780470316801<\/a><\/li>\n<li>Lance, G. N., &amp; Williams, W. T. (1967). A general theory of classificatory sorting strategies: 1. Hierarchical systems. <em>The Computer Journal<\/em>, 9(4), 373\u2013380. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/comjnl\/9.4.373\">https:\/\/doi.org\/10.1093\/comjnl\/9.4.373<\/a><\/li>\n<li>Murtagh, F., &amp; Legendre, P. (2014). Ward\u2019s hierarchical agglomerative clustering method: Which algorithms implement Ward\u2019s criterion? <em>Journal of Classification<\/em>, 31(3), 274\u2013295. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/s00357-014-9161-z\">https:\/\/doi.org\/10.1007\/s00357-014-9161-z<\/a><\/li>\n<li>Rousseeuw, P. J. (1987). Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. <em>Journal of Computational and Applied Mathematics<\/em>, 20, 53\u201365. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/0377-0427(87)90125-7\">https:\/\/doi.org\/10.1016\/0377-0427(87)90125-7<\/a><\/li>\n<li>Suzuki, R., &amp; Shimodaira, H. (2006). Pvclust: An R package for assessing the uncertainty in hierarchical clustering. <em>Bioinformatics<\/em>, 22(12), 1540\u20131542. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/bioinformatics\/btl117\">https:\/\/doi.org\/10.1093\/bioinformatics\/btl117<\/a><\/li>\n<li>Tibshirani, R., Walther, G., &amp; Hastie, T. (2001). Estimating the number of clusters in a data set via the gap statistic. <em>Journal of the Royal Statistical Society: Series B (Statistical Methodology)<\/em>, 63(2), 411\u2013423. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/1467-9868.00293\">https:\/\/doi.org\/10.1111\/1467-9868.00293<\/a><\/li>\n<li>Tukey, J. W. (1977). <em>Exploratory data analysis<\/em>. Addison-Wesley.<\/li>\n<li>Ward, J. H., Jr. (1963). Hierarchical grouping to optimize an objective function. <em>Journal of the American Statistical Association<\/em>, 58(301), 236\u2013244. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1963.10500845\">https:\/\/doi.org\/10.1080\/01621459.1963.10500845<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet sur la classification hi\u00e9rarchique dans R : th\u00e9orie math\u00e9matique, dendrogrammes, validation statistique et application pas \u00e0 pas.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2937,2269],"tags":[7457,7458,16,2658,2344],"class_list":["post-5549","psych_stat","type-psych_stat","status-publish","hentry","category-analyse-de-donnees","category-methodes-quantitatives","tag-classification-hierarchique","tag-dendrogramme","tag-psychometrie","tag-r","tag-statistiques"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/5549","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=5549"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/5549\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=5549"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=5549"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=5549"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}