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.