N <- 10; shape <- 0.5; y <- 1:N
proby <- dzipf(y, N = N, shape = shape)
if (FALSE) plot(proby ~ y, type = "h", col = "blue",
ylim = c(0, 0.2), ylab = "Probability", lwd = 2, las = 1,
main = paste0("Zipf(N = ", N, ", shape = ", shape, ")"))
sum(proby) # Should be 1
max(abs(cumsum(proby) - pzipf(y, N = N, shape = shape))) # 0?
Run the code above in your browser using DataLab