0

给定一个矩阵 M 和具有不同可能值的较小矩阵，我试图列出将这些小矩阵的组合叠加到矩阵 M 中产生的所有可能矩阵。小矩阵将插入到 M 具有相同行/列的位置名字。

例如，说有：

 M <- matrix(rep(0, 49), nrow =7, ncol =7)
 rownames(M) <- colnames(M) <-seq(1,7)
 > M
 1 2 3 4 5 6 7
 1 0 0 0 0 0 0 0
 2 0 0 0 0 0 0 0
 3 0 0 0 0 0 0 0
 4 0 0 0 0 0 0 0
 5 0 0 0 0 0 0 0
 6 0 0 0 0 0 0 0
 7 0 0 0 0 0 0 0

# Generate first set of small matrices:
sub_mat_1_1 <- matrix(rep(1, 9), nrow =3, ncol =3)
rownames(sub_mat_1_1) <- colnames(sub_mat_1_1) <- c(2,3,5)
sub_mat_1_2 <- matrix(rep(2, 9), nrow =3, ncol =3)
rownames(sub_mat_1_2) <- colnames(sub_mat_1_2) <- c(2,3,5)
sub_mat_1_3 <- matrix(rep(3, 9), nrow =3, ncol =3)
rownames(sub_mat_1_3) <- colnames(sub_mat_1_3) <- c(2,3,5)
submatrix_1 <- list(sub_mat_1_1, sub_mat_1_2, sub_mat_1_3)

# Generate second set of small matrices:
submatrix_2 <- list()
sub_mat_2_1 <- matrix(rep(1, 4), nrow =2, ncol =2)
rownames(sub_mat_2_1) <- colnames(sub_mat_2_1) <- c(1,6)
sub_mat_2_2 <- matrix(rep(2, 4), nrow =2, ncol =2)
rownames(sub_mat_2_2) <- colnames(sub_mat_2_2) <- c(1,6)
submatrix_2 <- list(sub_mat_2_1, sub_mat_2_2)

# Generate list of small matrices:
submatrices <- list()
submatrices[[1]] <- submatrix_1
submatrices[[2]] <- submatrix_2

[[1]]
[[1]][[1]]
  2 3 5
2 1 1 1
3 1 1 1
5 1 1 1

[[1]][[2]]
  2 3 5
2 2 2 2
3 2 2 2
5 2 2 2

[[1]][[3]]
  2 3 5
2 3 3 3
3 3 3 3
5 3 3 3


[[2]]
[[2]][[1]]
  1 6
1 1 1
6 1 1

[[2]][[2]]
  1 6
1 2 2
6 2 2

由于第一个小矩阵集有 3 种可能性，第二个有 2 种可能性，我试图在不使用 for 循环的情况下输出所有 6 个可能的矩阵作为列表：

[[1]]
  1 2 3 4 5 6 7
1 1 0 0 0 0 1 0
2 0 1 1 0 1 0 0
3 0 1 1 0 1 0 0
4 0 0 0 0 0 0 0
5 0 1 1 0 1 0 0
6 1 0 0 0 0 1 0
7 0 0 0 0 0 0 0

[[2]]
  1 2 3 4 5 6 7
1 1 0 0 0 0 1 0
2 0 2 2 0 2 0 0
3 0 2 2 0 2 0 0
4 0 0 0 0 0 0 0
5 0 2 2 0 2 0 0
6 1 0 0 0 0 1 0
7 0 0 0 0 0 0 0

[[3]]
  1 2 3 4 5 6 7
1 1 0 0 0 0 1 0
2 0 3 3 0 3 0 0
3 0 3 3 0 3 0 0
4 0 0 0 0 0 0 0
5 0 3 3 0 3 0 0
6 1 0 0 0 0 1 0
7 0 0 0 0 0 0 0

[[4]]
  1 2 3 4 5 6 7
1 2 0 0 0 0 2 0
2 0 1 1 0 1 0 0
3 0 1 1 0 1 0 0
4 0 0 0 0 0 0 0
5 0 1 1 0 1 0 0
6 2 0 0 0 0 2 0
7 0 0 0 0 0 0 0

[[5]]
  1 2 3 4 5 6 7
1 2 0 0 0 0 2 0
2 0 2 2 0 2 0 0
3 0 2 2 0 2 0 0
4 0 0 0 0 0 0 0
5 0 2 2 0 2 0 0
6 2 0 0 0 0 2 0
7 0 0 0 0 0 0 0

[[6]]
  1 2 3 4 5 6 7
1 2 0 0 0 0 2 0
2 0 3 3 0 3 0 0
3 0 3 3 0 3 0 0
4 0 0 0 0 0 0 0
5 0 3 3 0 3 0 0
6 2 0 0 0 0 2 0
7 0 0 0 0 0 0 0

一般来说，我可能有 n 个给定的“小矩阵列表”，每个都有自己的矩阵数量。在这种情况下，我将如何使用应用类型函数？

4

1 回答 1

1

这给出了请求的输出，但我必须同意评论者的观点，即问题不是很好。您如何定义插入子矩阵的值的位置？我只是假设您希望将它们插入到特定的行和列子集中，如您请求的输出中所示...

将两个矩阵插入指定行/列的函数M（全局定义，可能是不好的做法）

tmpmatf <- function(m1,m2,rc1=c(2,3,5),rc2=c(1,6)) {
    pos1 <- as.matrix(expand.grid(rc1,rc1))
    pos2 <- as.matrix(expand.grid(rc2,rc2))
    M[pos1] <- m1
    M[pos2] <- m2
    M
}

现在使用expand.grid从每个子矩阵列表中创建一个包含所有索引组合的数据框，并使用alply(array-to-list) from在每个组合plyr上运行：tmpmatf

library(plyr)
alply(expand.grid(seq(length(submatrices[[1]])),
                  seq(length(submatrices[[2]]))),
      1,
      function(x) {
          tmpmatf(submatrices[[1]][[x[[1]]]],submatrices[[2]][[x[[2]]]])
      })

这应该适用于您的两个子矩阵列表中的任意数量的子矩阵，但是如果您确实有更多（> 2）个子矩阵列表，那么您没有给我们足够的信息来指定（例如）第三个列表子矩阵应该粘贴到更大的矩阵中......

请注意，此部分的第一部分（通过双列矩阵使用矩阵索引）比使用for循环将单个元素插入矩阵要快得多，但第二部分（alply）实际上并不比（任何？）快一对for迭代所有组合的嵌套循环——在这种情况下，后者可能更清晰/更容易调试......

于 2012-06-03T12:25:01.153 回答