R语言如何将变量名列表批量传入lm函数构建线性回归交互项

你要的效果完全可以通过原生R的公式动态构造实现，无需额外安装第三方包，实现步骤如下：

步骤1：提取所有X组变量的变量名 #

如果你的X变量命名规律（X1到X100），可以直接生成：

x_vars <- paste0("X", 1:100)

如果X变量命名不规律，可以从数据集列名中匹配提取：

# 示例：提取所有以"X"开头的变量作为分组变量
x_vars <- grep("^X", names(data), value = TRUE)

步骤2：动态构造回归公式 #

用reformulate函数批量生成所有交互项，自动拼接为合法的回归公式：

# 构造公式右侧项：主效应C + A与所有X的双向交互、B与所有X的双向交互
# 注：R中`A*X1`等价于`A + X1 + A:X1`，符合你原公式的项定义
rhs_terms <- c("C", paste0("A*", x_vars), paste0("B*", x_vars))
reg_formula <- reformulate(rhs_terms, response = "Y")

步骤3：拟合模型 #

直接把构造好的公式传入lm即可：

model <- lm(reg_formula, data = data)

补充说明 #

• R的公式语法会自动识别重复的主效应项（比如X1在A*X1和B*X1中各出现一次），不会重复计算，无需额外处理
• 如果仅需要交互项不需要X的主效应，可以把A*替换为A:，即可仅生成交互项
• 大量交互项会快速消耗自由度，拟合前请确认样本量足够，避免过拟合或多重共线性问题

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