我正在使用dplyr
+ broom
包为我的
my_data %>% group_by(my_groups) %>%
do(tidy(lm( attended_school ~ treated + I_post_treatmet, data=.)
))
(实际数据限制访问,不能在此分享)
我收到以下错误:
Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
0 (non-NA) cases
如果我理解正确,这意味着无法为某个子组运行回归,可能是由于某些变量中的所有NA。
我如何捕获这些错误,以便最终收到一个数据集,其中包含回归运行成功的子组的回归输出?
编辑:尝试了类似于评论中建议的链接的语法并遇到了同样的问题:
my_data %>% group_by(my_groups) %>%
do(aa=tryCatch(lm( attended_school ~ treated + I_post_treatmet, data=.)
))
lm.fit中的错误(x,y,offset = offset,singular.ok = singular.ok,...): 0(非NA)案件