Zero-inflated Negative Binomial Model in Rethinking

161 views Asked by At

I am looking for help/resources to guide me on how to change this zero-inflate Poisson model into zero-inflated negative binomial model in Rethinking. I don't see a function like dzipois that would easily take care of this issue. Could I somehow use dgampois in a zero-inflated negative binomial model? Any ideas? Thank you!!!

SitePoisson <- ulam(
  alist(
    SportfishCount ~ dzipois(pbar, lamba), 
    logit(pbar) <- z1 + z2*Dim2
    log(lamba)a + b*Dim1 + c*Dim2 + d*HTOY, 
    a ~ dnorm(1.85, 2.76), 
    b ~ dnorm(1, 100), 
    c ~ dnorm(1, 100), 
    d ~ dnorm(0, 1), 
    c(z1, z2) ~ dnorm(0, 1)
  ), 
  data=ModelData, 
  start=list(a=0, b=0, c=0, d=0, z1=0.1, z2=0.2), 
  iter=5000, warmup=2000, chains=3, cores=4, log_lik=TRUE
)

I've read as much as I could, but I am not finding any information on zero-inflated negative binomial models in bayesian techniques.

0

There are 0 answers