将线性模型限制为 R 中的第三级交互

Question

我有一个包含 14 个二进制变量的数据集。我已经测试了重要的单个变量，但我还想检查重要的交互作用。然而，我知道更高级别的交互不太可能很重要，而且只会混淆模型。是否有办法在 R 中运行线性模型，但告诉它仅测试最多 3 个变量之间的相互作用？

Answer 1

使用内置 anscombe 数据集的前 5 列：

lm(y1 ~ .^3, anscombe[1:5])

给予：

Call:
lm(formula = y1 ~ .^3, data = anscombe[1:5])

Coefficients:
(Intercept)           x1           x2           x3           x4        x1:x2  
   12.81992     -2.60371           NA           NA     -0.16258      0.36279  
      x1:x3        x1:x4        x2:x3        x2:x4        x3:x4     x1:x2:x3  
         NA           NA           NA           NA           NA     -0.01345  
   x1:x2:x4     x1:x3:x4     x2:x3:x4  
         NA           NA           NA