R：使用索引矩阵从数组中提取矩阵

Question

我在R中编码，并且我有一个包含数据的3维数组（在本例中为ab）。然后我有一个包含第三个数组维度（idx）索引的矩阵。该矩阵具有相同数量的行和列数组。我想使用idx中包含的索引从数组中提取数据，以获得具有相同维度idx的矩阵。请看下面的例子：

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] 

有没有人知道我怎么能得到它？

非常感谢！

萨拉

Answer 1

我们需要用于行/列和第三维度中的索引（从“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 

Answer 2

可能有点复杂，但它的作品至少： 首先我所有的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