ggplot2でヒストグラム＋密度関数のオーバープロット

良く使うのでメモ。正規分布（平均１）だとこんな感じ。

library(ggplot2)
x <- rnorm(1000)+1
#plot
qplot(x, geom = "blank") + 
  geom_histogram(aes(y = ..density..),fill="dodgerblue",colour="black") + 
  stat_function(fun=dnorm,color="brown4",size=2,arg=list(mean=1))