Articulo de referencia

Distribución hipoexponencial

\\lambda_{1},\\dots,\\lambda_{k} > 0\\, rates ([[real number|real]])"},"support":{"wt":" x \\in [0; \\infty)\\! "},"pdf":{"wt":"Expressed as a [[phase-type distribution]] -\\bol...

En teoría de la probabilidad, la distribución hipoexponencial o distribución de Erlang generalizada es una distribución continua que se utiliza en campos similares a la distribución de Erlang, como la teoría de colas , la ingeniería de tráfico telefónico y, en general, en procesos estocásticos . Se denomina distribución hipoexponencial porque tiene un coeficiente de variación menor que uno, en comparación con la distribución hiperexponencial , cuyo coeficiente de variación es mayor que uno, y la distribución exponencial, cuyo coeficiente de variación es igual a uno.

Descripción general

La distribución de Erlang es una serie de k distribuciones exponenciales, todas con tasaλ{\displaystyle \lambda }La hipoexponencial es una serie de k distribuciones exponenciales, cada una con su propia tasaλi{\displaystyle \lambda _{i}}, la tasa de laith{\displaystyle i^{th}}distribución exponencial. Si tenemos k variables aleatorias exponenciales distribuidas independientementeincógnitai{\displaystyle X_{i}}, entonces la variable aleatoria,

incógnita=i=1kincógnitai{\displaystyle X=\sum _{i=1}^{k}X_{i}}

tiene una distribución hipoexponencial. La distribución hipoexponencial tiene un coeficiente de variación mínimo de1/k{\displaystyle 1/k}.

Relación con la distribución de tipo fase

Como resultado de la definición, es más fácil considerar esta distribución como un caso especial de la distribución de tipo fase . [ 2 ] La distribución de tipo fase es el tiempo de absorción de un proceso de Markov de estado finito . Si tenemos un proceso de k+1 estados, donde los primeros k estados son transitorios y el estado k+1 es un estado absorbente, entonces la distribución del tiempo desde el inicio del proceso hasta que se alcanza el estado absorbente es de tipo fase. Esto se convierte en la hipoexponencial si comenzamos en el primer 1 y nos movemos sin saltos del estado i al i+1 con una tasaλi{\displaystyle \lambda _{i}}hasta que el estado k transiciona con tasaλk{\displaystyle \lambda _{k}}al estado absorbente k+1 . Esto se puede escribir en forma de matriz subgeneradora,

[λ1λ10000λ2λ20000λk2λk20000λk1λk10000λk].{\displaystyle \left[{\begin{matrix}-\lambda _{1}&\lambda _{1}&0&\dots &0&0\\0&-\lambda _{2}&\lambda _{2}&\ddots &0&0\\\vdots &\ddots &\ddots &\ddots &\ddots &\vdots \\0&0&\ddots &-\lambda _{k-2}&\lambda _{k-2}&0\\0&0&\dots &0&-\lambda _{k-1}&\lambda _{k-1}\\0&0&\dots &0&0&-\lambda _{k}\end{matrix}}\right]\;.}

Para simplificar, denotemos la matriz anterior como la matrizΘΘ(λ1,,λk){\displaystyle \Theta \equiv \Theta (\lambda _{1},\dots ,\lambda _{k})}. Si la probabilidad de comenzar en cada uno de los k estados es

α=(1,0,,0){\displaystyle {\boldsymbol {\alpha }}=(1,0,\dots ,0)}

entoncesHypago(λ1,,λk)=PAGH(α,Θ).{\displaystyle Hypo(\lambda _{1},\dots ,\lambda _{k})=PH({\boldsymbol {\alpha }},\Theta ).}

caso de dos parámetros

Donde la distribución tiene dos parámetros (λ1λ2{\displaystyle \lambda _{1}\neq \lambda _{2}}) Las formas explícitas de las funciones de probabilidad y las estadísticas asociadas son: [ 3 ]

CDF:F(incógnita)=1λ2λ2λ1miλ1incógnitaλ1λ1λ2miλ2incógnita{\displaystyle F(x)=1-{\frac {\lambda _{2}}{\lambda _{2}-\lambda _{1}}}e^{-\lambda _{1}x}-{\frac {\lambda _{1}}{\lambda _{1}-\lambda _{2}}}e^{-\lambda _{2}x}}

PDF:F(incógnita)=λ1λ2λ1λ2(miincógnitaλ2miincógnitaλ1){\displaystyle f(x)={\frac {\lambda _{1}\lambda _{2}}{\lambda _{1}-\lambda _{2}}}(e^{-x\lambda _{2}}-e^{-x\lambda _{1}})}

Significar:1λ1+1λ2{\displaystyle {\frac {1}{\lambda _{1}}}+{\frac {1}{\lambda _{2}}}}

