Cette distribution a-t-elle un nom? Ou quel est un processus stochastique qui pourrait le générer?

9

Une distribution discrète avec fonction de masse

p(x;k)=k(x+k)(x+k1),x=1,2,

se trouve à la page 9 de ce document .

Pour c'est une distribution de Yule-Simon avec ρ = 1 , mais je n'ai trouvé aucun autre exemple.k=1ρ=1

At-il un nom? Apparaît-il dans d'autres contextes? Existe-t-il un processus stochastique simple qui pourrait le générer?

Simon Byrne
la source

Réponses:

13

C'est une loi de pouvoir discrète.

(Il s'agit d'une description - dont le sens sera précisé ci-dessous - plutôt que d'un terme technique. L'expression "loi de puissance discrète" a une signification technique légèrement différente, comme indiqué par @Cardinal dans les commentaires de cette réponse.)

Pour voir cela, notez que la décomposition partielle des fractions peut s'écrire

p(x;k)=k(x+k)(x+k1)=11+(x1)/k11+x/k.

Les télescopes CDF sous une forme fermée:

CDF(i)=x=1ip(x;k)=[11+0/k11+1/k]+[11+1/k11+2/k]++[11+(i1)/k11+i/k]=11+0/k+[11+1/k+11+1/k]+[11+2/k++11+(i1)/k]11+i/k=1+0++011+i/k=ii+k.

( Soit dit en passant , car cela est facilement inversé, il fournit immédiatement un moyen efficace de générer des variables aléatoires à partir de cette distribution: il suffit de calculer uest uniformément réparti sur(0,1).)ku1uu(0,1)

La différenciation de cette expression par rapport à montre comment le CDF peut être écrit comme une intégrale,i

CDF(i)=ii+k=0idt/k(1+t/k)2=x=1ix1xdt/k(1+t/k)2,

D'où

p(x;k)=x1xdt/k(1+t/k)2.

k

f(ξ)dξ=(1+ξ)2dξ

p(x;k) fkx1x2

whuber
la source
p(x;k)kx2xp(x;k)x2/k1x
p
Je ne suis pas sûr de la distinction que vous essayez de faire. Malheureusement, je n'ai pas eu l'occasion d'y réfléchir attentivement, mais il semble que vous définissiez une distribution de loi de puissance discrète comme une version discrétisée d'une distribution de loi de puissance continue. Suis-je en train d'interpréter correctement votre commentaire? En tout cas, quand je vois des références à des lois de puissance discrètes dans la littérature, la définition habituelle semble être la plus faible (c'est-à-dire asymptotique) que j'ai utilisée. (suite)
cardinal
(Suite) D'un autre côté, une distribution Zipf semble être aussi pure que possible d'une loi de puissance discrète, mais je ne pense pas qu'elle puisse être générée comme une discrétisation d'une loi de puissance continue. Ai-je mal interprété votre intention? (Soit dit en passant, votre développement ci-dessus est assez agréable. La reconnaissance de la somme télescopique pour le cdf est excellente, tout comme la reconnaissance d'un schéma d'échantillonnage facile.)
Cardinal
10

D'accord, après un peu plus d'enquête, j'ai trouvé plus de détails.

PBeta(1,k)
X|PGeometric(P)
Y=X+1

Il a quelques autres propriétés intéressantes:

  • Il a une moyenne infinie
  • XkXt|X>tt+k
Simon Byrne
la source