如何对R中自变量之间的所有关系建模？

我有一个小数据集，有4个独立变量（称为a、b、c、d）和1个因变量。因为自变量很少，所以我想探讨这些变量的所有组合。只能有14种型号（a、b、c、d、a+b、a+c、a+d、b+c、b+d、c+d、a+b+c、a+b+d、b+c+d、a+b+c+d）。 我手工制作所有模型，这很耗时。所以我想把它自动化。在R中可能吗

glm（dep~a，family=“二项式”，data=data） glm（dep~b+c，family=“二项式”，data=data）

---

etc

使用内置的

anscombe

数据框架，这将找到最佳组合（尽管现在不鼓励逐步回归，支持套索、弹性网和其他基于惩罚的类似方法）

fm
fm <- glm(y1 ~ x1 + x2 + x3 + x4, data = anscombe)
step(fm)