如何在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