#!/usr/bin/R
#first execute in shell: links2 -dump http://www.reddit.com/r/Ingress/comments/25um5v/lets_create_capsule_hacks_statistics/|grep '[0-9]/[0-9]'|sed -e 's/^[^0-9]*//' -e 's/ .*//' -e 's/\([0-9/]\+\).*/\1/'|grep '[0-9]/[0-9]' >dat
read.table("dat",sep="/")->X
probb=function(x,p) {dbinom(x=x[[1]],size=x[[2]],prob=p)}
MLE=function(p) prod(apply(X,1,function(x) probb(x,p)))
q=(1:300)/100000
plot(q,sapply(q,MLE))
print(q[which.max(sapply(q,MLE))])