你剛做了些貝氏資料分析！

R 的 Bayesian 資料分析基礎

Rasmus Bååth

Data Scientist

先驗與後驗

先驗是一個機率分佈，表示模型在看到資料前的認知。
後驗是一個機率分佈，表示模型在看到資料後的認知。

擲 5 顆骰子時出現 6 的個數之機率分佈

$$\text{p(x)} = {5\choose x}(\frac{1}{6})^x(1-(\frac{1}{6}))^{5-x}$$

number_of_sixes

1 1 0 0 1 0 1 2 0 1 0 0 1 0 0 0 0 0 1 1 1 0 1 1 0 0 2 0
0 1 0 0 1 0 0 1 0 1 2 0 1 0 0 0 1 2 1 2 0 0 1 1 3 3 0 0
1 1 1 1 1 0 0 1 2 0 1 3 1 1 1 0 1 0 1 2 0 1 1 0 1 1 1 0
2 1 0 4 0 1 2 1 1 1 2 0 1 0 1 1 0 0 2 0 0 0 0 0 1 1 0 1
0 0 0 0 2 0 0 0 0 0 1 1 0 0 2 1 1 1 0 2 1 1 1 0 0 1 1 1
...

mean(number_of_sixes)

0.83

posterior <- prop_model(data)
print(posterior)

0.23 0.36 0.20 0.21 0.12
0.10 0.03 0.16 0.09 0.14
0.23 0.05 0.15 0.26 0.22
...

完成殭屍藥物分析！

R 的 Bayesian 資料分析基礎