估计 zoib beta 回归模型的 DIC

Question

我有下面的数据和代码，非常感谢您帮助估计

DIC

（或

AIC

）的 beta 回归模型，其中

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

（或

rd.md

），我将不胜感激。

Answer 1

你有没有弄清楚这一点？我遇到了同样的问题，不知道如何继续。