如何用ggplot制作Y轴为计数的基因测序数据堆叠柱状图

嘿，我来帮你搞定这个ggplot堆叠柱状图的问题！你的核心问题是没有让ggplot正确统计每个分组的数量，反而把类别变量直接放到了Y轴上，所以才会显示tRF.type的类别而不是计数。下面给你两种可行的解决方法：

问题根源 #

你原来的代码里把y = tRF.type放进了aes()映射中，这会让ggplot把Y轴识别为类别变量。而堆叠柱状图需要的是每个氨基酸分组下，不同tRF.type的出现次数，所以需要调整映射参数或者提前汇总数据。

方法1：让ggplot自动统计计数（最简单） #

geom_bar()默认的stat参数是"count"，它会自动帮你统计每个分组的数量，你只需要去掉Y轴的错误映射，把填充色绑定到tRF.type即可：

# 先加载你的示例数据
tRF <- c('tRF-26-OB1690PQR3E', 'tRF-27-OB1690PQR3P', 'tRF-30-MIF91SS2P46I')
tRF.type <- c('5-tRF', 'i-tRF', '3-tRF')
Amino.Acid <- c('Ser', 'Lys', 'Ser')
tRF.data <- data.frame(tRF, tRF.type, Amino.Acid)

library(ggplot2)

ggplot(tRF.data, aes(x = Amino.Acid, fill = tRF.type)) +
  geom_bar(position = "stack") +  # 默认stat="count"，自动计数
  ggtitle("LAN5 - 4 days post CNTF treatment") +
  xlab("Amino Acid") +
  ylab("Count of tRF Type")  # 把Y轴标签改成计数相关的描述

方法2：手动汇总数据（更灵活） #

如果之后需要对数据做更多自定义处理（比如过滤、计算比例），可以先用dplyr提前统计好每个组的计数，再传入ggplot：

library(dplyr)
library(ggplot2)

# 按氨基酸和tRF.type分组，统计每组的数量
summary_data <- tRF.data %>%
  group_by(Amino.Acid, tRF.type) %>%
  summarise(count = n(), .groups = "drop")  # .groups="drop"避免分组残留

# 用汇总后的数据画图
ggplot(summary_data, aes(x = Amino.Acid, y = count, fill = tRF.type)) +
  geom_bar(stat = "identity", position = "stack") +
  ggtitle("LAN5 - 4 days post CNTF treatment") +
  xlab("Amino Acid") +
  ylab("Count of tRF Type")

效果说明 #

修改后，Y轴会显示每个氨基酸类别下对应tRF.type的出现次数，堆叠的柱子会按tRF.type的颜色区分，完全符合你想要的可视化需求。

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