如何绘制横向条形图并在X轴外侧添加p值？解决拼接布局问题

问题背景 #

我有如下R语言数据集：

microbial_data <- data.frame(
  Taxonomy = c("Nocardioidaceae", "Nocardioidaceae", "Nocardioidaceae",
              "Odoribacteraceae", "Odoribacteraceae", "Odoribacteraceae",
              "Rikenellaceae", "Rikenellaceae", "Rikenellaceae"),
  
  KW_p_adj = c(0.0177116, 0.0177116, 0.0177116,
             0.0177116, 0.0177116, 0.0177116,
             0.0177116, 0.0177116, 0.0177116),
  
  Group = c("C", "T1", "T2",
           "C", "T1", "T2",
           "C", "T1", "T2"),
  
  Abundance = c(0.000537303, 0.011049230, 0.002599831,
               0.000072700, 0.000391100, 0.002194485,
               0.000313588, 0.002816217, 0.013430997)
)

我想要绘制横向条形图，并在X轴外侧展示p值。之前用patchwork拼接两个ggplot图形的方法效果很差，会出现大量空白区域，图例位置也不合理。当前代码如下：

p <- ggplot(plot_data, aes(y = Taxonomy, x = Abundance, fill = Group)) + 
  geom_bar(stat = "identity", position = "dodge") +
  scale_fill_brewer(palette = "Set1") #+


p1 <- ggplot(plot_data, aes(y = Taxonomy, x = Abundance, fill = Group)) +
geom_text(aes(y = Taxonomy, x = 0), label = plot_data$KW_p_adj  ,
hjust = 0, fontface = "bold", inherit.aes = FALSE,)+xlim(c(0,1)) +
theme(panel.background = element_blank(),
        panel.grid = element_blank(),
        axis.line = element_blank(),
        axis.ticks = element_blank(),
        axis.text = element_blank(),
        axis.title = element_blank())
library(patchwork)
pall <- p + p1

优化方案 #

不需要用patchwork拼接，直接在同一个ggplot图层里添加p值文本，同时调整坐标轴范围来预留显示空间，即可解决空白区域和图例问题。

完整代码 #

library(ggplot2)

# 提取唯一的分类单元与对应p值，避免重复绘制文本
p_value_data <- unique(microbial_data[, c("Taxonomy", "KW_p_adj")])

# 计算X轴最大范围，预留p值显示空间
max_abundance <- max(microbial_data$Abundance)
x_max <- max_abundance * 1.2  # 扩展20%的空间用于放置p值

# 绘制图形
ggplot(microbial_data, aes(y = Taxonomy, x = Abundance, fill = Group)) +
  geom_bar(stat = "identity", position = position_dodge(width = 0.8)) +
  # 添加p值文本到X轴外侧
  geom_text(data = p_value_data, 
            aes(y = Taxonomy, x = x_max, label = paste("p =", round(KW_p_adj, 6))),
            hjust = 0, fontface = "bold", inherit.aes = FALSE) +
  scale_fill_brewer(palette = "Set1") +
  # 设置X轴范围，包含预留空间
  xlim(0, x_max) +
  labs(x = "Abundance", y = "Taxonomy") +
  theme_minimal() +
  theme(
    legend.position = "top",  # 调整图例到顶部，优化布局
    plot.margin = margin(10, 60, 10, 10)  # 右侧留边，防止p值被截断
  )

核心改进点 #

• 直接在主图层添加p值文本，消除拼接带来的空白区域
• 用unique处理p值数据，避免同一分类单元重复绘制文本
• 动态计算X轴范围，预留足够空间显示p值
• 调整图例位置和绘图边距，优化整体布局

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