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")
plot
28
May 10
R: More plotting fun with Poission
3
May 10
First annual R plot replication prize
$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()“). - Including the “plot” statement, I created this with 3 short lines of code.
This is based on a random sampling of unstated size, so I don’t expect that your graph will be an absolute, exact match.
I’ll add $1 to the prize for every day that goes by without a winner until the end of the year. After that I’ll consider it an unsolved mystery and reveal the code I used.
Post your guesses for the code as a comment to this post. First correct answer wins. Good luck to all!
21
Apr 10
R: more plotting fun, this time with the Poisson
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")
8
Apr 10
R: another nifty graph
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)) r = list(x=x, y=y) return(r) } r = moxbuller(50000) par(bg="black") par(mar=c(0,0,0,0)) plot(r$x,r$y, pch=".", col="blue", cex=1.2)
25
Mar 10
R plotting fun
Not easy to produce cool looking graphs in R, but it can be done. The results of some messing around are above. Here is the code I used:
x = runif(1000) y = x/runif(1000) cexes = 10*y/max(y) # For circle size par(bg="black") # I see a white background and I want it painted black. par(mar=c(0,0,0,0)) # Margins? We don't kneed no stinkin' margins. plot(x,log(y), pch=20, col="white", cex=cexes)
