Quel est le nom de cette distribution discrète (équation de différence récursive) que j'ai dérivée?

11

Je suis tombé sur cette distribution dans un jeu vidéo et je voulais en savoir plus sur son comportement. Cela vient de la décision de savoir si un certain événement doit se produire après un certain nombre d'actions de joueurs. Les détails au-delà de cela ne sont pas pertinents. Il semble applicable à d'autres situations, et je l'ai trouvé intéressant car il est facile à calculer et crée une longue queue.

A chaque étape , le jeu génère un nombre aléatoire uniforme . Si , alors l'événement est déclenché. Une fois que l'événement s'est produit, le jeu réinitialise et exécute à nouveau la séquence. Je ne suis intéressé que par une occurrence de l'événement pour ce problème, car cela représente la distribution que le jeu utilise. (En outre, toutes les questions concernant plusieurs occurrences peuvent être répondues avec un seul modèle d'occurrence.)n0X<1X<p(n)n=0

La principale "anomalie" ici est que le paramètre de probabilité dans cette distribution augmente avec le temps, ou en d'autres termes, le seuil augmente avec le temps. Dans l'exemple, il change linéairement mais je suppose que d'autres règles pourraient s'appliquer. Après étapes ou actions de l'utilisateur,n

p(n)=kn

pour une constante . À un certain point , nous obtenons p (n _ {\ max}) \ geq 1 . L'événement est simplement garanti de se produire à cette étape.0<k<1nmaxp(nmax)1

J'ai pu déterminer que

F ( n ) = p ( n ) + F ( n - 1 ) [ 1 - p ( n ) ] f ( n ) F ( n ) n p ( n )

f(n)=p(n)[1F(n1)]
et pour PMF et CDF . En bref, la probabilité que l'événement se produise à la ème étape est égale à la probabilité , moins la probabilité qu'il se soit déjà produit à une étape précédente.
F(n)=p(n)+F(n1)[1p(n)]
f(n)F(n)np(n)

Voici une intrigue de notre ami Monte Carlo, pour le plaisir, avec . La médiane s'établit à 21 et la moyenne à 22. k0.003entrez la description de l'image ici

Cela équivaut à peu près à une équation de différence de premier ordre du traitement du signal numérique, qui est mon expérience, et j'ai donc trouvé cela assez nouveau. Je suis également intrigué par l'idée que pourrait varier selon n'importe quelle formule arbitraire.p(n)

Mes questions:

  1. Quel est le nom de cette distribution, si elle en a une?
  2. Existe-t-il un moyen de dériver une expression pour sans référence à ?F ( n )f(n)F(n)
  3. Existe-t-il d'autres exemples de distributions récursives discrètes comme celle-ci?

Modifie processus Clarification sur la génération de nombres aléatoires.

jbarlow
la source
1
Une raison pour laquelle vous avez choisi des crochets au lieu de ()?
Cam.Davidson.Pilon
1
@ Cam.Davidson.Pilon: Mon arrière-plan DSP s'est glissé dedans. Nous avons tendance à utiliser des crochets pour les fonctions temporelles discrètes. Je suppose que cela doit être choquant, donc je vais le changer.
jbarlow
1
Le processus que vous supposez n'apparaît pas clairement défini ici. Vous dites "A chaque étape , le jeu lance un nombre aléatoire Si , alors l'événement est déclenché." Mais, vous ne donnez aucune spécification sur la façon dont est dessiné. Je pense qu'il serait utile que le processus puisse être décrit un peu plus précisément. X X < p ( n ) XnXX<p(n)X
cardinal
2
@jbarlow: Désolé si ma remarque précédente n'était pas claire. Si pour quelque , il n'y a aucun moyen que votre processus puisse avoir plus de étapes, car un nombre aléatoire uniforme entre zéro et un serait certainement plus petit que pour tout . La quantité en fonction de est très étroitement liée à ce que l'on appelle la fonction de risque dans le sous-domaine des statistiques appelé analyse de survie . 0 < k < 1 k - 1p ( n ) n > 1 / k p ( n ) np(n)=kn0<k<1k1p(n)n>1/kp(n)n
Cardinal
1
Pour les petits , l'utilisation de l'analogue différentiel de cette équation de différence montre que ( pas !) Est proche de la gaussienne. (De cela, nous déduisons immédiatement, par exemple, que la moyenne doit être proche de ). Veuillez également noter qu'il existe des restrictions (fortes) sur , sinon, une fois que dépasse (ce qu'il finit par faire), rien ne garantit que reste inférieur ou égal à . F f kF fkp(n)1F11/k=33318kp(n)1F1
whuber

