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

Question

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

例如：

输入= c（年龄，性别，国籍）

输出=重要互动列表= c（年龄*性别，性别*国籍，年龄*国籍*性别）

Answer 1

试用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值是显着性的指标。）