PERMANOVA-adonis2报错排查及后续统计分析咨询

背景概述 #

我有9个样地的生态学丰富度（richness）数据，分为3组（每组3个样地）：前2组包含5年的丰富度数据，第3组仅有3年数据。目标是检验时间（年份）、样地（组内）、组的主效应及所有交互作用。

遇到的问题 #

使用vegan包的adonis2函数进行分析时，出现错误提示：contrasts can only be applied to factors with 2 or more levels，需要排查原因。

同时存在疑问：若任一因子检验出显著效应，是否应该采用t检验？之前尝试过用easynova包处理非平衡样本，结果显示无差异，转而使用adonis2函数。

附数据与代码 #

数据集 #

dados <- read.table(text = "
Linha Parcela Ano Riqueza
1 1 a1 8
1 1 a2 8
1 1 a3 12
1 1 a4 12
1 1 a5 12
1 2 a1 12
1 2 a2 12
1 2 a3 13
1 2 a4 13
1 2 a5 13
1 3 a1 8
1 3 a2 7
1 3 a3 9
1 3 a4 9
1 3 a5 9
2 1 a1 12
2 1 a2 12
2 1 a3 11
2 1 a4 17
2 1 a5 16
2 2 a1 13
2 2 a2 13
2 2 a3 13
2 2 a4 13
2 2 a5 12
2 3 a1 8
2 3 a2 7
2 3 a3 6
2 3 a4 10
2 3 a5 10
3 1 a1 13
3 1 a2 12
3 1 a3 11
3 1 a4 NA
3 1 a5 NA
3 2 a1 11
3 2 a2 11
3 2 a3 9
3 2 a4 NA
3 2 a5 NA
3 3 a1 9
3 3 a2 9
3 3 a3 9
3 3 a4 NA
3 3 a5 NA
", header = TRUE)

出错的分析代码 #

# 填充缺失值
data$Wealth[is.na(data$Wealth)] <- mean(data$Wealth, na.rm = TRUE)

# PERMANOVA分析
permanova_resultado <- adonis2(data$Wealth ~ Line + Installment + Year +
Line*Parcel + Line*Year + Year*Parcel, data = data,
permutations = 1000, strata = Line)

print(permanova_resultado)

后续尝试的代码 #

# 检验样地效应
riq_parc<-adonis2(data$Riqueza ~ Parcel, data = data, permutations = 1000, strata = data$Line)

print(riq_parc)

# 样地间成对t检验
Parcel <- factor(data$Parcel, labels = data$Parcel[1:3])

pairwise.t.test(data$Wealth, data$Parcel)

问题解答 #

1. 错误原因排查 #

• 变量名完全不匹配：定义的数据集是dados，但后续代码全使用data；且数据列名为Linha/Parcela/Ano/Riqueza，代码里却用Line/Parcel/Year/Wealth，导致程序无法正确读取分组变量，部分因子被错误识别为无有效水平。
• 分组变量未转为因子：Linha、Parcela、Ano默认是字符/数值类型，adonis2要求分组变量为因子类型，若未转换，可能被误判为只有1个水平，触发对比错误。
• 无效变量混入模型：公式中的Installment在数据中不存在，属于冗余项，会干扰模型对因子水平的识别。
• 缺失值填充失效：填充代码用了data$Wealth，但实际数据是dados$Riqueza，填充操作未执行，后续分析仍受缺失值影响。

2. 显著效应后的检验建议 #

• 不建议直接使用t检验：PERMANOVA是非参数检验，t检验是参数检验，二者假设前提差异大；且你的数据是非平衡设计，t检验无法控制样地嵌套在组内的结构误差，会增大假阳性风险。
• 替代方案：
  • 保持方法一致性：若PERMANOVA结果显著，使用pairwise.adonis（需额外安装）做组间成对PERMANOVA，检验具体组间差异。
  • 参数检验路径：若要使用参数方法，先验证数据是否满足正态性、方差齐性，再用混合效应模型（如lme4包的lmer函数）处理嵌套结构，非平衡设计的参数检验优先选择这类模型而非t检验。

内容的提问来源于stack exchange，提问作者Thielly Schmidt furtado