######################################################################
# Plot 3 normal distributions #
# #
######################################################################
curve(expr = dnorm(x = x, mean = 24.3, sd = 0.6),
xlim = c(20,30), col = "darkgreen", lwd = 2, ylab = "f(y)", xlab = "y")
curve(expr = dnorm(x = x, mean = 24.3, sd = 1.3),
xlim = c(20,30), col = "blue", add = TRUE, lwd = 2)
curve(expr = dnorm(x = x, mean = 23.1, sd = 0.6),
xlim = c(20,30), col = "red", add = TRUE, lwd = 2)
abline(h = 0)
legend(x = 25, y = 0.6, legend = c("mu = 24.3, sigma = 0.6", "mu = 24.3, sigma = 1.3", "mu = 23.1, sigma = 0.6"),
lty = c(1,1,1), col = c("darkgreen", "blue","red"), bty = "n", lwd = c(2,2,2))
#