公式就像这样。我不是很理解符号“/”的用法。它似乎通常用于虚拟变量。但是我对它的用法不确定。
lm(y~x/z,data = data.frame(x = rnorm(6), y = rnorm(6), z = rep(0:1,each=3)))
lm(y ~ x/z, data)
只是lm(y ~ x + x:z, data)
的一种缩写形式。
这两者给出相同的结果:
lm(mpg ~ disp/hp,data = mtcars)
Call:
lm(formula = mpg ~ disp/hp, data = df)
Coefficients:
(Intercept) disp disp:hp
2.932e+01 -3.751e-02 -1.433e-05
lm(mpg ~ disp + disp:hp, data = mtcars)
Call:
lm(formula = mpg ~ disp + disp:hp, data = mtcars)
Coefficients:
(Intercept) disp disp:hp
2.932e+01 -3.751e-02 -1.433e-05
所以,你正在基于 disp
本身和 disp
和 hp
之间的交互来建模 mpg
。