관측된 데이터를 PyMC의 observed 변수에 합쳐 추론 알고리즘을 실행한다.

 

---

# 베르누이 관측치를 포함한다.
obs = pm.Bernoulli(“obs”, p, value=occurrences, observed=True)
 
# 3장에서 설명할 예정
mcmc = pm.MCMC([p, obs])
mcmc.sample(20000, 1000)

---

 

 

결과

[—————–100%—————–] 20000 of 20000 complete in 0.8 sec

---