背景
最常用的弱先验方差之一是带参数的反伽马 (格尔曼 2006 年)。
但是,此分布的 90%CI 约为.
library(pscl)
sapply(c(0.05, 0.95), function(x) qigamma(x, 0.001, 0.001))
[1] 3.362941e+19 Inf
由此,我解释为给出方差很大的概率很低,方差小于 1 的概率很低.
pigamma(1, 0.001, 0.001)
[1] 0.006312353
问题
我错过了什么还是这实际上是一个信息丰富的先验?
更新澄清,我考虑这个“信息性”的原因是因为它非常强烈地声称方差是巨大的,并且远远超出了几乎任何测量过的方差的规模。
后续对大量方差估计的荟萃分析会提供更合理的先验吗?
参考
Gelman 2006。层次模型中方差参数的先验分布。贝叶斯分析 1(3):515–533