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