我按组/级别运行数据回归。所以每个级别的注册表，即我的代码看起来非常相似：
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个级别，我必须做更多。非常感谢帮助。