在R语言中用apply函数实现列表对象间的运算（含对角项）

我来帮你优化这个实现，既能覆盖所有需要的项（包括对角项和非对角项），又不会做重复计算——毕竟C_i.j和C_j.i是完全一样的，没必要算两次。下面是具体的实现步骤和代码：

第一步：准备示例数据 #

先复现你提供的测试数据：

set.seed(123)
A <- matrix(c(5.2,6.9,32,40,8.3), ncol = 5 )

B <- list(matrix(rnorm(100 * 100, mean = 0, sd = 1), 100, 100), 
          matrix(rnorm(100 * 100, mean = 0, sd = 1), 100, 100), 
          matrix(rnorm(100 * 100, mean = 0, sd = 1), 100, 100), 
          matrix(rnorm(100 * 100, mean = 0, sd = 1), 100, 100), 
          matrix(rnorm(100 * 100, mean = 0, sd = 1), 100, 100))

第二步：计算对角项 C_i.i #

按照你指定的公式单独计算每个对角项，这部分没有重复的可能：

# 生成所有对角项
diag_list <- lapply(1:ncol(A), function(i) {
  (A[, i] * B[[i]]^2) / 2
})
# 给对角项命名为 C_1.1、C_2.2...
names(diag_list) <- paste0("C_", 1:ncol(A), ".", 1:ncol(A))

第三步：计算非对角项（仅保留i<j的组合） #

用combn生成所有i<j的无序组合，这样只会计算一次C_i.j，完全避免重复计算C_j.i：

# 生成所有i<j的索引组合
non_diag_combs <- combn(1:ncol(A), 2)

# 计算非对角项
non_diag_list <- lapply(1:ncol(non_diag_combs), function(idx) {
  i <- non_diag_combs[1, idx]
  j <- non_diag_combs[2, idx]
  (max(A[, i], A[, j]) * B[[i]] * B[[j]]) / 2
})
# 给非对角项命名为 C_1.2、C_1.3...
names(non_diag_list) <- paste0("C_", non_diag_combs[1, ], ".", non_diag_combs[2, ])

第四步：合并结果列表 #

把对角项和非对角项的列表合并，得到完整的结果：

result_list <- c(diag_list, non_diag_list)

方案优势 #

• 无重复计算：非对角项只计算一次i<j的组合，相比生成所有i和j的全组合，节省了一半的计算量，尤其当B中的矩阵很大或者数量很多时，这个优化效果会更明显。
• 完整覆盖所有需要的项：最终列表包含5个对角项 + 10个非对角项（C(5,2)=10），总共15个元素，没有冗余。
• 命名清晰：所有元素的命名完全符合你的要求，后续调用非常方便（比如需要C_2.1的话，直接取result_list[["C_1.2"]]即可）。

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