2

我在 R 中有一个包含 3 列的数组,如下所示:

>x
[,1]     [,2] [,3]
V1 7.803550 0.000000    0
V2 7.842899 0.000000    0
V3 7.310273 0.000000    0
V4 7.230438 0.000000    0
V5 6.582147 1.097024    0
V6 9.364152 0.000000    0

对于我的数组切片,我想要每一行的不同列,所以我想要的索引是:

rows<-1:6
cols<-c(1,1,1,2,3,1)

当我尝试在 R 中获取适当的数组切片时,而不是获取具有适当值的向量,

> t
      V1       V2       V3       V4       V5       V6 
7.803550 7.842899 7.310273 0.000000 0.000000 9.364152 

我得到一个 6x6 的输出矩阵,其中为所有选定的行打印出给定的列

> t
       [,1]     [,2]     [,3]     [,4] [,5]     [,6]
V1 7.803550 7.803550 7.803550 0.000000    0 7.803550
V2 7.842899 7.842899 7.842899 0.000000    0 7.842899
V3 7.310273 7.310273 7.310273 0.000000    0 7.310273
V4 7.230438 7.230438 7.230438 0.000000    0 7.230438
V5 6.582147 6.582147 6.582147 1.097024    0 6.582147
V6 9.364152 9.364152 9.364152 0.000000    0 9.364152

如何在没有 for 循环的向量中获得所需的索引对?

根据大卫的建议,我希望能工作的代码是

t<-x[rows,cols]
4

2 回答 2

10

怎么样:

x[cbind(rows, cols)]

读入?“[”为什么这应该有效。(我认为它比罗宾逊的 sapply 策略更快。)

于 2012-08-30T04:14:58.443 回答
1

可能不能做得更好

sapply(1:length(rows), function(i) x[rows[i], cols[i]])
于 2012-08-30T03:59:34.440 回答