为r中的索引对子集数组

赛义德·奥米德采购

尽管已经进行了搜索,但是找不到我的问题的简单答案。假设我有一个数组:

    vector1 <- c(5,9,3)
vector2 <- c(10,11,12,13,14,15)
result <- array(c(vector1,vector2),dim = c(3,3,2))

现在,我想以某种方式对这个数组进行子集化,以便从某些行和列中获取元素。例如:

result[1,3,1:2]
result[3,1,1:2]

因为我有很多索引,所以它们深陷其中

rowind=c(1,3)
colind=c(3,1)

对于子集,我尝试以这种方式使用行和列的向量

dim(result[rowind,colind,1:2])

[1] 2 2 2

这不是我想要的。我希望我的输出是每对索引的一个值。实际上,我希望提取的矩阵的尺寸为2x2(而不是2x2x2)。让您知道我尝试了“循环”。但是问题是我有很多数组,这将非常耗时。

result_c=matrix(NA,nrow=2,ncol=2)

for (i in 1:2){
  result_c[i,]=result[rowind[i],colind[i],] 
}

谢谢你的帮助。

莫里斯·埃弗斯(Maurits Evers)

只要我对您的理解正确,我们就可以使用 Map

Map(function(i, j) result[i, j, 1:2], rowind, colind)
#[[1]]
#[1] 13 13
#
#[[2]]
#[1] 3 3

或简化为matrix使用mapply

mapply(function(i, j) result[i, j, 1:2], rowind, colind)
#     [,1] [,2]
#[1,]   13    3
#[2,]   13    3

Mapmapply不简化结果的包装器,从而产生一个list(而不是matrix)。

本文收集自互联网,转载请注明来源。

如有侵权,请联系 [email protected] 删除。

编辑于
0

我来说两句

0 条评论
登录 后参与评论

相关文章