R语言MuMIn包dredge函数子集设置：限定交互项包含规则

要实现你需要的模型筛选逻辑，关键是利用dredge()的subset参数定义自定义的模型包含/排除规则。以下是具体的实现步骤和代码：

核心思路 #

对于每个环境变量（Forest、Scrub、Meadow、RiverLgthm、DailySunHrs），我们需要确保：只有当该变量与GISalt的交互项存在时，才允许包含它与I(GISalt^2)的交互项。用逻辑表达式描述就是：如果模型包含I(GISalt^2):X（X为对应环境变量），那么必须同时包含GISalt:X；反之，如果模型不包含I(GISalt^2):X，则对GISalt:X的存在与否没有限制。

同时，从你给出的示例模型来看，MLE200405是所有模型都必须保留的变量，所以我们也要把这一点加入规则。

具体代码 #

library(MuMIn)

# 定义你的全局模型
globmod <- glm(MLE2017 ~ MLE200405 + GISalt * Forest + GISalt * Scrub + GISalt * Meadow + 
                 GISalt * RiverLgthm + GISalt * DailySunHrs + I(GISalt^2) * Forest + 
                 I(GISalt^2) * Scrub + I(GISalt^2) * Meadow + I(GISalt^2) * RiverLgthm + 
                 I(GISalt^2) * DailySunHrs, data = GLMdata, family = quasipoisson(link = "log"))

# 构建子集筛选规则
subset_conditions <- expression(
  # 强制保留MLE200405在所有模型中
  MLE200405 &
  # Forest的交互项规则：I(GISalt^2):Forest 仅在 GISalt:Forest 存在时允许
  (!(I(GISalt^2):Forest) | (GISalt:Forest)) &
  # Scrub的交互项规则
  (!(I(GISalt^2):Scrub) | (GISalt:Scrub)) &
  # Meadow的交互项规则
  (!(I(GISalt^2):Meadow) | (GISalt:Meadow)) &
  # RiverLgthm的交互项规则
  (!(I(GISalt^2):RiverLgthm) | (GISalt:RiverLgthm)) &
  # DailySunHrs的交互项规则
  (!(I(GISalt^2):DailySunHrs) | (GISalt:DailySunHrs))
)

# 运行dredge并应用规则
dredged_result <- dredge(globmod, subset = subset_conditions)

规则验证 #

• 你的mod1（同时包含GISalt*Forest和I(GISalt^2)*Forest）会被保留：因为GISalt:Forest存在，所以允许I(GISalt^2):Forest存在。
• 你的mod2（仅包含GISalt*Forest）会被保留：没有包含I(GISalt^2):Forest，规则不限制GISalt:Forest的存在。
• 你的mod3（仅包含I(GISalt^2)*Forest）会被排除：此时I(GISalt^2):Forest存在，但GISalt:Forest不存在，违反了规则中的逻辑条件。

额外说明 #

• dredge()默认开启marginality = TRUE，这会确保交互项的主效应（比如Forest、GISalt本身）会在交互项存在时自动被包含，符合统计建模的边际性原则。如果你不需要这个行为，可以手动设置marginality = FALSE，但一般不建议这么做。
• 如果你需要调整规则（比如对主效应也做限制），可以在subset_conditions中添加对应的逻辑表达式即可。

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