如何在R中将Welch ANOVA事后pairwise.t.test结果转为指定DataFrame

场景说明 #

执行Welch ANOVA后使用pairwise.t.test做事后检验，得到的结果是类似矩阵的结构，需要将其转换为包含group1、group2、P、Sig列的DataFrame，其中Sig列根据P值标注显著性（P>0.05为"ns"，P<0.05为"*"，P<0.01为"**"），且支持分组数大于3的情况。

实现步骤 #

先执行检验代码：

# Welch ANOVA分析
welch <- oneway.test(Petal.Width ~ Species, data = iris)
# 事后检验（Bonferroni校正）
ph <- pairwise.t.test(iris$Petal.Width, iris$Species, p.adjust.method = "bonferroni")

接下来处理结果矩阵：

1. 提取并转换矩阵为长格式数据框
   使用tidyr包的pivot_longer函数将宽格式的p值矩阵转为长格式，自动提取分组名称并过滤无效的NA值：

   library(tidyr)
   
   # 提取p值矩阵并转为数据框，添加行名为group1
   p_matrix <- as.data.frame(ph$p.value)
   p_matrix$group1 <- rownames(p_matrix)
   
   # 转为长格式，保留有效配对结果
   result_df <- pivot_longer(p_matrix,
                             cols = -group1,
                             names_to = "group2",
                             values_to = "P",
                             values_drop_na = TRUE)
   

2. 添加显著性标记列
   使用cut函数根据P值区间生成对应的显著性标记：

   result_df$Sig <- cut(result_df$P,
                        breaks = c(-Inf, 0.01, 0.05, Inf),
                        labels = c("**", "*", "ns"))
   

3. 格式化P值显示（可选）
   如果需要将极小的P值统一显示为0.000这类格式，用sprintf处理：

   result_df$P <- sprintf("%.3f", result_df$P)
   

最终结果 #

执行上述代码后，result_df就是符合要求的DataFrame，示例输出：

group1     group2     P Sig
1 versicolor  setosa   0.000 **
2 virginica   setosa   0.000 **
3 virginica versicolor 0.000 **

该方法支持任意数量的分组，分组名称从原结果矩阵的行/列名自动提取，无需手动指定分组数量。

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