使用R语言Reduce函数求矩阵和结果错误，请求排查代码问题

主要错误点 #

• 循环内重复计算求和：原代码在每次循环中都调用Reduce求和，虽不影响最终结果，但效率低下且逻辑冗余。
• 需求与实现可能不符：你提到矩阵由多元矩阵的列生成，但原代码处理的是矩阵的行（x[i,]），若需处理列需调整索引。
• 代码冗余易出错：原代码用繁琐的矩阵构造方式计算外积，增加出错概率。
• 未初始化列表：动态创建列表可能导致意外行为，提前初始化更稳妥。

修正后的代码（处理行的情况） #

library(MASS)
# 生成5行5列的多元正态数据
x <- mvrnorm(5, mu = rep(0,5), diag(5))
p <- ncol(x) # 自动获取列数，替代硬编码
n <- nrow(x)

# 初始化空列表存储每个生成的矩阵
matrices <- list()

# 循环生成矩阵
for(i in 1:n){
  xi <- x[i,]
  norm_xi <- norm(xi, type = "2")
  # 简化外积计算
  proj_matrix <- tcrossprod(xi) / (norm_xi^2)
  # 构造目标矩阵并加入列表
  matrices[[i]] <- (diag(p) - proj_matrix) / norm_xi
}

# 一次性对所有矩阵求和
sum_matrices <- Reduce("+", matrices)
sum_matrices

若需处理列的情况 #

只需调整索引和维度定义：

library(MASS)
x <- mvrnorm(5, mu = rep(0,5), diag(5))
p <- nrow(x) 
n <- ncol(x)

matrices <- list()

for(i in 1:n){
  xi <- x[,i] # 改为取列向量
  norm_xi <- norm(xi, type = "2")
  proj_matrix <- tcrossprod(xi) / (norm_xi^2)
  matrices[[i]] <- (diag(p) - proj_matrix) / norm_xi
}

sum_matrices <- Reduce("+", matrices)
sum_matrices

修正说明 #

• 将Reduce求和移至循环外，仅计算一次，逻辑更清晰。
• 用tcrossprod(xi)替代繁琐的矩阵转置相乘，简化外积计算。
• 单独定义xi、norm_xi变量，提升代码可读性与调试便利性。
• 用ncol(x)/nrow(x)替代硬编码数值，让代码适配不同维度的输入矩阵。

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