我正在使用 R 中的zoib包来构建零膨胀的 beta 回归模型。我正在寻找一种简单的方法来使用 zoib 生成的模型来计算新数据集的预测响应。“新数据集”是指未用于构建原始 zoib 模型的数据。
我知道我可以只取zoib模型参数并在 R 中手动编写一个函数来预测，但我想利用 zoib 模型是贝叶斯的事实，所以我可以获得可能响应值的后验分布。我的计划是使用后验分布来计算每个预测的置信区间。
因为 zoib 在RJags中使用 MCMC 方法，所以我研究了这两种解决方案：
在 RJags 中操作代码
第一个解决方案我不知道如何实现，因为 zoib 在内部运行 RJags，而且它运行的零膨胀模型非常复杂。我尝试了第二种解决方案，但它只是忽略了我附加“NA”响应值的数据行。