permute中dbrda/PERMANOVA的三因子嵌套重复测量设计置换限制

实验设计背景 #

• 3个样地（区组1、2、3），每个样地包含3个不同火烧强度（高、中、低）的样方
• 每个样方内设置10个亚样方：5个处理组、5个未处理组；处理水平覆盖所有火烧强度，火烧强度水平覆盖所有样地
• 每个亚样方在两年内被测量4次，核心关注火烧强度×处理的交互作用随时间对物种组成的影响，需控制样地区组差异

核心问题解答 #

1. 单个时间点的置换结构是否恰当？ #

你当前的置换设置存在逻辑漏洞，仅用blocks = dat$Site会允许同一样地内不同火烧强度组的亚样方互相置换，这不符合实验的随机化逻辑——Treatment是嵌套在Burn_Severity组内的，置换必须限制在同一样地+同一火烧强度的亚样方之间进行。

修正后的代码如下：

# 构建精准分组变量：样地与火烧强度的组合
dat$Site_Burn <- interaction(dat$Site, dat$Burn_Severity)

# 设置置换规则：仅在同一Site_Burn组内置换亚样方
ctrl <- how(
  nperm = 999,
  blocks = dat$Site_Burn
)

# 运行dbRDA，Condition(Site)控制区组效应
dbrda(sp.mat ~ Treatment*Burn_Severity + Condition(Site),
      data = dat,
      permutations = ctrl,
      method = "bray")

该设置能保证置换仅发生在匹配实验设计的亚组内，避免无效跨组置换导致的假阳性结果。

2. 整合重复测量（Time因素）的置换设置 #

重复测量数据的核心要求是：同一亚样方的4个时间点属于绑定观测，不能拆分置换，同时保留组内（Site+Burn_Severity）的亚样方置换逻辑，以下是两种可行方案：

方案1：用strata锁定亚样方的重复观测

# 构建分组与亚样方标识变量
dat$Site_Burn <- interaction(dat$Site, dat$Burn_Severity)
dat$Subplot_ID <- factor(paste(dat$Site, dat$Burn_Severity, dat$Subplot, sep = "_"))

# 设置置换规则：Site_Burn内置换，且同一亚样方的时间点绑定
ctrl <- how(
  nperm = 999,
  blocks = dat$Site_Burn,
  strata = dat$Subplot_ID,
  within = Within(type = "free", mirror = FALSE)
)

# 包含时间交互项的dbRDA模型
dbrda(sp.mat ~ Treatment*Burn_Severity*Time + Condition(Site),
      data = dat,
      permutations = ctrl,
      method = "bray")

方案2：结合混合效应框架（adonis2）

如果更强调亚样方的随机效应，可采用adonis2配合置换限制：

dat$Site_Burn <- interaction(dat$Site, dat$Burn_Severity)

ctrl <- how(
  nperm = 999,
  blocks = dat$Site_Burn,
  strata = dat$Subplot_ID
)

# 将Subplot_ID设为随机效应，控制重复测量的非独立性
adonis2(sp.mat ~ Treatment*Burn_Severity*Time + Condition(Site) + (1|Subplot_ID),
        data = dat,
        permutations = ctrl,
        method = "bray")

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