¿Cómo muestrear desde una distribución discreta en los enteros no negativos?


9

Tengo la siguiente distribución discreta, donde son constantes conocidas:α,β

pag(X;α,β)=Beta(α+1,β+X)Beta(α,β)para X=0 0,1,2,...

¿Cuáles son algunos enfoques para tomar muestras de manera eficiente de esta distribución?

Respuestas:


9

Esta es una distribución binomial Beta negativa , con el parámetro en su caso, utilizando la notación de Wikipedia. También llamó distribución Beta-Pascal cuando rr=1r es un número entero. Como señaló en un comentario, esta es una distribución predictiva en el modelo binomial negativo bayesiano con un Beta conjugado antes de la probabilidad de éxito.

Por lo tanto, puede muestrearlo muestreando una variable u y luego muestreando una variable binomial negativa NB ( r , u ) (con r = 1 en su caso, es decir, una distribución geométrica).Beta(α,β)tunótese bien(r,tu)r=1

Esta distribución se implementa en el paquete R brr. El muestreador tiene nombre rbeta_nbinom, el pmf tiene nombre dbeta_nbinom, etc. Las anotaciones son , c = α , d = β . Cheque:una=rC=αre=β

> Alpha <- 2; Beta <- 3
> a <- 1
> all.equal(brr::dbeta_nbinom(0:10, a, Alpha, Beta), beta(Alpha+a, Beta+0:10)/beta(Alpha,Beta))
[1] TRUE

Mirando el código, uno puede ver que en realidad llama a la ghyperfamilia (distribuida hipergeométrica) de distribuciones del SuppDistspaquete:

brr::rbeta_nbinom
function(n, a, c, d){
  rghyper(n, -d, -a, c-1)
}

Ineed, la distribución BNB se conoce como distribución hipergeométrica generalizada de tipo IV . Vea la ayuda de ghyperen el SuppDistspaquete. Creo que esto también se puede encontrar en el libro de Johnson & al . Distribuciones discretas univariantes .


Esta respuesta es excelente, pero sería aún mejor si probaras que la densidad OP publicada es la misma que la densidad binomial negativa.
Sycorax dice Reinstate Monica el

1
@ user777 Creo que el autor de la OP lo demostró él mismo, en vista de su comentario a la respuesta de Xian (distribución predictiva posterior en el modelo binomial negativo con un Beta conjugado anterior).
Stéphane Laurent

10

Beta(α+1,β+X)Beta(α,β)=αα+β+Xβ+X-1α+β+X-1βα+β
XtuU(0 0,1)
Sk=X=0 0kBeta(α+1,β+X)Beta(α,β)
Sk>tu
k
RX=Beta(α+1,β+X)Beta(α,β)=αα+β+Xβ+X-1α+β+X-1βα+β=α+β+X-1α+β+Xβ+X-1α+β+X-1RX-1=β+X-1α+β+XRX-1
Sk=Sk-1+Rk

1
Sk=1-Γ(una+si)Γ(si+k+1)Γ(si)Γ(una+si+k+1)

1
ktuαβtuΓ(si+k+1)Γ(una+si+k+1)αβson integrales, entonces la solución es la raíz de un polinomio, pero aun así, usar Gamma puede ser el camino a seguir.
whuber

1
pag
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.