set.seed(1234)
x <- rbinom(20, 1, 0.4) # Structural zeroes
y <- rbinom(20, 1, 0.7*x) # Non-structural zeroes
fit <- bayesZIB(y~1|1, priors=list(c(0,0.5), c(0.5,1)))
print(fit$fit, pars=c("theta", "beta"))
Run the code above in your browser using DataLab