如何在ggplot2分面标签中显示希腊字母phi=数值与英文字符sd=数值
问题说明
你需要修改ggplot2分面网格图的标签,具体要求如下:
- 列标签显示为
φ = 0.8、φ = 0.9、φ = 0.95,其中phi要渲染为希腊字母符号 - 行标签显示为
sd = 1、sd = 3、sd = 5、sd = 10,sd保留原英文字母格式
现有模拟ARIMA(1,0,0)序列并绘图的代码如下:
## 模拟ARIMA(1, 0, 0)序列 set.seed(289805) x1 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 1) set.seed(671086) x2 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 1) set.seed(799837) x3 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 1) set.seed(289805) x4 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 3) set.seed(671086) x5 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 3) set.seed(799837) x6 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 3) set.seed(289805) x7 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 5) set.seed(671086) x8 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 5) set.seed(799837) x9 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 5) set.seed(289805) x10 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 10) set.seed(671086) x11 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 10) set.seed(799837) x12 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 10) xx <- 1:10 df <- data.frame(xx, x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12) reshapp <- reshape2::melt(df, id = "xx") NEWDAT <- data.frame(y = reshapp$value, x = reshapp$xx, sd = rep(rep(c(sd=1, sd=3, sd=5, sd=10), each = 10), each = 3),phi = rep(rep(c(.8, .9, .95), each = 10), 4)) ggplot(NEWDAT, aes(x = x, y = y)) + geom_line() + geom_point() + labs(x = 'lb', y = 'RMSE') + facet_grid(sd ~ phi, scales = "free_y") + theme_bw() + ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
原有输出效果:
解决方法
核心是通过factor设置分面变量的标签为plotmath表达式格式,再通过facet_grid的labeller参数调用label_parsed解析表达式,自动渲染希腊字母。
修改后的完整代码
## 模拟ARIMA(1, 0, 0)序列 set.seed(289805) x1 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 1) set.seed(671086) x2 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 1) set.seed(799837) x3 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 1) set.seed(289805) x4 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 3) set.seed(671086) x5 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 3) set.seed(799837) x6 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 3) set.seed(289805) x7 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 5) set.seed(671086) x8 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 5) set.seed(799837) x9 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 5) set.seed(289805) x10 <- arima.sim(n = 10, model = list(ar = 0.8, order = c(1, 0, 0)), sd = 10) set.seed(671086) x11 <- arima.sim(n = 10, model = list(ar = 0.9, order = c(1, 0, 0)), sd = 10) set.seed(799837) x12 <- arima.sim(n = 10, model = list(ar = 0.95, order = c(1, 0, 0)), sd = 10) xx <- 1:10 df <- data.frame(xx, x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12) reshapp <- reshape2::melt(df, id = "xx") NEWDAT <- data.frame(y = reshapp$value, x = reshapp$xx, sd = rep(rep(c(sd=1, sd=3, sd=5, sd=10), each = 10), each = 3),phi = rep(rep(c(.8, .9, .95), each = 10), 4)) # 新增代码:设置分面变量的标签为plotmath表达式格式 NEWDAT$sd <- factor(NEWDAT$sd, levels = c(1,3,5,10), labels = paste0("sd == ", c(1,3,5,10))) NEWDAT$phi <- factor(NEWDAT$phi, levels = c(0.8,0.9,0.95), labels = paste0("phi == ", c(0.8,0.9,0.95))) ggplot(NEWDAT, aes(x = x, y = y)) + geom_line() + geom_point() + labs(x = 'lb', y = 'RMSE') + # 修改facet_grid,增加labeller参数 facet_grid(sd ~ phi, scales = "free_y", labeller = label_parsed) + theme_bw() + ggplot2::scale_y_continuous(expand = c(0.0, 0.00))
修改说明
- 将
sd和phi转为因子时,标签使用==连接变量名和取值,符合R的plotmath语法要求 labeller = label_parsed会自动将因子标签解析为表达式,其中phi会自动渲染为希腊字母φ,sd保持原有英文字母格式- 不需要额外安装其他包,仅用ggplot2内置功能即可实现需求
内容的提问来源于stack exchange,提问作者Daniel James
相关产品推荐
相关产品推荐

