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