5

我想从矩阵列表中创建交错矩阵。

示例输入：

> l <- list(a=matrix(1:4,2),b=matrix(5:8,2))
> l
$a
     [,1] [,2]
[1,]    1    3
[2,]    2    4

$b
     [,1] [,2]
[1,]    5    7
[2,]    6    8

预期输出：

1    3
5    7
2    4
6    8

我已经检查了 gdata 中的交错函数，但它没有显示列表的这种行为。任何帮助表示赞赏。

4

2 回答 2

5

这是一个单行：

do.call(rbind, l)[order(sequence(sapply(l, nrow))), ]
#      [,1] [,2]
# [1,]    1    3
# [2,]    5    7
# [3,]    2    4
# [4,]    6    8

为了帮助理解，矩阵首先用 相互堆叠do.call(rbind, l)，然后以正确的顺序提取行：

sequence(sapply(l, nrow))
# a1 a2 b1 b2 
#  1  2  1  2 

order(sequence(sapply(l, nrow)))
# [1] 1 3 2 4

它适用于任意数量的矩阵，即使它们的行数不同，它也会做“正确的事情”（主观的）。

于 2013-11-05T12:02:08.007 回答
4

您无需重新发明轮子，只需对其进行修改即可到达目的地。

interleave“gdata”中的函数首先...允许您指定data.frame要放在一起的多个 s 或矩阵。该函数的前几行如下所示：

head(interleave)
# 
# 1 function (..., append.source = TRUE, sep = ": ", drop = FALSE) 
# 2 {                                                              
# 3     sources <- list(...)                                       
# 4     sources[sapply(sources, is.null)] <- NULL                  
# 5     sources <- lapply(sources, function(x) if (is.matrix(x) || 
# 6         is.data.frame(x))

你可以像我在这个 Gist中所做的那样重写第 1 行和第 3 行来创建一个list版本interleave（在这里，我称之为Interleave）

head(Interleave)
#                                                                     
# 1 function (myList, append.source = TRUE, sep = ": ", drop = FALSE) 
# 2 {                                                                 
# 3     sources <- myList                                             
# 4     sources[sapply(sources, is.null)] <- NULL                     
# 5     sources <- lapply(sources, function(x) if (is.matrix(x) ||    
# 6         is.data.frame(x))

它有效吗？

l <- list(a=matrix(1:4,2),b=matrix(5:8,2), c=matrix(9:12,2))
Interleave(l)
#   [,1] [,2]
# a    1    3
# b    5    7
# c    9   11
# a    2    4
# b    6    8
# c   10   12
于 2013-11-05T04:30:40.353 回答