2

我需要制作一个矩阵R,其元素来自我之前定义的矩阵。例如,我有 4 个矩阵,

w <- matrix(c(1,2,3,4),2,2)
x <- matrix(c(5,6,7,8),2,2)
y <- matrix(c(9,10,11,12),2,2)
z <- matrix(c(13,14,15,16),2,2)

那么,新矩阵应该是一个4X4矩阵,其中w是一个[1:2,1:2]元素,x一个[1:2,3:4]元素,y一个[3:4,1:2]元素,z一个[3:4,3:4]元素。

我怎样才能快速做到这一点?

4

1 回答 1

4

我们可以创建一个array然后循环遍历第三维,然后rbind它。

ar1 <- array(c(w, x, y,z), dim=c(2, 4,2))
do.call(rbind,lapply(seq(dim(ar1)[3]), function(i) ar1[,,i]))
#     [,1] [,2] [,3] [,4]
#[1,]    1    3    5    7
#[2,]    2    4    6    8
#[3,]    9   11   13   15
#[4,]   10   12   14   16

或者正如评论中提到的@thelatemail

apply(array(c(w,x,y,z), dim=c(2,4,2)), 2, I)

whereI代表inhibit interpretation或使用identity代替I

于 2016-05-05T06:17:47.653 回答