定義、編譯與模擬

在你的選戰情境中，設 $p$ 為支持你的潛在選民比例。根據先前的民調與選舉資料，你對 $p$ 的「先驗」模型為 Beta($a$,$b$)，其形狀參數為 $a = 45$、$b = 55$。為了進一步了解 $p$，你也對 $n$ 位可能投票者進行了調查。$X$（支持你的受訪者人數）對 $p$ 的依賴關係以 Bin($n$,$p$) 分配建模。

在最終的民調中，$n = 10$ 位受訪者中有 $X = 6$ 位支持你。接下來的目標是根據這些觀察到的民調資料來更新你對 $p$ 的模型！為此，你將使用 rjags 套件來近似 $p$ 的「後驗」模型。我們把這個練習分成 3 個 rjags 步驟：define（定義）、compile（編譯）、simulate（模擬）。

本練習屬於課程

使用 RJAGS 的貝氏建模

檢視課程

動手互動練習

試著完成這個範例程式碼，體驗一下這個練習。

# DEFINE the model
vote_model <- "model{
    # Likelihood model for X
    X ~ ___
    
    # Prior model for p
    p ~ ___
}"

編輯並執行程式碼