如何用ggplot绘制物种配对p值的三角相关矩阵风格图？

1. 锁死物种顺序（核心） #

把长格式数据里的物种列转成因子，指定levels为原数据里的物种顺序——这是让图中物种不自动排序的关键：

# 假设长格式数据叫df_long，物种列是species1和species2
# 提取原数据里的物种顺序（这里取species1的唯一值顺序，按需调整）
species_order <- unique(df_long$species1)
# 转成因子，锁定顺序
df_long$species1 <- factor(df_long$species1, levels = species_order)
df_long$species2 <- factor(df_long$species2, levels = species_order)

2. 给p值做二元分组 #

新增一列标记显著性，直接区分两种情况：

df_long$sig_status <- ifelse(df_long$p_value <= 0.05, "显著 (p≤0.05)", "不显著 (p>0.05)")

3. 实现三角布局 #

用索引过滤的方式只保留三角区域的数据，比如要下三角（含对角线）就这么做：

# 给每个物种分配索引
df_long$idx1 <- as.integer(df_long$species1)
df_long$idx2 <- as.integer(df_long$species2)
# 过滤出下三角数据（要上三角就改成idx2 <= idx1）
df_tri <- df_long[df_long$idx2 >= df_long$idx1, ]

4. 绘制二元颜色矩阵图 #

用geom_tile搭配手动颜色映射，避免渐变效果：

library(ggplot2)

ggplot(df_tri, aes(x = species1, y = species2, fill = sig_status)) +
  geom_tile(color = "white", size = 0.5) + # 加白色边框区分单元格
  # 自定义两种颜色，按需替换
  scale_fill_manual(values = c("显著 (p≤0.05)" = "#dc322f", "不显著 (p>0.05)" = "#268bd2")) +
  # 旋转x轴标签防止重叠
  theme(axis.text.x = element_text(angle = 45, hjust = 1)) +
  # 强制用我们锁死的物种顺序
  scale_x_discrete(limits = species_order) +
  scale_y_discrete(limits = species_order) +
  # 去掉多余的轴标签，修改图例标题
  labs(x = "", y = "", fill = "显著性") +
  theme_minimal() +
  # 去掉网格线，让图更清爽
  theme(panel.grid = element_blank())

如果你是从宽格式数据开始 #

先把宽格式转成长格式，再按上面的步骤操作：

library(tidyr)

# 假设宽格式数据叫df_wide，第一列是物种名，其他列是配对物种的p值
df_long <- df_wide %>%
  pivot_longer(cols = -1, names_to = "species2", values_to = "p_value") %>%
  rename(species1 = 1)

内容的提问来源于stack exchange，提问作者N.F.F.