permute包how()函数能否支持不平衡数据适配vegan::adonis2()置换限制？

研究背景与问题 #

我正在开展昆虫群落对比研究，设计为分层不平衡结构：

• 3个林区；
• 每个林区包含2个林窗；
• 每个林窗对应3个样方：1个位于林窗内，2个位于邻近郁闭林（此为数据不平衡的来源）。

目标是用vegan::adonis2()实现PERMANOVA，对比林窗与郁闭林的群落差异。由于设计具有分层性，必须限制置换范围，但看到资料提及permute包要求“数据集必须平衡”，因此有两个核心问题：

1. permute包真的完全不支持不平衡数据吗？
2. 若支持，如何在vegan::adonis2()中用permute定义适配该不平衡设计的置换限制？

模拟数据集 #

# 加载包
library(tibble)
library(vegan)
library(permute)

# 构建模拟数据：行对应样方
data <- tibble(
  # 3个林区，每个林区2个林窗
  forestarea = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3),
  # 6个唯一林窗编号
  gapNo = c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6),
  # 样方生境类型：g=林窗，f=郁闭林（每个林窗对应1g2f）
  gapOrforest = rep(c("g", "f", "f"), 6),
  # 10种昆虫的计数数据
  species1 = c(50, 0, 0, 65, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 42, 0, 0),
  species2 = c(0, 10, 8, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 5, 9, 9),
  species3 = c(30, 0, 10, 0, 0, 23, 0, 0, 0, 45, 0, 0, 23, 10, 0, 43, 21, 60),
  species4 = c(0, 5, 12, 0, 0, 3, 20, 0, 0, 0, 15, 0, 0, 0, 18, 0, 0, 0),
  species5 = c(25, 0, 0, 0, 13, 0, 0, 0, 22, 0, 0, 0, 0, 0, 0, 14, 0, 0),
  species6 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 17, 4, 5, 0, 0, 0, 6, 0),
  species7 = c(14, 0, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 4),
  species8 = c(0, 7, 0, 0, 0, 5, 0, 0, 0, 12, 0, 0, 0, 0, 0, 0, 0, 0),
  species9 = c(0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 5, 0, 0, 0, 0, 0),
  species10 = c(8, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 7, 0, 0, 0)
)

平衡数据下的PERMANOVA示例代码 #

dist <- vegdist(data[,4:ncol(data)])

control <- how(within = Within(type = "free"),
             plots = Plots(strata = data$gapNo, type = "none"),
             nperm = 999,
             observed = TRUE)

adonis2(dist ~ gapNo, 
        data=data,
        permutations = control)

解答 #

1. permute包是否支持不平衡数据？ #

并非完全不支持。所谓“平衡要求”仅针对部分特定置换场景（如完全交叉设计的全局置换）。对于分层嵌套的不平衡设计，只要正确定义置换策略，permute完全可以适配。

2. 适配不平衡设计的置换限制设置 #

针对你的设计（每个林窗内固定1个林窗样方、2个郁闭林样方），核心是限制置换仅在单个林窗内部进行，且保持每个林窗内的生境类型样方数量不变，这样既能控制林窗间的固有差异，又适配了数据的不平衡结构。

修改后的代码如下：

# 计算群落距离矩阵（常用Bray-Curtis距离）
dist <- vegdist(data[,4:ncol(data)], method = "bray")

# 定义置换控制策略
control <- how(
  within = Within(type = "free"),
  # 指定分层单元为每个林窗（gapNo）
  plots = Plots(strata = data$gapNo, type = "none"),
  nperm = 999,
  observed = TRUE,
  # 关键：保持每个分层内的生境类型样方数量，适配不平衡结构
  permControl = permControl(balance = "strata")
)

# 运行PERMANOVA：检验生境类型（gapOrforest）的效应
adonis2(dist ~ gapOrforest, 
        data = data,
        permutations = control,
        method = "bray")

代码说明：

• strata = data$gapNo：将每个林窗设置为独立的置换分层，确保置换仅在同一林窗的样方间进行；
• permControl(balance = "strata")：强制每个分层内的生境类型样方数量与原数据一致（即每个林窗始终保留1个g、2个f样方），完美适配你的不平衡设计；
• 模型公式改为dist ~ gapOrforest，直接检验生境类型对昆虫群落的影响。

如果需要进一步控制更高层级的林区（forestarea）差异，可将分层设置为strata = interaction(data$forestarea, data$gapNo)，但由于你的gapNo是全局唯一的，此步骤非必需。

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