在ggplot2中为数据及轴标签添加上标与负号的方案求助

完整代码 #

library(ggplot2)

# 原始数据
dt <- data.frame(outcome = c("A", "B", "C", "D"), 
                 beta = c(-0.007, -0.008, -0.009, -0.009),
                 LCI = c(-0.005, -0.006, -0.007, -0.008), 
                 UCI = c(-0.008, -0.009, -0.01, -0.01), 
                 P = c(0.000001, 0.00002, 0.000003, 0.00004))

# 自定义函数：格式化p值为"系数 × 10^指数"的上标形式
format_p <- function(p_val) {
  # 转成科学计数字符串
  sci_str <- format(p_val, scientific = TRUE)
  # 拆分系数和指数
  parts <- strsplit(sci_str, "e")[[1]]
  coeff <- as.numeric(parts[1])
  exp <- as.integer(parts[2])
  # 生成带解析标记的字符串，用于ggplot的parse=TRUE
  paste0(coeff, " %*% 10^", exp)
}

# 自定义函数：替换短横线为Unicode负号
format_neg <- function(num) {
  gsub("-", "\u2212", format(num, nsmall = 3))
}

# 绘图
pt <- ggplot(data = dt, aes(
  y = reorder(outcome, beta),
  x = beta,
  xmin = LCI,
  xmax = UCI
)) +
  geom_point(color = "black", shape = 18, size = 5) +
  geom_errorbarh(height = 0.1, linewidth = 1, color = "black") +
  geom_text(
    aes(
      label = paste0(
        format_neg(beta), " (", format_neg(LCI), ";", format_neg(UCI), "), p = ", format_p(P)
      ),
      x = beta,
      y = outcome
    ),
    size = 4,
    color = "black",
    vjust = -1.5,
    show.legend = FALSE,
    parse = TRUE  # 启用解析，让p值的上标生效
  ) +
  labs(title = 'My plot', x = 'Beta', y = 'Outcome') +
  theme_classic() +
  theme(
    plot.title = element_text(hjust = 0.05, size = 10, face = "bold"),
    axis.text = element_text(size = 10, face = "bold"),
    axis.title = element_text(size = 10, face = "bold")
  ) +
  # 修正x轴刻度的负号
  scale_x_continuous(
    breaks = c(-0.01, -0.009, -0.008, -0.007, -0.006, -0.005),
    labels = function(x) gsub("-", "\u2212", format(x, nsmall = 3))
  )

print(pt)

关键问题解决说明 #

• 负号替换：

  1. 你之前的gsub语法错误，正确写法是gsub("-", "\u2212", 目标字符串)，参数顺序为pattern, replacement, x。
  2. 定义format_neg函数，先将数值格式化为固定小数位的字符串，再替换短横线为Unicode负号（\u2212），确保beta、LCI、UCI的标签和x轴刻度都使用美观的负号。
  3. x轴刻度通过匿名函数直接替换刻度文本中的短横线。

• p值格式化：

  1. 定义format_p函数，先将p值转为科学计数字符串，拆分出系数和指数，生成符合ggplot解析语法的字符串（paste0(coeff, " %*% 10^", exp)）。
  2. 在geom_text中设置parse = TRUE，让ggplot自动解析字符串中的数学表达式，把指数显示为上标，最终得到类似1 × 10^-6的美观格式。

额外说明 #

• 如果需要更精细的系数格式化（比如保留1位小数），可以修改format_p函数中的coeff部分，例如coeff <- sprintf("%.1f", as.numeric(parts[1]))。
• 若遇到标签重叠问题，可以调整vjust或hjust参数，或者使用ggrepel包的geom_text_repel替代geom_text。

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