pobs0
.dzabinom(x, size, prob, pobs0 = 0, log = FALSE)
pzabinom(q, size, prob, pobs0 = 0)
qzabinom(p, size, prob, pobs0 = 0)
rzabinom(n, size, prob, pobs0 = 0)
length(n) > 1
then the length is taken to be the number required.dbinom
).pobs0 = 0
corresponds
to the response having a positive binomial distribution.dzabinom
gives the density and
pzabinom
gives the distribution function,
qzabinom
gives the quantile function, and
rzabinom
generates random deviates.pobs0
,
else a positive
binomial(size, prob)
distribution.zibinomial
,
rposbinom
.size <- 10; prob = 0.15; pobs0 <- 0.05; x <- (-1):7
dzabinom(x, size = size, prob = prob, pobs0 = pobs0)
table(rzabinom(100, size = size, prob = prob, pobs0 = pobs0))
x = 0:10
barplot(rbind(dzabinom(x, size = size, prob = prob, pobs0 = pobs0),
dbinom(x, size = size, prob = prob)),
beside = TRUE, col = c("blue", "orange"), cex.main = 0.7, las = 1,
ylab = "Probability", names.arg = as.character(x),
main = paste("ZAB(size = ", size, ", prob = ", prob, ", pobs0 = ", pobs0,
") [blue] vs", " Binom(size = ", size, ", prob = ", prob,
") [orange] densities", sep = ""))
Run the code above in your browser using DataLab