我正在按组/级别对数据进行回归。所以每个级别都有一个 reg，即我的代码看起来与此非常相似：
library('dplyr')
regressions <- mtcars %>%
group_by(cyl) %>%
do(mod = lm(mpg ~ wt, .))
如何在具有相同组的较小新数据集上提取预测？即水平相同但样本较小。我想一次对所有级别的新数据集进行预测。我知道扫帚中的“增强”可以为您提供预测。但我不知道如何同时为所有级别做到这一点。
我使用的代码如下所示。
library('broom')
aa <- as.data.frame(augment(regressions, newdata=test))
我也试过
aa <- as.data.frame(augment(regressions, mod, newdata=test))
它给出了拟合值，但复制了我的新数据以匹配原始数据，但每个观察的预测都不同。我的数据有 44 个级别，我必须做更多。非常感谢您的帮助。