R 中多個 glm 對象的變量列表循環

Question

我無法終生明白為什么這在 R 中對我不起作用:(

## Unadjusted model 
noadj_var <- list("a",
                  "b",
                  "c",
                  "d",
                  )

for (i in 1:length(noadj_var)) {
  paste0("unadjusted_model_",i) <- glm(
    df$dependent_variable ~ i,
    data = df[df$present_2018 == "Yes" , ],
    family = binomial(link=logit),
    na.action = na.exclude
  )
}

我收到錯誤：

Error in model.frame.default(formula = df$dependent_variable ~ i, data = df[df$present_2018 ==  : 
 variable lengths differ (found for 'i')

我不明白為什么當我有 na.action = na.exclude 時會出現此錯誤

還有其他東西被搞砸了嗎？

在此之前，變量名稱“a”、“b”、“c”不是真實名稱，我不能只將它們更改為序列，它需要是一個字符串列表。

Answer 1

咕嚕聲purrr的另一種可能性：

 list("a",
       "b",
       "c",
       "d",
  ) %>% 
    purrr::set_names(paste0("unadjusted_model_", .)) %>% 
    purrr::map(~glm(reformulate(.x, 'dependent_variable'),data = df),
               family = binomial(link=logit),
               na.action = na.exclude ) %>% 
    purrr::map(broom::tidy)