具有尊重顺序的向量的可能组合

Question

x <- letters[1:4]
x
# [1] "a" "b" "c" "d"

t(combn(x, 2))
#   [,1] [,2]
# [1,] "a"  "b" 
# [2,] "a"  "c" 
# [3,] "a"  "d" 
# [4,] "b"  "c" 
# [5,] "b"  "d" 
# [6,] "c"  "d"

如果我还有b-a,c-a...d-c的逆组合，我应该怎么写代码呢。一共12种组合。

Answer 1

您可以使用基数 R 中的 expand.grid 来获取向量的所有可能组合（即 16 种组合），然后使用 subset（或 [.data.frame）以便两列中的值永远不会连续相等（导致预期的 12 种组合）：

x <- letters[1:4]
subset(expand.grid(rep(list(x),2)), Var1 != Var2)
#   Var1 Var2
#2     b    a
#3     c    a
#4     d    a
#5     a    b
#7     c    b
#8     d    b
#9     a    c
#10    b    c
#12    d    c
#13    a    d
#14    b    d
#15    c    d

data.table 的 cross-join (CJ) 函数的替代方法：

libraray(data.table)
CJ(x, x)[V1 != V2]
#    V1 V2
# 1:  a  b
# 2:  a  c
# 3:  a  d
# 4:  b  a
# 5:  b  c
# 6:  b  d
# 7:  c  a
# 8:  c  b
# 9:  c  d
#10:  d  a
#11:  d  b
#12:  d  c

Answer 2

程序包 gtools 具有 permutations 功能。

x <- letters[1:4]
library(gtools)
permutations(length(x), 2, x, repeats = FALSE)
#       [,1] [,2]
#  [1,] "a"  "b" 
#  [2,] "a"  "c" 
#  [3,] "a"  "d" 
#  [4,] "b"  "a" 
#  [5,] "b"  "c" 
#  [6,] "b"  "d" 
#  [7,] "c"  "a" 
#  [8,] "c"  "b" 
#  [9,] "c"  "d" 
# [10,] "d"  "a" 
# [11,] "d"  "b" 
# [12,] "d"  "c"

Answer 3

另一个选项dplyr/tidyr

library(dplyr)
library(tidyr)
data_frame(x1=x,x2=x) %>% 
        expand(x1, x2) %>%
        filter(x1!=x2)
#      x1    x2
#    (chr) (chr)
#1      a     b
#2      a     c
#3      a     d
#4      b     a
#5      b     c
#6      b     d
#7      c     a
#8      c     b
#9      c     d
#10     d     a
#11     d     b
#12     d     c

具有尊重顺序的向量的可能组合

Possible combinations of a vector with respect order

r

combn