按矩陣中的行獲取所有可能的組合

Question

我有一個m x n矩陣，如下所示：

1 2 3
4 5 6

按行獲取所有可能組合的最快方法是什么？ 在這種情況下，那將是c(1,4), c(1,5), c(1,6), c(2,4), c(2,5) ... c(3,5), c(3,6)

如何使用矢量化方法解決這個問題？ 通常， m × n矩陣將具有n^m這樣的組合。

Answer 1

您可以使用expand.grid函數來獲取每一行中的元素的所有組合來說，使用行的列表split如圖所示這里並通過該列表中的每個元素expand.grid與do.call功能：

(m <- rbind(1:3, 4:6))
#      [,1] [,2] [,3]
# [1,]    1    2    3
# [2,]    4    5    6
do.call(expand.grid, split(m, rep(1:nrow(m), ncol(m))))
#   1 2
# 1 1 4
# 2 2 4
# 3 3 4
# 4 1 5
# 5 2 5
# 6 3 5
# 7 1 6
# 8 2 6
# 9 3 6

這是一個使用3 x 2矩陣而不是2 x 3矩陣的示例：

(m <- matrix(1:6, nrow=3))
#      [,1] [,2]
# [1,]    1    4
# [2,]    2    5
# [3,]    3    6
do.call(expand.grid, split(m, rep(1:nrow(m), ncol(m))))
#   1 2 3
# 1 1 2 3
# 2 4 2 3
# 3 1 5 3
# 4 4 5 3
# 5 1 2 6
# 6 4 2 6
# 7 1 5 6
# 8 4 5 6