R语言线性回归：如何循环双变量构建成对回归模型？

我来帮你搞定这个循环线性回归的问题——你想要的是自动生成所有Test类变量的两两配对，然后对每一组配对运行Y ~ X的线性回归，不用手动逐个替换变量，对吧？其实核心思路是先拿到所有目标变量，生成所有可能的组合，再循环构建模型就行，下面一步步来：

第一步：提取所有Test类变量 #

首先先把数据里所有以Test开头的变量筛选出来（如果你的变量命名规则不是Test开头，也可以用位置索引，比如你之前用的names(df)[3:7]，不过用名称匹配更灵活）：

df <- read.csv("test.csv")
# 提取所有Test开头的变量，返回变量名向量
test_vars <- grep("^Test", names(df), value = TRUE)
# 如果是按位置选，比如Test变量在第2到Control变量之前的列，就用：
# test_vars <- names(df)[2:(which(names(df)=="Control1")-1)]

第二步：生成所有两两变量组合 #

接下来要生成所有可能的Y和X配对，分两种情况：

情况1：需要所有有序组合（比如同时跑Test1_{Test2和Test2}Test1） #

用expand.grid生成所有有序配对，再去掉Y和X相同的无效组合：

all_pairs <- expand.grid(y_var = test_vars, x_var = test_vars, stringsAsFactors = FALSE)
# 剔除Y和X是同一个变量的情况
all_pairs <- all_pairs[all_pairs$y_var != all_pairs$x_var, ]

情况2：只需要无序的唯一组合（比如只保留Test1_{Test2，不重复跑Test2}Test1） #

用combn生成无序的两两组合，更高效：

# 生成2个变量的所有组合，转成数据框方便后续处理
all_pairs <- t(combn(test_vars, 2))
colnames(all_pairs) <- c("y_var", "x_var")
all_pairs <- as.data.frame(all_pairs, stringsAsFactors = FALSE)

第三步：循环构建回归模型 #

现在就可以循环每一组配对，动态生成公式并运行回归了：

# 用lapply循环每一行组合
models <- lapply(1:nrow(all_pairs), function(i) {
  # 取出当前组的Y和X变量名
  current_y <- all_pairs$y_var[i]
  current_x <- all_pairs$x_var[i]
  # 动态构建公式：把变量名拼接成"Y ~ X"的字符串，转成formula对象
  reg_formula <- as.formula(paste(current_y, "~", current_x))
  # 运行线性回归
  lm(reg_formula, data = df)
})

# 给每个模型命名，方便后续快速定位（比如"Test1 ~ Test2"）
names(models) <- paste(all_pairs$y_var, "~", all_pairs$x_var)

# 查看单个模型结果，比如Test1和Test2的回归
models[["Test1 ~ Test2"]]
# 查看所有模型的摘要统计
lapply(models, summary)

额外小技巧：把结果整理成数据框 #

如果想把所有模型的关键结果（比如系数、p值、R²）整理成结构化的数据框，推荐用broom包，非常方便：

install.packages("broom") # 第一次用先安装
library(broom)

# 提取每个模型的系数细节（截距、自变量系数、p值等）
coef_results <- lapply(models, tidy)
coef_df <- do.call(rbind, coef_results)
# 添加模型名称列，方便对应
coef_df$model <- rep(names(models), each = 2) # 每个模型有2行结果：截距和自变量

# 提取模型的整体统计量（R²、调整R²、F值等）
model_stats <- lapply(models, glance)
stats_df <- do.call(rbind, model_stats)
stats_df$model <- names(models)

这样哪怕你的数据集有上千个Test变量，也能自动完成所有两两组合的回归，完全不用手动替换变量啦。

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