我正在使用runjags,我正在尝试为我的模型获取DIC.

我知道如果我打印runjags对象,它将返回包含DIC的模型摘要,但我的模型包含很多变量,因此计算摘要统计数据需要很长时间.当我使用summarize = FALSE时,不返回DIC.

我尝试过使用：extract(runJagsObj,what =’dic’)

但它返回此错误消息：
“错误：请提供列名”

有没有办法在摘要信息中查看DIC而无需计算每个变量的汇总统计信息？或者为什么提取物不起作用？

如果这是有用的信息我正在使用method = parallel.

最佳答案 我终于找到了“错误：请提供列名”错误.

我正在使用ggmcmc进行绘图,这需要tidyr. Tidyr与提取物有命名重叠.我能够通过在解决错误的ggmcmc之后加载runjags来解决命名冲突.