Question

我需要将dplyr / broom代码的一部分（有效的）转换为一个函数，因为我会多次（多次）调用它。

我陷入困境-这很可能与将非标准评估与标准评估混为一谈。

这里我直接从the vignette 'broom and dplyr'中获取代码

library(tidyverse)
library(broom)

data(Orange)

Orange %>% 
  nest(-Tree) %>% 
  mutate(
    test = map(data, ~ cor.test(.x$age, .x$circumference)), 
    tidied = map(test, tidy)
  ) %>% 
  unnest(tidied, .drop = TRUE)

这有效：

    Tree  estimate statistic      p.value parameter  conf.low conf.high                               method alternative
1    1 0.9854675  12.97258 4.851902e-05         5 0.9012111 0.9979400 Pearson's product-moment correlation   two.sided
2    2 0.9873624  13.93129 3.425041e-05         5 0.9136142 0.9982101 Pearson's product-moment correlation   two.sided
3    3 0.9881766  14.41188 2.901046e-05         5 0.9189858 0.9983260 Pearson's product-moment correlation   two.sided
4    4 0.9844610  12.53575 5.733090e-05         5 0.8946782 0.9977964 Pearson's product-moment correlation   two.sided
5    5 0.9877376  14.14686 3.177093e-05         5 0.9160865 0.9982635 Pearson's product-moment correlation   two.sided

现在的重点是，我想利用它做一个功能。

所以，如果我尝试这样做：

afunction <- function(data, var) {
  data %>% 
    nest(-Tree) %>% 
    mutate(
      test = map(data, ~ cor.test(.x$age, .x$var)), # S3 list-col
      tidied = map(test, tidy)
    ) %>% 
    unnest(tidied, .drop = TRUE)
}

它失败了。

 Error in cor.test.default(.x$age, .x$var) : 'x' and 'y' must have the same length

我尝试使用NSE，引号，半引号。我承认我随机尝试了一下，因为我找不到如何让NSE和SE与$运算符一起很好地发挥作用的适当教程。

任何解决方案-尤其是一种可以扩展并教会我如何一劳永逸地解决这些问题的解决方案？我也很高兴获得相关书籍/教程的指导。

功能编程：使用broom nest-> tidy-> unnest并在函数内映射

0 个答案: