ggplot2中如何将颜色标签格式化为带表达式的专业样式？

问题背景 #

需要用ggplot2绘制不同V_gs分组下value随vds变化的曲线，要求：

• 不同V_gs对应不同颜色和形状
• 图例标签采用U[GS] = 数值这类带表达式的格式（类似坐标轴标签的数学表达式样式）

尝试以下代码时出现类型错误：

scale_colour_discrete(
  labels = {function(x) {
      value <- gsub(".*\\.\\.(\+*)", "\\1", x)
      bquote('U[GS]' ~ ' = ' ~ .(value))}
    }
) +

报错信息：

Error in rep(value, length.out = nrows) : 
  attempt to replicate an object of type 'language'

后续尝试label_bquote()也出现同类的closure类型错误。

错误原因 #

1. 类型不匹配：bquote()返回单个语言对象，而ggplot的标签函数需要返回一个表达式向量，直接返回单个bquote结果会导致类型校验失败。
2. 正则表达式错误：原正则.*\\.\\.(\+*)针对的是用..分隔的字符串，但实际输入的列名是V_"gs"=3V这类格式，无法匹配出正确的数值。

可行解决方案 #

完整实现代码 #

修正正则表达式，并将标签函数改为返回表达式向量，同时补全原代码中缺失的坐标轴标签定义：

library(ggplot2)
library(reshape2)
library(dplyr)

# 读取数据
data <- read.csv("Input.txt", header=TRUE, sep="\t", quote = "")

# 转换为长格式
data_long <- melt(data, id.vars = "vds", variable.name = "V_gs", value.name = "value")

# 定义图例标签函数：提取数值并生成表达式向量
color_label <- function(x) {
  # 从列名中提取V_gs的数值（匹配"=数值V"部分）
  value <- gsub('.*=(.*)V', '\\1', x)
  # 生成表达式向量
  sapply(value, function(val) bquote(U[GS] ~ '=' ~ .(val) ~ 'V'))
}

# 自定义坐标轴标签函数
axis_label <- function(x) paste(x, "V")

ggplot(data_long, aes(
  x = vds,
  y = value,
  color = V_gs,
  shape = V_gs
)) +
  geom_line() +
  geom_point() +
  labs(
    x = expression(U[DS]),
    y = expression(I[DS]), # 修正原代码中y轴标签的笔误，改为电流更合理
    color = expression(U[GS]),
    shape = expression(U[GS])
  ) +
  scale_y_continuous(label = axis_label) +
  scale_x_continuous(label = axis_label) +
  guides(colour = guide_legend(reverse = TRUE), shape = guide_legend(reverse = TRUE)) +
  scale_colour_discrete(labels = color_label) +
  scale_shape_discrete(labels = color_label) +
  theme_minimal()

数据示例（Input.txt） #

vds V_"gs"=3V   V_"gs"=3.25V    V_"gs"=3.5V V_"gs"=3.75V    V_"gs"=4V   V_"gs"=4.25V    V_"gs"=4.5V V_"gs"=4.75V    V_"gs"=5V
0.000000000000000e+00   8.923415e-34    8.923415e-34    8.923415e-34    8.923415e-34    8.923415e-34    8.923415e-34    8.923415e-34    8.923415e-34
1.000000000000000e+00   6.629273e-02    8.460001e-02    1.005602e-01    1.145959e-01    1.270345e-01    1.381333e-01    1.480970e-01    1.570910e-01    1.652500e-01
2.000000000000000e+00   7.004236e-02    9.939607e-02    1.327373e-01    1.684299e-01    2.003164e-01    2.283812e-01    2.532684e-01    2.754860e-01    2.954397e-01
3.000000000000000e+00   7.004274e-02    9.940270e-02    1.327725e-01    1.697193e-01    2.100430e-01    2.534481e-01    2.992511e-01    3.413388e-01    3.786857e-01
4.000000000000000e+00   7.004237e-02    9.939607e-02    1.327373e-01    1.697194e-01    2.100469e-01    2.534753e-01    2.998012e-01    3.486163e-01    4.001946e-01
5.000000000000000e+00   7.004237e-02    9.939608e-02    1.327373e-01    1.697193e-01    2.100430e-01    2.534481e-01    2.997063e-01    3.486189e-01    4.000209e-01

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