批量使用buildglmmTMB建模时如何保留指定随机效应？

核心结论 #

默认使用direction = "backward"或direction = "order"参数时，buildglmmTMB不会移除随机效应。它的向后逐步消除逻辑仅针对固定效应，随机效应结构会全程保留在模型中，不会被纳入筛选范围。

确保随机效应不被移除的实操方法 #

• 维持默认方向参数：继续使用direction = "backward"或direction = "order"，这两种模式下随机效应不会被改动。
• 明确固定随机效应结构：在模型公式中显式写出所有需要保留的随机效应，例如：

  model <- buildglmmTMB(
    y ~ var1 + var2 + var3 + ... + (1|group1) + (1|group2),
    data = group_data,
    direction = "backward"
  )
  

  这种写法下，逐步选择过程只会对var1、var2等固定效应进行筛选，随机效应(1|group1)和(1|group2)会被强制保留。
• 验证模型一致性：循环处理完每个类群后，可通过summary(model)或ranef(model)检查随机效应是否存在，确保所有模型的随机效应结构统一，满足跨类群比较的要求。

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