R语言绘图中的数学公式标注：从基础语法到ggplot2与latex2exp实战

R语言数据可视化数学公式

于 2026-08-04 07:12:21 修改

·本内容遵循CC 4.0 BY-SA版权协议

1. 项目概述：为什么R绘图中的公式标注是个技术活？

如果你用R语言做过科研图表，尤其是需要在坐标轴标签、图例或者图形内部标注一些数学符号、上下标或者复杂公式时，大概率会和我一样，经历过一段“抓狂”的时期。明明在LaTeX或者Word里信手拈来的表达式，一到R的plot()函数里就变得支离破碎，要么是上标显示不出来，要么是希腊字母变成了乱码，更别提复杂的分式或者求和符号了。这不仅仅是美观问题，它直接关系到图表信息的准确传达和论文、报告的专业性。

这个项目的核心，就是彻底解决在R基础绘图系统（base R graphics）以及ggplot2等高级绘图包中，如何正确、优雅地输出包含上下标、数学公式和特殊符号的文本。这绝不仅仅是调用一个expression()函数那么简单，它涉及到R图形设备对数学注释（mathematical annotation）的解析规则、不同函数（如text(), mtext(), axis()）的细微差别，以及在ggplot2语法体系下的特殊实现方式。掌握这些技巧，意味着你能让数据可视化成果在学术和专业场景下，真正达到出版级的水准。

2. 核心原理：R如何理解你的数学表达式？

在开始实操之前，我们必须先理解背后的原理。R本身并不直接渲染数学公式，它依赖的是底层图形设备（如屏幕的windows()、quartz()，或用于输出的pdf()、png()）对一种特殊标记语言的处理能力。这种标记语言，我们可以粗略地理解为R版本的“简易LaTeX”。

2.1 `expression()`、`bquote()` 与 `~` 的作用域

最核心的函数是expression()。它的作用是将一段文本“标记”为需要被数学解析的表达式对象，而不是普通的字符串。

# 错误示范：直接使用字符串

plot(1, main = “y = x^2”) # “^2”会被当作普通字符显示，不会变成上标

# 正确示范：使用expression()

plot(1, main = expression(y == x^2)) # “^”被正确解析为上标运算符

这里有几个关键点：

== 表示等号：在expression()内部，单个等号=通常有赋值或其他含义，因此显示数学等号需要使用双等号==。
^ 表示上标：这是固定的语法。
expression() 返回一个对象：main = expression(...) 是将这个表达式对象传递给main参数，图形设备在绘制标题时，会识别并渲染这个对象。

但expression()有一个局限：它内部的文本是“写死”的，很难方便地插入变量值。这时就需要bquote()出场。bquote()实现了“部分引用”，允许你在表达式中用 .() 来包裹需要求值的R变量。

x_power <- 3

plot(1, main = bquote(y == x^.(x_power))) # 动态插入变量x_power的值3

另一个重要的符号是波浪号 ~，它在表达式中代表一个空格。在数学排版中，空格对于可读性至关重要，但直接打空格可能会被忽略。~ 能插入一个适当宽度的空格。

plot(1, main = expression(Alpha~”vs.”~Beta)) # 在希腊字母和单词间插入空格

2.2 常用数学符号的语法速查

R的数学注释支持一套相对固定的符号集，以下是最常用的部分：

上下标：x^2 表示x的平方（上标），x[1] 表示x下标1。
希腊字母：直接拼写其英文名称，如 alpha, beta, gamma, Omega, Sigma。首字母大写表示大写希腊字母。
数学运算符：
- 加减乘除：+, -, *, /
- 乘号：%*%
- 无穷大：infinity
- 部分微分：partial
- 正负号：+-
- 角度符号：degree
括号与分隔符：圆括号()可直接使用。对于可伸缩的括号（如用于矩阵或分段函数），使用 group(“(“, list, “)”) 结构，但更复杂的排版建议依赖latex2exp包。
分式与根式：frac(a, b) 表示分数a/b。sqrt(x) 表示平方根，sqrt[n](x) 表示n次根号下x。
修饰符：bold() 加粗，italic() 斜体，plain() 正常字体。

注意：expression() 中的语法是大小写敏感且对空格极其敏感的。alpha 和 Alpha 是天壤之别。符号名称必须完全正确，拼写错误会导致其被当作普通文本显示。

3. 实战演练：在图形的各个位置添加公式

理解了基本语法后，我们将其应用到绘图的具体环节。我将以基础绘图系统为例，因为其原理完全适用于ggplot2。

3.1 主标题、坐标轴标签与图例

这是最常用的场景。你需要修改main（主标题）、xlab/ylab（坐标轴标签）、sub（副标题）等参数。

# 创建一个示例图形框架

