Question

我有一个相关矩阵，我试图保持每对（行/列）的最大值（考虑绝对值）。我想问一下，如果我有具有特定最大值的位置索引，如何提取值。值。这是我的样本：

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，但我不知道该怎么做。或者如果还有其他方法可以做到这一点。

Answer 1

由于您的索引位于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

Answer 2

这可以为你做到：

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

的结果名称

使用数据框中的索引选择列和行值

2 个答案: