我正在尝试summarise()通过任意组编写一个简单的包装到任意变量的包装,并且现在已经取得了进展我已经加载了正确的库版本,但我很困惑(再次)关于如何使用多个值取消引用参数.
我目前有以下功能......
table_summary <- function(df = .,
id = individual_id,
select = c(),
group = site,
...){
## Quote all arguments (see http://dplyr.tidyverse.org/articles/programming.html)
quo_id <- enquo(id)
quo_select <- enquo(select)
quo_group <- enquo(group)
## Subset the data
df <- df %>%
dplyr::select(!!quo_id, !!quo_select, !!quo_group) %>%
unique()
## gather() data, just in case there is > 1 variable selected to be summarised
df <- df %>%
gather(key = variable, value = value, !!quo_select)
## Summarise selected variables by specified groups
results <- df %>%
group_by(!!quo_group, variable) %>%
summarise(n = n(),
mean = mean(value, na.rm = TRUE))
return(results)
}
Run Code Online (Sandbox Code Playgroud)
如果我指定一个分组变量,那么大部分都可以使用并且有效...
> table_summary(df = mtcars, id = model, select = c(mpg), group = gear)
# A tibble: 3 x 4
# Groups: c(gear) [?]
gear variable n mean
<dbl> <chr> <int> <dbl>
1 3 mpg 15 16.10667
2 4 mpg 12 24.53333
3 5 mpg 5 21.38000
Run Code Online (Sandbox Code Playgroud)
...但是group_by(!!quo_group, variable)当我指定多个时失败了group = c(gear, hp)......
> mtcars$model <- rownames(mtcars)
> table_summary(df = mtcars, id = model, select = c(mpg), group = c(gear, hp))
Error in mutate_impl(.data, dots) :
Column `c(gear, hp)` must be length 32 (the group size) or one, not 64
Run Code Online (Sandbox Code Playgroud)
我回去和重新阅读编程dplyr文件,我看你能捕捉多个变量使用quos()代替enquo(),然后用所享有,它们接起来!!!,所以尝试...
table_summary <- function(df = .,
id = individual_id,
select = c(),
group = c(),
digits = 3,
...){
## Quote all arguments (see http://dplyr.tidyverse.org/articles/programming.html)
quo_id <- enquo(id)
quo_select <- enquo(select)
quo_group <- quos(group) ## Use quos() rather than enquo()
UQS(quo_group) %>% print() ## Check to see what quo_group holds
## Subset the data
df <- df %>%
dplyr::select(!!quo_id, !!quo_select, !!!quo_group)) %>%
unique()
## gather() data, just in case there is > 1 variable selected to be summarised
df <- df %>%
gather(key = variable, value = value, !!quo_select)
## Summarise selected variables by specified groups
results <- df %>%
group_by(!!!quo_group, variable) %>%
summarise(n = n(),
mean = mean(value, na.rm = TRUE))
return(results)
}
Run Code Online (Sandbox Code Playgroud)
...现在在第一次引用!!!quo_group``withindplyr :: select()regardless of how many variables are specified undergroup =`时失败了
> table_summary(df = mtcars, id = model, select = c(mpg), group = c(gear))
[[1]]
<quosure: frame>
~group
attr(,"class")
[1] "quosures"
Error in overscope_eval_next(overscope, expr) : object 'gear' not found
> traceback()
17: .Call(rlang_eval, f_rhs(quo), overscope)
16: overscope_eval_next(overscope, expr)
15: FUN(X[[i]], ...)
14: lapply(.x, .f, ...)
13: map(.x[matches], .f, ...)
12: map_if(ind_list, !is_helper, eval_tidy, data = names_list)
11: select_vars(names(.data), !(!(!quos(...))))
10: select.data.frame(., !(!quo_id), !(!quo_select), !(!(!quo_group)))
9: dplyr::select(., !(!quo_id), !(!quo_select), !(!(!quo_group)))
8: function_list[[i]](value)
7: freduce(value, `_function_list`)
6: `_fseq`(`_lhs`)
5: eval(quote(`_fseq`(`_lhs`)), env, env)
4: eval(quote(`_fseq`(`_lhs`)), env, env)
3: withVisible(eval(quote(`_fseq`(`_lhs`)), env, env))
2: df %>% dplyr::select(!(!quo_id), !(!quo_select), !(!(!quo_group))) %>%
unique()
1: table_summary(df = mtcars, id = model, select = c(mpg), group = c(gear))
Run Code Online (Sandbox Code Playgroud)
看起来很奇怪,我认为这个问题的根源在于!!!quo_group(即UQS(quo_group))打印出来~gear而不是列表,因为print()在工作示例中添加了一个显示...
> my_summarise <- function(df, ...) {
group_by <- quos(...)
UQS(group_by) %>% print()
df %>%
group_by(!!!group_by) %>%
summarise(a = mean(a))
}
> df <- tibble(
g1 = c(1, 1, 2, 2, 2),
g2 = c(1, 2, 1, 2, 1),
a = sample(5),
b = sample(5)
)
> my_summarise(df, g1, g2)
[[1]]
<quosure: global>
~g1
[[2]]
<quosure: global>
~g2
attr(,"class")
[1] "quosures"
# A tibble: 4 x 3
# Groups: g1 [?]
g1 g2 a
<dbl> <dbl> <dbl>
1 1 1 1.0
2 1 2 5.0
3 2 1 2.5
4 2 2 4.0
Run Code Online (Sandbox Code Playgroud)
我想明确地提供我希望分组的变量作为我的参数的参数,但是如果我将它们指定为是可行的...但是我决定在提供分组变量时测试我的函数是否有效...
table_summary <- function(df = .,
id = individual_id,
select = c(),
group = c(),
digits = 3,
...){
## Quote all arguments (see http://dplyr.tidyverse.org/articles/programming.html)
quo_id <- enquo(id)
quo_select <- enquo(select)
## quo_group <- quos(group)
quo_group <- quos(...)
UQS(quo_group) %>% print()
## Subset the data
df <- df %>%
dplyr::select(!!quo_id, !!quo_select, !!!quo_group) %>%
unique()
## gather() data, just in case there is > 1 variable selected to be summarised
df <- df %>%
gather(key = variable, value = value, !!quo_select)
## Summarise selected variables by specified groups
results <- df %>%
group_by(!!!quo_group, variable) %>%
summarise(n = n(),
mean = mean(value, na.rm = TRUE))
return(results)
}
Run Code Online (Sandbox Code Playgroud)
...但它没有,quos()再次取消引用拼接,NULL因此变量既不被选择也不被分组......
> table_summary(df = mtcars, id = model, select = c(mpg), gear, hp)
NULL
# A tibble: 1 x 3
variable n mean
<chr> <int> <dbl>
1 mpg 32 20.09062
> table_summary(df = mtcars, id = model, select = c(mpg), gear)
NULL
# A tibble: 1 x 3
variable n mean
<chr> <int> <dbl>
1 mpg 32 20.09062
Run Code Online (Sandbox Code Playgroud)
我已经通过这个周期,现在去了几次检查使用的每种方法enquo()和quos(),但不能看到我错了,尽管看了编程dplyr文档几次.
IIUC 你的帖子,你想提供c(col1, col2)给group_by(). 该动词不支持这一点:
group_by(mtcars, c(cyl, am))
#> Error in mutate_impl(.data, dots) :
#> Column `c(cyl, am)` must be length 32 (the number of rows) or one, not 64
Run Code Online (Sandbox Code Playgroud)
那是因为group_by()具有变异语义,而不是选择语义。这意味着您提供的表达式group_by()是转换表达式。这是一个令人惊讶但非常方便的功能。例如,您可以disp按如下方式划分为三个间隔进行分组:
group_by(mtcars, cut3 = cut(disp, 3))
Run Code Online (Sandbox Code Playgroud)
这也意味着,如果您提供c(cyl, am),它将把两列连接在一起并返回长度为 64 的向量,而它预期的长度为 32(行数)。
所以你的问题是你想要一个group_by()具有选择语义的包装器。使用 可以很容易地做到这一点dplyr::select_vars(),它将很快被提取到新的 tidyselect 包中:
library("dplyr")
group_wrapper <- function(df, groups = rlang::chr()) {
groups <- select_vars(tbl_vars(df), !! enquo(groups))
group_by(df, !!! rlang::syms(groups))
}
Run Code Online (Sandbox Code Playgroud)
group_by_at()或者,您可以包装具有选择语义的新动词:
group_wrapper <- function(df, groups = rlang::chr()) {
group_by_at(df, vars(!! enquo(groups)))
}
Run Code Online (Sandbox Code Playgroud)
让我们尝试一下:
group_wrapper(mtcars, c(disp, am))
#> # A tibble: 32 x 11
#> # Groups: disp, am [27]
#> mpg cyl disp hp drat wt qsec vs am gear carb
#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 21.0 6 160 110 3.90 2.62 16.5 0 1 4 4
#> # ... with 22 more rows
Run Code Online (Sandbox Code Playgroud)
该接口的优点是支持select()选择要分组的列的所有操作。
请注意,我使用rlang::chr()作为默认参数,因为选择函数不支持c()返回NULL(我们将来可能希望更改它)。chr()不带参数调用会返回长度为 0 的字符向量。