plot(1:10, 1:10, type = “n”,

main = expression(“Relationship between” ~ alpha ~ “and” ~ beta^2),

xlab = expression(“Input Variable” ~ (mu*mol~m^{-2}~s^{-1})),

ylab = expression(“Response” ~ (frac(Delta*Yield, “%”)))

)

这段代码实现了：

主标题：将单词和希腊字母、上标混合排版。
X轴标签：包含单位，其中 mu 渲染为微符号μ，m^{-2} 表示平方米的负二次方，s^{-1} 表示每秒。~ 确保了单位各部分间的空格。
Y轴标签：使用 frac 创建了一个分式，Delta 渲染为增量符号Δ，”%” 用引号包裹表示普通文本百分号。

为图例添加公式稍微特殊一些，因为legend()函数的legend参数可以直接接受表达式向量。

x <- 1:10

y1 <- x^2

y2 <- log(x)

plot(x, y1, type = “l”, col = “blue”, ylim = c(0, max(y1, y2)))

lines(x, y2, col = “red”)

legend(“topleft”,

legend = c(expression(y == x^2),

expression(y == log(x))),

col = c(“blue”, “red”), lty = 1)

3.2 在图形内部添加文本：`text()` 与 `mtext()`

当需要在数据区域内部特定坐标位置添加注释时，使用text()。

plot(1:10, (1:10)^2)

# 在点(5, 25)附近添加一个带公式的注释

text(5, 25,

expression(hat(y) == beta[0] + beta[1]*x),

cex = 1.2, col = “darkred”)

# 在点(8, 64)附近添加另一个注释

text(8, 64,

bquote(R^2 == .(round(cor(1:10, (1:10)^2)^2, 3))),

adj = 0)

这里展示了两个技巧：

hat(y) 渲染为y的估计值（y上加帽）。
在text()中使用bquote()动态插入计算得到的R平方值，.()内的cor()函数会先被计算，结果再嵌入表达式。

mtext()用于在图形区域的**边空（margin）**添加文本，常用于添加统一的单位说明或脚注。

plot(rnorm(100))

mtext(side = 1, # 底部边空

expression(“Data:” ~ N(mu == 0, sigma^2 == 1)),

line = 3, cex = 0.8)

3.3 自定义坐标轴刻度标签

有时，默认的数值刻度需要被替换为包含符号的标签，比如你想把X轴的1,2,3显示为alpha, beta, gamma。

plot(1:3, 1:3, xaxt = “n”, xlab = “Phase”) # 先取消默认X轴

axis(side = 1, at = 1:3,

labels = expression(alpha, beta, gamma))

axis()函数的labels参数可以接受一个表达式向量，每个表达式对应一个刻度位置（at参数）。

4. 进阶应用：在ggplot2中驾驭公式与使用latex2exp

ggplot2的语法哲学与基础绘图不同，但它处理数学表达式的方式一脉相承，主要使用labeller函数和plotmath语法（即前述的expression语法）。

4.1 ggplot2中的公式标注

在ggplot2中，为图形属性添加标签通常使用labs()函数，而数学表达式需要借助expression()或bquote()。

library(ggplot2)

df <- data.frame(x = 1:10, y = (1:10)^2)

ggplot(df, aes(x, y)) +

geom_point() +

labs(

title = expression(“Quadratic Relationship:” ~ y == x^2),

x = expression(“Independent Variable” ~ (xi)),

y = expression(“Dependent Variable” ~ (frac(mm^2, s)))

) +

annotate(“text”, x = 5, y = 80,

label = expression(hat(y) == 1.2*x^2 - 0.5),

parse = TRUE)

这里有一个至关重要的参数：parse = TRUE。 在ggplot2中，除了labs()等直接接受表达式的函数外，在annotate()或geom_text()中使用字符串标签时，必须设置parse = TRUE，告诉ggplot2将此字符串解析为数学表达式。如果你直接传递一个expression()对象，则不需要parse=TRUE。

4.2 突破局限：`latex2exp`包的力量

尽管R原生的plotmath语法功能强大，但对于用惯LaTeX的用户来说，它的学习成本高，且表达能力有限，无法处理非常复杂的公式（如多行公式、大型矩阵、特定符号等）。

这时，latex2exp包就是救星。它允许你直接在R图形中使用LaTeX语法！

# 安装并加载包

# install.packages(“latex2exp”)

library(latex2exp)

plot(1, 1,

main = LatexExpr(r”($\int_{-\infty}^{\infty} e^{-x^2} dx = \sqrt{\pi}$)”),

xlab = LatexExpr(r”($\frac{\partial u}{\partial t} = \alpha \nabla^2 u$)”),

cex.main = 1.5)

使用技巧与注意事项：

原始字符串 r”()：在R 4.0及以上版本，使用原始字符串可以避免对LaTeX中的反斜杠\进行转义，让语法更清晰。低版本R需使用双反斜杠\\。
LatexExpr() 函数：这是latex2exp包的核心函数，它将LaTeX字符串转换为R图形设备可以理解的表达式对象。
适用范围：它同样适用于ggplot2。

