如何在R中的线性模型中包含所有可能的双向交互项？

Question

是否有一种简单的方法可以在R？模型中包含所有可能的双向交互？

鉴于此模型：

lm(a~b+c+d)

将使用什么语法使模型包含b，c，d，bc，bd和cd作为解释变量，bc是主效应b和c的交互项。

Answer 1

您可以写下以下内容：

lm(a ~ (b + c + d)^2)

这会创建b，c和d之间双向互动的所有组合

例如：

lm(mpg ~ (cyl+disp+hp)^2, data = mtcars)

给出：

Call:
lm(formula = mpg ~ (cyl + disp + hp)^2, data = mtcars)

Coefficients:
(Intercept)          cyl         disp           hp     cyl:disp       cyl:hp      disp:hp  
  5.601e+01   -4.427e+00   -1.184e-01   -1.142e-01    1.439e-02    1.556e-02   -8.567e-05