r - 新样本均值的 R 预测区间

Question

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

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

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

Answer 1

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

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

然后可能只是：

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

我说“可能”是因为你没有提供数据集来测试这个，我认为提供这样的数据是你的责任。我不知道这是否复制了 JMP 方法或 WinBUGS 方法。