###################################################################### # 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") #