Pinned, entropy augmented, digitally normal distribution, of no particular work-related use and thus perfectly suitable for today. Code in R: iters = 1000 sd = 2 precision = 20 results = rep(0,it...
https://statisticsblog.com/2014/09/01/labor-day-distribution-fun/
A couple days ago I noticed a fun piece of R code by Allan Roberts, which lets you create a digital snowflake by cutting out virtual triangles. Go give it a try. Roberts inspired me to create a w...
Fantastic map which shows the location of every death on US roads from 2001 to 2009. Go take a look then come back and tell me…. What location did you zoom in on first, and why?
https://statisticsblog.com/2011/11/25/interactive-map-of-us-road-fatalities/
Imagine a unit square. Every side has length 1, perfectly square. Now imagine this square was really a fence, and you picked two spots at random along the fence, with uniform probability over the...
https://statisticsblog.com/2010/07/19/r-clash-of-the-cannon-cycles/
Last night, while playing Boggle, we ended up with a board without a single vowel. Not even a “Y” or “Qu”. This seemed fairly unusual, so I wondered what the chances were of such an occur...
https://statisticsblog.com/2010/07/02/word-games-in-probability-and-r/
What if you wanted to randomly place objects into a field, and the more objects you had, the more they rejected newcomers placed nearby? To find out, I setup a simulation. The code, shown at the ...
https://statisticsblog.com/2010/06/14/repulsive-dots-pattern-the-difference-of-distance/
Coded as follows: x = seq(.001,50,.001) par(bg="black") par(mar=c(0,0,0,0)) plot(x,sin(1/x)*rpois(length(x),x),pch=20,col="blue")
https://statisticsblog.com/2010/05/28/r-more-plotting-fun-with-poission/
$100 to the first person who can figure out how I created this plot and replicate it. Some hints: It was done in R. There is only one underlying probability distribution involved (one “rdist()�...
https://statisticsblog.com/2010/05/03/first-annual-r-plot-replication-prize/
Click on image for a larger version. Here is the code: par(bg="black") par(mar=c(0,0,0,0)) plot(sort(rpois(10000,100))/rpois(10000,100),frame.plot=F,pch=20,col="blue")
https://statisticsblog.com/2010/04/21/r-more-plotting-fun-this-time-with-the-poisson/
Make sure to click on the image to see the large version. Code for this graph: moxbuller = function(n) { u = runif(n) v = runif(n) x = cos(2*pi*u)*sqrt(-2*log(v)) y = sin(2*pi*v)*sqrt(-2*log(u)) ...
https://statisticsblog.com/2010/04/08/r-another-nifty-graph/