{"id":4871,"date":"2026-09-11T13:56:59","date_gmt":"2026-09-11T13:56:59","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/introduction-analyse-discriminante-lineaire\/"},"modified":"2026-09-11T13:56:59","modified_gmt":"2026-09-11T13:56:59","slug":"introduction-analyse-discriminante-lineaire","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=introduction-analyse-discriminante-lineaire","title":{"rendered":"Introduction \u00e0 l&rsquo;analyse discriminante lin\u00e9aire"},"content":{"rendered":"<p>L&rsquo;analyse statistique multivari\u00e9e constitue le socle fondamental sur lequel reposent la mod\u00e9lisation pr\u00e9dictive contemporaine, la biom\u00e9trie computationnelle et la psychom\u00e9trie quantitative. Parmi l&rsquo;arsenal m\u00e9thodologique d\u00e9velopp\u00e9 au cours du XXe si\u00e8cle pour d\u00e9m\u00ealer les structures latentes au sein de donn\u00e9es multidimensionnelles, l&rsquo;analyse discriminante lin\u00e9aire \u2014 fr\u00e9quemment d\u00e9sign\u00e9e par le sigle ADL ou LDA pour <em>Linear Discriminant Analysis<\/em> \u2014 occupe une place \u00e9pist\u00e9mologique et pratique singuli\u00e8re. Con\u00e7ue initialement pour r\u00e9pondre \u00e0 des imp\u00e9ratifs de classification morphologique et taxonomique, cette m\u00e9thode combine avec une remarquable \u00e9l\u00e9gance la rigueur de l&rsquo;alg\u00e8bre lin\u00e9aire, la g\u00e9om\u00e9trie des projections vectorielles et la puissance inf\u00e9rentielle de la th\u00e9orie bay\u00e9sienne de la d\u00e9cision. Elle offre ainsi un cadre unifi\u00e9 permettant simultan\u00e9ment de r\u00e9duire la dimensionnalit\u00e9 d&rsquo;un ensemble complexe de pr\u00e9dicteurs continus et d&rsquo;assigner des entit\u00e9s \u00e0 des cat\u00e9gories discr\u00e8tes mutuellement exclusives.<\/p>\n<p>Dans le champ des sciences du comportement, de la neuropsychologie et de la psychiatrie computationnelle, l&rsquo;investigation empirique se heurte continuellement \u00e0 la haute dimensionnalit\u00e9 des batteries d&rsquo;\u00e9valuation. Qu&rsquo;il s&rsquo;agisse de profils cognitifs composites, de faisceaux de biomarqueurs neurophysiologiques ou de scores d&rsquo;\u00e9chelles psychopathologiques, le chercheur doit naviguer entre deux \u00e9cueils m\u00e9thodologiques majeurs : la perte d&rsquo;interpr\u00e9tabilit\u00e9 inh\u00e9rente aux mod\u00e8les opaques dits de \u00ab bo\u00eete noire \u00bb et l&rsquo;instabilit\u00e9 statistique induite par le fl\u00e9au de la dimensionnalit\u00e9. L&rsquo;analyse discriminante lin\u00e9aire r\u00e9sout ce dilemme avec une parcimonie math\u00e9matique exemplaire, en identifiant les combinaisons lin\u00e9aires de variables qui maximisent l&rsquo;\u00e9cartement inter-groupes tout en minimisant la dispersion intra-groupe. Ce faisant, elle ne se contente pas de proposer un algorithme de classification supervis\u00e9e hautement performant : elle r\u00e9v\u00e8le les axes fondamentaux selon lesquels des sous-populations cliniques ou comportementales se diff\u00e9rencient structurellement.<\/p>\n<p>Le pr\u00e9sent trait\u00e9 propose une exploration exhaustive, rigoureuse et didactique de l&rsquo;analyse discriminante lin\u00e9aire. Con\u00e7u \u00e0 l&rsquo;intention des chercheurs, des statisticiens et des praticiens de la mod\u00e9lisation des donn\u00e9es psychom\u00e9triques, cet article d\u00e9cortique les fondations th\u00e9oriques, les postulats probabilistes, les d\u00e9rivations alg\u00e9briques matricielles, ainsi que les protocoles de validation empirique qui r\u00e9gissent cette discipline. Du crit\u00e8re de s\u00e9paration de Fisher aux extensions r\u00e9gularis\u00e9es de haute dimension, en passant par les comparaisons formelles avec la r\u00e9gression logistique multinomiale et l&rsquo;analyse en composantes principales, ce document fournit les rep\u00e8res conceptuels et m\u00e9thodologiques indispensables pour ma\u00eetriser et d\u00e9ployer l&rsquo;ADL avec une rigueur irr\u00e9prochable.<\/p>\n<h2>1. Fondements conceptuels et \u00e9mergence de l&rsquo;analyse discriminante lin\u00e9aire<\/h2>\n<h3>1.1 Origines historiques et formulation par Ronald Fisher<\/h3>\n<p>L&rsquo;acte de naissance formel de l&rsquo;analyse discriminante lin\u00e9aire remonte \u00e0 l&rsquo;ann\u00e9e 1936, sous la plume du g\u00e9n\u00e9ticien et statisticien britannique <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Ronald_Aylmer_Fisher\">Ronald Aylmer Fisher<\/a>, dans son article s\u00e9minal intitul\u00e9 <em>The Use of Multiple Measurements in Taxonomic Problems<\/em> publi\u00e9 au sein des <em>Annals of Eugenics<\/em>. Confront\u00e9 aux donn\u00e9es morphom\u00e9triques collect\u00e9es par le botaniste Edgar Anderson sur trois esp\u00e8ces sympatriques d&rsquo;iris (<em>Iris setosa<\/em>, <em>Iris versicolor<\/em> et <em>Iris virginica<\/em>), Fisher cherchait \u00e0 r\u00e9soudre un probl\u00e8me fondamental de taxonomie : comment combiner quatre mesures ph\u00e9notypiques continues (la longueur et la largeur des s\u00e9pales, ainsi que la longueur et la largeur des p\u00e9tales) afin d&rsquo;\u00e9tablir une r\u00e8gle de d\u00e9marcation objective et reproductible capable d&rsquo;assigner sans \u00e9quivoque chaque sp\u00e9cimen v\u00e9g\u00e9tal \u00e0 son esp\u00e8ce botanique d&rsquo;origine.<\/p>\n<p>\u00c0 cette \u00e9poque, les statistiques appliqu\u00e9es demeuraient largement confin\u00e9es \u00e0 l&rsquo;analyse univari\u00e9e ou \u00e0 des d\u00e9marches descriptives fragmentaires. L&rsquo;innovation conceptuelle de Fisher a consist\u00e9 \u00e0 op\u00e9rer une transition paradigmatique majeure vers l&rsquo;analyse multivari\u00e9e inf\u00e9rentielle. Plut\u00f4t que d&rsquo;examiner chaque variable morphologique isol\u00e9ment \u2014 d\u00e9marche qui \u00e9chouait en raison des chevauchements consid\u00e9rables des distributions marginales \u2014, Fisher a postul\u00e9 l&rsquo;existence d&rsquo;un axe g\u00e9om\u00e9trique optimal. Cet axe correspond \u00e0 une combinaison lin\u00e9aire des variables pr\u00e9dictives le long de laquelle les centres de gravit\u00e9 des populations sont les plus distants possible, comparativement \u00e0 l&rsquo;\u00e9talement des donn\u00e9es autour de ces m\u00eames centres. Cette intuition a scell\u00e9 l&rsquo;\u00e9mergence d&rsquo;une mod\u00e9lisation conjointe de la variance et de la covariance dans l&rsquo;espace multidimensionnel.<\/p>\n<p>L&rsquo;approche initiale de Fisher \u00e9tait purement d\u00e9terministe et g\u00e9om\u00e9trique, reposant sur l&rsquo;optimisation d&rsquo;un rapport de variances sans postulat probabiliste explicite quant \u00e0 la distribution sous-jacente des populations. Il a fallu attendre les d\u00e9veloppements ult\u00e9rieurs de statisticiens tels que Harold Hotelling, P. C. Mahalanobis et, plus tard, C. R. Rao pour que cette formulation g\u00e9om\u00e9trique soit int\u00e9gr\u00e9e dans le cadre probabiliste de la th\u00e9orie de la d\u00e9cision et du formalisme bay\u00e9sien. Cette \u00e9volution a transform\u00e9 l&rsquo;outil de s\u00e9paration en un classificateur probabiliste \u00e0 part enti\u00e8re, capable de calculer des probabilit\u00e9s <em>a posteriori<\/em> d&rsquo;appartenance \u00e0 des groupes tout en pr\u00e9servant l&rsquo;\u00e9l\u00e9gance analytique de la solution originale de Fisher.<\/p>\n<h3>1.2 D\u00e9finition formelle et objectifs en psychom\u00e9trie<\/h3>\n<p>Dans le domaine de la psychom\u00e9trie et des sciences quantitatives du comportement, l&rsquo;analyse discriminante lin\u00e9aire se d\u00e9finit comme une m\u00e9thode statistique multivari\u00e9e supervis\u00e9e ayant pour vocation d&rsquo;\u00e9tudier les relations entre une variable d\u00e9pendante qualitative polytomique (repr\u00e9sentant <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> classes, groupes diagnostiques ou typologies comportementales distinctes) et un vecteur de <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> variables ind\u00e9pendantes quantitatives continues. Ces pr\u00e9dicteurs sont g\u00e9n\u00e9ralement issus de batteries de tests psychologiques, de questionnaires standardis\u00e9s, de bilans neuropsychologiques ou d&rsquo;enregistrements physiologiques. L&rsquo;objectif m\u00e9thodologique premier consiste \u00e0 d\u00e9terminer si ces sous-groupes d&rsquo;individus, pr\u00e9alablement d\u00e9finis sur des crit\u00e8res th\u00e9oriques ou cliniques ind\u00e9pendants, pr\u00e9sentent des profils distincts sur l&rsquo;ensemble de ces \u00e9preuves cognitives ou \u00e9motionnelles.<\/p>\n<p>Cette m\u00e9thodologie r\u00e9pond \u00e0 une double vocation th\u00e9orique et appliqu\u00e9e, indissociable de la recherche psychom\u00e9trique contemporaine :<\/p>\n<ul>\n<li><strong>La s\u00e9paration g\u00e9om\u00e9trique des groupes (vocation structurelle et descriptive) :<\/strong> Il s&rsquo;agit d&rsquo;identifier le nombre et la nature des dimensions sous-jacentes orthogonales selon lesquelles les populations se diff\u00e9rencient. L&rsquo;ADL condense l&rsquo;espace des descripteurs initiaux en un sous-espace discriminant de dimension r\u00e9duite, r\u00e9v\u00e9lant comment les fonctions psychologiques mesur\u00e9es s&rsquo;articulent pour isoler une trajectoire clinique d&rsquo;une trajectoire neurotypique.<\/li>\n<li><strong>La classification et la pr\u00e9diction d&rsquo;appartenance (vocation d\u00e9cisionnelle) :<\/strong> Le mod\u00e8le produit des r\u00e8gles math\u00e9matiques d&rsquo;assignation permettant d&rsquo;attribuer un nouvel individu, dont le profil psychom\u00e9trique est connu mais dont l&rsquo;appartenance de groupe est inconnue, \u00e0 l&rsquo;une des populations de r\u00e9f\u00e9rence, tout en quantifiant l&rsquo;incertitude li\u00e9e \u00e0 cette pr\u00e9diction.<\/li>\n<\/ul>\n<p>Cette approche s&rsquo;av\u00e8re particuli\u00e8rement f\u00e9conde pour l&rsquo;identification de signatures cognitives sp\u00e9cifiques. Par exemple, au sein des troubles neurod\u00e9veloppementaux, l&rsquo;ADL permet de d\u00e9terminer si le trouble du spectre de l&rsquo;autisme, le trouble du d\u00e9ficit de l&rsquo;attention avec ou sans hyperactivit\u00e9 (TDAH) et les troubles sp\u00e9cifiques du langage forment des constellations cognitives distinctes dans l&rsquo;espace des fonctions ex\u00e9cutives. En condensant les scores bruts en <em>scores discriminants individuels<\/em>, elle fournit aux cliniciens une \u00e9chelle standardis\u00e9e unique indiquant le positionnement pr\u00e9cis d&rsquo;un patient le long d&rsquo;un gradient de diff\u00e9renciation pathologique, facilitant ainsi le diagnostic diff\u00e9rentiel objectif.<\/p>\n<h3>1.3 Distinction entre inf\u00e9rence explicative et classification pr\u00e9dictive<\/h3>\n<p>Un aspect \u00e9pist\u00e9mologique fondamental de l&rsquo;analyse discriminante lin\u00e9aire r\u00e9side dans son statut charni\u00e8re entre la mod\u00e9lisation statistique inf\u00e9rentielle et l&rsquo;apprentissage automatique (<em>machine learning<\/em>). Contrairement \u00e0 de nombreux algorithmes modernes purement pr\u00e9dictifs \u2014 tels que les r\u00e9seaux de neurones profonds ou les for\u00eats al\u00e9atoires, qui fonctionnent souvent comme des estimateurs non param\u00e9triques opaques \u2014, l&rsquo;ADL maintient un \u00e9quilibre structurel entre explicabilit\u00e9 des coefficients et minimisation du taux d&rsquo;erreur de classification.<\/p>\n<p>Sur le plan explicatif, l&rsquo;analyste d\u00e9cortique les coefficients des fonctions discriminantes canoniques, assimilables aux coefficients de r\u00e9gression ou aux saturations factorielles. L&rsquo;\u00e9valuation de ces poids statistiques permet d&rsquo;ordonner les pr\u00e9dicteurs selon leur contribution relative \u00e0 la d\u00e9marcation des groupes, neutralisant les intercorr\u00e9lations entre variables. L&rsquo;analyste peut ainsi affirmer, par exemple, que dans la s\u00e9paration entre un vieillissement cognitif sain et une phase prodromique de d\u00e9mence, la composante de rappel diff\u00e9r\u00e9 de la m\u00e9moire \u00e9pisodique exerce un pouvoir discriminant trois fois sup\u00e9rieur \u00e0 celui de la vitesse de traitement de l&rsquo;information, m\u00eame si les deux tests pr\u00e9sentent une corr\u00e9lation bivari\u00e9e \u00e9lev\u00e9e.<\/p>\n<p>Sur le versant pr\u00e9dictif, l&rsquo;accent est mis sur l&rsquo;optimisation des hyperplans s\u00e9parateurs afin de minimiser le risque global de mauvaise classification sur de futures observations. Cet arbitrage entre inf\u00e9rence et pr\u00e9diction s&rsquo;inscrit au c\u0153ur du compromis biais-variance. L&rsquo;ADL impose une structure param\u00e9trique rigide (lin\u00e9arit\u00e9 des fronti\u00e8res d\u00e9cisionnelles), ce qui induit un biais de mod\u00e8le sup\u00e9rieur \u00e0 celui des classificateurs non lin\u00e9aires hautement flexibles. En contrepartie, cette parcimonie conf\u00e8re \u00e0 l&rsquo;ADL une variance d&rsquo;\u00e9chantillonnage exceptionnellement faible, pr\u00e9venant le surapprentissage (<em>overfitting<\/em>) sur des donn\u00e9es psychom\u00e9triques o\u00f9 la taille des cohortes est structurellement restreinte par des contraintes de recrutement clinique.<\/p>\n<p>Il importe enfin de distinguer rigoureusement l&rsquo;ADL de l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Analyse_en_composantes_principales\">analyse en composantes principales (ACP)<\/a>. Bien que toutes deux proc\u00e8dent par transformation lin\u00e9aire et r\u00e9duction de dimensionnalit\u00e9 via la diagonalisation matricielle, l&rsquo;ACP est intrins\u00e8quement <em>non supervis\u00e9e<\/em> : elle maximise la variance globale totale des donn\u00e9es sans \u00e9gard pour une quelconque labellisation des observations. \u00c0 l&rsquo;inverse, l&rsquo;ADL est <em>supervis\u00e9e<\/em> : elle recherche sp\u00e9cifiquement la variance inter-classes en la calibrant sur la dispersion r\u00e9siduelle intra-classe. Il est fr\u00e9quent que la premi\u00e8re composante principale d&rsquo;une batterie de tests capture une variance g\u00e9n\u00e9rale non informative (telle que le facteur g\u00e9n\u00e9ral d&rsquo;intelligence ou le niveau socio-\u00e9ducatif), tandis que l&rsquo;axe discriminant de Fisher isole pr\u00e9cis\u00e9ment la variation pathognomonique subtile permettant d&rsquo;individualiser une pathologie.<\/p>\n<h2>2. L&rsquo;analyse discriminante lin\u00e9aire face \u00e0 la r\u00e9gression logistique<\/h2>\n<h3>2.1 Limites de la r\u00e9gression logistique pour les variables multiclasses<\/h3>\n<p>Dans la litt\u00e9rature biom\u00e9dicale et psychosociale, la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/R%C3%A9gression_logistique\">r\u00e9gression logistique<\/a> est fr\u00e9quemment \u00e9rig\u00e9e en m\u00e9thode de r\u00e9f\u00e9rence pour mod\u00e9liser des issues cat\u00e9gorielles binaires. Cependant, d\u00e8s lors que la variable d\u00e9pendante comporte trois modalit\u00e9s ou plus (<span class=\"arabpsy-math-inline\" data-latex=\"k ge 3\">$k ge 3$<\/span>) \u2014 configuration omnipr\u00e9sente en clinique o\u00f9 l&rsquo;on oppose par exemple des sous-types cliniques (trouble d\u00e9pressif caract\u00e9ris\u00e9, trouble bipolaire en phase d\u00e9pressive, trouble schizo-affectif et t\u00e9moins sains) \u2014, la r\u00e9gression logistique multinomiale r\u00e9v\u00e8le des contraintes structurelles et computationnelles majeures.<\/p>\n<p>Le mod\u00e8le logistique multinomial impose le choix arbitraire d&rsquo;une cat\u00e9gorie de r\u00e9f\u00e9rence par rapport \u00e0 laquelle <span class=\"arabpsy-math-inline\" data-latex=\"k - 1\">$k &#8211; 1$<\/span> \u00e9quations logit distinctes doivent \u00eatre estim\u00e9es simultan\u00e9ment. Cette sp\u00e9cification asym\u00e9trique entra\u00eene une prolif\u00e9ration exponentielle du nombre de param\u00e8tres libres : pour <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> pr\u00e9dicteurs et <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> classes, le mod\u00e8le requiert l&rsquo;estimation de <span class=\"arabpsy-math-inline\" data-latex=\"(k - 1) \\times (p + 1)\">$(k &#8211; 1) \\times (p + 1)$<\/span> coefficients. Cette complexit\u00e9 param\u00e9trique accrue d\u00e9t\u00e9riore rapidement la puissance statistique d\u00e8s lors que le ratio entre le nombre de sujets par groupe et le nombre de variables diminue. L&rsquo;estimation, reposant sur l&rsquo;optimisation it\u00e9rative de la fonction de vraisemblance via l&rsquo;algorithme de Newton-Raphson ou des m\u00e9thodes de quasi-Newton, devient sujette \u00e0 des instabilit\u00e9s d&rsquo;optimisation num\u00e9rique, telles que la non-convergence ou l&rsquo;arr\u00eat sur des points d&rsquo;inflexion locaux en pr\u00e9sence de corr\u00e9lations mod\u00e9r\u00e9es.<\/p>\n<p>L&rsquo;analyse discriminante lin\u00e9aire contourne ces \u00e9cueils algorithmiques gr\u00e2ce \u00e0 une formulation analytique ferm\u00e9e directe. Aucune it\u00e9ration num\u00e9rique n&rsquo;est requise pour converger vers la solution optimale : les param\u00e8tres du mod\u00e8le sont obtenus directement par des op\u00e9rations matricielles exactes d&rsquo;inversion et de d\u00e9composition spectrale. De surcro\u00eet, l&rsquo;ADL traite l&rsquo;ensemble des <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> groupes de fa\u00e7on rigoureusement sym\u00e9trique au sein d&rsquo;un m\u00eame espace vectoriel, \u00e9vitant ainsi le recours \u00e0 une modalit\u00e9 de r\u00e9f\u00e9rence artificielle et fournissant une g\u00e9om\u00e9trie des fronti\u00e8res de d\u00e9cision intrins\u00e8quement plus parcimonieuse et \u00e9l\u00e9gante.<\/p>\n<h3>2.2 Stabilit\u00e9 computationnelle et pr\u00e9cision sur petits \u00e9chantillons<\/h3>\n<p>La recherche clinique, neuropsychologique et d&rsquo;imagerie c\u00e9r\u00e9brale fonctionnelle se caract\u00e9rise presque invariablement par des tailles d&rsquo;\u00e9chantillons modestes, dict\u00e9es par le co\u00fbt financier, la lourdeur des protocoles d&rsquo;acquisition ou la raret\u00e9 \u00e9pid\u00e9miologique des affections cibl\u00e9es. Dans ces contextes o\u00f9 la taille totale de cohorte <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span> ne d\u00e9passe que rarement quelques dizaines ou centaines de participants pour un vecteur de mesures \u00e9tendu, le comportement asymptotique de l&rsquo;ADL surpasse syst\u00e9matiquement celui de la r\u00e9gression logistique.<\/p>\n<p>Ce ph\u00e9nom\u00e8ne s&rsquo;explique par la nature des paradigmes d&rsquo;estimation :<\/p>\n<ul>\n<li><strong>R\u00e9gression logistique (mod\u00e8le discriminatif conditionnel) :<\/strong> Elle mod\u00e9lise uniquement la distribution conditionnelle <span class=\"arabpsy-math-inline\" data-latex=\"P(Y|X)\">$P(Y|X)$<\/span>, ignorant d\u00e9lib\u00e9r\u00e9ment la distribution marginale <span class=\"arabpsy-math-inline\" data-latex=\"P(X)\">$P(X)$<\/span> des pr\u00e9dicteurs. Ce choix \u00e9vite d&rsquo;\u00e9mettre des hypoth\u00e8ses distributionnelles sur les variables pr\u00e9dictives, mais sacrifie une quantit\u00e9 substantielle d&rsquo;information statistique contenue dans la structure spatiale globale des donn\u00e9es observ\u00e9es.<\/li>\n<li><strong>Analyse discriminante lin\u00e9aire (mod\u00e8le g\u00e9n\u00e9ratif conjoint) :<\/strong> Elle mod\u00e9lise la distribution conjointe <span class=\"arabpsy-math-inline\" data-latex=\"P(X, Y) = P(X|Y)P(Y)\">$P(X, Y) = P(X|Y)P(Y)$<\/span>. En supposant que les pr\u00e9dicteurs suivent une distribution gaussienne multivari\u00e9e au sein de chaque classe, l&rsquo;ADL exploite la totalit\u00e9 des informations g\u00e9om\u00e9triques des centro\u00efdes et des matrices de dispersion empiriques.<\/li>\n<\/ul>\n<p>Conform\u00e9ment aux travaux th\u00e9oriques classiques d&rsquo;Efron (1975), lorsque les postulats de normalit\u00e9 multivari\u00e9e et d&rsquo;homosc\u00e9dasticit\u00e9 sont raisonnablement respect\u00e9s, l&rsquo;ADL est asymptotiquement plus efficiente que la r\u00e9gression logistique. Cela signifie qu&rsquo;elle atteint une erreur d&rsquo;estimation plus faible et des intervalles de confiance plus \u00e9troits pour une m\u00eame taille d&rsquo;\u00e9chantillon. En pratique, sur de petites cohortes, l&rsquo;ADL pr\u00e9sente une variance d&rsquo;\u00e9chantillonnage remarquablement r\u00e9duite, l\u00e0 o\u00f9 la r\u00e9gression logistique souffre d&rsquo;un surajustement prononc\u00e9 de ses coefficients, conduisant \u00e0 des pr\u00e9dictions erratiques lors de l&rsquo;application du mod\u00e8le \u00e0 de nouveaux patients ind\u00e9pendants.<\/p>\n<h3>2.3 Le probl\u00e8me de la s\u00e9paration parfaite<\/h3>\n<p>L&rsquo;une des vuln\u00e9rabilit\u00e9s m\u00e9thodologiques les plus handicapantes de la r\u00e9gression logistique r\u00e9side dans ce que la litt\u00e9rature statistique nomme le probl\u00e8me de la <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Separation_(statistics)\">s\u00e9paration quasi-compl\u00e8te ou parfaite<\/a> (ph\u00e9nom\u00e8ne identifi\u00e9 formellement par Albert et Anderson en 1984). Ce ph\u00e9nom\u00e8ne survient lorsqu&rsquo;un pr\u00e9dicteur univari\u00e9 ou une combinaison lin\u00e9aire de plusieurs variables permet de scinder impeccablement deux cat\u00e9gories sans la moindre zone de chevauchement. Dans ce cas de figure, qui correspond intuitivement \u00e0 une discrimination id\u00e9ale, le maximum de vraisemblance n&rsquo;existe pas au sens strict :<\/p>\n<ul>\n<li>Les coefficients de r\u00e9gression logistique associ\u00e9s tendent math\u00e9matiquement vers l&rsquo;infini (<span class=\"arabpsy-math-inline\" data-latex=\"+\\infty\">$+\\infty$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"-\\infty\">$-\\infty$<\/span>).<\/li>\n<li>Les erreurs-types asymptotiques explosent, atteignant des valeurs aberrantes de plusieurs milliers d&rsquo;unit\u00e9s.<\/li>\n<li>Les tests d&rsquo;hypoth\u00e8ses conventionnels (tests de Wald) deviennent caducs, concluant paradoxalement \u00e0 la non-significativit\u00e9 des pr\u00e9dicteurs les plus puissants du mod\u00e8le.<\/li>\n<li>Les algorithmes d&rsquo;optimisation num\u00e9rique \u00e9chouent \u00e0 converger, renvoyant des messages d&rsquo;avertissement critiques quant \u00e0 l&rsquo;instabilit\u00e9 de la matrice hessienne.<\/li>\n<\/ul>\n<p>L&rsquo;analyse discriminante lin\u00e9aire est totalement insensible \u00e0 cette pathologie computationnelle. D\u00e8s lors que deux groupes sont parfaitement disjoints dans l&rsquo;espace des pr\u00e9dicteurs, l&rsquo;ADL calcule sans difficult\u00e9 les centro\u00efdes respectifs <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_1\">$boldsymbol{\\mu}_1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_2\">$boldsymbol{\\mu}_2$<\/span>, d\u00e9termine la matrice de covariance combin\u00e9e intra-classe <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span>, et \u00e9tablit l&rsquo;hyperplan s\u00e9parateur optimal par une inversion matricielle standard. Les coefficients des fonctions discriminantes demeurent parfaitement stables, born\u00e9s et imm\u00e9diatement interpr\u00e9tables. L&rsquo;ADL constitue par cons\u00e9quent la m\u00e9thode d&rsquo;\u00e9lection face \u00e0 des variables psychom\u00e9triques hautement r\u00e9solutives g\u00e9n\u00e9rant des fronti\u00e8res de s\u00e9paration franches entre populations cliniques et populations t\u00e9moins.<\/p>\n<h2>3. Les postulats th\u00e9oriques et conditions d&rsquo;application<\/h2>\n<h3>3.1 L&rsquo;hypoth\u00e8se de normalit\u00e9 multivari\u00e9e<\/h3>\n<p>Pour que la r\u00e8gle de classification d\u00e9riv\u00e9e de l&rsquo;ADL soit th\u00e9oriquement optimale au sens de la minimisation de l&rsquo;erreur bay\u00e9sienne, le mod\u00e8le repose sur le postulat d&rsquo;une distribution normale multivari\u00e9e des pr\u00e9dicteurs au sein de chacune des populations parentes. Formellement, pour chaque classe <span class=\"arabpsy-math-inline\" data-latex=\"k in {1, dots, K}\">$k in {1, dots, K}$<\/span>, le vecteur al\u00e9atoire des pr\u00e9dicteurs <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{X} = (X_1, X_2, dots, X_p)^T\">$\\mathbf{X} = (X_1, X_2, dots, X_p)^T$<\/span> doit suivre une loi gaussienne \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> dimensions not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{N}_p(boldsymbol{\\mu}_k, boldsymbol{\\Sigma})\">$\\mathcal{N}_p(boldsymbol{\\mu}_k, boldsymbol{\\Sigma})$<\/span>.<\/p>\n<p>Ce postulat implique que non seulement chaque variable individuelle pr\u00e9sente une distribution normale univari\u00e9e, mais \u00e9galement que toutes les combinaisons lin\u00e9aires possibles de ces variables suivent des distributions normales univari\u00e9es, avec des surfaces de niveau de densit\u00e9 ellipso\u00efdales dans l&rsquo;espace <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^p\">$\\mathbb{R}^p$<\/span>. Dans la r\u00e9alit\u00e9 de l&rsquo;\u00e9valuation psychom\u00e9trique, ce postulat est fr\u00e9quemment mis \u00e0 l&rsquo;\u00e9preuve par des donn\u00e9es pr\u00e9sentant une forte asym\u00e9trie (<em>skewness<\/em>) ou un aplatissement marqu\u00e9 (<em>kurtosis<\/em>), r\u00e9sultant d&rsquo;effets de plancher ou de plafond sur des tests de performance cognitive, ou de distributions hautement biais\u00e9es sur des \u00e9chelles \u00e9valuant des comportements d\u00e9viants rares.<\/p>\n<p>L&rsquo;impact de la violation de ce postulat d\u00e9pend de son amplitude :<\/p>\n<ul>\n<li><strong>D\u00e9viations mineures \u00e0 mod\u00e9r\u00e9es :<\/strong> L&rsquo;ADL fait preuve d&rsquo;une robustesse substantielle. Fisher l&rsquo;ayant formul\u00e9e sans hypoth\u00e8se distributionnelle, la projection g\u00e9om\u00e9trique demeure op\u00e9rante m\u00eame si les densit\u00e9s ne sont pas strictement gaussiennes.<\/li>\n<li><strong>D\u00e9viations extr\u00eames ou pr\u00e9sence d&rsquo;asym\u00e9tries divergentes :<\/strong> Si une classe pr\u00e9sente une asym\u00e9trie positive tandis qu&rsquo;une autre pr\u00e9sente une asym\u00e9trie n\u00e9gative, les centro\u00efdes empiriques et l&rsquo;hyperplan lin\u00e9aire sont syst\u00e9matiquement d\u00e9cal\u00e9s, augmentant les taux de mauvaise classification au d\u00e9triment des sous-groupes minoritaires.<\/li>\n<\/ul>\n<p>L&rsquo;\u00e9valuation rigoureuse de la normalit\u00e9 multivari\u00e9e requiert le recours \u00e0 des instruments diagnostiques sp\u00e9cialis\u00e9s. Les coefficients d&rsquo;asym\u00e9trie et de kurtosis multivari\u00e9s de Mardia constituent la r\u00e9f\u00e9rence inf\u00e9rentielle : un test statistiquement significatif indique un rejet de l&rsquo;hypoth\u00e8se de gaussianit\u00e9 conjointe. Ces statistiques doivent \u00eatre compl\u00e9t\u00e9es par des approches graphiques, notamment le trac\u00e9 des diagrammes quantiles-quantiles (Q-Q plots) comparant le carr\u00e9 des <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Distance_de_Mahalanobis\">distances de Mahalanobis<\/a> des individus par rapport au centre de gravit\u00e9 de leur classe \u00e0 une distribution th\u00e9orique du Chi-deux \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> degr\u00e9s de libert\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"\\chi^2_p\">$\\chi^2_p$<\/span>). En cas d&rsquo;\u00e9cart manifeste, l&rsquo;application de transformations math\u00e9matiques univari\u00e9es pr\u00e9liminaires \u2014 telles que la famille de transformations de Box-Cox ou les transformations logarithmiques et racines carr\u00e9es \u2014 permet fr\u00e9quemment de restaurer la normalit\u00e9 multivari\u00e9e avant d&rsquo;ajuster le mod\u00e8le discriminant.<\/p>\n<h3>3.2 L&rsquo;homog\u00e9n\u00e9it\u00e9 des matrices de covariance<\/h3>\n<p>Le second pilier axiomatique de l&rsquo;analyse discriminante lin\u00e9aire, indispensable pour garantir la parfaite lin\u00e9arit\u00e9 des fronti\u00e8res de d\u00e9cision, r\u00e9side dans l&rsquo;hypoth\u00e8se d&rsquo;homog\u00e9n\u00e9it\u00e9 des matrices de variance-covariance inter-groupes, couramment d\u00e9sign\u00e9e sous le terme d&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Homosc%C3%A9dasticit%C3%A9\">homosc\u00e9dasticit\u00e9 multivari\u00e9e<\/a>. Ce postulat stipule que bien que les groupes poss\u00e8dent des vecteurs de moyennes distincts (<span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_1 ne boldsymbol{\\mu}_2 ne dots ne boldsymbol{\\mu}_K\">$boldsymbol{\\mu}_1 ne boldsymbol{\\mu}_2 ne dots ne boldsymbol{\\mu}_K$<\/span>), la structure de dispersion, de variance interne et de covariation entre les variables pr\u00e9dictives demeure strictement identique au sein de chaque population parente :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"boldsymbol{\\Sigma}_1 = boldsymbol{\\Sigma}_2 = dots = boldsymbol{\\Sigma}_K = boldsymbol{\\Sigma}\">$$boldsymbol{\\Sigma}_1 = boldsymbol{\\Sigma}_2 = dots = boldsymbol{\\Sigma}_K = boldsymbol{\\Sigma}$$<\/span><\/p>\n<p>Cette hypoth\u00e8se permet d&rsquo;estimer une matrice unique de covariance intra-groupe combin\u00e9e (<em>pooled covariance matrix<\/em>), not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}\">$\\mathbf{S}_{\\text{intra}}$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}_{\\text{pooled}}\">$boldsymbol{\\Sigma}_{\\text{pooled}}$<\/span>, obtenue par la moyenne pond\u00e9r\u00e9e des matrices empiriques de variance-covariance de chaque classe. Sur le plan g\u00e9om\u00e9trique, l&rsquo;homosc\u00e9dasticit\u00e9 garantit que les ellipso\u00efdes de dispersion repr\u00e9sentant les densit\u00e9s de probabilit\u00e9 de chaque classe poss\u00e8dent exactement la m\u00eame taille, la m\u00eame forme g\u00e9om\u00e9trique et la m\u00eame orientation spatiale dans l&rsquo;espace factoriel. Les groupes ne diff\u00e8rent d\u00e8s lors que par la translation de leur centre de gravit\u00e9 respectif.<\/p>\n<p>Lorsque cette condition n&rsquo;est pas remplie \u2014 c&rsquo;est-\u00e0-dire en cas d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 structurelle \u2014, les surfaces de densit\u00e9s conditionnelles \u00e9quivalentes ne se coupent plus selon des hyperplans plats, mais selon des hyperquadriques (courbures paraboliques, hyperboliques ou ellipso\u00efdales). Le maintien forc\u00e9 d&rsquo;une fronti\u00e8re lin\u00e9aire simplificatrice engendre alors une distorsion asym\u00e9trique de la r\u00e8gle d\u00e9cisionnelle : l&rsquo;hyperplan tend \u00e0 \u00eatre repouss\u00e9 vers le groupe pr\u00e9sentant la plus forte dispersion interne, provoquant une augmentation artificielle du taux de faux n\u00e9gatifs dans la population la plus variable. En psychologie clinique, cette situation est fr\u00e9quente : les cohortes saines pr\u00e9sentent typiquement une variance cognitive homog\u00e8ne et resserr\u00e9e, tandis que les groupes pathologiques manifestent une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 d&rsquo;atteinte consid\u00e9rable g\u00e9n\u00e9rant des variances substantiellement plus amples.<\/p>\n<p>Le contr\u00f4le formel de ce postulat s&rsquo;op\u00e8re canoniquement par le test <span class=\"arabpsy-math-inline\" data-latex=\"M\">$M$<\/span> de Box. Ce test inf\u00e9rentiel \u00e9value l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9galit\u00e9 conjointe des matrices de variance-covariance :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"H_0 : boldsymbol{\\Sigma}_1 = boldsymbol{\\Sigma}_2 = dots = boldsymbol{\\Sigma}_K\">$$H_0 : boldsymbol{\\Sigma}_1 = boldsymbol{\\Sigma}_2 = dots = boldsymbol{\\Sigma}_K$$<\/span><\/p>\n<p>Toutefois, le test <span class=\"arabpsy-math-inline\" data-latex=\"M\">$M$<\/span> de Box souffre d&rsquo;une sensibilit\u00e9 extr\u00eame \u00e0 la non-normalit\u00e9 des donn\u00e9es : une violation mineure de la normalit\u00e9 multivari\u00e9e suffit \u00e0 provoquer le rejet de l&rsquo;hypoth\u00e8se d&rsquo;homosc\u00e9dasticit\u00e9 \u00e0 tort. D\u00e8s lors, il est d&rsquo;usage en m\u00e9thodologie quantitative d&rsquo;adopter un seuil de significativit\u00e9 conservateur (par exemple <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0{,}001\">$\\alpha = 0{,}001$<\/span>) et d&rsquo;examiner attentivement les d\u00e9terminants logarithmiques des matrices individuelles <span class=\"arabpsy-math-inline\" data-latex=\"|\\mathbf{S}_k|\">$|\\mathbf{S}_k|$<\/span>. Si l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 est substantielle et irr\u00e9ductible, l&rsquo;analyste doit imp\u00e9rativement renoncer \u00e0 la lin\u00e9arit\u00e9 pour s&rsquo;orienter vers une analyse discriminante quadratique (QDA).<\/p>\n<h3>3.3 Ind\u00e9pendance et absence de multicolin\u00e9arit\u00e9 extr\u00eame<\/h3>\n<p>L&rsquo;analyse discriminante lin\u00e9aire s&rsquo;articule autour de deux hypoth\u00e8ses structurelles r\u00e9gissant la relation entre les unit\u00e9s statistiques et les pr\u00e9dicteurs :<\/p>\n<ul>\n<li><strong>L&rsquo;ind\u00e9pendance stricte des observations :<\/strong> Chaque sujet inclus dans l&rsquo;\u00e9chantillon doit avoir \u00e9t\u00e9 s\u00e9lectionn\u00e9 de mani\u00e8re ind\u00e9pendante, sans structure de grappes non mod\u00e9lis\u00e9e, sans mesures r\u00e9p\u00e9t\u00e9es longitudinales directes sur le m\u00eame individu au sein du m\u00eame jeu d&rsquo;apprentissage, et sans d\u00e9pendance spatio-temporelle. Les corr\u00e9lations r\u00e9siduelles intra-sujets violent le calcul des degr\u00e9s de libert\u00e9 et biaisent substantiellement les tests d&rsquo;inf\u00e9rence globale tels que le lambda de Wilks.<\/li>\n<li><strong>L&rsquo;absence de multicolin\u00e9arit\u00e9 parfaite ou extr\u00eame :<\/strong> L&rsquo;ADL requiert l&rsquo;inversion analytique de la matrice de dispersion intra-classe combin\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}\">$\\mathbf{S}_{\\text{intra}}$<\/span>. Si une variable psychom\u00e9trique est une combinaison lin\u00e9aire quasi-parfaite d&rsquo;autres \u00e9preuves de la batterie (par exemple, si le score total d&rsquo;un test est inclus simultan\u00e9ment avec l&rsquo;ensemble de ses sous-scores \u00e9l\u00e9mentaires), la matrice devient singuli\u00e8re ou mal conditionn\u00e9e.<\/li>\n<\/ul>\n<p>Sur le plan math\u00e9matique, une singularit\u00e9 matricielle se traduit par un d\u00e9terminant nul (<span class=\"arabpsy-math-inline\" data-latex=\"|\\mathbf{S}_{\\text{intra}}| = 0\">$|\\mathbf{S}_{\\text{intra}}| = 0$<\/span>), rendant l&rsquo;inversion matricielle rigoureusement impossible. Dans les cas de multicolin\u00e9arit\u00e9 extr\u00eame non parfaite mais s\u00e9v\u00e8re, l&rsquo;inversion num\u00e9rique produit des coefficients discriminants instables, caract\u00e9ris\u00e9s par des variances d&rsquo;\u00e9chantillonnage d\u00e9mesur\u00e9es : des variations imperceptibles dans le jeu de donn\u00e9es d&rsquo;apprentissage peuvent inverser radicalement le signe et l&rsquo;amplitude des poids discriminants canoniques.<\/p>\n<p>Avant d&rsquo;ajuster une analyse discriminante, le statisticien doit imp\u00e9rativement mener un diagnostic de colin\u00e9arit\u00e9 multivari\u00e9. Deux indicateurs statistiques s&rsquo;imposent :<\/p>\n<ul>\n<li><strong>Le facteur d&rsquo;inflation de la variance (VIF &#8211; <em>Variance Inflation Factor<\/em>) :<\/strong> Calcul\u00e9 pour chaque pr\u00e9dicteur r\u00e9gress\u00e9 sur tous les autres. Des valeurs de VIF exc\u00e9dant un seuil conventionnel de 5 ou 10 signalent une redondance critique n\u00e9cessitant l&rsquo;\u00e9lagage des variables redondantes.<\/li>\n<li><strong>L&rsquo;indice de conditionnement (<em>Condition Index<\/em>) de la matrice de donn\u00e9es :<\/strong> Un indice sup\u00e9rieur \u00e0 30, associ\u00e9 \u00e0 des proportions de variance \u00e9lev\u00e9es partag\u00e9es sur deux variables ou plus, d\u00e9montre une instabilit\u00e9 num\u00e9rique de la d\u00e9composition spectrale. L&rsquo;\u00e9limination des scores composites d\u00e9riv\u00e9s, la s\u00e9lection guid\u00e9e par la th\u00e9orie psychom\u00e9trique ou l&rsquo;application pr\u00e9alable d&rsquo;une extraction de facteurs orthogonaux constituent les rem\u00e8des m\u00e9thodologiques appropri\u00e9s.<\/li>\n<\/ul>\n<h2>4. Cadre math\u00e9matique et mod\u00e9lisation probabiliste bay\u00e9sienne<\/h2>\n<h3>4.1 Formulation probabiliste bay\u00e9sienne<\/h3>\n<p>L&rsquo;analyse discriminante lin\u00e9aire contemporaine s&rsquo;enracine rigoureusement dans la th\u00e9orie statistique de la d\u00e9cision formalis\u00e9e par la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Inference_bayesienne\">r\u00e8gle d&rsquo;inversion bay\u00e9sienne<\/a>. Dans cette perspective, la t\u00e2che de classification consiste \u00e0 d\u00e9duire la probabilit\u00e9 conditionnelle qu&rsquo;un individu pr\u00e9sentant un profil psychom\u00e9trique pr\u00e9cis <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x} in \\mathbb{R}^p\">$\\mathbf{x} in \\mathbb{R}^p$<\/span> appartienne \u00e0 la cat\u00e9gorie clinique <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"P(Y = k | \\mathbf{X} = \\mathbf{x})\">$P(Y = k | \\mathbf{X} = \\mathbf{x})$<\/span>. Cette quantit\u00e9 repr\u00e9sente la <em>probabilit\u00e9 a posteriori<\/em> de la classe.<\/p>\n<p>En vertu du th\u00e9or\u00e8me de Bayes, cette probabilit\u00e9 s&rsquo;\u00e9nonce comme le produit normalis\u00e9 de la vraisemblance conditionnelle de la classe par sa probabilit\u00e9 d&rsquo;occurrence globale dans la population parente :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"P(Y = k | \\mathbf{X} = \\mathbf{x}) = \\frac{f_k(\\mathbf{x}) , \\pi_k}{\\sum_{j=1}^{K} f_j(\\mathbf{x}) , \\pi_j}\">$$P(Y = k | \\mathbf{X} = \\mathbf{x}) = \\frac{f_k(\\mathbf{x}) , \\pi_k}{\\sum_{j=1}^{K} f_j(\\mathbf{x}) , \\pi_j}$$<\/span><\/p>\n<p>Cette d\u00e9composition structurelle fait intervenir trois termes d\u00e9cisionnels fondamentaux :<\/p>\n<ul>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"\\pi_k = P(Y = k)\">$\\pi_k = P(Y = k)$<\/span> repr\u00e9sente la <strong>probabilit\u00e9 a priori<\/strong> d&rsquo;appartenir \u00e0 la classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>. Elle quantifie le niveau de croyance ou la pr\u00e9valence \u00e9pid\u00e9miologique objective de l&rsquo;\u00e9tat pathologique au sein de la population g\u00e9n\u00e9rale avant toute administration de test psychom\u00e9trique.<\/li>\n<li><span class=\"arabpsy-math-inline\" data-latex=\"f_k(\\mathbf{x}) = P(\\mathbf{X} = \\mathbf{x} | Y = k)\">$f_k(\\mathbf{x}) = P(\\mathbf{X} = \\mathbf{x} | Y = k)$<\/span> d\u00e9signe la <strong>densit\u00e9 de probabilit\u00e9 conditionnelle<\/strong> de classe (la vraisemblance). Elle \u00e9value la densit\u00e9 de probabilit\u00e9 d&rsquo;observer pr\u00e9cis\u00e9ment la configuration de scores <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}\">$\\mathbf{x}$<\/span> chez un patient issu de la population clinique <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>.<\/li>\n<li>Le d\u00e9nominateur <span class=\"arabpsy-math-inline\" data-latex=\"\\sum_{j=1}^{K} f_j(\\mathbf{x}) , \\pi_j\">$\\sum_{j=1}^{K} f_j(\\mathbf{x}) , \\pi_j$<\/span> constitue la <strong>densit\u00e9 marginale totale<\/strong> des pr\u00e9dicteurs <span class=\"arabpsy-math-inline\" data-latex=\"P(\\mathbf{X} = \\mathbf{x})\">$P(\\mathbf{X} = \\mathbf{x})$<\/span>, agissant comme une constante de normalisation veillant \u00e0 ce que la somme des probabilit\u00e9s <em>a posteriori<\/em> sur l&rsquo;ensemble des <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span> classes soit rigoureusement \u00e9gale \u00e0 l&rsquo;unit\u00e9 : <span class=\"arabpsy-math-inline\" data-latex=\"\\sum_{k=1}^K P(Y = k | \\mathbf{X} = \\mathbf{x}) = 1\">$\\sum_{k=1}^K P(Y = k | \\mathbf{X} = \\mathbf{x}) = 1$<\/span>.<\/li>\n<\/ul>\n<h3>4.2 Densit\u00e9 conditionnelle gaussienne multivari\u00e9e<\/h3>\n<p>L&rsquo;ADL sp\u00e9cifie la densit\u00e9 de probabilit\u00e9 multivari\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"f_k(\\mathbf{x})\">$f_k(\\mathbf{x})$<\/span> au moyen de la fonction de densit\u00e9 gaussienne \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> dimensions. Pour une classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> caract\u00e9ris\u00e9e par un vecteur d&rsquo;esp\u00e9rances ou centro\u00efde <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_k in \\mathbb{R}^p\">$boldsymbol{\\mu}_k in \\mathbb{R}^p$<\/span> et partageant avec toutes les autres classes une matrice de variance-covariance unique commune <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma} in \\mathbb{R}^{p \\times p}\">$boldsymbol{\\Sigma} in \\mathbb{R}^{p \\times p}$<\/span> sym\u00e9trique d\u00e9finie positive, l&rsquo;expression analytique s&rsquo;\u00e9crit :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"f_k(\\mathbf{x}) = \\frac{1}{(2\\pi)^{p\/2} |boldsymbol{\\Sigma}|^{1\/2}} \\exp\\left( -\\frac{1}{2} (\\mathbf{x} - boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} - boldsymbol{\\mu}_k) \\right)\">$$f_k(\\mathbf{x}) = \\frac{1}{(2\\pi)^{p\/2} |boldsymbol{\\Sigma}|^{1\/2}} \\exp\\left( -\\frac{1}{2} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k) \\right)$$<\/span><\/p>\n<p>Au c\u0153ur de cette expression math\u00e9matique se trouve un op\u00e9rateur m\u00e9trique : la forme quadratique <span class=\"arabpsy-math-inline\" data-latex=\"(\\mathbf{x} - boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} - boldsymbol{\\mu}_k)\">$(\\mathbf{x} &#8211; boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k)$<\/span>. Cette grandeur scalaire correspond au carr\u00e9 de la <strong>distance statistique de Mahalanobis<\/strong>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"D_M^2(\\mathbf{x}, boldsymbol{\\mu}_k)\">$D_M^2(\\mathbf{x}, boldsymbol{\\mu}_k)$<\/span>, s\u00e9parant l&rsquo;observation <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}\">$\\mathbf{x}$<\/span> du centro\u00efde <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_k\">$boldsymbol{\\mu}_k$<\/span> dans la m\u00e9trique induite par la covariance intra-groupe.<\/p>\n<p>Contrairement \u00e0 la distance euclidienne classique \u2014 qui postule implicitement des axes isotropes non corr\u00e9l\u00e9s d&rsquo;\u00e9gale variance \u2014, la distance de Mahalanobis d\u00e9forme l&rsquo;espace g\u00e9om\u00e9trique en pond\u00e9rant chaque variable par l&rsquo;inverse de sa variance r\u00e9siduelle et en orthogonalisant les corr\u00e9lations mutuelles :<\/p>\n<ul>\n<li>Elle standardise les composantes h\u00e9t\u00e9rog\u00e8nes des batteries de tests (annulant les disparit\u00e9s d&rsquo;unit\u00e9s de mesure entre temps de r\u00e9action millisecond\u00e9s et scores d&rsquo;inventaires d&rsquo;humeur born\u00e9s de 0 \u00e0 100).<\/li>\n<li>Elle neutralise la redondance informationnelle : deux variables hautement corr\u00e9l\u00e9es re\u00e7oivent un poids effectif ajust\u00e9 de fa\u00e7on \u00e0 ne pas compter double dans l&rsquo;estimation de l&rsquo;\u00e9cart au centro\u00efde.<\/li>\n<li>Sur le plan g\u00e9om\u00e9trique, les hypersurfaces d&rsquo;isodensit\u00e9 associ\u00e9es \u00e0 chaque classe forment des ellipso\u00efdes de confiance multidimensionnels concentriques, centr\u00e9s en <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_k\">$boldsymbol{\\mu}_k$<\/span>, dont les axes principaux sont orient\u00e9s le long des vecteurs propres de <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span>.<\/li>\n<\/ul>\n<h3>4.3 R\u00e8gle de d\u00e9cision bay\u00e9sienne optimale<\/h3>\n<p>Dans un sch\u00e9ma d\u00e9cisionnel bay\u00e9sien sym\u00e9trique o\u00f9 l&rsquo;on postule une fonction de co\u00fbt de mauvaise classification unitaire dite \u00ab 0-1 \u00bb (qui sanctionne uniform\u00e9ment toute affectation erron\u00e9e sans privil\u00e9gier la pr\u00e9vention d&rsquo;un faux n\u00e9gatif sur un faux positif), l&rsquo;assignation optimale minimise le risque de Bayes en allouant l&rsquo;individu \u00e0 la classe pour laquelle sa probabilit\u00e9 conditionnelle <em>a posteriori<\/em> est maximale :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{Y} = arg\\max_{k in {1, dots, K}} P(Y = k | \\mathbf{X} = \\mathbf{x})\">$$\\hat{Y} = arg\\max_{k in {1, dots, K}} P(Y = k | \\mathbf{X} = \\mathbf{x})$$<\/span><\/p>\n<p>Compte tenu de la stricte monotonie croissante de la fonction logarithme n\u00e9p\u00e9rien, maximiser <span class=\"arabpsy-math-inline\" data-latex=\"P(Y = k | \\mathbf{X} = \\mathbf{x})\">$P(Y = k | \\mathbf{X} = \\mathbf{x})$<\/span> \u00e9quivaut \u00e0 maximiser son log-ratio :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\ln P(Y = k | \\mathbf{X} = \\mathbf{x}) = \\ln f_k(\\mathbf{x}) + \\ln \\pi_k - \\ln P(\\mathbf{X} = \\mathbf{x})\">$$\\ln P(Y = k | \\mathbf{X} = \\mathbf{x}) = \\ln f_k(\\mathbf{x}) + \\ln \\pi_k &#8211; \\ln P(\\mathbf{X} = \\mathbf{x})$$<\/span><\/p>\n<p>Le terme <span class=\"arabpsy-math-inline\" data-latex=\"-\\ln P(\\mathbf{X} = \\mathbf{x})\">$-\\ln P(\\mathbf{X} = \\mathbf{x})$<\/span> \u00e9tant strictement invariant par rapport \u00e0 l&rsquo;indice de classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, il agit comme une constante additive commune et peut \u00eatre ignor\u00e9 dans l&rsquo;algorithme d&rsquo;optimisation. En substituant l&rsquo;expression de la densit\u00e9 normale multivari\u00e9e au sein de cette \u00e9quation, on obtient :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\ln f_k(\\mathbf{x}) + \\ln \\pi_k = -\\frac{p}{2} \\ln(2\\pi) - \\frac{1}{2}ln|boldsymbol{\\Sigma}| - \\frac{1}{2} (\\mathbf{x} - boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} - boldsymbol{\\mu}_k) + \\ln \\pi_k\">$$\\ln f_k(\\mathbf{x}) + \\ln \\pi_k = -\\frac{p}{2} \\ln(2\\pi) &#8211; \\frac{1}{2}ln|boldsymbol{\\Sigma}| &#8211; \\frac{1}{2} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k) + \\ln \\pi_k$$<\/span><\/p>\n<p>Les deux premiers termes, <span class=\"arabpsy-math-inline\" data-latex=\"-\\frac{p}{2} \\ln(2\\pi)\">$-\\frac{p}{2} \\ln(2\\pi)$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"-\\frac{1}{2}ln|boldsymbol{\\Sigma}|\">$-\\frac{1}{2}ln|boldsymbol{\\Sigma}|$<\/span>, sont rigoureusement identiques pour toutes les cat\u00e9gories en vertu de l&rsquo;homosc\u00e9dasticit\u00e9 multivari\u00e9e postulant une matrice <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span> unique et partag\u00e9e. Ils peuvent donc \u00eatre omis. D\u00e9veloppons maintenant la forme quadratique centrale :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"-\\frac{1}{2} (\\mathbf{x} - boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} - boldsymbol{\\mu}_k) = -\\frac{1}{2} \\left[ \\mathbf{x}^T boldsymbol{\\Sigma}^{-1} \\mathbf{x} - 2 boldsymbol{\\mu}_k^T boldsymbol{\\Sigma}^{-1} \\mathbf{x} + boldsymbol{\\mu}_k^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_k \\right]\">$$-\\frac{1}{2} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}^{-1} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k) = -\\frac{1}{2} \\left[ \\mathbf{x}^T boldsymbol{\\Sigma}^{-1} \\mathbf{x} &#8211; 2 boldsymbol{\\mu}_k^T boldsymbol{\\Sigma}^{-1} \\mathbf{x} + boldsymbol{\\mu}_k^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_k \\right]$$<\/span><\/p>\n<p>Le terme quadratique r\u00e9siduel <span class=\"arabpsy-math-inline\" data-latex=\"-\\frac{1}{2} \\mathbf{x}^T boldsymbol{\\Sigma}^{-1} \\mathbf{x}\">$-\\frac{1}{2} \\mathbf{x}^T boldsymbol{\\Sigma}^{-1} \\mathbf{x}$<\/span> ne d\u00e9pend nullement de l&rsquo;indice de classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> : il est rigoureusement invariable quelle que soit la cat\u00e9gorie \u00e9valu\u00e9e pour un profil individuel <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}\">$\\mathbf{x}$<\/span> donn\u00e9. En le retranchant syst\u00e9matiquement, il ne subsiste qu&rsquo;une \u00e9quation strictement affine en <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x}\">$\\mathbf{x}$<\/span>, formellement d\u00e9sign\u00e9e sous le nom de <strong>fonction discriminante lin\u00e9aire de classification<\/strong> pour la classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\delta_k(\\mathbf{x})\">$\\delta_k(\\mathbf{x})$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\delta_k(\\mathbf{x}) = \\mathbf{x}^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_k - \\frac{1}{2} boldsymbol{\\mu}_k^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_k + \\ln \\pi_k\">$$\\delta_k(\\mathbf{x}) = \\mathbf{x}^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_k &#8211; \\frac{1}{2} boldsymbol{\\mu}_k^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_k + \\ln \\pi_k$$<\/span><\/p>\n<p>Cette formulation d\u00e9montre l&rsquo;essence math\u00e9matique de l&rsquo;analyse discriminante lin\u00e9aire : l&rsquo;\u00e9limination alg\u00e9brique exacte des composantes quadratiques induite par l&rsquo;homog\u00e9n\u00e9it\u00e9 de covariance garantit que la r\u00e8gle de classement final est une fonction lin\u00e9aire des pr\u00e9dicteurs. La fronti\u00e8re de d\u00e9cision s\u00e9parant deux classes arbitraires <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"l\">$l$<\/span> correspond \u00e0 l&rsquo;ensemble des points d&rsquo;ind\u00e9cision g\u00e9om\u00e9trique v\u00e9rifiant <span class=\"arabpsy-math-inline\" data-latex=\"\\delta_k(\\mathbf{x}) = \\delta_l(\\mathbf{x})\">$\\delta_k(\\mathbf{x}) = \\delta_l(\\mathbf{x})$<\/span>, ce qui d\u00e9finit de mani\u00e8re rigoureuse un hyperplan affine plat dans l&rsquo;espace multidimensionnel <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^p\">$\\mathbb{R}^p$<\/span>.<\/p>\n<h2>5. D\u00e9rivation analytique des fonctions discriminantes lin\u00e9aires<\/h2>\n<h3>5.1 Maximisation du crit\u00e8re de Fisher<\/h3>\n<p>Dans sa d\u00e9rivation originelle de 1936, Ronald Fisher a con\u00e7u le probl\u00e8me sous un angle purement variationnel, ind\u00e9pendamment de toute formulation bay\u00e9sienne. Consid\u00e9rons le cas fondamental de la discrimination entre deux sous-populations (<span class=\"arabpsy-math-inline\" data-latex=\"K = 2\">$K = 2$<\/span>). Fisher cherchait \u00e0 projeter le vecteur des pr\u00e9dicteurs <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{x} in \\mathbb{R}^p\">$\\mathbf{x} in \\mathbb{R}^p$<\/span> sur un vecteur de projection unidimensionnel <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w} in \\mathbb{R}^p\">$\\mathbf{w} in \\mathbb{R}^p$<\/span> pour g\u00e9n\u00e9rer un score scalaire composite unique <span class=\"arabpsy-math-inline\" data-latex=\"y = \\mathbf{w}^T \\mathbf{x}\">$y = \\mathbf{w}^T \\mathbf{x}$<\/span>, de telle sorte que les scores des deux groupes soient mutuellement les plus s\u00e9par\u00e9s possible.<\/p>\n<p>Pour formaliser cet objectif, d\u00e9finissons au pr\u00e9alable les deux matrices structurales de dispersion :<\/p>\n<ul>\n<li>La <strong>matrice de dispersion intra-classe (<em>Within-class scatter matrix<\/em>)<\/strong>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W\">$\\mathbf{S}_W$<\/span> :\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{S}_W = \\sum_{k=1}^{K} \\sum_{i in C_k} (\\mathbf{x}_i - \\bar{\\mathbf{x}}_k)(\\mathbf{x}_i - \\bar{\\mathbf{x}}_k)^T\">$$\\mathbf{S}_W = \\sum_{k=1}^{K} \\sum_{i in C_k} (\\mathbf{x}_i &#8211; \\bar{\\mathbf{x}}_k)(\\mathbf{x}_i &#8211; \\bar{\\mathbf{x}}_k)^T$$<\/span><\/p>\n<p>    o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{\\mathbf{x}}_k\">$\\bar{\\mathbf{x}}_k$<\/span> repr\u00e9sente le vecteur moyen empirique de la classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>. Cette matrice mesure la variabilit\u00e9 g\u00e9om\u00e9trique r\u00e9siduelle interne des individus autour de leurs centro\u00efdes respectifs. Divis\u00e9e par les degr\u00e9s de libert\u00e9 appropri\u00e9s, elle constitue l&rsquo;estimateur sans biais de <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span>.<\/li>\n<li>La <strong>matrice de dispersion inter-classe (<em>Between-class scatter matrix<\/em>)<\/strong>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_B\">$\\mathbf{S}_B$<\/span> :\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{S}_B = \\sum_{k=1}^{K} N_k (\\bar{\\mathbf{x}}_k - \\bar{\\mathbf{x}})(\\bar{\\mathbf{x}}_k - \\bar{\\mathbf{x}})^T\">$$\\mathbf{S}_B = \\sum_{k=1}^{K} N_k (\\bar{\\mathbf{x}}_k &#8211; \\bar{\\mathbf{x}})(\\bar{\\mathbf{x}}_k &#8211; \\bar{\\mathbf{x}})^T$$<\/span><\/p>\n<p>    o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"N_k\">$N_k$<\/span> est l&rsquo;effectif du groupe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{\\mathbf{x}}\">$\\bar{\\mathbf{x}}$<\/span> le centre de gravit\u00e9 global de l&rsquo;ensemble de l&rsquo;\u00e9chantillon. Pour deux classes, cette matrice se simplifie selon la formule directe :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{S}_B = (\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2)(\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2)^T\">$$\\mathbf{S}_B = (\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2)(\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2)^T$$<\/span><\/p>\n<\/li>\n<\/ul>\n<p>Apr\u00e8s projection vectorielle le long de l&rsquo;axe directionnel <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}\">$\\mathbf{w}$<\/span>, la variance inter-classes projet\u00e9e scalaire s&rsquo;exprime par <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}^T \\mathbf{S}_B \\mathbf{w}\">$\\mathbf{w}^T \\mathbf{S}_B \\mathbf{w}$<\/span>, tandis que la variance intra-classe projet\u00e9e s&rsquo;\u00e9crit <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w}\">$\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w}$<\/span>. Le g\u00e9nie m\u00e9thodologique de Fisher a \u00e9t\u00e9 de formaliser la recherche du vecteur optimal <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}^*\">$\\mathbf{w}^*$<\/span> par l&rsquo;optimisation d&rsquo;un ratio de formes quadratiques, universellement connu en analyse spectrale sous le nom de <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Rayleigh_quotient\">quotient de Rayleigh<\/a> g\u00e9n\u00e9ralis\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"J(\\mathbf{w}) = \\frac{\\mathbf{w}^T \\mathbf{S}_B \\mathbf{w}}{\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w}}\">$$J(\\mathbf{w}) = \\frac{\\mathbf{w}^T \\mathbf{S}_B \\mathbf{w}}{\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w}}$$<\/span><\/p>\n<p>L&rsquo;objectif analytique consiste \u00e0 trouver le vecteur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}\">$\\mathbf{w}$<\/span> qui rend ce quotient maximal, assurant un contraste maximal des moyennes pond\u00e9r\u00e9 par la compacit\u00e9 interne des nuages de points.<\/p>\n<h3>5.2 D\u00e9composition spectrale et calcul des vecteurs propres<\/h3>\n<p>Puisque l&rsquo;amplitude scalaire de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}\">$\\mathbf{w}$<\/span> peut \u00eatre mise \u00e0 l&rsquo;\u00e9chelle arbitrairement sans modifier la valeur du quotient <span class=\"arabpsy-math-inline\" data-latex=\"J(\\mathbf{w})\">$J(\\mathbf{w})$<\/span> (la direction vectorielle spatiale \u00e9tant seule signifiante), le probl\u00e8me d&rsquo;optimisation sans contrainte peut \u00eatre reformul\u00e9 comme un probl\u00e8me d&rsquo;optimisation quadratique sous contrainte d&rsquo;\u00e9galit\u00e9 :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\max_{\\mathbf{w}} ; \\mathbf{w}^T \\mathbf{S}_B \\mathbf{w} \\quad \\text{sujet \u00e0} \\quad \\mathbf{w}^T \\mathbf{S}_W \\mathbf{w} = 1\">$$\\max_{\\mathbf{w}} ; \\mathbf{w}^T \\mathbf{S}_B \\mathbf{w} \\quad \\text{sujet \u00e0} \\quad \\mathbf{w}^T \\mathbf{S}_W \\mathbf{w} = 1$$<\/span><\/p>\n<p>La r\u00e9solution de ce programme d&rsquo;optimisation sous contrainte mobilise canoniquement la m\u00e9thode des multiplicateurs de Lagrange. D\u00e9finissons la fonction lagrangienne <span class=\"arabpsy-math-inline\" data-latex=\"\\mathcal{L}(\\mathbf{w}, \\lambda)\">$\\mathcal{L}(\\mathbf{w}, \\lambda)$<\/span> associ\u00e9e :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathcal{L}(\\mathbf{w}, \\lambda) = \\mathbf{w}^T \\mathbf{S}_B \\mathbf{w} - \\lambda (\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w} - 1)\">$$\\mathcal{L}(\\mathbf{w}, \\lambda) = \\mathbf{w}^T \\mathbf{S}_B \\mathbf{w} &#8211; \\lambda (\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w} &#8211; 1)$$<\/span><\/p>\n<p>En annulant la d\u00e9riv\u00e9e vectorielle partielle du lagrangien par rapport \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}\">$\\mathbf{w}$<\/span>, on obtient la condition de stationnarit\u00e9 du premier ordre :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\frac{\\partial \\mathcal{L}}{\\partial \\mathbf{w}} = 2 \\mathbf{S}_B \\mathbf{w} - 2 \\lambda \\mathbf{S}_W \\mathbf{w} = \\mathbf{0} implies \\mathbf{S}_B \\mathbf{w} = \\lambda \\mathbf{S}_W \\mathbf{w}\">$$\\frac{\\partial \\mathcal{L}}{\\partial \\mathbf{w}} = 2 \\mathbf{S}_B \\mathbf{w} &#8211; 2 \\lambda \\mathbf{S}_W \\mathbf{w} = \\mathbf{0} implies \\mathbf{S}_B \\mathbf{w} = \\lambda \\mathbf{S}_W \\mathbf{w}$$<\/span><\/p>\n<p>Cette formulation constitue l&rsquo;expression math\u00e9matique classique d&rsquo;un <strong>probl\u00e8me de valeurs propres g\u00e9n\u00e9ralis\u00e9<\/strong>. Sous r\u00e9serve que la matrice de dispersion intra-classe <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W\">$\\mathbf{S}_W$<\/span> soit non singuli\u00e8re et donc inversible (ce qui est v\u00e9rifi\u00e9 lorsque la taille de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span> exc\u00e8de substantiellement le nombre de pr\u00e9dicteurs <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> et en l&rsquo;absence de colin\u00e9arit\u00e9 parfaite), nous pouvons pr\u00e9multiplier les deux membres par l&rsquo;inverse <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W^{-1}\">$\\mathbf{S}_W^{-1}$<\/span> :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{S}_W^{-1} \\mathbf{S}_B \\mathbf{w} = \\lambda \\mathbf{w}\">$$\\mathbf{S}_W^{-1} \\mathbf{S}_B \\mathbf{w} = \\lambda \\mathbf{w}$$<\/span><\/p>\n<p>Le probl\u00e8me se r\u00e9duit ainsi \u00e0 la d\u00e9composition spectrale standard de la matrice non sym\u00e9trique <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W^{-1} \\mathbf{S}_B\">$\\mathbf{S}_W^{-1} \\mathbf{S}_B$<\/span>. Dans le cas particulier \u00e0 deux classes (<span class=\"arabpsy-math-inline\" data-latex=\"K = 2\">$K = 2$<\/span>), la matrice de dispersion inter-classes <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_B\">$\\mathbf{S}_B$<\/span> poss\u00e8de un rang matriciel \u00e9gal \u00e0 1, car elle d\u00e9rive d&rsquo;un produit tensoriel de dimension 1 : <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_B \\mathbf{w} = (\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2) \\left[ (\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2)^T \\mathbf{w} \\right]\">$\\mathbf{S}_B \\mathbf{w} = (\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2) \\left[ (\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2)^T \\mathbf{w} \\right]$<\/span>. Le terme entre crochets \u00e9tant un pur scalaire, le vecteur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_B \\mathbf{w}\">$\\mathbf{S}_B \\mathbf{w}$<\/span> est obligatoirement colin\u00e9aire au vecteur de diff\u00e9rence des centro\u00efdes <span class=\"arabpsy-math-inline\" data-latex=\"(\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2)\">$(\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2)$<\/span>. En injectant cette propri\u00e9t\u00e9 dans l&rsquo;\u00e9quation aux valeurs propres, la direction du vecteur discriminant optimal de Fisher se r\u00e9sout explicitement sans d\u00e9composition it\u00e9rative :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{w}^* propto \\mathbf{S}_W^{-1} (\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2)\">$$\\mathbf{w}^* propto \\mathbf{S}_W^{-1} (\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2)$$<\/span><\/p>\n<p>Dans le cas polytomique g\u00e9n\u00e9ral impliquant <span class=\"arabpsy-math-inline\" data-latex=\"K ge 3\">$K ge 3$<\/span> groupes distincts, le rang de la matrice <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_B\">$\\mathbf{S}_B$<\/span> est math\u00e9matiquement born\u00e9 par <span class=\"arabpsy-math-inline\" data-latex=\"min(p, K - 1)\">$min(p, K &#8211; 1)$<\/span>. Il existe donc au maximum <span class=\"arabpsy-math-inline\" data-latex=\"m = min(p, K - 1)\">$m = min(p, K &#8211; 1)$<\/span> valeurs propres r\u00e9elles positives non nulles <span class=\"arabpsy-math-inline\" data-latex=\"\\lambda_1 ge \\lambda_2 ge dots ge \\lambda_m &gt; 0\">$\\lambda_1 ge \\lambda_2 ge dots ge \\lambda_m > 0$<\/span> associ\u00e9es \u00e0 une famille de vecteurs propres lin\u00e9airement ind\u00e9pendants <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}_1, \\mathbf{w}_2, dots, \\mathbf{w}_m\">$\\mathbf{w}_1, \\mathbf{w}_2, dots, \\mathbf{w}_m$<\/span>. Ces axes constituent les <em>fonctions discriminantes canoniques successives<\/em>. Chaque valeur propre <span class=\"arabpsy-math-inline\" data-latex=\"\\lambda_r\">$\\lambda_r$<\/span> traduit la proportion de variance inter-classes isol\u00e9e par l&rsquo;axe factoriel <span class=\"arabpsy-math-inline\" data-latex=\"r\">$r$<\/span> correspondant, quantifiant directement son efficacit\u00e9 diff\u00e9rentielle.<\/p>\n<h3>5.3 \u00c9quations des hyperplans s\u00e9parateurs<\/h3>\n<p>Une fois le vecteur directeur <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}^*\">$\\mathbf{w}^*$<\/span> d\u00e9termin\u00e9 pour s\u00e9parer deux classes, la partition effective de l&rsquo;espace g\u00e9om\u00e9trique requiert le calcul d&rsquo;une constante scalaire de coupure, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"c\">$c$<\/span>, d\u00e9finissant la fronti\u00e8re o\u00f9 l&rsquo;assignation bascule d&rsquo;une cat\u00e9gorie \u00e0 l&rsquo;autre selon que <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}^{*T} \\mathbf{x} - c &gt; 0\">$\\mathbf{w}^{*T} \\mathbf{x} &#8211; c > 0$<\/span> ou $< 0$.<\/p>\n<p>Sous la perspective g\u00e9om\u00e9trique \u00e9quilibr\u00e9e de Fisher, le seuil de d\u00e9cision se positionne \u00e0 mi-distance exacte entre les projections scalaires des centro\u00efdes des deux groupes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"c = \\frac{1}{2} \\mathbf{w}^{*T} (\\bar{\\mathbf{x}}_1 + \\bar{\\mathbf{x}}_2) = \\frac{1}{2} (\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2)^T \\mathbf{S}_W^{-1} (\\bar{\\mathbf{x}}_1 + \\bar{\\mathbf{x}}_2)\">$$c = \\frac{1}{2} \\mathbf{w}^{*T} (\\bar{\\mathbf{x}}_1 + \\bar{\\mathbf{x}}_2) = \\frac{1}{2} (\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2)^T \\mathbf{S}_W^{-1} (\\bar{\\mathbf{x}}_1 + \\bar{\\mathbf{x}}_2)$$<\/span><\/p>\n<p>L&rsquo;\u00e9quation de l&rsquo;hyperplan s\u00e9parateur \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"(p-1)\">$(p-1)$<\/span> dimensions est formellement donn\u00e9e par :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathcal{H} = \\left{ \\mathbf{x} in \\mathbb{R}^p ;middle|; (\\bar{\\mathbf{x}}_1 - \\bar{\\mathbf{x}}_2)^T \\mathbf{S}_W^{-1} \\left( \\mathbf{x} - \\frac{\\bar{\\mathbf{x}}_1 + \\bar{\\mathbf{x}}_2}{2} \\right) = 0 \\right}\">$$\\mathcal{H} = \\left{ \\mathbf{x} in \\mathbb{R}^p ;middle|; (\\bar{\\mathbf{x}}_1 &#8211; \\bar{\\mathbf{x}}_2)^T \\mathbf{S}_W^{-1} \\left( \\mathbf{x} &#8211; \\frac{\\bar{\\mathbf{x}}_1 + \\bar{\\mathbf{x}}_2}{2} \\right) = 0 \\right}$$<\/span><\/p>\n<p>Lorsque le statisticien r\u00e9int\u00e8gre le cadre probabiliste bay\u00e9sien afin de mod\u00e9liser des probabilit\u00e9s <em>a priori<\/em> in\u00e9gales (<span class=\"arabpsy-math-inline\" data-latex=\"\\pi_1 ne \\pi_2\">$\\pi_1 ne \\pi_2$<\/span>), le seuil g\u00e9om\u00e9trique subit une translation parall\u00e8le proportionnelle au logarithme n\u00e9p\u00e9rien du rapport des pr\u00e9valences :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"c_{\\text{Bayes}} = \\frac{1}{2} boldsymbol{\\mu}_1^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_1 - \\frac{1}{2} boldsymbol{\\mu}_2^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_2 - \\ln\\left(\\frac{\\pi_1}{\\pi_2}\\right)\">$$c_{\\text{Bayes}} = \\frac{1}{2} boldsymbol{\\mu}_1^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_1 &#8211; \\frac{1}{2} boldsymbol{\\mu}_2^T boldsymbol{\\Sigma}^{-1} boldsymbol{\\mu}_2 &#8211; \\ln\\left(\\frac{\\pi_1}{\\pi_2}\\right)$$<\/span><\/p>\n<p>L&rsquo;hyperplan s\u00e9parateur affine conserve rigoureusement la m\u00eame orientation angulaire dans l&rsquo;espace multidimensionnel \u2014 son vecteur normal orthogonal demeurant <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}^* = boldsymbol{\\Sigma}^{-1}(boldsymbol{\\mu}_1 - boldsymbol{\\mu}_2)\">$\\mathbf{w}^* = boldsymbol{\\Sigma}^{-1}(boldsymbol{\\mu}_1 &#8211; boldsymbol{\\mu}_2)$<\/span> \u2014, mais sa position spatiale est translat\u00e9e le long de ce vecteur vers le centre de gravit\u00e9 de la classe la moins probable. Cette translation ajuste math\u00e9matiquement la fronti\u00e8re d\u00e9cisionnelle afin de restreindre le volume de l&rsquo;espace affect\u00e9 \u00e0 la classe rare, compensant exactement la dissym\u00e9trie des pr\u00e9valences d\u00e9mographiques.<\/p>\n<h2>6. R\u00e9duction de dimensionnalit\u00e9 et g\u00e9om\u00e9trie de l&rsquo;espace latent<\/h2>\n<h3>6.1 Projection canonique et sous-espaces discriminants<\/h3>\n<p>Au-del\u00e0 de sa fonction de classification supervis\u00e9e, l&rsquo;analyse discriminante lin\u00e9aire s&rsquo;affirme comme une m\u00e9thode de r\u00e9duction de la dimensionnalit\u00e9 guid\u00e9e par les \u00e9tiquettes de classes. Lorsqu&rsquo;un protocole neuropsychologique ou psychom\u00e9trique mobilise plusieurs dizaines de variables (<span class=\"arabpsy-math-inline\" data-latex=\"p gg 1\">$p gg 1$<\/span>) pour diff\u00e9rencier un nombre restreint de groupes <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span>, la visualisation directe et la mod\u00e9lisation parcimonieuse deviennent intenables au sein de l&rsquo;espace g\u00e9om\u00e9trique originel.<\/p>\n<p>L&rsquo;ADL r\u00e9sout ce probl\u00e8me en projetant les donn\u00e9es depuis l&rsquo;espace de d\u00e9part <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^p\">$\\mathbb{R}^p$<\/span> vers un sous-espace latent affine <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^m\">$\\mathbb{R}^m$<\/span>, o\u00f9 la dimension maximale <span class=\"arabpsy-math-inline\" data-latex=\"m\">$m$<\/span> est strictement contrainte par la g\u00e9om\u00e9trie inter-groupes :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"m = \\min(p, K - 1)\">$$m = \\min(p, K &#8211; 1)$$<\/span><\/p>\n<p>Cette propri\u00e9t\u00e9 fondamentale d\u00e9coule du fait que les centro\u00efdes des <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span> populations <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_1, dots, boldsymbol{\\mu}_K\">$boldsymbol{\\mu}_1, dots, boldsymbol{\\mu}_K$<\/span> sous-tendent au maximum un sous-espace vectoriel de dimension <span class=\"arabpsy-math-inline\" data-latex=\"K - 1\">$K &#8211; 1$<\/span>. Par exemple, pour distinguer trois entit\u00e9s nosologiques (<span class=\"arabpsy-math-inline\" data-latex=\"K = 3\">$K = 3$<\/span>) au moyen d&rsquo;une batterie de vingt tests cognitifs (<span class=\"arabpsy-math-inline\" data-latex=\"p = 20\">$p = 20$<\/span>), la totalit\u00e9 de l&rsquo;information diff\u00e9rentielle inter-classes est int\u00e9gralement encapsul\u00e9e dans un plan factoriel bidimensionnel (<span class=\"arabpsy-math-inline\" data-latex=\"m = min(20, 2) = 2\">$m = min(20, 2) = 2$<\/span>).<\/p>\n<p>La transformation s&rsquo;op\u00e8re par le produit matriciel :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{z}_i = \\mathbf{W}^T \\mathbf{x}_i\">$$\\mathbf{z}_i = \\mathbf{W}^T \\mathbf{x}_i$$<\/span><\/p>\n<p>o\u00f9 la matrice de projection <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{W} in \\mathbb{R}^{p \\times m}\">$\\mathbf{W} in \\mathbb{R}^{p \\times m}$<\/span> r\u00e9unit en colonnes les <span class=\"arabpsy-math-inline\" data-latex=\"m\">$m$<\/span> vecteurs propres canoniques de la matrice <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W^{-1} \\mathbf{S}_B\">$\\mathbf{S}_W^{-1} \\mathbf{S}_B$<\/span>. Les coordonn\u00e9es <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{z}_i\">$\\mathbf{z}_i$<\/span> obtenues sont appel\u00e9es <em>variables canoniques discriminantes<\/em>. Ces nouveaux axes synth\u00e9tiques poss\u00e8dent la propri\u00e9t\u00e9 d&rsquo;\u00eatre mutuellement non corr\u00e9l\u00e9s au sens de la m\u00e9trique intra-classe :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{w}_i^T \\mathbf{S}_W \\mathbf{w}_j = \\begin{\\cases} 1 &amp; \\text{si } i = j \\ 0 &amp; \\text{si } i ne j \\end{\\cases}\">$$\\mathbf{w}_i^T \\mathbf{S}_W \\mathbf{w}_j = \\begin{\\cases} 1 &#038; \\text{si } i = j \\ 0 &#038; \\text{si } i ne j \\end{\\cases}$$<\/span><\/p>\n<p>L&rsquo;espace g\u00e9om\u00e9trique r\u00e9sultant offre une compression optimale : il \u00e9limine toutes les dimensions de variations psychom\u00e9triques au sein desquelles les groupes ne manifestent aucun \u00e9cart de moyenne, concentrant la variance diff\u00e9rentielle sur une poign\u00e9e d&rsquo;axes interpr\u00e9tables.<\/p>\n<h3>6.2 Comparaison syst\u00e9matique entre ACP et ADL<\/h3>\n<p>La distinction conceptuelle entre l&rsquo;Analyse en Composantes Principales (m\u00e9thode de projection orthogonale non supervis\u00e9e) et l&rsquo;Analyse Discriminante Lin\u00e9aire (projection oblique canonique supervis\u00e9e) structure l&rsquo;\u00e9pist\u00e9mologie de l&rsquo;analyse multivari\u00e9e. Ces deux m\u00e9thodes poursuivent des optimisations variationnelles radicalement distinctes :<\/p>\n<table border=\"1\" cellspacing=\"0\" cellpadding=\"8\" style=\"width:100%; border-collapse:collapse; text-align:left; margin: 20px 0;\">\n<thead>\n<tr style=\"background-color: #f2f2f2;\">\n<th>Crit\u00e8re m\u00e9thodologique<\/th>\n<th>Analyse en composantes principales (ACP)<\/th>\n<th>Analyse discriminante lin\u00e9aire (ADL)<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td><strong>Paradigme d&rsquo;apprentissage<\/strong><\/td>\n<td>Non supervis\u00e9 (aucune \u00e9tiquette de classe exploit\u00e9e)<\/td>\n<td>Supervis\u00e9 (guid\u00e9 par les labels cat\u00e9goriels de groupes)<\/td>\n<\/tr>\n<tr>\n<td><strong>Objectif variationnel fondamental<\/strong><\/td>\n<td>Maximisation de la variance totale globale projet\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}^T \\mathbf{S}_T \\mathbf{w}\">$\\mathbf{w}^T \\mathbf{S}_T \\mathbf{w}$<\/span>)<\/td>\n<td>Maximisation du rapport de Rayleigh (<span class=\"arabpsy-math-inline\" data-latex=\"\\frac{\\mathbf{w}^T \\mathbf{S}_B \\mathbf{w}}{\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w}}\">$\\frac{\\mathbf{w}^T \\mathbf{S}_B \\mathbf{w}}{\\mathbf{w}^T \\mathbf{S}_W \\mathbf{w}}$<\/span>)<\/td>\n<\/tr>\n<tr>\n<td><strong>Nombre maximal d&rsquo;axes factoriels<\/strong><\/td>\n<td><span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> (nombre total de pr\u00e9dicteurs initiaux)<\/td>\n<td><span class=\"arabpsy-math-inline\" data-latex=\"min(p, K - 1)\">$min(p, K &#8211; 1)$<\/span> (contraint par le nombre de cat\u00e9gories)<\/td>\n<\/tr>\n<tr>\n<td><strong>Sensibilit\u00e9 aux gradients intra-groupes<\/strong><\/td>\n<td>Domin\u00e9e par la variabilit\u00e9 globale, m\u00eame si non s\u00e9paratrice<\/td>\n<td>Neutralise la variabilit\u00e9 intra-groupe pour isoler l&rsquo;\u00e9cartement<\/td>\n<\/tr>\n<tr>\n<td><strong>Sensibilit\u00e9 d&rsquo;\u00e9chelle<\/strong><\/td>\n<td>\u00c9lev\u00e9e (d\u00e9pend de la matrice de covariance vs corr\u00e9lation)<\/td>\n<td>Invariante aux transformations affines inversibles des pr\u00e9dicteurs<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Il est courant qu&rsquo;une ACP \u00e9choue totalement \u00e0 s\u00e9parer des sous-populations distinctes lorsque la direction de discrimination inter-groupes ne co\u00efncide pas avec l&rsquo;axe d&rsquo;\u00e9talement principal de l&rsquo;ensemble des donn\u00e9es. Si, par exemple, la variabilit\u00e9 physiologique normale entre individus est nettement plus vaste que la variation pathologique subtile induite par un stade pr\u00e9coce de maladie neurod\u00e9g\u00e9n\u00e9rative, les composantes principales 1 et 2 de l&rsquo;ACP captureront exclusivement la diversit\u00e9 morphologique saine, superposant totalement les sujets sains et les malades. En revanche, l&rsquo;ADL comprime cette variabilit\u00e9 normale partag\u00e9e via l&rsquo;inversion de la matrice intra-groupe <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W^{-1}\">$\\mathbf{S}_W^{-1}$<\/span>, orientant le regard analytique sur l&rsquo;axe subtil o\u00f9 s&rsquo;effectue la d\u00e9marcation clinique.<\/p>\n<p>Ces deux m\u00e9thodes s&rsquo;av\u00e8rent hautement compl\u00e9mentaires au sein des protocoles m\u00e9thodologiques exploratoires pr\u00e9alables. Dans les situations dites de grande dimensionnalit\u00e9 o\u00f9 le nombre de mesures psychom\u00e9triques exc\u00e8de la taille de la cohorte (<span class=\"arabpsy-math-inline\" data-latex=\"p &gt; N\">$p > N$<\/span>), la matrice <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W\">$\\mathbf{S}_W$<\/span> devient singuli\u00e8re, rendant l&rsquo;application directe de l&rsquo;ADL impossible. Une strat\u00e9gie consiste \u00e0 impl\u00e9menter un pipeline s\u00e9quentiel d&rsquo;ACP-ADL : l&rsquo;ACP est initialement exploit\u00e9e pour d\u00e9bruiter le signal et r\u00e9duire la dimensionalit\u00e9 de <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> \u00e0 $k&rsquo; < N$ composantes principales, sur lesquelles une ADL supervis\u00e9e est ensuite ajust\u00e9e sans instabilit\u00e9 d'inversion matricielle.<\/p>\n<h3>6.3 Visualisation cartographique des profils psychom\u00e9triques<\/h3>\n<p>La projection des participants au sein de l&rsquo;espace factoriel issu de l&rsquo;ADL fournit un outil de cartographie clinique puissant. En positionnant chaque patient sur un rep\u00e8re orthogonal bivari\u00e9 structur\u00e9 par la premi\u00e8re variable canonique (Axe 1 en abscisse) et la seconde variable canonique (Axe 2 en ordonn\u00e9e), l&rsquo;analyste g\u00e9n\u00e8re un portrait spatial intelligible des proximit\u00e9s ph\u00e9notypiques.<\/p>\n<p>Pour conf\u00e9rer une rigueur probabiliste \u00e0 cette cartographie, la repr\u00e9sentation spatiale s&rsquo;enrichit de plusieurs couches d&rsquo;information vectorielle :<\/p>\n<ul>\n<li><strong>Les centro\u00efdes canoniques des groupes :<\/strong> Repr\u00e9sent\u00e9s par des barycentres d&rsquo;attraction g\u00e9om\u00e9trique indiquant le comportement typique moyen de chaque classe dans le sous-espace discriminant.<\/li>\n<li><strong>Les ellipses de concentration \u00e9quiprobables :<\/strong> D\u00e9finies pour des seuils statistiques formalis\u00e9s (par exemple \u00e0 68 % ou 95 % de confiance). Ces contours ellipso\u00efdaux, mod\u00e9lis\u00e9s \u00e0 partir de la distribution bivari\u00e9e du Chi-deux, objectivent visuellement la dispersion r\u00e9siduelle de chaque pathologie. Le degr\u00e9 d&rsquo;interp\u00e9n\u00e9tration ou de disjonction des ellipses mat\u00e9rialise imm\u00e9diatement l&rsquo;ambigu\u00eft\u00e9 diagnostique ou, \u00e0 l&rsquo;inverse, l&rsquo;individualisation nette des entit\u00e9s cliniques.<\/li>\n<li><strong>Les graphiques de biplot factoriel (fl\u00e8ches de projection des pr\u00e9dicteurs) :<\/strong> La superposition des vecteurs structurels repr\u00e9sentant les variables cognitives de d\u00e9part \u00e9claire la s\u00e9mantique clinique des axes factoriels. La longueur d&rsquo;une fl\u00e8che refl\u00e8te son pouvoir d&rsquo;entra\u00eenement discriminant, tandis que son angle renseigne sur ses corr\u00e9lations directes avec les fonctions canoniques latentes.<\/li>\n<\/ul>\n<p>Cette approche permet d&rsquo;identifier visuellement les profils cognitifs atypiques, d&rsquo;explorer les chevauchements nosologiques aux zones de fronti\u00e8res inter-classes, et de rep\u00e9rer les \u00e9ventuels cas d&rsquo;observations aberrantes multidimensionnelles (<em>outliers<\/em>) dont l&rsquo;excentricit\u00e9 par rapport \u00e0 l&rsquo;ellipso\u00efde de leur propre groupe pourrait fausser les inf\u00e9rences statistiques de l&rsquo;\u00e9tude.<\/p>\n<h2>7. Estimation empirique et calibrage des probabilit\u00e9s a priori<\/h2>\n<h3>7.1 Param\u00e9trisation par maximum de vraisemblance<\/h3>\n<p>En pratique empirique, les param\u00e8tres distributionnels th\u00e9oriques r\u00e9gissant la population globale \u2014 \u00e0 savoir les centro\u00efdes de classes <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_k\">$boldsymbol{\\mu}_k$<\/span> et la matrice de variance-covariance intra-groupe commune <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span> \u2014 sont rigoureusement inconnus du chercheur. Ils doivent \u00eatre estim\u00e9s \u00e0 partir des observations r\u00e9elles d&rsquo;un \u00e9chantillon d&rsquo;apprentissage constitu\u00e9 de <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span> individus r\u00e9partis au sein de <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span> groupes de tailles <span class=\"arabpsy-math-inline\" data-latex=\"N_1, N_2, dots, N_K\">$N_1, N_2, dots, N_K$<\/span> telles que <span class=\"arabpsy-math-inline\" data-latex=\"\\sum_{k=1}^K N_k = N\">$\\sum_{k=1}^K N_k = N$<\/span>.<\/p>\n<p>L&rsquo;estimation conjointe sous l&rsquo;hypoth\u00e8se de normalit\u00e9 multivari\u00e9e mobilise les estimateurs du <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Maximum_de_vraisemblance\">maximum de vraisemblance<\/a>, calibr\u00e9s par les corrections usuelles de degr\u00e9s de libert\u00e9 afin de garantir des propri\u00e9t\u00e9s asymptotiques sans biais :<\/p>\n<ul>\n<li>Le vecteur des moyennes empiriques d&rsquo;\u00e9chantillon pour la classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{\\mathbf{x}}_k\">$\\bar{\\mathbf{x}}_k$<\/span>, constitue l&rsquo;estimateur sans biais du centro\u00efde th\u00e9orique <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\mu}_k\">$boldsymbol{\\mu}_k$<\/span> :\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{boldsymbol{\\mu}}_k = \\bar{\\mathbf{x}}_k = \\frac{1}{N_k} \\sum_{i in C_k} \\mathbf{x}_i\">$$\\hat{boldsymbol{\\mu}}_k = \\bar{\\mathbf{x}}_k = \\frac{1}{N_k} \\sum_{i in C_k} \\mathbf{x}_i$$<\/span><\/p>\n<\/li>\n<li>La matrice de covariance empirique individuelle du groupe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_k\">$\\mathbf{S}_k$<\/span>, s&rsquo;\u00e9crit :\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\mathbf{S}_k = \\frac{1}{N_k - 1} \\sum_{i in C_k} (\\mathbf{x}_i - \\bar{\\mathbf{x}}_k)(\\mathbf{x}_i - \\bar{\\mathbf{x}}_k)^T\">$$\\mathbf{S}_k = \\frac{1}{N_k &#8211; 1} \\sum_{i in C_k} (\\mathbf{x}_i &#8211; \\bar{\\mathbf{x}}_k)(\\mathbf{x}_i &#8211; \\bar{\\mathbf{x}}_k)^T$$<\/span><\/p>\n<\/li>\n<li>L&rsquo;estimateur sans biais de la matrice de variance-covariance combin\u00e9e intra-groupe <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}\">$boldsymbol{\\Sigma}$<\/span>, not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}\">$\\mathbf{S}_{\\text{intra}}$<\/span>, s&rsquo;obtient par la moyenne pond\u00e9r\u00e9e des matrices de covariance individuelles, ajust\u00e9e sur les degr\u00e9s de libert\u00e9 r\u00e9siduels totaux :\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{boldsymbol{\\Sigma}} = \\mathbf{S}_{\\text{intra}} = \\frac{1}{N - K} \\sum_{k=1}^{K} (N_k - 1) \\mathbf{S}_k = \\frac{1}{N - K} \\mathbf{S}_W\">$$\\hat{boldsymbol{\\Sigma}} = \\mathbf{S}_{\\text{intra}} = \\frac{1}{N &#8211; K} \\sum_{k=1}^{K} (N_k &#8211; 1) \\mathbf{S}_k = \\frac{1}{N &#8211; K} \\mathbf{S}_W$$<\/span><\/p>\n<\/li>\n<\/ul>\n<p>Ces estimateurs empiriques sont asymptotiquement convergents et d&rsquo;une variance minimale sous r\u00e9serve que la normalit\u00e9 soit respect\u00e9e. Cependant, lorsque la taille de l&rsquo;\u00e9chantillon est restreinte ou que le nombre de pr\u00e9dicteurs <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> augmente, l&rsquo;incertitude d&rsquo;\u00e9chantillonnage grevant chaque terme de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}\">$\\mathbf{S}_{\\text{intra}}$<\/span> se propage dans l&rsquo;inversion matricielle <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}^{-1}\">$\\mathbf{S}_{\\text{intra}}^{-1}$<\/span>, introduisant un biais optimiste syst\u00e9matique dans la performance de s\u00e9paration des fonctions calcul\u00e9es.<\/p>\n<h3>7.2 Int\u00e9gration des pr\u00e9valences r\u00e9elles en sant\u00e9 mentale<\/h3>\n<p>L&rsquo;un des \u00e9cueils m\u00e9thodologiques les plus fr\u00e9quents en psychologie quantitative r\u00e9side dans la confusion entre les fr\u00e9quences marginales observ\u00e9es au sein du jeu d&rsquo;\u00e9chantillonnage de recherche et les probabilit\u00e9s <em>a priori<\/em> <span class=\"arabpsy-math-inline\" data-latex=\"\\pi_k\">$\\pi_k$<\/span> r\u00e9gissant la survenue des affections en population g\u00e9n\u00e9rale. Dans les protocoles de recherche cas-t\u00e9moins, les cohortes sont typiquement compos\u00e9es de fa\u00e7on paritaire (par exemple, 50 patients pr\u00e9sentant un premier \u00e9pisode psychotique appari\u00e9s \u00e0 50 t\u00e9moins sains, fixant arbitrairement une fr\u00e9quence relative d&rsquo;\u00e9chantillon de 50 % \/ 50 %).<\/p>\n<p>Si l&rsquo;analyste laisse le mod\u00e8le discriminant d\u00e9river ses probabilit\u00e9s <em>a priori<\/em> directement \u00e0 partir des effectifs de l&rsquo;\u00e9chantillon (<span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\pi}_k = N_k \/ N\">$\\hat{\\pi}_k = N_k \/ N$<\/span>), le classificateur appliquera une r\u00e8gle bay\u00e9sienne postulant que l&rsquo;apparition d&rsquo;un \u00e9pisode psychotique en vie r\u00e9elle poss\u00e8de une probabilit\u00e9 \u00e9quivalente \u00e0 une chance sur deux. Or, dans la population g\u00e9n\u00e9rale, la pr\u00e9valence \u00e9pid\u00e9miologique r\u00e9elle de cette pathologie est de l&rsquo;ordre de 1 % \u00e0 3 %. Appliquer un tel mod\u00e8le sans correction pour des fins de d\u00e9pistage en population g\u00e9n\u00e9rale engendrera une cascade d\u00e9sastreuse de faux positifs cliniques :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\ln\\left(\\frac{\\hat{\\pi}_{\\text{cas}}}{\\hat{\\pi}_{\\text{t\u00e9moin}}}\\right)_{\\text{\u00e9chantillon}} = \\ln\\left(\\frac{0{,}50}{0{,}50}\\right) = 0 \\quad \\text{alors que} \\quad \\ln\\left(\\frac{\\pi_{\\text{r\u00e9elle}}}{\\pi_{\\text{t\u00e9moin}}}\\right)_{\\text{\u00e9pid\u00e9miologique}} = \\ln\\left(\\frac{0{,}01}{0{,}99}\\right) \\approx -4{,}595\">$$\\ln\\left(\\frac{\\hat{\\pi}_{\\text{cas}}}{\\hat{\\pi}_{\\text{t\u00e9moin}}}\\right)_{\\text{\u00e9chantillon}} = \\ln\\left(\\frac{0{,}50}{0{,}50}\\right) = 0 \\quad \\text{alors que} \\quad \\ln\\left(\\frac{\\pi_{\\text{r\u00e9elle}}}{\\pi_{\\text{t\u00e9moin}}}\\right)_{\\text{\u00e9pid\u00e9miologique}} = \\ln\\left(\\frac{0{,}01}{0{,}99}\\right) \\approx -4{,}595$$<\/span><\/p>\n<p>L&rsquo;omission de cette correction d\u00e9place massivement l&rsquo;hyperplan s\u00e9parateur au d\u00e9triment de la classe saine majoritaire. L&rsquo;analyse discriminante lin\u00e9aire permet d&rsquo;ajuster ce seuil d\u00e9cisionnel en injectant explicitement les pr\u00e9valences \u00e9pid\u00e9miologiques objectives document\u00e9es dans les registres de sant\u00e9 publique au titre des param\u00e8tres <span class=\"arabpsy-math-inline\" data-latex=\"\\pi_k\">$\\pi_k$<\/span>. Cette modification n&rsquo;alt\u00e8re en rien l&rsquo;orientation des fonctions discriminantes (d\u00e9finies par les barycentres et la matrice intra-groupe), mais r\u00e9\u00e9talonne rigoureusement les termes constants pour une pr\u00e9diction adapt\u00e9e au milieu \u00e9cologique de d\u00e9ploiement clinique.<\/p>\n<h3>7.3 Gestion des d\u00e9s\u00e9quilibres structurels de donn\u00e9es<\/h3>\n<p>\u00c0 l&rsquo;inverse des designs exp\u00e9rimentaux artificiellement paritaires, les \u00e9tudes \u00e9pid\u00e9miologiques observationnelles sont fr\u00e9quemment confront\u00e9es \u00e0 des d\u00e9s\u00e9quilibres de classes massifs (par exemple, 95 % d&rsquo;individus asymptomatiques contre 5 % de patients suicidaires). Dans cette configuration, l&rsquo;application pure et simple du crit\u00e8re bay\u00e9sien standard produit un ph\u00e9nom\u00e8ne d&rsquo;\u00e9crasement statistique : la cat\u00e9gorie minoritaire est \u00e9touff\u00e9e par le poids a priori disproportionn\u00e9 de la classe majoritaire, conduisant le mod\u00e8le \u00e0 maximiser son exactitude brute en affectant la quasi-totalit\u00e9 des individus \u00e0 la modalit\u00e9 dominante.<\/p>\n<p>Pour rem\u00e9dier \u00e0 cette asym\u00e9trie pr\u00e9dictive tout en pr\u00e9servant le cadre analytique de l&rsquo;ADL, plusieurs strat\u00e9gies m\u00e9thodologiques doivent \u00eatre envisag\u00e9es :<\/p>\n<ul>\n<li><strong>L&rsquo;ajustement asym\u00e9trique des co\u00fbts d&rsquo;erreur :<\/strong> La matrice de gain\/perte bay\u00e9sienne est modifi\u00e9e pour attribuer un co\u00fbt p\u00e9nalisant sup\u00e9rieur au fait de classer \u00e0 tort un individu vuln\u00e9rable dans la cat\u00e9gorie saine (faux n\u00e9gatif), par rapport au co\u00fbt induit par une fausse alerte (faux positif). Cette asym\u00e9trie compense math\u00e9matiquement la faiblesse du param\u00e8tre <span class=\"arabpsy-math-inline\" data-latex=\"\\pi_k\">$\\pi_k$<\/span>.<\/li>\n<li><strong>Le r\u00e9\u00e9chantillonnage \u00e9quilibr\u00e9 ou synth\u00e9tique en amont :<\/strong> L&rsquo;application d&rsquo;algorithmes de sur\u00e9chantillonnage de la classe minoritaire au sein du jeu d&rsquo;apprentissage \u2014 tels que <a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Oversampling_and_undersampling_in_data_analysis\">SMOTE<\/a> (<em>Synthetic Minority Over-sampling Technique<\/em>) adapt\u00e9 aux donn\u00e9es multivari\u00e9es continues \u2014 permet de r\u00e9\u00e9quilibrer la structure des donn\u00e9es avant l&rsquo;estimation de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}\">$\\mathbf{S}_{\\text{intra}}$<\/span>.<\/li>\n<li><strong>Le recalibrage empirique du seuil de coupure :<\/strong> L&rsquo;abandon de l&rsquo;affectation automatique au mode de la probabilit\u00e9 conditionnelle (<span class=\"arabpsy-math-inline\" data-latex=\"P &gt; 0{,}5\">$P > 0{,}5$<\/span>) au profit d&rsquo;une optimisation empirique du seuil de d\u00e9cision le long de la fonction discriminante, pilot\u00e9e par l&rsquo;indice de Youden sur une courbe ROC d\u00e9di\u00e9e.<\/li>\n<\/ul>\n<h2>8. \u00c9valuation de la performance pr\u00e9dictive et validation du mod\u00e8le<\/h2>\n<h3>8.1 Indicateurs m\u00e9triques et matrices de confusion<\/h3>\n<p>L&rsquo;\u00e9valuation rigoureuse de la qualit\u00e9 discriminante d&rsquo;un mod\u00e8le ajust\u00e9 ne saurait s&rsquo;appuyer sur la seule m\u00e9trique de l&rsquo;exactitude globale (<em>accuracy<\/em>), hautement trompeuse en pr\u00e9sence de groupes d&rsquo;effectifs in\u00e9gaux. L&rsquo;investigation empirique requiert la construction syst\u00e9matique d&rsquo;une <strong>matrice de confusion crois\u00e9e<\/strong>, confrontant les appartenances r\u00e9elles observ\u00e9es aux assignations cat\u00e9gorielles pr\u00e9dites par l&rsquo;algorithme sur un jeu de donn\u00e9es.<\/p>\n<p>\u00c0 partir de cette matrice, l&rsquo;analyste d\u00e9cline un ensemble articul\u00e9 de m\u00e9triques statistiques adapt\u00e9es \u00e0 la d\u00e9cision diagnostique :<\/p>\n<ul>\n<li><strong>La sensibilit\u00e9 clinique (taux de vrais positifs) :<\/strong> Capacit\u00e9 du mod\u00e8le \u00e0 identifier exhaustivement l&rsquo;ensemble des individus porteurs d&rsquo;un \u00e9tat pathologique donn\u00e9.<\/li>\n<li><strong>La sp\u00e9cificit\u00e9 diagnostique (taux de vrais n\u00e9gatifs) :<\/strong> Capacit\u00e9 du mod\u00e8le \u00e0 \u00e9carter scrupuleusement les sujets exempts du diagnostic cibl\u00e9, pr\u00e9venant la stigmatisation li\u00e9e au faux diagnostic.<\/li>\n<li><strong>La valeur pr\u00e9dictive positive (VPP) et la valeur pr\u00e9dictive n\u00e9gative (VPN) :<\/strong> Indicateurs critiques pour le clinicien, mesurant la probabilit\u00e9 conditionnelle qu&rsquo;un diagnostic pr\u00e9dit par l&rsquo;ADL soit av\u00e9r\u00e9 en pratique clinique, compte tenu des pr\u00e9valences appliqu\u00e9es.<\/li>\n<li><strong>Le coefficient Kappa de Cohen (<span class=\"arabpsy-math-inline\" data-latex=\"kappa\">$kappa$<\/span>) :<\/strong> M\u00e9trique inf\u00e9rentielle ajustant le taux de concordance observ\u00e9 entre le mod\u00e8le et la r\u00e9alit\u00e9 empirique en soustrayant le niveau d&rsquo;accord purement fortuit li\u00e9 aux marges de la table. Un <span class=\"arabpsy-math-inline\" data-latex=\"\\kappa &gt; 0{,}75\">$\\kappa > 0{,}75$<\/span> atteste d&rsquo;une concordance substantielle \u00e0 excellente.<\/li>\n<\/ul>\n<p>Dans les contextes multiclasses (<span class=\"arabpsy-math-inline\" data-latex=\"K ge 3\">$K ge 3$<\/span>), ces indicateurs doivent \u00eatre calcul\u00e9s selon des sch\u00e9mas d&rsquo;agr\u00e9gation macro-analytique (moyenne non pond\u00e9r\u00e9e des performances calcul\u00e9es au sein de chaque groupe) ou micro-analytique (pond\u00e9r\u00e9e par les effectifs). De surcro\u00eet, le recours aux <em>courbes ROC<\/em> (Receiver Operating Characteristic) \u00e9tendues aux situations multiclasses, accompagn\u00e9es du calcul de l&rsquo;aire sous la courbe (AUC de Hand et Till), permet d&rsquo;\u00e9valuer la discriminabilit\u00e9 intrins\u00e8que du mod\u00e8le ind\u00e9pendamment de tout seuil d\u00e9cisionnel binaire arbitraire.<\/p>\n<h3>8.2 Protocoles de validation crois\u00e9e rigoureux<\/h3>\n<p>L&rsquo;\u00e9valuation de la performance pr\u00e9dictive sur les observations m\u00eames qui ont servi \u00e0 estimer les centro\u00efdes et la matrice de variance-covariance engendre un <em>biais d&rsquo;optimisme<\/em> syst\u00e9matique, \u00e9galement d\u00e9sign\u00e9 sous le terme de surapprentissage de r\u00e9\u00e9valuation (<em>resubstitution bias<\/em>). Ce ph\u00e9nom\u00e8ne conduit \u00e0 une surestimation dramatique des capacit\u00e9s de g\u00e9n\u00e9ralisation r\u00e9elles de la batterie de tests psychom\u00e9triques.<\/p>\n<p>Pour neutraliser ce biais m\u00e9thodologique, le d\u00e9ploiement de protocoles de validation crois\u00e9e \u00e9tanches constitue un imp\u00e9ratif m\u00e9thodologique incontournable :<\/p>\n<ul>\n<li><strong>La validation crois\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>-fold stratifi\u00e9e :<\/strong> L&rsquo;\u00e9chantillon total est scind\u00e9 en <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> sous-ensembles (typiquement 5 ou 10) mutuellement exclusifs, en veillant scrupuleusement \u00e0 ce que chaque fold reproduise la distribution proportionnelle exacte des classes initiales. Le mod\u00e8le d&rsquo;ADL est entra\u00een\u00e9 sur <span class=\"arabpsy-math-inline\" data-latex=\"k - 1\">$k &#8211; 1$<\/span> sous-ensembles, puis sa performance est \u00e9valu\u00e9e sur le pli r\u00e9siduel gard\u00e9 totalement aveugle. Cette proc\u00e9dure est r\u00e9it\u00e9r\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> fois, puis les m\u00e9triques sont moyenn\u00e9es.<\/li>\n<li><strong>La validation par exclusion unitaire (<em>Leave-One-Out Cross-Validation &#8211; LOOCV<\/em>) :<\/strong> Particuli\u00e8rement indiqu\u00e9e en neuropsychologie face \u00e0 des cohortes tr\u00e8s restreintes ($N < 50<span class=\"arabpsy-math-inline\" data-latex=\"), cette m\u00e9thode consiste \u00e0 masquer un unique individu, \u00e0 estimer l&#039;ADL sur les\">$), cette m\u00e9thode consiste \u00e0 masquer un unique individu, \u00e0 estimer l&rsquo;ADL sur les$<\/span>N &#8211; 1$ patients restants, puis \u00e0 tester l&rsquo;assignation de l&rsquo;individu exclu. R\u00e9p\u00e9t\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span> fois, elle fournit un estimateur quasi sans biais de l&rsquo;erreur r\u00e9elle de g\u00e9n\u00e9ralisation.<\/li>\n<li><strong>La validation externe ind\u00e9pendante :<\/strong> L&rsquo;\u00e9talon-or ultime demeure l&rsquo;application du mod\u00e8le fig\u00e9 \u00e0 une cohorte de r\u00e9plication prospective totalement ind\u00e9pendante, recrut\u00e9e au sein d&rsquo;autres centres hospitaliers ou universitaires, afin d&rsquo;attester de la transf\u00e9rabilit\u00e9 inter-sites et de la robustesse psychom\u00e9trique du mod\u00e8le face aux variations de passation des \u00e9preuves.<\/li>\n<\/ul>\n<h3>8.3 Statistiques inf\u00e9rentielles de s\u00e9paration globale<\/h3>\n<p>L&rsquo;analyse discriminante lin\u00e9aire int\u00e8gre un cadre inf\u00e9rentiel multivari\u00e9 rigoureux permettant de tester formellement l&rsquo;hypoth\u00e8se nulle d&rsquo;\u00e9galit\u00e9 absolue des profils psychom\u00e9triques moyens \u00e0 travers les <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span> groupes exp\u00e9rimentaux :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"H_0 : boldsymbol{\\mu}_1 = boldsymbol{\\mu}_2 = dots = boldsymbol{\\mu}_K\">$$H_0 : boldsymbol{\\mu}_1 = boldsymbol{\\mu}_2 = dots = boldsymbol{\\mu}_K$$<\/span><\/p>\n<p>La statistique multivari\u00e9e la plus universellement rapport\u00e9e est le <strong>Lambda de Wilks (<span class=\"arabpsy-math-inline\" data-latex=\"Lambda\">$Lambda$<\/span>)<\/strong>. Sur le plan alg\u00e9brique, elle repr\u00e9sente le rapport du d\u00e9terminant de la matrice de dispersion intra-classe au d\u00e9terminant de la matrice de dispersion totale :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\Lambda = \\frac{|\\mathbf{S}_W|}{|\\mathbf{S}_T|} = \\frac{|\\mathbf{S}_W|}{|\\mathbf{S}_W + \\mathbf{S}_B|} = \\prod_{r=1}^{m} \\frac{1}{1 + \\lambda_r}\">$$\\Lambda = \\frac{|\\mathbf{S}_W|}{|\\mathbf{S}_T|} = \\frac{|\\mathbf{S}_W|}{|\\mathbf{S}_W + \\mathbf{S}_B|} = \\prod_{r=1}^{m} \\frac{1}{1 + \\lambda_r}$$<\/span><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"\\lambda_r\">$\\lambda_r$<\/span> correspond aux valeurs propres issues de la d\u00e9composition de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W^{-1} \\mathbf{S}_B\">$\\mathbf{S}_W^{-1} \\mathbf{S}_B$<\/span>. Le Lambda de Wilks est un indicateur inverse de taille d&rsquo;effet variant dans l&rsquo;intervalle <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span> :<\/p>\n<ul>\n<li>Une valeur proche de 1 signifie que la matrice intra-groupe absorbe la quasi-totalit\u00e9 de la dispersion totale, indiquant une absence d&rsquo;\u00e9cartement moyen des groupes (rejet impossible de <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span>).<\/li>\n<li>Une valeur proche de 0 traduit une dispersion inter-classes massivement dominante, signalant une discrimination multivari\u00e9e puissante des populations cliniques.<\/li>\n<\/ul>\n<p>Le Lambda de Wilks peut \u00eatre transform\u00e9 via l&rsquo;approximation de Rao en une statistique <span class=\"arabpsy-math-inline\" data-latex=\"F\">$F$<\/span> de Fisher-Snedecor permettant de d\u00e9river une <em>p<\/em>-valeur exacte. Si le test global s&rsquo;av\u00e8re significatif, une proc\u00e9dure d&rsquo;\u00e9limination s\u00e9quentielle permet d&rsquo;\u00e9valuer la significativit\u00e9 statistique de chaque racine canonique successive, d\u00e9terminant pr\u00e9cis\u00e9ment combien de dimensions factorielles discriminantes doivent \u00eatre conserv\u00e9es pour l&rsquo;interpr\u00e9tation. D&rsquo;autres crit\u00e8res multivari\u00e9s robustes \u2014 tels que la trace de Pillai-Bartlett, la trace d&rsquo;Hotelling-Lawley ou la plus grande racine de Roy \u2014 compl\u00e8tent avantageusement l&rsquo;analyse, la trace de Pillai \u00e9tant particuli\u00e8rement privil\u00e9gi\u00e9e en pr\u00e9sence de d\u00e9viations mod\u00e9r\u00e9es de l&rsquo;hypoth\u00e8se d&rsquo;homosc\u00e9dasticit\u00e9.<\/p>\n<p>Enfin, pour chaque axe discriminant canonique retenu, le <strong>coefficient de corr\u00e9lation canonique<\/strong> (<span class=\"arabpsy-math-inline\" data-latex=\"R_c\">$R_c$<\/span>) est calcul\u00e9. D\u00e9riv\u00e9 de la valeur propre associ\u00e9e selon la relation <span class=\"arabpsy-math-inline\" data-latex=\"R_c = \\sqrt{\\lambda_r \/ (1 + \\lambda_r)}\">$R_c = \\sqrt{\\lambda_r \/ (1 + \\lambda_r)}$<\/span>, il quantifie l&rsquo;intensit\u00e9 de l&rsquo;association statistique entre l&rsquo;axe psychom\u00e9trique latent et la partition des groupes cat\u00e9goriels, son carr\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"R_c^2\">$R_c^2$<\/span>) s&rsquo;interpr\u00e9tant directement comme la proportion de variance inter-classes totale expliqu\u00e9e par cette dimension.<\/p>\n<h2>9. Applications empiriques en psychologie et neurosciences cognitives<\/h2>\n<h3>9.1 Classification diff\u00e9rentielle des troubles psychopathologiques<\/h3>\n<p>Dans l&rsquo;exercice de la nosographie psychiatrique, le clinicien se heurte \u00e0 des recouvrements symptomatologiques consid\u00e9rables compliquant le diagnostic pr\u00e9coce. Par exemple, la diff\u00e9renciation entre le trouble d\u00e9pressif unipolaire caract\u00e9ris\u00e9, le trouble bipolaire en phase d\u00e9pressive inaugurale et les \u00e9tats d&rsquo;anxi\u00e9t\u00e9 g\u00e9n\u00e9ralis\u00e9e s\u00e9v\u00e8re pose des d\u00e9fis th\u00e9rapeutiques cruciaux, le recours inappropri\u00e9 aux antid\u00e9presseurs chez un patient bipolaire risquant d&rsquo;induire un virage maniaque iatrog\u00e8ne.<\/p>\n<p>L&rsquo;analyse discriminante lin\u00e9aire apporte une contribution d\u00e9cisive en exploitant simultan\u00e9ment des batteries h\u00e9t\u00e9rog\u00e8nes regroupant des inventaires d&rsquo;auto-\u00e9valuation de la r\u00e9gulation \u00e9motionnelle, des \u00e9preuves neuropsychologiques informatis\u00e9es (par exemple, des t\u00e2ches de flexibilit\u00e9 cognitive de type Wisconsin Card Sorting Test) et des biomarqueurs \u00e9lectrophysiologiques (tels que l&rsquo;amplitude de l&rsquo;onde P300 ou la variabilit\u00e9 du rythme cardiaque en condition de stress). L&rsquo;ADL comprime cette multidimensionnalit\u00e9 en isolant des signatures sp\u00e9cifiques :<\/p>\n<ul>\n<li>Le premier axe canonique discriminant s\u00e9pare typiquement le niveau d&rsquo;activation \u00e9motionnelle et d&rsquo;arousal physiologique global, dissociant les \u00e9tats anxieux purs des \u00e9tats de ralentissement psychomoteur d\u00e9pressif.<\/li>\n<li>Le second axe canonique capture les indices de labilit\u00e9 cognitive et d&rsquo;asym\u00e9trie frontale r\u00e9siduelle, op\u00e9rant la d\u00e9marcation fine entre les sous-types unipolaires et bipolaires.<\/li>\n<\/ul>\n<p>En pr\u00e9sence de comorbidit\u00e9s complexes, l&rsquo;ADL calcule les probabilit\u00e9s <em>a posteriori<\/em> d&rsquo;appartenance \u00e0 chaque entit\u00e9. Lorsqu&rsquo;un patient se positionne dans une zone d&rsquo;\u00e9qui-probabilit\u00e9 d\u00e9cisionnelle (par exemple <span class=\"arabpsy-math-inline\" data-latex=\"P(\\text{Unipolaire}) = 0{,}48\">$P(\\text{Unipolaire}) = 0{,}48$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"P(\\text{Bipolaire}) = 0{,}46\">$P(\\text{Bipolaire}) = 0{,}46$<\/span>), le mod\u00e8le formalise objectivement l&rsquo;ind\u00e9cision nosographique, incitant le coll\u00e8ge m\u00e9dical \u00e0 une surveillance clinique rapproch\u00e9e plut\u00f4t qu&rsquo;\u00e0 l&rsquo;application d&rsquo;un diagnostic unilat\u00e9ral rigide.<\/p>\n<h3>9.2 Orientation vocationnelle et profils de comp\u00e9tences<\/h3>\n<p>Dans le domaine de la psychologie du travail, de l&rsquo;ergonomie cognitive et de l&rsquo;orientation acad\u00e9mique, l&rsquo;\u00e9valuation des aptitudes et des int\u00e9r\u00eats vocationnels a historiquement constitu\u00e9 l&rsquo;un des premiers terrains d&rsquo;exp\u00e9rimentation de l&rsquo;ADL. L&rsquo;enjeu m\u00e9thodologique consiste \u00e0 mod\u00e9liser la mani\u00e8re dont une batterie de tests mesurant des dimensions cognitives diversifi\u00e9es \u2014 raisonnement logico-math\u00e9matique, fluidit\u00e9 verbale, vision spatiale tridimensionnelle, vitesse d&rsquo;int\u00e9gration perceptivo-motrice \u2014 pr\u00e9dit la trajectoire de r\u00e9ussite et d&rsquo;\u00e9panouissement au sein de fili\u00e8res d&rsquo;\u00e9tudes ou de secteurs professionnels distincts.<\/p>\n<p>L&rsquo;application d&rsquo;une ADL multiclasses sur des cohortes longitudinales de dipl\u00f4m\u00e9s permet d&rsquo;extraire la g\u00e9om\u00e9trie diff\u00e9rentielle sous-tendant les exigences cognitives sectorielles :<\/p>\n<ul>\n<li>Une premi\u00e8re fonction discriminante canonique oppose g\u00e9n\u00e9ralement un p\u00f4le de comp\u00e9tences d&rsquo;abstraction verbale et relationnelle (satur\u00e9 positivement par les fili\u00e8res juridiques, litt\u00e9raires et de psychologie clinique) \u00e0 un p\u00f4le de manipulation visuo-spatiale et de traitement formel quantitatif (satur\u00e9 par les cursus d&rsquo;ing\u00e9nierie m\u00e9canique et d&rsquo;architecture).<\/li>\n<li>Une seconde dimension factorielle discriminant isole le facteur de vitesse d&rsquo;ex\u00e9cution sous contrainte temporelle, diff\u00e9renciant par exemple les m\u00e9tiers de r\u00e9gulation du trafic a\u00e9rien des disciplines de recherche th\u00e9orique fondamentale.<\/li>\n<\/ul>\n<p>Les fonctions discriminantes ainsi standardis\u00e9es permettent de situer un jeune adulte sur cette carte des profils de comp\u00e9tences. Le calcul des distances statistiques le s\u00e9parant des centro\u00efdes professionnels renseigne sur ses correspondances aptitudinales objectives, fournissant un guide probabiliste dans le cadre de bilans d&rsquo;orientation sans jamais enfermer le sujet dans une trajectoire d\u00e9terministe.<\/p>\n<h3>9.3 D\u00e9tection pr\u00e9coce du d\u00e9clin cognitif en neuropsychologie<\/h3>\n<p>Face au vieillissement d\u00e9mographique plan\u00e9taire, le d\u00e9pistage pr\u00e9coce des processus neurod\u00e9g\u00e9n\u00e9ratifs constitue une priorit\u00e9 sanitaire majeure. La neuropsychologie du troisi\u00e8me \u00e2ge se donne pour mission de discriminer trois trajectoires cognitives aux fronti\u00e8res cliniques particuli\u00e8rement t\u00e9nues : le vieillissement cognitif normal (s\u00e9nescence saine), le d\u00e9ficit cognitif l\u00e9ger (<em>Mild Cognitive Impairment<\/em> &#8211; MCI, sous ses formes amn\u00e9sique et multi-domaine) et la maladie d&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Maladie_d%27Alzheimer\">Alzheimer<\/a> au stade prodromique ou l\u00e9ger.<\/p>\n<p>Une batterie d&rsquo;\u00e9valuation neuropsychologique standardis\u00e9e g\u00e9n\u00e8re couramment plus d&rsquo;une vingtaine d&rsquo;indices cliniques (scores de rappel imm\u00e9diat et diff\u00e9r\u00e9 au test de Grober-Buschke, interf\u00e9rences de Stroop, fluence phonologique, empans envers, Trail Making Test A et B). L&rsquo;ajustement d&rsquo;une ADL permet de hi\u00e9rarchiser math\u00e9matiquement le pouvoir d&rsquo;\u00e9vocation de chacun de ces marqueurs :<\/p>\n<ul>\n<li>Le mod\u00e8le d\u00e9montre de fa\u00e7on r\u00e9currente que le score de rappel diff\u00e9r\u00e9 sous indi\u00e7age s\u00e9mantique pr\u00e9sente le coefficient de fonction discriminante le plus puissant, suivi imm\u00e9diatement par le co\u00fbt d&rsquo;interf\u00e9rence attentionnelle.<\/li>\n<li>Les variables de m\u00e9moire \u00e0 court terme imm\u00e9diate (empans forwards), bien que d\u00e9grad\u00e9es chez les patients s\u00e9v\u00e8rement atteints, re\u00e7oivent des coefficients discriminants proches de z\u00e9ro, r\u00e9v\u00e9lant qu&rsquo;elles ne contribuent pas \u00e0 la d\u00e9marcation s\u00e9lective des stades prodromiques.<\/li>\n<\/ul>\n<p>De surcro\u00eet, le score discriminant individuel canonique issu de cette mod\u00e9lisation peut \u00eatre exploit\u00e9 longitudinalement comme un <em>biomarqueur cognitif composite continu<\/em>. En \u00e9valuant l&rsquo;\u00e9volution temporelle de la coordonn\u00e9e discriminante d&rsquo;un patient au fil de bilans r\u00e9p\u00e9t\u00e9s \u00e0 un an d&rsquo;intervalle, les neuropsychologues quantifient la vitesse de translation de son profil cognitif vers le centro\u00efde de la pathologie av\u00e9r\u00e9e, offrant un indicateur de suivi de la progression pathologique.<\/p>\n<h2>10. Variantes et extensions de l&rsquo;analyse discriminante<\/h2>\n<h3>10.1 L&rsquo;analyse discriminante quadratique<\/h3>\n<p>Lorsque le postulat fondamental d&rsquo;homosc\u00e9dasticit\u00e9 multivari\u00e9e est formellement viol\u00e9 \u2014 situation fr\u00e9quente o\u00f9 le test <span class=\"arabpsy-math-inline\" data-latex=\"M\">$M$<\/span> de Box rejette l&rsquo;\u00e9galit\u00e9 des matrices de covariance et o\u00f9 les dispersions intra-groupes s&rsquo;av\u00e8rent structurellement divergentes \u2014, l&rsquo;analyse discriminante lin\u00e9aire produit des partitions d&rsquo;espace biais\u00e9es. Il devient alors th\u00e9oriquement imp\u00e9ratif de rel\u00e2cher la contrainte d&rsquo;homog\u00e9n\u00e9it\u00e9 en adoptant l&rsquo;<strong>analyse discriminante quadratique<\/strong> (<em>Quadratic Discriminant Analysis<\/em> &#8211; QDA).<\/p>\n<p>Sous la formulation de la QDA, le mod\u00e8le estime une matrice de variance-covariance empirique sp\u00e9cifique et distincte <span class=\"arabpsy-math-inline\" data-latex=\"boldsymbol{\\Sigma}_k\">$boldsymbol{\\Sigma}_k$<\/span> pour chaque cat\u00e9gorie individuelle <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>. En r\u00e9introduisant ces matrices individualis\u00e9es au sein du d\u00e9veloppement bay\u00e9sien, le terme quadratique <span class=\"arabpsy-math-inline\" data-latex=\"-\\frac{1}{2} \\mathbf{x}^T boldsymbol{\\Sigma}_k^{-1} \\mathbf{x}\">$-\\frac{1}{2} \\mathbf{x}^T boldsymbol{\\Sigma}_k^{-1} \\mathbf{x}$<\/span> ne peut plus \u00eatre neutralis\u00e9 puisqu&rsquo;il d\u00e9pend d\u00e9sormais de l&rsquo;indice <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>. La fonction discriminante r\u00e9sultante pour la classe <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"Q_k(\\mathbf{x})\">$Q_k(\\mathbf{x})$<\/span>, conserve par cons\u00e9quent sa forme quadratique int\u00e9grale :<\/p>\n<p><span class=\"arabpsy-math-display\" data-latex=\"Q_k(\\mathbf{x}) = -\\frac{1}{2} \\ln |boldsymbol{\\Sigma}_k| - \\frac{1}{2} (\\mathbf{x} - boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}_k^{-1} (\\mathbf{x} - boldsymbol{\\mu}_k) + \\ln \\pi_k\">$$Q_k(\\mathbf{x}) = -\\frac{1}{2} \\ln |boldsymbol{\\Sigma}_k| &#8211; \\frac{1}{2} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k)^T boldsymbol{\\Sigma}_k^{-1} (\\mathbf{x} &#8211; boldsymbol{\\mu}_k) + \\ln \\pi_k$$<\/span><\/p>\n<p>Les fronti\u00e8res de d\u00e9cision d\u00e9limitant les groupes (<span class=\"arabpsy-math-inline\" data-latex=\"Q_k(\\mathbf{x}) = Q_l(\\mathbf{x})\">$Q_k(\\mathbf{x}) = Q_l(\\mathbf{x})$<\/span>) ne forment plus des hyperplans plats, mais des hypersurfaces courb\u00e9es du second degr\u00e9 (hyperparabolo\u00efdes, ellipses ferm\u00e9es ou selles hyperboliques). Cette plasticit\u00e9 g\u00e9om\u00e9trique permet \u00e0 la QDA de capturer des topologies complexes o\u00f9 un groupe clinique homog\u00e8ne est spatialement enclav\u00e9 au centre d&rsquo;une population de sujets sains pr\u00e9sentant une dispersion nettement plus large.<\/p>\n<p>Toutefois, cette flexibilit\u00e9 expressive a un co\u00fbt statistique lourd en mati\u00e8re de degr\u00e9s de libert\u00e9. Alors que l&rsquo;ADL n&rsquo;estime qu&rsquo;une seule matrice de covariance combin\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"p(p+1)\/2\">$p(p+1)\/2$<\/span> param\u00e8tres), la QDA exige l&rsquo;estimation distincte de <span class=\"arabpsy-math-inline\" data-latex=\"K\">$K$<\/span> matrices compl\u00e8tes, totalisant <span class=\"arabpsy-math-inline\" data-latex=\"K \\times [p(p+1)\/2]\">$K \\times [p(p+1)\/2]$<\/span> param\u00e8tres de dispersion. En pr\u00e9sence d&rsquo;\u00e9chantillons modestes, la QDA souffre d&rsquo;une variance d&rsquo;\u00e9chantillonnage \u00e9lev\u00e9e et d&rsquo;un risque majeur de surapprentissage, justifiant souvent de privil\u00e9gier la parcimonie robuste de l&rsquo;ADL malgr\u00e9 une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 mod\u00e9r\u00e9e.<\/p>\n<h3>10.2 L&rsquo;analyse discriminante r\u00e9gularis\u00e9e<\/h3>\n<p>Pour d\u00e9passer l&rsquo;opposition dichotomique rigide entre la parcimonie lin\u00e9aire de l&rsquo;ADL et la gourmandise param\u00e9trique quadratique de la QDA, Jerome Friedman a formalis\u00e9 en 1989 l&rsquo;<strong>analyse discriminante r\u00e9gularis\u00e9e<\/strong> (<em>Regularized Discriminant Analysis<\/em> &#8211; RDA). Cette extension computationnelle s&rsquo;av\u00e8re particuli\u00e8rement providentielle pour traiter des probl\u00e9matiques de haute dimensionnalit\u00e9 en neuro-imagerie ou en g\u00e9nomique comportementale, o\u00f9 le nombre de pr\u00e9dicteurs <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> approche ou d\u00e9passe la taille totale de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span>.<\/p>\n<p>La RDA s&rsquo;articule autour de l&rsquo;introduction de deux hyperparam\u00e8tres de r\u00e9gularisation et de r\u00e9tr\u00e9cissement statistique (<em>shrinkage<\/em>), continus et born\u00e9s sur l&rsquo;intervalle <span class=\"arabpsy-math-inline\" data-latex=\"[0, 1]\">$[0, 1]$<\/span>, not\u00e9s <span class=\"arabpsy-math-inline\" data-latex=\"lambda\">$lambda$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\gamma\">$\\gamma$<\/span> :<\/p>\n<ul>\n<li>Le premier param\u00e8tre, <span class=\"arabpsy-math-inline\" data-latex=\"lambda\">$lambda$<\/span>, orchestre un compromis param\u00e9trique direct entre la matrice de covariance individuelle d&rsquo;un groupe <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_k\">$\\mathbf{S}_k$<\/span> et la matrice combin\u00e9e intra-groupe globale <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}\">$\\mathbf{S}_{\\text{intra}}$<\/span> :\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{boldsymbol{\\Sigma}}_k(\\lambda) = (1 - \\lambda)\\mathbf{S}_k + \\lambda \\mathbf{S}_{\\text{intra}}\">$$\\hat{boldsymbol{\\Sigma}}_k(\\lambda) = (1 &#8211; \\lambda)\\mathbf{S}_k + \\lambda \\mathbf{S}_{\\text{intra}}$$<\/span><\/p>\n<p>    Lorsque <span class=\"arabpsy-math-inline\" data-latex=\"lambda = 0\">$lambda = 0$<\/span>, le mod\u00e8le \u00e9quivaut strictement \u00e0 la QDA ; lorsque <span class=\"arabpsy-math-inline\" data-latex=\"lambda = 1\">$lambda = 1$<\/span>, il converge identiquement vers l&rsquo;ADL classique.<\/li>\n<li>Le second param\u00e8tre, <span class=\"arabpsy-math-inline\" data-latex=\"\\gamma\">$\\gamma$<\/span>, induit un r\u00e9tr\u00e9cissement diagonal r\u00e9gularisant la matrice r\u00e9sultante vers un mod\u00e8le sph\u00e9rique isotrope, en la combinant avec une matrice identit\u00e9 pond\u00e9r\u00e9e par la variance univari\u00e9e moyenne <span class=\"arabpsy-math-inline\" data-latex=\"\\sigma_k^2 = \\frac{1}{p} \\text{Tr}(\\hat{boldsymbol{\\Sigma}}_k(\\lambda))\">$\\sigma_k^2 = \\frac{1}{p} \\text{Tr}(\\hat{boldsymbol{\\Sigma}}_k(\\lambda))$<\/span> :\n<p><span class=\"arabpsy-math-display\" data-latex=\"\\hat{boldsymbol{\\Sigma}}_k(\\lambda, \\gamma) = (1 - \\gamma)\\hat{boldsymbol{\\Sigma}}_k(\\lambda) + \\gamma \\sigma_k^2 \\mathbf{I}_p\">$$\\hat{boldsymbol{\\Sigma}}_k(\\lambda, \\gamma) = (1 &#8211; \\gamma)\\hat{boldsymbol{\\Sigma}}_k(\\lambda) + \\gamma \\sigma_k^2 \\mathbf{I}_p$$<\/span><\/p>\n<p>    Ce second palier de r\u00e9gularisation garantit l&rsquo;inversibilit\u00e9 stricte de la matrice de dispersion m\u00eame lorsque <span class=\"arabpsy-math-inline\" data-latex=\"p gg N\">$p gg N$<\/span>, \u00e9liminant le ph\u00e9nom\u00e8ne de singularit\u00e9 matricielle.<\/li>\n<\/ul>\n<p>Le calibrage optimal de ce doublet d&rsquo;hyperparam\u00e8tres <span class=\"arabpsy-math-inline\" data-latex=\"(\\lambda, \\gamma)\">$(\\lambda, \\gamma)$<\/span> est conduit empiriquement par recherche sur grille (<em>grid search<\/em>) crois\u00e9e optimisant le taux d&rsquo;exactitude ou le score de log-vraisemblance. La RDA fournit ainsi un estimateur adaptatif stabilis\u00e9, capitalisant sur la robustesse de l&rsquo;ADL tout en int\u00e9grant des degr\u00e9s contr\u00f4l\u00e9s de courbure quadratique.<\/p>\n<h3>10.3 Approches non lin\u00e9aires et m\u00e9thodes \u00e0 noyaux<\/h3>\n<p>Lorsque les structures sous-jacentes s\u00e9parant les groupes comportementaux adoptent des g\u00e9om\u00e9tries hautement non lin\u00e9aires \u2014 telles que des cercles concentriques ou des vari\u00e9t\u00e9s spiral\u00e9es \u2014, aucune fronti\u00e8re lin\u00e9aire ni m\u00eame quadratique ne peut op\u00e9rer une discrimination satisfaisante au sein de l&rsquo;espace d&rsquo;origine <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbb{R}^p\">$\\mathbb{R}^p$<\/span>. Pour surmonter cette limitation g\u00e9om\u00e9trique sans sacrifier l&rsquo;ossature analytique de Fisher, Mika et ses collaborateurs ont formalis\u00e9 l&rsquo;<strong>analyse discriminante \u00e0 noyaux<\/strong> (<em>Kernel Fisher Discriminant Analysis<\/em> &#8211; KFDA).<\/p>\n<p>Cette approche applique l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Astuce_du_noyau\">astuce du noyau<\/a> (<em>kernel trick<\/em>), concept popularis\u00e9 par les machines \u00e0 vecteurs de support (SVM). Les donn\u00e9es d&rsquo;origine sont projet\u00e9es de mani\u00e8re implicite via une transformation non lin\u00e9aire <span class=\"arabpsy-math-inline\" data-latex=\"\\phi(\\mathbf{x})\">$\\phi(\\mathbf{x})$<\/span> vers un espace de Hilbert de dimension potentiellement infinie, au sein duquel les classes redeviennent s\u00e9parables par un hyperplan plat. Gr\u00e2ce aux propri\u00e9t\u00e9s de reproducing kernel Hilbert spaces (RKHS), tous les produits scalaires dans cet espace d&rsquo;arriv\u00e9e s&rsquo;expriment au moyen d&rsquo;une fonction de noyau <span class=\"arabpsy-math-inline\" data-latex=\"k(\\mathbf{x}_i, \\mathbf{x}_j) = \\langle \\phi(\\mathbf{x}_i), \\phi(\\mathbf{x}_j) \\rangle\">$k(\\mathbf{x}_i, \\mathbf{x}_j) = \\langle \\phi(\\mathbf{x}_i), \\phi(\\mathbf{x}_j) \\rangle$<\/span>, sans jamais n\u00e9cessiter le calcul explicite des coordonn\u00e9es transform\u00e9es :<\/p>\n<ul>\n<li>Le noyau gaussien \u00e0 base radiale (<em>RBF &#8211; Radial Basis Function<\/em>) : <span class=\"arabpsy-math-inline\" data-latex=\"k(\\mathbf{x}_i, \\mathbf{x}_j) = \\exp(-\\gamma |\\mathbf{x}_i - \\mathbf{x}_j|^2)\">$k(\\mathbf{x}_i, \\mathbf{x}_j) = \\exp(-\\gamma |\\mathbf{x}_i &#8211; \\mathbf{x}_j|^2)$<\/span><\/li>\n<li>Le noyau polynomial d&rsquo;ordre <span class=\"arabpsy-math-inline\" data-latex=\"d\">$d$<\/span> : <span class=\"arabpsy-math-inline\" data-latex=\"k(\\mathbf{x}_i, \\mathbf{x}_j) = (\\mathbf{x}_i^T \\mathbf{x}_j + c)^d\">$k(\\mathbf{x}_i, \\mathbf{x}_j) = (\\mathbf{x}_i^T \\mathbf{x}_j + c)^d$<\/span><\/li>\n<\/ul>\n<p>Si la KFDA permet de mod\u00e9liser des fronti\u00e8res de s\u00e9paration d\u00e9cisionnelles arbitrairement complexes rivalisant avec les algorithmes d&rsquo;apprentissage profond, elle induit une rupture \u00e9pist\u00e9mologique avec l&rsquo;ADL classique. Les fonctions discriminantes n&rsquo;\u00e9tant plus des combinaisons lin\u00e9aires directes des variables d&rsquo;origine, l&rsquo;interpr\u00e9tabilit\u00e9 psychom\u00e9trique des coefficients individuels est totalement oblit\u00e9r\u00e9e. D\u00e8s lors, la KFDA rel\u00e8ve purement de l&rsquo;apprentissage automatique pr\u00e9dictif (IA de type bo\u00eete noire), s&rsquo;\u00e9loignant de la double mission explicative-descriptive qui fait la force de l&rsquo;ADL conventionnelle.<\/p>\n<h2>11. Impl\u00e9mentation technique et protocole d&rsquo;analyse sous R et Python<\/h2>\n<h3>11.1 Pr\u00e9paration et nettoyage des donn\u00e9es multivari\u00e9es<\/h3>\n<p>L&rsquo;application rigoureuse d&rsquo;une analyse discriminante lin\u00e9aire impose un pipeline pr\u00e9liminaire de pr\u00e9traitement des donn\u00e9es, garantissant que les estimateurs matriciels ne soient pas compromis par des artefacts de mesure. Ce protocole syst\u00e9matique comprend quatre \u00e9tapes coordonn\u00e9es :<\/p>\n<ul>\n<li><strong>Le diagnostic et le traitement des valeurs aberrantes multivari\u00e9es :<\/strong> Une observation aberrante isol\u00e9e sur une seule dimension peut passer inaper\u00e7ue lors d&rsquo;une inspection univari\u00e9e, tout en g\u00e9n\u00e9rant un levier d\u00e9stabilisateur majeur dans l&rsquo;espace multidimensionnel. L&rsquo;analyste doit calculer syst\u00e9matiquement le carr\u00e9 des distances de Mahalanobis pour chaque sujet par rapport au centro\u00efde de sa classe et identifier les participants dont la distance exc\u00e8de le quantile critique de la distribution th\u00e9orique du Chi-deux au seuil <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0{,}001\">$\\alpha = 0{,}001$<\/span>. L&rsquo;exclusion motiv\u00e9e ou la winsorisation multivari\u00e9e de ces observations est imp\u00e9rative.<\/li>\n<li><strong>L&rsquo;imputation des donn\u00e9es manquantes :<\/strong> Les batteries psychom\u00e9triques exhaustives comportent fr\u00e9quemment des omissions ponctuelles d&rsquo;items. L&rsquo;exclusion pure et simple de tout dossier incomplet (<em>listwise deletion<\/em>) induit une r\u00e9duction drastique de la puissance statistique et un biais de s\u00e9lection potentiel. Le recours \u00e0 des algorithmes d&rsquo;imputation multiple par \u00e9quations cha\u00een\u00e9es (MICE) ou par l&rsquo;algorithme des <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> plus proches voisins (k-NN) sur les pr\u00e9dicteurs continus pr\u00e9serve la repr\u00e9sentativit\u00e9 de la cohorte sans alt\u00e9rer la matrice de covariance.<\/li>\n<li><strong>La standardisation centr\u00e9e-r\u00e9duite :<\/strong> Bien que l&rsquo;ADL math\u00e9matique soit formellement invariante \u00e0 tout changement d&rsquo;\u00e9chelle univari\u00e9 lin\u00e9aire (les fonctions de projection s&rsquo;adaptant proportionnellement), la standardisation pr\u00e9alable de tous les pr\u00e9dicteurs (centrage sur z\u00e9ro et r\u00e9duction \u00e0 une variance unitaire : <span class=\"arabpsy-math-inline\" data-latex=\"z = (x - \\mu)\/\\sigma\">$z = (x &#8211; \\mu)\/\\sigma$<\/span>) constitue une bonne pratique universelle. Elle permet l&rsquo;interpr\u00e9tation comparative directe des coefficients de discrimination canoniques bruts sans distorsion induite par l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des amplitudes m\u00e9triques initiales.<\/li>\n<li><strong>Le partitionnement stratifi\u00e9 \u00e9tanche :<\/strong> Avant tout calcul statistique, le jeu de donn\u00e9es doit \u00eatre scind\u00e9 en un ensemble d&rsquo;entra\u00eenement (<em>training set<\/em>, typiquement 70 % \u00e0 80 % de la cohorte) d\u00e9di\u00e9 \u00e0 l&rsquo;ajustement du mod\u00e8le et de ses seuils, et un ensemble de test ind\u00e9pendant (<em>test set<\/em>, 20 % \u00e0 30 %) r\u00e9serv\u00e9 exclusivement \u00e0 l&rsquo;\u00e9valuation finale des performances de g\u00e9n\u00e9ralisation pr\u00e9dictive.<\/li>\n<\/ul>\n<h3>11.2 D\u00e9ploiement sous R avec la biblioth\u00e8que MASS<\/h3>\n<p>Au sein de l&rsquo;environnement de programmation statistique R, l&rsquo;impl\u00e9mentation de r\u00e9f\u00e9rence de l&rsquo;analyse discriminante lin\u00e9aire est encapsul\u00e9e dans le package <a rel=\"noopener noreferrer\" href=\"https:\/\/cran.r-project.org\/package=MASS\">MASS<\/a> (<em>Modern Applied Statistics with S<\/em>), d\u00e9velopp\u00e9 par Venables et Ripley. Ce package fournit la fonction cardinale <code>lda()<\/code>, qui impl\u00e9mente l&rsquo;estimation non par inversion directe de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_{\\text{intra}}\">$\\mathbf{S}_{\\text{intra}}$<\/span>, mais par d\u00e9composition en valeurs singuli\u00e8res (SVD) de la matrice des donn\u00e9es centr\u00e9es, une approche computationnelle offrant une stabilit\u00e9 num\u00e9rique sup\u00e9rieure.<\/p>\n<p>Le protocole analytique sous R s&rsquo;articule autour des syntaxes et d\u00e9marches suivantes :<\/p>\n<ul>\n<li><strong>Ajustement du mod\u00e8le de base :<\/strong> La syntaxe mod\u00e9lise la variable cat\u00e9gorielle d\u00e9pendante en fonction de la matrice des scores continus, en autorisant la sp\u00e9cification manuelle du vecteur des probabilit\u00e9s <em>a priori<\/em> via l&rsquo;argument <code>prior<\/code> (permettant d&rsquo;ajuster les pr\u00e9valences \u00e9pid\u00e9miologiques) :<br \/>\n    <br \/><code>modele_lda &lt;- lda(Diagnostic ~ Memoire + Attention + Vitesse + Flexibilite, data = donnees_train, prior = c(0.70, 0.20, 0.10))<\/code><\/li>\n<li><strong>Extraction des param\u00e8tres canoniques :<\/strong> L&rsquo;affichage du mod\u00e8le r\u00e9v\u00e8le les moyennes des groupes au sein de chaque variable (centro\u00efdes empiriques), ainsi que les coefficients des discriminants lin\u00e9aires (<em>coefficients of linear discriminants<\/em>). Ces coefficients repr\u00e9sentent les pond\u00e9rations vectorielles <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{w}\">$\\mathbf{w}$<\/span> associ\u00e9es \u00e0 chaque fonction canonique. La fonction <code>summary()<\/code> permet d&rsquo;extraire les proportions de s\u00e9paration inter-classes captur\u00e9es par chaque axe factoriel successif.<\/li>\n<li><strong>Visualisation graphique de l&rsquo;espace factoriel :<\/strong> L&rsquo;appel de la m\u00e9thode graphique g\u00e9n\u00e9rique <code>plot(modele_lda, col = as.numeric(donnees_train$Diagnostic))<\/code> produit une cartographie des observations projet\u00e9es sur les deux premi\u00e8res fonctions discriminantes canoniques <span class=\"arabpsy-math-inline\" data-latex=\"LD_1\">$LD_1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"LD_2\">$LD_2$<\/span>. Pour des visualisations enrichies d&rsquo;ellipses d&rsquo;\u00e9quiprobabilit\u00e9 \u00e0 95 %, l&rsquo;int\u00e9gration des coordonn\u00e9es transform\u00e9es au sein de l&rsquo;\u00e9cosyst\u00e8me <code>ggplot2<\/code> (via l&rsquo;instruction <code>predict(modele_lda)$x<\/code>) permet d&rsquo;exporter des cartographies pr\u00eates pour publication.<\/li>\n<li><strong>Inf\u00e9rence pr\u00e9dictive sur de nouvelles observations :<\/strong> L&rsquo;application de la fonction <code>predict(modele_lda, newdata = donnees_test)<\/code> renvoie une liste hi\u00e9rarchis\u00e9e contenant trois sorties cardinales :\n<ul>\n<li><code>$class<\/code> : l&rsquo;assignation cat\u00e9gorielle finale issue de la maximisation bay\u00e9sienne.<\/li>\n<li><code>$posterior<\/code> : la matrice des probabilit\u00e9s <em>a posteriori<\/em> individuelles pour chaque classe, permettant de quantifier le degr\u00e9 de certitude du mod\u00e8le.<\/li>\n<li><code>$x<\/code> : les coordonn\u00e9es g\u00e9om\u00e9triques de chaque nouveau patient dans le sous-espace factoriel canonique.<\/li>\n<\/ul>\n<\/li>\n<\/ul>\n<h3>11.3 Mise en \u0153uvre sous Python avec Scikit-Learn<\/h3>\n<p>Au sein de l&rsquo;\u00e9cosyst\u00e8me d&rsquo;apprentissage automatique en langage Python, le module <a rel=\"noopener noreferrer\" href=\"https:\/\/scikit-learn.org\/stable\/modules\/generated\/sklearn.discriminant_analysis.LinearDiscriminantAnalysis.html\">LinearDiscriminantAnalysis<\/a> de la biblioth\u00e8que <code>scikit-learn<\/code> (sous le sous-module <code>sklearn.discriminant_analysis<\/code>) constitue l&rsquo;infrastructure logicielle standard. Cette classe g\u00e8re de mani\u00e8re int\u00e9gr\u00e9e tant la r\u00e9duction de dimension supervis\u00e9e que la classification bay\u00e9sienne multiclasses.<\/p>\n<p>Le d\u00e9ploiement industriel d&rsquo;un pipeline d&rsquo;analyse sous Python requiert le param\u00e9trage rigoureux des composantes de la classe :<\/p>\n<ul>\n<li><strong>S\u00e9lection des solveurs math\u00e9matiques (param\u00e8tre <code>solver<\/code>) :<\/strong>\n<ul>\n<li><code>solver='svd'<\/code> : solveur par d\u00e9faut, fond\u00e9 sur la d\u00e9composition en valeurs singuli\u00e8res. Il ne requiert pas le calcul direct de la matrice de covariance, ce qui le rend exceptionnellement rapide et stable pour des jeux de donn\u00e9es comportant un grand nombre d&rsquo;individus. N\u00e9anmoins, il ne permet pas le r\u00e9glage du param\u00e8tre de shrinkage.<\/li>\n<li><code>solver='lsqr'<\/code> : solveur bas\u00e9 sur les moindres carr\u00e9s r\u00e9gularis\u00e9s. Il s&rsquo;av\u00e8re hautement efficient et supporte l&rsquo;introduction dynamique du r\u00e9tr\u00e9cissement statistique de Ledoit-Wolf via l&rsquo;argument <code>shrinkage='auto'<\/code>.<\/li>\n<li><code>solver='eigen'<\/code> : solveur mobilisant la d\u00e9composition spectrale classique des valeurs propres de <span class=\"arabpsy-math-inline\" data-latex=\"\\mathbf{S}_W^{-1}\\mathbf{S}_B\">$\\mathbf{S}_W^{-1}\\mathbf{S}_B$<\/span>. Id\u00e9al lorsque l&rsquo;on souhaite extraire et analyser explicitement les racines et valeurs propres associ\u00e9es aux axes canoniques de Fisher.<\/li>\n<\/ul>\n<\/li>\n<li><strong>Int\u00e9gration au sein d&rsquo;un <code>Pipeline<\/code> Scikit-Learn :<\/strong> Afin de pr\u00e9venir rigoureusement toute fuite d&rsquo;information statistique entre l&rsquo;apprentissage et le test (<em>data leakage<\/em>), l&rsquo;ADL est encapsul\u00e9e au sein d&rsquo;une cha\u00eene s\u00e9quentielle int\u00e9grant la standardisation :<br \/>\n    <br \/><code>pipeline = Pipeline([('scaler', StandardScaler()), ('lda', LinearDiscriminantAnalysis(solver='lsqr', shrinkage='auto'))])<\/code><\/li>\n<li><strong>Rapport d&rsquo;\u00e9valuation de la classification :<\/strong> L&rsquo;\u00e9valuation finale des performances de test mobilise les fonctions de <code>sklearn.metrics<\/code>, notamment <code>classification_report()<\/code> qui calcule instantan\u00e9ment pour chaque classe la pr\u00e9cision, le rappel (sensibilit\u00e9) et le score F1 macro et pond\u00e9r\u00e9, ainsi que <code>confusion_matrix()<\/code>. Pour les probl\u00e9matiques de r\u00e9duction de dimensionnalit\u00e9, la m\u00e9thode <code>lda.transform(X)<\/code> g\u00e9n\u00e8re directement le tableau des coordonn\u00e9es projet\u00e9es sur les composantes canoniques les plus informatives.<\/li>\n<\/ul>\n<h2>12. Bonnes pratiques m\u00e9thodologiques, \u00e9cueils et perspectives de recherche<\/h2>\n<h3>12.1 Robustesse aux violations mod\u00e9r\u00e9es des postulats<\/h3>\n<p>L&rsquo;application rigoureuse de l&rsquo;ADL en contexte scientifique requiert un discernement lucide quant \u00e0 la tol\u00e9rance r\u00e9elle de la m\u00e9thode face aux in\u00e9vitables imperfections des donn\u00e9es empiriques. La litt\u00e9rature statistique de simulation de Monte-Carlo a largement explor\u00e9 le comportement de l&rsquo;ADL lors de violations mod\u00e9r\u00e9es de ses postulats th\u00e9oriques :<\/p>\n<ul>\n<li><strong>Tol\u00e9rance \u00e0 la non-normalit\u00e9 univari\u00e9e et multivari\u00e9e :<\/strong> L&rsquo;ADL manifeste une remarquable robustesse face aux d\u00e9viations sym\u00e9triques de la normalit\u00e9 (telles que des distributions platykurtiques ou l\u00e9g\u00e8rement leptokurtiques). En vertu du th\u00e9or\u00e8me central limite multivari\u00e9, les estimateurs des moyennes de classes <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{\\mathbf{x}}_k\">$\\bar{\\mathbf{x}}_k$<\/span> convergent rapidement vers la normalit\u00e9 d\u00e8s lors que les effectifs de chaque groupe d\u00e9passent <span class=\"arabpsy-math-inline\" data-latex=\"N_k ge 20\">$N_k ge 20$<\/span> \u00e0 30 individus. En revanche, des asym\u00e9tries prononc\u00e9es (<em>skewness<\/em> marqu\u00e9) et divergentes entre groupes faussent les fronti\u00e8res d\u00e9cisionnelles, rendant obligatoire l&rsquo;application pr\u00e9alable d&rsquo;une transformation normalisatrice de Box-Cox.<\/li>\n<li><strong>Donn\u00e9es issues d&rsquo;\u00e9chelles de Likert ordinales :<\/strong> Bien que l&rsquo;ADL postule formellement des variables continues, son application sur des scores composites issus de batteries d&rsquo;\u00e9chelles de Likert (lorsque les scores totaux agr\u00e8gent au moins 5 \u00e0 7 niveaux et se rapprochent d&rsquo;un continuum) n&rsquo;alt\u00e8re que tr\u00e8s marginalement les taux d&rsquo;erreur de classification par rapport \u00e0 des mod\u00e8les polytomiques non lin\u00e9aires plus lourds.<\/li>\n<li><strong>Limites face \u00e0 l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 de covariance :<\/strong> La violation de l&rsquo;homosc\u00e9dasticit\u00e9 constitue la faiblesse la plus dommageable de l&rsquo;ADL. Les simulations d\u00e9montrent que si le rapport des plus grandes variances de classes sur les plus faibles n&rsquo;exc\u00e8de pas un facteur de 2 \u00e0 3, et sous r\u00e9serve que les tailles de groupes soient strictement \u00e9quilibr\u00e9es (<span class=\"arabpsy-math-inline\" data-latex=\"N_1 \\approx N_2 dots \\approx N_K\">$N_1 \\approx N_2 dots \\approx N_K$<\/span>), l&rsquo;impact de l&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 sur le taux d&rsquo;erreur global demeure n\u00e9gligeable. En revanche, la combinaison d\u00e9l\u00e9t\u00e8re d&rsquo;une h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 marqu\u00e9e et de tailles d&rsquo;\u00e9chantillons fortement in\u00e9gales produit des biais d&rsquo;assignation majeurs : l&rsquo;ADL favorise massivement la classe dont la variance est la plus ample au d\u00e9triment de la classe compacte.<\/li>\n<\/ul>\n<p>Face \u00e0 ces configurations limites, le statisticien avis\u00e9 consolidera ses inf\u00e9rences en impl\u00e9mentant des proc\u00e9dures de r\u00e9\u00e9chantillonnage par <strong>bootstrap non param\u00e9trique<\/strong> (typiquement 2 000 r\u00e9plications). En g\u00e9n\u00e9rant des distributions empiriques des coefficients discriminants sans postulat distributionnel formel, le bootstrap permet de calculer des intervalles de confiance \u00e0 95 % corrig\u00e9s des biais, attestant de la solidit\u00e9 des variables discriminantes identifi\u00e9es.<\/p>\n<h3>12.2 Pi\u00e8ges d&rsquo;interpr\u00e9tation des coefficients de fonction discriminante<\/h3>\n<p>L&rsquo;erreur m\u00e9thodologique la plus insidieuse lors de la restitution des r\u00e9sultats d&rsquo;une analyse discriminante consiste \u00e0 interpr\u00e9ter na\u00efvement la magnitude brute des coefficients des fonctions discriminantes comme une mesure directe de l&rsquo;importance clinique causale ou diff\u00e9rentielle d&rsquo;un pr\u00e9dicteur. Trois artefacts statistiques doivent imp\u00e9rativement \u00eatre ma\u00eetris\u00e9s :<\/p>\n<ul>\n<li><strong>La disparit\u00e9 des \u00e9chelles de mesure initiales :<\/strong> Les coefficients canoniques bruts d\u00e9pendent intrins\u00e8quement de l&rsquo;unit\u00e9 physique ou m\u00e9trique de chaque variable. Un score mesur\u00e9 sur une amplitude de 0 \u00e0 1 recevra artificiellement un coefficient brut dix fois sup\u00e9rieur \u00e0 celui d&rsquo;un test mesur\u00e9 de 0 \u00e0 100, pour un pouvoir discriminant strictement identique. Seuls les <strong>coefficients canoniques standardis\u00e9s<\/strong> (obtenus en multipliant le coefficient brut par l&rsquo;\u00e9cart-type intra-groupe de la variable correspondante) permettent une comparaison relative rigoureuse du pouvoir discriminant au sein du mod\u00e8le.<\/li>\n<li><strong>L&rsquo;illusion causale induite par la multicolin\u00e9arit\u00e9 :<\/strong> \u00c0 l&rsquo;instar de ce qui s&rsquo;observe dans les mod\u00e8les de r\u00e9gression multiple, la pr\u00e9sence de corr\u00e9lations mod\u00e9r\u00e9es \u00e0 fortes entre plusieurs \u00e9preuves psychom\u00e9triques masque la contribution de pr\u00e9dicteurs cruciaux. Si deux tests d&rsquo;attention partagent 70 % de variance commune, le mod\u00e8le peut attribuer un coefficient \u00e9lev\u00e9 \u00e0 l&rsquo;un et un coefficient proche de z\u00e9ro \u00e0 l&rsquo;autre, non pas parce que ce dernier est d\u00e9pourvu de valeur discriminante clinique, mais parce que son apport d&rsquo;information diff\u00e9rentielle est redondant avec le premier.<\/li>\n<li><strong>La distinction entre coefficients standardis\u00e9s et matrices de structure :<\/strong> Pour d\u00e9m\u00ealer cet artefact, l&rsquo;analyste doit imp\u00e9rativement rapporter et examiner la <em>matrice de structure<\/em> (\u00e9galement appel\u00e9e corr\u00e9lations canoniques intra-groupes ou <em>structure loadings<\/em>). Ces indicateurs quantifient la corr\u00e9lation bivari\u00e9e directe entre chaque pr\u00e9dicteur individuel et l&rsquo;axe latent discriminant, ind\u00e9pendamment des autres variables du mod\u00e8le. Une variable peut avoir un coefficient standardis\u00e9 quasi-nul dans la fonction pr\u00e9dictive en raison de sa redondance, tout en affichant une saturation structurelle tr\u00e8s \u00e9lev\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"r &gt; 0{,}60\">$r > 0{,}60$<\/span>) avec la dimension diff\u00e9rentielle globale, attestant de sa centralit\u00e9 clinique dans le tableau pathologique.<\/li>\n<\/ul>\n<p>Conform\u00e9ment aux normes r\u00e9dactionnelles de l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/apastyle.apa.org\/\">American Psychological Association<\/a> (APA 7th edition), un rapport d&rsquo;analyse discriminante lin\u00e9aire exemplaire doit syst\u00e9matiquement consigner : les tests de normalit\u00e9 multivari\u00e9e et d&rsquo;homog\u00e9n\u00e9it\u00e9 (test <span class=\"arabpsy-math-inline\" data-latex=\"M\">$M$<\/span> de Box), les corr\u00e9lations de la matrice de structure conjointement aux coefficients standardis\u00e9s, le Lambda de Wilks associ\u00e9 \u00e0 son test de significativit\u00e9 statistique <span class=\"arabpsy-math-inline\" data-latex=\"F\">$F$<\/span>, les valeurs propres de chaque racine canonique, les pourcentages de variance expliqu\u00e9e et le coefficient Kappa de Cohen issu d&rsquo;une validation crois\u00e9e rigoureuse.<\/p>\n<h3>12.3 Int\u00e9gration de l&rsquo;ADL dans l&rsquo;apprentissage automatique moderne<\/h3>\n<p>Loin d&rsquo;\u00eatre rel\u00e9gu\u00e9e au statut de relique historique par l&rsquo;\u00e9mergence des technologies d&rsquo;intelligence artificielle contemporaines, l&rsquo;analyse discriminante lin\u00e9aire conna\u00eet un renouveau remarquable au sein des architectures d&rsquo;apprentissage automatique avanc\u00e9es. Sa rapidit\u00e9 computationnelle, sa faible consommation en ressources de calcul et, par-dessus tout, sa formulation analytique rigoureusement transparente en font un composant architectural de premier plan dans le paradigme de l&rsquo;<strong>intelligence artificielle explicable<\/strong> (<em>XAI &#8211; Explainable Artificial Intelligence<\/em>).<\/p>\n<p>Dans les protocoles de mod\u00e9lisation contemporains, l&rsquo;ADL est fr\u00e9quemment mobilis\u00e9e selon des logiques d&rsquo;hybridation m\u00e9thodologique :<\/p>\n<ul>\n<li><strong>Extracteur d&rsquo;attributs supervis\u00e9 en amont de classificateurs non lin\u00e9aires :<\/strong> Plut\u00f4t que de soumettre directement des milliers de descripteurs non filtr\u00e9s \u00e0 un classificateur complexe (tel qu&rsquo;une machine \u00e0 vecteurs de support \u00e0 noyau RBF ou un algorithme de Gradient Boosting de type XGBoost), l&rsquo;ADL est d\u00e9ploy\u00e9e en phase amont pour projeter les donn\u00e9es vers l&rsquo;espace discriminant canonique optimal de Fisher. Ce sous-espace compact de dimension <span class=\"arabpsy-math-inline\" data-latex=\"K - 1\">$K &#8211; 1$<\/span> sert ensuite d&rsquo;entr\u00e9e hautement purifi\u00e9e pour le classificateur non lin\u00e9aire, conjuguant l&rsquo;acc\u00e9l\u00e9ration de convergence, la r\u00e9duction de variance et la capture des courbures r\u00e9siduelles.<\/li>\n<li><strong>R\u00e9gularisation parcimonieuse des mod\u00e8les d&rsquo;ensemble :<\/strong> Dans les architectures d&#8217;empilement (<em>stacking ensembles<\/em>), les pr\u00e9dictions probabilistes g\u00e9n\u00e9r\u00e9es par une ADL lin\u00e9aire robuste sont fusionn\u00e9es avec celles issues de mod\u00e8les non param\u00e9triques flexibles. L&rsquo;ADL agit comme une ancre bay\u00e9sienne stabilisatrice, emp\u00eachant les mod\u00e8les d&rsquo;ensemble de surajuster le bruit idiosyncrasique de l&rsquo;\u00e9chantillon d&rsquo;apprentissage.<\/li>\n<li><strong>Interpr\u00e9tation des espaces latents dans l&rsquo;apprentissage profond :<\/strong> Dans les r\u00e9seaux de neurones convolutionnels ou les auto-encodeurs variationnels appliqu\u00e9s \u00e0 la neuro-imagerie m\u00e9dicale, l&rsquo;ADL est appliqu\u00e9e sur l&rsquo;avant-derni\u00e8re couche dense du r\u00e9seau (l&rsquo;espace de repr\u00e9sentations latentes). Cela permet de v\u00e9rifier si les repr\u00e9sentations g\u00e9om\u00e9triques abstraites forg\u00e9es de mani\u00e8re non supervis\u00e9e par les neurones profonds s&rsquo;alignent naturellement sur des axes lin\u00e9aires s\u00e9parateurs discriminants conformes \u00e0 la clinique m\u00e9dicale.<\/li>\n<\/ul>\n<p>Ainsi, pr\u00e8s d&rsquo;un si\u00e8cle apr\u00e8s sa formulation princeps par Ronald Fisher, l&rsquo;analyse discriminante lin\u00e9aire s&rsquo;affirme plus que jamais comme un pilier m\u00e9thodologique indispensable : un pont entre la mod\u00e9lisation statistique inf\u00e9rentielle classique et la puissance pr\u00e9dictive de l&rsquo;apprentissage automatique contemporain, garantissant que l&rsquo;investigation des structures quantitatives du comportement humain demeure fermement ancr\u00e9e dans les principes de rigueur, de parcimonie et d&rsquo;interpr\u00e9tabilit\u00e9 math\u00e9matique.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Albert, A., &amp; Anderson, J. A. (1984). On the existence of maximum likelihood estimates in logistic regression models. <em>Biometrika<\/em>, 71(1), 1\u201310. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/71.1.1\">https:\/\/doi.org\/10.1093\/biomet\/71.1.1<\/a><\/li>\n<li>American Psychological Association. (2020). <em>Publication manual of the American Psychological Association<\/em> (7th ed.). American Psychological Association. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1037\/0000165-000\">https:\/\/doi.org\/10.1037\/0000165-000<\/a><\/li>\n<li>Efron, B. (1975). The efficiency of logistic regression compared to normal discriminant analysis. <em>Journal of the American Statistical Association<\/em>, 70(352), 892\u2013898. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1975.10480319\">https:\/\/doi.org\/10.1080\/01621459.1975.10480319<\/a><\/li>\n<li>Fisher, R. A. (1936). The use of multiple measurements in taxonomic problems. <em>Annals of Eugenics<\/em>, 7(2), 179\u2013188. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.1469-1809.1936.tb02137.x\">https:\/\/doi.org\/10.1111\/j.1469-1809.1936.tb02137.x<\/a><\/li>\n<li>Friedman, J. H. (1989). Regularized discriminant analysis. <em>Journal of the American Statistical Association<\/em>, 84(405), 165\u2013175. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1989.10478752\">https:\/\/doi.org\/10.1080\/01621459.1989.10478752<\/a><\/li>\n<li>Hand, D. J., &amp; Till, R. J. (2001). A simple generalisation of the area under the ROC curve for multiple class classification problems. <em>Machine Learning<\/em>, 45(2), 171\u2013186. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1023\/A:1010920819831\">https:\/\/doi.org\/10.1023\/A:1010920819831<\/a><\/li>\n<li>Hastie, T., Tibshirani, R., &amp; Friedman, J. (2009). <em>The elements of statistical learning: Data mining, inference, and prediction<\/em> (2nd ed.). Springer. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-0-387-84858-7\">https:\/\/doi.org\/10.1007\/978-0-387-84858-7<\/a><\/li>\n<li>Huberty, C. J., &amp; Olejnik, S. (2006). <em>Applied MANOVA and discriminant analysis<\/em> (2nd ed.). John Wiley &amp; Sons. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/047178947X\">https:\/\/doi.org\/10.1002\/047178947X<\/a><\/li>\n<li>Klecka, W. R. (1980). <em>Discriminant analysis<\/em>. SAGE Publications. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.4135\/9781412983990\">https:\/\/doi.org\/10.4135\/9781412983990<\/a><\/li>\n<li>Mardia, K. V. (1970). Measures of multivariate skewness and kurtosis with applications. <em>Biometrika<\/em>, 57(3), 519\u2013530. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/57.3.519\">https:\/\/doi.org\/10.1093\/biomet\/57.3.519<\/a><\/li>\n<li>Mika, S., Ratsch, G., Weston, J., Scholkopf, B., &amp; Muller, K. R. (1999). Fisher discriminant analysis with kernels. In <em>Neural Networks for Signal Processing IX: Proceedings of the 1999 IEEE Signal Processing Society Workshop<\/em> (pp. 41\u201348). IEEE. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1109\/NNSP.1999.788121\">https:\/\/doi.org\/10.1109\/NNSP.1999.788121<\/a><\/li>\n<li>Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., Blondel, M., Prettenhofer, P., Weiss, R., Dubourg, V., Vanderplas, J., Passos, A., Cournapeau, D., Brucher, M., Perrot, M., &amp; Duchesnay, E. (2011). Scikit-learn: Machine learning in Python. <em>Journal of Machine Learning Research<\/em>, 12, 2825\u20132830. <a rel=\"noopener noreferrer\" href=\"https:\/\/jmlr.org\/papers\/v12\/pedregosa11a.html\">https:\/\/jmlr.org\/papers\/v12\/pedregosa11a.html<\/a><\/li>\n<li>Tabachnick, B. G., &amp; Fidell, L. S. (2019). <em>Using multivariate statistics<\/em> (7th ed.). Pearson.<\/li>\n<li>Venables, W. N., &amp; Ripley, B. D. (2002). <em>Modern applied statistics with S<\/em> (4th ed.). Springer. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-0-387-21706-2\">https:\/\/doi.org\/10.1007\/978-0-387-21706-2<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique approfondi sur l&rsquo;analyse discriminante lin\u00e9aire (ADL) : fondements math\u00e9matiques, comparaison avec la r\u00e9gression logistique et applications.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2,5048],"tags":[6814,6834,6835,16,6829],"class_list":["post-4871","psych_stat","type-psych_stat","status-publish","hentry","category-psychometrie","category-statistiques-multivariees","tag-analyse-discriminante-lineaire","tag-classification-multiclasses","tag-methodes-statistiques","tag-psychometrie","tag-reduction-de-dimensionnalite"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/4871","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=4871"}],"version-history":[{"count":0,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/4871\/revisions"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=4871"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=4871"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=4871"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}