¿Cuál es el nombre de esta distribución discreta (ecuación de diferencia recursiva) que obtuve?

11

Encontré esta distribución en un juego de computadora y quería aprender más sobre su comportamiento. Viene de la decisión de si un determinado evento debe ocurrir después de un número determinado de acciones del jugador. Los detalles más allá de esto no son relevantes. Parece aplicable a otras situaciones, y lo encontré interesante porque es fácil de calcular y crea una larga cola.

Cada paso , el juego genera un número aleatorio uniforme 0 X < 1 . Si X < p ( n ) , entonces se desencadena el evento. Una vez que se produce el evento, el juego restablece n = 0 y vuelve a ejecutar la secuencia. Solo estoy interesado en una ocurrencia del evento para este problema, porque eso representa la distribución que está usando el juego. (Además, cualquier pregunta relacionada con múltiples ocurrencias se puede responder con un solo modelo de ocurrencia).n0X<1X<p(n)n=0

La "anormalidad" principal aquí es que el parámetro de probabilidad en esta distribución aumenta con el tiempo, o dicho de otra manera, el umbral aumenta con el tiempo. En el ejemplo, cambia linealmente, pero supongo que podrían aplicarse otras reglas. Después de pasos, o acciones del usuario,n

p(n)=kn

para alguna constante . En cierto punto n max , obtenemos p ( n max ) 1 . El evento simplemente se garantiza que ocurra en ese paso.0<k<1nmaxp(nmax)1

Pude determinar eso

y F ( n ) = p ( n ) + F ( n - 1 ) [ 1 - p ( n ) ] para PMF f ( n ) y CDF F ( n ) . En resumen, la probabilidad de que el evento ocurra en el n

f(n)=p(n)[1F(n1)]
F(n)=p(n)+F(n1)[1p(n)]
f(n)F(n)nEl paso es igual a la probabilidad , menos la probabilidad de que ya haya sucedido en cualquier paso anterior.p(n)

Aquí hay una trama de nuestro amigo Monte Carlo, por diversión, con . La mediana es de 21 y promedio de 22. k0.003ingrese la descripción de la imagen aquí

Esto es ampliamente equivalente a una ecuación de diferencia de primer orden del procesamiento de señal digital, que es mi experiencia, por lo que me pareció bastante novedoso. También me intriga la idea de que podría variar de acuerdo con cualquier fórmula arbitraria.p(n)

Mis preguntas:

  1. ¿Cuál es el nombre de esta distribución, si tiene una?
  2. f(n)F(n)
  3. ¿Hay otros ejemplos de distribuciones recursivas discretas como esta?

Ediciones Se aclaró el proceso sobre la generación de números aleatorios.

jbarlow
fuente
1
¿Alguna razón por la que elegiste corchetes en lugar de ()?
Cam.Davidson.Pilon
1
@ Cam.Davidson.Pilon: Mi fondo DSP se coló. Tendemos a usar corchetes para funciones de tiempo discreto. Supongo que esto debe ser discordante, así que lo cambiaré.
jbarlow
1
nXX<p(n)X
2
p(n)=kn0<k<1k1p(n)n>1/kp(n)nfunción de riesgo en el subcampo de estadística conocido comoAnálisis de supervivencia .
cardenal
1
kF f1/k=33318kp(n)1F1

Respuestas:

9

En cierto sentido, lo que ha hecho es caracterizar todas las distribuciones no negativas de valores enteros.

Dejemos de lado la descripción del proceso aleatorio por un momento y centrémonos en las recurrencias en la pregunta.

fn=pn(1Fn1)Fn=pn+(1pn)Fn1 Sn=1Fn=P(T>n)TF

Sn=1Fn=(1pn)Sn1,
Sn=k=0n(1pk).
(pn)[0,1]n

Más específicamente,

(pn)[0,1]

n=0log(1pn)=,

(pn)[0,1]

(pn)[0,1]0<pn<1nNpn=0n>N

¡Pero espera hay mas!

Ff

h(t)=f(t)S(t).

Λ(t)=0th(s)ds

S(t)=exp(Λ(t))=exp(0th(s)ds).
hh(t)0t0th(s)dst

t>t0

S(t)=et0th(s)dsS(t0).

h(t)S(t)

Conectando de nuevo al caso discreto

S(n)

h(t)=hn=log(1pn),
(n1,n](pn)

pnlog(1pn)pn=fn/Sn1

pn=knfnn=k1fn=0n>k1

cardenal
fuente
1
kk=1/2f=(0,1/2,1/2,0,)k=1/mf(m+1)=f(m+2)==0
2
fnFnfn
2
Gran respuesta. Esto es muy perspicaz. Estaba realmente interesado en ver este problema conectado con otras áreas y conceptos.
jbarlow
1
@jbarlow: gracias. ¡Me alegra que lo hayas encontrado útil! Disfruté pensar un poco al respecto, ya que es una buena pregunta.
cardenal
8

p(n)=p<1

F(n)=p+F(n1)(1p);F(0)=p

tiene la solución

F(n)=P(Nn)=1(1p)n+1

p(n)

Otros casos:

  1. p(n)=pn;p<1;F(0)=p
    F(n)=1(1p)Γ(n+1p)Γ(1p)Γ(n+1)
  2. S(n)=1F(n)
    S(n)=(1p(n))S(n1)
  3. p(n)np(n)=knp>1
    p(n)=1(1p)n+1p<1
    p(0)=pp(n)
    F(n)=1(1p)n+1n!
    S(n)=1F(n)=(1p)n+1n!
    i=0S(i)=E[N]
    E[N]=(1p)e(1p)
Cam.Davidson.Pilon
fuente
2
Cam, esa no es la función de peligro, sino la función de supervivencia . :-)
cardenal
1
ty, * edited to survival
Cam.Davidson.Pilon