我正在使用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)案例