自动从线性模型中去除不显着的术语

Question

我正在研究一个包含分类和连续自变量的数据集，并希望找出最小的适当模型。

这是首发模式：

mod1 <- lm(Richness ~ Distances*Flower*Veg*Canopy*factor(Vines), data = Data)
anova(mod1)

然后我通过手动过程删除不显着的交互术语和变量，如下所示：

mod2 <- update(mod1, ~.-Canopy:factor(Vines))
anova(mod2)

有没有办法让这个过程自动化？

感谢您的帮助。

Answer 1

根据@nongkrong关于dropEvery函数的建议，我认为我找到了一个解决方案：

step

到目前为止似乎有效，如果有人发现此代码存在问题，请发表评论。