Base R绘图中用expression()设置含上标与希腊字母的X轴标签报错求助

问题分析 #

你遇到的语法错误源于expression()函数内的字符串拼接不符合R的表达式规则：直接用"δ"^18" O"这种相邻字符串常量的写法会被R判定为语法错误，同时原代码里数据框的创建也存在问题（用<-而非=会导致列名异常）。

解决方案 #

要实现δ¹⁸O [‰] vs. PDB的X轴标签，需遵循expression()的语法规则，正确处理上标、特殊字符和文本拼接：

• 上标用^实现，比如δ^18（或用R内置的delta关键字）
• 用*连接表达式中的不同元素（比如上标后的符号和字母O）
• 用~添加空格，或直接用字符串包裹需要保留格式的文本

同时修正数据框创建错误：将data.frame(DO <- c(...))改为data.frame(DO = c(...))，确保列名正确。

修正后的完整代码 #

# 修正数据框创建方式：用=替代<-，避免列名异常
plankton_tow <- data.frame(DO = c(-1.2, -0.8, -0.5, -1.5), depth = c(10, 20, 30, 40))
core_tops <- data.frame(DO = c(-2.0, -1.6, -1.3, -2.5), depth = c(50, 60, 70, 80))

max_depth <- max(c(plankton_tow$depth, core_tops$depth))
plot(plankton_tow$DO, max_depth - plankton_tow$depth, 
     xlim = range(c(plankton_tow$DO, core_tops$DO)), 
     ylim = c(0, max_depth),
     # 正确的expression写法：*连接上标元素与文本，~添加空格
     xlab = expression(δ^18*O~"[‰] vs. PDB"), 
     ylab = "Depth [m BSL]", 
     pch = 16, col = "black", xaxt = "n", yaxt = "n")
points(core_tops$DO, max_depth - core_tops$depth, pch = 16, col = "red")

# 可选：添加坐标轴刻度（原代码xaxt="n"和yaxt="n"会隐藏刻度，如需显示可添加）
axis(1)
axis(2)

其他可选写法 #

如果更习惯用paste()拼接文本，也可以这样写：

xlab = expression(paste(δ^18, "O [‰] vs. PDB"))

或者用R内置的delta关键字替代Unicode的δ：

xlab = expression(delta^18*O~"[‰] vs. PDB")

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