在他被广泛引用的论文中,层次模型中方差参数的先验分布 (迄今为止在 Google Scholar 上引用了 916 次)Gelman 提出,分层贝叶斯模型中方差的良好非信息性先验分布是均匀分布和半 t 分布。如果我理解正确,那么当主要关注位置参数(例如平均值)时,它会很好地工作。然而,有时方差参数是主要关注点,例如在分析来自计时任务的人类响应数据时,意味着计时可变性通常是感兴趣的度量。在这些情况下,我不清楚如何使用例如均匀分布对可变性进行分层建模,因为我在分析后希望在参与者级别和组级别上获得均值方差的可信度。
那么我的问题是:当主要关注数据的方差时,在构建分层贝叶斯模型时推荐什么分布?
我知道可以将伽马分布重新参数化为由均值和 SD 指定。例如,下面的层次模型来自 Kruschke 的《做贝叶斯数据分析》一书。但是 Gelman 在他的文章中概述了 gamma 分布的一些问题,我将不胜感激替代方案的建议,最好是在 BUGS/JAGS 中不难工作的替代方案。