1. 学ぶ
  2. /
  3. コース
  4. /
  5. RJAGS によるベイズモデリング

Connected

演習

定義・コンパイル・シミュレーション

選挙戦において、\(p\) をあなたを支持する有権者の母集団における割合とします。過去の世論調査と選挙データをもとに構築された \(p\) の事前モデルは、形状パラメータ $a = 45$、\(b = 55\) を持つ Beta(\(a\),\(b\)) 分布です。\(p\) をさらに詳しく把握するため、\(n\) 人の有権者候補にも調査を行いました。あなたを支持する人数 \(X\) の \(p\) への依存関係は、Bin(\(n\),\(p\)) 分布でモデル化されています。

実施した調査では、\(n = 10\) 人中 \(X = 6\) 人があなたを支持しました。次のステップは、この調査結果をもとに \(p\) のモデルを更新することです。そのために rjags パッケージを使って \(p\) の事後モデルを近似します。この演習は rjags の 3 つのステップ、すなわち定義(define)・コンパイル(compile)・シミュレーション(simulate)に分けて進めましょう。

指示1 / 3

undefined XP
    1
    2
    3

ベイズモデルを定義します。

  • \(X\) の尤度モデルを Bin(\(n\),\(p\))、\(p\) の事前分布を Beta(\(a\),\(b\)) として指定します。rjags の構文では、それぞれ dbin(p, n) と dbeta(a, b) で記述します。
  • このモデル文字列を vote_model として保存しましょう。