Блог

Есть ли способ в R подогнать термин взаимодействия к каждой предыдущей переменной в спецификации модели? Я хотел бы сделать следующее, но в более краткой форме.

 
data("mtcars")
head(mtcars)
mod1<-lm(mpg~ cyl disp hp wt cyl:wt disp:wt hp:wt, data=mtcars)
summary(mod1)

  

Это похоже на арифметику.

 f1 <- lm(mpg ~ cyl   disp   hp   wt   cyl:wt   disp:wt   hp:wt, data=mtcars)
f2 <- lm(mpg ~ (cyl   disp   hp)*wt, data=mtcars)

stopifnot(all.equal(f1$coe, f2$coe))
  

Я думаю, это должно сработать.

 mod1<-lm(mpg~ (cyl disp hp wt)^2, data=mtcars)