2015-11-30 65 views
4

我有一个相关矩阵,我试图保持每对(行/列)的最大值(考虑绝对值)。我想问一下,如果我具有特定最大值的位置索引,如何提取值。值。 这是我的示例:在数据帧中选择具有索引的列和行值

mat <- structure(c(0, 0.428291512801413, 0.124436112431533, -0.345870125921382, 
      0.391613957773281, 0.428291512801413, 0, 0.341415068127906, -0.346724601510298, 
      0.486360835614514, 0.124436112431533, 0.341415068127906, 0, -0.496213980990412, 
      0.41819049956841, -0.345870125921382, -0.346724601510298, -0.496213980990412, 
      0, -0.80231408836218, 0.391613957773281, 0.486360835614514, 0.41819049956841, 
      -0.80231408836218, 0), .Dim = c(5L, 5L), .Dimnames = list(c("LO3","Tx", "Gh", "RH", "SR"), c("LO3", "Tx", "Gh", "RH", "SR"))) 

然后,我以最大价值的指标:这使我

ind <- apply(abs(mat), 2, which.max) 

LO3 Tx Gh RH SR 
2 5 4 5 4 

我现在想要的东西,它得到的这些位置的值为每列。 这将是:

LO3  Tx  Gh 
0.4282915 0.4863608 -0.4962140 ..... 

我试图使用apply,但我不知道该怎么做。或者如果还有其他方法可以做到这一点。

+3

只是'垫[cbind(1:nrow(mat),ind)]' –

回答

3

既然你有你的指数在ind一种方式可以是使用mapply

#the first argument is the function call 
#second argument is your matrix coerced to data.frame 
#third argument is your indices 
#each time an index will be used in conjunction to a column 
#and you get your result 
mapply(function(x,y) x[y], as.data.frame(mat), ind) 
#  LO3   Tx   Gh   RH   SR 
# 0.4282915 0.4863608 -0.4962140 -0.8023141 -0.8023141 
+0

非常感谢,这是完美的! – user3231352

+0

非常欢迎@ user3231352。很高兴我可以帮忙:) – LyzandeR

0

这能为你做到这一点:

mapply(function(i,j) sample[i,j], seq_len(ncol(sample)), ind) 

> mapply(function(i,j) sample[i,j], seq_len(ncol(sample)), ind) 
[1] 0.4282915 0.4863608 -0.4962140 -0.8023141 -0.8023141 

,如果你愿意,你可以设置的名称结果来自ind

相关问题