如何在 r 的列表中找到最大值的 i j k 索引？

Question

我在 R 中有一个列表，其中有数据 (Data)：

#create data
tmintest=array(1:100, c(12,256,512))

#create the list
Variable <- list(varName = c("tmin","tmin","tmin","tmin","tmin","tmin","tmin","tmin","tmin","tmin","tmin","tmin"),level = c(NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA))
Data     <- tmintest
xyCoords <- list(x = seq(-40.37,64.37,length.out=420), y = seq(25.37,72.37,length.out=189))
Dates <- list(start = seq(as.Date("2012-01-01"), as.Date("2015-12-31"), by="days"), end=seq(as.Date("2012-01-01"), as.Date("2015-12-31"), by="days"))
All <- list(Variable = Variable,Data=Data, xyCoords=xyCoords,Dates=Dates)

我怎样才能准确地找到ALL$Data 中最大数字的位置？例如，如果它是第一个“切片”或“网格”的第 4 行第 100 列，我想返回：[1,4,100]。

我试过 which.max(All$Data) 但只返回一个数字？

Answer 1

which 支持arr.ind，它返回条件的数组索引。不幸的是，which.max 没有这样的参数，因此我们可以将值与值的最大值进行比较。

head( which(All$Data == max(All$Data), arr.ind = TRUE) )
#      dim1 dim2 dim3
# [1,]    4    9    1
# [2,]    8   17    1
# [3,]   12   25    1
# [4,]    4   34    1
# [5,]    8   42    1
# [6,]   12   50    1

在这里我会有点谨慎：当精度涉及大量小数时，浮点数的严格相等测试可能会成为问题。请参阅 Why are these numbers not equal?、Is floating point math broken? 和 https://en.wikipedia.org/wiki/IEEE_754 以获得对此的良好讨论。

更好的测试是严格的不等式之一，寻找一个容忍度。在这里我将使用1e-5，因为我们知道它明显小于数字范围（1 到 100），但是如果您的实数具有更高的精度，您可能需要更适合您需要的东西。

head( which( (max(All$Data) - All$Data) < 1e-5, arr.ind = TRUE) )
#      dim1 dim2 dim3
# [1,]    4    9    1
# [2,]    8   17    1
# [3,]   12   25    1
# [4,]    4   34    1
# [5,]    8   42    1
# [6,]   12   50    1

请注意，如果您将公差值 1e-5过低，您可能会开始丢失值。它不会在这里发生（因为您的数据要大得多）。