---

# 베르누이 확률변수를 통해 확률과 관측치를 연결한다.
observed = pm.Bernoulli(“bernoulli_obs”, p, value=D, observed=True)
 
model = pm.Model([observed, beta, alpha])
 
# 3장에서 설명할 예정
map_ = pm.MAP(model)
map_.fit()
mcmc = pm.MCMC(model)
mcmc.sample(120000, 100000, 2)

---

 

 

결과

[—————–100%—————–] 120000 of 120000 complete in 6.9 sec

---