如何在R语言lmer中添加单个预测变量E与B、C、D的交互项？

包含主效应的简便写法 #

如果需要同时保留E、B、C、D的主效应，以及E分别与B、C、D的双向交互项，直接使用分组公式即可：

lmer(MyVar ~ (1|ID) + E*(B + C + D), data = data, REML = FALSE)

这个公式会自动展开为：

MyVar ~ (1|ID) + E + B + C + D + E:B + E:C + E:D

完全覆盖需求，无需逐个硬编码*或:。

仅保留交互项（不含主效应） #

如果不需要主效应，只需要E与B、C、D的交互项，用:运算符分组：

lmer(MyVar ~ (1|ID) + E:(B + C + D), data = data, REML = FALSE)

展开后仅包含三个交互项：E:B、E:C、E:D。

错误原因说明 #

你之前尝试的(B + C + D)^E写法无效，因为R公式中的幂次运算符^仅支持数字（比如^2表示所有双向交互），不能用变量作为幂次参数。

内容的提问来源于stack exchange，提问作者Larissa Cury