Purrr：哪个嵌套数据框遇到了错误？

Question

这是一个嵌套数据。

df1 <- tibble::tribble(
   ~A,        ~B,  ~group,
   4L,         1L,  "A",
   7L,         4L,  "A",
  NA_integer_, 1L,  "B",
  NA_integer_, 10L, "B")

df2 <- df1 %>% group_by(group) %>% nest()

我需要使用 purrr:map 运行 lm。

map(df2$data, ~lm(A~B, data=.x))

找出哪些嵌套数据遇到错误的最佳方法是什么，即我如何知道 B 组有问题。

你能建议一个可能或安全地使用 purrr 的解决方案吗？

Answer 1

你可以possibly()在这里使用。NA在我的示例中，如果模型有错误，我会返回它。

首先，我 make posslm， usingotherwise告诉它如果使用后出现错误则返回什么lm。

posslm = possibly(lm, otherwise = NA)

map然后您可以使用in创建一个新的模型列mutate。完成后，filter转到NA新变量的行，然后拉出group.

mutate(df2, mod = map(data, ~posslm(A~B, data=.x))) %>%
    filter( is.na(mod) ) %>%
    pull(group)

[1] "B"