## Not run: m <- 0.8; x <- seq(-1, 4, len = 501)
# plot(x, dposnorm(x, m = m), type = "l", las = 1, ylim = 0:1,
# ylab = paste("posnorm(m = ", m, ", sd = 1)"), col = "blue",
# main = "Blue is density, orange is cumulative distribution function",
# sub = "Purple lines are the 10,20,...,90 percentiles")
# abline(h = 0, col = "grey")
# lines(x, pposnorm(x, m = m), col = "orange", type = "l")
# 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
# ## End(Not run)
Run the code above in your browser using DataLab