如何生成要在GLM中使用的所有可能重要交互项的列表?

时间:2016-09-07 10:24:35

标签: r interaction

是否有可以返回R中所有重要交互项的命令?

例如:

输入= c(年龄,性别,国籍)

输出=重要互动列表= c(年龄*性别,性别*国籍,年龄*国籍*性别)

1 个答案:

答案 0 :(得分:0)

试用broom包:

library(broom)
m = lm(mpg ~ am * hp * disp, data = mtcars)
tidy.m = tidy(m)
tidy.m$term[tidy.m$p.value < 0.05]

这给了你:

[1] "(Intercept)" "am"          "am:disp"     "am:hp:disp" 

这就是你想要的吗?

(假设p值是显着性的指标。)