后验概率
您已经使用 RJAGS 的输出探索并量化了后验趋势与不确定性 $b$。您也可以用 RJAGS 的输出来检验具体假设。例如：在平均意义下，身高每增加 1 cm，体重是否有超过 1.1 kg 的增幅？也就是说，\(b > 1.1\) 的后验概率是多少？
您将通过超过 1.1 的 \(b\) 马尔可夫链取值所占比例来"近似"这个概率。包含 100,000 次迭代的马尔可夫链输出的数据框 weight_chains 已在您的工作区中。
本练习是课程的一部分
使用 RJAGS 的贝叶斯建模
练习说明
- 构造 \(b\) 马尔可夫链取值的密度图，并使用
geom_vline()在 1.1 处叠加一条竖直线。 - 使用
table()汇总超过 1.1 的 \(b\) 马尔可夫链取值的个数。 - 使用
mean()计算超过 1.1 的 \(b\) 马尔可夫链取值所占的"比例"。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Mark 1.1 on a posterior density plot for b
ggplot(___, aes(x = ___)) +
geom_density() +
geom_vline(xintercept = ___, color = "red")
# Summarize the number of b chain values that exceed 1.1
# Calculate the proportion of b chain values that exceed 1.1