如何在R中绘制带Y轴断点的分组条形图（含散点与P值）

核心问题分析 #

你遇到的报错non-numeric argument to binary operator，直接原因是data1$Value是字符型，无法进行数值运算；另外gap.barplot默认处理单组数据，不支持直接生成分组条形图，因此X轴会显示单个数据点而非A/B分组。

完整解决方案 #

推荐用ggplot2结合ggbreak包实现需求——这套组合更灵活，能轻松添加条形图、抖动散点、Y轴断点和P值。

步骤1：修正数据类型并准备统计量

# 构造数据并将Value转为数值型
Condition <- c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B")
Value <- c("32", "43", "33", "65", "78", "66", "3", "1.2", "0.8", "2", "1", "1.4")
data1 <- data.frame(Condition, Value = as.numeric(Value))

# 计算分组均值和标准差（用于条形图和误差棒）
library(dplyr)
summary_data <- data1 %>%
  group_by(Condition) %>%
  summarise(mean_val = mean(Value),
            sd_val = sd(Value))

步骤2：绘制带断点的完整图形

library(ggplot2)
library(ggbreak)
library(ggsignif)

ggplot() +
  # 绘制分组条形图（基于均值）
  geom_bar(data = summary_data, aes(x = Condition, y = mean_val), 
           stat = "identity", fill = "#6baed6", width = 0.6) +
  # 添加标准差误差棒
  geom_errorbar(data = summary_data, aes(x = Condition, ymin = mean_val - sd_val, ymax = mean_val + sd_val),
                width = 0.2, color = "#2171b5") +
  # 添加抖动散点（避免数据重叠）
  geom_jitter(data = data1, aes(x = Condition, y = Value), 
              color = "#d73027", size = 2, alpha = 0.7, width = 0.2) +
  # 添加Y轴断点，跳过44-64区间
  scale_y_break(breaks = c(44, 64), space = 0.1) +
  # 设置Y轴范围
  ylim(0, 80) +
  # 自动计算并添加组间t检验P值
  geom_signif(comparisons = list(c("A", "B")), 
              map_signif_level = FALSE,
              test = "t.test",
              y_position = 75, size = 1) +
  # 美化图形
  labs(title = "分组带Y轴断点的统计图",
       y = "数值", x = "分组") +
  theme_bw() +
  theme(plot.title = element_text(hjust = 0.5))

关键细节说明 #

1. 数据类型修正：必须将Value转为数值型，否则所有数值相关的绘图、统计都会报错。
2. Y轴断点：ggbreak::scale_y_break的space参数控制断点处的空白宽度，可根据需求调整。
3. P值添加：如果需要用其他检验方法（如 Wilcoxon 秩和检验），只需修改geom_signif的test参数为"wilcox.test"；也可以用annotate手动指定P值：

   annotate("text", x = 1.5, y = 75, label = "P = 1.2e-05", size = 4)
   

4. 散点抖动：geom_jitter的width参数控制散点在X轴方向的抖动幅度，避免同组数据点重叠。

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