# NOT RUN {
alpha <- 3; k <- exp(1); x <- seq(2.8, 8, len = 300)
# }
# NOT RUN {
plot(x, dpareto(x, scale = alpha, shape = k), type = "l",
main = "Pareto density split into 10 equal areas")
abline(h = 0, col = "blue", lty = 2)
qvec <- qpareto(seq(0.1, 0.9, by = 0.1), scale = alpha, shape = k)
lines(qvec, dpareto(qvec, scale = alpha, shape = k),
col = "purple", lty = 3, type = "h")
# }
# NOT RUN {
pvec <- seq(0.1, 0.9, by = 0.1)
qvec <- qpareto(pvec, scale = alpha, shape = k)
ppareto(qvec, scale = alpha, shape = k)
qpareto(ppareto(qvec, scale = alpha, shape = k),
scale = alpha, shape = k) - qvec # Should be 0
# }
Run the code above in your browser using DataLab