【ggplot2】改坐标轴breaks标签，加数学表达式，数学表达式legend对齐(align the text of legend)_ggplot2中坐标轴标注和坐标轴对齐

figure_display_batch <- function(res) {
  res_tmp <- res
  size <- dim(res_tmp)[3]
  res_list <- list()
  for (i in 1:size) {
    res_list[[i]] <- res_tmp[, , i]
  }
  res_list <- do.call("rbind", res_list)
  col_name <- c("Prediction Error", "Model Size")
  colnames(res_list) <- col_name
  res_list <- as.data.frame(res_list)
  rownames(res_list) <- NULL
  row_name <-  c(bquote("G"~L[2]~"PDAS"), as.expression(bquote("P"~L[2]~"PDAS")),
                 as.expression(bquote(L[0]~"-SPDAS")), as.expression(bquote(L[0]~"-GPDAS")),
                 as.expression(bquote("LASSO")), "SCAD",
                 "MCP", "Elastic-Net",
                 "Relaxed-LASSO", as.expression(bquote(L[0]~L[2]~"-CD")),
                 as.expression(bquote(L[0]~L[2]~"-CDPSI")),
                 as.expression(bquote(L[0]~"-CD")), as.expression(bquote(L[0]~"-CDPSI")))

  res_list[["method"]] <- rep(row_name, size)
  res_list[["method"]] <- factor(res_list[["method"]],
                                 levels =   c(bquote("G"~L[2]~"PDAS"), as.expression(bquote("P"~L[2]~"PDAS")),
                                              as.expression(bquote(L[0]~"-SPDAS")), as.expression(bquote(L[0]~"-GPDAS")),
                                              as.expression(bquote("LASSO")), "SCAD",
                                              "MCP", "Elastic-Net",
                                              "Relaxed-LASSO", as.expression(bquote(L[0]~L[2]~"-CD")),
                                              as.expression(bquote(L[0]~L[2]~"-CDPSI")),
                                              as.expression(bquote(L[0]~"-CD")), as.expression(bquote(L[0]~"-CDPSI")))
  )
  library(tidyr)
  plot_data = gather(res_list, metric, value, -method)
  plot_data$metric = factor(plot_data$metric , levels =  c("Prediction Error", "Model Size", "Infinity Norm"))
 plot_data$method = factor(plot_data$method, levels =  c(bquote("G"~L[2]~"PDAS"), as.expression(bquote("P"~L[2]~"PDAS")),
                                                         as.expression(bquote(L[0]~"-SPDAS")), as.expression(bquote(L[0]~"-GPDAS")),
                                                         as.expression(bquote("LASSO")), "SCAD",
                                                         "MCP", "Elastic-Net",
                                                         "Relaxed-LASSO", as.expression(bquote(L[0]~L[2]~"-CD")),
                                                         as.expression(bquote(L[0]~L[2]~"-CDPSI")),
                                                         as.expression(bquote(L[0]~"-CD")), as.expression(bquote(L[0]~"-CDPSI")))
 )

 calc_stat <- function(x) {
   coef <- 5
   n <- sum(!is.na(x))
   # calculate quantiles
   stats <- quantile(x, probs = c(0, 0.25, 0.5, 0.75, 1))
   names(stats) <- c("ymin", "lower", "middle", "upper", "ymax")
   return(stats)
 }
 # color = c('#b2182b', '#CD443E', '#E76F51', '#EE8959','#F4A261','#EFB366', '#E9C46A', '#F3C891','#BABB74', '#8AB17D',
 #           '#5AA786', '#2A9D8F', '#287271', '#264653'
 # )
 color = c('#B2182B','#BF0C49','#CC0066','#B34D33','#AD6027', '#A6731A', '#A0860D', '#999900',
           '#809940', '#669980', '#3399FF', '#2D70A9','#2D70A9', '#264653')
 p = ggplot(plot_data, aes(x = method, y = value, fill = method), coef = 5) +
   # geom_boxplot() + #coef = 5
    stat_summary(fun.data = calc_stat, geom="boxplot", width = 0.75, alpha = 0.8) +
    facet_wrap(~metric, scales = "free") +
   scale_fill_manual(values = color, labels = c(as.expression(bquote("G"~L[2]~"PDAS")), as.expression(bquote("P"~L[2]~"PDAS")),
                                                as.expression(bquote(L[0]~"-SPDAS")), as.expression(bquote(L[0]~"-GPDAS")),
                                                as.expression(bquote("LASSO")), "SCAD",
                                                "MCP", "Elastic-Net",
                                                "Relaxed-LASSO", as.expression(bquote(L[0]~L[2]~"-CD")),
                                                as.expression(bquote(L[0]~L[2]~"-CDPSI")),
                                                as.expression(bquote(L[0]~"-CD")), as.expression(bquote(L[0]~"-CDPSI")))
                     ) +
   # 更改坐标轴的breaks标签
   scale_x_discrete(labels = c(as.expression(bquote("G"~L[2]~"PDAS")), as.expression(bquote("P"~L[2]~"PDAS")),
                               as.expression(bquote(L[0]~"-SPDAS")), as.expression(bquote(L[0]~"-GPDAS")),
                               as.expression(bquote("LASSO")), "SCAD",
                               "MCP", "Elastic-Net",
                               "Relaxed-LASSO", as.expression(bquote(L[0]~L[2]~"-CD")),
                               as.expression(bquote(L[0]~L[2]~"-CDPSI")),
                               as.expression(bquote(L[0]~"-CD")), as.expression(bquote(L[0]~"-CDPSI")))
   ) +
    theme_bw()+
    theme(
      legend.position = "bottom",
          panel.grid = element_blank(),
          axis.text.x = element_text(angle = 45, hjust = 0.5, vjust = 0.5),
          axis.title = element_blank(),
      legend.text.align = 0)

  p
}
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86