m = 0.8; x = seq(-1, 4, len = 501)
plot(x, dposnorm(x, m = m), type = "l", ylim = 0:1, las = 1,
ylab = paste("posnorm(m = ", m, ", sd = 1)"), col = "blue",
main = "Blue is density, red is cumulative distribution function",
sub = "Purple lines are the 10,20,...,90 percentiles")
lines(x, pposnorm(x, m = m), col = "red")
abline(h = 0)
probs = seq(0.1, 0.9, by = 0.1)
Q = qposnorm(probs, m = m)
lines(Q, dposnorm(Q, m = m), col = "purple", lty = 3, type = "h")
lines(Q, pposnorm(Q, m = m), col = "purple", lty = 3, type = "h")
abline(h = probs, col = "purple", lty = 3)
max(abs(pposnorm(Q, m = m) - probs)) # Should be 0
Run the code above in your browser using DataLab