我在 WINBUGS 中使用以下模型来运行分层贝叶斯回归,其中beta
是我的协变量:
如果我通过添加以下代码来修改此模型:
# posterior probabilities of Positive beta's
p.beta0 <- step( beta0 )
p.beta1 <- step( beta1 )
然后我可以评估PP
beta 协变量(正或负)关联的后验概率 ( )。
我的 beta 值是:
beta0 = 0.23434
beta1 = -0.4582
使用此代码,PP
ofbeta0
是0.959033
,而PP
ofbeta1
是0.015043
。我的解释beta0
是,这种 covaraite 有 95.9033% 的正相关。但是,我不确定如何解释,beta1
因为这具有负关联和低后验概率。我不确定计算后验概率的代码是否存在问题。
欢迎任何见解。