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