Distribució de Yule-Simon

En Probabilitat i Estadística, la distribució de Yule-Simon és una distribució de probabilitat discreta anomenada així en honor d'Udny Yule i Herbert A. Simon. Simon originalment la va anomenar distribució de Yule.

La funció densitat de probabilitat (fdp) de la distribució de Yule-Simon (ρ) és:

$f(k;\rho )=\rho \mathrm{B}(k,\rho +1)$,

per un enter k ≥ 1 i ρ > 0 real, on B és la funció beta. De manera equivalent la fdp es pot escriure en termes de factorials com:

$f(k;\rho )=\frac{\rho \mathrm{\Gamma }(\rho +1)}{(k+\rho {)}^{\underset{\_}{\rho +1}}}$

on $\mathrm{\Gamma }$ es la funció gamma. Així si ρ és un enter:

$f(k;\rho )=\frac{\rho \rho !(k-1)!}{(k+\rho )!}$.

Gràfic de la distribució de Yule–Simon (vermell) i l'asimptòtica llei de Zipf (blau)

El paràmetre ρ es pot estimar utilitzant un algoritme de punt fix. La funció densitat de probabilitat $f(k;\rho )$ té la propietat que per a k prou gran tenim

$f(k;\rho )\approx \frac{\rho \mathrm{\Gamma }(\rho +1)}{k^{\rho +1}}\propto \frac{1}{k^{\rho +1}}$.

Això vol dir que la cua de la distribució de Yule-Simon compleix en realitat la llei de Zipf, i $f(k;\rho )$ es pot utilitzar com a model, per exemple, de la freqüència relativa de la k-èsima paraula més freqüent en un text llarg, que segons la llei de Zipf és inversament proporcional a una (generalment petita) potència de k.

Ocurrència

La distribució de Yule-Simon va sorgir originalment com la distribució limitant d'un procés estocàstic particular, estudiat per Yule com a model per a la distribució dels tàxons biològics i els subtàxons. Simon va nomenar aquest procés el "procés de Yule", però avui és més conegut com un procés de vinculació preferencial. El procés de vinculació preferencial és un procés urna en el que s'afegeixen les boles a un nombre creixent d'urnes, cada bola que s'assigna a una urna amb probabilitat lineal al nombre de l'urna que ja en conté.

La distribució també apareix com una distribució composta, en la que el paràmetre d'una distribució geomètrica es tracta com una funció de la variable aleatòria que té una distribució exponencial. Específicament, se suposa que W segueix una distribució exponencial amb l'escala 1/ρ o ratio ρ:

$W\sim \mathrm{Exponential}(\rho )$,

amb funció densitat

$h(w;\rho )=\rho \exp (-\rho w)$..

Llavors una variable K que segueixi la distribució de Yule-Simon té la següent distribució geomètrica condicional en W:

$K\sim \mathrm{Geometric}(\exp (-W)).$

La fdp de una distribució geomètrica és:

$g(k;p)=p(1-p{)}^{k-1}$

per a $k\in \{1,2,\dots \}$. La funció densitat de probabilitat és la següent distribució composta exponencial-geomètrica:

$f(k;\rho )=\underset{0}{\overset{\mathrm{\infty }}{\int }}g(k;\exp (-w))h(w;\rho )dw$.

amb la següent relació de recurrència:

$\left\{\begin{array}{c}kP(k)=(\alpha +k+1)P(k+1),\\ P(1)=\alpha B(\alpha +1,1)\end{array}\right\}$

Generalització

La generalització dels dos paràmetres de la distribució de Yule porten a reemplaçar la funció beta per una funció beta incompleta. La funció densitat de probabilitat de la distribució de Yule-Simon generalitzada (ρ, α) es defineix com:

$f(k;\rho ,\alpha )=\frac{\rho }{1-{\alpha }^{\rho }}{\mathrm{B}}_{1-\alpha }(k,\rho +1)$

amb 0 ⩽ α < 1. Per a α = 0 s'obté com un cas particular la distribució ordinària de Yule-Simon (ρ). L'ús de la funció beta incompleta té l'efecte d'introduir un límit exponencial de la cua superior.

Bibliografia

• Colin Rose and Murray D. Smith, "Mathematical Statistics with Mathematica". New York: Springer, 2002, Plantilla:ISBN