我正在使用 runjags 并试图为我的模型获取 DIC。

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

我试过使用:extract(runJagsObj, what='dic')
但它返回此错误消息:"Error: Please supply column name"
有没有办法在摘要信息中查看 DIC 而不必计算每个变量的摘要统计数据?或者为什么提取物不起作用?

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

最佳答案

我终于弄清楚了 "Error: Please supply column name" 错误。

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

关于r - 无法提取 DIC,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/31838807/

10-12 21:16