Shading values/areas under the normal curve is a quite frequent taks in eg educational contexts. Thanks to Hadley in this post, I found this easy solution. library(ggplot2) r ggplot(NULL, aes(c(-3,3))) + geom_area(stat = "function", fun = dnorm, fill = "#00998a", xlim = c(-3, 0)) + geom_area(stat = "function", fun = dnorm, fill = "grey80", xlim = c(0, 3)) plot of chunk unnamed-chunk-1 Simple, right? Some minor beautification: ```r ggplot(NULL, aes(c(-3,3))) + geom_area(stat = "function", fun …