【发布时间】:2018-02-01 22:53:19
【问题描述】:
我正在使用 R 中的 zoib 包来构建零膨胀的 beta 回归模型。我正在寻找一种简单的方法来使用 zoib 生成的模型来计算新数据集的预测响应。 “新数据集”是指未用于构建原始 zoib 模型的数据。
我知道我可以只取 zoib 模型参数并在 R 中手动编写一个函数来预测,但我想利用 zoib 模型是贝叶斯的事实,所以我可以获得后验分布可能的响应值。我的计划是使用后验分布来计算每个预测的置信区间。
因为 zoib 在RJags 中使用了 MCMC 方法,所以我研究了这两种解决方案:
manipulating the code 在 RJags 中
用"NA" response variable 附加新数据
第一个解决方案我不知道如何实现,因为 zoib 在内部运行 RJags 并且它运行的零膨胀模型非常复杂。我尝试了第二种解决方案,但它只是忽略了我附加“NA”响应值的数据行。
【问题讨论】: