R:使用索引矩阵从数组中提取矩阵
问题描述:
我在R中编码,并且我有一个包含数据的3维数组(在本例中为ab)。然后我有一个包含第三个数组维度(idx)索引的矩阵。该矩阵具有相同数量的行和列数组。我想使用idx中包含的索引从数组中提取数据,以获得具有相同维度idx的矩阵。请看下面的例子:R:使用索引矩阵从数组中提取矩阵
a <- c(1:9)
b <- rev(a)
#array of data
ab <- array(c(a,b), dim = c(3,3,2))
ab
, , 1
[,1] [,2] [,3]
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
, , 2
[,1] [,2] [,3]
[1,] 9 6 3
[2,] 8 5 2
[3,] 7 4 1
#matrix of indices
idx <- matrix(sample(1:2,9,replace=TRUE), nrow = 3)
idx
[,1] [,2] [,3]
[1,] 2 2 2
[2,] 2 1 1
[3,] 1 1 1
#now I want to get the following matrix:
[,1] [,2] [,3]
[1,] 9 6 3
[2,] 8 5 8
[3,] 3 6 9
#these two don´t do the job
ab[idx]
ab[ , ,idx]
有没有人知道我怎么能得到它?
非常感谢!
萨拉
答
我们需要用于行/列和第三维度中的索引(从“IDX”),以提取的元素。我们通过cbind
这个行索引,列索引和'idx'来做到这一点。
i1 <- dim(ab)[1]
j1 <- dim(ab)[2]
matrix(ab[cbind(rep(seq_len(i1), j1),rep(seq_len(j1), each = i1), c(idx))], ncol=3)
# [,1] [,2] [,3]
#[1,] 9 6 3
#[2,] 8 5 8
#[3,] 3 6 9
答
可能有点复杂,但它的作品至少: 首先我所有的IDX的:
[,1] [,2] [,3]
[1,] 2 2 1
[2,] 1 2 1
[3,] 2 2 2
首先我把它都为矢量:
idVec <- as.vector(idx)
d3mat <- apply(ab,3,as.vector)
然后我构建了一个矢量,以遍历每个元素:
len <- length(idVec)
len <- 1:len
然后我apllied一个函数来向量LEN的每个元素:
resultvec <- sapply(len, function(x, vec, mat){return(mat[x,vec[x]])}, mat=d3mat, vec=idVec)
然后变换结果返回给一个矩阵:
matrix(resultVec,3,3)
[,1] [,2] [,3]
[1,] 9 6 7
[2,] 2 5 8
[3,] 7 4 1