如何在ggplot自动化标注中使用变量值作为下标文本

问题说明 #

我需要在geom_density_ridges()绘制的密度脊图中添加标注，用来展示scale_x_continuous()定义的x轴限值。最终要实现的效果是用变量limit_quantile的值作为标注文本的下标（参考示例中grob1的标注样式），因为后续要通过apply族函数批量绘图，不同变量需要对应不同参数，必须实现自动化标注。

尝试过将变量转字符传入bquote()、使用get()函数，但都没成功，示例代码如下：

## 尝试用变量值生成下标文本
# 2022-DEC-02

df1 <- data.frame(v1 = c(1:100), class = 1, v2 = rnorm(100) )
df2 <- data.frame(v1 = c(1:100), class = 2, v2 = rnorm(100)*5+10 )
df3 <- data.frame(v1 = c(1:100), class = 3, v2 = rnorm(100)*10+25 )

df <- rbind(df1,df2,df3)


limit_quantile <- 0.99

ridge_plot <- ggplot(df, aes(x = v2, y = as.factor(class), group = as.factor(class))) + 
    geom_density_ridges(stat = "binline", rel_min_height = 0.005 , binwidth = 1,  alpha = 0.5,scale = 0.99) +
    scale_y_discrete(expand = c(0.01, 0), breaks = c(1:3)) +
    scale_x_continuous(expand = c(0.01, 0), limits = quantile(df$v2, c(0, limit_quantile))) +
    theme_ridges() +
    xlab("Values of v2") +
    ylab(NULL) +
    theme(axis.title.x = element_text(hjust = 0.5))


# 手动指定下标0.95的效果（目标样式）
grob1 <- grobTree(textGrob(bquote("X-Axis-Limit: Q"[0.95]), x=0.9,  y=0.1, hjust=1,
                           gp=gpar(fontsize = 10.5, fontface="italic")))

ridge_plot_text1 <- ridge_plot  + annotation_custom(grob1)
ridge_plot_text1

# 尝试用变量limit_quantile做下标（失败）
grob2 <- grobTree(textGrob(bquote("X-Axis-Limit: Q"[limit_quantile]), x=0.9,  y=0.1, hjust=1,
                           gp=gpar(fontsize = 10.5, fontface="italic")))
ridge_plot_text2 <- ridge_plot  + annotation_custom(grob2)
ridge_plot_text2

# 其他尝试
limit_quantile_text <- as.character(limit_quantile)
grobTree(textGrob(bquote("X-Axis-Limit: Q"[limit_quantile_text ]), x=0.9,  y=0.1, hjust=1,
                           gp=gpar(fontsize = 10.5, fontface="italic")))

# 第二次尝试
grob2 <- grobTree(textGrob(bquote("X-Axis-Limit: Q"[get(limit_quantile)]), x=0.9,  y=0.1, hjust=1,
                           gp=gpar(fontsize = 10.5, fontface="italic")))

解决方案 #

方法1：使用bquote()的.()语法嵌入变量 #

bquote()中可以用.()来引用当前环境中的变量值，这是最直接的解决方式：

limit_quantile <- 0.99

# 正确生成带变量下标的文本grob
grob3 <- grobTree(textGrob(bquote("X-Axis-Limit: Q"[.(limit_quantile)]), 
                           x=0.9, y=0.1, hjust=1,
                           gp=gpar(fontsize = 10.5, fontface="italic")))

ridge_plot_text3 <- ridge_plot + annotation_custom(grob3)
ridge_plot_text3

方法2：使用glue结合parse=TRUE（更直观的字符串拼接） #

如果觉得bquote的语法不够直观，可以用glue包生成带LaTeX式下标语法的字符串，再通过parse=TRUE解析：

library(glue)

annotation_text <- glue("X-Axis-Limit: Q[{limit_quantile}]")

# 用annotate直接添加，无需grobTree
ridge_plot + annotate("text", x = Inf, y = -Inf, 
                      label = annotation_text, 
                      hjust = 1.1, vjust = -0.5,
                      fontface = "italic", size = 3.5,
                      parse = TRUE)

方法3：直接提取x轴限值并标注（更精准） #

如果需要标注的是实际的x轴截断值（而非分位数本身），可以直接从绘图对象中提取限值，再生成标注：

# 提取x轴的上限值
x_limit <- layer_scales(ridge_plot)$x$limits[2]

# 生成包含实际限值的标注
ridge_plot + annotate("text", x = x_limit, y = -Inf,
                      label = glue("X-Axis Limit: {round(x_limit,2)}"),
                      hjust = 1.1, vjust = -0.5,
                      fontface = "italic", size = 3.5)

关键说明 #

• bquote(.())：用于在表达式中插入变量的实际值，是R中生成动态数学表达式的常用方法
• parse=TRUE：告诉ggplot将字符串解析为数学表达式，支持下标（[]）、上标（^）等LaTeX语法
• 批量绘图时，只需将limit_quantile作为参数传入绘图函数，上述方法都能自动生成对应标注

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