R语言 Tidy eval在函数中失败

disbfnqx 于 2023-10-13 发布在其他

关注(0)|答案(1)|浏览(102)

我试图自动化.tex表的生成，我想在stargazer的add.lines部分添加一个名为mean_outcome的值，但我在生成这个值时遇到了麻烦，我想我处理tidy eval时出错了。我已经试过了！”““但是什么都不起作用。

df <- tibble::tribble(~hhid, ~treatment_arm, ~round_pooled, ~tot_lstock_count, ~tot_animal_tlu,
                      "1020201023",1,1,10,1,
                      "1020201023",1,11,6,0.600000023841858,
                      "1020201036",1,1,10,0.100000001490116,
                      "1020201036",1,11,4,0.400000005960464,
                      "1020201039",1,1,9,0.0900000035762787,
                      "1020201057",1,1,12,0.480000019073486,
                      "1020201068",1,1,2,0.200000002980232,
                      "1020201095",1,1,6,0.330000013113022,
                      "1020201095",1,11,7,0.400000005960464,
                      "1020201103",1,1,2,0.200000002980232,
                      "1020201116",1,1,7,0.0599999986588955,
                      "1020201116",1,11,2,0.0199999995529652,
                      "1020201139",1,11,12,0.480000019073486,
                      "1020201144",1,1,5,0.5,
                      "1020201144",1,11,8,0.800000011920929,
                      "1020201146",1,1,6,0.0599999986588955,
                      "1020201146",1,11,6,0.0599999986588955,
                      "1020201159",1,11,15,1.5,
                      "1020201180",1,11,11,0.560000002384186,
                      "1020201205",1,1,2,0.0199999995529652)

reg_coef_tables <- function(df, y) {
  
  outcome = eval(`i`)  
  
  mean_outcome  <- df %>% 
    
    filter(round_pooled == 11 & treatment_arm == 1) %>% 
    
    summarise(mean = mean(!!outcome, na.rm = T)) %>% 
    
    pull(mean)
  
  reg1 <- felm(outcome ~ treatment_arm, data = df)

stargazer(reg1,
          align = TRUE, 
          dep.var.labels = c("(1)"),
          omit.stat = c("f", "adj.rsq", "ser"),
          no.space = TRUE,
          digits = 3, # number of decimals
          add.lines = list(c('Mean outcome', rep(`mean_outcome`, 1)),
                           dep.var.caption = "", # remove dep var header
                           out = paste0(outcome, ".tex"))
}          
          
liv_vars  = c("tot_lstock_count", "tot_animal_tlu")
liv_out_list = vector("list", 2)
for(i in liv_vars) {
            liv_out_list[[i]] <- reg_coef_tables(df, i)
}

r

来源：https://stackoverflow.com/questions/77222622/tidy-eval-failing-within-a-function

1条答案

按热度按时间

vuktfyat1#

撇开stargazer（因为，恕我直言，它很可怕）和felm（因为我不知道它来自哪个包），因为这两个都与tidyverse的问题无关，试试这个

reg_coef_tables <- function(df, y) {
  outcome = as.name(y) 
  mean_outcome  <- df %>% 
    filter(round_pooled == 11 & treatment_arm == 1) %>% 
    summarise(mean = mean(!!outcome, na.rm = T)) %>% 
    pull(mean)
}
liv_vars  = c("tot_lstock_count", "tot_animal_tlu")
liv_out_list = list()
for(i in liv_vars) {
  liv_out_list[[i]] <- reg_coef_tables(df, i)
}  
liv_out_list
$tot_lstock_count
[1] 7.888889

$tot_animal_tlu
[1] 0.5355556

或者，更简洁地说：

mean_outcome <- df %>% 
  filter(round_pooled == 11 & treatment_arm == 1) %>% 
  summarise(
    across(
      c(tot_lstock_count, tot_animal_tlu), 
      \(x) mean(x, na.rm = TRUE)
    )
  )
mean_outcome
# A tibble: 1 × 2
  tot_lstock_count tot_animal_tlu
             <dbl>          <dbl>
1             7.89          0.536

你可以在NSE here上找到一个有用的引物。

赞(0）回复(0）举报 2023-10-13

我来回答

R语言 Tidy eval在函数中失败

1条答案

相关问题

热门标签

最新问答