Draw normal distribution on histogram of any dataset

hist_and_norm.R

# Draw a normal distribution curve on top of a histogram
# of any vector
#
# @param data A numeric vector

hist_and_norm <- function(data){
  hist_min <- mean(data) - 4*sd(data)
  hist_max <- mean(data) + 4*sd(data)
  normalx <- seq(hist_min, hist_max, by = 1)
  normaly <- dnorm(normalx, mean = mean(data), sd = sd(data))
  hist(rvs, xlim = c(hist_min, hist_max), ylim =c(0, max(normaly)),
       freq = FALSE)
  lines(normalx, normaly, col = "red")
}

rvs <- sample(1:1000, replace = TRUE)
hist_and_norm(rvs)