R

ggplot(df, aes(x, y)) +

geom_point() +

labs(title = LatexExpr(r”($y = \beta_0 + \sum_{i=1}^n \beta_i x_i + \epsilon$)”))
性能：复杂公式的首次渲染可能会稍慢，因为包需要调用LaTeX引擎进行编译和转换。但对于最终输出的出版级图表，这点时间成本完全可以接受。

实操心得：对于95%的日常图表，R原生的expression()和bquote()语法完全够用，且与图形系统集成度更高。只有当遇到极其复杂的公式，或者你本人就是LaTeX高手，希望复用现有代码时，才优先考虑latex2exp。混合使用两者时要注意，不要在同一个标签字符串里混用两种语法。

5. 常见问题排查与字体陷阱

即使语法正确，你可能还是会遇到显示问题。以下是一些常见坑点及其解决方案。

5.1 公式显示为乱码或方框

这几乎总是字体问题。R的图形设备在渲染数学符号时，依赖于特定的符号字体（如Computer Modern数学字体）。在某些设备（特别是Windows的默认windows()屏幕设备）或输出格式（如png()）中，如果找不到这些字体，就会显示乱码。

解决方案：

使用PDF或PostScript输出：这是最可靠的方法。pdf()和postscript()设备内置了对数学字体的完美支持。

R

pdf(“my_plot.pdf”, width=8, height=6)

# 绘制你的带公式的图

dev.off()
在png()/jpeg()中指定字体族：如果必须输出为位图，可以尝试在打开设备时指定一个包含数学符号的字体族。

R

png(“my_plot.png”, width=800, height=600, family=“serif”) # 或 “Times”

# 绘制图形

dev.off()

你可以用names(pdfFonts())查看系统可用的PDF字体，但并非所有都支持数学符号。”Times”, ”Helvetica”, ”Courier” 等标准字体族通常支持较好。
在ggplot2中全局设置字体：

R

library(showtext) # 一个强大的字体管理包

# 可以加载系统字体或谷歌字体

theme_set(theme_gray(base_family = “Times New Roman”))

5.2 上下标位置不对或大小不协调

这通常是由于错误地混合了expression()模式和普通文本模式。

错误：expression(paste(“R^2 = “, 0.95))。paste()在expression()内部是无效的，它会被当作普通文本。
正确：使用 ~ 或 * 连接。* 表示无空格连接，~ 表示有空格连接。

R

expression(R^2 == 0.95) # 等号连接

expression(“Model Fit:” ~ R^2 == 0.95) # 空格连接文本和公式

expression(alpha*”%”) # 无空格连接希腊字母和百分号

5.3 如何在公式中换行？

这是一个常见需求，比如一个很长的公式需要分成两行显示。在基础绘图中，可以使用atop()或paste()在表达式中模拟换行，但效果不完美。

plot(1, main = expression(atop(“First Line with” ~ alpha,

“Second Line with” ~ beta^2)))

在ggplot2中，更优雅的方式是使用\n换行符，并结合parse = TRUE（但需要注意，换行符本身可能会破坏表达式解析，通常需要将多行拆分为多个独立的图形元素，如使用多个annotate）。

更稳健的解决方案是直接使用latex2exp，因为LaTeX本身对多行公式（如aligned环境）有完美的支持。

library(latex2exp)

latex_str <- r”($\begin{aligned}

y &= \beta_0 + \beta_1 x_1 + \\

&\quad \beta_2 x_2 + \epsilon

\end{aligned}$)”

plot(1, main = LatexExpr(latex_str))

5.4 表达式与因子（factor）或字符串拼接的冲突

当你需要动态生成大量包含变量的标签时，可能会想用paste()或sprintf()先构造字符串，再传给expression()，这是行不通的。

# 错误

r2_val <- 0.956

my_label <- paste(“R^2 =“, r2_val)

text(5, 5, expression(my_label)) # 只会显示字符串“my_label”

# 正确：使用bquote()

text(5, 5, bquote(R^2 == .(r2_val)))

# 正确：循环生成表达式列表

coef_names <- c(“alpha”, “beta”, “gamma”)

label_list <- lapply(coef_names, function(name) bquote(.(name)^2))

# label_list 就是一个表达式列表，可以用于legend()的labels参数

掌握R中公式绘制的技巧，是从“能画图”到“能画出专业图”的关键一步。它没有捷径，核心就是理解expression()、bquote()的语法规则，并在实践中反复调试。我的经验是，对于任何复杂的标签，先在图形窗口用简单的plot(1, main=expression(...))测试显示效果，确认无误后再整合到最终的绘图代码中。当原生语法力有不逮时，不要犹豫，请出latex2exp这个终极武器。最后，永远优先选择PDF作为最终输出格式，它能最大程度保证数学排版的保真度。