Probabilidade de um processo independente de Poisson ultrapassar outro

9

Já fiz essa pergunta de outra maneira em outras trocas de pilha, desculpe-me pela repostagem.

Eu perguntei ao meu professor e a alguns estudantes de doutorado sobre, sem uma resposta definitiva. Primeiro vou declarar o problema, depois minha solução em potencial e o problema com a minha solução, desculpe a parede de texto.

O problema:

Suponha dois Poisson independente processa M e R , com λR e λM para o mesmo intervalo de tempo, sujeito a λR>λM . Qual é a probabilidade de que, em qualquer ponto do tempo, à medida que o tempo tende ao infinito, a saída agregada do processo M seja maior que a saída agregada do processo R mais D , ou seja, P(M>R+D) . Para ilustrar com um exemplo, assuma duas pontes R e M , em média λR eλM carros passar por cima da ponteR eM , respectivamente, por cada intervalo, eλR>λM . D carros já têm impulsionado sobre a ponteR , qual é a probabilidade de que a qualquer ponto no tempo mais carros no total têm impulsionado sobre a ponteM deR .

Minha maneira de resolver esse problema:

Primeiro, definimos dois processos de Poisson:

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

O próximo passo é o de encontrar uma função que descreve após um determinado número de intervalos I . Isso acontecerá no caso de M ( I ) > k + D condicional à saída de R ( I ) = k , para todos os valores não negativos de k . Para ilustrar, se a saída do agregado de R é X , então a saída total de M tem de ser maior que X + D . Como mostrado abaixo.P(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

Devido à independência, isso pode ser reescrito como o produto dos dois elementos, onde o primeiro elemento é 1-CDF da distribuição Poisson e o segundo elemento é o Poisson pmf:

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

n

Para criar um exemplo, assuma , λ R = 0,6 e λ M = 0,4 , abaixo está o gráfico dessa função sobre I :D=6λR=0.6λM=0.4I

insira a descrição da imagem aqui

O próximo passo é encontrar a probabilidade de que isso aconteça em qualquer ponto no tempo, permite chamada que . Meu pensamento é que este é equivalente a encontrar um menos a probabilidade de M nunca estar acima de R + D . Isto é deixar N abordagem infinito que é P ( R ( N ) + D H ( N ) ) condicionada por este também ser verdadeiro para todos os valores anteriores de N .QMR+DNP(R(N)+DM(N))N

é o mesmo que 1 - P ( M ( I ) > R ( I ) + D ) , vamos definir isso como função g (I):P(R(I)+DM(I))1P(M(I)>R(I)+D)

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

Como tende ao infinito, isso também pode ser reescrito como a integral geométrica sobre a função g ( I ) .Ng(I)

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

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

N

Onde temos a função de de cima.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

Agora, para mim isso deve me dar o valor final da , para um dado D , λ R e λ M . No entanto, há um problema, devemos ser capazes de reescrever as lambdas como queremos, pois a única coisa que importa é a proporção entre si. Para desenvolver o exemplo anterior com D = 6 , λ R = 0,6 e λ M = 0,4 , é efetivamente o mesmo que D = 6 , λ R = 0,06 e λ M = 0,04QDλRλMD=6λR=0.6λM=0.4D=6λR=0.06λM=0.04, contanto que o intervalo seja dividido por 10. Ou seja, 10 carros a cada 10 minutos é o mesmo que 1 carro a cada minuto. No entanto, fazer isso produz um resultado diferente. , λ R = 0,6 e λ H = 0,4 rendimentos um Q de 0,5856116 e D = 6 , λ R = 0,06 e λ M = 0,04 rendimentos um Q de 0,9998507 . A realização imediata é que 1 - ( 1 -D=6λR=0.6λM=0.4Q0.5856116D=6λR=0.06λM=0.04Q0.9998507 , e o motivo é realmente bastante simples se compararmos os gráficos dos dois resultados, o gráfico abaixo mostra a função para D = 6 , λ R = 0,06 e λ M = 0,04 .1(10.5856116)10=0.9998507D=6λR=0.06λM=0.04

insira a descrição da imagem aqui

Q0.040.060.40.611.50.40.60.040.06

QttλM=0.4λR=λM1.5

insira a descrição da imagem aqui

É aqui que estou preso, para mim a abordagem parece correta e correta, mas o resultado está obviamente errado. Meu pensamento inicial é que estou perdendo uma redimensionamento fundamental em algum lugar, mas não consigo descobrir onde.

Obrigado pela leitura, toda e qualquer ajuda é muito apreciada.

Além disso, se alguém quiser meu código R, informe-me e eu o carregarei.

não nein
fonte
Fiz uma limpeza bastante extensa do seu código MathJax. Se você der uma olhada, verá algumas coisas sobre o uso padrão e adequado. (Mais trabalho poderia ser feito, talvez mais tarde.)
Michael Hardy
Impressionante! Muito obrigado, eu não sabia disso, existe um guia específico que devo seguir?
No nein
Eu editei algumas coisas extras de acordo com o que você fez.
No nein
@nonein Há um pouquinho na ajuda de edição, mas além disso há o tutorial básico e a referência rápida do math.SE's MathJax . Os guias sobre como escrever matemática no LaTeX (que são fáceis de pesquisar no Google) geralmente ajudam se você estiver tentando encontrar algo não abordado na referência rápida (embora agora tenha uma cobertura bastante abrangente do subconjunto do MathJax).
Glen_b -Reinstala Monica

Respostas:

3

T=(0=t0<t1<t2<).i>0,

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

B(i)W:W(0)=0W(i+1)=W(i)+B(i)i>0. W(i)RMti.

Figura: simulação

RM(ti,W(i))R(ti)M(ti)

b=0,1,2,,EbWibf(b)

f(D+1).

λ=λR+λM.W

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

Portanto,

A resposta é igual à chance de que esse passeio aleatório binomial encontre uma barreira absorvente emWD1.

A maneira mais elementar de encontrar essa chance é observar que

f(0)=1

porquee, para todo os dois próximos passos possíveis de produzem recursivamenteW(0)=0;b>0,±1

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

Supondo a solução exclusiva para éλRλM,b0

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

como você pode verificar, inserindo isso nas equações definidoras anteriores. Portanto,

A resposta é

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

whuber
fonte