使用lm()执行含交互项的多元线性回归报错，求助解决方法

你写的公式y ~ x1 + x2 + x1*x2存在语法冗余——在R的公式规则里，x1*x2会自动展开为x1 + x2 + x1:x2（包含两个变量的主效应+交互项），所以你等于重复添加了x1和x2的主效应。虽然这不会直接导致报错，但属于多余写法，推荐用以下两种标准写法：

• 简洁自动版：直接用*符号，R会自动生成主效应和交互项

  data <- data.frame(x1 = rnorm(100), x2 = rnorm(100), y = rnorm(100))
  model <- lm(y ~ x1*x2, data = data)
  summary(model)
  plot(model)
  

• 显式指定版：用:表示纯交互项，搭配主效应一起写，适合需要明确控制变量的场景

  data <- data.frame(x1 = rnorm(100), x2 = rnorm(100), y = rnorm(100))
  model <- lm(y ~ x1 + x2 + x1:x2, data = data)
  summary(model)
  plot(model)
  

如果你实际运行原代码时出现错误，大概率不是公式的x1*x2部分导致的，可能是变量名冲突、数据存在缺失值或其他环境问题，可以检查数据完整性和变量命名。

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