Last chance! 50% off unlimited learning
Sale ends in
pobs0
.dzapois(x, lambda, pobs0 = 0, log = FALSE)
pzapois(q, lambda, pobs0 = 0)
qzapois(p, lambda, pobs0 = 0)
rzapois(n, lambda, pobs0 = 0)
length(n) > 1
then the length is taken to be the number required.pobs0 = 0
corresponds
to the response having a positive Poisson distribution.dzapois
gives the density,
pzapois
gives the distribution function,
qzapois
gives the quantile function, and
rzapois
generates random deviates.pobs0
, else a positive $Poisson(\lambda)$.zapoisson
,
dzipois
.lambda <- 3; pobs0 <- 0.2; x <- (-1):7
(ii <- dzapois(x, lambda, pobs0))
max(abs(cumsum(ii) - pzapois(x, lambda, pobs0))) # Should be 0
table(rzapois(100, lambda, pobs0))
table(qzapois(runif(100), lambda, pobs0))
round(dzapois(0:10, lambda, pobs0) * 100) # Should be similar
x <- 0:10
barplot(rbind(dzapois(x, lambda, pobs0), dpois(x, lambda)),
beside = TRUE, col = c("blue", "green"), las = 1,
main = paste("ZAP(", lambda, ", pobs0 = ", pobs0, ") [blue] vs",
" Poisson(", lambda, ") [green] densities", sep = ""),
names.arg = as.character(x), ylab = "Probability")
Run the code above in your browser using DataLab