# Setup
source('https://myweb.uiowa.edu/pbreheny/7110/f25/notes/fun.R')

# Gamma vs Normal figure
col <- c("#FF4E37FF", "#008DFFFF")
x <- seq(0, 8, len = 101)
plot(x, dnorm(x, 4, 2), type = "l", ylim = c(0, 0.25), las = 1, bty = "n", ylab = "Density", lwd = 3, col = col[1])
lines(x, dgamma(x, 4, 1), lwd = 3, col = col[2])
rightlegend(legend = c("Model", "True"), lwd = 3, col = col)