Respuestas:
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 r 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).
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:
> 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 ghyper
familia (distribuida hipergeométrica) de distribuciones del SuppDists
paquete:
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 ghyper
en el SuppDists
paquete. Creo que esto también se puede encontrar en el libro de Johnson & al . Distribuciones discretas univariantes .