ddiffzeta(1:20, 0.5, start = 2)
rdiffzeta(20, 0.5)
if (FALSE) shape <- 0.8; x <- 1:10
plot(x, ddiffzeta(x, sh = shape), type = "h", ylim = 0:1, las = 1,
sub = "shape=0.8", col = "blue", ylab = "Probability",
main = "Differenced zeta distribution: blue=PMF; orange=CDF")
lines(x + 0.1, pdiffzeta(x, shape = shape), col = "orange",
lty = 3, type = "h")
Run the code above in your browser using DataLab