Probabilidad de que un proceso de Poisson independiente supere a otro

9

He hecho esta pregunta antes de otra manera en otros intercambios de pila, lo siento mucho por la nueva publicación.

Le pregunté a mi profesor y a un par de estudiantes de doctorado, sin una respuesta definitiva. Primero expondré el problema, luego mi solución potencial y el problema con mi solución, lo siento mucho por el muro de texto.

El problema:

Suponga que dos procesos independientes de Poisson y , con y para el mismo intervalo, sujeto a . ¿Cuál es la probabilidad de que en cualquier momento, como el tiempo tiende al infinito, que la salida agregada del proceso sea ​​mayor que la salida agregada del proceso más , es decir, ? Para ilustrar con un ejemplo, suponga que dos puentes y , en promedio, los y sobre el puente yR λ R λ M λ R > λ M M R D P ( M > R + D ) R M λ R λ M R MMRλRλMλR>λMMRDP(M>R+D)RMλRλMRMrespectivamente por intervalo, y . autos ya han conducido sobre el puente , ¿cuál es la probabilidad de que en algún momento más autos en total hayan manejado sobre el puente que ? D RλR>λMDRMR

Mi forma de resolver este problema:

Primero definimos dos procesos de Poisson:

M(I)Poisson(μMI)R(I)Poisson(μRI)

El siguiente paso es encontrar una función que describe después de un número determinado de intervalos . Esto sucederá en el caso de que esté condicionado a la salida de , para todos los valores no negativos de . Para ilustrar, si la producción total de es entonces la producción total de tiene que ser más grande que . Como se muestra abajo.I M ( I ) > k + D R ( I ) = k k R X M X + DP(M>R+D)IM(I)>k+DR(I)=kkRXMX+D

P(M(I))>R(I)+D)=k=0n[P(M(I)>k+DR(I)=k)]

n

Debido a la independencia, esto puede reescribirse como el producto de los dos elementos, donde el primer elemento es 1-CDF de la distribución de Poisson y el segundo elemento es el Poisson pmf:

P(M(I)>R(I)+D)=k=0n[P(M(I)>k+D)1Poisson CDFP(R(I)=k)Poisson pmf]

n

Para crear un ejemplo, suponga que , y , a continuación se muestra el gráfico de esa función sobre :λ R = 0.6 λ M = 0.4 ID=6λR=0.6λM=0.4I

ingrese la descripción de la imagen aquí

El siguiente paso es encontrar la probabilidad de que esto ocurra en cualquier punto en el tiempo, deja de llamadas que . Mi idea es que esto es equivalente a encontrar 1 menos la probabilidad de que no estar nunca por encima de . Es decir, dejar enfoque infinito lo que es condicionada a esto también ser cierto para todos los valores anteriores de .M R + D N P ( R ( N ) + D M ( N ) ) NQMR+DNP(R(N)+DM(N))N

1 - P ( M ( I ) > R ( I ) + D )P(R(I)+DM(I)) es lo mismo que , definamos eso como la función g (I):1P(M(I)>R(I)+D)

g(I)=1P(M(I)>R(I)+D)

Como tiende al infinito, esto también puede reescribirse como la integral geométrica sobre la función .Ng(I)

Q=1exp(0Nln(g(I))dI)

Q=1exp(0Nln(1P(M(I)>R(I)+D))dI)

N

Donde tenemos la función de desde arriba.P(M(I)>R(I)+D)

Q=1exp(0Nln(1k=0n[P(M(I)>k+D)1Poisson CDFP(R(I)=k)Poisson pmf])dI)

N

n

Ahora para mí esto debería darme el valor final de , para cualquier , y . Sin embargo, hay un problema, deberíamos poder reescribir las lambdas como queramos, ya que lo único que debería importar es su proporción entre sí. Para construir sobre el ejemplo anterior con , y , esto es efectivamente lo mismo que , y , siempre que su intervalo se divida entre 10. Es decir, 10 autos cada 10 minutos es lo mismo que 1 auto cada minuto. Sin embargo, hacer esto produce un resultado diferente. ,QDλRλMD=6λR=0.6λM=0.4D=6λR=0.06λM=0.04D=6λR=0.6 y produce un de y , y produce un de . La realización inmediata es que , y la razón es bastante simple si comparamos los gráficos de los dos resultados, el siguiente gráfico muestra la función para , y .λM=0.4Q0.5856116D=6λR=0.06λM=0.04Q0.99985071(10.5856116)10=0.9998507D=6λR=0.06λM=0.04

