r - 将堆叠矩阵转换为分块矩阵

Question

我当前的数据是一系列相互堆叠的 3x3 矩阵。结构如下所示(以下示例是 3 矩阵，尽管我的实际数据更大/更长):

0   1   1
0   0   1
1   1   0
0   1   1
1   0   1
1   0   0
0   1   1
0   0   1
1   1   0

使用R，我希望它看起来像这样:

0   1   1                       
0   0   1                       
1   1   0                       
            0   1   1           
            1   0   1           
            1   0   0           
                        0   1   1
                        0   0   1
                        1   1   0

Answer 1

我们通过使用 gl 创建的分组索引将矩阵拆分为矩阵列表，并应用bdiag(来自Matrix)获取 block 对角稀疏矩阵

library(Matrix)
m2 <- bdiag(lapply(split(d1, as.integer(gl(nrow(d1), 3, nrow(d1)))), as.matrix))
m2
#9 x 9 sparse Matrix of class "dgCMatrix"
#                       
# [1,] . 1 1 . . . . . .
# [2,] . . 1 . . . . . .
# [3,] 1 1 . . . . . . .
# [4,] . . . . 1 1 . . .
# [5,] . . . 1 . 1 . . .
# [6,] . . . 1 . . . . .
# [7,] . . . . . . . 1 1
# [8,] . . . . . . . . 1
# [9,] . . . . . . 1 1 .

可以通过使用 as.matrix 包装将其转换为常规矩阵

as.matrix(m2)

数据

m1 <- structure(c(0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 
       1, 1, 1, 0, 1, 1, 0, 1, 1, 0), .Dim = c(9L, 3L))
d1 <- as.data.frame(m1)