是否有可以返回R中所有重要交互项的命令?
例如:
输入= c(年龄,性别,国籍)
输出=重要互动列表= c(年龄*性别,性别*国籍,年龄*国籍*性别)
答案 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值是显着性的指标。)