python - PyMCモデルで状態はどこに保持されますか?

Question

次のモデルが与えられた場合、私の質問は、 、、およびSについて何かを知る方法です。別のファイル (つまり、Python モジュール) で指定されたモデルが指定されている例を見たことがありますが、これはある程度理にかなっています。しかし、ここでは明示的にデータを渡していません。これがどのように機能するかを理解しようとしています。alphabetathetaMCMCS

import pymc as pm
import numpy as np

n = 5 * np.ones(4, dtype=int)
x = np.array([-0.86, -0.3, -0.05, 0.73])

alpha = pm.Normal('alpha', mu=0, tau=0.01)
beta = pm.Normal('beta', mu=0, tau=0.01)

@pm.deterministic
def theta(a=alpha, b=beta):
    """theta = logit^{-1}(a+b)"""
    return pm.invlogit(a + b * x)

d = pm.Binomial('d',n=n,p=theta,value=np.array([0.,1.,3.,5.]),observed=True)

S = pm.MCMC(db='pickle')
S.sample(iter=10000, burn=5000, thin=2)
pm.Matplot.plot(S)