情報的事前分布を使用するようにモデルを変更する
表示されているコードは、第2章でゼロから構築した古いモデルです。
この演習はコースの一部です
RによるベイズデータAnalysisの基礎
演習の手順
- このモデルを変更して、
proportion_clicksに次の新しい情報的事前分布を使用しましょう。
rbeta(n_draws, shape1 = 5, shape2 = 95)
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
n_draws <- 100000
n_ads_shown <- 100
# Change the prior on proportion_clicks
proportion_clicks <-
runif(n_draws, min = 0.0, max = 0.2)
n_visitors <-
rbinom(n_draws, size = n_ads_shown,
prob = proportion_clicks)
prior <-
data.frame(proportion_clicks, n_visitors)
posterior <-
prior[prior$n_visitors == 13, ]
# This plots the prior and the posterior in the same plot
par(mfcol = c(2, 1))
hist(prior$proportion_clicks,
xlim = c(0, 0.25))
hist(posterior$proportion_clicks,
xlim = c(0, 0.25))