基于Group变量两类别的CTaa_beta共享情况冲积图绘制求助

ggplot2原生对复杂冲积图的支持有限，推荐使用专门的ggalluvial包来实现你的需求，以下是具体解决步骤：

1. 安装并加载必要工具包 #

install.packages("ggalluvial")
library(ggalluvial)
library(ggplot2)

2. 确认数据集结构 #

确保你的数据集a包含以下核心列：

• Group：取值为InEx和ReA SFMC的分类变量
• CTaa_beta：作为填充色的分组变量
• n：数值型的频率变量，用于控制冲积带厚度

3. 绘制标准冲积图 #

使用ggalluvial的专属图层实现符合预期的效果：

ggplot(a,
       aes(y = n, axis1 = Group, axis2 = CTaa_beta)) +
  # 绘制连接分组的冲积带，厚度由n控制
  geom_alluvium(aes(fill = CTaa_beta), width = 1/12) +
  # 绘制分组的矩形层，增强可读性
  geom_stratum(width = 1/12, fill = "white", color = "black") +
  # 为分组添加标签
  geom_text(stat = "stratum", aes(label = after_stat(stratum))) +
  # 设置X轴参数
  scale_x_discrete(limits = c("Group", "CTaa_beta"), expand = c(0.05, 0.05)) +
  # 自定义配色（可替换为你需要的方案）
  scale_fill_viridis_d(option = "plasma") +
  # 添加标题与轴标签
  labs(title = "Group与CTaa_beta的冲积频率图",
       y = "频率(n)",
       fill = "CTaa_beta") +
  theme_minimal()

4. 突出两类间共享的CTaa_beta #

如果需要仅展示两个Group中都存在的CTaa_beta，可先筛选数据再绘图：

# 筛选出在两个Group中都出现的CTaa_beta值
shared_cta <- unique(a$CTaa_beta[duplicated(a$CTaa_beta)])
a_shared <- subset(a, CTaa_beta %in% shared_cta)

# 绘制共享类别的冲积图
ggplot(a_shared,
       aes(y = n, axis1 = Group, axis2 = CTaa_beta)) +
  geom_alluvium(aes(fill = CTaa_beta), width = 1/12) +
  geom_stratum(width = 1/12, fill = "white", color = "black") +
  geom_text(stat = "stratum", aes(label = after_stat(stratum))) +
  scale_x_discrete(limits = c("Group", "CTaa_beta"), expand = c(0.05, 0.05)) +
  labs(title = "Group间共享CTaa_beta的冲积频率图",
       y = "频率(n)",
       fill = "CTaa_beta") +
  theme_minimal()

常见问题排查 #

• 若冲积带未正确连接：检查数据是否包含每个CTaa_beta对应两个Group的完整记录
• 若厚度显示异常：确认n列是数值型变量，无缺失或非数值内容
• 避免使用ggplot2原生的geom_flow：该图层在多类别连接场景下易出现逻辑错误，ggalluvial是更专业的替代工具

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