import pymc as mc

mc.__version__

alpha = 2
beta = 2
n = 20
yes = 15

with mc.Model() as model:
    p = mc.Beta('p',alpha,beta)
    surv = mc.Binomial('surv',n=n,p=p,observed=yes)
    survsim = mc.Binomial('survsim',n=n,p=p)

with model:
    start = mc.find_MAP()

start

with model:
    step = mc.NUTS()

with model:
    trace = mc.sample(10000,step,start)

mc.traceplot(trace);