如何根据 2 个不同矩阵的元素对矩阵进行排序

Question

我有两个矩阵。一个（per）是一个排列，另一个矩阵（wh）基于每个矩阵的每个元素包含一些计算。

    > per=permutations(3,3,c(1,2,3))
    > per
         [,1] [,2] [,3]
    [1,]    1    2    3
    [2,]    1    3    2
    [3,]    2    1    3
    [4,]    2    3    1
    [5,]    3    1    2
    [6,]    3    2    1

> wh
     [,1] [,2] [,3]
[1,] "w1" "w2" "w3"
[2,] "w1" "w3" "w2"
[3,] "w2" "w1" "w3"
[4,] "w2" "w3" "w1"
[5,] "w3" "w1" "w2"
[6,] "w3" "w2" "w1"

我想要做的是创建一个新矩阵 (whr)，它将以 w1,w2,w3 的形式对矩阵 wh 进行排序

我现在正在做的是尝试用 for 循环填充矩阵。对于第 i 行，我检测每个矩阵中的元素所在的位置（例如 4）并给定这个位置，我检测它在相应 w4 的位置，然后我相应地填充矩阵 whr。

我用循环来做这个

for (i in 1:720){
    whr[i,per[i,1]]=wh[i,1]
    whr[i,per[i,2]]=wh[i,2]
    whr[i,per[i,3]]=wh[i,3]
    whr[i,per[i,4]]=wh[i,4]
    whr[i,per[i,5]]=wh[i,5]
    whr[i,per[i,6]]=wh[i,6]

它可以完成这项工作，但速度非常慢。我尝试对其进行矢量化，但我不确定如何。

我给出一般命令

whr[,per[,1]]=wh[,1]

对所有元素 1 和相应的 2-6 执行相同的操作，但它不起作用。 任何建议都更受欢迎。

Answer 1

你给我们留下了一些空白来尝试重现你的代码（没有提到包以及混淆了 3 或 6），但假设你想要类似的东西

library(gtools)
n   = 3
per = permutations(n, n, 1:n)
wh  = permutations(n, n, paste0("wh", 1:n))
whr = matrix(NA, nrow=factorial(n), ncol=n)
for (i in 1:factorial(n)){ for (j in 1:n){ whr[i,per[i,j]] = wh[i,j] } }

如果您不希望在最后一行出现双循环，则将其替换为以下内容将减少为单循环并产生相同的结果

for (i in 1:factorial(n)){ whr[i,] <- wh[i, order(per[i,])] }