Question

我有几个矩阵，让我们简单说我有3个矩阵。我想创建一个列表然后使用rbind将一个放在另一个上面。如果我手工完成，使用以下代码，它可以工作：

list<-list(matrix1,matrix2,matrix3)
test<-do.call("rbind",list)

我得到一个97947行4列的矩阵，这就是我想要的。但如果我做一个循环，它就不起作用了：

list2<-list()
for (i in 1:3)
{
y<-paste0("matrix",x)
list2[[x]] <- y
}
test2<-do.call("rbind",list2)

我得到一个3x1字符矩阵??? 有人可以指点我错误吗？任何评论将不胜感激。谢谢！！！！

Answer 1

考虑使用类似mget的函数来获取来自globalenvironment（默认环境）的所有矩阵对象，并将它们放在列表中。然后，您可以使用do.call方法并避免循环。这是一个玩具示例：

#  Some data
m1 <- matrix( 1:4 , 2 , byrow = TRUE )
m2 <- matrix( 1:4 , 2 , byrow = TRUE )
m3 <- matrix( 1:4 , 2 , byrow = TRUE )

#  Use mget to put them in a list. mget searches the .GlobalEnvironment (by default) for the object names in it's first argument
list <- mget( paste0( "m" , 1:3 ) )
list
#$m1
#     [,1] [,2]
#[1,]    1    2
#[2,]    3    4

#$m2
#     [,1] [,2]
#[1,]    1    2
#[2,]    3    4

#$m3
#     [,1] [,2]
#[1,]    1    2
#[2,]    3    4

#  rbind them
do.call( rbind , list )
#     [,1] [,2]
#[1,]    1    2
#[2,]    3    4
#[3,]    1    2
#[4,]    3    4
#[5,]    1    2
#[6,]    3    4

使用循环创建矩阵列表并将它们合并到R中

1 个答案: