Question

我正在使用R包'zoib'在R中执行beta回归。我试图复制包装作者在The R Journal中发表的论文第41页中包含的例子：

Lui F and Kong Y. 2015. zoib: An R Package for Bayesian Inference for Beta Regression and Zero/One Inflated Beta Regression. The R Journal 7(2)

我相信我使用的是他们使用的完全相同的数据和代码：

library(zoib)
data("GasolineYield", package="zoib")
GasolineYield$batch <- as.factor(GasolineYield$batch)
d <- GasolineYield
eg1.fixed <- zoib(yield ~ temp + as.factor(batch) | 1, data=GasolineYield, joint=FALSE,
         random=0, EUID=1:nrow(d), zero.inflation=F, one.inflation=F, 
         n.iter=1050, n.thin=5, n.burn=50)
sample1 <- eg1$coeff
traceplot(sample1)
autocorr.plot(sample1)
gelman.diag(sample1)

但是，当我尝试对样本进行诊断图时，我收到错误。这是错误消息：

Error in ts(seq(from = start(x), to = end(x), by = thin(x)), start = start(x),  : 
  invalid time series parameters specified

我无法理解为什么代码无效或我可以做些什么来解决问题。我可以将错误跟踪到由time调用的zoib函数，似乎问题是示例对象没有tsp属性，但是{ {1}}包作者明确表示他们的模型输出意图与zoib一起使用，所以我很困惑。我没有太多使用MCMC或时间序列对象的经验，所以也许我只是遗漏了一些明显的东西。任何人都可以解释为什么包作者提供的示例失败了，解决方案是什么？

Answer 1

我通过电子邮件发送了包裹作者（刘芳），她告诉我实际上我的软件包版本中有一个错误，但错误是在{{1}的最新版本中修复的（版本1.4.2）。使用最新版本，代码现在可以使用。

使用'zoib'R包进行beta回归时'ts'函数出错

1 个答案: