我正在使用dplyr中生成的模型列表
library(dplyr)
library(magrittr)
mod.list <- iris %>%
group_by(Species) %>%
do(mod = lm(Petal.Length ~ Petal.Width, data = .))
如果我绘制每个模型，则所有工程都能按预期进行
par(mfrow = c(2,2))
mod.list %>%
do(.$mod %>% plot)
但是，如果引入过滤器，则会出现错误
mod.list %>%
filter(row_number() <= 1) %>%
do(.$mod %>% plot)
Error in xy.coords(x, y, xlabel, ylabel, log) :
'x' is a list, but does not have components 'x' and 'y'
我也尝试了另一种类型的过滤器，但错误是相同的
mod.list %>%
filter(Species == "setosa") %>%
do(.$mod %>% plot)
有人知道为什么会这样吗？
要进行调试，您可以使用
mod.list %>%
filter(row_number() <= 1) %>%
do(.$mod %>% (function(x) browser()))
然后，您会看到这class(x)是一个list。您需要第一个（唯一的）元素，所以
mod.list %>%
filter(row_number() <= 1) %>%
do(.$mod %>% `[[`(i=1) %>% plot)
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句