在using predict with a list of lm() objects中,@ JD Long要求解决我经常遇到的问题;通常,当模型适合并且要生产图表时,它会出现。在@Joshua Ulrich编辑的JD Long的解决方案很优雅,除了paste
所需要的。
ddply
帝国的首领@hadley发表了mdply
的解决方案,我发现这个解决方案难以理解,远非优雅。
现在dplyr
(可能还有tidyr
)正在接管:是否有更直观的方式来处理这项常见工作?
有关数据和示例代码,请参阅上面引用的JD Long示例。
答案 0 :(得分:1)
使用引用的question中的数据,您可以对dplyr使用以下方法。
library(dplyr)
newData %>%
group_by(state) %>%
mutate(pred = predict(modelList[[as.character(state[1])]],
as.data.frame(year)))
但是,我不确定这种解决方案是否比现有解决方案更优雅。