根据标题,有没有一种方法可以将一个数据切片(一个n乘n矩阵)附加到R中现有的n维矩阵上?
例如,我有以下内容:
one <- array(1, dim = c(3, 3))
two <- array(2, dim = c(3, 3))
three <- array(3, dim = c(6, 6))
然后我想把它转换成一个6x6x3矩阵,我可以使用它,它看起来如下:
[[1]]
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 1 1 NA NA NA
[2,] 1 1 1 NA NA NA
[3,] 1 1 1 NA NA NA
[4,] NA NA NA NA NA NA
[5,] NA NA NA NA NA NA
[6,] NA NA NA NA NA NA
[[2]]
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 2 2 2 NA NA NA
[2,] 2 2 2 NA NA NA
[3,] 2 2 2 NA NA NA
[4,] NA NA NA NA NA NA
[5,] NA NA NA NA NA NA
[6,] NA NA NA NA NA NA
[[3]]
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 3 3 3 3 3 3
[2,] 3 3 3 3 3 3
[3,] 3 3 3 3 3 3
[4,] 3 3 3 3 3 3
[5,] 3 3 3 3 3 3
[6,] 3 3 3 3 3 3
我知道如何通过自己的代码做到这一点,所以我更感兴趣的是是否有现有的库函数支持这一点。
除了abind
,我想您还需要计算出矩阵的最大大小,并创建用适当数量的NA填充的矩阵?
padmat <- function(X,m,n) {
Y <- matrix(NA,m,n)
Y[1:nrow(X),1:ncol(X)] <- X
Y
}
one <- array(1, dim = c(3, 3))
two <- array(2, dim = c(3, 3))
three <- array(3, dim = c(6, 6))
mlist <- list(one,two,three)
maxrows <- max(sapply(mlist,nrow))
maxcols <- max(sapply(mlist,ncol))
mlist2 <- lapply(mlist,padmat,m=maxrows,n=maxcols)
library(abind)
abind(mlist2,along=3)