{"id":7220,"date":"2026-09-16T18:23:27","date_gmt":"2026-09-16T18:23:27","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/qu-est-ce-que-la-distribution-d-erlang\/"},"modified":"2026-09-16T18:23:29","modified_gmt":"2026-09-16T18:23:29","slug":"qu-est-ce-que-la-distribution-d-erlang","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=qu-est-ce-que-la-distribution-d-erlang","title":{"rendered":"Qu&rsquo;est-ce que la distribution d&rsquo;Erlang ?"},"content":{"rendered":"<p>Au c\u0153ur de la th\u00e9orie des probabilit\u00e9s appliqu\u00e9es et de la mod\u00e9lisation des processus stochastiques, la <strong>distribution d&rsquo;Erlang<\/strong> occupe une place \u00e9minente et singuli\u00e8re. Con\u00e7ue initialement pour r\u00e9pondre aux d\u00e9fis techniques et logistiques pos\u00e9s par l&rsquo;essor des r\u00e9seaux t\u00e9l\u00e9phoniques au d\u00e9but du XXe si\u00e8cle, cette loi de probabilit\u00e9 continue s&rsquo;est progressivement \u00e9mancip\u00e9e de son champ d&rsquo;application originel pour devenir un instrument analytique universel. En reliant de mani\u00e8re \u00e9l\u00e9gante les lois continues et les processus de comptage discrets, elle offre un cadre rigoureux pour formaliser le temps requis pour accomplir une suite ordonn\u00e9e d&rsquo;\u00e9tapes ind\u00e9pendantes, chacune r\u00e9gie par un m\u00e9canisme exponentiel sans m\u00e9moire.<\/p>\n<p>D&rsquo;un point de vue structural, la loi d&rsquo;Erlang repr\u00e9sente un pont fondamental entre la distribution exponentielle \u2014 dont elle constitue la somme de variables ind\u00e9pendantes et identiquement distribu\u00e9es \u2014 et la distribution Gamma, dont elle est la restriction exacte aux param\u00e8tres de forme entiers. Cette sp\u00e9cificit\u00e9 alg\u00e9brique conf\u00e8re \u00e0 la loi d&rsquo;Erlang une tractabilit\u00e9 analytique remarquable : l\u00e0 o\u00f9 les distributions continues arbitraires imposent des approximations num\u00e9riques complexes, la distribution d&rsquo;Erlang permet des r\u00e9solutions explicites et exactes gr\u00e2ce \u00e0 des sommes finies de termes poissoniens. De ce fait, elle demeure une pierre angulaire dans la mod\u00e9lisation des ph\u00e9nom\u00e8nes de flux, d&rsquo;attente et de d\u00e9gradation.<\/p>\n<p>Le pr\u00e9sent article propose une exploration exhaustive et approfondie de la distribution d&rsquo;Erlang, depuis sa gen\u00e8se historique au sein des centraux t\u00e9l\u00e9phoniques de Copenhague jusqu&rsquo;\u00e0 ses utilisations les plus contemporaines en biologie math\u00e9matique, en informatique distribu\u00e9e et en sciences cognitives. \u00c0 travers une analyse d\u00e9taill\u00e9e de ses fondements analytiques, de ses propri\u00e9t\u00e9s statistiques, de ses dynamiques de transition et de ses protocoles d&rsquo;estimation, cette synth\u00e8se acad\u00e9mique vise \u00e0 offrir une compr\u00e9hension globale de ce mod\u00e8le probabiliste indispensable \u00e0 l&rsquo;ing\u00e9nieur, au statisticien et au chercheur en sciences appliqu\u00e9es.<\/p>\n<h2>1. Introduction et gen\u00e8se historique de la distribution d&rsquo;Erlang<\/h2>\n<h3>1.1 Les travaux fondateurs d&rsquo;Agner Krarup Erlang<\/h3>\n<p>L&rsquo;\u00e9mergence de la distribution d&rsquo;Erlang est indissociable du parcours intellectuel du math\u00e9maticien et ing\u00e9nieur danois <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Agner_Krarup_Erlang\">Agner Krarup Erlang<\/a> (1878\u20131929). Embauch\u00e9 en 1908 par la <em>Copenhagen Telephone Company<\/em> (KTAS) sous l&rsquo;impulsion de son directeur Frantz Johannsen, Erlang avait pour mission de rationaliser le fonctionnement d&rsquo;un r\u00e9seau alors en pleine explosion technologique. Au tournant du XXe si\u00e8cle, la t\u00e9l\u00e9phonie subissait une mutation radicale : le passage progressif de la commutation enti\u00e8rement manuelle, op\u00e9r\u00e9e par des standardistes, vers des dispositifs \u00e9lectrom\u00e9caniques automatis\u00e9s imposait une gestion rigoureuse des goulets d&rsquo;\u00e9tranglement.<\/p>\n<p>Confront\u00e9 \u00e0 l&rsquo;incapacit\u00e9 des ing\u00e9nieurs de l&rsquo;\u00e9poque \u00e0 quantifier scientifiquement la capacit\u00e9 optimale des faisceaux de c\u00e2bles, Erlang aborda le probl\u00e8me non pas sous un angle purement empirique, mais par le prisme de la physique statistique et du calcul des probabilit\u00e9s. Ses observations m\u00e9ticuleuses au sein des centraux de Copenhague lui permirent de comprendre que les appels t\u00e9l\u00e9phoniques ne se pr\u00e9sentaient pas de mani\u00e8re d\u00e9terministe ou synchrone, mais constituaient un flux d&rsquo;\u00e9v\u00e9nements al\u00e9atoires discontinus. C&rsquo;est dans ce laboratoire industriel qu&rsquo;il jeta les bases de ce qui allait devenir la discipline fondamentale de l&rsquo;ing\u00e9nierie du trafic et de la stochastique appliqu\u00e9e.<\/p>\n<p>En formulant l&rsquo;hypoth\u00e8se audacieuse que les dur\u00e9es de communication et les intervalles inter-arriv\u00e9es pouvaient \u00eatre d\u00e9compos\u00e9s en segments \u00e9l\u00e9mentaires, Erlang a ouvert la voie \u00e0 la formalisation math\u00e9matique de l&rsquo;encombrement des lignes. Ses travaux, publi\u00e9s notamment en 1909 dans son m\u00e9moire inaugural <em>The Theory of Probabilities and Telephone Conversations<\/em>, d\u00e9montr\u00e8rent pour la premi\u00e8re fois comment mod\u00e9liser math\u00e9matiquement les blocages de lignes d&rsquo;abonn\u00e9s et d\u00e9terminer avec une pr\u00e9cision rigoureuse le compromis \u00e9conomique entre le co\u00fbt des infrastructures et la qualit\u00e9 de service ressentie par les usagers.<\/p>\n<h3>1.2 Le probl\u00e8me initial de commutation et du trafic t\u00e9l\u00e9phonique<\/h3>\n<p>Le probl\u00e8me physique auquel Erlang \u00e9tait confront\u00e9 r\u00e9sidait dans l&rsquo;asym\u00e9trie structurelle entre la demande des usagers et la capacit\u00e9 finie des centraux t\u00e9l\u00e9phoniques. D&rsquo;une part, les appels entrants survenaient de fa\u00e7on purement impr\u00e9visible au fil de la journ\u00e9e ; d&rsquo;autre part, la dur\u00e9e d&rsquo;occupation des circuits par les conversations variait de mani\u00e8re sensible d&rsquo;un abonn\u00e9 \u00e0 l&rsquo;autre. Dans un tel syst\u00e8me, si un abonn\u00e9 tentait d&rsquo;\u00e9tablir une liaison alors que l&rsquo;ensemble des commutateurs ou des lignes de transmission disponibles \u00e9tait d\u00e9j\u00e0 sollicit\u00e9, l&rsquo;appel \u00e9tait irr\u00e9m\u00e9diablement rejet\u00e9 ou perdu, g\u00e9n\u00e9rant une friction technique et commerciale majeure.<\/p>\n<p>Le dimensionnement optimal des r\u00e9seaux n\u00e9cessitait une approche novatrice. Construire un r\u00e9seau capable d&rsquo;absorber la demande maximale absolue sans aucun blocage aurait exig\u00e9 des investissements mat\u00e9riels colossaux, \u00e9conomiquement insoutenables pour la compagnie de t\u00e9l\u00e9phone. \u00c0 l&rsquo;inverse, sous-dimensionner les installations entra\u00eenait un taux de saturation inacceptable et un engorgement perp\u00e9tuel du syst\u00e8me. La r\u00e9solution de ce dilemme r\u00e9sidait dans le calcul de la probabilit\u00e9 exacte qu&rsquo;un syst\u00e8me compos\u00e9 d&rsquo;un nombre d\u00e9termin\u00e9 de circuits tombe en \u00e9tat de blocage sous une charge de trafic donn\u00e9e.<\/p>\n<p>Erlang a formalis\u00e9 ce trafic en introduisant une m\u00e9trique adimensionnelle, qui portera plus tard son nom : l&rsquo;erlang, d\u00e9fini comme le produit du taux moyen d&rsquo;arriv\u00e9e des appels par la dur\u00e9e moyenne de maintien de la communication. En traduisant les flux physiques en variables stochastiques, il a op\u00e9r\u00e9 une rupture \u00e9pist\u00e9mologique majeure, substituant aux r\u00e8gles de pouce approximatives et \u00e0 l&#8217;empirisme d&rsquo;atelier une v\u00e9ritable th\u00e9orie probabiliste du trafic des t\u00e9l\u00e9communications, universellement applicable \u00e0 tout syst\u00e8me d&rsquo;acheminement de requ\u00eates.<\/p>\n<h3>1.3 Naissance de la th\u00e9orie des files d&rsquo;attente<\/h3>\n<p>En mod\u00e9lisant math\u00e9matiquement les transitions d&rsquo;\u00e9tat d&rsquo;un r\u00e9seau t\u00e9l\u00e9phonique, Agner Krarup Erlang ne s&rsquo;est pas content\u00e9 de r\u00e9soudre un probl\u00e8me d&rsquo;ing\u00e9nierie locale : il a fond\u00e9 une branche enti\u00e8rement nouvelle des math\u00e9matiques appliqu\u00e9es, la <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Queueing_theory\">th\u00e9orie des files d&rsquo;attente<\/a>. Le c\u0153ur de sa d\u00e9marche consistait \u00e0 formaliser le concept d&rsquo;\u00e9tat d&rsquo;\u00e9quilibre statistique, selon lequel, en r\u00e9gime permanent, le flux de probabilit\u00e9 entrant dans un \u00e9tat donn\u00e9 du syst\u00e8me compense rigoureusement le flux de probabilit\u00e9 qui en sort. Ce principe fondamental de balance globale a permis d&rsquo;\u00e9crire des syst\u00e8mes diff\u00e9rentiels solubles d\u00e9crivant l&rsquo;\u00e9volution temporelle des congestions.<\/p>\n<p>De ces recherches ont \u00e9merg\u00e9 deux \u00e9quations c\u00e9l\u00e8bres qui dominent encore aujourd&rsquo;hui l&rsquo;ing\u00e9nierie des r\u00e9seaux : la formule d&rsquo;Erlang B, d\u00e9di\u00e9e aux syst\u00e8mes sans attente o\u00f9 tout appel rejet\u00e9 est perdu (mod\u00e8le \u00e0 pertes), et la formule d&rsquo;Erlang C, applicable aux syst\u00e8mes o\u00f9 les requ\u00eates exc\u00e9dentaires sont plac\u00e9es dans un registre d&rsquo;attente jusqu&rsquo;\u00e0 la lib\u00e9ration d&rsquo;un serveur (mod\u00e8le \u00e0 retard). La d\u00e9rivation de ces lois exigeait une mod\u00e9lisation fine du temps de service, conduisant Erlang \u00e0 envisager des dur\u00e9es qui n&rsquo;\u00e9taient plus seulement distribu\u00e9es selon une loi exponentielle simple, mais compos\u00e9es de plusieurs phases de service cons\u00e9cutives.<\/p>\n<p>Cette conceptualisation des processus stochastiques par phases a \u00e9tabli les bases th\u00e9oriques de la recherche op\u00e9rationnelle contemporaine. L&rsquo;influence d&rsquo;Erlang s&rsquo;est \u00e9tendue bien au-del\u00e0 des t\u00e9l\u00e9phones danois : ses concepts ont inspir\u00e9 les travaux de math\u00e9maticiens comme David George Kendall, F\u00e9lix Pollaczek et Andre\u00ef Kolmogorov, conf\u00e9rant aux structures stochastiques une port\u00e9e universelle qui innerve aujourd&rsquo;hui l&rsquo;analyse des r\u00e9seaux informatiques, les cha\u00eenes logistiques mondiales et l&rsquo;optimisation des services de sant\u00e9.<\/p>\n<h2>2. D\u00e9finition math\u00e9matique et fonction de densit\u00e9 de probabilit\u00e9<\/h2>\n<h3>2.1 Formulation analytique de la densit\u00e9 de probabilit\u00e9<\/h3>\n<p>La distribution d&rsquo;Erlang est une distribution de probabilit\u00e9 continue d\u00e9finie sur le demi-axe des r\u00e9els positifs. Elle est caract\u00e9ris\u00e9e par deux param\u00e8tres fondamentaux : un param\u00e8tre de forme, traditionnellement not\u00e9 <em>k<\/em>, qui doit imp\u00e9rativement appartenir \u00e0 l&rsquo;ensemble des entiers strictement positifs (<em>k<\/em> &isin; &amp;naturals;<sup>*<\/sup> = {1, 2, 3, &#8230;}), et un param\u00e8tre d&rsquo;intensit\u00e9 ou de taux, not\u00e9 &lambda;, qui est un nombre r\u00e9el strictement sup\u00e9rieur \u00e0 z\u00e9ro (&lambda; &isin; &amp;reals;<sub>+<\/sub><sup>*<\/sup>). Dans certaines formulations d&rsquo;ing\u00e9nierie, on substitue au taux &lambda; un param\u00e8tre d&rsquo;\u00e9chelle not\u00e9 &mu; ou &theta;, avec la relation d&rsquo;inversion directe &mu; = 1 \/ &lambda;.<\/p>\n<p>Pour toute variable al\u00e9atoire continue <em>X<\/em> suivant une loi d&rsquo;Erlang de param\u00e8tres <em>k<\/em> et &lambda;, que l&rsquo;on note de mani\u00e8re canonique <em>X<\/em> ~ Erlang(<em>k<\/em>, &lambda;), la fonction de densit\u00e9 de probabilit\u00e9 (PDF, pour <em>Probability Density Function<\/em>) s&rsquo;\u00e9nonce rigoureusement pour tout <em>x<\/em> &ge; 0 selon l&rsquo;expression analytique suivante :<\/p>\n<p><em>f(x; k, &lambda;) = [ &lambda;<sup>k<\/sup> &middot; x<sup>k-1<\/sup> &middot; e<sup>-&lambda;x<\/sup> ] \/ (k &#8211; 1)!<\/em><\/p>\n<p>Pour toute valeur strictement n\u00e9gative de la variable, c&rsquo;est-\u00e0-dire pour <em>x<\/em> &lt; 0, la fonction de densit\u00e9 est identiquement nulle : <em>f(x; k, &lambda;) = 0<\/em>. Le d\u00e9nominateur de cette fraction fait intervenir la factorielle <em>(k &#8211; 1)!<\/em>, qui agit comme la constante de normalisation indispensable pour garantir que l&rsquo;int\u00e9grale de la densit\u00e9 sur son domaine de d\u00e9finition soit rigoureusement \u00e9gale \u00e0 l&rsquo;unit\u00e9 :<\/p>\n<p>&int;<sub>0<\/sub><sup>&infin;<\/sup> <em>f(x; k, &lambda;) dx = 1<\/em><\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/erlang1-1.png\" alt=\"Erlang distribution\" \/><figcaption>Erlang distribution<\/figcaption><\/figure>\n<h3>2.2 D\u00e9composition structurelle de l&rsquo;expression math\u00e9matique<\/h3>\n<p>Une analyse attentive de la fonction de densit\u00e9 de probabilit\u00e9 d&rsquo;Erlang r\u00e9v\u00e8le qu&rsquo;elle proc\u00e8de de la multiplication synergique de trois composantes distinctes mais compl\u00e9mentaires. La premi\u00e8re composante est le facteur de normalisation et d&rsquo;intensit\u00e9 globale, repr\u00e9sent\u00e9 par la quantit\u00e9 <em>&lambda;<sup>k<\/sup> \/ (k &#8211; 1)!<\/em>. Ce terme scalaire assure la coh\u00e9rence dimensionnelle de la distribution, garantissant que l&rsquo;aire totale d\u00e9limit\u00e9e par la courbe au-dessus de l&rsquo;axe des abscisses demeure invariante et \u00e9gale \u00e0 un, quelles que soient les valeurs attribu\u00e9es \u00e0 <em>k<\/em> et &lambda;.<\/p>\n<p>La deuxi\u00e8me composante est le terme monomial ou polynomial <em>x<sup>k-1<\/sup><\/em>. Ce facteur g\u00e9om\u00e9trique r\u00e9git le comportement de la distribution au voisinage imm\u00e9diat de l&rsquo;origine (lorsque <em>x<\/em> tend vers 0). Lorsque <em>k = 1<\/em>, ce terme se r\u00e9duit \u00e0 <em>x<sup>0<\/sup> = 1<\/em>, conf\u00e9rant \u00e0 la densit\u00e9 sa valeur maximale d\u00e8s l&rsquo;origine. En revanche, d\u00e8s que <em>k<\/em> &ge; 2, le mon\u00f4me s&rsquo;annule en <em>x = 0<\/em> avec un ordre de multiplicit\u00e9 <em>k &#8211; 1<\/em>. Cela conf\u00e8re \u00e0 la courbe une pente nulle \u00e0 l&rsquo;origine d\u00e8s que <em>k<\/em> &ge; 3, traduisant une latence initiale o\u00f9 la probabilit\u00e9 qu&rsquo;un \u00e9v\u00e9nement survienne de mani\u00e8re instantan\u00e9e est infinit\u00e9simale.<\/p>\n<p>La troisi\u00e8me composante est le facteur exponentiel d\u00e9croissant <em>e<sup>-&lambda;x<\/sup><\/em>. Alors que le terme polynomial tire la distribution vers le haut au fur et \u00e0 mesure que <em>x<\/em> augmente, l&rsquo;exponentielle n\u00e9gative impose un amortissement asymptotique puissant d\u00e8s que <em>x<\/em> devient grand. C&rsquo;est pr\u00e9cis\u00e9ment l&rsquo;affrontement entre la croissance polynomiale et l&rsquo;extinction exponentielle qui donne \u00e0 la densit\u00e9 d&rsquo;Erlang son profil unimodal asym\u00e9trique caract\u00e9ristique, culminant en un sommet unique avant de d\u00e9cro\u00eetre doucement vers z\u00e9ro \u00e0 l&rsquo;infini.<\/p>\n<h3>2.3 Interpr\u00e9tation stochastique : somme de variables exponentielles<\/h3>\n<p>L&rsquo;attrait conceptuel majeur de la distribution d&rsquo;Erlang r\u00e9side dans son interpr\u00e9tation probabiliste concr\u00e8te. Soit une suite de <em>k<\/em> variables al\u00e9atoires ind\u00e9pendantes et identiquement distribu\u00e9es (i.i.d.), not\u00e9es <em>Y<sub>1<\/sub>, Y<sub>2<\/sub>, &#8230;, Y<sub>k<\/sub><\/em>, o\u00f9 chaque variable <em>Y<sub>i<\/sub><\/em> suit une distribution exponentielle de param\u00e8tre d&rsquo;intensit\u00e9 &lambda;, ce que l&rsquo;on note <em>Y<sub>i<\/sub> ~ Exp(&lambda;)<\/em>. Si l&rsquo;on d\u00e9finit la variable al\u00e9atoire <em>X<\/em> comme la somme arithm\u00e9tique exacte de ces <em>k<\/em> composantes :<\/p>\n<p><em>X = &sum;<sub>i=1<\/sub><sup>k<\/sup> Y<sub>i<\/sub> = Y<sub>1<\/sub> + Y<sub>2<\/sub> + &#8230; + Y<sub>k<\/sub><\/em><\/p>\n<p>Alors, en vertu des propri\u00e9t\u00e9s de cl\u00f4ture des distributions sous l&rsquo;op\u00e9ration de convolution, la variable al\u00e9atoire <em>X<\/em> suit exactement une distribution d&rsquo;Erlang de param\u00e8tres <em>k<\/em> et &lambda;. La d\u00e9monstration de cette propri\u00e9t\u00e9 s&rsquo;\u00e9tablit par r\u00e9currence au moyen du produit de convolution des densit\u00e9s. Pour deux variables exponentielles ind\u00e9pendantes <em>Y<sub>1<\/sub><\/em> et <em>Y<sub>2<\/sub><\/em> de param\u00e8tre &lambda;, la densit\u00e9 de leur somme <em>S<sub>2<\/sub> = Y<sub>1<\/sub> + Y<sub>2<\/sub><\/em> s&rsquo;\u00e9crit :<\/p>\n<p><em>f<sub>S<sub>2<\/sub><\/sub>(x) = &int;<sub>0<\/sub><sup>x<\/sup> f<sub>Y<sub>1<\/sub><\/sub>(t) &middot; f<sub>Y<sub>2<\/sub><\/sub>(x &#8211; t) dt = &int;<sub>0<\/sub><sup>x<\/sup> &lambda; e<sup>-&lambda;t<\/sup> &middot; &lambda; e<sup>-&lambda;(x &#8211; t)<\/sup> dt = &lambda;<sup>2<\/sup> e<sup>-&lambda;x<\/sup> &int;<sub>0<\/sub><sup>x<\/sup> dt = &lambda;<sup>2<\/sup> x e<sup>-&lambda;x<\/sup><\/em><\/p>\n<p>Cette expression co\u00efncide rigoureusement avec la formule d&rsquo;Erlang pour <em>k = 2<\/em> puisque <em>(2 &#8211; 1)! = 1! = 1<\/em>. En r\u00e9it\u00e9rant cette convolution pas \u00e0 pas jusqu&rsquo;\u00e0 l&rsquo;ordre <em>k<\/em>, l&rsquo;int\u00e9gration successive des puissances de <em>t<\/em> fait \u00e9merger naturellement le terme factoriel <em>(k &#8211; 1)!<\/em> au d\u00e9nominateur. Ce r\u00e9sultat offre une intuition physique imm\u00e9diate : la loi d&rsquo;Erlang mod\u00e9lise le temps d&rsquo;attente total n\u00e9cessaire pour observer l&rsquo;accomplissement s\u00e9quentiel de <em>k<\/em> phases ind\u00e9pendantes et cons\u00e9cutives, chacune pr\u00e9sentant un taux d&rsquo;avancement constant &lambda;.<\/p>\n<h2>3. Param\u00e9trisation : r\u00f4les des param\u00e8tres de forme et d&rsquo;\u00e9chelle<\/h2>\n<h3>3.1 Le param\u00e8tre de forme k et sa contrainte d&rsquo;int\u00e9gralit\u00e9<\/h3>\n<p>Le param\u00e8tre <em>k<\/em>, couramment d\u00e9sign\u00e9 sous le terme de <em>param\u00e8tre de forme<\/em> (ou <em>shape parameter<\/em>), gouverne la morphologie intrins\u00e8que de la courbe de densit\u00e9 et d\u00e9termine le nombre d&rsquo;\u00e9tapes sous-jacentes du processus physique mod\u00e9lis\u00e9. La condition d&rsquo;appartenance de <em>k<\/em> au corps des entiers naturels non nuls constitue la fronti\u00e8re s\u00e9paratrice formelle entre la distribution d&rsquo;Erlang et la distribution Gamma plus g\u00e9n\u00e9rale. Dans le cadre d&rsquo;Erlang, <em>k<\/em> ne saurait \u00eatre fractionnaire ou n\u00e9gatif : il quantifie obligatoirement une collection discr\u00e8te d&rsquo;\u00e9tapes stochastiques identiquement distribu\u00e9es.<\/p>\n<p>L&rsquo;effet g\u00e9om\u00e9trique d&rsquo;une incr\u00e9mentation de la valeur de <em>k<\/em> est spectaculaire sur la morphologie de la distribution. Lorsque <em>k = 1<\/em>, la courbe pr\u00e9sente une allure strictement d\u00e9croissante avec une convexit\u00e9 continue et un pic maximal situ\u00e9 \u00e0 l&rsquo;origine. D\u00e8s que <em>k = 2<\/em>, le profil bascule vers une forme unimodale asym\u00e9trique, s&rsquo;\u00e9levant depuis l&rsquo;origine jusqu&rsquo;\u00e0 atteindre un sommet avant de redescendre. Lorsque <em>k<\/em> cro\u00eet vers des valeurs \u00e9lev\u00e9es (par exemple <em>k = 10, 20, 50<\/em>), la distribution se sym\u00e9trise progressivement, se resserre autour de sa valeur centrale et perd son \u00e9tirement initial pour adopter un profil en cloche tr\u00e8s marqu\u00e9.<\/p>\n<p>D&rsquo;un point de vue syst\u00e9mique, l&rsquo;augmentation de <em>k<\/em> traduit une r\u00e9gularisation du processus global. Plus le nombre d&rsquo;\u00e9tapes n\u00e9cessaires pour atteindre l&rsquo;\u00e9v\u00e9nement final est grand, plus la variance relative du temps global se r\u00e9duit par compensation mutuelle des fluctuations individuelles. Ce ph\u00e9nom\u00e8ne d&rsquo;agr\u00e9gation d&rsquo;\u00e9tapes \u00e9l\u00e9mentaires illustre de mani\u00e8re concr\u00e8te la loi des grands nombres au sein d&rsquo;un m\u00e9canisme de transmission s\u00e9quentiel.<\/p>\n<h3>3.2 Param\u00e8tre d&rsquo;\u00e9chelle (&mu;) versus param\u00e8tre de taux (&lambda;)<\/h3>\n<p>La litt\u00e9rature probabiliste et statistique utilise deux param\u00e9trisations duales pour r\u00e9gir l&rsquo;\u00e9talement horizontal de la distribution d&rsquo;Erlang, une dualit\u00e9 qui exige une vigilance rigoureuse lors de la manipulation des \u00e9quations ou de l&rsquo;utilisation de progiciels statistiques. La premi\u00e8re convention repose sur le <strong>param\u00e8tre de taux<\/strong> &lambda; (souvent d\u00e9sign\u00e9 par le terme anglais <em>rate parameter<\/em>). Ce param\u00e8tre quantifie la fr\u00e9quence ou l&rsquo;intensit\u00e9 temporelle \u00e0 laquelle chaque phase \u00e9l\u00e9mentaire est franchie ; son unit\u00e9 physique est l&rsquo;inverse d&rsquo;un temps (par exemple, nombre de transitions par seconde, s<sup>-1<\/sup>).<\/p>\n<p>La seconde convention recourt au <strong>param\u00e8tre d&rsquo;\u00e9chelle<\/strong> &mu; ou &theta; (<em>scale parameter<\/em>), d\u00e9fini comme l&rsquo;inverse strict du taux :<\/p>\n<p><em>&mu; = 1 \/ &lambda;<\/em><\/p>\n<p>Le param\u00e8tre d&rsquo;\u00e9chelle poss\u00e8de quant \u00e0 lui la dimension physique d&rsquo;une dur\u00e9e (secondes, heures, jours) et correspond exactement au temps moyen n\u00e9cessaire pour achever une seule des <em>k<\/em> phases \u00e9l\u00e9mentaires. L&rsquo;expression de la fonction de densit\u00e9 d&rsquo;Erlang sous sa param\u00e9trisation d&rsquo;\u00e9chelle s&rsquo;\u00e9crit alors :<\/p>\n<p><em>f(x; k, &mu;) = [ x<sup>k-1<\/sup> &middot; e<sup>-x \/ &mu;<\/sup> ] \/ [ &mu;<sup>k<\/sup> &middot; (k &#8211; 1)! ]<\/em><\/p>\n<p>L&rsquo;impact de ces param\u00e8tres sur la repr\u00e9sentation g\u00e9om\u00e9trique est parfaitement pr\u00e9visible. Une \u00e9l\u00e9vation du taux &lambda; (ou une diminution sym\u00e9trique de l&rsquo;\u00e9chelle &mu;) comprime horizontalement la distribution le long de l&rsquo;axe temporel et en augmente corr\u00e9lativement la hauteur maximale en raison de la contrainte d&rsquo;int\u00e9grale unitaire. \u00c0 l&rsquo;inverse, une augmentation de l&rsquo;\u00e9chelle &mu; \u00e9tire la courbe vers la droite, \u00e9talant la dispersion temporelle des observations et abaissant la valeur modale.<\/p>\n<h3>3.3 Dynamique comparative des courbes selon les variations param\u00e9triques<\/h3>\n<p>L&rsquo;analyse visuelle conjointe des variations de <em>k<\/em> et de &lambda; met en lumi\u00e8re la remarquable plasticit\u00e9 de la famille des distributions d&rsquo;Erlang. En ajustant ad\u00e9quatement ce tandem de param\u00e8tres, l&rsquo;analyste peut reproduire une gamme \u00e9tendue de comportements probabilistes, allant de syst\u00e8mes \u00e0 forte impr\u00e9visibilit\u00e9 temporelle jusqu&rsquo;\u00e0 des m\u00e9canismes temporels quasi d\u00e9terministes. Le tableau ci-dessous synth\u00e9tise les comportements morphologiques en fonction des r\u00e9gimes param\u00e9triques :<\/p>\n<ul>\n<li><strong>R\u00e9gime k = 1 (Loi exponentielle) :<\/strong> Dispersion maximale, absence totale de m\u00e9moire, valeur modale situ\u00e9e en <em>x = 0<\/em>, pente initiale strictement n\u00e9gative \u00e9gale \u00e0 <em>-&lambda;<sup>2<\/sup><\/em>.<\/li>\n<li><strong>R\u00e9gime pour k mod\u00e9r\u00e9 (k &isin; [2, 5]) :<\/strong> Apparition d&rsquo;un maximum int\u00e9rieur bien marqu\u00e9, asym\u00e9trie positive forte (queue de distribution \u00e9tir\u00e9e vers la droite), probabilit\u00e9 nulle en <em>x = 0<\/em>.<\/li>\n<li><strong>R\u00e9gime pour k grand (k &gt; 20) :<\/strong> Amortissement prononc\u00e9 de l&rsquo;asym\u00e9trie, resserrement de l&rsquo;intervalle interquartile, convergence progressive de la forme globale vers une cloche gaussienne sym\u00e9trique.<\/li>\n<\/ul>\n<p>Le choix de privil\u00e9gier l&rsquo;une ou l&rsquo;autre param\u00e9trisation d\u00e9pend principalement du domaine d&rsquo;application. En ing\u00e9nierie des t\u00e9l\u00e9communications, en th\u00e9orie des r\u00e9seaux et dans l&rsquo;\u00e9tude des files d&rsquo;attente markoviennes, l&rsquo;usage du taux d&rsquo;intensit\u00e9 &lambda; est quasi syst\u00e9matique, car il s&rsquo;int\u00e8gre harmonieusement dans les matrices de g\u00e9n\u00e9rateurs infinit\u00e9simaux (matrices <em>Q<\/em>). En revanche, en fiabilit\u00e9 industrielle, en chronom\u00e9trie exp\u00e9rimentale et en m\u00e9trologie, l&#8217;emploi de l&rsquo;\u00e9chelle temporelle &mu; est souvent privil\u00e9gi\u00e9 en raison de son interpr\u00e9tation directe comme d\u00e9lai moyen de r\u00e9alisation.<\/p>\n<h2>4. Propri\u00e9t\u00e9s statistiques fondamentales et moments d&rsquo;ordre sup\u00e9rieur<\/h2>\n<h3>4.1 Esp\u00e9rance math\u00e9matique et mode<\/h3>\n<p>L&rsquo;esp\u00e9rance math\u00e9matique, qui caract\u00e9rise le centre de gravit\u00e9 probabiliste de la variable al\u00e9atoire, se calcule directement \u00e0 partir de la propri\u00e9t\u00e9 de lin\u00e9arit\u00e9 de l&rsquo;esp\u00e9rance appliqu\u00e9e \u00e0 la somme des <em>k<\/em> variables exponentielles composantes. Sachant que l&rsquo;esp\u00e9rance d&rsquo;une variable exponentielle <em>Y<sub>i<\/sub> ~ Exp(&lambda;)<\/em> est rigoureusement \u00e9gale \u00e0 <em>1 \/ &lambda;<\/em>, l&rsquo;esp\u00e9rance d&rsquo;une variable <em>X ~ Erlang(k, &lambda;)<\/em> est donn\u00e9e par :<\/p>\n<p><em>E[X] = E[&sum;<sub>i=1<\/sub><sup>k<\/sup> Y<sub>i<\/sub>] = &sum;<sub>i=1<\/sub><sup>k<\/sup> E[Y<sub>i<\/sub>] = k &middot; (1 \/ &lambda;) = k \/ &lambda;<\/em><\/p>\n<p>Sous la param\u00e9trisation d&rsquo;\u00e9chelle, cette relation s&rsquo;exprime simplement comme <em>E[X] = k &middot; &mu;<\/em>. La d\u00e9rivation du mode, d\u00e9fini comme la valeur de <em>x<\/em> en laquelle la fonction de densit\u00e9 de probabilit\u00e9 atteint son maximum absolu, s&rsquo;obtient en annulant la d\u00e9riv\u00e9e premi\u00e8re de <em>f(x; k, &lambda;)<\/em> par rapport \u00e0 <em>x<\/em> sur l&rsquo;intervalle ouvert ]0, &infin;[. En diff\u00e9renciant l&rsquo;expression analytique :<\/p>\n<p><em>d\/dx [ f(x; k, &lambda;) ] = [ &lambda;<sup>k<\/sup> \/ (k &#8211; 1)! ] &middot; d\/dx [ x<sup>k-1<\/sup> e<sup>-&lambda;x<\/sup> ] = [ &lambda;<sup>k<\/sup> \/ (k &#8211; 1)! ] &middot; [ (k &#8211; 1) x<sup>k-2<\/sup> e<sup>-&lambda;x<\/sup> &#8211; &lambda; x<sup>k-1<\/sup> e<sup>-&lambda;x<\/sup> ]<\/em><\/p>\n<p>En \u00e9galisant cette expression \u00e0 z\u00e9ro pour <em>x &gt; 0<\/em>, on factorise par <em>x<sup>k-2<\/sup> e<sup>-&lambda;x<\/sup><\/em>, ce qui m\u00e8ne directement \u00e0 l&rsquo;\u00e9quation lin\u00e9aire <em>(k &#8211; 1) &#8211; &lambda;x = 0<\/em>. On en d\u00e9duit l&rsquo;expression exacte de la valeur modale :<\/p>\n<p><em>Mode(X) = (k &#8211; 1) \/ &lambda; = (k &#8211; 1) &middot; &mu; &nbsp; (pour k &ge; 1)<\/em><\/p>\n<p>Lorsque <em>k = 1<\/em>, le mode vaut rigoureusement 0, ce qui concorde avec la densit\u00e9 exponentielle d\u00e9croissante. Pour toute valeur <em>k &ge; 2<\/em>, le mode est strictement positif mais demeure syst\u00e9matiquement inf\u00e9rieur \u00e0 la moyenne arithm\u00e9tique <em>k \/ &lambda;<\/em>, ce qui d\u00e9montre analytiquement l&rsquo;asym\u00e9trie droite inh\u00e9rente \u00e0 la distribution.<\/p>\n<h3>4.2 Variance et indicateurs de dispersion<\/h3>\n<p>L&rsquo;ind\u00e9pendance mutuelle des variables exponentielles somm\u00e9 permet d&rsquo;exploiter la propri\u00e9t\u00e9 d&rsquo;additivit\u00e9 de la variance pour les op\u00e9rateurs stochastiques non corr\u00e9l\u00e9s. Pour une variable exponentielle \u00e9l\u00e9mentaire <em>Y<sub>i<\/sub><\/em>, la variance th\u00e9orique est donn\u00e9e par <em>Var(Y<sub>i<\/sub>) = 1 \/ &lambda;<sup>2<\/sup><\/em>. Par cons\u00e9quent, la variance de la distribution d&rsquo;Erlang s&rsquo;exprime selon la formule \u00e9l\u00e9gante :<\/p>\n<p><em>Var(X) = Var(&sum;<sub>i=1<\/sub><sup>k<\/sup> Y<sub>i<\/sub>) = &sum;<sub>i=1<\/sub><sup>k<\/sup> Var(Y<sub>i<\/sub>) = k \/ &lambda;<sup>2<\/sup> = k &middot; &mu;<sup>2<\/sup><\/em><\/p>\n<p>L&rsquo;\u00e9cart-type &sigma;(X), qui mesure l&rsquo;incertitude absolue ou la dispersion des temps autour de la moyenne, est obtenu par la racine carr\u00e9e de la variance, soit :<\/p>\n<p><em>&sigma;(X) = &radic;k \/ &lambda; = &mu; &middot; &radic;k<\/em><\/p>\n<p>Un indicateur synth\u00e9tique crucial en th\u00e9orie des files d&rsquo;attente et en fiabilit\u00e9 est le <strong>coefficient de variation<\/strong> (not\u00e9 <em>c<sub>v<\/sub><\/em>), d\u00e9fini comme le rapport sans dimension de l&rsquo;\u00e9cart-type sur l&rsquo;esp\u00e9rance math\u00e9matique :<\/p>\n<p><em>c<sub>v<\/sub> = &sigma;(X) \/ E[X] = (&radic;k \/ &lambda;) \/ (k \/ &lambda;) = 1 \/ &radic;k<\/em><\/p>\n<p>Cette relation d\u00e9montre un principe d&rsquo;ing\u00e9nierie fondamental : le coefficient de variation de la distribution d&rsquo;Erlang est strictement inf\u00e9rieur ou \u00e9gal \u00e0 1 pour tout <em>k &ge; 1<\/em>, et il tend asymptotiquement vers z\u00e9ro lorsque le nombre d&rsquo;\u00e9tapes <em>k<\/em> tend vers l&rsquo;infini. Cela prouve math\u00e9matiquement que la distribution d&rsquo;Erlang est structurellement \u00ab\u00a0plus r\u00e9guli\u00e8re\u00a0\u00bb que la distribution exponentielle (pour laquelle <em>c<sub>v<\/sub> = 1<\/em>), ce qui en fait un outil de mod\u00e9lisation de choix pour les processus industriels ma\u00eetris\u00e9s pr\u00e9sentant une variabilit\u00e9 r\u00e9siduelle mod\u00e9r\u00e9e.<\/p>\n<h3>4.3 Asym\u00e9trie et aplatissement (Skewness et Kurtosis)<\/h3>\n<p>Les propri\u00e9t\u00e9s morphologiques fines de la distribution d&rsquo;Erlang sont d\u00e9crites par ses moments centr\u00e9s normalis\u00e9s d&rsquo;ordre trois et quatre. Le coefficient d&rsquo;asym\u00e9trie de Fisher (<em>skewness<\/em>), not\u00e9 &gamma;<sub>1<\/sub>, permet de quantifier le degr\u00e9 de non-sym\u00e9trie de la densit\u00e9 par rapport \u00e0 sa moyenne :<\/p>\n<p><em>&gamma;<sub>1<\/sub> = E[(X &#8211; E[X])<sup>3<\/sup>] \/ &sigma;<sup>3<\/sup>(X) = (2 \/ &radic;k)<\/em><\/p>\n<p>Cette valeur \u00e9tant strictement positive pour tout entier <em>k<\/em> fini, la distribution d&rsquo;Erlang est invariablement asym\u00e9trique \u00e0 droite, pr\u00e9sentant une queue de distribution plus \u00e9tir\u00e9e vers les valeurs temporelles \u00e9lev\u00e9es. Toutefois, lorsque <em>k<\/em> devient arbitrairement grand, &gamma;<sub>1<\/sub> d\u00e9cro\u00eet vers z\u00e9ro \u00e0 un taux proportionnel \u00e0 <em>1 \/ &radic;k<\/em>, traduisant une sym\u00e9trisation progressive de la cloche.<\/p>\n<p>Le coefficient d&rsquo;aplatissement de Pearson ou kurtosis (&beta;<sub>2<\/sub>) et l&rsquo;exc\u00e8s de kurtosis (&gamma;<sub>2<\/sub> = &beta;<sub>2<\/sub> &#8211; 3) mesurent la lourdeur des queues de probabilit\u00e9 et l&rsquo;acuit\u00e9 de son sommet par comparaison avec la loi normale :<\/p>\n<p><em>&beta;<sub>2<\/sub> = E[(X &#8211; E[X])<sup>4<\/sup>] \/ &sigma;<sup>4<\/sup>(X) = 3 + (6 \/ k)<\/em><\/p>\n<p><em>&gamma;<sub>2<\/sub> = &beta;<sub>2<\/sub> &#8211; 3 = 6 \/ k<\/em><\/p>\n<p>L&rsquo;exc\u00e8s de kurtosis est strictement positif pour toute valeur finie de <em>k<\/em>, ce qui classe formellement la distribution d&rsquo;Erlang parmi les lois leptokurtiques. La convergence conjointe de l&rsquo;asym\u00e9trie vers 0 et de l&rsquo;exc\u00e8s d&rsquo;aplatissement vers 0 lorsque <em>k &rarr; &infin;<\/em> constitue la signature alg\u00e9brique du th\u00e9or\u00e8me central limite : convenablement centr\u00e9e et r\u00e9duite, la distribution d&rsquo;Erlang converge en loi vers la distribution normale standard <em>N(0, 1)<\/em> lorsque son nombre de phases devient infini.<\/p>\n<h3>4.4 Fonction g\u00e9n\u00e9ratrice des moments et fonction caract\u00e9ristique<\/h3>\n<p>La <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Moment-generating_function\">fonction g\u00e9n\u00e9ratrice des moments<\/a> (MGF), not\u00e9e <em>M<sub>X<\/sub>(t)<\/em>, offre un m\u00e9canisme puissant pour caract\u00e9riser l&rsquo;ensemble des moments d&rsquo;une variable stochastique par simple d\u00e9rivation analytique. Pour une variable <em>X ~ Erlang(k, &lambda;)<\/em>, elle est d\u00e9finie par l&rsquo;esp\u00e9rance <em>E[e<sup>tX<\/sup>]<\/em> sous r\u00e9serve de convergence :<\/p>\n<p><em>M<sub>X<\/sub>(t) = &int;<sub>0<\/sub><sup>&infin;<\/sup> e<sup>tx<\/sup> &middot; [ &lambda;<sup>k<\/sup> x<sup>k-1<\/sup> e<sup>-&lambda;x<\/sup> \/ (k &#8211; 1)! ] dx = [ &lambda;<sup>k<\/sup> \/ (k &#8211; 1)! ] &int;<sub>0<\/sub><sup>&infin;<\/sup> x<sup>k-1<\/sup> e<sup>-(&lambda; &#8211; t)x<\/sup> dx<\/em><\/p>\n<p>Cette int\u00e9grale converge de mani\u00e8re absolue si et seulement si la partie r\u00e9elle de <em>t<\/em> v\u00e9rifie la condition stricte <em>t &lt; &lambda;<\/em>. En effectuant le changement de variable <em>u = (&lambda; &#8211; t)x<\/em>, on fait r\u00e9appara\u00eetre la fonction gamma d&rsquo;Euler &Gamma;(k) = (k &#8211; 1)!, ce qui donne la formulation compacte :<\/p>\n<p><em>M<sub>X<\/sub>(t) = ( 1 &#8211; t \/ &lambda; )<sup>-k<\/sup> = [ &lambda; \/ (&lambda; &#8211; t) ]<sup>k<\/sup><\/em><\/p>\n<p>Cette factorisation sous la forme d&rsquo;une puissance <em>k<\/em>-i\u00e8me confirme imm\u00e9diatement la structure convolutive de la loi, puisque la MGF d&rsquo;une variable exponentielle isol\u00e9e vaut <em>(1 &#8211; t \/ &lambda;)<sup>-1<\/sup><\/em>. La fonction caract\u00e9ristique &phi;<sub>X<\/sub>(t), quant \u00e0 elle, s&rsquo;obtient par prolongement analytique sur l&rsquo;axe imaginaire pour tout <em>t &isin; &amp;reals;<\/em> :<\/p>\n<p><em>&phi;<sub>X<\/sub>(t) = E[e<sup>itX<\/sup>] = ( 1 &#8211; it \/ &lambda; )<sup>-k<\/sup> = [ &lambda; \/ (&lambda; &#8211; it) ]<sup>k<\/sup><\/em><\/p>\n<p>La fonction caract\u00e9ristique pr\u00e9sente l&rsquo;avantage majeur d&rsquo;\u00eatre d\u00e9finie sur l&rsquo;ensemble de la droite r\u00e9elle sans restriction de convergence. Elle simplifie consid\u00e9rablement les calculs th\u00e9oriques impliquant des combinaisons lin\u00e9aires de dur\u00e9es d&rsquo;Erlang ind\u00e9pendantes.<\/p>\n<h2>5. Fonction de r\u00e9partition et calcul des probabilit\u00e9s cumul\u00e9es<\/h2>\n<h3>5.1 Expression analytique de la fonction de r\u00e9partition<\/h3>\n<p>La fonction de r\u00e9partition (CDF, pour <em>Cumulative Distribution Function<\/em>), not\u00e9e <em>F(x; k, &lambda;)<\/em>, quantifie la probabilit\u00e9 cumul\u00e9e que la variable al\u00e9atoire <em>X<\/em> prenne une valeur inf\u00e9rieure ou \u00e9gale \u00e0 un seuil temporel donn\u00e9 <em>x<\/em>. Par d\u00e9finition, elle r\u00e9sulte de l&rsquo;int\u00e9gration de la fonction de densit\u00e9 depuis l&rsquo;origine jusqu&rsquo;\u00e0 la borne sup\u00e9rieure <em>x<\/em> :<\/p>\n<p><em>F(x; k, &lambda;) = P(X &le; x) = &int;<sub>0<\/sub><sup>x<\/sup> [ &lambda;<sup>k<\/sup> u<sup>k-1<\/sup> e<sup>-&lambda;u<\/sup> \/ (k &#8211; 1)! ] du<\/em><\/p>\n<p>Pour r\u00e9soudre cette int\u00e9grale sans recourir \u00e0 des quadratures num\u00e9riques, on applique une succession de <em>k &#8211; 1<\/em> int\u00e9grations par parties cons\u00e9cutives. \u00c0 chaque \u00e9tape, la d\u00e9rivation abaisse la puissance de <em>u<\/em> tout en int\u00e9grant le terme exponentiel. Au terme de cette r\u00e9duction r\u00e9cursive, la fonction de r\u00e9partition d&rsquo;Erlang se formalise sous la forme \u00e9l\u00e9gante d&rsquo;un compl\u00e9ment \u00e0 un d&rsquo;une somme finie de termes pond\u00e9r\u00e9s :<\/p>\n<p><em>F(x; k, &lambda;) = 1 &#8211; &sum;<sub>n=0<\/sub><sup>k-1<\/sup> [ (&lambda;x)<sup>n<\/sup> \/ n! ] &middot; e<sup>-&lambda;x<\/sup> &nbsp; (pour x &ge; 0)<\/em><\/p>\n<p>Cette formule ferm\u00e9e rev\u00eat une importance capitale : elle \u00e9limine le besoin d&rsquo;\u00e9valuer des int\u00e9grales transcendantes continues en r\u00e9duisant le calcul probabiliste \u00e0 l&rsquo;\u00e9valuation d&rsquo;une somme discr\u00e8te de <em>k<\/em> termes \u00e9l\u00e9mentaires. Pour toute valeur n\u00e9gative <em>x &lt; 0<\/em>, la fonction de r\u00e9partition est identiquement nulle : <em>F(x; k, &lambda;) = 0<\/em>.<\/p>\n<h3>5.2 Lien avec la fonction gamma incompl\u00e8te<\/h3>\n<p>Sur le plan de l&rsquo;analyse fonctionnelle avanc\u00e9e, l&rsquo;int\u00e9grale d\u00e9finissant la fonction de r\u00e9partition de la loi d&rsquo;Erlang peut \u00eatre exprim\u00e9e \u00e0 l&rsquo;aide des fonctions eul\u00e9riennes g\u00e9n\u00e9ralis\u00e9es. En proc\u00e9dant au changement de variable adimensionnel <em>s = &lambda;u<\/em> dans l&rsquo;int\u00e9grale de d\u00e9finition, on obtient :<\/p>\n<p><em>F(x; k, &lambda;) = [ 1 \/ (k &#8211; 1)! ] &int;<sub>0<\/sub><sup>&lambda;x<\/sup> s<sup>k-1<\/sup> e<sup>-s<\/sup> ds = &gamma;(k, &lambda;x) \/ &Gamma;(k) = P(k, &lambda;x)<\/em><\/p>\n<p>Dans cette formulation, &gamma;(k, &lambda;x) d\u00e9signe la <strong>fonction gamma incompl\u00e8te inf\u00e9rieure<\/strong>, d\u00e9finie analytiquement par l&rsquo;int\u00e9grale param\u00e9tr\u00e9e :<\/p>\n<p><em>&gamma;(a, z) = &int;<sub>0<\/sub><sup>z<\/sup> t<sup>a-1<\/sup> e<sup>-t<\/sup> dt<\/em><\/p>\n<p>Le terme <em>P(k, &lambda;x) = &gamma;(k, &lambda;x) \/ &Gamma;(k)<\/em> correspond \u00e0 la fonction gamma incompl\u00e8te r\u00e9gularis\u00e9e. Cette correspondance math\u00e9matique rigoureuse d\u00e9montre que la somme finie de Poisson \u00e9tablie au paragraphe pr\u00e9c\u00e9dent n&rsquo;est autre que l&rsquo;expression d\u00e9velopp\u00e9e de la fonction gamma incompl\u00e8te lorsque son premier argument est un entier naturel non nul.<\/p>\n<p>D&rsquo;un point de vue algorithmique et informatique, cette passerelle vers les fonctions eul\u00e9riennes est particuli\u00e8rement pr\u00e9cieuse. Elle permet aux biblioth\u00e8ques de calcul num\u00e9rique haute performance d&rsquo;\u00e9valuer la fonction de r\u00e9partition de distributions d&rsquo;Erlang d&rsquo;ordre tr\u00e8s \u00e9lev\u00e9 (par exemple pour <em>k &gt; 500<\/em>) sans subir les instabilit\u00e9s num\u00e9riques ou les d\u00e9passements de capacit\u00e9 arithm\u00e9tique (<em>overflow<\/em>) inh\u00e9rents au calcul direct des factorielles ou des termes de Poisson exponentiels.<\/p>\n<h3>5.3 Fonction de survie et taux de d\u00e9faillance<\/h3>\n<p>En ing\u00e9nierie de la fiabilit\u00e9, en actuariat et en analyse de survie, la fonction d&rsquo;int\u00e9r\u00eat primaire est souvent la <strong>fonction de survie<\/strong> ou de fiabilit\u00e9 (not\u00e9e <em>R(x)<\/em> ou <em>S(x)<\/em>), qui quantifie la probabilit\u00e9 qu&rsquo;un syst\u00e8me ne connaisse aucune d\u00e9faillance avant l&rsquo;instant <em>x<\/em>. Par compl\u00e9mentarit\u00e9 avec la fonction de r\u00e9partition, elle s&rsquo;\u00e9crit :<\/p>\n<p><em>R(x; k, &lambda;) = P(X &gt; x) = 1 &#8211; F(x; k, &lambda;) = &sum;<sub>n=0<\/sub><sup>k-1<\/sup> [ (&lambda;x)<sup>n<\/sup> \/ n! ] &middot; e<sup>-&lambda;x<\/sup><\/em><\/p>\n<p>Le <strong>taux de d\u00e9faillance instantan\u00e9<\/strong> (ou taux de hasard, d\u00e9sign\u00e9 par <em>h(x)<\/em>), qui mesure le risque conditionnel d&rsquo;extinction dans l&rsquo;intervalle infinit\u00e9simal suivant sachant la survie jusqu&rsquo;\u00e0 l&rsquo;instant <em>x<\/em>, est d\u00e9fini par le ratio entre la densit\u00e9 de probabilit\u00e9 et la fonction de survie :<\/p>\n<p><em>h(x; k, &lambda;) = f(x; k, &lambda;) \/ R(x; k, &lambda;) = [ &lambda;<sup>k<\/sup> x<sup>k-1<\/sup> \/ (k &#8211; 1)! ] \/ [ &sum;<sub>n=0<\/sub><sup>k-1<\/sup> (&lambda;x)<sup>n<\/sup> \/ n! ]<\/em><\/p>\n<p>L&rsquo;analyse dynamique de cette fonction <em>h(x)<\/em> r\u00e9v\u00e8le une rupture conceptuelle majeure :<\/p>\n<ul>\n<li>Lorsque <em>k = 1<\/em>, le taux de d\u00e9faillance est rigoureusement constant : <em>h(x) = &lambda;<\/em>. C&rsquo;est l&rsquo;expression exacte de la propri\u00e9t\u00e9 d&rsquo;amn\u00e9sie (ou absence de m\u00e9moire) de la loi exponentielle.<\/li>\n<li>D\u00e8s que <em>k &ge; 2<\/em>, le taux de d\u00e9faillance instantan\u00e9 v\u00e9rifie <em>h(0) = 0<\/em> et constitue une fonction <strong>strictement croissante et monotone<\/strong> de la variable temporelle <em>x<\/em>, convergeant asymptotiquement vers une valeur plafond \u00e9gale \u00e0 &lambda; lorsque <em>x &rarr; &infin;<\/em> :<\/li>\n<\/ul>\n<p><em>lim<sub>x&rarr;&infin;<\/sub> h(x; k, &lambda;) = &lambda;<\/em><\/p>\n<p>Ce profil croissant caract\u00e9rise un ph\u00e9nom\u00e8ne d&rsquo;usure ou de vieillissement positif (classe <em>IFR<\/em> pour <em>Increasing Failure Rate<\/em>). Contrairement \u00e0 un composant exponentiel qui ne s&rsquo;use jamais, un syst\u00e8me mod\u00e9lis\u00e9 par une loi d&rsquo;Erlang d&rsquo;ordre sup\u00e9rieur \u00e0 un accumule de la fatigue stochastique \u00e0 mesure que les phases interm\u00e9diaires sont successivement franchies.<\/p>\n<h2>6. Relations structurales avec d&rsquo;autres lois de probabilit\u00e9<\/h2>\n<h3>6.1 La distribution Gamma comme g\u00e9n\u00e9ralisation continue<\/h3>\n<p>La relation reliant la distribution d&rsquo;Erlang \u00e0 la <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Gamma_distribution\">distribution Gamma<\/a> est analogue \u00e0 celle reliant la factorielle standard \u00e0 la fonction Gamma d&rsquo;Euler. La loi Gamma est une distribution continue caract\u00e9ris\u00e9e par une densit\u00e9 identique dans sa forme analytique, mais dont le param\u00e8tre de forme \u2014 traditionnellement not\u00e9 &alpha; au lieu de <em>k<\/em> \u2014 n&rsquo;est plus assujetti \u00e0 la contrainte d&rsquo;int\u00e9gralit\u00e9 et peut prendre n&rsquo;importe quelle valeur r\u00e9elle strictement positive (&alpha; &isin; &amp;reals;<sub>+<\/sub><sup>*<\/sup>) :<\/p>\n<p><em>f<sub>Gamma<\/sub>(x; &alpha;, &beta;) = [ &beta;<sup>&alpha;<\/sup> &middot; x<sup>&alpha;-1<\/sup> &middot; e<sup>-&beta;x<\/sup> ] \/ &Gamma;(&alpha;)<\/em><\/p>\n<p>La distribution d&rsquo;Erlang est donc le cas particulier de la loi Gamma obtenu lorsque &alpha; = <em>k<\/em> &isin; &amp;naturals;<sup>*<\/sup>. Cette restriction aux entiers conf\u00e8re \u00e0 la distribution d&rsquo;Erlang une signification physique sp\u00e9cifique que la loi Gamma g\u00e9n\u00e9rale ne poss\u00e8de pas : celle d&rsquo;une somme finie d&rsquo;\u00e9tapes markoviennes. Si la loi Gamma permet un ajustement statistique empirique parfois plus souple en acceptant des param\u00e8tres continus, elle perd en contrepartie la possibilit\u00e9 d&rsquo;\u00eatre interpr\u00e9t\u00e9e directement comme une cha\u00eene d&rsquo;\u00e9tapes exponentielles homog\u00e8nes dans un espace d&rsquo;\u00e9tats discret.<\/p>\n<h3>6.2 Identit\u00e9 avec la loi exponentielle<\/h3>\n<p>La distribution exponentielle constitue l&rsquo;\u00e9l\u00e9ment atomique de la distribution d&rsquo;Erlang. Lorsque le param\u00e8tre de forme est fix\u00e9 \u00e0 l&rsquo;unit\u00e9 (<em>k = 1<\/em>), la formule g\u00e9n\u00e9rale de la densit\u00e9 se simplifie instantan\u00e9ment :<\/p>\n<p><em>f(x; 1, &lambda;) = [ &lambda;<sup>1<\/sup> &middot; x<sup>1-1<\/sup> &middot; e<sup>-&lambda;x<\/sup> ] \/ (1 &#8211; 1)! = [ &lambda; &middot; 1 &middot; e<sup>-&lambda;x<\/sup> ] \/ 0! = &lambda; e<sup>-&lambda;x<\/sup><\/em><\/p>\n<p>Cette identit\u00e9 alg\u00e9brique r\u00e9tablit la loi exponentielle standard. Ce n&rsquo;est qu&rsquo;exclusivement dans ce sous-cas d\u00e9g\u00e9n\u00e9r\u00e9 que la distribution d&rsquo;Erlang b\u00e9n\u00e9ficie de la propri\u00e9t\u00e9 d&rsquo;absence de m\u00e9moire (<em>memoryless property<\/em>), formalis\u00e9e par l&rsquo;\u00e9galit\u00e9 <em>P(X &gt; s + t | X &gt; s) = P(X &gt; t)<\/em>. D\u00e8s que <em>k &ge; 2<\/em>, le passage d&rsquo;une phase \u00e0 la suivante garde la m\u00e9moire du temps \u00e9coul\u00e9, car le syst\u00e8me s&rsquo;est rapproch\u00e9 de l&rsquo;ach\u00e8vement de la derni\u00e8re phase. La distribution d&rsquo;Erlang peut ainsi \u00eatre comprise comme un m\u00e9canisme math\u00e9matique visant \u00e0 d\u00e9passer et contr\u00f4ler la contrainte d&rsquo;amn\u00e9sie de l&rsquo;exponentielle tout en conservant sa tractabilit\u00e9 alg\u00e9brique.<\/p>\n<h3>6.3 Interactions avec le processus de Poisson<\/h3>\n<p>Il existe une dualit\u00e9 stochastique profonde et rigoureuse entre la distribution d&rsquo;Erlang et le <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Poisson_point_process\">processus de Poisson<\/a>. Consid\u00e9rons un processus ponctuel de Poisson homog\u00e8ne d&rsquo;intensit\u00e9 constante &lambda;, o\u00f9 les \u00e9v\u00e9nements successifs surviennent de fa\u00e7on totalement al\u00e9atoire et ind\u00e9pendante dans le temps. Par construction m\u00eame du processus, les intervalles temporels s\u00e9parant deux \u00e9v\u00e9nements cons\u00e9cutifs sont des variables al\u00e9atoires exponentielles ind\u00e9pendantes de param\u00e8tre &lambda;.<\/p>\n<p>En cons\u00e9quence directe, le temps d&rsquo;attente <em>T<sub>k<\/sub><\/em> n\u00e9cessaire pour observer pr\u00e9cis\u00e9ment l&rsquo;occurrence du <em>k<\/em>-i\u00e8me \u00e9v\u00e9nement du processus de Poisson est la somme de <em>k<\/em> intervalles inter-arriv\u00e9es exponentiels. Par d\u00e9finition, <em>T<sub>k<\/sub><\/em> suit donc rigoureusement une distribution d&rsquo;Erlang de param\u00e8tres <em>k<\/em> et &lambda; :<\/p>\n<p><em>T<sub>k<\/sub> ~ Erlang(k, &lambda;)<\/em><\/p>\n<p>Cette relation de r\u00e9ciprocit\u00e9 permet de d\u00e9montrer la fonction de r\u00e9partition d&rsquo;Erlang par un simple raisonnement logique de comptage. Affirmer que le temps d&rsquo;attente pour voir survenir le <em>k<\/em>-i\u00e8me \u00e9v\u00e9nement d\u00e9passe un instant donn\u00e9 <em>t<\/em> (c&rsquo;est-\u00e0-dire <em>T<sub>k<\/sub> &gt; t<\/em>) \u00e9quivaut rigoureusement \u00e0 affirmer que le nombre total d&rsquo;\u00e9v\u00e9nements <em>N(t)<\/em> survenus durant l&rsquo;intervalle [0, t] est strictement inf\u00e9rieur \u00e0 <em>k<\/em> (soit au plus <em>k &#8211; 1<\/em> \u00e9v\u00e9nements). En exploitant la loi de Poisson de param\u00e8tre &lambda;t pour la variable de comptage <em>N(t)<\/em>, on \u00e9crit imm\u00e9diatement :<\/p>\n<p><em>P(T<sub>k<\/sub> &gt; t) = P(N(t) &le; k &#8211; 1) = &sum;<sub>n=0<\/sub><sup>k-1<\/sup> [ (&lambda;t)<sup>n<\/sup> \/ n! ] &middot; e<sup>-&lambda;t<\/sup><\/em><\/p>\n<p>Cette d\u00e9monstration \u00e9l\u00e9gante \u00e9tablit l&rsquo;isomorphisme parfait reliant la dynamique continue du temps d&rsquo;attente d&rsquo;Erlang \u00e0 l&rsquo;alg\u00e8bre discr\u00e8te du processus de comptage de Poisson.<\/p>\n<h3>6.4 Connexion avec la loi du Chi-deux<\/h3>\n<p>La distribution d&rsquo;Erlang entretient \u00e9galement une parent\u00e9 structurelle \u00e9troite avec la distribution du Chi-deux (&chi;<sup>2<\/sup>), un pivot m\u00e9thodologique de la statistique inf\u00e9rentielle. La distribution du Chi-deux \u00e0 &nu; degr\u00e9s de libert\u00e9, not\u00e9e &chi;<sup>2<\/sup>(&nu;), est un cas particulier de la loi Gamma o\u00f9 le param\u00e8tre de forme vaut &alpha; = &nu; \/ 2 et le param\u00e8tre de taux vaut &beta; = 1 \/ 2 (ou \u00e9chelle &mu; = 2).<\/p>\n<p>Si l&rsquo;on consid\u00e8re le cas o\u00f9 le nombre de degr\u00e9s de libert\u00e9 de la loi du Chi-deux est un entier pair, c&rsquo;est-\u00e0-dire &nu; = 2k (avec <em>k &isin; &amp;naturals;<sup>*<\/sup><\/em>), la distribution &chi;<sup>2<\/sup>(2k) co\u00efncide exactement avec une distribution d&rsquo;Erlang de param\u00e8tres <em>k<\/em> et d&rsquo;intensit\u00e9 &lambda; = 1\/2 :<\/p>\n<p><em>&chi;<sup>2<\/sup>(2k) &equiv; Erlang(k, &lambda; = 1\/2)<\/em><\/p>\n<p>De mani\u00e8re plus g\u00e9n\u00e9rale, si une variable al\u00e9atoire <em>X<\/em> suit une loi d&rsquo;Erlang quelconque <em>X ~ Erlang(k, &lambda;)<\/em>, la transformation par changement d&rsquo;\u00e9chelle scalaire <em>Y = 2&lambda;X<\/em> produit une nouvelle variable al\u00e9atoire qui suit rigoureusement une distribution du Chi-deux \u00e0 2k degr\u00e9s de libert\u00e9 :<\/p>\n<p><em>2&lambda;X ~ &chi;<sup>2<\/sup>(2k)<\/em><\/p>\n<p>Cette \u00e9quivalence est d&rsquo;une grande port\u00e9e pratique en statistique appliqu\u00e9e. Elle permet d&rsquo;exploiter imm\u00e9diatement les tables num\u00e9riques hautement optimis\u00e9es et les routines informatiques consacr\u00e9es \u00e0 la loi du Chi-deux pour construire des intervalles de confiance exacts ou r\u00e9aliser des tests d&rsquo;ad\u00e9quation sur des param\u00e8tres de mod\u00e8les erlangiens.<\/p>\n<h2>7. R\u00f4le central dans la mod\u00e9lisation des files d&rsquo;attente<\/h2>\n<h3>7.1 La m\u00e9thode des phases fictives d&rsquo;Erlang<\/h3>\n<p>L&rsquo;une des contributions m\u00e9thodologiques les plus novatrices apport\u00e9es par Agner Krarup Erlang \u00e0 l&rsquo;analyse stochastique est sans conteste la <strong>m\u00e9thode des phases fictives<\/strong> (parfois qualifi\u00e9e de m\u00e9thode des \u00e9tapes). Dans la th\u00e9orie des files d&rsquo;attente, l&rsquo;hypoth\u00e8se selon laquelle les temps de service ou les intervalles inter-arriv\u00e9es suivent une distribution exponentielle est commode, car elle conf\u00e8re au syst\u00e8me la propri\u00e9t\u00e9 markovienne (processus sans m\u00e9moire) permettant une r\u00e9solution analytique ais\u00e9e via les cha\u00eenes de Markov \u00e0 temps continu (CTMC).<\/p>\n<p>Cependant, dans la quasi-totalit\u00e9 des contextes r\u00e9els (t\u00e9l\u00e9coms, industrie, logistique), supposer un temps de service exponentiel est faux : la dispersion est trop forte et le mode en z\u00e9ro est physiquement absurde. C&rsquo;est ici qu&rsquo;intervient le g\u00e9nie d&rsquo;Erlang : pour mod\u00e9liser une t\u00e2che dont la dur\u00e9e globale n&rsquo;est pas exponentielle mais pr\u00e9sente une variance plus restreinte, il propose de concevoir cette t\u00e2che comme le franchissement successif de <em>k<\/em> sous-\u00e9tapes exponentielles s\u00e9quentielles invisibles de l&rsquo;ext\u00e9rieur. D\u00e8s qu&rsquo;une phase s&rsquo;ach\u00e8ve, le client bascule instantan\u00e9ment dans la suivante ; le service complet n&rsquo;est valid\u00e9 que lorsque la <em>k<\/em>-i\u00e8me phase est purg\u00e9e.<\/p>\n<p>En enrichissant l&rsquo;espace des \u00e9tats du syst\u00e8me \u2014 c&rsquo;est-\u00e0-dire en ne m\u00e9morisant plus seulement le nombre total de clients pr\u00e9sents mais \u00e9galement l&rsquo;indice exact de la phase en cours d&rsquo;ex\u00e9cution \u2014 le processus stochastique global redevient strictement markovien. Cette technique d&rsquo;expansion de l&rsquo;espace d&rsquo;\u00e9tats permet de convertir des syst\u00e8mes non markoviens insolubles en syst\u00e8mes d&rsquo;\u00e9quations de bilan lin\u00e9aires \u00e0 flux stationnaires solubles par l&rsquo;alg\u00e8bre matricielle classique.<\/p>\n<h3>7.2 L&rsquo;architecture des syst\u00e8mes de files d&rsquo;attente M\/E_k\/1<\/h3>\n<p>Dans la notation standard de Kendall, le syst\u00e8me <strong>M\/E<sub>k<\/sub>\/1<\/strong> d\u00e9signe une file d&rsquo;attente mono-serveur dot\u00e9e d&rsquo;un processus d&rsquo;arriv\u00e9e poissonnien de taux moyen &Lambda; (la lettre <em>M<\/em> signifiant markovien), d&rsquo;un temps de service mod\u00e9lis\u00e9 selon une distribution d&rsquo;Erlang \u00e0 <em>k<\/em> phases de taux total &mu; (not\u00e9e <em>E<sub>k<\/sub><\/em>), et d&rsquo;un guichet unique de traitement. La dur\u00e9e moyenne d&rsquo;un service complet est <em>1 \/ &mu;<\/em>, et l&rsquo;intensit\u00e9 de trafic est caract\u00e9ris\u00e9e par le taux d&rsquo;occupation &rho; = &Lambda; \/ &mu; &lt; 1 pour garantir la stabilit\u00e9 de la file.<\/p>\n<p>L&rsquo;utilisation de la distribution d&rsquo;Erlang \u00e0 <em>k<\/em> phases pour mod\u00e9liser le temps de service exerce un impact b\u00e9n\u00e9fique mesurable sur la performance du syst\u00e8me. En vertu de la c\u00e9l\u00e8bre formule de Pollaczek-Khintchine, la longueur moyenne de la file d&rsquo;attente <em>L<sub>q<\/sub><\/em> (nombre moyen de clients en attente avant acc\u00e8s au serveur) s&rsquo;exprime en fonction du coefficient de variation <em>c<sub>v<\/sub><\/em> du temps de service :<\/p>\n<p><em>L<sub>q<\/sub> = [ &rho;<sup>2<\/sup> &middot; (1 + c<sub>v<\/sub><sup>2<\/sup>) ] \/ [ 2 &middot; (1 &#8211; &rho;) ]<\/em><\/p>\n<p>Comme d\u00e9montr\u00e9 pr\u00e9c\u00e9demment, pour une distribution d&rsquo;Erlang, nous avons <em>c<sub>v<\/sub><sup>2<\/sup> = 1 \/ k<\/em>. En injectant cette relation dans la formule de Pollaczek-Khintchine, la formule devient :<\/p>\n<p><em>L<sub>q<\/sub> = [ &rho;<sup>2<\/sup> \/ (2(1 &#8211; &rho;)) ] &middot; [ 1 + 1\/k ] = [ &rho;<sup>2<\/sup> &middot; (k + 1) ] \/ [ 2k &middot; (1 &#8211; &rho;) ]<\/em><\/p>\n<p>Cette \u00e9quation illustre math\u00e9matiquement l&rsquo;influence de la r\u00e9gularit\u00e9 du service sur la congestion du syst\u00e8me. Dans une file purement exponentielle <em>M\/M\/1<\/em> (\u00e9quivalente \u00e0 <em>k = 1<\/em>), le terme multiplicatif vaut <em>(1 + 1)\/2 = 1<\/em>. En augmentant la valeur de <em>k<\/em>, le terme <em>(k + 1) \/ (2k)<\/em> d\u00e9cro\u00eet de mani\u00e8re monotone pour converger vers <em>1\/2<\/em> lorsque <em>k &rarr; &infin;<\/em> (correspondant au syst\u00e8me d\u00e9terministe <em>M\/D\/1<\/em>). Ainsi, en substituant un service d&rsquo;Erlang d&rsquo;ordre \u00e9lev\u00e9 \u00e0 un service exponentiel, on parvient \u00e0 <strong>r\u00e9duire de moiti\u00e9<\/strong> la taille moyenne de la file d&rsquo;attente et le temps d&rsquo;attente des requ\u00eates sans modifier le d\u00e9bit nominal du serveur.<\/p>\n<h3>7.3 Syst\u00e8mes \u00e0 arriv\u00e9es erlangiennes E_k\/M\/1<\/h3>\n<p>\u00c0 l&rsquo;inverse des syst\u00e8mes pr\u00e9c\u00e9dents, la file <strong>E<sub>k<\/sub>\/M\/1<\/strong> mod\u00e9lise une configuration o\u00f9 les temps s\u00e9parant deux arriv\u00e9es successives de clients ob\u00e9issent \u00e0 une distribution d&rsquo;Erlang \u00e0 <em>k<\/em> phases, tandis que la dur\u00e9e de service au guichet demeure purement exponentielle. Ce mod\u00e8le intervient fr\u00e9quemment dans les processus manufacturiers ou les syst\u00e8mes logistiques o\u00f9 les pi\u00e8ces ou usagers n&rsquo;arrivent pas de mani\u00e8re totalement erratique, mais ont pr\u00e9alablement franchi une succession de <em>k<\/em> postes de conditionnement ou de contr\u00f4le qui r\u00e9gulent et lissent leur flux.<\/p>\n<p>L&rsquo;analyse stochastique du syst\u00e8me E<sub>k<\/sub>\/M\/1 repose sur l&rsquo;introduction d&rsquo;un processus de Markov semi-r\u00e9g\u00e9n\u00e9ratif. L&rsquo;\u00e9criture des \u00e9quations aux diff\u00e9rences r\u00e9gissant les probabilit\u00e9s d&rsquo;\u00e9tat stationnaires imm\u00e9diatement avant l&rsquo;arriv\u00e9e d&rsquo;un client montre que la distribution g\u00e9om\u00e9trique stationnaire d\u00e9pend d&rsquo;une racine unique <em>r<sub>0<\/sub><\/em> comprise strictement entre 0 et 1, solution de l&rsquo;\u00e9quation implicite caract\u00e9ristique :<\/p>\n<p><em>r = [ &lambda; \/ (&lambda; + &mu;(1 &#8211; r)) ]<sup>k<\/sup><\/em><\/p>\n<p>La pr\u00e9sence du m\u00e9canisme d&rsquo;arriv\u00e9e erlangien att\u00e9nue consid\u00e9rablement l&rsquo;apparition de grappes d&rsquo;arriv\u00e9es (ou <em>bursts<\/em>). Par rapport \u00e0 un processus de Poisson classique, les arriv\u00e9es erlangiennes sont plus espac\u00e9es de fa\u00e7on r\u00e9guli\u00e8re dans le temps, ce qui diminue drastiquement la probabilit\u00e9 qu&rsquo;une surtension instantan\u00e9e sature le serveur d&rsquo;attente, am\u00e9liorant ainsi la r\u00e9silience globale du r\u00e9seau.<\/p>\n<h2>8. Applications en psychologie cognitive et mod\u00e9lisation comportementale<\/h2>\n<h3>8.1 Chronom\u00e9trie mentale et temps de r\u00e9action humaine<\/h3>\n<p>L&rsquo;une des transpositions disciplinaires les plus fascinantes de la distribution d&rsquo;Erlang concerne le champ de la psychologie exp\u00e9rimentale et de la chronom\u00e9trie cognitive. D\u00e8s les ann\u00e9es 1960, des chercheurs en neurosciences cognitives ont cherch\u00e9 \u00e0 formaliser la distribution statistique des <strong>temps de r\u00e9action<\/strong> (TR) humains lors de t\u00e2ches perceptives ou d\u00e9cisionnelles simples et complexes. Les enregistrements chronom\u00e9triques montrent invariablement une courbe asym\u00e9trique \u00e0 queue droite que la loi normale \u00e9choue syst\u00e9matiquement \u00e0 d\u00e9crire.<\/p>\n<p>Dans ce contexte, la distribution d&rsquo;Erlang fournit une analogie th\u00e9orique puissante : l&rsquo;hypoth\u00e8se du traitement s\u00e9quentiel par \u00e9tapes. Selon ce mod\u00e8le neuronal, la r\u00e9ponse motrice finale cons\u00e9cutive \u00e0 l&rsquo;apparition d&rsquo;un stimulus sensoriel n&rsquo;est pas un processus continu et indivisible, mais requiert l&rsquo;ex\u00e9cution ordonn\u00e9e d&rsquo;un nombre entier <em>k<\/em> d&rsquo;\u00e9tapes de traitement neurocognitif ind\u00e9pendantes. Ces \u00e9tapes englobent successivement :<\/p>\n<ul>\n<li>La transduction sensorielle p\u00e9riph\u00e9rique (r\u00e9tinienne ou cochl\u00e9aire).<\/li>\n<li>Le relais thalamique et l&rsquo;activation des aires corticales primaires.<\/li>\n<li>La cat\u00e9gorisation cognitive du stimulus au sein du cortex associatif.<\/li>\n<li>La s\u00e9lection programmatique de la r\u00e9ponse comportementale ad\u00e9quate.<\/li>\n<li>La commande motrice descendante et la latence neuromusculaire p\u00e9riph\u00e9rique.<\/li>\n<\/ul>\n<p>Si chacune de ces micro-\u00e9tapes c\u00e9r\u00e9brales poss\u00e8de une dur\u00e9e al\u00e9atoire distribu\u00e9e selon une loi exponentielle de taux moyen homog\u00e8ne &lambda;, le temps de r\u00e9action global r\u00e9sultant correspond avec une remarquable fid\u00e9lit\u00e9 \u00e0 une distribution d&rsquo;Erlang \u00e0 <em>k<\/em> phases. L&rsquo;ajustement empirique de cette loi sur des cohortes de sujets permet ainsi d&rsquo;inf\u00e9rer statistiquement le nombre probable de relais neuronaux impliqu\u00e9s dans la cha\u00eene de traitement.<\/p>\n<h3>8.2 Mod\u00e9lisation de l&rsquo;attention et des processus d\u00e9cisionnels<\/h3>\n<p>Dans les paradigmes contemporains de prise de d\u00e9cision rapide, comme le mod\u00e8le de diffusion ou les mod\u00e8les d&rsquo;accumulation de preuves stochastiques (<em>race models<\/em>), la distribution d&rsquo;Erlang est fr\u00e9quemment mobilis\u00e9e pour repr\u00e9senter la dynamique temporelle de franchissement des barri\u00e8res attentionnelles. Lorsqu&rsquo;un sujet humain doit discriminer deux alternatives perceptives conflictuelles, son cortex fronto-pari\u00e9tal accumule s\u00e9quentiellement des paquets d&rsquo;informations neuronales (quanta d&rsquo;\u00e9vidence) jusqu&rsquo;\u00e0 d\u00e9passer un seuil critique d&rsquo;activation.<\/p>\n<p>Si l&rsquo;on postule que chaque incr\u00e9ment informationnel \u00e9l\u00e9mentaire survient selon un processus de Poisson ind\u00e9pendant, le d\u00e9lai requis pour collecter le quota de <em>k<\/em> indices n\u00e9cessaires pour arr\u00eater d\u00e9finitivement la d\u00e9cision suit une distribution d&rsquo;Erlang. Cette approche permet de mod\u00e9liser avec \u00e9l\u00e9gance la variabilit\u00e9 intra-individuelle des temps de r\u00e9ponse : sous contrainte temporelle forte ou en situation de fatigue cognitive, le param\u00e8tre <em>k<\/em> refl\u00e8te le niveau d&rsquo;exigence probatoire fix\u00e9 par le sujet, permettant de quantifier le compromis classique entre vitesse d&rsquo;ex\u00e9cution et exactitude de la d\u00e9cision (<em>speed-accuracy tradeoff<\/em>).<\/p>\n<h3>8.3 Psychologie ergonomique et perception de l&rsquo;attente<\/h3>\n<p>Au-del\u00e0 de la chronom\u00e9trie mentale fondamentale, la distribution d&rsquo;Erlang trouve une application op\u00e9rationnelle d\u00e9terminante en psychologie ergonomique et dans la conception de l&rsquo;exp\u00e9rience utilisateur (<em>UX design<\/em>). Des \u00e9tudes comportementales approfondies ont d\u00e9montr\u00e9 que l&rsquo;anxi\u00e9t\u00e9, le stress cognitif et l&rsquo;insatisfaction subjective ressentis par un utilisateur ou un usager lors d&rsquo;une attente ne sont pas corr\u00e9l\u00e9s uniquement \u00e0 la dur\u00e9e brute moyenne de cette attente, mais d\u00e9pendent de fa\u00e7on pr\u00e9pond\u00e9rante de sa <strong>variance per\u00e7ue<\/strong>.<\/p>\n<p>Un service caract\u00e9ris\u00e9 par un temps distribu\u00e9 selon une loi exponentielle (<em>k = 1<\/em>) g\u00e9n\u00e8re un sentiment d&rsquo;impr\u00e9visibilit\u00e9 et d&rsquo;injustice temporelle en raison de sa queue lourde et de ses fluctuations chaotiques. En revanche, un service articul\u00e9 selon une distribution d&rsquo;Erlang de param\u00e8tre <em>k<\/em> \u00e9lev\u00e9 (par exemple <em>k &ge; 8<\/em>) pr\u00e9sente une variance contract\u00e9e qui conf\u00e8re \u00e0 l&rsquo;utilisateur un sentiment de pr\u00e9visibilit\u00e9 rassurante. En ergonomie logicielle et dans l&rsquo;ing\u00e9nierie des services marchands, segmenter une attente globale en sous-phases visibles (barres de progression d\u00e9coup\u00e9es en \u00e9tapes logiques) permet de modifier la perception temporelle de l&rsquo;op\u00e9rateur en alignant l&rsquo;exp\u00e9rience sur la dynamique d&rsquo;une distribution d&rsquo;Erlang ma\u00eetris\u00e9e.<\/p>\n<h2>9. M\u00e9thodes d&rsquo;estimation statistique des param\u00e8tres<\/h2>\n<h3>9.1 Estimation par la m\u00e9thode des moments<\/h3>\n<p>L&rsquo;estimation des param\u00e8tres <em>k<\/em> et &lambda; \u00e0 partir d&rsquo;un \u00e9chantillon empirique d&rsquo;observations ind\u00e9pendantes <em>x<sub>1<\/sub>, x<sub>2<\/sub>, &#8230;, x<sub>n<\/sub><\/em> constitue une probl\u00e9matique statistique classique mais sp\u00e9cifique en raison de la contrainte d&rsquo;int\u00e9gralit\u00e9 pesant sur <em>k<\/em>. La m\u00e9thode des moments repose sur l&rsquo;\u00e9galisation des moments th\u00e9oriques de premier et second ordre avec leurs homologues empiriques calcul\u00e9s sur l&rsquo;\u00e9chantillon :<\/p>\n<p>Moyenne empirique : <em>x\u0304 = (1 \/ n) &sum;<sub>i=1<\/sub><sup>n<\/sup> x<sub>i<\/sub><\/em><\/p>\n<p>Variance empirique : <em>s<sup>2<\/sup> = [ 1 \/ (n &#8211; 1) ] &sum;<sub>i=1<\/sub><sup>n<\/sup> (x<sub>i<\/sub> &#8211; x\u0304)<sup>2<\/sup><\/em><\/p>\n<p>En exploitant les \u00e9quations th\u00e9oriques d\u00e9montr\u00e9es au paragraphe 4, nous \u00e9tablissons le syst\u00e8me d&rsquo;\u00e9galit\u00e9s :<\/p>\n<p><em>E[X] = k \/ &lambda; = x\u0304<\/em> &nbsp;&nbsp; et &nbsp;&nbsp; <em>Var(X) = k \/ &lambda;<sup>2<\/sup> = s<sup>2<\/sup><\/em><\/p>\n<p>En divisant le carr\u00e9 de l&rsquo;esp\u00e9rance par la variance, on isole le param\u00e8tre de forme <em>k<\/em>, puis on en d\u00e9duit le taux &lambda; :<\/p>\n<p><em>k\u0303 = (x\u0304)<sup>2<\/sup> \/ s<sup>2<\/sup><\/em><\/p>\n<p><em>\u03bb\u0303 = x\u0304 \/ s<sup>2<\/sup><\/em><\/p>\n<p>Cependant, dans la quasi-totalit\u00e9 des contextes r\u00e9els, le ratio empirique <em>(x\u0304)<sup>2<\/sup> \/ s<sup>2<\/sup><\/em> donne un nombre r\u00e9el fractionnaire non entier. Pour respecter le cadre formel de la distribution d&rsquo;Erlang, il est indispensable de contraindre l&rsquo;estimateur de <em>k<\/em> \u00e0 l&rsquo;ensemble des entiers naturels. La pratique canonique consiste \u00e0 arrondir la valeur brute \u00e0 l&rsquo;entier le plus proche :<\/p>\n<p><em>k\u0302 = max( 1, round( (x\u0304)<sup>2<\/sup> \/ s<sup>2<\/sup> ) )<\/em><\/p>\n<p>Une fois cet entier <em>k\u0302<\/em> fix\u00e9, on r\u00e9ajuste l&rsquo;estimateur du param\u00e8tre de taux pour conserver l&rsquo;exactitude de la moyenne observ\u00e9e : <em>\u03bb\u0302 = k\u0302 \/ x\u0304<\/em>. Bien que simple et rapide \u00e0 impl\u00e9menter, cette m\u00e9thode pr\u00e9sente un biais pour les \u00e9chantillons de petite taille.<\/p>\n<h3>9.2 Estimation par le maximum de vraisemblance<\/h3>\n<p>La <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Maximum_likelihood_estimation\">m\u00e9thode du maximum de vraisemblance<\/a> (MLE) fournit des estimateurs asymptotiquement efficaces et non biais\u00e9s. Pour un \u00e9chantillon de <em>n<\/em> r\u00e9alisations ind\u00e9pendantes <em>X = {x<sub>1<\/sub>, &#8230;, x<sub>n<\/sub>}<\/em>, la fonction de vraisemblance conjointe s&rsquo;\u00e9crit comme le produit des densit\u00e9s individuelles :<\/p>\n<p><em>L(k, &lambda; | X) = &prod;<sub>i=1<\/sub><sup>n<\/sup> [ &lambda;<sup>k<\/sup> x<sub>i<\/sub><sup>k-1<\/sup> e<sup>-&lambda;x<sub>i<\/sub><\/sup> \/ (k &#8211; 1)! ] = [ &lambda;<sup>nk<\/sup> \/ ((k &#8211; 1)!)<sup>n<\/sup> ] &middot; ( &prod;<sub>i=1<\/sub><sup>n<\/sup> x<sub>i<\/sub> )<sup>k-1<\/sup> &middot; exp( -&lambda; &sum;<sub>i=1<\/sub><sup>n<\/sup> x<sub>i<\/sub> )<\/em><\/p>\n<p>Le passage au logarithme n\u00e9p\u00e9rien permet d&rsquo;obtenir la fonction de log-vraisemblance (<em>ln L<\/em>) :<\/p>\n<p><em>ln L(k, &lambda;) = n &middot; k &middot; ln(&lambda;) &#8211; n &middot; ln((k &#8211; 1)!) + (k &#8211; 1) &sum;<sub>i=1<\/sub><sup>n<\/sup> ln(x<sub>i<\/sub>) &#8211; &lambda; &sum;<sub>i=1<\/sub><sup>n<\/sup> x<sub>i<\/sub><\/em><\/p>\n<p>Pour maximiser cette fonction, on examine en premier lieu la d\u00e9riv\u00e9e partielle par rapport au param\u00e8tre continu &lambda; en supposant <em>k<\/em> temporairement fix\u00e9 :<\/p>\n<p><em>&part; ln L \/ &part;&lambda; = (n &middot; k \/ &lambda;) &#8211; &sum;<sub>i=1<\/sub><sup>n<\/sup> x<sub>i<\/sub> = 0 &rArr; &lambda;\u0302(k) = n &middot; k \/ &sum;<sub>i=1<\/sub><sup>n<\/sup> x<sub>i<\/sub> = k \/ x\u0304<\/em><\/p>\n<p>Ce r\u00e9sultat remarquable montre que pour toute valeur discr\u00e8te de <em>k<\/em>, l&rsquo;estimateur optimal du taux &lambda; est strictement proportionnel \u00e0 <em>k<\/em> divis\u00e9 par la moyenne arithm\u00e9tique. En r\u00e9injectant cette expression conditionnelle <em>\u03bb\u0302(k)<\/em> dans la log-vraisemblance, on obtient la fonction de log-vraisemblance profil\u00e9e ne d\u00e9pendant plus que du param\u00e8tre entier <em>k<\/em> :<\/p>\n<p><em>ln L<sub>prof<\/sub>(k) = n &middot; k &middot; ln(k \/ x\u0304) &#8211; n &middot; ln((k &#8211; 1)!) + (k &#8211; 1) &sum;<sub>i=1<\/sub><sup>n<\/sup> ln(x<sub>i<\/sub>) &#8211; n &middot; k<\/em><\/p>\n<p>En raison de la nature discr\u00e8te de <em>k<\/em> &isin; {1, 2, 3, &#8230;}, la maximisation finale ne peut s&rsquo;effectuer par annulation d&rsquo;un gradient continu. On recourt donc \u00e0 un <strong>algorithme de recherche par balayage s\u00e9quentiel<\/strong> (<em>grid search<\/em> ou recherche lin\u00e9aire) sur les entiers successifs <em>k = 1, 2, 3, &#8230;<\/em> jusqu&rsquo;\u00e0 observer une stricte diminution de la log-vraisemblance profil\u00e9e, garantissant ainsi l&rsquo;obtention du maximum global.<\/p>\n<h3>9.3 Approches bay\u00e9siennes et inf\u00e9rence moderne<\/h3>\n<p>L&rsquo;inf\u00e9rence statistique bay\u00e9sienne aborde l&rsquo;estimation des param\u00e8tres de la distribution d&rsquo;Erlang en traitant <em>k<\/em> et &lambda; comme des variables al\u00e9atoires dot\u00e9es de distributions <em>a priori<\/em>. Pour le param\u00e8tre d&rsquo;intensit\u00e9 &lambda;, la famille conjugu\u00e9e naturelle (\u00e0 <em>k<\/em> fix\u00e9) est constitu\u00e9e par la distribution Gamma :<\/p>\n<p><em>&lambda; ~ Gamma(&alpha;<sub>0<\/sub>, &beta;<sub>0<\/sub>)<\/em><\/p>\n<p>Gr\u00e2ce \u00e0 la propri\u00e9t\u00e9 de conjugaison alg\u00e9brique, la distribution <em>a posteriori<\/em> conditionnelle de &lambda; sachant les donn\u00e9es <em>X<\/em> et sachant <em>k<\/em> est \u00e9galement une distribution Gamma actualis\u00e9e :<\/p>\n<p><em>p(&lambda; | X, k) ~ Gamma( &alpha;<sub>0<\/sub> + n &middot; k, &beta;<sub>0<\/sub> + &sum;<sub>i=1<\/sub><sup>n<\/sup> x<sub>i<\/sub> )<\/em><\/p>\n<p>Pour le param\u00e8tre de forme discret <em>k<\/em>, on lui assigne g\u00e9n\u00e9ralement une loi <em>a priori<\/em> non informative ou discr\u00e8te uniforme born\u00e9e sur un intervalle r\u00e9aliste [1, <em>K<sub>max<\/sub><\/em>], ou encore une loi de Poisson tronqu\u00e9e. L&rsquo;\u00e9chantillonnage de la distribution conjointe <em>a posteriori<\/em> <em>p(k, &lambda; | X)<\/em> ne pouvant \u00eatre r\u00e9solu analytiquement sous forme ferm\u00e9e, on d\u00e9ploie des algorithmes d&rsquo;\u00e9chantillonnage par cha\u00eenes de Markov Monte-Carlo (MCMC), tels que l&rsquo;\u00e9chantillonneur de Gibbs combin\u00e9 \u00e0 un pas de Metropolis-Hastings pour le param\u00e8tre discret <em>k<\/em>.<\/p>\n<p>Cette m\u00e9thodologie bay\u00e9sienne pr\u00e9sente des atouts consid\u00e9rables lorsque l&rsquo;on manipule des donn\u00e9es \u00e9parses ou des \u00e9chantillons de tr\u00e8s faible effectif (par exemple lors de l&rsquo;\u00e9tude d&rsquo;accidents industriels rares), car elle int\u00e8gre l&rsquo;incertitude \u00e9pist\u00e9mique sur les param\u00e8tres et permet d&rsquo;\u00e9viter les surajustements fr\u00e9quents des m\u00e9thodes fr\u00e9quentistes.<\/p>\n<h2>10. Simulation num\u00e9rique et g\u00e9n\u00e9ration de variables al\u00e9atoires<\/h2>\n<h3>10.1 M\u00e9thode par sommation de variables exponentielles<\/h3>\n<p>La g\u00e9n\u00e9ration num\u00e9rique de nombres pseudo-al\u00e9atoires distribu\u00e9s selon une loi d&rsquo;Erlang s&rsquo;appuie directement sur son th\u00e9or\u00e8me constructif fondamental. Pour simuler une variable al\u00e9atoire exponentielle de param\u00e8tre &lambda;, l&rsquo;algorithme standard utilise la m\u00e9thode de la transform\u00e9e inverse. Si <em>U<\/em> est une variable al\u00e9atoire uniform\u00e9ment distribu\u00e9e sur l&rsquo;intervalle ouvert ]0, 1[, alors la variable <em>Y = -ln(U) \/ &lambda;<\/em> suit exactement une distribution exponentielle <em>Exp(&lambda;)<\/em>.<\/p>\n<p>Pour engendrer une r\u00e9alisation stochastique <em>X<\/em> d&rsquo;une distribution d&rsquo;Erlang de param\u00e8tres <em>k<\/em> et &lambda;, il suffit donc d&rsquo;additionner <em>k<\/em> r\u00e9alisations exponentielles ind\u00e9pendantes :<\/p>\n<p><em>X = &sum;<sub>i=1<\/sub><sup>k<\/sup> [ -ln(U<sub>i<\/sub>) \/ &lambda; ] = &#8211; (1 \/ &lambda;) &middot; &sum;<sub>i=1<\/sub><sup>k<\/sup> ln(U<sub>i<\/sub>)<\/em><\/p>\n<p>En exploitant les propri\u00e9t\u00e9s fondamentales du logarithme selon lesquelles la somme de logarithmes est \u00e9gale au logarithme du produit, on optimise la formulation math\u00e9matique pour n&rsquo;\u00e9valuer qu&rsquo;un seul appel \u00e0 la fonction transcendante <em>ln()<\/em> :<\/p>\n<p><em>X = &#8211; (1 \/ &lambda;) &middot; ln( &prod;<sub>i=1<\/sub><sup>k<\/sup> U<sub>i<\/sub> )<\/em><\/p>\n<p>Cet algorithme direct de multiplication de <em>k<\/em> variables uniformes ind\u00e9pendantes est \u00e9l\u00e9gant, exact et d&rsquo;une remarquable simplicit\u00e9 logique. N\u00e9anmoins, sa complexit\u00e9 temporelle est proportionnelle \u00e0 <em>O(k)<\/em>, ce qui signifie que le temps de calcul augmente de mani\u00e8re lin\u00e9aire avec <em>k<\/em>. En outre, pour des valeurs de <em>k<\/em> particuli\u00e8rement \u00e9lev\u00e9es (par exemple <em>k &gt; 100<\/em>), le produit r\u00e9p\u00e9t\u00e9 de nombres strictement compris entre 0 et 1 peut g\u00e9n\u00e9rer un sous-d\u00e9passement arithm\u00e9tique (<em>underflow<\/em>), imposant de scinder le calcul ou de basculer vers des g\u00e9n\u00e9rateurs g\u00e9n\u00e9raux de lois Gamma bas\u00e9s sur la m\u00e9thode du rejet de Marsaglia et Tsang.<\/p>\n<h3>10.2 Impl\u00e9mentation algorithmique dans les environnements scientifiques<\/h3>\n<p>Les environnements modernes de calcul scientifique int\u00e8grent des modules hautement optimis\u00e9s pour manipuler et g\u00e9n\u00e9rer des variables de distribution d&rsquo;Erlang. Dans l&rsquo;\u00e9cosyst\u00e8me <strong>Python<\/strong>, la biblioth\u00e8que scientifique <em>SciPy<\/em> propose la classe sp\u00e9cialis\u00e9e <code>scipy.stats.erlang<\/code>. Il convient de souligner que dans cette impl\u00e9mentation, la classe h\u00e9rite directement de la loi Gamma, avec l&rsquo;exigence explicite que le param\u00e8tre de forme soit un entier strict. Le param\u00e9trage utilise l&rsquo;argument <code>a = k<\/code> pour la forme et <code>scale = 1 \/ &lambda;<\/code> pour l&rsquo;\u00e9chelle temporelle.<\/p>\n<p>Dans l&rsquo;environnement de programmation <strong>R<\/strong>, tr\u00e8s pris\u00e9 en bio-statistique et en actuariat, la loi d&rsquo;Erlang ne dispose pas syst\u00e9matiquement d&rsquo;un pr\u00e9fixe d\u00e9di\u00e9 propre (comme <code>rerlang<\/code>) mais s&rsquo;ex\u00e9cute nativement et rigoureusement via les fonctions canoniques de la famille Gamma : <code>dgamma(x, shape = k, rate = lambda)<\/code> pour la densit\u00e9 et <code>rgamma(n, shape = k, scale = 1 \/ lambda)<\/code> pour la g\u00e9n\u00e9ration d&rsquo;\u00e9chantillons stochastiques, en contraignant explicitement la valeur attribu\u00e9e au param\u00e8tre <code>shape<\/code> \u00e0 un entier naturel.<\/p>\n<p>Sous la plateforme logicielle <strong>MATLAB<\/strong>, la fonction <code>gamrnd(k, 1\/lambda)<\/code> est couramment sollicit\u00e9e pour produire des r\u00e9alisations stochastiques vectoris\u00e9es. Sur le plan des benchmarks de performance, les g\u00e9n\u00e9rateurs modernes impl\u00e9mentent des optimisations de bas niveau en C ou Fortran exploitant des registres vectoriels SIMD, permettant la g\u00e9n\u00e9ration simultan\u00e9e de millions de r\u00e9alisations erlangiennes en quelques fractions de seconde, m\u00eame pour des structures de r\u00e9seaux complexes.<\/p>\n<h3>10.3 Protocoles de validation et tests d&rsquo;ad\u00e9quation<\/h3>\n<p>Une fois les variables stochastiques g\u00e9n\u00e9r\u00e9es ou ajust\u00e9es sur des jeux de donn\u00e9es r\u00e9els, il est scientifiquement imp\u00e9ratif de valider la conformit\u00e9 distributionnelle par des tests d&rsquo;ad\u00e9quation statistique rigoureux. Le test non param\u00e9trique de <strong>Kolmogorov-Smirnov<\/strong> (K-S) est fr\u00e9quemment d\u00e9ploy\u00e9 pour comparer la fonction de r\u00e9partition empirique <em>F<sub>n<\/sub>(x)<\/em> \u00e0 la fonction de r\u00e9partition th\u00e9orique <em>F(x; k, &lambda;)<\/em>. La statistique de test <em>D<sub>n<\/sub><\/em> est d\u00e9finie par le supremum absolu des \u00e9carts :<\/p>\n<p><em>D<sub>n<\/sub> = sup<sub>x<\/sub> | F<sub>n<\/sub>(x) &#8211; F(x; k, &lambda;) |<\/em><\/p>\n<p>Si la valeur critique calcul\u00e9e est inf\u00e9rieure au seuil tabul\u00e9 pour un niveau de signification &alpha; donn\u00e9 (habituellement &alpha; = 0,05), l&rsquo;hypoth\u00e8se nulle d&rsquo;ad\u00e9quation \u00e0 la loi d&rsquo;Erlang ne peut \u00eatre rejet\u00e9e. Le <strong>test d&rsquo;ad\u00e9quation du Chi-deux de Pearson<\/strong> offre une m\u00e9thode alternative particuli\u00e8rement robuste pour les effectifs importants, consistant \u00e0 partitionner le domaine temporel en classes disjointes et \u00e0 comparer les fr\u00e9quences observ\u00e9es aux effectifs attendus calcul\u00e9s \u00e0 partir de la fonction de r\u00e9partition de la loi d&rsquo;Erlang.<\/p>\n<p>Enfin, sur le plan heuristique et qualitatif, le protocole diagnostique repose sur le trac\u00e9 du <strong>diagramme quantile-quantile<\/strong> (Q-Q plot). Ce graphique projette les quantiles empiriques de l&rsquo;\u00e9chantillon observ\u00e9 en fonction des quantiles th\u00e9oriques calcul\u00e9s pour la loi d&rsquo;Erlang de param\u00e8tres estim\u00e9s. Un alignement lin\u00e9aire rectiligne des points le long de la bissectrice confirme l&rsquo;ad\u00e9quation du mod\u00e8le, tandis que toute d\u00e9viation syst\u00e9matique au niveau des extr\u00e9mit\u00e9s graphiques signale des queues de distribution divergentes n\u00e9cessitant une r\u00e9vision param\u00e9trique.<\/p>\n<h2>11. Applications interdisciplinaires avanc\u00e9es<\/h2>\n<h3>11.1 Ing\u00e9nierie moderne des r\u00e9seaux informatiques et t\u00e9l\u00e9coms<\/h3>\n<p>Bien que n\u00e9e dans les infrastructures t\u00e9l\u00e9phoniques \u00e0 circuits commut\u00e9s du XXe si\u00e8cle, la distribution d&rsquo;Erlang s&rsquo;est impos\u00e9e comme un concept fondamental au sein des architectures num\u00e9riques contemporaines. Dans les r\u00e9seaux de commutation par paquets, la transmission d&rsquo;un datagramme IP \u00e0 travers une cha\u00eene s\u00e9quentielle de commutateurs et de routeurs Internet implique la travers\u00e9e ordonn\u00e9e de plusieurs tampons m\u00e9moire (<em>buffers<\/em>) et op\u00e9rations de routage. Chaque passage dans un n\u0153ud interm\u00e9diaire g\u00e9n\u00e9rant un temps de latence, la distribution du <strong>d\u00e9lai d&rsquo;acheminement de bout en bout<\/strong> (<em>jitter<\/em> et latence totale) s&rsquo;ajuste parfaitement \u00e0 une loi d&rsquo;Erlang d&rsquo;ordre <em>k<\/em>, o\u00f9 <em>k<\/em> mat\u00e9rialise le nombre de sauts de routage (<em>hop count<\/em>).<\/p>\n<p>Dans le domaine des t\u00e9l\u00e9communications mobiles de cinqui\u00e8me et sixi\u00e8me g\u00e9n\u00e9ration (5G et 6G), la gestion dynamique de la bande passante sous contrainte stricte de <em>Quality of Service<\/em> (QoS) et d&rsquo;ultra-fiabilit\u00e9 (protocoles URLLC pour <em>Ultra-Reliable and Low-Latency Communications<\/em>) exploite des mod\u00e8les d&rsquo;Erlang pour dimensionner le partitionnement virtuel des r\u00e9seaux (<em>network slicing<\/em>). Les algorithmes d&rsquo;admission d&rsquo;appels et d&rsquo;allocation spectrale calculent en temps r\u00e9el les probabilit\u00e9s de blocage instantan\u00e9 et de mise en m\u00e9moire tampon des trames de donn\u00e9es \u00e0 l&rsquo;aide des structures erlangiennes \u00e9tendues.<\/p>\n<p>De m\u00eame, au sein des architectures en grappe (<em>clusters<\/em> de serveurs, architectures micro-services et informatique en nuage), le traitement d&rsquo;une requ\u00eate web complexe est couramment d\u00e9coup\u00e9 en un pipeline s\u00e9quentiel d&rsquo;appels RPC (<em>Remote Procedure Calls<\/em>) non parall\u00e8les. La distribution d&rsquo;Erlang sert alors \u00e0 caract\u00e9riser le temps global de compl\u00e9tion des requ\u00eates distribu\u00e9es et \u00e0 d\u00e9tecter les micro-pannes induites par la d\u00e9rive temporelle de serveurs applicatifs sp\u00e9cifiques.<\/p>\n<h3>11.2 Biologie math\u00e9matique et mod\u00e9lisation \u00e9pid\u00e9miologique<\/h3>\n<p>En \u00e9pid\u00e9miologie th\u00e9orique et en sant\u00e9 publique computationnelle, la distribution d&rsquo;Erlang joue un r\u00f4le strat\u00e9gique pour d\u00e9passer les simplifications irr\u00e9alistes des mod\u00e8les compartimentaux classiques de type <strong>SEIR<\/strong> (Sensible &ndash; Expos\u00e9 &ndash; Infect\u00e9 &ndash; R\u00e9tabli). Dans le mod\u00e8le historique de Kermack et McKendrick, le transfert des individus d&rsquo;un compartiment \u00e0 l&rsquo;autre repose sur des \u00e9quations diff\u00e9rentielles lin\u00e9aires impliquant des temps de s\u00e9jour distribu\u00e9s selon une loi exponentielle.<\/p>\n<p>Or, en r\u00e9alit\u00e9 biologique, la p\u00e9riode d&rsquo;incubation d&rsquo;un agent pathog\u00e8ne (virus ou bact\u00e9rie) et la dur\u00e9e de la p\u00e9riode infectieuse ne pr\u00e9sentent jamais un taux d&rsquo;\u00e9limination maximal au premier jour de la contamination : elles exigent une phase de prolif\u00e9ration intra-h\u00f4te minimale avant l&rsquo;apparition des sympt\u00f4mes ou de la contagion. En rempla\u00e7ant un compartiment d&rsquo;exposition unique par une cha\u00eene de <em>k<\/em> sous-compartiments cons\u00e9cutifs <em>E<sub>1<\/sub>, E<sub>2<\/sub>, &#8230;, E<sub>k<\/sub><\/em>, le temps d&rsquo;incubation global suit automatiquement une distribution d&rsquo;Erlang \u00e0 <em>k<\/em> phases :<\/p>\n<p><em>dE<sub>1<\/sub>\/dt = &beta; S I &#8211; k &sigma; E<sub>1<\/sub><\/em><\/p>\n<p><em>dE<sub>j<\/sub>\/dt = k &sigma; E<sub>j-1<\/sub> &#8211; k &sigma; E<sub>j<\/sub> &nbsp;&nbsp; (pour j = 2, &#8230;, k)<\/em><\/p>\n<p>Cette approche, valid\u00e9e exp\u00e9rimentalement lors de l&rsquo;analyse des \u00e9pid\u00e9mies d&rsquo;Ebola, de la grippe H1N1 et de la COVID-19, permet d&rsquo;ajuster avec une fid\u00e9lit\u00e9 nettement sup\u00e9rieure le temps d&rsquo;incubation moyen et son \u00e9talement statistique sans quitter le cadre commode des syst\u00e8mes d&rsquo;\u00e9quations diff\u00e9rentielles ordinaires. Les projections sanitaires et l&rsquo;\u00e9valuation de l&rsquo;efficacit\u00e9 des mesures de mise en quarantaine gagnent ainsi en rigueur et en repr\u00e9sentativit\u00e9 terrain.<\/p>\n<p>\u00c0 l&rsquo;\u00e9chelle cellulaire et mol\u00e9culaire, la distribution d&rsquo;Erlang est \u00e9galement appliqu\u00e9e \u00e0 la mod\u00e9lisation de la cin\u00e9tique enzymatique et du cycle de division cellulaire. Le passage d&rsquo;une cellule par les phases successives de l&rsquo;interphase (G1, S, G2) avant la mitose M n\u00e9cessite l&rsquo;accumulation d&rsquo;un ensemble de d\u00e9clencheurs biochimiques, chaque transition enzymatique \u00e9l\u00e9mentaire agissant comme une \u00e9tape exponentielle au sein d&rsquo;une distribution d&rsquo;Erlang macroscopique.<\/p>\n<h3>11.3 Gestion des flux de patients en milieu hospitalier<\/h3>\n<p>L&rsquo;optimisation des parcours de soins et de l&rsquo;ing\u00e9nierie hospitali\u00e8re constitue un terrain d&rsquo;\u00e9lection pour la distribution d&rsquo;Erlang. Les flux de patients se pr\u00e9sentant dans les services d&rsquo;urgence ou programm\u00e9s pour des interventions chirurgicales au sein des blocs op\u00e9ratoires traversent des fili\u00e8res hautement norm\u00e9es comprenant un protocole m\u00e9dical strict d\u00e9coup\u00e9 en phases :<\/p>\n<ul>\n<li>Accueil, enregistrement m\u00e9dico-administratif et tri initial par l&rsquo;infirmi\u00e8re organisatrice (IOA).<\/li>\n<li>Examen m\u00e9dical et anamn\u00e8se initiale par le m\u00e9decin urgentiste.<\/li>\n<li>Examens compl\u00e9mentaires d&rsquo;imagerie m\u00e9dicale (scanner, radiographie, \u00e9chographie) et analyses biologiques.<\/li>\n<li>Intervention th\u00e9rapeutique, s\u00e9dation ou pose de pl\u00e2tre.<\/li>\n<li>P\u00e9riode d&rsquo;observation clinique et organisation de la sortie ou du transfert vers un lit d&rsquo;hospitalisation.<\/li>\n<\/ul>\n<p>La dur\u00e9e totale de prise en charge d&rsquo;un patient pr\u00e9sente par cons\u00e9quent une variance r\u00e9duite par rapport \u00e0 un simple processus exponentiel et s&rsquo;ajuste avec une pertinence math\u00e9matique \u00e9prouv\u00e9e \u00e0 une loi d&rsquo;Erlang \u00e0 plusieurs phases. La mod\u00e9lisation conjointe de ces flux via des r\u00e9seaux de files d&rsquo;attente \u00e0 serveurs multiples (mod\u00e8les d&rsquo;Erlang avec files d&rsquo;attente \u00e9tendues) permet aux gestionnaires de sant\u00e9 de calculer le taux de rotation optimal des brancards, de minimiser les temps de souffrance des patients en salle d&rsquo;attente et de dimensionner de fa\u00e7on optimale les plannings du personnel m\u00e9dical soignant face aux variations de l&rsquo;afflux quotidien.<\/p>\n<h2>12. Limites m\u00e9thodologiques et distributions alternatives<\/h2>\n<h3>12.1 Contraintes li\u00e9es \u00e0 la rigidit\u00e9 du param\u00e8tre entier<\/h3>\n<p>Malgr\u00e9 sa robustesse conceptuelle et son \u00e9l\u00e9gance analytique, la distribution d&rsquo;Erlang souffre de limitations m\u00e9thodologiques intrins\u00e8ques qu&rsquo;un statisticien se doit de reconna\u00eetre. La premi\u00e8re contrainte r\u00e9side dans la discontinuit\u00e9 stricte du param\u00e8tre de forme <em>k<\/em> &isin; &amp;naturals;<sup>*<\/sup>. Alors que la variabilit\u00e9 d&rsquo;un jeu de donn\u00e9es empiriques issues du monde r\u00e9el peut sugg\u00e9rer un param\u00e8tre de forme optimal de <em>k = 2,43<\/em>, la loi d&rsquo;Erlang oblige l&rsquo;analyste \u00e0 un arbitrage artificiel en choisissant soit <em>k = 2<\/em>, soit <em>k = 3<\/em>.<\/p>\n<p>Cette discontinuit\u00e9 pose un verrou math\u00e9matique majeur lors de la mise en \u0153uvre d&rsquo;algorithmes d&rsquo;optimisation bas\u00e9s sur les gradients (comme la descente de gradient stochastique ou les algorithmes de quasi-Newton de type BFGS). L&rsquo;absence de d\u00e9rivabilit\u00e9 de la fonction de vraisemblance par rapport \u00e0 l&rsquo;entier <em>k<\/em> interdit toute optimisation conjointe diff\u00e9rentiable directe, compliquant l&rsquo;int\u00e9gration de la distribution d&rsquo;Erlang au sein d&rsquo;architectures d&rsquo;apprentissage automatique moderne ou de r\u00e9seaux de neurones probabilistes.<\/p>\n<p>De surcro\u00eet, le couplage structurel impos\u00e9 par la distribution entre l&rsquo;esp\u00e9rance <em>E[X] = k \/ &lambda;<\/em> et la variance <em>Var(X) = k \/ &lambda;<sup>2<\/sup><\/em> restreint le champ de ses applications possibles. Dans la loi d&rsquo;Erlang, le coefficient de variation au carr\u00e9 vaut rigoureusement <em>c<sub>v<\/sub><sup>2<\/sup> = 1 \/ k<\/em>, ce qui signifie que le ratio de la variance sur le carr\u00e9 de la moyenne ne peut prendre que des fractions discr\u00e8tes bien pr\u00e9cises : 1, 1\/2, 1\/3, 1\/4, etc. D\u00e8s lors qu&rsquo;un syst\u00e8me physique pr\u00e9sente un comportement sur-dispers\u00e9 (<em>overdispersion<\/em>), caract\u00e9ris\u00e9 par un coefficient de variation sup\u00e9rieur \u00e0 un (<em>c<sub>v<\/sub> &gt; 1<\/em>), la distribution d&rsquo;Erlang est structurellement incapable de mod\u00e9liser le ph\u00e9nom\u00e8ne.<\/p>\n<h3>12.2 Comparaison avec la loi Gamma g\u00e9n\u00e9rale et la loi de Weibull<\/h3>\n<p>Face aux contraintes inh\u00e9rentes \u00e0 la loi d&rsquo;Erlang, les statisticiens recourent fr\u00e9quemment \u00e0 deux familles param\u00e9triques concurrentes : la distribution Gamma g\u00e9n\u00e9rale et la <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Weibull_distribution\">distribution de Weibull<\/a>. La distribution Gamma, en s&rsquo;affranchissant de l&rsquo;obligation d&rsquo;un param\u00e8tre entier, \u00e9limine l&rsquo;art\u00e9fact de troncature et autorise un lissage parfait des donn\u00e9es tout en autorisant des m\u00e9thodes d&rsquo;ajustement diff\u00e9rentiables globales.<\/p>\n<p>La distribution de Weibull, quant \u00e0 elle, s&rsquo;appuie sur une formulation analytique diff\u00e9rente de sa fonction de r\u00e9partition :<\/p>\n<p><em>F<sub>Weibull<\/sub>(x; k, &lambda;) = 1 &#8211; e<sup>-(x \/ &lambda;)<sup>k<\/sup><\/sup><\/em><\/p>\n<p>La force de la loi de Weibull r\u00e9side dans la flexibilit\u00e9 exceptionnelle de sa fonction de taux de d\u00e9faillance instantan\u00e9 <em>h(x)<\/em>. Alors que le taux de d\u00e9faillance de la distribution d&rsquo;Erlang est contraint de d\u00e9buter \u00e0 0 et d&rsquo;atteindre asymptotiquement un palier constant &lambda;, la loi de Weibull permet de mod\u00e9liser des taux de d\u00e9faillance en loi de puissance pure : strictement d\u00e9croissants pour <em>k &lt; 1<\/em> (mortalit\u00e9 infantile pr\u00e9coce), strictement constants pour <em>k = 1<\/em>, ou croissant de mani\u00e8re exponentielle sans plafonnement pour <em>k &gt; 1<\/em> (vieillissement acc\u00e9l\u00e9r\u00e9 ou fatigue m\u00e9canique structurelle).<\/p>\n<p>L&rsquo;arbitrage entre ces lois repose donc sur un compromis explicite : si le ph\u00e9nom\u00e8ne \u00e9tudi\u00e9 proc\u00e8de manifestement d&rsquo;une cha\u00eene causale d&rsquo;\u00e9v\u00e9nements discrets successifs et doit \u00eatre incorpor\u00e9 dans une cha\u00eene de Markov, la distribution d&rsquo;Erlang s&rsquo;impose par sa coh\u00e9rence m\u00e9canique. En revanche, pour des \u00e9tudes purement ph\u00e9nomenologiques de survie de mat\u00e9riaux o\u00f9 le comportement aux extr\u00eames pr\u00e9vaut, la distribution de Weibull surpasse fr\u00e9quemment le mod\u00e8le erlangien.<\/p>\n<h3>12.3 Crit\u00e8res rigoureux de s\u00e9lection de mod\u00e8le<\/h3>\n<p>Afin de d\u00e9terminer scientifiquement si l&rsquo;utilisation d&rsquo;une distribution d&rsquo;Erlang est pr\u00e9f\u00e9rable \u00e0 celle d&rsquo;une loi Gamma g\u00e9n\u00e9rale, d&rsquo;une loi de Weibull ou d&rsquo;une loi log-normale, le statisticien applique des protocoles d&rsquo;\u00e9valuation formels fond\u00e9s sur la th\u00e9orie de l&rsquo;information. Les deux m\u00e9triques les plus couramment mobilis\u00e9es sont le <strong>crit\u00e8re d&rsquo;information d&rsquo;Akaike<\/strong> (AIC) et le <strong>crit\u00e8re d&rsquo;information bay\u00e9sien<\/strong> (BIC) :<\/p>\n<p><em>AIC = 2p &#8211; 2 ln(L\u0302)<\/em><\/p>\n<p><em>BIC = p &middot; ln(n) &#8211; 2 ln(L\u0302)<\/em><\/p>\n<p>Dans ces expressions math\u00e9matiques, <em>p<\/em> repr\u00e9sente le nombre de param\u00e8tres libres estim\u00e9s dans le mod\u00e8le, <em>n<\/em> la taille effective de l&rsquo;\u00e9chantillon empirique, et <em>ln(L\u0302)<\/em> la valeur maximale de la log-vraisemblance obtenue \u00e0 convergence. Ces crit\u00e8res op\u00e8rent un arbitrage quantitatif objectif entre la qualit\u00e9 de l&rsquo;ajustement statistique sur les donn\u00e9es (mesur\u00e9e par la vraisemblance) et la complexit\u00e9 structurelle du mod\u00e8le (p\u00e9nalis\u00e9e par le nombre de param\u00e8tres <em>p<\/em>).<\/p>\n<p>Dans le cas d&rsquo;une comparaison entre une distribution d&rsquo;Erlang et une distribution Gamma continue, la loi d&rsquo;Erlang b\u00e9n\u00e9ficie souvent d&rsquo;un traitement avantageux : si la valeur estim\u00e9e du param\u00e8tre continu de la loi Gamma est statistiquement indiscernable d&rsquo;un entier naturel, le principe de parcimonie scientifique (ou rasoir d&rsquo;Ockham) conduit \u00e0 retenir le mod\u00e8le d&rsquo;Erlang en raison de sa plus grande simplicit\u00e9 th\u00e9orique et de son interpr\u00e9tabilit\u00e9 physique directe dans l&rsquo;analyse des syst\u00e8mes complexes.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<p>Aalen, O. O., &amp; Gjessing, H. K. (2001). Understanding the shape of the hazard rate: A process point of view. <em>Statistical Science<\/em>, 16(1), 1-14. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/ss\/1009213285\">https:\/\/doi.org\/10.1214\/ss\/1009213285<\/a><\/p>\n<p>Asmussen, S. (2003). <em>Applied probability and queues<\/em> (2nd ed.). Springer. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/b97237\">https:\/\/doi.org\/10.1007\/b97237<\/a><\/p>\n<p>Brockmeyer, E., Halstr\u00f8m, H. L., &amp; Jensen, A. (1948). <em>The life and works of A. K. Erlang<\/em>. Copenhagen Telephone Company.<\/p>\n<p>Cox, D. R. (1955). A use of complex probabilities in the theory of stochastic processes. <em>Mathematical Proceedings of the Cambridge Philosophical Society<\/em>, 51(2), 313-319. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1017\/S030500410003022X\">https:\/\/doi.org\/10.1017\/S030500410003022X<\/a><\/p>\n<p>Erlang, A. K. (1909). The theory of probabilities and telephone conversations. <em>Nyt Tidsskrift for Matematik B<\/em>, 20, 33-39.<\/p>\n<p>Erlang, A. K. (1917). Solution of some problems in the theory of probabilities of significance in automatic telephone exchanges. <em>Elektroteknikeren<\/em>, 13, 5-13.<\/p>\n<p>Gross, D., Shortle, J. F., Thompson, J. M., &amp; Harris, C. M. (2008). <em>Fundamentals of queueing theory<\/em> (4th ed.). John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/9780470377048\">https:\/\/doi.org\/10.1002\/9780470377048<\/a><\/p>\n<p>Johnson, N. L., Kemp, A. W., &amp; Kotz, S. (2005). <em>Univariate discrete distributions<\/em> (3rd ed.). John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/0471715816\">https:\/\/doi.org\/10.1002\/0471715816<\/a><\/p>\n<p>Kendall, D. G. (1953). Stochastic processes occurring in the theory of queues and their analysis by the method of the imbedded Markov chain. <em>The Annals of Mathematical Statistics<\/em>, 24(3), 338-354. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aoms\/1177728975\">https:\/\/doi.org\/10.1214\/aoms\/1177728975<\/a><\/p>\n<p>Kleinrock, L. (1975). <em>Queueing systems, Volume 1: Theory<\/em>. John Wiley &amp; Sons.<\/p>\n<p>Luce, R. D. (1986). <em>Response times: Their role in inferring elementary mental organization<\/em>. Oxford University Press.<\/p>\n<p>Nelson, R. (1995). <em>Probability, stochastic processes, and queueing theory: The mathematics of computer performance modeling<\/em>. Springer-Verlag. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/978-1-4757-2430-1\">https:\/\/doi.org\/10.1007\/978-1-4757-2430-1<\/a><\/p>\n<p>Ross, S. M. (2014). <em>Introduction to probability models<\/em> (11th ed.). Academic Press. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1016\/C2012-0-06283-0\">https:\/\/doi.org\/10.1016\/C2012-0-06283-0<\/a><\/p>\n<p>Tijms, H. C. (2003). <em>A first course in stochastic models<\/em>. John Wiley &amp; Sons. <a target=\"_blank\" rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1002\/047001363X\">https:\/\/doi.org\/10.1002\/047001363X<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<p>D\u00e9couvrez la distribution d&rsquo;Erlang : fondements math\u00e9matiques, propri\u00e9t\u00e9s statistiques, liens avec d&rsquo;autres lois et applications en th\u00e9orie des files d&rsquo;attente.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[647,359],"tags":[4269,8693,5133,4452,5426],"class_list":["post-7220","psych_stat","type-psych_stat","status-publish","hentry","category-psychologie-cognitive","category-statistiques-appliquees","tag-chronometrie-mentale","tag-distribution-derlang","tag-loi-de-probabilite","tag-modelisation-stochastique","tag-theorie-des-files-dattente"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7220","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=7220"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7220\/revisions"}],"predecessor-version":[{"id":7222,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/7220\/revisions\/7222"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=7220"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=7220"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=7220"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}