Question

鉴于从一个数据集创建的回归模型，我一直在使用WinBUGS来构建围绕第二个数据集的平均值的预测间隔（PI）。我刚刚在R中发现了“预测”功能，但它在第二个数据集中围绕每个预测值提供了PI。我在这里和网上搜索了R帮助，只找到了单独成员的间隔。

这些区间的平均值显然与预测样本均值周围的PI不同（我已根据WinBUGS得到的值进行了测试）。

如何让R给我平均值的PI？

Answer 1

曾经有一个R mean.data.frame函数，但它已被弃用，然后被删除。您可以通过以下方式获得相同的结果：

  mean.vec <-  lapply(na.omit(dfrm), mean)

然后可能只是：

 predict(fit, newdata=data.frame(mean.vec) )

我说＆＃39;可能＆＃39;因为您没有提供任何数据集来测试这一点，因此我认为这是您的责任。我不知道这是复制JMP方法还是WinBUGS方法。