如何无需循环对矩阵列表执行逐元素取最大值操作？

你拥有一个同维度矩阵组成的列表，想要获取各矩阵对应位置的元素最大值，这里提供几种无需显式循环的实现方式：

方法1：do.call() 配合 pmax() #

pmax 支持对多个同维度矩阵逐元素取最大值，do.call 可以将列表中的所有矩阵作为参数传递给 pmax，直接得到结果矩阵：

mat <- list(matrix(1:6,2,3),matrix(5:10,2,3),matrix(4:9,2,3))
result <- do.call(pmax, mat)
result

输出结果：

[,1] [,2] [,3]
[1,]    5    7    9
[2,]    6    8   10

方法2：Reduce() 迭代应用 pmax() #

Reduce 会迭代地将 pmax 应用于列表元素，依次将当前的逐元素最大值结果与下一个矩阵取最大值，最终得到全局结果：

result <- Reduce(pmax, mat)
result

该方法逻辑直观，适合处理元素较多的矩阵列表。

方法3：转三维数组后用 apply() #

先将矩阵列表转换为三维数组，再针对每个位置（行+列维度）计算所有矩阵的最大值：

# 构造三维数组，维度为：原矩阵行 × 原矩阵列 × 列表长度
arr <- array(unlist(mat), dim = c(dim(mat[[1]]), length(mat)))
# 对行、列维度取最大值
result <- apply(arr, c(1, 2), max)
result

以上三种方法均无需编写显式循环，且在处理大规模数据时效率优于嵌套循环实现。

内容的提问来源于stack exchange，提问作者Triparna Poddar