估算zoib beta回归模型的DIC

Question

我有下面的数据和代码，非常感谢您帮助估算DIC的{​​{1}}（或AIC）与zoib的回归模型：

library(zoib)
data("GasolineYield", package = "zoib")
re.md <- zoib(yield ~ temp | 1 | 1, data=GasolineYield, 
joint = FALSE, random=1, EUID=GasolineYield$batch, 
zero.inflation = FALSE, one.inflation = FALSE, 
n.iter=3200, n.thin=15, n.burn=200)
sample2 <- re.md$coeff
summary(sample2)

根据作者发布的软件包vignette和article，我实施了以下代码以获取DIC，但收到错误消息：

dic.samples(sample2, n.iter=201, thin = 1, type="pD")
Error in nchain(model) : Invalid JAGS model object in nchain

我使用coda.samples {rjags}文档中的对象运行相同的代码并且它有效：

data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
summary(LINE.out)
dic.samples(LINE, n.iter=201, thin = 1, type="pD")

然后我检查了LINE，sample2，re.md和LINE.out上方的每个对象的类，输出为：

class(LINE)
"jags"
class(re.md)
"list"
class(sample2)
"mcmc.list"
class(LINE.out)
"mcmc.list"

这表明我的错误可能是因为对象sample2不属于jags类。

因此，如果我能够sample2以dic.samples可以接受的形式获取DIC以便能够获得AIC（或{{1我的rd.md模型。