curve

Think I got it: set.seed(1) df <- data.frame(PF = 10*rnorm(1000)) ggplot(df, aes(x = PF)) + geom_histogram(aes(y =..density..), breaks = seq(-50, 50, by = 10), colour = “black”, fill = “white”) + stat_function(fun = dnorm, args = list(mean = mean(df$PF), sd = sd(df$PF)))