这是使用dplyr::bind_cols 和purrr::map2 的一种方法,它在大量列时似乎比base 快得多。我对分析的了解不足,无法猜测原因,因为它比其他选项感觉更复杂。我不确定使用mutate_ 动词是否容易做到这一点,尽管可以更正。
编辑:添加了一个带有dplyr::do 的选项,这是进行计算的“预期”方式,它不能完全适合mutate 函数。 mutate 的问题在于它希望只创建一列。我认为除了使用 map 来构造单独的 mutate 调用(我无法想象它会更快)之外,这是最好的选择。
library(tidyverse)
set.seed(4321)
df <- matrix(rnorm(1000000), ncol = 1000) %>%
as_tibble()
microbenchmark::microbenchmark(
base = df[, 1001:1500] <- df[, 1:500] - df[, 501:1000],
base2 = df %>% magrittr::inset(, 1001:1500, .[, 1:500] - .[, 501:1000]),
map = df %>% bind_cols(map2(.x = .[, 1:500], .y = .[, 501:1000], .f = ~.x - .y)),
nomap = df %>% bind_cols(.[, 1:500] - .[, 501:1000]),
do = df %>% do(.[, 1:500] - .[, 501:1000])
)
#> Unit: milliseconds
#> expr min lq mean median uq max neval
#> base 32.928171 36.394238 39.362308 37.361149 39.454822 112.76356 100
#> base2 33.302556 35.500491 38.888530 37.433863 40.207799 84.08674 100
#> map 4.693637 5.139985 5.967655 5.468398 6.264793 12.20658 100
#> nomap 23.061348 25.016053 28.598282 26.973913 29.574478 79.97451 100
#> do 21.906042 23.460822 27.049262 25.135640 26.596373 80.01928 100
#> cld
#> c
#> c
#> a
#> b
#> b
由reprex package (v0.2.0) 于 2018 年 5 月 11 日创建。