始める無料で始める

事後分布の更新

選挙支持率 \(p\) の事後モデルは、\(p\) の事前モデルと世論調査データ \(X\) の両方をもとに構築されます。右のスクリプトを実行して、元の事前分布(Beta(45, 55))と元の調査データ(\(n = 10\) 人中 \(X = 6\) 人が支持)から得られた事後分布を確認しましょう。定義済みの vote_model はワークスペースに用意されています。

3ステップの演習を通じて、異なる事前モデルを使用した場合や新しいデータを観測した場合(あるいはその組み合わせ)に、事後分布がどのように変化するかを探っていきましょう。

この演習はコースの一部です

RJAGS によるベイズモデリング

コースを見る

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# COMPILE the model    
vote_jags <- jags.model(textConnection(vote_model), 
    data = list(a = 45, b = 55, X = 6, n = 10),
    inits = list(.RNG.name = "base::Wichmann-Hill", .RNG.seed = 100))

# SIMULATE the posterior
vote_sim <- coda.samples(model = vote_jags, variable.names = c("p"), n.iter = 10000)

# PLOT the posterior
plot(vote_sim, trace = FALSE, xlim = c(0,1), ylim = c(0,18))
コードを編集して実行