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

Question

我有一个小型数据集，其中包含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 =“ binomial”，数据=数据） glm（dep〜b + c，family =“ binomial”，数据=数据）  等

Answer 1

使用内置的anscombe数据框，它将找到最佳的组合（尽管如今不鼓励逐步回归，而倾向于套索，弹性网和其他基于惩罚的类似方法）。

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

另请参阅jumps和glmnet软件包。