r - 如何在一个循环中拟合多个交互模型？

Question

假设我有 3 个响应变量 A、C 和 M，我想为所有可能的模型拟合一个模型，即拟合 Y ~ A、Y ~ C、Y ~ M、Y ~ A * C、Y ~ A * M、Y ~ C * M 等。有没有一种快速的方法可以做到这一点，而无需每次都手动指定交互？

我不想写

M1 = glm(Y ~ A , data = subs, family = "poisson")
M2 = glm(Y ~ C , data = subs, family = "poisson")
M3 = glm(Y ~ M , data = subs, family = "poisson")
M4 = glm(Y ~ A*C , data = subs, family = "poisson")
...

实际上我有 3 个以上的变量并且想要某种循环，这是否可能。 谢谢

Answer 1

这应该有效:

glmulti::glmulti(
  Y = "Y", 
  xr = c("A", "C", "M"),
  data = subs,
  filename = "my_results",
  family = "poisson"
) 

它将创建一个文本文件 my_results.txt，其中包含有关每个模型的信息。

您也可以将它与其他包一起使用，leaps、bestglm，可能还有其他包。