Как смоделировать все отношения между независимыми переменными в 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", data = data) glm (dep ~ b + c, family = "binomial", data = data) и т. Д.

Answer 1

Используя встроенный фрейм данных anscombe, он найдет наилучшую комбинацию (хотя в наши дни не рекомендуется пошаговая регрессия в пользу лассо, эластик и других подобных подходов, основанных на штрафах).

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

Также смотрите пакеты leaps и glmnet.