Diferencia:1λ12+1λ22{\displaystyle {\frac {1}{\lambda _{1}^{2}}}+{\frac {1}{\lambda _{2}^{2}}}}

Coeficiente de variación:λ12+λ22λ1+λ2{\displaystyle {\frac {\sqrt {\lambda _{1}^{2}+\lambda _{2}^{2}}}{\lambda _{1}+\lambda _{2}}}}

El coeficiente de variación siempre es menor que 1.

Dado el promedio de la muestra (incógnita¯{\displaystyle {\bar {x}}}) y coeficiente de variación de la muestra (do{\displaystyle c}), los parámetrosλ1{\displaystyle \lambda _{1}}yλ2{\displaystyle \lambda _{2}}se puede estimar de la siguiente manera:

λ1=2incógnita¯[1+1+2(do21)]1{\displaystyle \lambda _{1}={\frac {2}{\bar {x}}}\left[1+{\sqrt {1+2(c^{2}-1)}}\right]^{-1}}

λ2=2incógnita¯[11+2(do21)]1{\displaystyle \lambda _{2}={\frac {2}{\bar {x}}}\left[1-{\sqrt {1+2(c^{2}-1)}}\right]^{-1}}

Estos estimadores se pueden derivar de los métodos de momentos estableciendo1λ1+1λ2=incógnita¯{\displaystyle {\frac {1}{\lambda _{1}}}+{\frac {1}{\lambda _{2}}}={\bar {x}}}yλ12+λ22λ1+λ2=do{\displaystyle {\frac {\sqrt {\lambda _{1}^{2}+\lambda _{2}^{2}}}{\lambda _{1}+\lambda _{2}}}=c}.

Los parámetros resultantesλ1{\displaystyle \lambda _{1}}yλ2{\displaystyle \lambda _{2}}son valores reales sido2[0,5,1]{\displaystyle c^{2}\in [0.5,1]}.

Caracterización

Una variable aleatoriaincógnitaHypago(λ1,,λk){\displaystyle X\sim Hypo(\lambda _{1},\dots ,\lambda _{k})}tiene una función de distribución acumulativa dada por,

F(incógnita)=1αmiincógnitaΘ1{\displaystyle F(x)=1-{\boldsymbol {\alpha }}e^{x\Theta }{\boldsymbol {1}}}

y función de densidad ,

F(incógnita)=αmiincógnitaΘΘ1,{\displaystyle f(x)=-{\boldsymbol {\alpha }}e^{x\Theta }\Theta {\boldsymbol {1}}\;,}

dónde1{\displaystyle {\boldsymbol {1}}}es un vector columna de unos de tamaño k ymiA{\displaystyle e^{A}}es la exponencial matricial de A. Cuandoλiλj{\displaystyle \lambda _{i}\neq \lambda _{j}}a pesar deij{\displaystyle i\neq j}, la función de densidad se puede escribir como

F(incógnita)=i=1kλimiincógnitaλi(j=1,jikλjλjλi)=i=1ki(0)λimiincógnitaλi{\displaystyle f(x)=\sum _{i=1}^{k}\lambda _{i}e^{-x\lambda _{i}}\left(\prod _{j=1,j\neq i}^{k}{\frac {\lambda _{j}}{\lambda _{j}-\lambda _{i}}}\right)=\sum _{i=1}^{k}\ell _{i}(0)\lambda _{i}e^{-x\lambda _{i}}}

dónde1(incógnita),,k(incógnita){\displaystyle \ell _{1}(x),\dots ,\ell _{k}(x)}son los polinomios de la base de Lagrange asociados con los puntosλ1,,λk{\displaystyle \lambda _{1},\dots ,\lambda _{k}}.

La distribución tiene transformada de Laplace de

L{F(incógnita)}=α(sIΘ)1Θ1{\displaystyle {\mathcal {L}}\{f(x)\}=-{\boldsymbol {\alpha }}(sI-\Theta )^{-1}\Theta {\boldsymbol {1}}}

Que se puede utilizar para encontrar momentos,

mi[incógnitanorte]=(1)nortenorte¡αΘnorte1.{\displaystyle E[X^{n}]=(-1)^{n}n!{\boldsymbol {\alpha }}\Theta ^{-n}{\boldsymbol {1}}\;.}

Caso general

En el caso general donde haya{\displaystyle a}sumas distintas de distribuciones exponenciales con tasasλ1,λ2,,λa{\displaystyle \lambda _{1},\lambda _{2},\cdots ,\lambda _{a}}y un número de términos en cada suma es igual ar1,r2,,ra{\displaystyle r_{1},r_{2},\cdots ,r_{a}}respectivamente. La función de distribución acumulativa parat0{\displaystyle t\geq 0}es dado por

F(t)=1(j=1aλjrj)k=1al=1rkΨk,l(λk)trklexp(λkt)(rkl)¡(l1)¡,{\displaystyle F(t)=1-\left(\prod _{j=1}^{a}\lambda _{j}^{r_{j}}\right)\sum _{k=1}^{a}\sum _{l=1}^{r_{k}}{\frac {\Psi _{k,l}(-\lambda _{k})t^{r_{k}-l}\exp(-\lambda _{k}t)}{(r_{k}-l)!(l-1)!}},}

con

Ψk,l(incógnita)=l1incógnital1(j=0,jka(λj+incógnita)rj).{\displaystyle \Psi _{k,l}(x)=-{\frac {\partial ^{l-1}}{\partial x^{l-1}}}\left(\prod _{j=0,j\neq k}^{a}\left(\lambda _{j}+x\right)^{-r_{j}}\right).}

con la convención adicionalλ0=0,r0=1{\displaystyle \lambda _{0}=0,r_{0}=1}. [ 4 ]

Usos

Esta distribución se ha utilizado en genética de poblaciones, [ 5 ] biología celular, [ 6 ] [ 7 ] y teoría de colas. [ 8 ] [ 9 ]

Véase también

Referencias

  1. "Distribución hipoexponencial" . Centro de documentación de lenguajes y sistemas de Wolfram . Wolfram. 2012. Consultado el 27 de febrero de 2024 .
  2. Legros, Benjamin; Jouini, Oualid (2015). "Un enfoque algebraico lineal para el cálculo de sumas de variables aleatorias de Erlang" . Modelado Matemático Aplicado . 39 (16): 4971– 4977. doi : 10.1016/j.apm.2015.04.013 .
  3. Bolch, Gunter; Greiner, Stefan; de Meer, Hermann; Trivedi, Kishor S. (2006). Redes de colas y cadenas de Markov: modelado y evaluación del rendimiento con aplicaciones en informática (2.ª ed.). Wiley. pp. 24–25 . doi : 10.1002/0471791571 . ISBN   978-0-471-79157-7.
  4. Amari, Suprasad V.; Misra, Ravindra B. (1997). "Expresiones de forma cerrada para la distribución de la suma de variables aleatorias exponenciales". IEEE Transactions on Reliability . 46 (4): 519– 522. doi : 10.1109/24.693785 .
  5. Strimmer, Korbinian; Pybus, Oliver G. (2001). "Explorando la historia demográfica de las secuencias de ADN usando el gráfico de horizonte generalizado" . Biología molecular y evolución . 18 (12): 2298– 2305. doi : 10.1093/oxfordjournals.molbev.a003776 . PMID 11719579 . 
  6. Yates, Christian A.; Ford, Matthew J.; Mort, Richard L. (2017). "Una representación multietapa de la proliferación celular como un proceso de Markov" . Bulletin of Mathematical Biology . 79 (12): 2905– 2928. arXiv : 1705.09718 . doi : 10.1007/ s11538-017-0356-4 . PMC 5709504. PMID 29030804 .  
  7. Gavagnin, Enrico; Ford, Matthew J.; Mort, Richard L.; Rogers, Tim; Yates, Christian A. (2019). "La velocidad de invasión de los modelos de migración celular con distribuciones realistas del tiempo del ciclo celular". Journal of Theoretical Biology . 481 : 91–99 . arXiv : 1806.03140 . doi : 10.1016/j.jtbi.2018.09.010 . PMID 30219568 . 
  8. Călinescu, Malenia (agosto de 2009). "Pronóstico y planificación de capacidad para los servicios de ambulancia" (PDF) . Facultad de Ciencias . Universidad Libre de Ámsterdam . Archivado del original (PDF) el 15 de febrero de 2010.
  9. Bekker, René; Koeleman, Paulien M. (2011). "Programación de admisiones y reducción de la variabilidad en la demanda de camas" . Health Care Management Science . 14 (3): 237– 249. doi : 10.1007/s10729-011-9163-x . PMC 3158339. PMID 21667090 .  

Lecturas adicionales

  • MF Neuts . (1981) Soluciones matriciales geométricas en modelos estocásticos: un enfoque algorítmico, Capítulo 2: Distribuciones de probabilidad de tipo fase; Dover Publications Inc.
  • G. Latouche, V. Ramaswami. (1999) Introducción a los métodos analíticos matriciales en el modelado estocástico, 1.ª edición. Capítulo 2: Distribuciones PH; ASA SIAM,
  • Colm A. O'Cinneide (1999). Distribución de tipo fase: problemas abiertos y algunas propiedades , Communication in Statistic - Stochastic Models, 15(4), 731 757.
  • L. Leemis y J. McQueston (2008). Relaciones de distribución univariadas , The American Statistician, 62(1), 45—53.
  • S. Ross. (2007) Introducción a los modelos de probabilidad, 9.ª edición, Nueva York: Academic Press