######################################################################
# One-parameter beta PDF #
# #
######################################################################
par(xaxs = "i", yaxs = "i")
curve(expr = dbeta(x = x, shape1 = 1, shape2 = 1), xlim = c(0,1), ylim = c(0,10), col = "red",
ylab = "f(y)", xlab = "y")
curve(expr = dbeta(x = x, shape1 = 1, shape2 = 2), col = "green", add = TRUE)
curve(expr = dbeta(x = x, shape1 = 1, shape2 = 5), col = "yellow", add = TRUE)
curve(expr = dbeta(x = x, shape1 = 1, shape2 = 10), col = "blue", add = TRUE)
legend(x = 0.4, y = 6, legend = c("beta = 1", "beta = 2", "beta = 5", "beta = 10"), lty = 1,
col = c("red", "green", "yellow", "blue"), bty = "n")
#