JAGS版本更新后贝叶斯分层模型结果不同的原因问询
贝叶斯分层建模结果复现疑问
- 使用R语言
R2jags工具开展贝叶斯分层建模,运行经Windows资源管理器“修改日期”验证的、一个月前的代码与数据集,得到的结果与之前不同。唯一变量为更换新工作电脑并安装JAGS 4.3.0,此前使用JAGS 4.2.0。 - 仅因JAGS版本更新就导致结果不同是否有可能?所有条件均一致,无需故障排查协助。
- 补充信息:模型收敛情况良好——Gewekes检验、自相关图及迹图表现正常,与之前无差异。
- 已通过
set.seed()和jags.seed设置随机种子,这是否足够?此前从未遇到结果无法复现的问题。 - 结果差异足以影响推断结论:正在评估336名研究对象中30种化学暴露与健康结局的关联,其中Chemical B的可信区间偏移问题最为突出,Chemical A是典型例子(附图:化学物关联结果对比)。
- 将迭代次数从5万增加至10万,仅产生微小且无意义的差异。
- 补充2:已在JAGS官方讨论区发帖询问版本间默认随机数生成器(RNGs)的差异。
内容的提问来源于stack exchange,提问作者pakalla
相关产品推荐
相关产品推荐

