为什么我的 dplyr 百分位数计算不适用于 tidy 评估？

Question

我对学生的测试数据很感兴趣，我希望使用 dplyr 将这些数据转换为百分位数。为了有一个最小的例子，想象一下三个学生的以下设置。

require(tidyverse)

tbl <- tibble(Name = c("Alice", "Bob", "Cat"), Test = c(16, 13, 15))

以下代码有效并产生所需的输出。

tbl %>% mutate(TestPercentile = cume_dist(Test) * 100)

# A tibble: 3 x 3
  Name   Test TestPercentile
  <chr> <dbl>          <dbl>
1 Alice    16          100  
2 Bob      13           33.3
3 Cat      15           66.7

但是，我实际上想以编程方式进行，因为有很多这样的列。

colname <- "Test"
percname <- str_c(colname, "Percentile")
tbl %>% mutate({{percname}} := cume_dist({{colname}}) * 100)

# A tibble: 3 x 3
  Name   Test TestPercentile
  <chr> <dbl>          <dbl>
1 Alice    16            100
2 Bob      13            100
3 Cat      15            100

为什么 cume_dist 当我尝试使用这样的 tidy 评估时，所有学生的百分位数都为 100？（理想情况下，如果允许我提出第二个问题，我该如何解决？）

Answer 1

如果以编程方式表示您想编写自己的函数，您可以这样做：

calculate_percentile <- function(data, colname) {

   data %>% 
    mutate("{{colname}}Percentile" := cume_dist({{colname}} * 100))

}

tbl %>% 
  calculate_percentile(Test)

 # A tibble: 3 x 3
  Name   Test TestPercentile
  <chr> <dbl>          <dbl>
1 Alice    16          1    
2 Bob      13          0.333
3 Cat      15          0.667

编辑多列 新数据

tbl <- tibble(Name = c("Alice", "Bob", "Cat"), Test = c(16, 13, 15), Test_math = c(16, 30, 55), Test_music = c(3, 78, 34))

calculate_percentile <- function(data, colnames) {

  data %>% 

    mutate(across({{colnames}}, ~cume_dist(.) * 100, .names = "{col}Percentile"))

}

test_columns <- c("Test_math", "Test_music")
tbl %>% 
  calculate_percentile(test_columns) 

# A tibble: 3 x 6
  Name   Test Test_math Test_music Test_mathPercentile Test_musicPercentile
  <chr> <dbl>     <dbl>      <dbl>               <dbl>                <dbl>
1 Alice    16        16          3                33.3                 33.3
2 Bob      13        30         78                66.7                100  
3 Cat      15        55         34               100                   66.7

为什么您的解决方案不起作用？因为您的解决方案按字面意义将 cume_dist 应用于字符串 "test":

tbl %>% mutate({{percname}} := print({{colname}}))

[1] "Test"
# A tibble: 3 x 5
  Name   Test Test_math Test_music TestPercentile
  <chr> <dbl>     <dbl>      <dbl> <chr>         
1 Alice    16        16          3 Test          
2 Bob      13        30         78 Test          
3 Cat      15        55         34 Test

为什么 TestPercentile 的值为 100？因为 "test" 的 cume_dist 是 1:

cume_dist("test")
#[1] 1

所以我们需要 R 告诉它不要对字符串 "test" 本身求值，而是要查找具有该名称的变量，我们可以这样做：

tbl %>% mutate({{percname}} := cume_dist(!!parse_quo(colname, env = global_env())) * 100)

# A tibble: 3 x 5
  Name   Test Test_math Test_music TestPercentile
  <chr> <dbl>     <dbl>      <dbl>          <dbl>
1 Alice    16        16          3          100  
2 Bob      13        30         78           33.3
3 Cat      15        55         34           66.7

#Check that this uses the values of "Test" and not "Test" per se:
tbl %>% mutate({{percname}} := print(!!parse_quo(colname, env = global_env())))

[1] 16 13 15
# A tibble: 3 x 5
  Name   Test Test_math Test_music TestPercentile
  <chr> <dbl>     <dbl>      <dbl>          <dbl>
1 Alice    16        16          3             16
2 Bob      13        30         78             13
3 Cat      15        55         34             15

Answer 2

将列名作为字符串传递：

library(dplyr)
library(rlang)

return_percentile <- function(data, colname) {
   percname <- paste0(colname, "Percentile")
   data %>% mutate({{percname}} := cume_dist(!!sym(colname)) * 100)
}

tbl %>% return_percentile("Test")

# A tibble: 3 x 3
#  Name   Test TestPercentile
#  <chr> <dbl>          <dbl>
#1 Alice    16          100  
#2 Bob      13           33.3
#3 Cat      15           66.7

传递不带引号的列名：

return_percentile <- function(data, colname) {
  percname <- paste0(deparse(substitute(colname)), "Percentile")
  data %>% mutate({{percname}} := cume_dist({{colname}}) * 100)
}

tbl %>% return_percentile(Test)

# A tibble: 3 x 3
#  Name   Test TestPercentile
#  <chr> <dbl>          <dbl>
#1 Alice    16          100  
#2 Bob      13           33.3
#3 Cat      15           66.7

为什么我的 dplyr 百分位数计算不适用于 tidy 评估？

Why does my dplyr percentile calculation not work with tidy evaluation?

r

dplyr

tidyeval