Articulo de referencia

Teoría de la renovación

La teoría de la renovación es la rama de la teoría de la probabilidad que generaliza el proceso de Poisson para tiempos de retención arbitrarios. En lugar de tiempos de retenció...

La teoría de la renovación es la rama de la teoría de la probabilidad que generaliza el proceso de Poisson para tiempos de retención arbitrarios. En lugar de tiempos de retención con distribución exponencial , un proceso de renovación puede tener tiempos de retención independientes e idénticamente distribuidos (IID) con esperanza finita. Un proceso de renovación-recompensa también tiene una secuencia aleatoria de recompensas incurridas en cada tiempo de retención, las cuales son IID pero no necesariamente independientes de los tiempos de retención.

Un proceso de renovación tiene propiedades asintóticas análogas a la ley fuerte de los grandes números y al teorema del límite central . La función de renovaciónmetro(t){\displaystyle m(t)}(número esperado de llegadas) y función de recompensagramo(t){\displaystyle g(t)}(valor de recompensa esperado) son de vital importancia en la teoría de renovación. La función de renovación satisface una ecuación integral recursiva, la ecuación de renovación. La ecuación de renovación clave proporciona el valor límite de la convolución demetro(t){\displaystyle m'(t)}con una función no negativa adecuada. La superposición de procesos de renovación puede estudiarse como un caso especial de procesos de renovación de Markov .

Las aplicaciones incluyen calcular la mejor estrategia para reemplazar maquinaria desgastada en una fábrica; comparar los beneficios a largo plazo de diferentes pólizas de seguro; y modelar la transmisión de enfermedades infecciosas, donde "Uno de los medios más ampliamente adoptados para inferir el número de reproducción es a través de la ecuación de renovación". [ 1 ] La paradoja de la inspección se relaciona con el hecho de que observar un intervalo de renovación en el tiempo t da un intervalo con un valor promedio mayor que el de un intervalo de renovación promedio.

Procesos de renovación

Introducción

El proceso de renovación es una generalización del proceso de Poisson . En esencia, el proceso de Poisson es un proceso de Markov de tiempo continuo sobre los enteros positivos (generalmente comenzando en cero) que tiene tiempos de retención independientes distribuidos exponencialmente en cada entero.i{\displaystyle i}antes de avanzar al siguiente entero,i+1{\displaystyle i+1}En un proceso de renovación, los tiempos de retención no necesitan tener una distribución exponencial; más bien, los tiempos de retención pueden tener cualquier distribución en los números positivos, siempre que los tiempos de retención sean independientes e idénticamente distribuidos ( IID ) y tengan una media finita.

Definición formal

Ejemplo de evolución de un proceso de renovación con tiempos de retención S i y tiempos de salto J n .

Dejar(Si)i1{\displaystyle (S_{i})_{i\geq 1}}sea ​​una secuencia de variables aleatorias positivas, independientes e idénticamente distribuidas con un valor esperado finito.

0<mi[Si]<.{\displaystyle 0<\operatorname {E} [S_{i}]<\infty .}

Nos referimos a la variable aleatoriaSi{\displaystyle S_{i}}como el "i{\displaystyle i}-tiempo de espera".

Definir para cada n > 0  :

Jnorte=i=1norteSi,{\displaystyle J_{n}=\sum _{i=1}^{n}S_{i},}

cadaJnorte{\displaystyle J_{n}}se le conoce como el "norte{\displaystyle n}-tiempo de salto" y los intervalos[Jnorte,Jnorte+1]{\displaystyle [J_{n},J_{n+1}]}se denominan "intervalos de renovación".

Entonces(incógnitat)t0{\displaystyle (X_{t})_{t\geq 0}}viene dada por una variable aleatoria

incógnitat=norte=1I{Jnortet}=sorber{norte:Jnortet}{\displaystyle X_{t}=\sum _{n=1}^{\infty }\operatorname {\mathbb {I} } _{\{J_{n}\leq t\}}=\sup \left\{\,n:J_{n}\leq t\,\right\}}

dóndeI{Jnortet}{\displaystyle \operatorname {\mathbb {I} } _{\{J_{n}\leq t\}}}es la función indicadora

I{Jnortet}={1,si Jnortet0,de lo contrario{\displaystyle \operatorname {\mathbb {I} } _{\{J_{n}\leq t\}}={\begin{cases}1,&{\text{si }}J_{n}\leq t\\0,&{\text{en otro caso}}\end{cases}}}

(incógnitat)t0{\displaystyle (X_{t})_{t\geq 0}}representa el número de saltos que han ocurrido hasta el tiempo t , y se denomina proceso de renovación.

Interpretación

Si se consideran eventos que ocurren en momentos aleatorios, se puede optar por pensar en los tiempos de retención.{Si:i1}{\displaystyle \{S_{i}:i\geq 1\}}como el tiempo aleatorio transcurrido entre dos eventos consecutivos. Por ejemplo, si el proceso de renovación modela la cantidad de averías de diferentes máquinas, entonces el tiempo de espera representa el tiempo entre la avería de una máquina y la de otra.

El proceso de Poisson es el único proceso de renovación con la propiedad de Markov , [ 2 ] así como la distribución exponencial es la única variable aleatoria continua con la propiedad de falta de memoria.

Procesos de renovación y recompensa

Ejemplo de evolución de un proceso de renovación-recompensa con tiempos de retención S i , tiempos de salto J n y recompensas W i

DejarW1,W2,{\displaystyle W_{1},W_{2},\ldots }ser una secuencia de variables aleatorias IID ( recompensas ) que satisfacen

mi|Wi|<.{\displaystyle \operatorname {E} |W_{i}|<\infty .\,}

Entonces la variable aleatoria

Yt=i=1incógnitatWi{\displaystyle Y_{t}=\sum _{i=1}^{X_{t}}W_{i}}

se denomina proceso de renovación-recompensa . Tenga en cuenta que, a diferencia delSi{\displaystyle S_{i}}, cadaWi{\displaystyle W_{i}}puede tomar valores negativos además de valores positivos.

La variable aleatoriaYt{\displaystyle Y_{t}}depende de dos secuencias: los tiempos de retenciónS1,S2,{\displaystyle S_{1},S_{2},\ldots }y las recompensas W1,W2,{\displaystyle W_{1},W_{2},\ldots }Estas dos secuencias no tienen por qué ser independientes. En particular,Wi{\displaystyle W_{i}}puede ser una función deSi{\displaystyle S_{i}}.

Interpretación

En el contexto de la interpretación anterior de los tiempos de espera como el tiempo entre fallas sucesivas de una máquina, las "recompensas"W1,W2,{\displaystyle W_{1},W_{2},\ldots }(que en este caso resultan ser negativos) pueden considerarse como los sucesivos costos de reparación incurridos como resultado de las sucesivas fallas.

Una analogía alternativa es que tenemos una gallina mágica que pone huevos a intervalos (tiempos de espera) distribuidos comoSi{\displaystyle S_{i}}A veces pone huevos de oro de peso aleatorio, y a veces pone huevos tóxicos (también de peso aleatorio) que requieren una eliminación responsable (y costosa). Las "recompensas"Wi{\displaystyle W_{i}}son las sucesivas pérdidas/ganancias financieras (aleatorias) resultantes de huevos sucesivos ( i = 1,2,3,...) yYt{\displaystyle Y_{t}}registra la "recompensa" financiera total en el momento t .

Función de renovación

Definimos la función de renovación como el valor esperado del número de saltos observados hasta cierto tiempo.t{\displaystyle t}:

metro(t)=mi[incógnitat].{\displaystyle m(t)=\operatorname {E} [X_{t}].\,}

Teorema elemental de renovación

La función de renovación satisface

límitet1tmetro(t)=1mi[S1].{\displaystyle \lim _{t\to \infty }{\frac {1}{t}}m(t)={\frac {1}{\operatorname {E} [S_{1}]}}.}

Teorema elemental de renovación para procesos de recompensa por renovación

Definimos la función de recompensa :

gramo(t)=mi[Yt].{\displaystyle g(t)=\operatorname {E} [Y_{t}].\,}

La función de recompensa satisface

límitet1tgramo(t)=mi[W1]mi[S1].{\displaystyle \lim _{t\to \infty }{\frac {1}{t}}g(t)={\frac {\operatorname {E} [W_{1}]}{\operatorname {E} [S_{1}]}}.}

Ecuación de renovación

La función de renovación satisface

metro(t)=FS(t)+0tmetro(ts)FS(s)ds{\displaystyle m(t)=F_{S}(t)+\int _{0}^{t}m(ts)f_{S}(s)\,ds}

dóndeFS{\displaystyle F_{S}}es la función de distribución acumulativa deS1{\displaystyle S_{1}}yFS{\displaystyle f_{S}}es la función de densidad de probabilidad correspondiente.

Teorema clave de renovación

Sea X un proceso de renovación con función de renovaciónmetro(t){\displaystyle m(t)}y la renovación intermedia significaμ{\displaystyle \mu }. Dejargramo:[0,)[0,){\displaystyle g:[0,\infty )\rightarrow [0,\infty )}ser una función que satisfaga:

  • 0gramo(t)dt<{\displaystyle \int _{0}^{\infty }g(t)\,dt<\infty }
  • g es monótona y no creciente.

El teorema clave de renovación establece que, comot{\displaystyle t\rightarrow \infty }: [ 4 ]

0tgramo(tincógnita)metro(incógnita)dincógnita1μ0gramo(incógnita)dincógnita{\displaystyle \int _{0}^{t}g(t-x)m'(x)\,dx\rightarrow {\frac {1}{\mu }}\int _{0}^{\infty }g(x)\,dx}

Teorema de renovación

En vista degramo(incógnita)=I[0,h](incógnita){\displaystyle g(x)=\mathbb {I} _{[0,h]}(x)}para cualquierh>0{\displaystyle h>0}da como caso especial el teorema de renovación: [ 5 ]

metro(t+h)metro(t)hμ{\displaystyle m(t+h)-m(t)\rightarrow {\frac {h}{\mu }}}comot{\displaystyle t\rightarrow \infty }

El resultado puede demostrarse mediante ecuaciones integrales o mediante un argumento de acoplamiento . [ 6 ] Aunque es un caso especial del teorema de renovación clave, puede utilizarse para deducir el teorema completo, considerando funciones escalonadas y luego secuencias crecientes de funciones escalonadas. [ 4 ]

Propiedades asintóticas

Los procesos de renovación y los procesos de renovación-recompensa tienen propiedades análogas a la ley fuerte de los grandes números , que se puede derivar del mismo teorema. Si(incógnitat)t0{\displaystyle (X_{t})_{t\geq 0}}es un proceso de renovación y(Yt)t0{\displaystyle (Y_{t})_{t\geq 0}}Se trata entonces de un proceso de renovación y recompensa:

límitet1tincógnitat=1mi[S1]{\displaystyle \lim _{t\to \infty }{\frac {1}{t}}X_{t}={\frac {1}{\operatorname {E} [S_{1}]}}}[ 7 ]
límitet1tYt=1mi[S1]mi[W1]{\displaystyle \lim _{t\to \infty }{\frac {1}{t}}Y_{t}={\frac {1}{\operatorname {E} [S_{1}]}}\operatorname {E} [W_{1}]}

casi con seguridad.

Los procesos de renovación tienen además una propiedad análoga al teorema del límite central : [ 7 ]

incógnitatt/μtσ2/μ3norte(0,1){\displaystyle {\frac {X_{t}-t/\mu }{\sqrt {t\sigma ^{2}/\mu ^{3}}}}\to {\mathcal {N}}(0,1)}

paradoja de la inspección

El intervalo de renovación determinado por el punto aleatorio t (mostrado en rojo) es estocásticamente mayor que el primer intervalo de renovación.

Una característica curiosa de los procesos de renovación es que si esperamos un tiempo predeterminado t y luego observamos cuán grande es el intervalo de renovación que contiene t , deberíamos esperar que sea típicamente mayor que un intervalo de renovación de tamaño promedio.

Matemáticamente, la paradoja de la inspección establece que, para cualquier t > 0, el intervalo de renovación que contiene a t es estocásticamente mayor que el primer intervalo de renovación. Es decir, para todo x > 0 y para todo t > 0:

PAG(Sincógnitat+1>incógnita)PAG(S1>incógnita)=1FS(incógnita){\displaystyle \operatorname {P} (S_{X_{t}+1}>x)\geq \operatorname {P} (S_{1}>x)=1-F_{S}(x)}

donde F S es la función de distribución acumulativa de los tiempos de espera IID S i . Un ejemplo ilustrativo es la paradoja del tiempo de espera del autobús : para una distribución aleatoria dada de llegadas de autobuses, el pasajero promedio en una parada de autobús observa más retrasos que el operador promedio de los autobuses.

La solución a la paradoja reside en que nuestra distribución muestreada en el instante t está sesgada por el tamaño (véase sesgo de muestreo ), ya que la probabilidad de que se elija un intervalo es proporcional a su tamaño. Sin embargo, un intervalo de renovación de tamaño medio no está sesgado por el tamaño.

Superposición

A menos que el proceso de renovación sea un proceso de Poisson, la superposición (suma) de dos procesos de renovación independientes no es un proceso de renovación. [ 8 ] Sin embargo, tales procesos pueden describirse dentro de una clase más amplia de procesos llamados procesos de renovación de Markov . [ 9 ] Sin embargo, la función de distribución acumulativa del primer tiempo entre eventos en el proceso de superposición viene dada por [ 10 ]

R(t)=1k=1Kαkl=1Kαl(1Rk(t))j=1,jkKαjt(1Rj())d{\displaystyle R(t)=1-\sum _{k=1}^{K}{\frac {\alpha _{k}}{\sum _{l=1}^{K}\alpha _{l}}}(1-R_{k}(t))\prod _{j=1,j\neq k}^{K}\alpha _{j}\int _{t}^{\infty }(1-R_{j}(u))\,{\text{d}}u}

donde R k ( t ) y α k  >  0 son la función de distribución acumulada de los tiempos entre eventos y la tasa de llegada del proceso k . [ 11 ]

Ejemplo de aplicación

Eric, el empresario, posee n máquinas, cada una con una vida útil uniformemente distribuida entre cero y dos años. Eric puede dejar que cada máquina funcione hasta que falle, con un coste de reposición de 2600 €; o bien, puede reemplazarla en cualquier momento mientras aún funcione, con un coste de 200 €.

¿Cuál es su política de reemplazo óptima?

Véase también

Notas

  1. Green, William D.; Ferguson, Neil M.; Cori, Anne (2022). "Inferencia del número de reproducción mediante la ecuación de renovación en epidemias heterogéneas" . Journal of the Royal Society Interface . 19 (188) 20210429. doi : 10.1098/rsif.2021.0429 . PMC 8965414. PMID 35350879 .  
  2. Grimmett y Stirzaker (1992) , pág. 393.
  3. Grimmett y Stirzaker (1992) , pág. 390.
  4. 1 2 3 Grimmett y Stirzaker (1992) , pág. 395.
  5. Feller (1971) , págs. 347–351.
  6. Grimmett y Stirzaker (1992) , págs. 394-395.
  7. 1 2 Grimmett y Stirzaker (1992) , pág. 394.
  8. Grimmett y Stirzaker (1992) , pág. 405.
  9. Çinlar, Erhan (1969). "Teoría de renovación de Markov". Avances en probabilidad aplicada . 1 (2). Applied Probability Trust: 123– 187. doi : 10.2307/1426216 . JSTOR 1426216 . 
  10. Lawrence, AJ (1973). "Dependencia de los intervalos entre eventos en procesos de superposición". Journal of the Royal Statistical Society. Serie B (Metodológica) . 35 (2): 306–315 . doi : 10.1111/j.2517-6161.1973.tb00960.x . JSTOR 2984914 . fórmula 4.1
  11. Choungmo Fofack, Nicaise; Nain, Philippe; Neglia, Giovanni; Towsley, Don (6 de marzo de 2012). Análisis de redes de caché basadas en TTL . Actas de la 6.ª Conferencia Internacional sobre Metodologías y Herramientas de Evaluación del Rendimiento (informe) . Recuperado el 15 de noviembre de 2012 .

Referencias

  • Cox, David (1970). Teoría de la renovación . Londres: Methuen & Co. pág.  142. ISBN 0-412-20570-X.
  • Doob, JL (1948). "Teoría de la renovación desde el punto de vista de la teoría de la probabilidad" (PDF) . Transactions of the American Mathematical Society . 63 (3): 422– 438. doi : 10.2307/1990567 . JSTOR 1990567 . 
  • Feller, William (1971). Introducción a la teoría de la probabilidad y sus aplicaciones . Vol.  2 (segunda  edición). Wiley.
  • Grimmett, GR ; Stirzaker, DR (1992). Probabilidad y procesos aleatorios (segunda  ed.). Oxford University Press. ISBN 0-19-857222-0.
  • Smith, Walter L. (1958). "Teoría de la renovación y sus ramificaciones". Journal of the Royal Statistical Society, Serie B. 20 ( 2): 243– 302. doi : 10.1111/j.2517-6161.1958.tb00294.x . JSTOR 2983891 . 
  • Wanli Wang, Johannes HP Schulz, Weihua Deng y Eli Barkai (2018). "Teoría de renovación con tiempos de permanencia distribuidos de cola pesada: típico versus raro". Phys. Rev. E . 98 (4) 042139. arXiv : 1809.05856 . Bibcode : 2018PhRvE..98d2139W . doi : 10.1103/PhysRevE.98.042139 . S2CID 54727926 . {{cite journal}}: CS1 maint: varios nombres: lista de autores ( enlace )