x = seq(-0.2, 4, by=0.01)
loc = 0; Scale = 1; ineq = 1; shape = 1.0;
plot(x, dparetoIV(x, loc, Scale, ineq, shape), type="l", col="blue",
main="Blue is density, red is cumulative distribution function",
sub="Purple are 5,10,...,95 percentiles", ylim=0:1, las=1, ylab="")
abline(h=0, col="blue", lty=2)
Q = qparetoIV(seq(0.05,0.95,by=0.05), loc, Scale, ineq, shape)
lines(Q, dparetoIV(Q, loc, Scale, ineq, shape), col="purple", lty=3, type="h")
lines(x, pparetoIV(x, loc, Scale, ineq, shape), col="red")
abline(h=0, lty=2)
Run the code above in your browser using DataLab