我有一个数字类型的矩阵,dim 10000 * 50.现在我想按照它们的值的顺序找到每行中前5个元素的索引。例如样本可能如下所示:
set.seed(2)
v1 <- matrix(runif(20 , 0 ,20 ) , 2 ,10)
v1
# [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
#[1,] 3.697645 11.466527 18.87679 2.58318 9.36037 11.053481 15.210266 8.105644 19.527970 8.896185
#[2,] 14.047481 3.361038 18.86950 16.66898 10.99967 4.777895 3.616402 17.070969 4.516509 1.499588
然后我希望输出看起来像:
#[1,] 9 3 7 2 6
#[2,] 3 8 4 1 5
我只能找到this question,它解释了如何找到前n个元素,但不是按值的顺序。
答案 0 :(得分:3)
apply()
非常适合矩阵上的逐行操作。你可以做到
t(apply(v1, 1, function(x) order(-x)[1:5]))
# [,1] [,2] [,3] [,4] [,5]
# [1,] 9 3 7 2 6
# [2,] 3 8 4 1 5
这会在矩阵order()
下逐行运行v1
函数,然后为每一个值取前五个值,转换结果,因为您指定的行不是列。
答案 1 :(得分:1)
这也可以在data.table
进入&#39; long&#39;之后melting
完成。格式,按&#39; Var1&#39;分组,我们得到order
&#39;值&#39;
library(reshape2)
library(data.table)
setDT(melt(v1))[, head(order(-value),5), Var1]
# Var1 V1
#1: 1 9
#2: 1 3
#3: 1 7
#4: 1 2
#5: 1 6
#6: 2 3
#7: 2 8
#8: 2 4
#9: 2 1
#10: 2 5
或使用base R
ave(-v1, row(v1), FUN = order)[,1:5]
# [,1] [,2] [,3] [,4] [,5]
#[1,] 9 3 7 2 6
#[2,] 3 8 4 1 5