用zoib模型预测(MCMC / RJags)

时间:2018-02-01 22:53:19

标签: r bayesian predict rjags

我正在使用R中的zoib包来构建零膨胀的β回归模型。我正在寻找一种简单的方法来使用zoib生成的模型来计算新数据集的预测响应。通过"新数据集"我的意思是不用于构建原始zoib模型的数据。

我知道我可以只使用 zoib 模型参数并在R中手动编写一个函数进行预测,但我想利用zoib模型是贝叶斯的事实,这样我就可以得到后验分布可能的响应值。我的计划是使用后验分布来计算每个预测周围的置信区间。

因为zoib在RJags内使用了MCMC方法,所以我研究了这两个解决方案:

    在RJags中
  1. manipulating the code

  2. 使用"NA" response variable

  3. 附加新数据

    第一个解决方案我不知道如何实现,因为zoib在内部运行RJag并且它运行的零膨胀模型非常复杂。我尝试了第二种解决方案,但它忽略了我追加的数据行" NA"响应值。

1 个答案:

答案 0 :(得分:0)

我通过电子邮件发送了zoib软件包开发人员,这就是响应。

  

目前,zoib函数只能输出应用了zoib回归的数据集中的X的Y的后验预测样本,但不能输出新的X集合。您的建议可以很容易地纳入新版本的软件包中,预计将在几周后推出。