x = seq(0, 14, by = 0.01); d = 1.5; Scale = 2; k = 6
plot(x, dgengamma(x, Scale, d, k), type = "l", col = "blue", ylim = 0:1,
main = "Blue is density, red is cumulative distribution function",
sub = "Purple are 5,10,...,95 percentiles", las = 1, ylab = "")
abline(h = 0, col = "blue", lty = 2)
lines(qgengamma(seq(0.05,0.95,by = 0.05), Scale, d, k),
dgengamma(qgengamma(seq(0.05,0.95,by = 0.05), Scale, d, k),
Scale, d, k), col = "purple", lty = 3, type = "h")
lines(x, pgengamma(x, Scale, d, k), type = "l", col = "red")
abline(h = 0, lty = 2)
Run the code above in your browser using DataLab