Questa distribuzione ha un nome? O cos'è un processo stocastico che potrebbe generarlo?


9

Una distribuzione discreta con funzione di massa

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

sorge a pagina 9 di questo documento .

Per k=1 è una distribuzione di Yule-Simon con ρ=1 , ma non ho trovato altri esempi.

ha un nome? Appare in altri contesti? Esiste un semplice processo stocastico che potrebbe generarlo?

Risposte:


13

È una legge di potere discreta.

(Questa è una descrizione - il cui significato sarà precisato di seguito - piuttosto che un termine tecnico. La frase "legge del potere discreto" ha un significato tecnico leggermente diverso, come indicato da @Cardinal nei commenti a questa risposta.)

Per vedere questo, osserva che la decomposizione della frazione parziale può essere scritta

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

Il CDF si telescopica in una forma chiusa:

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.

(Per inciso, perché questo è facilmente invertita, fornisce immediatamente un modo efficiente per generare variabili aleatorie di questa distribuzione: è sufficiente calcolare dove è distribuito in modo uniforme su .)ku1uu(0,1)

Differenziare questa espressione rispetto a mostra come il CDF possa essere scritto come integrale,i

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

da dove

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

Questa forma di scrittura mostra come parametro di scala per la famiglia di distribuzioni (continue) determinata dalla densitàk

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

e mostra come è la versione discretizzata di (ridimensionata da ) ottenuta integrando la probabilità continua nell'intervallo da a . Questa è ovviamente una legge di potere con esponente . Questa osservazione ti dà accesso a una vasta letteratura sulle leggi di potere e su come si presentano in scienza, ingegneria e statistica, che può suggerire molte risposte alle tue ultime due domande.p(x;k) fkx1x2


(+1) Dalla funzione di massa della probabilità, è chiaro che come , il che sembra essere sufficiente per concludere che si tratta di una distribuzione di potere-legge. In effetti, come . p(x;k)kx2xp(x;k)x2/k1x
cardinale il

@cardinal Hai ragione, ma c'è un limite a questo argomento: mostra solo che è asintoticamente una legge di potere. I calcoli mostrano che è esattamente una versione discretizzata di una legge di potere. p
whuber

Non sono del tutto sicuro della distinzione che stai cercando di disegnare. Sfortunatamente, non ho avuto la possibilità di pensarci attentamente, ma sembra che tu stia definendo una distribuzione discreta della legge del potere come una versione discretizzata di una distribuzione continua della legge del potere. Sto interpretando correttamente il tuo commento? Ad ogni modo, quando vedo riferimenti a leggi sul potere discreto in letteratura, la definizione abituale sembra essere quella più debole (cioè asintotica) che ho usato. (cont.)
cardinale il

(Cont.) D'altra parte, una distribuzione Zipf sembrerebbe la più pura possibile di una legge di potere discreta, ma non credo che possa essere generata come una discretizzazione di una legge di potere continua. Ho frainteso il tuo intento? (A proposito, il tuo sviluppo sopra è abbastanza bello. Il riconoscimento della somma telescopica per il cdf è ottimo, così come il riconoscimento di un semplice schema di campionamento.)
Cardinale

10

Bene, dopo un po 'più di indagine, ho trovato qualche dettaglio in più.

È un caso speciale di una miscela continua di una distribuzione geometrica con una Beta, quindi potrebbe essere definita una distribuzione geometrica Beta . In particolare, se: e: la distribuzione marginale di ha questa distribuzione. Come tale, è un caso speciale di una distribuzione binomiale beta-negativa .

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

Ha un paio di altre proprietà interessanti:

  • Ha una media infinita
  • Descrive la propria distribuzione di coda: se ha questa distribuzione con il parametro , allora ha il parametro .XkXt|X>tt+k
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.