munb = 3; phi = 0.2; k = 10
x = 0:10
(i = dzinb(x, phi=phi, mu=munb, size=k))
cumsum(i) - pzinb(x, phi=phi, mu=munb, size=k) # Should be 0s
table(rzinb(100, phi=phi, mu=munb, size=k))
table(qzinb(runif(1000), phi=phi, mu=munb, size=k))
round(dzinb(x, phi=phi, mu=munb, size=k) * 1000) # Should be similar
barplot(rbind(dzinb(x, phi=phi, mu=munb, size=k),
dnbinom(x, mu=munb, size=k)),
beside = TRUE, col = c("blue","green"),
main=paste("ZINB(mu=", munb, ", k=", k, ", phi=", phi,
") (blue) vs negative binomial(mu=", munb,
", k=", k, ") (green)", sep=""),
names.arg = as.character(x))
Run the code above in your browser using DataLab