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

给定此模型:

lm(a~b+c+d)


如果bc是主要效果b和c的相互作用项,将使用什么语法使模型包括b,c,d,bc,bd和cd作为解释变量。

最佳答案

您可以编写以下内容:

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


这将创建bcd之间双向交互的所有组合

例如:

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

关于r - 如何在R的线性模型中包括所有可能的双向相互作用项?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/47144532/

10-12 20:07