微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何仅在rstanarm中的一个特定预测器上指定信息先验

如何解决如何仅在rstanarm中的一个特定预测器上指定信息先验

我目前正在尝试使用 rstanarm 拟合贝叶斯多级模型。 我预计 - 并且有文献证据 - 我的主要预测变量的系数介于 0.15 和 0.65 之间。因此,我想设置一个信息先验,但仅针对此变量,并为其他变量保留弱信息认值。 到目前为止,我有

   stan_glmer(isei_r ~ 1 + maxisei_cntr + agea + as.factor(gender) + as.factor(emp_status) + (1 + maxisei_cntr | country),data = ess,seed = 349,prior = normal(0.40,0.25,autoscale=F))  

但通过这种方式,它在我所有的协变量之前应用了信息。是否可以仅为一个预测变量指定信息先验?

谢谢

解决方法

是的,但您需要传递大小等于系数数量(不包括截距)的先前位置和/或比例的向量。因此,在您的示例中,它可能类似于:

stan_glmer(...,prior = normal(location = c(0,0.4,0),scale = c(10,0.25,10,10))

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。