馬可夫鏈的密度圖
追蹤圖會呈現馬可夫鏈的縱向行為，而密度圖則顯示鏈值的最終分佈。密度圖進而提供後驗模型的近似。你將為下方的 \(m\) 馬可夫鏈建立並檢視密度圖。mcmc.list 物件 sleep_sim 與 sleep_chains 資料框已在你的工作區：
sleep_sim <- coda.samples(model = sleep_jags, variable.names = c("m", "s"), n.iter = 10000)
sleep_chains <- data.frame(sleep_sim[[1]], iter = 1:10000)
本練習屬於課程
使用 RJAGS 的貝氏建模
練習說明
對
sleep_sim使用plot()，並設定trace = FALSE，以建立 \(m\) 與 \(s\) 鏈的密度圖。對
sleep_chains使用ggplot()，重新建立 \(m\) 鏈的密度圖。
動手互動練習
試著完成這個範例程式碼，體驗一下這個練習。
# Use plot() to construct density plots of the m and s chains
# Use ggplot() to construct a density plot of the m chain
ggplot(___, aes(x = ___)) +
___()