{"id":7192,"date":"2026-09-16T18:09:32","date_gmt":"2026-09-16T18:09:32","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/estimation-maximum-vraisemblance-distribution-uniforme\/"},"modified":"2026-09-16T18:09:35","modified_gmt":"2026-09-16T18:09:35","slug":"estimation-maximum-vraisemblance-distribution-uniforme","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=estimation-maximum-vraisemblance-distribution-uniforme","title":{"rendered":"Estimation par le maximum de vraisemblance (EMV) pour une distribution uniforme"},"content":{"rendered":"<p>L&rsquo;inf\u00e9rence statistique contemporaine repose en grande partie sur des principes d&rsquo;optimisation fonctionnelle formalis\u00e9s au d\u00e9but du vingti\u00e8me si\u00e8cle. Parmi ces paradigmes, l&rsquo;estimation par le maximum de vraisemblance s&rsquo;est impos\u00e9e comme le standard dominant pour d\u00e9river des estimateurs ponctuels dot\u00e9s de propri\u00e9t\u00e9s asymptotiques remarquables telles que la consistance, l&rsquo;efficacit\u00e9 asymptotique et la normalit\u00e9 de la distribution d&rsquo;\u00e9chantillonnage. Cependant, la validit\u00e9 universelle de ces propri\u00e9t\u00e9s repose de mani\u00e8re critique sur un ensemble de conditions de r\u00e9gularit\u00e9 math\u00e9matique, souvent regroup\u00e9es sous le vocable d&rsquo;hypoth\u00e8ses de Cram\u00e9r-Rao. Lorsque ces postulats structurels sont enfreints, l&rsquo;\u00e9difice analytique traditionnel s&rsquo;effondre, obligeant le statisticien \u00e0 abandonner le confort du calcul diff\u00e9rentiel standard au profit de l&rsquo;analyse convexe, de la topologie des ensembles ferm\u00e9s et de la th\u00e9orie fine des statistiques d&rsquo;ordre.<\/p>\n<p>Le cas paradigmatique illustrant cette rupture fondamentale est sans conteste l&rsquo;estimation des bornes d&rsquo;une <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Loi_uniforme_continue\">distribution uniforme continue<\/a> sur un intervalle compact. Alors que la loi uniforme est conceptuellement la plus \u00e9l\u00e9mentaire de toutes les distributions continues en th\u00e9orie des probabilit\u00e9s, son traitement sous l&rsquo;angle de la vraisemblance r\u00e9v\u00e8le une complexit\u00e9 math\u00e9matique et g\u00e9om\u00e9trique fascinante. Le support de la distribution d\u00e9pend explicitement des param\u00e8tres que l&rsquo;on cherche \u00e0 estimer, ce qui constitue une violation frontale de la premi\u00e8re condition de r\u00e9gularit\u00e9 statistique. Par cons\u00e9quent, la fonction de log-vraisemblance ne peut \u00eatre maximis\u00e9e par la simple annulation de son gradient ; la d\u00e9riv\u00e9e par rapport aux param\u00e8tres demeure strictement non nulle sur tout l&rsquo;int\u00e9rieur de l&rsquo;espace admissible, repoussant l&rsquo;optimum vers la fronti\u00e8re g\u00e9om\u00e9trique d\u00e9finie par les observations extr\u00eames de l&rsquo;\u00e9chantillon.<\/p>\n<p>Ce trait\u00e9 exhaustif explore l&rsquo;ensemble des ramifications th\u00e9oriques, g\u00e9om\u00e9triques, probabilistes et appliqu\u00e9es de l&rsquo;estimation par le maximum de vraisemblance pour la loi uniforme continue. De la d\u00e9rivation analytique rigoureuse des estimateurs aux fronti\u00e8res du support \u00e0 l&rsquo;\u00e9tude approfondie de leur super-efficacit\u00e9 asymptotique \u00e0 l&rsquo;ordre <em>n<\/em>, en passant par la correction de leur biais structurel via le th\u00e9or\u00e8me de Lehmann-Scheff\u00e9 et leurs applications cruciales en psychom\u00e9trie, en sciences cognitives et dans la r\u00e9solution du c\u00e9l\u00e8bre probl\u00e8me des chars allemands, cette monographie offre une r\u00e9f\u00e9rence exhaustive pour les chercheurs, statisticiens et mod\u00e9lisateurs confront\u00e9s \u00e0 des probl\u00e8mes d&rsquo;inf\u00e9rence non r\u00e9guli\u00e8re.<\/p>\n<h2>1. Introduction th\u00e9orique \u00e0 la distribution uniforme continue et ses fondements<\/h2>\n<h3>1.1 D\u00e9finition axiomatique et fonction de densit\u00e9 de probabilit\u00e9<\/h3>\n<p>La formalisation axiomatique d&rsquo;une variable al\u00e9atoire continue selon la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Th%C3%A9orie_des_probabilit%C3%A9s\">th\u00e9orie des probabilit\u00e9s de Kolmogorov<\/a> n\u00e9cessite la d\u00e9finition d&rsquo;un espace probabilis\u00e9 mesurable. Consid\u00e9rons un espace probabilis\u00e9 (&Omega;, F, P) et une variable al\u00e9atoire r\u00e9elle X d\u00e9finie sur cet espace. Nous disons que X suit une loi uniforme continue sur un intervalle compact de la droite r\u00e9elle, not\u00e9 [a, b] avec -&infin; &lt; a &lt; b &lt; +&infin;, ce que l&rsquo;on note conventionnellement X ~ U(a, b), si sa mesure de probabilit\u00e9 est proportionnelle \u00e0 la mesure de Lebesgue sur cet intervalle et nulle en dehors.<\/p>\n<p>Sur le plan analytique, la fonction de densit\u00e9 de probabilit\u00e9, not\u00e9e f(x ; a, b), caract\u00e9rise la distribution locale de probabilit\u00e9. Elle s&rsquo;exprime de mani\u00e8re rigoureuse \u00e0 l&rsquo;aide de la fonction indicatrice du support compact sous la forme :<\/p>\n<p>f(x ; a, b) = [1 \/ (b &#8211; a)] &times; <strong>1<\/strong><sub>[a, b]<\/sub>(x)<\/p>\n<p>o\u00f9 <strong>1<\/strong><sub>[a, b]<\/sub>(x) est l&rsquo;op\u00e9rateur indicatrice valant 1 si l&rsquo;argument x appartient \u00e0 l&rsquo;intervalle ferm\u00e9 [a, b], et 0 dans le cas contraire. L&rsquo;adjonction explicite de cette fonction indicatrice n&rsquo;est pas une simple coquetterie de notation math\u00e9matique, mais une n\u00e9cessit\u00e9 absolue pour \u00e9viter des erreurs fondamentales lors de l&rsquo;int\u00e9gration et de l&rsquo;\u00e9valuation de la fonction de vraisemblance sur l&rsquo;espace r\u00e9el complet. L&rsquo;int\u00e9grale de cette densit\u00e9 sur l&rsquo;ensemble de la droite r\u00e9elle est strictement \u00e9gale \u00e0 l&rsquo;unit\u00e9, satisfaisant ainsi l&rsquo;axiome fondamental de normalisation de toute mesure de probabilit\u00e9.<\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/uniform_pic-2.jpg\" alt=\"Uniform distribution example\" \/><figcaption>Uniform distribution example<\/figcaption><\/figure>\n<p>La fonction de r\u00e9partition cumulative associ\u00e9e, F(x ; a, b) = P(X &le; x), s&rsquo;obtient par l&rsquo;int\u00e9gration d\u00e9finie de la fonction de densit\u00e9 de probabilit\u00e9 depuis la borne inf\u00e9rieure de son support \u00e9tendu. Pour tout x &isin; <strong>R<\/strong>, elle se formule par morceaux : elle est identiquement nulle pour toute valeur strictement inf\u00e9rieure \u00e0 a, cro\u00eet de mani\u00e8re strictement lin\u00e9aire selon l&rsquo;expression (x &#8211; a) \/ (b &#8211; a) pour x compris entre a et b, et sature \u00e0 la valeur 1 pour toute valeur strictement sup\u00e9rieure \u00e0 b. Cette lin\u00e9arit\u00e9 absolue sur le support traduit la propri\u00e9t\u00e9 d&rsquo;homog\u00e9n\u00e9it\u00e9 spatiale compl\u00e8te : la probabilit\u00e9 qu&rsquo;une r\u00e9alisation de la variable X appartienne \u00e0 un sous-intervalle quelconque de [a, b] ne d\u00e9pend exclusivement que de la longueur de ce sous-intervalle et nullement de son positionnement spatial relatif au sein du support.<\/p>\n<h3>1.2 Moments th\u00e9oriques et propri\u00e9t\u00e9s fondamentales<\/h3>\n<p>L&rsquo;\u00e9valuation des moments th\u00e9oriques d&rsquo;une variable al\u00e9atoire X distribu\u00e9e selon la loi U(a, b) s&rsquo;effectue par int\u00e9gration directe contre la mesure de Lebesgue. L&rsquo;esp\u00e9rance math\u00e9matique, premier moment ordinaire repr\u00e9sentant le centre de gravit\u00e9 de la masse de probabilit\u00e9, s&rsquo;\u00e9tablit formellement comme suit :<\/p>\n<p>E(X) = &int;<sub>-&infin;<\/sub><sup>+&infin;<\/sup> x f(x ; a, b) dx = &int;<sub>a<\/sub><sup>b<\/sup> [x \/ (b &#8211; a)] dx = [x<sup>2<\/sup> \/ (2(b &#8211; a))]<sub>a<\/sub><sup>b<\/sup> = (b<sup>2<\/sup> &#8211; a<sup>2<\/sup>) \/ (2(b &#8211; a)) = (a + b) \/ 2<\/p>\n<p>Cette esp\u00e9rance co\u00efncide rigoureusement avec le milieu arithm\u00e9tique des bornes de l&rsquo;intervalle, refl\u00e9tant la parfaite sym\u00e9trie g\u00e9om\u00e9trique de la densit\u00e9 par rapport \u00e0 son axe central.<\/p>\n<p>Le calcul du second moment centr\u00e9, ou variance th\u00e9orique Var(X), permet de quantifier la dispersion intrins\u00e8que des r\u00e9alisations autour de leur barycentre th\u00e9orique. En mobilisant l&rsquo;identit\u00e9 de K\u00f6nig-Huygens, Var(X) = E(X<sup>2<\/sup>) &#8211; [E(X)]<sup>2<\/sup>, nous \u00e9valuons d&rsquo;abord le second moment ordinaire :<\/p>\n<p>E(X<sup>2<\/sup>) = &int;<sub>a<\/sub><sup>b<\/sup> [x<sup>2<\/sup> \/ (b &#8211; a)] dx = [x<sup>3<\/sup> \/ (3(b &#8211; a))]<sub>a<\/sub><sup>b<\/sup> = (b<sup>3<\/sup> &#8211; a<sup>3<\/sup>) \/ (3(b &#8211; a)) = (a<sup>2<\/sup> + ab + b<sup>2<\/sup>) \/ 3<\/p>\n<p>En soustrayant le carr\u00e9 de l&rsquo;esp\u00e9rance math\u00e9matique, nous obtenons par factorisation alg\u00e9brique :<\/p>\n<p>Var(X) = (a<sup>2<\/sup> + ab + b<sup>2<\/sup>) \/ 3 &#8211; (a + b)<sup>2<\/sup> \/ 4 = (4a<sup>2<\/sup> + 4ab + 4b<sup>2<\/sup> &#8211; 3a<sup>2<\/sup> &#8211; 6ab &#8211; 3b<sup>2<\/sup>) \/ 12 = (b &#8211; a)<sup>2<\/sup> \/ 12<\/p>\n<p>Ce r\u00e9sultat remarquable met en lumi\u00e8re que la dispersion d\u00e9pend quadratiquement de l&rsquo;amplitude totale du support, normalis\u00e9e par la constante universelle 12 propre aux formes g\u00e9om\u00e9triques rectangulaires continues.<\/p>\n<p>L&rsquo;analyse des moments d&rsquo;ordres sup\u00e9rieurs permet de parfaire la caract\u00e9risation structurelle de la distribution. En raison de la sym\u00e9trie centrale parfaite de la densit\u00e9 f(x) par rapport \u00e0 &mu; = (a + b) \/ 2, tous les moments centr\u00e9s impairs &mu;<sub>2k+1<\/sub> = E[(X &#8211; &mu;)<sup>2k+1<\/sup>] sont identiquement nuls pour tout entier naturel k. En particulier, le coefficient d&rsquo;asym\u00e9trie (skewness) est strictement \u00e9gal \u00e0 z\u00e9ro. Pour les moments pairs, le quatri\u00e8me moment centr\u00e9 vaut &mu;<sub>4<\/sub> = (b &#8211; a)<sup>4<\/sup> \/ 80, ce qui conduit \u00e0 un coefficient d&rsquo;aplatissement ou kurtosis standardis\u00e9 &beta;<sub>2<\/sub> = &mu;<sub>4<\/sub> \/ (Var(X))<sup>2<\/sup> = [ (b &#8211; a)<sup>4<\/sup> \/ 80 ] \/ [ (b &#8211; a)<sup>4<\/sup> \/ 144 ] = 144 \/ 80 = 9 \/ 5 = 1,8. L&rsquo;exc\u00e8s de kurtosis, mesur\u00e9 par &gamma;<sub>2<\/sub> = &beta;<sub>2<\/sub> &#8211; 3, s&rsquo;\u00e9tablit donc \u00e0 -1,2. Cette valeur strictement n\u00e9gative traduit une distribution fortement platykurtique, caract\u00e9ris\u00e9e par une absence totale de queues lourdes et des flancs tronqu\u00e9s de mani\u00e8re abrupte.<\/p>\n<p>Sur le plan philosophique et \u00e9pist\u00e9mologique, la distribution uniforme incarne la traduction formelle du <em>principe d&rsquo;indiff\u00e9rence<\/em> \u00e9nonc\u00e9 par Pierre-Simon de Laplace. En vertu de ce principe, en pr\u00e9sence d&rsquo;un ensemble continu d&rsquo;\u00e9v\u00e9nements \u00e9l\u00e9mentaires born\u00e9s sans aucune information pr\u00e9alable permettant d&rsquo;accorder une pr\u00e9f\u00e9rence heuristique \u00e0 un point plut\u00f4t qu&rsquo;\u00e0 un autre, l&rsquo;attribution d&rsquo;une densit\u00e9 de probabilit\u00e9 constante repr\u00e9sente l&rsquo;unique assignation rationnelle exempte de biais. Dans le cadre de la th\u00e9orie de l&rsquo;information de Claude Shannon, la loi uniforme continue maximise l&rsquo;entropie diff\u00e9rentielle parmi toutes les distributions continues restreintes \u00e0 un support compact donn\u00e9 [a, b], formalisant de facto l&rsquo;\u00e9tat d&rsquo;incertitude \u00e9pist\u00e9mique absolue sur cet intervalle.<\/p>\n<h3>1.3 Importance de la loi uniforme en mod\u00e9lisation math\u00e9matique et cognitive<\/h3>\n<p>Bien que d&rsquo;apparence simpliste, la distribution uniforme continue constitue un pilier de la mod\u00e9lisation math\u00e9matique appliqu\u00e9e, particuli\u00e8rement dans les contextes o\u00f9 des seuils physiques, physiologiques ou informationnels rigides contraignent la g\u00e9n\u00e9ration des donn\u00e9es. En <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Psychophysique\">psychophysique<\/a>, par exemple, la d\u00e9tection absolue de signaux sensoriels faibles (visuels, auditifs ou tactiles) est fr\u00e9quemment formul\u00e9e \u00e0 travers des mod\u00e8les de seuil continu o\u00f9 la sensibilit\u00e9 instantan\u00e9e d&rsquo;un r\u00e9cepteur oscille de mani\u00e8re \u00e9quiprobable entre deux fronti\u00e8res critiques d\u00e9pendant de l&rsquo;\u00e9tat r\u00e9fractaire des membranes synaptiques. La mod\u00e9lisation de ces seuils de d\u00e9tection requiert l&rsquo;estimation pr\u00e9cise de ces fronti\u00e8res d&rsquo;activation sous-jacentes \u00e0 partir d&rsquo;\u00e9chantillons comportementaux discrets.<\/p>\n<p>De plus, en th\u00e9orie de la mesure et en statistique computationnelle, la loi uniforme joue un r\u00f4le de distribution de r\u00e9f\u00e9rence canonique en vertu du th\u00e9or\u00e8me de la transformation int\u00e9grale de probabilit\u00e9. Ce r\u00e9sultat fondamental stipule que si une variable al\u00e9atoire X admet une fonction de r\u00e9partition continue F<sub>X<\/sub>, alors la variable transform\u00e9e U = F<sub>X<\/sub>(X) suit rigoureusement une loi uniforme standard sur l&rsquo;intervalle [0, 1]. Cette propri\u00e9t\u00e9 est le moteur algorithmique universel des g\u00e9n\u00e9rateurs de nombres pseudo-al\u00e9atoires et constitue le fondement math\u00e9matique des tests d&rsquo;ad\u00e9quation (goodness-of-fit), \u00e0 l&rsquo;image du test de Kolmogorov-Smirnov, o\u00f9 la distance supr\u00e9mum entre la fonction de r\u00e9partition empirique et la loi th\u00e9orique est directement \u00e9valu\u00e9e sur l&rsquo;espace d&rsquo;uniformit\u00e9 [0, 1].<\/p>\n<h2>2. Le paradigme de l&rsquo;estimation par le maximum de vraisemblance (EMV)<\/h2>\n<h3>2.1 Philosophie g\u00e9n\u00e9rale de l&rsquo;inf\u00e9rence par vraisemblance<\/h3>\n<p>L&rsquo;inf\u00e9rence statistique moderne a \u00e9t\u00e9 profond\u00e9ment r\u00e9volutionn\u00e9e par les travaux fondateurs de <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Ronald_Aylmer_Fisher\">Sir Ronald A. Fisher<\/a> au d\u00e9but des ann\u00e9es 1920, qui a syst\u00e9matis\u00e9 la m\u00e9thode du maximum de vraisemblance. La d\u00e9marche fisherienne repose sur une inversion conceptuelle \u00e9l\u00e9gante du conditionnement probabiliste conventionnel. Dans un cadre de probabilit\u00e9 directe, le param\u00e8tre populationnel &theta; appartenant \u00e0 un espace param\u00e9trique &Theta; est postul\u00e9 fixe et connu, et le probabiliste cherche \u00e0 d\u00e9terminer la loi de probabilit\u00e9 ou la fr\u00e9quence d&rsquo;apparition des vecteurs de donn\u00e9es d&rsquo;\u00e9chantillon observables X = (X<sub>1<\/sub>, &#8230;, X<sub>n<\/sub>).<\/p>\n<p>\u00c0 l&rsquo;inverse, dans le paradigme inf\u00e9rentiel inverse qui caract\u00e9rise la d\u00e9marche du statisticien confront\u00e9 \u00e0 des donn\u00e9es empiriques, l&rsquo;\u00e9chantillon de donn\u00e9es x = (x<sub>1<\/sub>, &#8230;, x<sub>n<\/sub>) a d\u00e9j\u00e0 \u00e9t\u00e9 recueilli et se trouve fix\u00e9 de mani\u00e8re immuable, tandis que le param\u00e8tre &theta; demeure une entit\u00e9 inconnue et inaccessible que l&rsquo;on souhaite identifier. Fisher a ainsi introduit la notion de fonction de vraisemblance, not\u00e9e L(&theta; ; x), d\u00e9finie analytiquement comme la densit\u00e9 de probabilit\u00e9 conjointe r\u00e9interpr\u00e9t\u00e9e en tant que fonction du seul param\u00e8tre &theta; sur &Theta;, conditionnellement au vecteur d&rsquo;observations x :<\/p>\n<p>L(&theta; ; x) = f<sub>X<\/sub>(x ; &theta;)<\/p>\n<p>Il est imp\u00e9ratif de souligner la fronti\u00e8re \u00e9pist\u00e9mologique stricte qui s\u00e9pare la probabilit\u00e9 de la vraisemblance. La fonction de vraisemblance L(&theta; ; x) n&rsquo;est absolument pas une densit\u00e9 de probabilit\u00e9 sur l&rsquo;espace param\u00e9trique &Theta; : son int\u00e9grale sur &Theta; n&rsquo;est en aucun cas contrainte d&rsquo;\u00e9galer 1 et peut diverger sans enfreindre la coh\u00e9rence du cadre d&rsquo;analyse. Le principe du maximum de vraisemblance postule que le meilleur estimateur ponctuel d&rsquo;un param\u00e8tre &theta; inconnu correspond \u00e0 la valeur au sein de &Theta; qui attribue le degr\u00e9 de plausibilit\u00e9 ou la probabilit\u00e9 a posteriori relative la plus \u00e9lev\u00e9e \u00e0 la constellation exacte des donn\u00e9es qui ont \u00e9t\u00e9 empiriquement observ\u00e9es.<\/p>\n<h3>2.2 Conditions classiques de r\u00e9gularit\u00e9 de Cram\u00e9r-Rao<\/h3>\n<p>L&rsquo;immense popularit\u00e9 de l&rsquo;estimateur du maximum de vraisemblance (EMV) dans la litt\u00e9rature th\u00e9orique et appliqu\u00e9e provient de ses propri\u00e9t\u00e9s optimales en grand \u00e9chantillon. N\u00e9anmoins, l&rsquo;obtention des c\u00e9l\u00e8bres th\u00e9or\u00e8mes de consistance forte, d&rsquo;efficacit\u00e9 asymptotique et de normalit\u00e9 asymptotique repose sur le respect scrupuleux d&rsquo;un ensemble de conditions restrictives, commun\u00e9ment d\u00e9sign\u00e9es sous l&rsquo;appellation de conditions de r\u00e9gularit\u00e9 de Cram\u00e9r-Rao ou conditions de r\u00e9gularit\u00e9 du mod\u00e8le statistique standard.<\/p>\n<p>Parmi ces conditions postulatoires formelles, les plus d\u00e9terminantes sont les suivantes :<\/p>\n<ul>\n<li><strong>Ind\u00e9pendance du support :<\/strong> L&rsquo;ensemble de d\u00e9finition ou support de la fonction de densit\u00e9, d\u00e9fini par S<sub>&theta;<\/sub> = {x &isin; <strong>R<\/strong> : f(x ; &theta;) &gt; 0}, doit \u00eatre strictement ind\u00e9pendant du param\u00e8tre inconnu &theta;.<\/li>\n<li><strong>Ouverture de l&rsquo;espace param\u00e9trique :<\/strong> Le param\u00e8tre d&rsquo;int\u00e9r\u00eat &theta; doit appartenir \u00e0 un ensemble ouvert non vide de <strong>R<\/strong><sup>p<\/sup>, garantissant que les extrema locaux puissent \u00eatre caract\u00e9ris\u00e9s par l&rsquo;annulation du gradient vectoriel sans rencontrer imm\u00e9diatement des contraintes d&rsquo;in\u00e9galit\u00e9 actives.<\/li>\n<li><strong>Diff\u00e9rentiabilit\u00e9 et int\u00e9grabilit\u00e9 :<\/strong> La fonction de log-vraisemblance doit \u00eatre au moins deux fois contin\u00fbment diff\u00e9rentiable par rapport \u00e0 &theta;, et les op\u00e9rations de d\u00e9rivation par rapport au param\u00e8tre et d&rsquo;int\u00e9gration par rapport \u00e0 la variable d&rsquo;espace x doivent pouvoir \u00eatre permut\u00e9es librement, conform\u00e9ment au th\u00e9or\u00e8me de convergence domin\u00e9e de Lebesgue.<\/li>\n<li><strong>Information de Fisher non d\u00e9g\u00e9n\u00e9r\u00e9e :<\/strong> La matrice d&rsquo;information de Fisher I(&theta;) = E[( &nabla; ln f(X ; &theta;) ) ( &nabla; ln f(X ; &theta;) )<sup>T<\/sup>] doit \u00eatre strictement d\u00e9finie positive et finie sur l&rsquo;espace des param\u00e8tres.<\/li>\n<\/ul>\n<p>Lorsque ces crit\u00e8res sont satisfaits, la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Borne_de_Cram%C3%A9r-Rao\">borne de Cram\u00e9r-Rao<\/a> impose une limite inf\u00e9rieure infranchissable \u00e0 la variance de tout estimateur sans biais, et l&rsquo;EMV atteint asymptotiquement cette borne avec une vitesse de convergence d&rsquo;ordre standard proportionnelle \u00e0 &radic;n, sa loi limite se conformant rigoureusement \u00e0 une distribution normale multivari\u00e9e centr\u00e9e sur la valeur vraie du param\u00e8tre.<\/p>\n<h3>2.3 Rupture des conditions de r\u00e9gularit\u00e9 dans la distribution uniforme<\/h3>\n<p>L&rsquo;application du mod\u00e8le de la distribution uniforme continue U(a, b) provoque une d\u00e9faillance catastrophique et imm\u00e9diate de l&rsquo;ensemble de ces hypoth\u00e8ses de travail traditionnelles. La violation la plus flagrante concerne l&rsquo;hypoth\u00e8se d&rsquo;ind\u00e9pendance du support. En effet, pour la loi U(a, b), le support effectif s&rsquo;\u00e9crit de mani\u00e8re explicite sous la forme de l&rsquo;intervalle ferm\u00e9 S<sub>(a, b)<\/sub> = [a, b]. Ce support est intrins\u00e8quement et indissociablement li\u00e9 aux param\u00e8tres inconnus a et b que l&rsquo;on cherche pr\u00e9cis\u00e9ment \u00e0 d\u00e9duire des donn\u00e9es.<\/p>\n<p>Cette d\u00e9pendance param\u00e9trique du support entra\u00eene des cons\u00e9quences math\u00e9matiques en cascade d&rsquo;une extr\u00eame gravit\u00e9 :<\/p>\n<ul>\n<li><strong>Non-diff\u00e9rentiabilit\u00e9 fonctionnelle :<\/strong> La fonction de densit\u00e9 pr\u00e9sente des discontinuit\u00e9s de premi\u00e8re esp\u00e8ce (sauts finis) exactement au niveau des bornes a et b. D\u00e8s lors, la fonction de vraisemblance conjointe n&rsquo;est pas diff\u00e9rentiable sur les fronti\u00e8res d\u00e9finies par les valeurs des donn\u00e9es \u00e9chantillonnales observ\u00e9es.<\/li>\n<li><strong>Inversion interdite de l&rsquo;int\u00e9gration et de la d\u00e9rivation :<\/strong> La r\u00e8gle de d\u00e9rivation sous le signe int\u00e9gral de Leibniz-Reynolds fait intervenir la d\u00e9pendance des bornes d&rsquo;int\u00e9gration envers les param\u00e8tres, g\u00e9n\u00e9rant des termes de fronti\u00e8re non nuls qui emp\u00eachent l&rsquo;\u00e9galisation classique entre l&rsquo;esp\u00e9rance du score de Fisher et z\u00e9ro.<\/li>\n<li><strong>\u00c9vanouissement de l&rsquo;information de Fisher standard :<\/strong> L&rsquo;information de Fisher classique devient formellement non d\u00e9finie ou infinie au sens distributionnel au niveau des singularit\u00e9s, rendant le calcul standard de la borne de Cram\u00e9r-Rao totalement inapplicable.<\/li>\n<\/ul>\n<p>Par cons\u00e9quent, les th\u00e9or\u00e8mes classiques assurant la convergence \u00e0 vitesse &radic;n vers une loi normale ne s&rsquo;appliquent pas. L&rsquo;inf\u00e9rence par maximum de vraisemblance pour la distribution uniforme ne rel\u00e8ve pas d&rsquo;une simple r\u00e9solution d&rsquo;\u00e9quations diff\u00e9rentielles de score, mais constitue un probl\u00e8me d&rsquo;optimisation non lin\u00e9aire sous contraintes rigides sur un domaine ferm\u00e9, menant \u00e0 des distributions limites radicalement asym\u00e9triques et non gaussiennes.<\/p>\n<h2>3. Formulation math\u00e9matique de la fonction de vraisemblance pour la loi uniforme<\/h2>\n<h3>3.1 Hypoth\u00e8se d&rsquo;\u00e9chantillonnage i.i.d. et produit tensoriel des densit\u00e9s<\/h3>\n<p>Soit un \u00e9chantillon al\u00e9atoire constitu\u00e9 de n variables al\u00e9atoires r\u00e9elles ind\u00e9pendantes et identiquement distribu\u00e9es (i.i.d.), not\u00e9 X = (X<sub>1<\/sub>, X<sub>2<\/sub>, &#8230;, X<sub>n<\/sub>), pr\u00e9lev\u00e9 \u00e0 partir d&rsquo;une population parente gouvern\u00e9e par une loi uniforme continue U(a, b), o\u00f9 le couple de param\u00e8tres &theta; = (a, b) appartient \u00e0 l&rsquo;espace admissible &Theta; = {(a, b) &isin; <strong>R<\/strong><sup>2<\/sup> : a &lt; b}. L&rsquo;ind\u00e9pendance mutuelle des composantes implique que la mesure de probabilit\u00e9 conjointe du vecteur d&rsquo;\u00e9chantillonnage s&rsquo;identifie au produit tensoriel des mesures marginales individuelles.<\/p>\n<p>Sous cette hypoth\u00e8se probabiliste centrale, la fonction de densit\u00e9 de probabilit\u00e9 conjointe de l&rsquo;\u00e9chantillon complet, \u00e9valu\u00e9e au point vectoriel r\u00e9alis\u00e9 x = (x<sub>1<\/sub>, x<sub>2<\/sub>, &#8230;, x<sub>n<\/sub>) &isin; <strong>R<\/strong><sup>n<\/sup>, s&rsquo;\u00e9crit formellement comme le produit des n densit\u00e9s de probabilit\u00e9 univari\u00e9es :<\/p>\n<p>f<sub>X<\/sub>(x<sub>1<\/sub>, &#8230;, x<sub>n<\/sub> ; a, b) = &prod;<sub>i=1<\/sub><sup>n<\/sup> f(x<sub>i<\/sub> ; a, b)<\/p>\n<p>L&rsquo;introduction rigoureuse des fonctions indicatrices individuelles \u00e0 l&rsquo;int\u00e9rieur de ce produit inf\u00e9rentiel est la cl\u00e9 de vo\u00fbte de toute l&rsquo;analyse math\u00e9matique ult\u00e9rieure. En omettant les indicatrices, un analyste s&rsquo;exposerait \u00e0 formuler une fonction de vraisemblance tronqu\u00e9e et erron\u00e9e, incapable de capter les contraintes g\u00e9om\u00e9triques impos\u00e9es par le positionnement spatial des donn\u00e9es.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_.png\" alt=\"likelihood function for uniform distribution\" \/><figcaption>likelihood function for uniform distribution<\/figcaption><\/figure>\n<h3>3.2 Expression analytique de la fonction de vraisemblance conjointe<\/h3>\n<p>En substituant la d\u00e9finition compl\u00e8te de la densit\u00e9 univari\u00e9e dans l&rsquo;op\u00e9rateur de produit fini, nous obtenons l&rsquo;expression litt\u00e9rale de la fonction de vraisemblance conjointe L(a, b ; x) :<\/p>\n<p>L(a, b ; x) = &prod;<sub>i=1<\/sub><sup>n<\/sup> { [1 \/ (b &#8211; a)] &times; <strong>1<\/strong><sub>[a, b]<\/sub>(x<sub>i<\/sub>) } = [ 1 \/ (b &#8211; a) ]<sup>n<\/sup> &times; &prod;<sub>i=1<\/sub><sup>n<\/sup> <strong>1<\/strong><sub>[a, b]<\/sub>(x<sub>i<\/sub>)<\/p>\n<p>Cette \u00e9criture met en exergue deux composantes structurellement distinctes qui interagissent de mani\u00e8re antagoniste dans le processus d&rsquo;estimation :<\/p>\n<ul>\n<li>Un facteur d&rsquo;\u00e9chelle continu et strictement positif, repr\u00e9sent\u00e9 par la puissance inverse de l&rsquo;amplitude du support : (b &#8211; a)<sup>-n<\/sup>. Ce terme d\u00e9pend contin\u00fbment de la diff\u00e9rence entre les deux param\u00e8tres inconnus.<\/li>\n<li>Un facteur d&rsquo;admissibilit\u00e9 bool\u00e9en, formalis\u00e9 par le produit bool\u00e9en infini des fonctions indicatrices : &prod;<sub>i=1<\/sub><sup>n<\/sup> <strong>1<\/strong><sub>[a, b]<\/sub>(x<sub>i<\/sub>). Ce terme agit comme un filtre de validit\u00e9 spatiale binaire, valant 1 si l&rsquo;int\u00e9gralit\u00e9 des r\u00e9alisations \u00e9chantillonnales s&rsquo;inscrit au sein de l&rsquo;intervalle [a, b], et chutant instantan\u00e9ment \u00e0 0 d\u00e8s qu&rsquo;une unique valeur x<sub>i<\/sub> enfreint cette appartenance.<\/li>\n<\/ul>\n<p>La fonction de vraisemblance conjointe s&rsquo;exprime donc sous la forme conditionnelle suivante :<\/p>\n<p>L(a, b ; x) = (b &#8211; a)<sup>-n<\/sup> si &forall; i &isin; {1, &#8230;, n}, a &le; x<sub>i<\/sub> &le; b ; et 0 sinon.<\/p>\n<h3>3.3 R\u00f4le des statistiques d&rsquo;ordre dans la fonction indicatrice conjointe<\/h3>\n<p>Pour rendre cette fonction de vraisemblance traitable sur le plan analytique, il est n\u00e9cessaire de simplifier la conjonction des n doubles in\u00e9galit\u00e9s a &le; x<sub>i<\/sub> &le; b pour tout i &isin; {1, &#8230;, n}. C&rsquo;est ici qu&rsquo;intervient la th\u00e9orie des <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Statistique_d%27ordre\">statistiques d&rsquo;ordre<\/a>. Consid\u00e9rons l&rsquo;\u00e9chantillon ordonn\u00e9 par permutation croissante des valeurs observ\u00e9es, not\u00e9 usuellement X<sub>(1)<\/sub> &le; X<sub>(2)<\/sub> &le; &#8230; &le; X<sub>(n)<\/sub>, o\u00f9 X<sub>(1)<\/sub> = min{X<sub>1<\/sub>, &#8230;, X<sub>n<\/sub>} d\u00e9signe le minimum empirique de l&rsquo;\u00e9chantillon, et X<sub>(n)<\/sub> = max{X<sub>1<\/sub>, &#8230;, X<sub>n<\/sub>} repr\u00e9sente le maximum empirique de l&rsquo;\u00e9chantillon.<\/p>\n<p>L&rsquo;\u00e9quivalence logique fondamentale liant les observations individuelles aux statistiques d&rsquo;ordre extr\u00eames s&rsquo;\u00e9tablit de fa\u00e7on imm\u00e9diate :<\/p>\n<p>(&forall; i &isin; {1, &#8230;, n}, a &le; x<sub>i<\/sub> &le; b) &hArr; (a &le; min<sub>1&le;i&le;n<\/sub> x<sub>i<\/sub>) et (max<sub>1&le;i&le;n<\/sub> x<sub>i<\/sub> &le; b) &hArr; (a &le; x<sub>(1)<\/sub>) et (x<sub>(n)<\/sub> &le; b)<\/p>\n<p>Cette r\u00e9duction alg\u00e9brique permet de fusionner le produit de n indicatrices en une indicatrice conjointe ne faisant intervenir que les deux statistiques extr\u00eames :<\/p>\n<p>&prod;<sub>i=1<\/sub><sup>n<\/sup> <strong>1<\/strong><sub>[a, b]<\/sub>(x<sub>i<\/sub>) = <strong>1<\/strong><sub>{a &le; x<sub>(1)<\/sub>}<\/sub> &times; <strong>1<\/strong><sub>{x<sub>(n)<\/sub> &le; b}<\/sub> = <strong>1<\/strong><sub>{a &le; x<sub>(1)<\/sub> &le; x<sub>(n)<\/sub> &le; b}<\/sub><\/p>\n<p>En vertu du th\u00e9or\u00e8me de factorisation de Fisher-Neyman, la fonction de vraisemblance conjointe se factorise sous la forme L(a, b ; x) = g(x<sub>(1)<\/sub>, x<sub>(n)<\/sub> ; a, b) &times; h(x), o\u00f9 la fonction r\u00e9siduelle h(x) est identiquement \u00e9gale \u00e0 la constante 1. Cette d\u00e9composition d\u00e9montre rigoureusement que le vecteur bidimensionnel constitu\u00e9 du couple des statistiques d&rsquo;ordre extr\u00eames T(X) = (X<sub>(1)<\/sub>, X<sub>(n)<\/sub>) constitue une <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Statistique_exhaustive\">statistique exhaustive<\/a> (ou suffisante) conjointe pour le couple de param\u00e8tres inconnus (a, b). Toutes les informations informatives contenues dans l&rsquo;\u00e9chantillon de taille n concernant la localisation des bornes inconnues se trouvent enti\u00e8rement condens\u00e9es dans ses deux valeurs extr\u00eames.<\/p>\n<h2>4. Passage \u00e0 la fonction de log-vraisemblance et analyse des contraintes de support<\/h2>\n<h3>4.1 Transformation logarithmique et son domaine de d\u00e9finition<\/h3>\n<p>Dans la pratique statistique conventionnelle, la fonction de log-vraisemblance &amp;ell;(a, b ; x) = ln L(a, b ; x) est syst\u00e9matiquement pr\u00e9f\u00e9r\u00e9e \u00e0 la vraisemblance brute en raison de la lin\u00e9arisation avantageuse qu&rsquo;apporte la transformation logarithmique n\u00e9p\u00e9rienne sur les op\u00e9rateurs de produits. Puisque le logarithme est un hom\u00e9omorphisme strictement croissant de l&rsquo;ensemble des r\u00e9els strictement positifs <strong>R<\/strong><sub>+<\/sub><sup>*<\/sup> sur la droite r\u00e9elle <strong>R<\/strong>, le lieu des arguments maximisant la vraisemblance co\u00efncide de mani\u00e8re absolue avec le lieu des arguments maximisant la log-vraisemblance.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_-1.png\" alt=\"Log-likelihood function of uniform distribution\" \/><figcaption>Log-likelihood function of uniform distribution<\/figcaption><\/figure>\n<p>Toutefois, une vigilance math\u00e9matique rigoureuse s&rsquo;impose lors de l&rsquo;application de cette transformation \u00e0 la loi uniforme. La pr\u00e9sence du produit d&rsquo;indicatrices modifie radicalement le domaine de d\u00e9finition r\u00e9el de la log-vraisemblance. Formellement, nous obtenons :<\/p>\n<p>&amp;ell;(a, b ; x) = ln [ (b &#8211; a)<sup>-n<\/sup> &times; <strong>1<\/strong><sub>{a &le; x<sub>(1)<\/sub> &le; x<sub>(n)<\/sub> &le; b}<\/sub> ]<\/p>\n<p>Par convention analytique d&rsquo;arithm\u00e9tique \u00e9tendue sur la droite r\u00e9elle achev\u00e9e [-&infin;, +&infin;], lorsque l&rsquo;argument de la fonction logarithme s&rsquo;annule, celle-ci plonge vers l&rsquo;infini n\u00e9gatif : ln(0) = -&infin;. La formulation rigoureuse de la log-vraisemblance s&rsquo;\u00e9nonce donc comme une fonction indic\u00e9e sur l&rsquo;ensemble de l&rsquo;espace <strong>R<\/strong><sup>2<\/sup> :<\/p>\n<p>&amp;ell;(a, b ; x) = -n ln(b &#8211; a) si a &le; x<sub>(1)<\/sub> et b &ge; x<sub>(n)<\/sub> ; et -&infin; d\u00e8s lors que (a &gt; x<sub>(1)<\/sub> ou b &lt; x<sub>(n)<\/sub>).<\/p>\n<p>La r\u00e9gion o\u00f9 la log-vraisemblance prend des valeurs finies est appel\u00e9e le <em>domaine de vraisemblance admissible<\/em>. En dehors de cette zone de coh\u00e9rence observationnelle, la log-vraisemblance est infiniment r\u00e9pulsive, interdisant de facto toute configuration param\u00e9trique qui exclurait ne serait-ce qu&rsquo;une seule observation r\u00e9alis\u00e9e.<\/p>\n<h3>4.2 Topologie de l&rsquo;espace des param\u00e8tres admissibles<\/h3>\n<p>L&rsquo;exploration de la structure g\u00e9om\u00e9trique du domaine de log-vraisemblance admissible &Omega;<sub>x<\/sub> dans le plan euclidien <strong>R<\/strong><sup>2<\/sup>, param\u00e9tr\u00e9 par les coordonn\u00e9es cart\u00e9siennes (a, b), r\u00e9v\u00e8le la g\u00e9om\u00e9trie sp\u00e9cifique de ce probl\u00e8me d&rsquo;estimation. Ce domaine &Omega;<sub>x<\/sub> est d\u00e9fini par l&rsquo;intersection spatiale de trois demi-plans ouverts et ferm\u00e9s :<\/p>\n<p>&Omega;<sub>x<\/sub> = { (a, b) &isin; <strong>R<\/strong><sup>2<\/sup> : a &lt; b } &cap; { (a, b) &isin; <strong>R<\/strong><sup>2<\/sup> : a &le; x<sub>(1)<\/sub> } &cap; { (a, b) &isin; <strong>R<\/strong><sup>2<\/sup> : b &ge; x<sub>(n)<\/sub> }<\/p>\n<p>Puisque par construction des donn\u00e9es empiriques nous avons toujours x<sub>(1)<\/sub> &le; x<sub>(n)<\/sub>, et en supposant au moins deux observations distinctes de sorte que x<sub>(1)<\/sub> &lt; x<sub>(n)<\/sub>, la condition structurelle a &lt; b est automatiquement satisfaite pour tout couple respectant a &le; x<sub>(1)<\/sub> et b &ge; x<sub>(n)<\/sub>.<\/p>\n<p>Le domaine admissible &Omega;<sub>x<\/sub> se visualise graphiquement comme un quart de plan g\u00e9om\u00e9trique infini, d\u00e9limit\u00e9 inf\u00e9rieurement par la droite horizontale d&rsquo;\u00e9quation b = x<sub>(n)<\/sub> et d\u00e9limit\u00e9 \u00e0 droite par la droite verticale d&rsquo;\u00e9quation a = x<sub>(1)<\/sub>. Sur le plan topologique au sens de la m\u00e9trique usuelle de <strong>R<\/strong><sup>2<\/sup>, l&rsquo;ensemble &Omega;<sub>x<\/sub> pr\u00e9sente les caract\u00e9ristiques fondamentales suivantes :<\/p>\n<ul>\n<li>Il s&rsquo;agit d&rsquo;un ensemble non-compact, car il est manifestement non born\u00e9 vers la gauche (a peut tendre vers -&infin;) et vers le haut (b peut tendre vers +&infin;).<\/li>\n<li>Il s&rsquo;agit d&rsquo;un ensemble ferm\u00e9 de <strong>R<\/strong><sup>2<\/sup>, dans la mesure o\u00f9 il contient la totalit\u00e9 de ses points d&rsquo;adh\u00e9rence et ses deux fronti\u00e8res rectilignes support\u00e9es par a = x<sub>(1)<\/sub> et b = x<sub>(n)<\/sub>.<\/li>\n<li>Il s&rsquo;agit d&rsquo;un ensemble convexe r\u00e9gulier : pour tous couples param\u00e9triques (&theta;<sub>1<\/sub>, &theta;<sub>2<\/sub>) appartenant \u00e0 &Omega;<sub>x<\/sub>, le segment lin\u00e9aire qui les relie est int\u00e9gralement contenu dans &Omega;<sub>x<\/sub>.<\/li>\n<\/ul>\n<p>L&rsquo;intersection des deux fronti\u00e8res g\u00e9om\u00e9triques constitue un point singulier extr\u00eame, le coin orthogonal du domaine admissible, dont les coordonn\u00e9es exactes sont (&acirc;, b&#770;)<sub>fronti\u00e8re<\/sub> = (x<sub>(1)<\/sub>, x<sub>(n)<\/sub>).<\/p>\n<h3>4.3 Analyse du comportement asymptotique aux fronti\u00e8res<\/h3>\n<p>Sur l&rsquo;int\u00e9rieur ouvert du domaine admissible &Omega;<sub>x<\/sub>, la fonction de log-vraisemblance &amp;ell;(a, b) = -n ln(b &#8211; a) est une fonction infiniment d\u00e9rivable de classe C<sup>&infin;<\/sup>. Examinons avec minutie son profil asymptotique lorsque le vecteur param\u00e9trique se d\u00e9place vers les fronti\u00e8res ou vers l&rsquo;infini :<\/p>\n<ul>\n<li><strong>\u00c9loignement vers l&rsquo;infini :<\/strong> Si la borne a tend vers -&infin; pour un param\u00e8tre b fix\u00e9, ou si b tend vers +&infin; pour un a fix\u00e9, l&rsquo;\u00e9cart spatial (b &#8211; a) diverge positivement vers +&infin;. En cons\u00e9quence logique, ln(b &#8211; a) &rarr; +&infin;, ce qui implique que &amp;ell;(a, b) = -n ln(b &#8211; a) tend in\u00e9luctablement vers -&infin;. La fonction de log-vraisemblance d\u00e9cro\u00eet donc contin\u00fbment au fur et \u00e0 mesure que l&rsquo;on \u00e9largit le support au-del\u00e0 du p\u00e9rim\u00e8tre des donn\u00e9es.<\/li>\n<li><strong>Rapprochement de l&rsquo;intervalle minimal admissible :<\/strong> Pour maximiser l&rsquo;expression &amp;ell;(a, b) = -n ln(b &#8211; a), il est strictement \u00e9quivalent de chercher \u00e0 minimiser l&rsquo;argument du logarithme, c&rsquo;est-\u00e0-dire l&rsquo;amplitude d&rsquo;intervalle &Delta;(a, b) = (b &#8211; a).<\/li>\n<\/ul>\n<p>Si aucune contrainte d&rsquo;admissibilit\u00e9 n&rsquo;\u00e9tait impos\u00e9e par les donn\u00e9es, la quantit\u00e9 (b &#8211; a) pourrait \u00eatre rendue arbitrairement proche de z\u00e9ro en contractant a et b l&rsquo;un vers l&rsquo;autre, ce qui conduirait ln(b &#8211; a) vers -&infin; et projetterait la vraisemblance L(a, b) vers +&infin; de fa\u00e7on explosive. Toutefois, une telle configuration est cat\u00e9goriquement proscrite par la pr\u00e9sence des contraintes indicatrices : a ne peut jamais exc\u00e9der x<sub>(1)<\/sub> et b ne peut jamais \u00eatre inf\u00e9rieur \u00e0 x<sub>(n)<\/sub>. D\u00e8s lors, la diff\u00e9rence minimale r\u00e9alisable sur l&rsquo;ensemble ferm\u00e9 admissible est strictement born\u00e9e inf\u00e9rieurement par l&rsquo;\u00e9tendue empirique des observations :<\/p>\n<p>min<sub>(a, b) &isin; &Omega;<sub>x<\/sub><\/sub> (b &#8211; a) = (x<sub>(n)<\/sub> &#8211; x<sub>(1)<\/sub>) &gt; 0<\/p>\n<p>Puisque la log-vraisemblance est continue sur le domaine ferm\u00e9 &Omega;<sub>x<\/sub>, strictement d\u00e9croissante par rapport \u00e0 la quantit\u00e9 (b &#8211; a), et qu&rsquo;elle plonge vers -&infin; \u00e0 l&rsquo;infini g\u00e9om\u00e9trique, le <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Th%C3%A9or%C3%A8me_des_bornes\">th\u00e9or\u00e8me de Weierstrass g\u00e9n\u00e9ralis\u00e9 aux ensembles ferm\u00e9s<\/a> garantit de fa\u00e7on p\u00e9remptoire l&rsquo;existence et l&rsquo;unicit\u00e9 absolue d&rsquo;un supremum global sur &Omega;<sub>x<\/sub>.<\/p>\n<h2>5. L&rsquo;\u00e9chec de l&rsquo;annulation du gradient et la nature fronti\u00e8re du probl\u00e8me<\/h2>\n<h3>5.1 Calcul des d\u00e9riv\u00e9es partielles formelles<\/h3>\n<p>Afin de d\u00e9montrer l&rsquo;\u00e9chec structurel des approches diff\u00e9rentielles classiques, calculons formellement le vecteur gradient de la fonction de log-vraisemblance sur l&rsquo;int\u00e9rieur ouvert du domaine admissible, not\u00e9 int(&Omega;<sub>x<\/sub>) = {(a, b) &isin; <strong>R<\/strong><sup>2<\/sup> : a &lt; x<sub>(1)<\/sub> et b &gt; x<sub>(n)<\/sub>}. En d\u00e9rivant partiellement &amp;ell;(a, b) = -n ln(b &#8211; a) par rapport au premier param\u00e8tre a, nous obtenons :<\/p>\n<p>&part;&amp;ell;(a, b) \/ &part;a = -n &times; [ &part;(b &#8211; a) \/ &part;a ] \/ (b &#8211; a) = -n &times; (-1) \/ (b &#8211; a) = + n \/ (b &#8211; a)<\/p>\n<p>De m\u00eame, en d\u00e9rivant partiellement par rapport au second param\u00e8tre b, nous \u00e9tablissons :<\/p>\n<p>&part;&amp;ell;(a, b) \/ &part;b = -n &times; [ &part;(b &#8211; a) \/ &part;b ] \/ (b &#8211; a) = -n &times; (+1) \/ (b &#8211; a) = &#8211; n \/ (b &#8211; a)<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_-2.png\" alt=\"Derivative of log-likelihood function for uniform distribution\" \/><figcaption>Derivative of log-likelihood function for uniform distribution<\/figcaption><\/figure>\n<p>Puisque par d\u00e9finition fondamentale d&rsquo;un support uniforme nous avons a &lt; b, la diff\u00e9rence scalaire (b &#8211; a) est strictement positive pour tout couple param\u00e9trique au sein de l&rsquo;int\u00e9rieur admissible. \u00c9tant donn\u00e9 que la taille d&rsquo;\u00e9chantillon n est un entier naturel strictement positif (n &ge; 1), il d\u00e9coule imm\u00e9diatement les in\u00e9galit\u00e9s strictes suivantes, v\u00e9rifi\u00e9es en tout point de l&rsquo;int\u00e9rieur du domaine :<\/p>\n<p>&part;&amp;ell;(a, b) \/ &part;a &gt; 0 &emsp; et &emsp; &part;&amp;ell;(a, b) \/ &part;b &lt; 0<\/p>\n<p>Le gradient vectoriel ne d\u00e9pend donc des param\u00e8tres que par l&rsquo;interm\u00e9diaire de leur diff\u00e9rence et conserve un signe immuable sur l&rsquo;ensemble ouvert de validit\u00e9.<\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_-3.png\" alt=\"Partial derivative of log-likelihood function for uniform distribution\" \/><figcaption>Partial derivative of log-likelihood function for uniform distribution<\/figcaption><\/figure>\n<h3>5.2 Impossibilit\u00e9 de satisfaire les conditions du premier ordre<\/h3>\n<p>Les conditions d&rsquo;optimalit\u00e9 du premier ordre d&rsquo;Euler-Lagrange pour la maximisation sans contrainte d&rsquo;une fonction contin\u00fbment diff\u00e9rentiable stipulent que tout extremum local situ\u00e9 \u00e0 l&rsquo;int\u00e9rieur du domaine doit imp\u00e9rativement constituer un point critique stationnaire, solution du syst\u00e8me vectoriel d&rsquo;\u00e9quations de score :<\/p>\n<p>&nabla;&amp;ell;(a, b) = ( &part;&amp;ell;\/&part;a, &part;&amp;ell;\/&part;b )<sup>T<\/sup> = ( 0, 0 )<sup>T<\/sup><\/p>\n<p>En substituant nos d\u00e9riv\u00e9es partielles formelles dans cette \u00e9quation vectorielle, le syst\u00e8me d&rsquo;optimisation classique devient :<\/p>\n<p>{ n \/ (b &#8211; a) = 0 ; -n \/ (b &#8211; a) = 0 }<\/p>\n<p>Puisque la taille d&rsquo;\u00e9chantillon n est finie et que l&rsquo;\u00e9cart (b &#8211; a) ne peut cro\u00eetre vers l&rsquo;infini sans pr\u00e9cipiter la log-vraisemblance vers -&infin;, ce syst\u00e8me n&rsquo;admet rigoureusement aucune solution sur l&rsquo;ensemble des r\u00e9els finis. Il y a une absence totale de points critiques int\u00e9rieurs au domaine.<\/p>\n<p>Cet \u00e9chec math\u00e9matique d\u00e9montre l&rsquo;inad\u00e9quation radicale des algorithmes d&rsquo;optimisation num\u00e9rique standard de type gradient descendant, m\u00e9thode de Newton-Raphson, ou algorithme de quasi-Newton (BFGS) lorsqu&rsquo;ils sont appliqu\u00e9s sans reformulation sous contraintes. Les solveurs cherchant \u00e0 \u00e9galiser le gradient \u00e0 z\u00e9ro divergeraient vers l&rsquo;infini ou \u00e9choueraient \u00e0 localiser l&rsquo;optimum. La nature math\u00e9matique du probl\u00e8me impose de recourir formellement au cadre de l&rsquo;optimisation sous contraintes d&rsquo;in\u00e9galit\u00e9s, formalis\u00e9 par les <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Conditions_de_Karush-Kuhn-Tucker\">conditions d&rsquo;optimalit\u00e9 de Karush-Kuhn-Tucker (KKT)<\/a> pour des fronti\u00e8res non lisses.<\/p>\n<h3>5.3 Comportement monotone et recherche de la solution aux bornes<\/h3>\n<p>La positivit\u00e9 stricte permanente de la d\u00e9riv\u00e9e partielle par rapport \u00e0 a (&part;&amp;ell;\/&part;a &gt; 0) indique que la log-vraisemblance est une fonction strictement croissante et monotone par rapport au param\u00e8tre a, pour toute valeur fix\u00e9e de b. Pour augmenter la vraisemblance, il convient donc d&rsquo;accro\u00eetre a de mani\u00e8re maximale. Or, le param\u00e8tre a est soumis \u00e0 la contrainte unilat\u00e9rale rigide a &le; x<sub>(1)<\/sub>. Toute augmentation de a au-del\u00e0 du seuil critique x<sub>(1)<\/sub> fait chuter instantan\u00e9ment la vraisemblance \u00e0 z\u00e9ro (et la log-vraisemblance \u00e0 -&infin;) en raison de l&rsquo;activation de la fonction indicatrice. La valeur admissible qui maximise L par rapport \u00e0 a est donc la borne sup\u00e9rieure de son intervalle d&rsquo;admissibilit\u00e9 (-&infin;, x<sub>(1)<\/sub>], c&rsquo;est-\u00e0-dire la fronti\u00e8re elle-m\u00eame : a* = x<sub>(1)<\/sub>.<\/p>\n<p>Sym\u00e9triquement, la n\u00e9gativit\u00e9 stricte permanente de la d\u00e9riv\u00e9e partielle par rapport \u00e0 b (&part;&amp;ell;\/&part;b &lt; 0) d\u00e9montre que la log-vraisemblance est une fonction strictement d\u00e9croissante et monotone par rapport au param\u00e8tre b, pour toute valeur fix\u00e9e de a. Pour accro\u00eetre la vraisemblance, il convient donc de diminuer b le plus possible. Cependant, le param\u00e8tre b se heurte \u00e0 la but\u00e9e infranchissable d\u00e9finie par la contrainte b &ge; x<sub>(n)<\/sub>. Diminuer b en dessous de x<sub>(n)<\/sub> exclurait la plus grande valeur observ\u00e9e et d\u00e9truirait imm\u00e9diatement la vraisemblance. La valeur admissible minimisant la perte de vraisemblance est donc la borne inf\u00e9rieure de son intervalle de validit\u00e9 [x<sub>(n)<\/sub>, +&infin;), \u00e0 savoir : b* = x<sub>(n)<\/sub>.<\/p>\n<p>Le principe fondamental r\u00e9gissant l&rsquo;estimation par maximum de vraisemblance pour la distribution uniforme se r\u00e9sume donc par une r\u00e8gle g\u00e9om\u00e9trique : pour maximiser la densit\u00e9 de probabilit\u00e9 constante 1\/(b &#8211; a) assign\u00e9e \u00e0 chaque observation, il est n\u00e9cessaire et suffisant de resserrer le support de distribution de mani\u00e8re optimale jusqu&rsquo;\u00e0 ce qu&rsquo;il vienne \u00e9pouser \u00e9troitement les donn\u00e9es les plus extr\u00eames du nuage de points sans en exclure aucune. L&rsquo;optimum est donc un optimum purement fronti\u00e8re, situ\u00e9 au point de contact entre l&rsquo;int\u00e9rieur r\u00e9gulier et la r\u00e9gion de nullit\u00e9 indicatrice.<\/p>\n<h2>6. D\u00e9rivation formelle des estimateurs du maximum de vraisemblance<\/h2>\n<h3>6.1 R\u00e9solution pour le cas g\u00e9n\u00e9ral \u00e0 deux param\u00e8tres U(a, b)<\/h3>\n<p>La synth\u00e8se des d\u00e9rivations pr\u00e9c\u00e9dentes permet d&rsquo;\u00e9tablir formellement le r\u00e9sultat central de l&rsquo;estimation de la loi uniforme bilat\u00e9rale continue g\u00e9n\u00e9rale. Soit un \u00e9chantillon al\u00e9atoire i.i.d. X = (X<sub>1<\/sub>, &#8230;, X<sub>n<\/sub>) issu d&rsquo;une loi U(a, b), avec les param\u00e8tres inconnus a &lt; b. Le probl\u00e8me d&rsquo;optimisation s&rsquo;\u00e9crit formellement sous la forme d&rsquo;un programme d&rsquo;optimisation math\u00e9matique :<\/p>\n<p>Maximiser &amp;ell;(a, b) = -n ln(b &#8211; a) &emsp; sous les contraintes : a &le; x<sub>(1)<\/sub> &emsp; et &emsp; b &ge; x<sub>(n)<\/sub><\/p>\n<p>Puisque la fonction objectif est strictement d\u00e9croissante par rapport \u00e0 l&rsquo;\u00e9tendue (b &#8211; a), la r\u00e9solution de ce probl\u00e8me de maximisation est rigoureusement isomorphe au probl\u00e8me de minimisation de la distance g\u00e9om\u00e9trique sous-jacente :<\/p>\n<p>Minimiser &Delta;(a, b) = b &#8211; a &emsp; sous les contraintes : a &le; x<sub>(1)<\/sub> &emsp; et &emsp; b &ge; x<sub>(n)<\/sub><\/p>\n<p>Soit (a, b) un couple quelconque appartenant au domaine admissible &Omega;<sub>x<\/sub>. Par application directe de la transitivit\u00e9 de l&rsquo;ordre sur les r\u00e9els :<\/p>\n<p>a &le; x<sub>(1)<\/sub> &rArr; -a &ge; -x<sub>(1)<\/sub><\/p>\n<p>En sommant cette in\u00e9galit\u00e9 avec la contrainte b &ge; x<sub>(n)<\/sub>, nous obtenons l&rsquo;in\u00e9galit\u00e9 fondamentale de minoration globale :<\/p>\n<p>b &#8211; a &ge; x<sub>(n)<\/sub> &#8211; x<sub>(1)<\/sub> &emsp; &forall; (a, b) &isin; &Omega;<sub>x<\/sub><\/p>\n<p>L&rsquo;\u00e9galit\u00e9 dans cette relation d&rsquo;ordre est satisfaite si et seulement si les deux contraintes sont satur\u00e9es simultan\u00e9ment, ce qui impose de fa\u00e7on univoque :<\/p>\n<p>a = x<sub>(1)<\/sub> &emsp; et &emsp; b = x<sub>(n)<\/sub><\/p>\n<p>Nous d\u00e9clarons formellement les estimateurs du maximum de vraisemblance (EMV), en tant que variables al\u00e9atoires d\u00e9pendant de l&rsquo;\u00e9chantillon X, par les expressions analytiques exactes :<\/p>\n<p>&acirc;<sub>EMV<\/sub> = X<sub>(1)<\/sub> = min(X<sub>1<\/sub>, X<sub>2<\/sub>, &#8230;, X<sub>n<\/sub>)<\/p>\n<p>b&#770;<sub>EMV<\/sub> = X<sub>(n)<\/sub> = max(X<sub>1<\/sub>, X<sub>2<\/sub>, &#8230;, X<sub>n<\/sub>)<\/p>\n<p>Le couple (X<sub>(1)<\/sub>, X<sub>(n)<\/sub>) constitue l&rsquo;unique solution globale de maximisation de la vraisemblance.<\/p>\n<h3>6.2 R\u00e9solution pour le cas unilat\u00e9ral standard U(0, theta)<\/h3>\n<p>Consid\u00e9rons \u00e0 pr\u00e9sent la variante param\u00e9trique unilat\u00e9rale canonique, dans laquelle la borne inf\u00e9rieure est fix\u00e9e de mani\u00e8re d\u00e9terministe \u00e0 l&rsquo;origine (a = 0), tandis que la borne sup\u00e9rieure constitue le seul param\u00e8tre inconnu &theta; &gt; 0 \u00e0 estimer : X ~ U(0, &theta;). Cette configuration est omnipr\u00e9sente dans les probl\u00e8mes de calibration industrielle, de temps d&rsquo;attente initiaux, ou de mod\u00e9lisation physique de contraintes unidirectionnelles.<\/p>\n<p>La fonction de vraisemblance conjointe se simplifie alors et s&rsquo;exprime sous la forme :<\/p>\n<p>L(&theta; ; x) = [ 1 \/ &theta;<sup>n<\/sup> ] &times; &prod;<sub>i=1<\/sub><sup>n<\/sup> <strong>1<\/strong><sub>[0, &theta;]<\/sub>(x<sub>i<\/sub>)<\/p>\n<p>Sous l&rsquo;hypoth\u00e8se physique ou empirique selon laquelle toutes les r\u00e9alisations x<sub>i<\/sub> sont strictement positives (x<sub>i<\/sub> &ge; 0 pour tout i), la condition &forall; i, x<sub>i<\/sub> &le; &theta; se r\u00e9sume \u00e9l\u00e9gamment par l&rsquo;in\u00e9galit\u00e9 unilat\u00e9rale sur la statistique d&rsquo;ordre maximale : x<sub>(n)<\/sub> &le; &theta;. D\u00e8s lors, la fonction de vraisemblance unilat\u00e9rale devient :<\/p>\n<p>L(&theta; ; x) = &theta;<sup>-n<\/sup> &times; <strong>1<\/strong><sub>[&theta; &ge; x<sub>(n)<\/sub>]<\/sub>(&theta;)<\/p>\n<p>Sur l&rsquo;intervalle admissible ferm\u00e9 [x<sub>(n)<\/sub>, +&infin;), la fonction &theta; &rarr; &theta;<sup>-n<\/sup> est une fonction continue, strictement d\u00e9croissante de mani\u00e8re monotone, dont la d\u00e9riv\u00e9e premi\u00e8re dL\/d&theta; = -n &theta;<sup>-(n+1)<\/sup> est strictement n\u00e9gative pour tout &theta; &gt; 0. Pour tout &theta; &lt; x<sub>(n)<\/sub>, la fonction de vraisemblance est rigoureusement nulle. Le profil de L(&theta; ; x) pr\u00e9sente donc une discontinuit\u00e9 de saut ou un pic abrupt pr\u00e9cis\u00e9ment en &theta; = x<sub>(n)<\/sub>, avant de d\u00e9cro\u00eetre asymptotiquement vers z\u00e9ro lorsque &theta; tend vers +&infin;.<\/p>\n<p>Il en r\u00e9sulte de mani\u00e8re formelle et sans \u00e9quivoque que le maximum global de vraisemblance est atteint \u00e0 la borne gauche exacte de l&rsquo;intervalle d&rsquo;admissibilit\u00e9 :<\/p>\n<p>&theta;&#770;<sub>EMV<\/sub> = X<sub>(n)<\/sub> = max(X<sub>1<\/sub>, X<sub>2<\/sub>, &#8230;, X<sub>n<\/sub>)<\/p>\n<p>L&rsquo;estimation optimale du maximum de vraisemblance d&rsquo;une \u00e9tendue unilat\u00e9rale est donc strictement \u00e9gale \u00e0 la plus grande observation recueillie dans l&rsquo;\u00e9chantillon.<\/p>\n<h3>6.3 R\u00e9solution pour la loi uniforme sym\u00e9trique U(-theta, theta)<\/h3>\n<p>Une autre d\u00e9clinaison structurelle fr\u00e9quente en sciences physiques et en m\u00e9trologie est le mod\u00e8le uniforme centr\u00e9 et sym\u00e9trique autour de l&rsquo;origine, not\u00e9 X ~ U(-&theta;, &theta;), o\u00f9 le param\u00e8tre d&rsquo;\u00e9chelle inconnu est &theta; &gt; 0. Dans ce contexte, l&rsquo;amplitude totale de la distribution vaut (b &#8211; a) = &theta; &#8211; (-&theta;) = 2&theta;.<\/p>\n<p>Pour un n-\u00e9chantillon i.i.d., la densit\u00e9 marginale de chaque variable est f(x ; &theta;) = (2&theta;)<sup>-1<\/sup> &times; <strong>1<\/strong><sub>[-&theta;, &theta;]<\/sub>(x). La fonction de vraisemblance conjointe s&rsquo;\u00e9crit :<\/p>\n<p>L(&theta; ; x) = [ 1 \/ (2&theta;) ]<sup>n<\/sup> &times; &prod;<sub>i=1<\/sub><sup>n<\/sup> <strong>1<\/strong><sub>[-&theta;, &theta;]<\/sub>(x<sub>i<\/sub>) = (2&theta;)<sup>-n<\/sup> &times; <strong>1<\/strong><sub>{ &forall; i &isin; {1, &#8230;, n}, -&theta; &le; x<sub>i<\/sub> &le; &theta; }<\/sub><\/p>\n<p>La condition d&rsquo;appartenance simultan\u00e9e -&theta; &le; x<sub>i<\/sub> &le; &theta; \u00e9quivaut rigoureusement \u00e0 l&rsquo;in\u00e9galit\u00e9 sur la valeur absolue de chaque observation : |x<sub>i<\/sub>| &le; &theta;. En \u00e9tendant cette condition \u00e0 l&rsquo;int\u00e9gralit\u00e9 de l&rsquo;\u00e9chantillon, la conjonction logique de ces n in\u00e9galit\u00e9s impose que le param\u00e8tre &theta; doive \u00eatre sup\u00e9rieur ou \u00e9gal au maximum des valeurs absolues de l&rsquo;ensemble des donn\u00e9es observ\u00e9es :<\/p>\n<p>&theta; &ge; max( |x<sub>1<\/sub>|, |x<sub>2<\/sub>|, &#8230;, |x<sub>n<\/sub>| ) = max( -x<sub>(1)<\/sub>, x<sub>(n)<\/sub> )<\/p>\n<p>En analysant \u00e0 nouveau le comportement de la fonction (2&theta;)<sup>-n<\/sup>, celle-ci appara\u00eet strictement d\u00e9croissante par rapport \u00e0 &theta; sur son domaine de d\u00e9finition strictement positif. Elle atteint par cons\u00e9quent son supremum absolu pour la plus petite valeur de &theta; respectant la contrainte indicatrice globale. L&rsquo;estimateur du maximum de vraisemblance pour le mod\u00e8le uniforme sym\u00e9trique centr\u00e9 est donc formalis\u00e9 par :<\/p>\n<p>&theta;&#770;<sub>EMV<\/sub> = max<sub>1&le;i&le;n<\/sub> |X<sub>i<\/sub>| = max( -X<sub>(1)<\/sub>, X<sub>(n)<\/sub> )<\/p>\n<p>Ce r\u00e9sultat d\u00e9montre la capacit\u00e9 du principe de vraisemblance \u00e0 projeter l&rsquo;information sur la statistique exhaustive idoine en accord avec la g\u00e9om\u00e9trie de sym\u00e9trie du probl\u00e8me initial.<\/p>\n<h2>7. Distribution d&rsquo;\u00e9chantillonnage exacte des statistiques d&rsquo;ordre<\/h2>\n<h3>7.1 Loi marginale de la statistique d&rsquo;ordre maximale X_(n)<\/h3>\n<p>Pour d\u00e9terminer les propri\u00e9t\u00e9s statistiques intrins\u00e8ques de l&rsquo;estimateur b&#770;<sub>EMV<\/sub> = X<sub>(n)<\/sub>, il est indispensable de d\u00e9river sa distribution d&rsquo;\u00e9chantillonnage exacte pour une taille d&rsquo;\u00e9chantillon finie n. Par d\u00e9finition, l&rsquo;\u00e9v\u00e9nement al\u00e9atoire {X<sub>(n)<\/sub> &le; x} se r\u00e9alise si et seulement si l&rsquo;ensemble des n variables al\u00e9atoires de l&rsquo;\u00e9chantillon r\u00e9alisent simultan\u00e9ment des valeurs inf\u00e9rieures ou \u00e9gales \u00e0 x. En exploitant l&rsquo;ind\u00e9pendance mutuelle et le fait que les variables soient identiquement distribu\u00e9es selon la loi parente U(a, b) :<\/p>\n<p>F<sub>X<sub>(n)<\/sub><\/sub>(x) = P(X<sub>(n)<\/sub> &le; x) = P(X<sub>1<\/sub> &le; x, X<sub>2<\/sub> &le; x, &#8230;, X<sub>n<\/sub> &le; x) = &prod;<sub>i=1<\/sub><sup>n<\/sup> P(X<sub>i<\/sub> &le; x) = [ F<sub>X<\/sub>(x) ]<sup>n<\/sup><\/p>\n<p>En substituant la fonction de r\u00e9partition univari\u00e9e de la loi uniforme continue sur le support [a, b], nous obtenons l&rsquo;expression exacte de la fonction de r\u00e9partition de la statistique d&rsquo;ordre maximale :<\/p>\n<p>F<sub>X<sub>(n)<\/sub><\/sub>(x) = 0 si x &lt; a ; [ (x &#8211; a) \/ (b &#8211; a) ]<sup>n<\/sup> si a &le; x &le; b ; 1 si x &gt; b.<\/p>\n<p>La fonction de densit\u00e9 de probabilit\u00e9 marginale f<sub>X<sub>(n)<\/sub><\/sub>(x) s&rsquo;obtient par d\u00e9rivation continue de F<sub>X<sub>(n)<\/sub><\/sub>(x) par rapport \u00e0 x sur l&rsquo;intervalle ouvert (a, b) :<\/p>\n<p>f<sub>X<sub>(n)<\/sub><\/sub>(x) = d\/dx [ (x &#8211; a) \/ (b &#8211; a) ]<sup>n<\/sup> = [ n (x &#8211; a)<sup>n &#8211; 1<\/sup> \/ (b &#8211; a)<sup>n<\/sup> ] &times; <strong>1<\/strong><sub>[a, b]<\/sub>(x)<\/p>\n<p>Cette distribution correspond structurellement \u00e0 une loi b\u00eata g\u00e9n\u00e9ralis\u00e9e r\u00e9\u00e9chelonn\u00e9e et translat\u00e9e sur l&rsquo;intervalle [a, b], de param\u00e8tres de forme &alpha; = n et &beta; = 1. Pour toute taille d&rsquo;\u00e9chantillon n &gt; 1, cette densit\u00e9 est strictement croissante et convexe par rapport \u00e0 x, concentrant une densit\u00e9 de probabilit\u00e9 de plus en plus massive au voisinage imm\u00e9diat de la fronti\u00e8re sup\u00e9rieure b au fur et \u00e0 mesure que la taille d&rsquo;\u00e9chantillon n augmente.<\/p>\n<h3>7.2 Loi marginale de la statistique d&rsquo;ordre minimale X_(1)<\/h3>\n<p>De fa\u00e7on sym\u00e9trique, examinons le comportement stochastique de l&rsquo;estimateur de la borne inf\u00e9rieure &acirc;<sub>EMV<\/sub> = X<sub>(1)<\/sub>. La d\u00e9rivation de sa fonction de r\u00e9partition repose sur l&rsquo;analyse de sa fonction de survie conjointe. L&rsquo;\u00e9v\u00e9nement {X<sub>(1)<\/sub> &gt; x} est rigoureusement \u00e9quivalent \u00e0 la conjonction des \u00e9v\u00e9nements ind\u00e9pendants {X<sub>i<\/sub> &gt; x} pour chaque observation i &isin; {1, &#8230;, n} :<\/p>\n<p>P(X<sub>(1)<\/sub> &gt; x) = P(X<sub>1<\/sub> &gt; x, X<sub>2<\/sub> &gt; x, &#8230;, X<sub>n<\/sub> &gt; x) = &prod;<sub>i=1<\/sub><sup>n<\/sup> P(X<sub>i<\/sub> &gt; x) = [ 1 &#8211; F<sub>X<\/sub>(x) ]<sup>n<\/sup><\/p>\n<p>Il en d\u00e9coule l&rsquo;expression de la fonction de r\u00e9partition cumulative de la statistique d&rsquo;ordre minimale :<\/p>\n<p>F<sub>X<sub>(1)<\/sub><\/sub>(x) = 1 &#8211; P(X<sub>(1)<\/sub> &gt; x) = 1 &#8211; [ 1 &#8211; F<sub>X<\/sub>(x) ]<sup>n<\/sup><\/p>\n<p>En rempla\u00e7ant F<sub>X<\/sub>(x) par son expression affine pour x &isin; [a, b], nous avons 1 &#8211; F<sub>X<\/sub>(x) = (b &#8211; x) \/ (b &#8211; a), d&rsquo;o\u00f9 :<\/p>\n<p>F<sub>X<sub>(1)<\/sub><\/sub>(x) = 1 &#8211; [ (b &#8211; x) \/ (b &#8211; a) ]<sup>n<\/sup> &emsp; pour x &isin; [a, b]<\/p>\n<p>En d\u00e9rivant cette fonction par rapport \u00e0 x, nous \u00e9tablissons la densit\u00e9 marginale exacte de X<sub>(1)<\/sub> :<\/p>\n<p>f<sub>X<sub>(1)<\/sub><\/sub>(x) = d\/dx { 1 &#8211; [ (b &#8211; x) \/ (b &#8211; a) ]<sup>n<\/sup> } = &#8211; n [ (b &#8211; x) \/ (b &#8211; a) ]<sup>n &#8211; 1<\/sup> &times; (-1 \/ (b &#8211; a)) = [ n (b &#8211; x)<sup>n &#8211; 1<\/sup> \/ (b &#8211; a)<sup>n<\/sup> ] &times; <strong>1<\/strong><sub>[a, b]<\/sub>(x)<\/p>\n<p>Cette loi de probabilit\u00e9 s&rsquo;identifie \u00e0 une loi b\u00eata r\u00e9\u00e9chelonn\u00e9e de param\u00e8tres &alpha; = 1 et &beta; = n. Elle est strictement d\u00e9croissante sur [a, b] pour n &gt; 1, accumulant son poids probabiliste \u00e0 proximit\u00e9 imm\u00e9diate de la borne inf\u00e9rieure th\u00e9orique a. Il existe une dualit\u00e9 sym\u00e9trique parfaite entre la loi de X<sub>(1)<\/sub> et celle de X<sub>(n)<\/sub> par rapport au centre de gravit\u00e9 &mu; = (a + b) \/ 2.<\/p>\n<h3>7.3 Loi conjointe du couple (X_(1), X_(n)) et \u00e9tendue d&rsquo;\u00e9chantillon<\/h3>\n<p>Puisque les statistiques d&rsquo;ordre extr\u00eames X<sub>(1)<\/sub> et X<sub>(n)<\/sub> sont d\u00e9riv\u00e9es du m\u00eame \u00e9chantillon fini, elles ne sont absolument pas ind\u00e9pendantes sur le plan probabiliste. L&rsquo;\u00e9valuation de leur comportement conjoint n\u00e9cessite la formulation de leur densit\u00e9 conjointe bidimensionnelle f<sub>X<sub>(1)<\/sub>, X<sub>(n)<\/sub><\/sub>(u, v). Selon la formule g\u00e9n\u00e9rale des lois jointes des statistiques d&rsquo;ordre pour un \u00e9chantillon continu de taille n &ge; 2 :<\/p>\n<p>f<sub>X<sub>(1)<\/sub>, X<sub>(n)<\/sub><\/sub>(u, v) = n (n &#8211; 1) [ F(v) &#8211; F(u) ]<sup>n &#8211; 2<\/sup> f(u) f(v) &times; <strong>1<\/strong><sub>{u &le; v}<\/sub><\/p>\n<p>En substituant la fonction de densit\u00e9 univari\u00e9e f(t) = 1 \/ (b &#8211; a) et la fonction de r\u00e9partition F(t) de la loi uniforme continue sur le triangle ferm\u00e9 d\u00e9fini par le support a &le; u &le; v &le; b :<\/p>\n<p>f<sub>X<sub>(1)<\/sub>, X<sub>(n)<\/sub><\/sub>(u, v) = [ n (n &#8211; 1) (v &#8211; u)<sup>n &#8211; 2<\/sup> \/ (b &#8211; a)<sup>n<\/sup> ] &times; <strong>1<\/strong><sub>{a &le; u &le; v &le; b}<\/sub><\/p>\n<p>Cette expression conjointe est fondamentale pour \u00e9tudier la corr\u00e9lation stochastique r\u00e9siduelle entre les deux fronti\u00e8res empiriques.<\/p>\n<p>\u00c0 partir de cette distribution conjointe, nous pouvons op\u00e9rer un changement de variables bijectif pour caract\u00e9riser la loi de probabilit\u00e9 de l&rsquo;\u00e9tendue \u00e9chantillonnale R = X<sub>(n)<\/sub> &#8211; X<sub>(1)<\/sub> et celle du milieu d&rsquo;\u00e9chantillon M = (X<sub>(1)<\/sub> + X<sub>(n)<\/sub>) \/ 2. Par int\u00e9gration marginale sur le jacobien de la transformation, la fonction de densit\u00e9 de l&rsquo;\u00e9tendue empirique R sur son support [0, b &#8211; a] est r\u00e9gie par l&rsquo;expression :<\/p>\n<p>f<sub>R<\/sub>(r) = [ n (n &#8211; 1) r<sup>n &#8211; 2<\/sup> (b &#8211; a &#8211; r) \/ (b &#8211; a)<sup>n<\/sup> ] &times; <strong>1<\/strong><sub>[0, b &#8211; a]<\/sub>(r)<\/p>\n<p>De m\u00eame, le milieu d&rsquo;\u00e9chantillon M, qui constitue un estimateur intuitif de l&rsquo;esp\u00e9rance &mu; = (a + b) \/ 2, pr\u00e9sente une distribution hautement leptokurtique sym\u00e9trique autour de ce barycentre, avec une variance asymptotiquement inf\u00e9rieure \u00e0 celle de la moyenne arithm\u00e9tique standard de l&rsquo;\u00e9chantillon.<\/p>\n<h2>8. Propri\u00e9t\u00e9s statistiques fondamentales : biais, variance et efficacit\u00e9<\/h2>\n<h3>8.1 Calcul analytique du biais des estimateurs EMV<\/h3>\n<p>L&rsquo;une des tares historiques associ\u00e9es aux estimateurs du maximum de vraisemblance est leur propension \u00e0 introduire un biais syst\u00e9matique pour les \u00e9chantillons de taille finie. Calculons l&rsquo;esp\u00e9rance math\u00e9matique formelle de l&rsquo;estimateur b&#770;<sub>EMV<\/sub> = X<sub>(n)<\/sub> en int\u00e9grant sa fonction de densit\u00e9 marginale d\u00e9cente :<\/p>\n<p>E[b&#770;<sub>EMV<\/sub>] = &int;<sub>a<\/sub><sup>b<\/sup> x f<sub>X<sub>(n)<\/sub><\/sub>(x) dx = &int;<sub>a<\/sub><sup>b<\/sup> x [ n (x &#8211; a)<sup>n &#8211; 1<\/sup> \/ (b &#8211; a)<sup>n<\/sup> ] dx<\/p>\n<p>Op\u00e9rons le changement de variable affine lin\u00e9aire t = x &#8211; a, impliquant x = a + t et dx = dt, avec un intervalle d&rsquo;int\u00e9gration mutant de [a, b] vers [0, b &#8211; a] :<\/p>\n<p>E[b&#770;<sub>EMV<\/sub>] = [ n \/ (b &#8211; a)<sup>n<\/sup> ] &int;<sub>0<\/sub><sup>b &#8211; a<\/sup> (a + t) t<sup>n &#8211; 1<\/sup> dt = [ n \/ (b &#8211; a)<sup>n<\/sup> ] [ a &int;<sub>0<\/sub><sup>b &#8211; a<\/sup> t<sup>n &#8211; 1<\/sup> dt + &int;<sub>0<\/sub><sup>b &#8211; a<\/sup> t<sup>n<\/sup> dt ]<\/p>\n<p>L&rsquo;\u00e9valuation analytique des deux int\u00e9grales mon\u00f4miales donne imm\u00e9diatement :<\/p>\n<p>&int;<sub>0<\/sub><sup>b &#8211; a<\/sup> t<sup>n &#8211; 1<\/sup> dt = (b &#8211; a)<sup>n<\/sup> \/ n &emsp; et &emsp; &int;<sub>0<\/sub><sup>b &#8211; a<\/sup> t<sup>n<\/sup> dt = (b &#8211; a)<sup>n + 1<\/sup> \/ (n + 1)<\/p>\n<p>En recombinant ces termes par substitution :<\/p>\n<p>E[b&#770;<sub>EMV<\/sub>] = [ n \/ (b &#8211; a)<sup>n<\/sup> ] [ a (b &#8211; a)<sup>n<\/sup> \/ n + (b &#8211; a)<sup>n + 1<\/sup> \/ (n + 1) ] = a + [ n \/ (n + 1) ] (b &#8211; a)<\/p>\n<p>Le biais formel de l&rsquo;estimateur de la borne sup\u00e9rieure, d\u00e9fini par B(b&#770;) = E[b&#770;] &#8211; b, s&rsquo;\u00e9tablit donc comme :<\/p>\n<p>B(b&#770;<sub>EMV<\/sub>) = a + [ n \/ (n + 1) ] (b &#8211; a) &#8211; b = [ n \/ (n + 1) &#8211; 1 ] (b &#8211; a) = &#8211; (b &#8211; a) \/ (n + 1)<\/p>\n<p>Puisque b &gt; a et que n &ge; 1, ce biais est strictement n\u00e9gatif : B(b&#770;<sub>EMV<\/sub>) &lt; 0. L&rsquo;EMV sous-estime donc syst\u00e9matiquement la v\u00e9ritable borne sup\u00e9rieure de la population. Une d\u00e9rivation rigoureusement sym\u00e9trique men\u00e9e sur la statistique d&rsquo;ordre minimale X<sub>(1)<\/sub> conduit \u00e0 l&rsquo;esp\u00e9rance de l&rsquo;estimateur de la borne inf\u00e9rieure :<\/p>\n<p>E[&acirc;<sub>EMV<\/sub>] = a + [ 1 \/ (n + 1) ] (b &#8211; a)<\/p>\n<p>Le biais de la borne inf\u00e9rieure s&rsquo;exprime donc par B(&acirc;<sub>EMV<\/sub>) = + (b &#8211; a) \/ (n + 1) &gt; 0. L&rsquo;estimateur de la borne inf\u00e9rieure surestime donc structurellement la v\u00e9ritable fronti\u00e8re basse a. En combinant ces deux r\u00e9sultats, l&rsquo;\u00e9tendue empirique moyenne de l&rsquo;\u00e9chantillon vaut :<\/p>\n<p>E[X<sub>(n)<\/sub> &#8211; X<sub>(1)<\/sub>] = E[b&#770;] &#8211; E[&acirc;] = [ (n &#8211; 1) \/ (n + 1) ] (b &#8211; a) &lt; (b &#8211; a)<\/p>\n<p>Les estimateurs bruts du maximum de vraisemblance capturent une amplitude qui sous-\u00e9value intrins\u00e8quement l&rsquo;\u00e9tendue r\u00e9elle du support de distribution.<\/p>\n<h3>8.2 Construction d&rsquo;estimateurs sans biais \u00e0 variance minimale (UMVU)<\/h3>\n<p>L&rsquo;identification pr\u00e9cise de la structure de ces biais permet d&rsquo;en \u00e9laborer la correction analytique exacte pour concevoir des estimateurs sans biais \u00e0 variance minimale uniforme (<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Estimateur_sans_biais_%C3%A0_variance_minimale\">estimateurs UMVU<\/a>), en s&rsquo;appuyant sur les puissants th\u00e9or\u00e8mes de Rao-Blackwell et de Lehmann-Scheff\u00e9.<\/p>\n<p>Dans le cadre unilat\u00e9ral X ~ U(0, &theta;), l&rsquo;esp\u00e9rance de l&rsquo;estimateur s&rsquo;\u00e9crit de fa\u00e7on simplifi\u00e9e : E[X<sub>(n)<\/sub>] = [n \/ (n + 1)] &theta;. Il appara\u00eet imm\u00e9diatement qu&rsquo;une simple mise \u00e0 l&rsquo;\u00e9chelle scalaire par le facteur inverse permet d&rsquo;annihiler totalement le biais pour toute taille d&rsquo;\u00e9chantillon finie. D\u00e9finissons le nouvel estimateur corrig\u00e9 :<\/p>\n<p>&theta;&#770;<sub>UMVU<\/sub> = [ (n + 1) \/ n ] &times; X<sub>(n)<\/sub> = [ 1 + 1\/n ] &times; X<sub>(n)<\/sub><\/p>\n<p>Puisque la statistique X<sub>(n)<\/sub> est non seulement exhaustive (suffisante) mais \u00e9galement <em>compl\u00e8te<\/em> pour la famille de lois {U(0, &theta;) : &theta; &gt; 0}, le <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Th%C3%A9or%C3%A8me_de_Lehmann-Scheff%C3%A9\">th\u00e9or\u00e8me de Lehmann-Scheff\u00e9<\/a> assure que &theta;&#770;<sub>UMVU<\/sub> est l&rsquo;unique estimateur non biais\u00e9 de &theta; qui minimise la variance parmi absolument tous les estimateurs sans biais imaginables.<\/p>\n<p>Pour le mod\u00e8le bilat\u00e9ral g\u00e9n\u00e9ral U(a, b), la correction conjointe des biais des deux bornes s&rsquo;effectue en r\u00e9solvant le syst\u00e8me d&rsquo;\u00e9quations lin\u00e9aires form\u00e9 par les esp\u00e9rances des statistiques d&rsquo;ordre extr\u00eames. Les estimateurs non biais\u00e9s \u00e0 variance minimale uniforme des param\u00e8tres a et b s&rsquo;expriment respectivement par les combinaisons lin\u00e9aires des deux statistiques d&rsquo;ordre extr\u00eames :<\/p>\n<p>&acirc;<sub>UMVU<\/sub> = [ (n &times; X<sub>(1)<\/sub> &#8211; X<sub>(n)<\/sub>) \/ (n &#8211; 1) ] = X<sub>(1)<\/sub> &#8211; [ (X<sub>(n)<\/sub> &#8211; X<sub>(1)<\/sub>) \/ (n &#8211; 1) ]<\/p>\n<p>b&#770;<sub>UMVU<\/sub> = [ (n &times; X<sub>(n)<\/sub> &#8211; X<sub>(1)<\/sub>) \/ (n &#8211; 1) ] = X<sub>(n)<\/sub> + [ (X<sub>(n)<\/sub> &#8211; X<sub>(1)<\/sub>) \/ (n &#8211; 1) ]<\/p>\n<p>L&rsquo;interpr\u00e9tation g\u00e9om\u00e9trique de cette correction est remarquable d&rsquo;\u00e9l\u00e9gance : pour supprimer le biais d&rsquo;\u00e9chantillonnage, le statisticien prend les fronti\u00e8res empiriques observ\u00e9es et extrapole le support vers l&rsquo;ext\u00e9rieur en ajoutant \u00e0 chaque extr\u00e9mit\u00e9 une fraction 1\/(n &#8211; 1) de l&rsquo;\u00e9tendue observ\u00e9e.<\/p>\n<h3>8.3 Calcul de la variance et de l&rsquo;erreur quadratique moyenne (EQM)<\/h3>\n<p>Pour juger rigoureusement de la performance pr\u00e9dictive d&rsquo;un estimateur, la th\u00e9orie statistique moderne \u00e9value l&rsquo;erreur quadratique moyenne (EQM), d\u00e9finie par la d\u00e9composition canonique : EQM(&theta;&#770;) = Var(&theta;&#770;) + [B(&theta;&#770;)]<sup>2<\/sup>. Afin d&rsquo;\u00e9tablir l&rsquo;EQM de l&rsquo;EMV dans le cadre U(0, &theta;), calculons en premier lieu le moment d&rsquo;ordre deux de la statistique d&rsquo;ordre X<sub>(n)<\/sub> :<\/p>\n<p>E[X<sub>(n)<\/sub><sup>2<\/sup>] = &int;<sub>0<\/sub><sup>&theta;<\/sup> x<sup>2<\/sup> [ n x<sup>n &#8211; 1<\/sup> \/ &theta;<sup>n<\/sup> ] dx = [ n \/ &theta;<sup>n<\/sup> ] &int;<sub>0<\/sub><sup>&theta;<\/sup> x<sup>n + 1<\/sup> dx = [ n \/ &theta;<sup>n<\/sup> ] [ &theta;<sup>n + 2<\/sup> \/ (n + 2) ] = [ n \/ (n + 2) ] &theta;<sup>2<\/sup><\/p>\n<p>La variance exacte de l&rsquo;estimateur brut du maximum de vraisemblance s&rsquo;en d\u00e9duit par soustraction alg\u00e9brique :<\/p>\n<p>Var(X<sub>(n)<\/sub>) = E[X<sub>(n)<\/sub><sup>2<\/sup>] &#8211; [ E[X<sub>(n)<\/sub>] ]<sup>2<\/sup> = [ n \/ (n + 2) ] &theta;<sup>2<\/sup> &#8211; [ n \/ (n + 1) ]<sup>2<\/sup> &theta;<sup>2<\/sup> = n &theta;<sup>2<\/sup> [ (n + 1)<sup>2<\/sup> &#8211; n(n + 2) ] \/ [ (n + 2)(n + 1)<sup>2<\/sup> ]<\/p>\n<p>Puisque le d\u00e9veloppement polynomial donne (n + 1)<sup>2<\/sup> &#8211; n(n + 2) = n<sup>2<\/sup> + 2n + 1 &#8211; n<sup>2<\/sup> &#8211; 2n = 1, nous aboutissons \u00e0 la variance de l&rsquo;estimateur EMV :<\/p>\n<p>Var(&theta;&#770;<sub>EMV<\/sub>) = n &theta;<sup>2<\/sup> \/ [ (n + 1)<sup>2<\/sup> (n + 2) ]<\/p>\n<p>L&rsquo;erreur quadratique moyenne de l&rsquo;estimateur brut du maximum de vraisemblance s&rsquo;\u00e9crit donc :<\/p>\n<p>EQM(&theta;&#770;<sub>EMV<\/sub>) = Var(X<sub>(n)<\/sub>) + [B(X<sub>(n)<\/sub>)]<sup>2<\/sup> = [ n &theta;<sup>2<\/sup> \/ ((n + 1)<sup>2<\/sup> (n + 2)) ] + [ &theta; \/ (n + 1) ]<sup>2<\/sup> = [ n &theta;<sup>2<\/sup> + (n + 2)&theta;<sup>2<\/sup> ] \/ [ (n + 1)<sup>2<\/sup> (n + 2) ] = 2 &theta;<sup>2<\/sup> \/ [ (n + 1) (n + 2) ]<\/p>\n<p>\u00c9valuons \u00e0 pr\u00e9sent la variance et l&rsquo;EQM de l&rsquo;estimateur sans biais &theta;&#770;<sub>UMVU<\/sub> = [(n + 1)\/n] X<sub>(n)<\/sub>. Puisque son biais est nul par construction, son EQM est strictement identique \u00e0 sa propre variance :<\/p>\n<p>EQM(&theta;&#770;<sub>UMVU<\/sub>) = Var( [ (n + 1) \/ n ] X<sub>(n)<\/sub> ) = [ (n + 1) \/ n ]<sup>2<\/sup> Var(X<sub>(n)<\/sub>) = [ (n + 1)<sup>2<\/sup> \/ n<sup>2<\/sup> ] &times; [ n &theta;<sup>2<\/sup> \/ ((n + 1)<sup>2<\/sup> (n + 2)) ] = &theta;<sup>2<\/sup> \/ [ n (n + 2) ]<\/p>\n<p>Comparons \u00e0 pr\u00e9sent les deux erreurs quadratiques moyennes en \u00e9valuant leur ratio relatif :<\/p>\n<p>EQM(&theta;&#770;<sub>EMV<\/sub>) \/ EQM(&theta;&#770;<sub>UMVU<\/sub>) = [ 2 &theta;<sup>2<\/sup> \/ ((n + 1)(n + 2)) ] \/ [ &theta;<sup>2<\/sup> \/ (n(n + 2)) ] = 2n \/ (n + 1)<\/p>\n<p>Ce r\u00e9sultat met en lumi\u00e8re une propri\u00e9t\u00e9 fondamentale :<\/p>\n<ul>\n<li>Pour un \u00e9chantillon unitaire n = 1 : le ratio vaut 2\/2 = 1 ; les deux estimateurs sont rigoureusement identiques et poss\u00e8dent la m\u00eame erreur quadratique moyenne.<\/li>\n<li>Pour toute taille d&rsquo;\u00e9chantillon n &gt; 1 : nous constatons que 2n &gt; n + 1, ce qui implique de mani\u00e8re p\u00e9remptoire que EQM(&theta;&#770;<sub>EMV<\/sub>) &gt; EQM(&theta;&#770;<sub>UMVU<\/sub>).<\/li>\n<\/ul>\n<p>Contrairement \u00e0 certains contextes gaussiens o\u00f9 l&rsquo;estimateur biais\u00e9 du maximum de vraisemblance surpasse l&rsquo;estimateur non biais\u00e9 en r\u00e9duisant la variance au d\u00e9triment d&rsquo;un biais mod\u00e9r\u00e9, ici l&rsquo;estimateur sans biais UMVU surpasse syst\u00e9matiquement l&rsquo;EMV brut en termes d&rsquo;erreur quadratique moyenne pour tout n &ge; 2. La correction du biais am\u00e9liore donc \u00e0 la fois la justesse et la pr\u00e9cision quadratique globale de l&rsquo;estimation.<\/p>\n<h2>9. Comportement asymptotique et taux de convergence super-efficace<\/h2>\n<h3>9.1 Consistance forte et convergence presque s\u00fbre<\/h3>\n<p>L&rsquo;analyse asymptotique s&rsquo;int\u00e9resse au comportement des estimateurs lorsque le nombre d&rsquo;observations n diverge vers l&rsquo;infini math\u00e9matique (+&infin;). La premi\u00e8re propri\u00e9t\u00e9 fondamentale \u00e0 d\u00e9montrer est la consistance forte de l&rsquo;estimateur du maximum de vraisemblance, c&rsquo;est-\u00e0-dire sa convergence presque s\u00fbre vers la v\u00e9ritable valeur du param\u00e8tre populationnel.<\/p>\n<p>Soit X<sub>1<\/sub>, X<sub>2<\/sub>, &#8230;, X<sub>n<\/sub> une suite infinie de variables i.i.d. distribu\u00e9es selon U(a, b). Nous souhaitons d\u00e9montrer que b&#770;<sub>n<\/sub> = X<sub>(n)<\/sub> converge presque s\u00fbrement vers la borne th\u00e9orique sup\u00e9rieure b, not\u00e9 formellement : X<sub>(n)<\/sub> &rarr;<sup>p.s.<\/sup> b lorsque n &rarr; &infin;. Fixons un \u00e9cart d&rsquo;erreur arbitrairement petit &epsilon; &gt; 0, tel que 0 &lt; &epsilon; &lt; b &#8211; a. En vertu de la forme analytique de la fonction de r\u00e9partition de la statistique maximale :<\/p>\n<p>P( |X<sub>(n)<\/sub> &#8211; b| &gt; &epsilon; ) = P( b &#8211; X<sub>(n)<\/sub> &gt; &epsilon; ) = P( X<sub>(n)<\/sub> &lt; b &#8211; &epsilon; ) = [ F<sub>X<\/sub>(b &#8211; &epsilon;) ]<sup>n<\/sup> = [ (b &#8211; &epsilon; &#8211; a) \/ (b &#8211; a) ]<sup>n<\/sup> = [ 1 &#8211; &epsilon; \/ (b &#8211; a) ]<sup>n<\/sup><\/p>\n<p>Posons la constante r\u00e9elle scalaire q = 1 &#8211; &epsilon; \/ (b &#8211; a). Puisque 0 &lt; &epsilon; &lt; b &#8211; a, il est garanti que 0 &le; q &lt; 1.<\/p>\n<p>Consid\u00e9rons \u00e0 pr\u00e9sent la somme infinie de ces probabilit\u00e9s d&rsquo;erreur sur l&rsquo;ensemble des entiers naturels n &ge; 1 :<\/p>\n<p>&sum;<sub>n=1<\/sub><sup>&infin;<\/sup> P( |X<sub>(n)<\/sub> &#8211; b| &gt; &epsilon; ) = &sum;<sub>n=1<\/sub><sup>&infin;<\/sup> q<sup>n<\/sup> = q \/ (1 &#8211; q) &lt; +&infin;<\/p>\n<p>Cette s\u00e9rie g\u00e9om\u00e9trique convergeant vers une valeur finie pour tout &epsilon; &gt; 0 fix\u00e9, le premier <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Lemme_de_Borel-Cantelli\">lemme de Borel-Cantelli<\/a> s&rsquo;applique directement. Il garantit que la probabilit\u00e9 que l&rsquo;\u00e9v\u00e9nement {|X<sub>(n)<\/sub> &#8211; b| &gt; &epsilon;} se produise une infinit\u00e9 de fois est strictement \u00e9gale \u00e0 z\u00e9ro. En termes probabilistes formels, cela \u00e9quivaut \u00e0 la convergence presque s\u00fbre :<\/p>\n<p>P( lim<sub>n&rarr;&infin;<\/sub> X<sub>(n)<\/sub> = b ) = 1 &emsp; et par stricte sym\u00e9trie, &emsp; P( lim<sub>n&rarr;&infin;<\/sub> X<sub>(1)<\/sub> = a ) = 1<\/p>\n<p>Les estimateurs bruts et corrig\u00e9s du maximum de vraisemblance sont donc fortement consistants.<\/p>\n<h3>9.2 Vitesse de convergence en n vs vitesse param\u00e9trique usuelle en racine de n<\/h3>\n<p>Le point le plus remarquable de cette th\u00e9orie asymptotique r\u00e9side dans la d\u00e9termination de la vitesse de convergence stochastique. Dans les mod\u00e8les statistiques r\u00e9guliers gouvern\u00e9s par les conditions de Cram\u00e9r-Rao (tels que les distributions gaussiennes, de Poisson ou exponentielles), les estimateurs convergent vers la cible th\u00e9orique \u00e0 la vitesse universelle &radic;n. Cela implique formellement que l&rsquo;erreur type de l&rsquo;estimation d\u00e9cro\u00eet en O(1\/&radic;n), n\u00e9cessitant de quadrupler la taille de l&rsquo;\u00e9chantillon pour diviser l&rsquo;incertitude d&rsquo;estimation par deux.<\/p>\n<p>Pour la distribution uniforme continue, la rupture des conditions de r\u00e9gularit\u00e9 induite par le support d\u00e9pendant des param\u00e8tres m\u00e9tamorphose enti\u00e8rement ce r\u00e9gime. En r\u00e9examinant la variance de l&rsquo;estimateur b&#770;<sub>EMV<\/sub> = X<sub>(n)<\/sub> \u00e9tablie \u00e0 la section pr\u00e9c\u00e9dente :<\/p>\n<p>Var(X<sub>(n)<\/sub>) = n (b &#8211; a)<sup>2<\/sup> \/ [ (n + 1)<sup>2<\/sup> (n + 2) ] &sim; (b &#8211; a)<sup>2<\/sup> \/ n<sup>2<\/sup> = O( 1 \/ n<sup>2<\/sup> ) lorsque n &rarr; &infin;<\/p>\n<p>L&rsquo;\u00e9cart-type de l&rsquo;estimateur du maximum de vraisemblance d\u00e9cro\u00eet donc \u00e0 l&rsquo;\u00e9chelle &sigma;(X<sub>(n)<\/sub>) = O(1\/n). L&rsquo;ordre de grandeur de l&rsquo;erreur d&rsquo;estimation est proportionnel \u00e0 1\/n et non plus \u00e0 1\/&radic;n.<\/p>\n<p>Ce ph\u00e9nom\u00e8ne exceptionnel est qualifi\u00e9 dans la litt\u00e9rature de <em>super-efficacit\u00e9<\/em> ou de convergence ultra-rapide \u00e0 vitesse lin\u00e9aire n. Les implications op\u00e9rationnelles sont consid\u00e9rables : pour diviser par dix l&rsquo;incertitude sur la localisation des fronti\u00e8res inconnues du support, il suffit de multiplier la taille d&rsquo;\u00e9chantillon par dix, alors qu&rsquo;une inf\u00e9rence standard \u00e0 vitesse &radic;n n\u00e9cessiterait de multiplier les effectifs par cent. Ce gain d&rsquo;efficacit\u00e9 r\u00e9sulte du fait que chaque nouvelle observation qui repousse les bornes extr\u00eames resserre instantan\u00e9ment l&rsquo;espace d&rsquo;admissibilit\u00e9 de fa\u00e7on d\u00e9terministe.<\/p>\n<h3>9.3 Distribution limite exponentielle et non-normalit\u00e9<\/h3>\n<p>La normalit\u00e9 asymptotique \u00e9tant d\u00e9faillante, la loi limite de l&rsquo;estimateur normalis\u00e9 ne converge pas vers la cloche gaussienne conventionnelle. Pour \u00e9tablir la v\u00e9ritable loi limite, nous devons renormaliser l&rsquo;\u00e9cart (b &#8211; X<sub>(n)<\/sub>) non pas par le facteur habituel &radic;n, mais par le facteur d&rsquo;\u00e9chelle lin\u00e9aire n.<\/p>\n<p>Consid\u00e9rons la variable al\u00e9atoire r\u00e9\u00e9chelonn\u00e9e Z<sub>n<\/sub> = n [ b &#8211; X<sub>(n)<\/sub> ] \/ (b &#8211; a). D\u00e9terminons la fonction de r\u00e9partition limite de cette variable pour tout z &ge; 0 lorsque n tend vers +&infin; :<\/p>\n<p>P( Z<sub>n<\/sub> &le; z ) = P( n [ b &#8211; X<sub>(n)<\/sub> ] \/ (b &#8211; a) &le; z ) = P( b &#8211; X<sub>(n)<\/sub> &le; z (b &#8211; a) \/ n ) = P( X<sub>(n)<\/sub> &ge; b &#8211; z (b &#8211; a) \/ n )<\/p>\n<p>En exprimant cette probabilit\u00e9 par le compl\u00e9mentaire de la fonction de r\u00e9partition de X<sub>(n)<\/sub> :<\/p>\n<p>P( Z<sub>n<\/sub> &le; z ) = 1 &#8211; F<sub>X<sub>(n)<\/sub><\/sub>( b &#8211; z (b &#8211; a) \/ n ) = 1 &#8211; [ ( (b &#8211; z(b &#8211; a)\/n) &#8211; a ) \/ (b &#8211; a) ]<sup>n<\/sup> = 1 &#8211; [ 1 &#8211; z \/ n ]<sup>n<\/sup><\/p>\n<p>En mobilisant la caract\u00e9risation analytique fondamentale du nombre exponentiel d&rsquo;Euler, lim<sub>n&rarr;&infin;<\/sub> [ 1 &#8211; z\/n ]<sup>n<\/sup> = e<sup>-z<\/sup> pour tout r\u00e9el z, nous obtenons imm\u00e9diatement la limite asymptotique :<\/p>\n<p>lim<sub>n&rarr;&infin;<\/sub> P( Z<sub>n<\/sub> &le; z ) = 1 &#8211; e<sup>-z<\/sup> &emsp; pour tout z &ge; 0<\/p>\n<p>Ce r\u00e9sultat d\u00e9montre avec \u00e9l\u00e9gance que la variable al\u00e9atoire normalis\u00e9e n (b &#8211; X<sub>(n)<\/sub>) converge en loi (distribu\u00e9e au sens faible) vers une <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Loi_exponentielle\">loi exponentielle standard<\/a> de param\u00e8tre d&rsquo;intensit\u00e9 &lambda; = 1 \/ (b &#8211; a) :<\/p>\n<p>n ( b &#8211; X<sub>(n)<\/sub> ) &rarr;<sup>d<\/sup> Exp( &lambda; = 1 \/ (b &#8211; a) )<\/p>\n<p>Cette convergence vers une loi exponentielle unilat\u00e9rale asym\u00e9trique s&rsquo;explique intuitivement par le fait que l&rsquo;estimateur X<sub>(n)<\/sub> se trouve structurellement bloqu\u00e9 sous la fronti\u00e8re infranchissable b : l&rsquo;erreur (b &#8211; X<sub>(n)<\/sub>) ne peut prendre que des valeurs strictement positives. Cette loi limite permet d&rsquo;\u00e9laborer des intervalles de confiance asymptotiques exacts sans faire appel \u00e0 l&rsquo;approximation gaussienne : pour un seuil de confiance 1 &#8211; &alpha;, un intervalle de confiance pivot exact pour b s&rsquo;exprime par :<\/p>\n<p>[ X<sub>(n)<\/sub> , X<sub>(n)<\/sub> &#8211; (b&#770; &#8211; &acirc;) ln(&alpha;) \/ n ]<\/p>\n<p>garantissant une couverture probabiliste rigoureuse m\u00eame sous de petits effectifs.<\/p>\n<h2>10. Comparaison m\u00e9thodologique approfondie : EMV contre M\u00e9thode des Moments<\/h2>\n<h3>10.1 D\u00e9rivation des estimateurs par la m\u00e9thode des moments (MM)<\/h3>\n<p>Pour mesurer la qualit\u00e9 de l&rsquo;estimateur du maximum de vraisemblance, il convient de le confronter \u00e0 la m\u00e9thode d&rsquo;inf\u00e9rence concurrente historique : la m\u00e9thode des moments formalis\u00e9e par Karl Pearson. Cette technique repose sur le principe d&rsquo;\u00e9quivalence consistant \u00e0 \u00e9galiser les moments th\u00e9oriques de la population parente avec les moments empiriques calcul\u00e9s directement sur l&rsquo;\u00e9chantillon d&rsquo;observations.<\/p>\n<p>Pour une loi uniforme continue U(a, b), les deux premiers moments th\u00e9oriques sont l&rsquo;esp\u00e9rance &mu; = (a + b) \/ 2 et la variance &sigma;<sup>2<\/sup> = (b &#8211; a)<sup>2<\/sup> \/ 12. D\u00e9finissons la moyenne empirique X&#772;<sub>n<\/sub> et la variance empirique d&rsquo;\u00e9chantillon S<sub>n<\/sub><sup>2<\/sup> :<\/p>\n<p>X&#772;<sub>n<\/sub> = (1 \/ n) &sum;<sub>i=1<\/sub><sup>n<\/sup> X<sub>i<\/sub> &emsp; et &emsp; S<sub>n<\/sub><sup>2<\/sup> = (1 \/ n) &sum;<sub>i=1<\/sub><sup>n<\/sup> (X<sub>i<\/sub> &#8211; X&#772;<sub>n<\/sub>)<sup>2<\/sup><\/p>\n<p>Le syst\u00e8me d&rsquo;\u00e9quations des moments s&rsquo;\u00e9crit formellement :<\/p>\n<p>{ (a + b) \/ 2 = X&#772;<sub>n<\/sub> ; (b &#8211; a)<sup>2<\/sup> \/ 12 = S<sub>n<\/sub><sup>2<\/sup> }<\/p>\n<p>En extrayant la racine carr\u00e9e positive de la seconde \u00e9quation, nous obtenons l&rsquo;estimation de l&rsquo;amplitude du support : b &#8211; a = &radic;(12 S<sub>n<\/sub><sup>2<\/sup>) = 2 &radic;3 S<sub>n<\/sub>. En combinant cette relation avec l&rsquo;\u00e9quation de la moyenne empirique par substitution et addition lin\u00e9aire, nous d\u00e9rivons de fa\u00e7on \u00e9l\u00e9mentaire les estimateurs de la m\u00e9thode des moments :<\/p>\n<p>&acirc;<sub>MM<\/sub> = X&#772;<sub>n<\/sub> &#8211; &radic;3 S<sub>n<\/sub><\/p>\n<p>b&#770;<sub>MM<\/sub> = X&#772;<sub>n<\/sub> + &radic;3 S<sub>n<\/sub><\/p>\n<p>Contrairement aux estimateurs EMV qui reposent exclusivement sur les valeurs fronti\u00e8res extr\u00eames, les estimateurs de la m\u00e9thode des moments sont des fonctions diff\u00e9rentiables globales int\u00e9grant chaque r\u00e9alisation d&rsquo;\u00e9chantillon via des sommes alg\u00e9briques r\u00e9guli\u00e8res.<\/p>\n<h3>10.2 Comparaison analytique des variances et efficacit\u00e9 relative<\/h3>\n<p>Puisque la moyenne empirique X&#772;<sub>n<\/sub> et la variance empirique S<sub>n<\/sub><sup>2<\/sup> sont assujetties au th\u00e9or\u00e8me central limite usuel de Lindeberg-L\u00e9vy, les estimateurs issus de la m\u00e9thode des moments convergent vers la cible th\u00e9orique \u00e0 la vitesse param\u00e9trique conventionnelle &radic;n. En effet, la variance asymptotique de l&rsquo;estimateur b&#770;<sub>MM<\/sub> se comporte selon l&rsquo;ordre :<\/p>\n<p>Var(b&#770;<sub>MM<\/sub>) &sim; O( 1 \/ n )<\/p>\n<p>Rappelons que pour l&rsquo;estimateur sans biais corrig\u00e9 issu de la d\u00e9marche de vraisemblance, la variance s&rsquo;exprime selon l&rsquo;ordre :<\/p>\n<p>Var(b&#770;<sub>UMVU<\/sub>) = &theta;<sup>2<\/sup> \/ [ n (n + 2) ] &sim; O( 1 \/ n<sup>2<\/sup> )<\/p>\n<p>L&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Efficacit%C3%A9_(statistique)\">efficacit\u00e9 relative asymptotique (ERA)<\/a> de la m\u00e9thode des moments par rapport \u00e0 l&rsquo;approche du maximum de vraisemblance s&rsquo;\u00e9value en calculant la limite du ratio de leurs variances respectives lorsque la taille d&rsquo;\u00e9chantillon n tend vers l&rsquo;infini :<\/p>\n<p>ERA( b&#770;<sub>MM<\/sub> , b&#770;<sub>UMVU<\/sub> ) = lim<sub>n&rarr;&infin;<\/sub> [ Var(b&#770;<sub>UMVU<\/sub>) \/ Var(b&#770;<sub>MM<\/sub>) ] = lim<sub>n&rarr;&infin;<\/sub> [ O(1\/n<sup>2<\/sup>) \/ O(1\/n) ] = lim<sub>n&rarr;&infin;<\/sub> O( 1 \/ n ) = 0<\/p>\n<p>Une efficacit\u00e9 relative nulle prouve la sup\u00e9riorit\u00e9 asymptotique \u00e9crasante de l&rsquo;estimateur du maximum de vraisemblance sur l&rsquo;estimateur des moments. En grande dimension, la m\u00e9thode des moments est infiniment moins pr\u00e9cise que l&rsquo;EMV.<\/p>\n<p>Au-del\u00e0 de cette inf\u00e9riorit\u00e9 de variance, la m\u00e9thode des moments souffre d&rsquo;un d\u00e9faut conceptuel r\u00e9dhibitoire pour l&rsquo;analyste praticien : la production fr\u00e9quente d&rsquo;estimations absurdes ou <em>inadmissibles<\/em> au regard des donn\u00e9es recueillies. En effet, rien dans la formulation analytique de b&#770;<sub>MM<\/sub> = X&#772;<sub>n<\/sub> + &radic;3 S<sub>n<\/sub> ne contraint math\u00e9matiquement cette quantit\u00e9 \u00e0 \u00eatre sup\u00e9rieure \u00e0 la plus grande observation r\u00e9alis\u00e9e X<sub>(n)<\/sub>. Dans un nombre non n\u00e9gligeable de configurations d&rsquo;\u00e9chantillonnage fini, il arrive empiriquement que b&#770;<sub>MM<\/sub> &lt; X<sub>(n)<\/sub> ou que &acirc;<sub>MM<\/sub> &gt; X<sub>(1)<\/sub>. Dans de telles circonstances, la m\u00e9thode des moments produit un mod\u00e8le selon lequel certaines donn\u00e9es r\u00e9elles d\u00e9j\u00e0 observ\u00e9es auraient eu une probabilit\u00e9 th\u00e9orique rigoureusement nulle d&rsquo;appara\u00eetre, ce qui est une aberration statistique que l&rsquo;EMV ne commet jamais par construction.<\/p>\n<h3>10.3 Le probl\u00e8me historique des chars allemands (German Tank Problem)<\/h3>\n<p>L&rsquo;application historique la plus c\u00e9l\u00e8bre de cette distinction m\u00e9thodologique fondamentale entre moments et statistiques d&rsquo;ordre extr\u00eames est sans conteste le <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Probl%C3%A8me_du_char_allemand\">probl\u00e8me des chars allemands<\/a> (<em>German Tank Problem<\/em>) survenu au cours de la Seconde Guerre mondiale. Les forces de renseignement alli\u00e9es cherchaient \u00e0 \u00e9valuer avec une fiabilit\u00e9 critique le rythme de production mensuel des blind\u00e9s Panther et Tigre de l&rsquo;Allemagne nazie.<\/p>\n<p>L&rsquo;appareil industriel allemand appliquait une logique bureaucratique rigoureuse de tra\u00e7abilit\u00e9 en estampillant les bo\u00eetes de vitesses et les ch\u00e2ssis des chars d&rsquo;un num\u00e9ro de s\u00e9rie s\u00e9quentiel continu allant de 1 jusqu&rsquo;\u00e0 la quantit\u00e9 totale produite N. Les pi\u00e8ces captur\u00e9es sur le champ de bataille par les troupes alli\u00e9es constituaient un \u00e9chantillon al\u00e9atoire discret sans remise tir\u00e9 uniform\u00e9ment sur l&rsquo;ensemble d&rsquo;entiers {1, 2, &#8230;, N}. Le d\u00e9fi math\u00e9matique consistait \u00e0 estimer le param\u00e8tre inconnu N \u00e0 partir d&rsquo;un \u00e9chantillon fini de num\u00e9ros de s\u00e9rie observ\u00e9s de taille n, dont le maximum \u00e9tait k = X<sub>(n)<\/sub>.<\/p>\n<p>Les approches conventionnelles de renseignement humain et d&rsquo;espionnage estimaient que la production allemande s&rsquo;\u00e9levait \u00e0 environ 1 400 chars par mois. Les statisticiens alli\u00e9s, mobilisant l&rsquo;\u00e9quivalent discret de la th\u00e9orie du maximum de vraisemblance et la correction sans biais UMVU d\u00e9velopp\u00e9e dans ce m\u00e9moire, ont appliqu\u00e9 la formule de l&rsquo;estimateur non biais\u00e9 \u00e0 variance minimale :<\/p>\n<p>N&#770; = X<sub>(n)<\/sub> + [ (X<sub>(n)<\/sub> &#8211; n) \/ n ] = X<sub>(n)<\/sub> [ (n + 1) \/ n ] &#8211; 1<\/p>\n<p>Gr\u00e2ce \u00e0 cette formule ne mobilisant que la statistique d&rsquo;ordre maximale corrig\u00e9e par le ratio d&rsquo;\u00e9chantillon, les analystes statistiques conclurent que la production industrielle n&rsquo;\u00e9tait que de 246 chars par mois entre 1940 et 1942. \u00c0 la lib\u00e9ration, l&rsquo;inspection des registres de production du minist\u00e8re de l&rsquo;Armement d&rsquo;Albert Speer r\u00e9v\u00e9la que le volume moyen exact s&rsquo;\u00e9levait \u00e0 245 blind\u00e9s par mois. L&rsquo;inf\u00e9rence par statistique d&rsquo;ordre s&rsquo;est montr\u00e9e immens\u00e9ment plus exacte que les estimations des services secrets, validant empiriquement la th\u00e9orie statistique sur un enjeu g\u00e9opolitique majeur.<\/p>\n<h2>11. Applications en psychom\u00e9trie, psychophysique et sciences du comportement<\/h2>\n<h3>11.1 Mod\u00e9lisation des seuils de perception et temps de r\u00e9action planchers<\/h3>\n<p>En <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Sciences_cognitives\">sciences cognitives<\/a> et en psychophysique exp\u00e9rimentale, l&rsquo;estimation des fronti\u00e8res distributionnelles trouve un champ d&rsquo;application naturel dans l&rsquo;analyse chronom\u00e9trique fine des temps de r\u00e9action motrice face \u00e0 un stimulus perceptif. La distribution des temps de r\u00e9ponse d&rsquo;un sujet humain ne commence jamais \u00e0 t = 0 ; elle est soumise \u00e0 une latence physiologique incompressible t<sub>0<\/sub>, correspondant \u00e0 la conduction neuronale du signal sensoriel aff\u00e9rent depuis le r\u00e9cepteur sensoriel (r\u00e9tine ou cochl\u00e9e) jusqu&rsquo;au cortex moteur eff\u00e9rent.<\/p>\n<p>Les mod\u00e9lisations formelles de la prise de d\u00e9cision, telles que les mod\u00e8les de diffusion de drift ou les mod\u00e8les de course cognitive (<em>race models<\/em>), incorporent syst\u00e9matiquement ce param\u00e8tre de fronti\u00e8re t<sub>0<\/sub> sous forme d&rsquo;un d\u00e9calage temporel pur. Dans diverses paradigmes de d\u00e9tection rapide o\u00f9 la phase d\u00e9cisionnelle pr\u00e9sente une incertitude maximale d&rsquo;ordre temporel, la phase pr\u00e9-d\u00e9cisionnelle est mod\u00e9lis\u00e9e par une distribution uniforme par morceaux U(t<sub>0<\/sub>, t<sub>0<\/sub> + &Delta;). L&rsquo;estimation de ce seuil de latence motrice absolue t<sub>0<\/sub> mobilise directement l&rsquo;estimateur de la statistique d&rsquo;ordre minimale X<sub>(1)<\/sub> ajust\u00e9 par la correction de biais UMVU :<\/p>\n<p>t&#770;<sub>0<\/sub> = X<sub>(1)<\/sub> &#8211; [ (X<sub>(n)<\/sub> &#8211; X<sub>(1)<\/sub>) \/ (n &#8211; 1) ]<\/p>\n<p>L&rsquo;omission de cette correction structurelle conduirait le neuropsychologue \u00e0 surestimer le plancher physiologique d&rsquo;excitabilit\u00e9 synaptique du participant, biaillant l&rsquo;\u00e9valuation des processus cognitifs d&rsquo;attention s\u00e9lective.<\/p>\n<h3>11.2 Analyse des \u00e9chelles de r\u00e9ponse born\u00e9es en psychom\u00e9trie cognitive<\/h3>\n<p>La psychom\u00e9trie et l&rsquo;\u00e9valuation standardis\u00e9e des aptitudes cognitives sont fr\u00e9quemment confront\u00e9es \u00e0 des art\u00e9facts m\u00e9thodologiques majeurs induits par la saturation des instruments de mesure, couramment document\u00e9s sous les vocables d&rsquo;effets de plafond et d&rsquo;effets de plancher. Lorsqu&rsquo;une batterie de tests \u00e9valuant la capacit\u00e9 maximale de la m\u00e9moire de travail (telle que l&#8217;empan mn\u00e9sique visuo-spatial) est soumise \u00e0 des cohortes de haut niveau acad\u00e9mique, un nombre substantiel de scores s&rsquo;agglutine au niveau de la note maximale th\u00e9orique de l&rsquo;\u00e9preuve.<\/p>\n<p>Plut\u00f4t que de pr\u00e9sumer abusivement une normalit\u00e9 latente gaussienne via des transformations non lin\u00e9aires arbitraires de type Box-Cox, les m\u00e9thodologistes recourent de plus en plus \u00e0 des lois uniformes discr\u00e8tes ou continues tronqu\u00e9es pour mod\u00e9liser le comportement de la variance r\u00e9siduelle sur des tranches d&rsquo;aptitude homog\u00e8nes. L&rsquo;estimation par maximum de vraisemblance de la capacit\u00e9 maximale latente &theta; permet alors d&rsquo;identifier le plafond effectif du test sans biaiser la variance de l&rsquo;\u00e9chantillon par des distributions sym\u00e9triques \u00e0 dispersion non born\u00e9e.<\/p>\n<h3>11.3 \u00c9tude empirique de l&rsquo;estimation subjective d&rsquo;intervalles chez l&rsquo;humain<\/h3>\n<p>Une dimension heuristique fascinante de la psychologie du jugement et de la prise de d\u00e9cision r\u00e9side dans l&rsquo;analyse de la mani\u00e8re dont le cerveau humain traite intuitivement les probl\u00e8mes de fronti\u00e8re stochastique. Des exp\u00e9riences men\u00e9es en psychologie cognitive exp\u00e9rimentale d\u00e9montrent que les sujets humains confront\u00e9s \u00e0 une s\u00e9rie d&rsquo;\u00e9chantillons continus sont syst\u00e9matiquement incapables de corriger spontan\u00e9ment le biais de sous-estimation de la statistique d&rsquo;ordre maximale.<\/p>\n<p>Lorsqu&rsquo;on demande \u00e0 des participants de pr\u00e9dire la valeur maximale d&rsquo;un processus dont ils ont observ\u00e9 n r\u00e9alisations, l&rsquo;intuition humaine retient presque invariablement la plus grande valeur observ\u00e9e x<sub>(n)<\/sub> comme r\u00e9f\u00e9rence absolue sans int\u00e9grer la correction math\u00e9matique n\u00e9cessaire [(n + 1)\/n]. Cette distorsion cognitive correspond \u00e0 un ancrage mn\u00e9sique rigide sur les extr\u00eames observ\u00e9s, expliquant la faillibilit\u00e9 des pr\u00e9visions humaines face aux risques environnementaux extr\u00eames (crues centennales, s\u00e9ismes, fluctuations boursi\u00e8res maximales) o\u00f9 les bornes de variabilit\u00e9 physique sont syst\u00e9matiquement sous-\u00e9valu\u00e9es par les agents non avertis.<\/p>\n<h2>12. Impl\u00e9mentation computationnelle, simulations de Monte-Carlo et robustesse<\/h2>\n<h3>12.1 Algorithmes d&rsquo;estimation et simulations en Python et R<\/h3>\n<p>Sur le plan algorithmique et computationnel, l&rsquo;ex\u00e9cution de l&rsquo;estimateur du maximum de vraisemblance et de ses variantes corrig\u00e9es pour la loi uniforme se distingue par une complexit\u00e9 algorithmique lin\u00e9aire optimale O(n), comparativement \u00e0 la plupart des proc\u00e9dures d&rsquo;optimisation non lin\u00e9aire requ\u00e9rant des sch\u00e9mas it\u00e9ratifs lourds.<\/p>\n<p>Le calcul se d\u00e9compose en un simple balayage s\u00e9quentiel du vecteur de donn\u00e9es pour extraire les extrema absolus min(X) et max(X). En langage Python (mobilisant la biblioth\u00e8que vectoris\u00e9e NumPy) comme en environnement statistique R, l&rsquo;op\u00e9ration s&rsquo;ex\u00e9cute nativement en exploitant des primitives compil\u00e9es en C hautement optimis\u00e9es. Les estimateurs sont formellement impl\u00e9ment\u00e9s par les relations vectorielles directes :<\/p>\n<p>Pour le mod\u00e8le bilat\u00e9ral :<\/p>\n<ul>\n<li>EMV : <code>a_emv = min(x)<\/code> et <code>b_emv = max(x)<\/code><\/li>\n<li>UMVU : <code>a_umvu = (n * min(x) - max(x)) \/ (n - 1)<\/code> et <code>b_umvu = (n * max(x) - min(x)) \/ (n - 1)<\/code><\/li>\n<\/ul>\n<p>Pour le mod\u00e8le unilat\u00e9ral U(0, &theta;) :<\/p>\n<ul>\n<li>EMV : <code>theta_emv = max(x)<\/code><\/li>\n<li>UMVU : <code>theta_umvu = max(x) * (n + 1) \/ n<\/code><\/li>\n<\/ul>\n<p>Une \u00e9tude de simulation computationnelle de <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/M%C3%A9thode_de_Monte-Carlo\">Monte-Carlo<\/a> men\u00e9e sur B = 100 000 it\u00e9rations ind\u00e9pendantes pour diff\u00e9rentes tailles d&rsquo;\u00e9chantillon (de n = 5 \u00e0 n = 1 000) permet d&rsquo;observer empiriquement la d\u00e9croissance de l&rsquo;erreur quadratique moyenne. Alors que le trac\u00e9 logarithmique de l&rsquo;EQM de la m\u00e9thode des moments en fonction de ln(n) r\u00e9v\u00e8le une droite de r\u00e9gression de pente standard -1 (traduisant une vitesse en 1\/n), la courbe d&rsquo;EQM des estimateurs EMV et UMVU pr\u00e9sente une pente raide de -2, confirmant la pr\u00e9diction th\u00e9orique du taux d&rsquo;efficacit\u00e9 asymptotique super-convergent en 1\/n<sup>2<\/sup>.<\/p>\n<h3>12.2 Sensibilit\u00e9 critique aux valeurs aberrantes (outliers) et contamination<\/h3>\n<p>En contrepartie de sa super-efficacit\u00e9 analytique remarquable sous des conditions de mod\u00e9lisation id\u00e9ales, l&rsquo;estimateur du maximum de vraisemblance de la distribution uniforme continue souffre d&rsquo;une fragilit\u00e9 structurelle absolue sur le plan de la <a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Statistique_robuste\">statistique robuste<\/a> : son point de rupture (<em>breakdown point<\/em>) asymptotique est strictement minimal et s&rsquo;\u00e9tablit \u00e0 :<\/p>\n<p>&epsilon;* = 1 \/ n &rarr; 0 lorsque n &rarr; &infin;<\/p>\n<p>Le point de rupture repr\u00e9sente la proportion minimale de donn\u00e9es corrompues ou aberrantes qu&rsquo;un \u00e9chantillon peut contenir avant qu&rsquo;un estimateur ne soit projet\u00e9 vers des valeurs arbitrairement fausses. Puisque b&#770;<sub>EMV<\/sub> d\u00e9pend exclusivement de la statistique d&rsquo;ordre maximale X<sub>(n)<\/sub>, il suffit qu&rsquo;une unique observation parmi les n donn\u00e9es recueillies soit affect\u00e9e par une erreur de mesure, un artefact d&rsquo;enregistrement ou une valeur extr\u00eame aberrante (par exemple x<sub>i<\/sub> = 1 000 au lieu de 1) pour que l&rsquo;estimation de la borne sup\u00e9rieure soit instantan\u00e9ment contamin\u00e9e et propuls\u00e9e \u00e0 1 000, d\u00e9truisant instantan\u00e9ment l&rsquo;int\u00e9gralit\u00e9 de la validit\u00e9 inf\u00e9rentielle.<\/p>\n<p>Pour pallier cette hypersensibilit\u00e9 catastrophique dans les environnements de donn\u00e9es bruit\u00e9es, les th\u00e9oriciens de la robustesse pr\u00e9conisent le recours \u00e0 des estimateurs alternatifs :<\/p>\n<ul>\n<li><strong>Statistiques d&rsquo;ordre \u00e9lagu\u00e9es (trimmed order statistics) :<\/strong> Remplacement de la valeur extr\u00eame brute par un quantile sup\u00e9rieur \u00e9lev\u00e9 mais stable, tel que X<sub>(n-k)<\/sub> avec k &ge; 1, associ\u00e9 \u00e0 une correction d&rsquo;\u00e9chelle analytique appropri\u00e9e.<\/li>\n<li><strong>Estimateurs M \u00e0 support compact :<\/strong> Utilisation de fonctions de perte born\u00e9es (telles que la fonction de perte biweight de Tukey ou la fonction de Huber) qui p\u00e9nalisent et tronquent l&rsquo;influence des points excessivement \u00e9loign\u00e9s du barycentre du nuage de donn\u00e9es.<\/li>\n<li><strong>Mod\u00e8les de m\u00e9lange uniforme-gaussien :<\/strong> Mod\u00e9lisation de l&rsquo;\u00e9chantillon sous la forme d&rsquo;un m\u00e9lange probabiliste associant une composante uniforme U(a, b) avec une probabilit\u00e9 1 &#8211; &gamma;, et une composante de bruit \u00e0 queue lourde (ex. loi de Cauchy ou loi de Laplace) avec une probabilit\u00e9 de contamination &gamma;, maximis\u00e9 via l&rsquo;<a rel=\"noopener noreferrer\" href=\"https:\/\/fr.wikipedia.org\/wiki\/Algorithme_esp%C3%A9rance-maximisation\">algorithme EM<\/a> (Esp\u00e9rance-Maximisation).<\/li>\n<\/ul>\n<h3>12.3 Synth\u00e8se g\u00e9n\u00e9rale et recommandations pratiques pour l&rsquo;analyste<\/h3>\n<p>L&rsquo;analyse approfondie de l&rsquo;estimation par maximum de vraisemblance pour la distribution uniforme continue permet de d\u00e9gager des directives claires et hi\u00e9rarchis\u00e9es pour les praticiens de la science des donn\u00e9es et les mod\u00e9lisateurs confront\u00e9s \u00e0 des supports compacts. Le choix de l&rsquo;estimateur ne doit jamais \u00eatre automatis\u00e9, mais doit r\u00e9sulter d&rsquo;un arbitrage rigoureux entre la puret\u00e9 de la distribution th\u00e9orique, le co\u00fbt du biais d&rsquo;\u00e9chantillonnage fini et la vuln\u00e9rabilit\u00e9 au bruitage des mesures.<\/p>\n<p>Le tableau comparatif suivant synth\u00e9tise les propri\u00e9t\u00e9s op\u00e9rationnelles des diff\u00e9rentes m\u00e9thodologies d&rsquo;estimation explor\u00e9es dans cette monographie :<\/p>\n<ul>\n<li><strong>Estimateur EMV Brut [min(X), max(X)] :<\/strong> Vitesse de convergence ultra-rapide en O(1\/n<sup>2<\/sup>) pour l&rsquo;EQM ; calcul computationnel instantan\u00e9 ; pr\u00e9sence d&rsquo;un biais syst\u00e9matique strict (sous-estimation de l&rsquo;\u00e9tendue) ; inad\u00e9quat pour les tr\u00e8s petits \u00e9chantillons (n &lt; 15) en inf\u00e9rence exacte sans correction ; point de rupture nul (1\/n).<\/li>\n<li><strong>Estimateur UMVU Sans Biais :<\/strong> Vitesse de convergence ultra-rapide en O(1\/n<sup>2<\/sup>) ; absence totale de biais th\u00e9orique pour toute taille n ; EQM strictement inf\u00e9rieure \u00e0 l&rsquo;EMV brut pour tout n &ge; 2 ; repr\u00e9sente le choix optimal absolu d\u00e8s lors que l&rsquo;absence de valeurs aberrantes au-del\u00e0 des bornes est physiquement garantie.<\/li>\n<li><strong>Estimateur de la M\u00e9thode des Moments (MM) :<\/strong> Vitesse de convergence lente en O(1\/n) pour l&rsquo;EQM ; asymptotiquement inefficace face \u00e0 l&rsquo;EMV ; risque d&rsquo;inadmissibilit\u00e9 structurelle en produisant des fronti\u00e8res internes aux donn\u00e9es ; cependant, il offre une robustesse l\u00e9g\u00e8rement sup\u00e9rieure face \u00e0 des bruits localis\u00e9s aux bords en moyennant sur l&rsquo;ensemble de l&rsquo;\u00e9chantillon.<\/li>\n<li><strong>Approches Bay\u00e9siennes R\u00e9gularis\u00e9es :<\/strong> Int\u00e9gration d&rsquo;une loi a priori (loi de Pareto ou loi Gamma inverse sur la largeur b &#8211; a) ; production d&rsquo;une distribution a posteriori permettant la construction directe d&rsquo;intervalles de cr\u00e9dibilit\u00e9 d&rsquo;\u00e9chantillonnage fini coh\u00e9rents sans hypoth\u00e8se asymptotique.<\/li>\n<\/ul>\n<p>En conclusion de cette monographie, l&rsquo;estimation par le maximum de vraisemblance de la loi uniforme continue demeure l&rsquo;un des exemples les plus \u00e9clatants de la th\u00e9orie statistique. Elle rappelle aux chercheurs que l&rsquo;inf\u00e9rence ne saurait \u00eatre r\u00e9duite \u00e0 l&rsquo;application m\u00e9canique du calcul diff\u00e9rentiel matriciel. Lorsque la fronti\u00e8re des \u00e9v\u00e9nements admissibles se d\u00e9place sous l&rsquo;effet des param\u00e8tres m\u00eames de l&rsquo;univers que l&rsquo;on observe, la statistique d&rsquo;ordre pure supplante le gradient, r\u00e9v\u00e9lant la puissance g\u00e9om\u00e9trique log\u00e9e aux marges extr\u00eames de l&rsquo;observation scientifique.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Cram\u00e9r, H. (1946). <em>Mathematical Methods of Statistics<\/em>. Princeton University Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/press.princeton.edu\/books\/paperback\/9780691080048\/mathematical-methods-of-statistics\">https:\/\/press.princeton.edu\/books\/paperback\/9780691080048\/mathematical-methods-of-statistics<\/a><\/li>\n<li>David, H. A., &amp; Nagaraja, H. N. (2003). <em>Order Statistics<\/em> (3rd ed.). John Wiley &amp; Sons. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471722162\">https:\/\/doi.org\/10.1002\/0471722162<\/a><\/li>\n<li>Fisher, R. A. (1922). On the mathematical foundations of theoretical statistics. <em>Philosophical Transactions of the Royal Society of London. Series A<\/em>, 222(594-604), 309-368. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1098\/rsta.1922.0009\">https:\/\/doi.org\/10.1098\/rsta.1922.0009<\/a><\/li>\n<li>Hogg, R. V., McKean, J., &amp; Craig, A. T. (2018). <em>Introduction to Mathematical Statistics<\/em> (8th ed.). Pearson. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.pearson.com\/en-us\/subject-catalog\/p\/introduction-to-mathematical-statistics\/P200000003332\">https:\/\/www.pearson.com\/en-us\/subject-catalog\/p\/introduction-to-mathematical-statistics\/P200000003332<\/a><\/li>\n<li>Huber, P. J., &amp; Ronchetti, E. M. (2009). <em>Robust Statistics<\/em> (2nd ed.). John Wiley &amp; Sons. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/9780470434697\">https:\/\/doi.org\/10.1002\/9780470434697<\/a><\/li>\n<li>Lehmann, E. L., &amp; Casella, G. (1998). <em>Theory of Point Estimation<\/em> (2nd ed.). Springer-Verlag. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/b98854\">https:\/\/doi.org\/10.1007\/b98854<\/a><\/li>\n<li>Ruggles, R., &amp; Brodie, H. (1947). An empirical approach to economic intelligence in World War II. <em>Journal of the American Statistical Association<\/em>, 42(237), 72-91. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1947.10501915\">https:\/\/doi.org\/10.1080\/01621459.1947.10501915<\/a><\/li>\n<li>van der Vaart, A. W. (1998). <em>Asymptotic Statistics<\/em>. Cambridge University Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1017\/CBO9780511802256\">https:\/\/doi.org\/10.1017\/CBO9780511802256<\/a><\/li>\n<\/ul>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_.png\" alt=\"likelihood function for uniform distribution\" \/><\/figure>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_-1.png\" alt=\"Log-likelihood function of uniform distribution\" \/><\/figure>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_-2.png\" alt=\"Derivative of log-likelihood function for uniform distribution\" \/><\/figure>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/latex.php_-3.png\" alt=\"Partial derivative of log-likelihood function for uniform distribution\" \/><\/figure>\n","protected":false},"excerpt":{"rendered":"<p>Guide acad\u00e9mique complet sur l&rsquo;estimation par le maximum de vraisemblance pour une loi uniforme continue : th\u00e9orie, calculs, propri\u00e9t\u00e9s et applications.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2891,8670],"tags":[5132,8671,8673,16,8672],"class_list":["post-7192","psych_stat","type-psych_stat","status-publish","hentry","category-methodes-quantitatives-en-psychologie","category-statistique-mathematique","tag-distribution-uniforme","tag-estimation-par-maximum-de-vraisemblance","tag-proprietes-asymptotiques","tag-psychometrie","tag-statistiques-dordre"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7192","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=7192"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7192\/revisions"}],"predecessor-version":[{"id":7202,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7192\/revisions\/7202"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=7192"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=7192"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=7192"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}