ingrese la descripción de la imagen aquí

Como se puede ver, la probabilidad no cambia, sin embargo, ahora se necesitan diez veces más intervalos para llegar a la misma probabilidad. Como depende del intervalo de la función, esto naturalmente tiene una implicación. Obviamente, esto significa que algo está mal, ya que el resultado no debería depender de mi lambda inicial, especialmente porque no hay una lambda inicial que sea correcta y es tan correcta como y o y , etc., siempre que el intervalo se escala en consecuencia. Por lo tanto, si bien puedo escalar fácilmente la probabilidad, es decir, pasar de y a yQ0.040.060.40.611.50.40.60.040.06 es lo mismo que escalar la probabilidad con un factor de 10. Esto obviamente produce el mismo resultado, pero como todas estas lambdas son puntos de partida igualmente válidos, entonces esto obviamente no es correcto.

Para mostrar este impacto, graficé en función de , donde es un factor de escala de las lambdas, con lambdas iniciales de y . La salida se puede ver en el siguiente gráfico:QttλM=0.4λR=λM1.5

ingrese la descripción de la imagen aquí

Aquí es donde estoy atrapado, para mí el enfoque se ve bien y correcto, pero el resultado es obviamente incorrecto. Mi pensamiento inicial es que me falta una nueva escala fundamental en alguna parte, pero por mi vida no puedo averiguar dónde.

Gracias por leer, toda ayuda es muy apreciada.

Además, si alguien quiere mi código R, avíseme y lo subiré.

no nein
fuente
Hice una limpieza bastante extensa de su código MathJax. Si echas un vistazo, verás algunas cosas sobre el uso estándar y adecuado. (Se podría hacer más trabajo; tal vez más tarde.)
Michael Hardy
¡Increíble! Muchas gracias, no lo sabía, ¿hay alguna guía específica que deba seguir?
no nein
Edité algunas cosas adicionales en línea con lo que hiciste.
no nein
@nonein Hay un poquito en la ayuda de edición, pero más allá de eso hay matemáticas. Tutorial básico MathJax de SE y referencia rápida . Las guías sobre cómo escribir matemáticas en LaTeX (que son fáciles de google) a menudo ayudan si está tratando de encontrar algo que no esté cubierto en la referencia rápida allí (aunque ahora tiene una cobertura bastante completa del subconjunto de MathJax).
Glen_b

Respuestas:

3

Deje que los tiempos colectivos de los procesos sean Debido a que estos son procesos independientes de Poisson, casi con seguridad se observa exactamente uno de ellos en cada uno de estos momentos. Para definaT=(0=t0<t1<t2<).i>0,

B(i)={+1if R(ti)=11if M(ti)=1

y acumule en el proceso es decir, y para todo cuenta cuántas veces más ha aparecido que justo después del tiempoB(i)W:W(0)=0W(i+1)=W(i)+B(i)i>0. W(i)RMti.

Figura: simulación

Esta figura muestra realizaciones de (en rojo) y (en azul medio) como "parcelas de alfombra" en la parte superior. Los puntos trazan los valores de . Cada punto rojo representa un aumento en el exceso de mientras que cada punto azul muestra una disminución en el exceso.RM(ti,W(i))R(ti)M(ti)

Para deje que sea ​​casualidad de que al menos uno de los sea ​​menor o igual que y sea su probabilidad.b=0,1,2,,EbWibf(b)

La pregunta pidef(D+1).

Deje Esta es la tasa de los procesos combinados. es una caminata aleatoria binomial, porqueλ=λR+λM.W

Pr(B(i)=1)=λRλ and Pr(B(i)=1)=λMλ.

Así,

La respuesta es igual a la posibilidad de que esta caminata aleatoria binomial encuentre una barrera absorbente enWD1.

La forma más elemental de encontrar esta oportunidad observa que

f(0)=1

porquey, para todo los dos siguientes pasos posibles de producen de forma recursivaW(0)=0;b>0,±1

f(b)=λRλf(b+1)+λMλf(b1).

Suponiendo la solución única para esλRλM,b0

f(b)=(λMλR)b,

como puede verificar al conectar esto en las ecuaciones definitorias anteriores. Así,

La respuesta es

Pr(ED+1)=f(D+1)=(λMλR)D+1.

whuber
fuente