geom_signif：对数y轴上绘制p值标签的问题排查

云原生阿岚

2026-4-29

解决对数Y轴下ggsignif添加p值的警告与显示问题

嘿，我看了你的代码和遇到的警告，问题根源其实很明确：你的数据里有大量0值，对数转换后会变成负无穷，导致这些行被自动过滤，同时你设置的p值标签位置和对数轴的刻度完全不匹配，所以才会出现一堆警告，标签也显示不出来。我给你整理了修正方案，一步步来：

首先分析警告原因

你用了scale_y_log10()，但log10(0)是负无穷，这就触发了“Transformation introduced infinite values”的警告，而且ggplot会自动移除这些非有限值的行，也就是那29个0值的样本。
你设置的y_pos是3.5、3.25这种正数，但你的数据最大值才约0.0004，对数轴上0.0004对应log10(-3.397)，你的y_pos完全不在轴的范围内，所以stat_signif找不到有效位置，只能移除这些标签。

修正方案与代码

1. 处理0值：添加极小偏移量

对数转换无法处理0，我们可以给所有0值加一个极小的常数（比如数据里最小非0值的1/10），这样既不会改变数据的相对分布，又能避免无穷值：

library(magrittr)
library(ggsignif)
library(ggpubr)

# 原始数据
value<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.000273200055869626, 0, 0, 0, 0, 0, 0.000391899951754252, 7.79999691577276e-06, 0, 0, 0.000341800044604173, 0, 0.000160399989889521, 0, 0, 0.00014739991923847, 0, 0, 0.000200399960107021, 0, 0, 0, 0, 6.35000054462966e-05)
condition<-c("E","E","E","E","E","E","E","E","E","E","E","E","E", "R","R","R","R","R","R","R","R","R","R","R","R","R", "S","S","S","S","S","S","S","S","S","S","S")
species<-rep('B_theta',length(condition))
df <- data.frame(value,condition,species)

# 给0值添加极小偏移量，避免对数转换出错
min_non_zero <- min(df$value[df$value > 0])
df$value_adjusted <- ifelse(df$value == 0, min_non_zero/10, df$value)

2. 重新计算统计量并设置正确的标签位置

这次用调整后的数据计算p值，同时基于数据实际范围设置p值标签的y轴位置：

# 计算组间差异检验
bact <- compare_means(value_adjusted~condition, data = df, method = 't.test') %>% 
  mutate(
    # 取调整后数据的最大值，按倍数设置标签位置（确保在图形上方）
    max_val = max(df$value_adjusted),
    y_pos = c(max_val * 2, max_val * 1.7, max_val * 1.4),
    # 按规则设置标签内容
    labels = ifelse(p < 0.1, p.format, p.signif)
  )

3. 绘制图形并修正geom_signif参数

注意x轴是因子变量，xmin和xmax需要用因子的索引位置，而不是字符：

p <- ggboxplot(df, x = "condition", y = "value_adjusted", color = "condition", palette = "jco", add = "jitter")+
  geom_signif(
    textsize=6, 
    tip_length = 0.01, 
    # 用match获取因子的位置索引
    xmin = match(bact$group1, unique(df$condition)),
    xmax = match(bact$group2, unique(df$condition)),
    annotations = bact$labels, 
    y_position = bact$y_pos
  )+
  scale_y_log10(
    # 调整y轴范围，包含所有数据和标签位置
    limits = c(min(df$value_adjusted)/2, max(bact$y_pos)*1.1),
    # 可以添加刻度标签，让对数轴更清晰
    breaks = c(min_non_zero/10, min_non_zero, 1e-4, 3e-4)
  )+
  annotation_logticks(sides = "l")+
  labs(y = "Value (log10 scale)")

print(p)

为什么这样修改有效？

偏移0值后，所有数据都能正常进行对数转换，不会产生无穷值，警告自然消失。
标签位置基于数据实际最大值设置，确保在对数轴的可见范围内，stat_signif能正确渲染标签。
用match()获取x轴因子的索引，解决了字符型x轴位置匹配错误的问题。

运行这段代码后，你就能得到一个没有警告、p值标签正确显示的对数轴箱线图啦！

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