我正在使用R中的佐伊布软件包来建立零膨胀的β回归模型。我正在寻找一种简单的方法来使用zoib生成的模型来计算新数据集的预测响应。所谓“新数据集”,我指的是不用于构建原始zoib模型的数据。
我知道我只需取zoib模型参数,然后用R手动编写一个函数来预测,但我想利用zoib模型是贝叶斯的事实,这样我就可以得到可能的响应值的后验分布。我的计划是使用后验分布来计算每个预测的置信区间。
由于zoib在RJags中使用MCMC方法,所以我研究了以下两种解决方案:
我不知道如何实现第一个解决方案,因为zoib在内部运行RJags,而它运行的零膨胀模型非常复杂。我尝试了第二个解决方案,但它忽略了我附加了"NA“响应值的数据行。
发布于 2018-02-03 00:05:02
我给zoib软件包开发人员发了电子邮件,这是他们的回应。
目前,zoib函数只能输出适用于zoib回归的数据集中的X,而不能输出Y的后验预测样本,但不能输出新的X集。您的建议可以很容易地纳入软件包的新版本,预计将在几周内发布。
https://stackoverflow.com/questions/48573063
复制相似问题