R语言如何对数据框二元变量两两求和并批量生成新列

完全可以用for循环或者apply系列函数实现，两种方案的代码都很简单，适配你说的50176组全配对需求（即所有i、j从1到224的有序两两组合，包含i=j、i>j的场景）：

方案1：for循环实现 #

优先推荐这个写法，逻辑直观，不容易出错，而且提前预分配了空间，效率足够应对你的数据规模：

# 先确认所有变量都是数值类型，如果是因子先转
df[] <- lapply(df, as.numeric)

n <- ncol(df)
# 预分配新列空间，避免动态扩充拉低效率
df[, (n+1):(n + n*n)] <- NA
current_col <- n + 1

for (i in 1:n) {
  for (j in 1:n) {
    df[, current_col] <- df[, i] + df[, j]
    # 给新列命名，方便后续识别对应哪两个变量的和，不需要可以删掉
    colnames(df)[current_col] <- paste0("sum_", colnames(df)[i], "_", colnames(df)[j])
    current_col <- current_col + 1
  }
}

方案2：apply系列函数实现 #

用expand.grid生成所有变量配对，再用sapply批量计算：

df[] <- lapply(df, as.numeric)
n <- ncol(df)
# 生成所有有序两两组合
var_pairs <- expand.grid(x = 1:n, y = 1:n)
# 批量计算所有配对的和
sum_result <- sapply(1:nrow(var_pairs), function(k) {
  df[, var_pairs$x[k]] + df[, var_pairs$y[k]]
})
# 列名绑定，不需要可以删掉
colnames(sum_result) <- paste0("sum_", var_pairs$x, "_", var_pairs$y)
# 合并到原数据框
df <- cbind(df, sum_result)

补充说明 #

如果你的需求不需要保留重复配对（比如1+2和2+1视为同一组，仅保留一次），只需要在生成配对时加过滤条件var_pairs <- var_pairs[var_pairs$x <= var_pairs$y, ]即可，最终生成的新列数量会减少到25200组。你的数据规模很小，88行观测、5万多列新变量，两种方案跑完全程都不会超过1秒，不用额外做性能优化。

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