## Not run:
# shape <- 0.5; Scale <- 1; nn <- 501
# x <- seq(-0.10, 3.0, len = nn)
# plot(x, dgenray(x, shape, scale = Scale), type = "l", las = 1, ylim = c(0, 1.2),
# ylab = paste("[dp]genray(shape = ", shape, ", scale = ", Scale, ")"),
# col = "blue", cex.main = 0.8,
# main = "Blue is density, orange is cumulative distribution function",
# sub = "Purple lines are the 10,20,...,90 percentiles")
# lines(x, pgenray(x, shape, scale = Scale), col = "orange")
# probs <- seq(0.1, 0.9, by = 0.1)
# Q <- qgenray(probs, shape, scale = Scale)
# lines(Q, dgenray(Q, shape, scale = Scale), col = "purple", lty = 3, type = "h")
# lines(Q, pgenray(Q, shape, scale = Scale), col = "purple", lty = 3, type = "h")
# abline(h = probs, col = "purple", lty = 3)
# max(abs(pgenray(Q, shape, scale = Scale) - probs)) # Should be 0
# ## End(Not run)
Run the code above in your browser using DataLab