Zanegbin {VGAM} | R Documentation |
Density, distribution function, quantile function and random generation
for the zero-altered negative binomial distribution with parameter p0
.
dzanegbin(x, p0, size, prob=NULL, munb=NULL, log = FALSE) pzanegbin(q, p0, size, prob=NULL, munb=NULL) qzanegbin(p, p0, size, prob=NULL, munb=NULL) rzanegbin(n, p0, size, prob=NULL, munb=NULL)
x, q |
vector of quantiles. |
p |
vector of probabilities. |
n |
number of observations. Must be a single positive integer. |
size, prob, munb, log |
Parameters from the ordinary negative binomial distribution
(see dnbinom ).
Some arguments have been renamed slightly.
|
p0 |
Probability of zero, called p0.
The default value of p0=0 corresponds
to the response having a positive negative binomial distribution.
|
The probability function of Y is 0 with probability p0
,
else a positive negative binomial(munb, size)
distribution.
dzanegbin
gives the density and
rzanegbin
generates random deviates.
The argument p0
is recycled to the required length, and
must have values which lie in the interval [0,1].
Thomas W. Yee
munb = 3 size = 4 p0 = 0.3 x = (-1):7 (i = dzanegbin(x, p0=p0, munb=munb, size=size)) table(rzanegbin(100, p0=p0, munb=munb, size=size)) ## Not run: x = 0:10 barplot(rbind(dzanegbin(x, p0=p0, munb=munb, size=size), dnbinom(x, mu=munb, size=size)), beside = TRUE, col = c("blue","green"), main=paste("ZANB(p0=", p0, ", munb=", munb, ", size=", size, ") (blue) vs", " NB(mu=", munb, ", size=", size, ") (green)", sep=""), names.arg = as.character(x)) ## End(Not run)