計算聯合分配
這裡顯示的是我們在上一支影片中撰寫的部分程式碼。它在 pars 中定義了點擊比例(proportion_clicks)與可能結果(n_visitors)的網格。接著為每個參數組合加入其 prior 機率,以及各 proportion_clicks 產生對應 n_visitors 的概似(likelihood)。
本練習屬於課程
R 的 Bayesian 資料分析基礎
練習說明
- 新增欄位
pars$probability:每個proportion_clicks與n_visitors組合的機率。和影片中一樣,作法是將likelihood乘上prior。 - 確認
pars$probability的總和為1.0,也就是將其「正規化」:用pars$probability除以pars$probability的總和。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
n_ads_shown <- 100
proportion_clicks <- seq(0, 1, by = 0.01)
n_visitors <- seq(0, 100, by = 1)
pars <- expand.grid(proportion_clicks = proportion_clicks,
n_visitors = n_visitors)
pars$prior <- dunif(pars$proportion_clicks, min = 0, max = 0.2)
pars$likelihood <- dbinom(pars$n_visitors,
size = n_ads_shown, prob = pars$proportion_clicks)
# Add the column pars$probability and normalize it
pars$probability <- ___