https://rud.is/b/2012/03/04/thinkstats-in-r-including-example-1-2/