在嵌套数据框架中使用扫描：：整齐将置信区间添加到线性模型中

Question

我一直在尝试按照Hadley提出的方法在嵌套数据框架中按照https://r4ds.had.co.nz/many-models.html运行多个模型

我在下面编写了下面的代码来创建多个线性模型：

#create nested data frame
by_clin_grp <- df_long %>%
  group_by(clin_risk) %>%
  nest()

#create a function to run the linear models
model <- function(df){
  lm(outcome ~ age + sex + clin_sub_grp, data =df)
}

#run the models
by_clin_grp <- by_clin_grp %>%
  mutate(model = map(data,model))


#unnest the models using the broom package 
by_clin_grp_unnest <- by_clin_grp %>%
  mutate(tidy = map(model, broom::tidy))%>%
  unnest(tidy)

但是，我需要在回归估计的前后增加置信区间。似乎我应该能够按照broom::tidy按照帮助页添加它们，但是我无法在上面的代码中正确地指定这一点？

Answer 1

您必须在map中指定相关参数。有两种可能性：

by_clin_grp_unnest <- by_clin_grp %>%
  mutate(tidy = map(model, ~broom::tidy(.x, conf.int = TRUE, conf.level = 0.99)))%>%
  unnest(tidy)

或

by_clin_grp_unnest <- by_clin_grp %>%
  mutate(tidy = map(model, broom::tidy, conf.int = TRUE, conf.level = 0.99)) %>%
  unnest(tidy)

两者都是等价的，只是语法不同。