R - aov() - 只做双向交互

Question

如果我在 R 中使用 aov() 运行 ANOVA，如何将交互项限制为仅双向交互？例如这段代码：

model <- aov(y ~ a * b * c, data=data)
summary(model)

将创建一个包含 a、b、c、ab、ac、bc 和 abc 的模型。如何限制这一点，使模型中出现的唯一交互是 ab、ac 和 bc？我试过了

aov(y ~ a .*. b .*. c, data = data)

但这引发了错误。

Answer 1

在公式中，^ 运算符将交互项限制为指定的顺序：

aov(y ~ (a + b + c)^2, ...)

Answer 2

您可以通过以下方式明确指定要包含在模型中的术语

model <- aov(y ~ a + b + c + a:b + b:c + a:c, data=data)