In [17]:
import slatkin
In [10]:
slatkin.montecarlo(100000, [10, 2, 3, 1,1,1])
Out[10]:
In [11]:
counts = [8, 4, 3, 2, 2, 1, 1, 1, 1, 1, 1]
In [18]:
(prob, theta) = slatkin.montecarlo(100000, counts
)
print "prob: %s theta: %s" % (prob, theta)
In [16]:
%reset
In [ ]: