## Not run: rate <- 2; shape <- 0.5; nn <- 201
# x <- seq(-0.05, 1.05, len = nn)
# plot(x, dexppois(x, rate = rate, shape), type = "l", las = 1, ylim = c(0, 3),
# ylab = paste("fexppoisson(rate = ", rate, ", shape = ", shape, ")"),
# col = "blue", cex.main = 0.8,
# main = "Blue is the density, orange the cumulative distribution function",
# sub = "Purple lines are the 10,20,...,90 percentiles")
# lines(x, pexppois(x, rate = rate, shape), col = "orange")
# probs <- seq(0.1, 0.9, by = 0.1)
# Q <- qexppois(probs, rate = rate, shape)
# lines(Q, dexppois(Q, rate = rate, shape), col = "purple", lty = 3, type = "h")
# lines(Q, pexppois(Q, rate = rate, shape), col = "purple", lty = 3, type = "h")
# abline(h = probs, col = "purple", lty = 3); abline(h = 0, col = "gray50")
# max(abs(pexppois(Q, rate = rate, shape) - probs)) # Should be 0
# ## End(Not run)
Run the code above in your browser using DataLab