# NOT RUN {
ddiffzeta(1:20, 0.5, start = 2)
rdiffzeta(20, 0.5)
# }
# NOT RUN {
shape <- 0.8; x <- 1:10
plot(x, ddiffzeta(x, shape = shape), type = "h", ylim = 0:1,
sub = "shape=0.8", las = 1, 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