如何在R中更新回归公式以排除交互项（适配多模型）

方法1：动态构造公式字符串（简洁适配多模型） #

一行代码即可移除所有交互项并保留主效应和随机效应，完美适配lm/glm/lmer/glmer：

update(out, as.formula(sprintf(". ~ . - %s", paste(grep(":", attr(terms(out), "term.labels"), value = TRUE), collapse = " - "))))

代码解析：

• grep(":", attr(terms(out), "term.labels"), value = TRUE)：从模型固定效应项中筛选出所有含:的交互项；
• paste(..., collapse = " - ")：将多个交互项拼接成"项1 - 项2 - ..."的字符串格式；
• sprintf(". ~ . - %s", ...)：构造更新公式，. ~ .会保留原模型的响应变量、主效应和随机效应，仅移除指定交互项；
• update(out, ...)：用新公式更新原模型，自动继承原模型的所有参数（如数据、分布族等）。

方法2：使用reformulate构造公式（更直观控制结构） #

如果需要更清晰地掌控公式组成，可以用reformulate直接生成仅含主效应和随机效应的公式：

update(out, reformulate(c(grep(":", attr(terms(out), "term.labels"), invert = TRUE, value = TRUE), sapply(lme4::findbars(formula(out)), deparse)), response = formula(out)[[2]]))

代码解析：

• grep(":", ..., invert = TRUE)：筛选出所有不含:的主效应项；
• sapply(lme4::findbars(formula(out)), deparse)：提取原模型中的随机效应部分（如(1|site)）并转为字符串；
• reformulate(...)：将主效应、随机效应与响应变量组合成新公式；
• update(out, ...)：更新模型，最终效果与方法1完全一致。

为什么你之前的方法失效？ #

1. drop.terms方法：该函数仅处理固定效应的terms对象，会直接丢失原模型的随机效应部分，因此在lmer/glmer中使用必然报错；
2. 直接在公式中使用字符向量：你尝试的update(out, . ~ . -out.fixvars[grep('[:]', out.fixvars)])中，out.fixvars[grep(...)]是一个字符向量，R会将其解析为单个变量名而非多个要移除的交互项，自然无法生效。

测试示例 #

library(lme4)
set.seed(123)
# 模拟测试数据
TXdata <- data.frame(
  abundance = rnorm(100),
  prop_seminatural = rnorm(100),
  flower_richness = rnorm(100),
  bare_ground = rnorm(100),
  site = rep(1:10, each=10)
)
# 原模型
out <- lmer(abundance~prop_seminatural*flower_richness+bare_ground+(1|site), data=TXdata)
# 方法1更新模型
new_out1 <- update(out, as.formula(sprintf(". ~ . - %s", paste(grep(":", attr(terms(out), "term.labels"), value = TRUE), collapse = " - "))))
# 方法2更新模型
new_out2 <- update(out, reformulate(c(grep(":", attr(terms(out), "term.labels"), invert = TRUE, value = TRUE), sapply(lme4::findbars(formula(out)), deparse)), response = formula(out)[[2]]))
# 查看新模型公式
formula(new_out1) # abundance ~ prop_seminatural + flower_richness + bare_ground + (1 | site)
formula(new_out2) # 与new_out1完全一致

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