R语言中对三维数组指定维度应用quantile分位数函数的方法

你已构造的combi.prod为维度(100, 150, 100)的三维数组，其中第二维度的150层分别对应150个100×100的目标矩阵，以下是两种常用实现方式：

方法1：基础R向量化实现（效率更高，无需额外依赖） #

直接使用apply函数对指定维度批量应用分位数函数即可：

# 定义需要计算的分位点
target_quant <- c(0.5, 0.9)
# 对第二维度的每一层矩阵计算分位数
quant_res <- apply(
  X = combi.prod,
  MARGIN = 2, # 指定沿第二维度切分
  FUN = function(mat) {
    quantile(mat, probs = target_quant, na.rm = TRUE) # 加na.rm避免存在缺失值时报错
  }
)

返回结果quant_res为2行150列的矩阵，每行对应一个分位点，每列对应一个原矩阵的计算结果。如果需要按矩阵为行展示，执行t(quant_res)转置即可。

方法2：显式循环实现（适合需要自定义中间处理逻辑的场景） #

如果需要在计算分位数前对单个矩阵做额外处理，可以用循环实现：

target_quant <- c(0.5, 0.9)
mat_num <- dim(combi.prod)[2]
# 初始化结果矩阵，行对应150个矩阵，列对应指定分位点
quant_res <- matrix(NA, nrow = mat_num, ncol = length(target_quant))
colnames(quant_res) <- paste0("Q_", target_quant)

for (k in 1:mat_num) {
  # 提取第k个100×100矩阵
  current_mat <- combi.prod[, k, ]
  # 计算分位数
  quant_res[k, ] <- quantile(current_mat, probs = target_quant, na.rm = TRUE)
}

补充说明：如果需要对每个矩阵的行/列单独计算分位数，只需调整内层函数的维度参数即可，比如对每个矩阵的行计算分位数，将quantile(mat, ...)替换为apply(mat, 1, quantile, probs = target_quant)即可。

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