Réponses:

9

Dans un sens, ce que vous avez fait est de caractériser toutes les distributions à valeur entière non négatives.

Mettons de côté la description du processus aléatoire pendant un moment et concentrons-nous sur les récursions de la question.

Si , alors certainement . Si nous réécrivons cette deuxième récursivité en termes de fonction de survie (où a la distribution ), nous obtenons quelque chose de très suggestif et facile à manipuler. Clairement, et donc Ainsi, tant que notre séquence prend des valeurs dansfn=pn(1Fn1) S n = 1 - F n = P ( T > n ) T F S n = 1 - F n = ( 1 - p n ) S n - 1Fn=pn+(1pn)Fn1 Sn=1Fn=P(T>n)TFS n = n k = 0 ( 1 - p k )

Sn=1-Fn=(1-pn)Sn-1,
( p n ) [ 0 , 1 ] n
Sn=k=0n(1-pk).
(pn)[0,1]et ne converge pas trop rapidement vers zéro, nous obtiendrons alors une fonction de survie valide (c'est-à-dire décroissant monotone à zéro en tant que ).n

Plus précisement,

Proposition : Une séquence prenant des valeurs dans [ 0 , 1 ] détermine une distribution sur les entiers non négatifs si et seulement si - n = 0 log ( 1 - p n ) = (pn)[0,1] Et toutes ces distributions ont une séquence correspondant (même si elle peut ne pas être unique).

-n=0Journal(1-pn)=,

(pn)[0,1]

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

Mais attendez, il y a plus!

FF

h(t)=F(t)S(t).

Λ(t)=0th(s)s

S(t)=exp(-Λ(t))=exp(-0th(s)s).
hh(t)0t0th(s)st

t>t0

S(t)=e-t0th(s)sS(t0).

h(t)S(t)

Se reconnecter au boîtier discret

S(n)

h(t)=hn=-Journal(1-pn),
(n-1,n](pn)

pn-Journal(1-pn)pn=Fn/Sn-1

pn=knFnn=k-1Fn=0n>k-1

cardinal
la source
1
kk=1/2F=(0,1/2,1/2,0,)k=1/mF(m+1)=F(m+2)==0
2
fnFnfn
2
Très bonne réponse. C'est très perspicace. J'étais vraiment intéressé de voir ce problème lié à d'autres domaines et concepts.
jbarlow
1
@jbarlow: Merci. Je suis content que vous l'ayez trouvé utile! J'ai aimé y réfléchir un peu, car c'est une belle question.
cardinal
8

p(n)=p<1

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

a la solution

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

p(n)

Autres cas:

  1. p(n)=pn;p<1;F(0)=p
    F(n)=1-(1-p)Γ(n+1-p)Γ(1-p)Γ(n+1)
  2. S(n)=1-F(n)
    S(n)=(1-p(n))S(n-1)
  3. p(n)np(n)=knp>1
    p(n)=1-(1-p)n+1p<1
    p(0)=pp(n)
    F(n)=1-(1-p)n+1n!
    S(n)=1F(n)=(1p)n+1n!
    i=0S(i)=E[N]
    E[N]=(1p)e(1p)
Cam.Davidson.Pilon
la source
2
Cam, ce n'est pas la fonction de hasard, mais plutôt la fonction de survie . :-)
Cardinal
1
ty, * modifié pour la survie
Cam.Davidson.Pilon