如何在R中构建交互回归模型：自动档线性/手动档1/x响应

数据预处理与观察 #

首先将mtcars处理为mtcars2：

mtcars2 <- within(mtcars, {
  vs <- factor(vs, labels = c("V", "S"))
  am <- factor(am, labels = c("automatic", "manual"))
  cyl  <- ordered(cyl)
  gear <- ordered(gear)
  carb <- ordered(carb)
})

通过分面散点图观察响应趋势：

ggplot(mtcars2, aes(x=disp, y=mpg))+
  geom_point()+
  facet_grid(~am)

观察得到：

• 自动档（am=automatic）：mpg与disp呈线性关系
• 手动档（am=manual）：mpg与disp呈1/x形式的响应

将automatic设为基准水平：

mtcars2 <- mtcars2 %>%
  mutate(am = relevel(am, ref = "automatic"))

遇到的问题 #

尝试构建模型时写出如下代码，但生成了不需要的单独1/disp项：

model4disp<-lm(data=mtcars2, mpg~am + disp +disp*am +I(1/disp)*am)

期望的数学模型为：
y = b0 + b1x1 + b2x2 + b3x1x2 + b4(1/x2)x1
其中：

• x1：虚拟变量，manual时取1，automatic时取0
• x2：disp，y：mpg

需求是仅手动档组使用1/disp的响应关系，自动档组保持线性响应。

正确的模型公式写法如下，可精准匹配你的需求：

model4disp <- lm(data = mtcars2, mpg ~ disp + am:(disp + I(1/disp)))

公式解析 #

• disp：对应自动档组的线性项，即模型中的b0 + b2x2
• am:(disp + I(1/disp))：表示手动档组（am的非基准水平）分别与disp、I(1/disp)构建交互项，对应模型中的b1x1 + b3x1x2 + b4(1/x2)x1

这种写法不会生成单独的1/disp项，只有手动档组才会包含disp交互项和1/disp交互项，完全符合你期望的分组响应逻辑。

验证模型结构 #

可以通过以下代码查看设计矩阵，确认变量结构是否符合预期：

head(model.matrix(model4disp))

输出会显示：

• 自动档组仅包含截距项和disp
• 手动档组额外包含ammanual:disp和ammanual:I(1/disp)两个交互项

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