推荐很好的R语言ggplot2科研数据可视化参考资料

data<-read.delim("data/20220622/example_data.csv",sep=",",header=TRUE)
library(ggplot2)
my_pal <- rcartocolor::carto_pal(n = 8, name = "Bold")[c(1, 3, 7, 2)]
g <- ggplot(data, aes(x = group, y = value, color = group, fill = group)) +
  scale_y_continuous(breaks = 1:9) +
  scale_color_manual(values = my_pal, guide = "none") +
  scale_fill_manual(values = my_pal, guide = "none")
g + 
  geom_boxplot(
    width = .2, fill = "white",
    size = 1.5, outlier.shape = NA
  ) +
  ggdist::stat_halfeye(
    adjust = .33, ## bandwidth
    width = .67, 
    color = NA, ## remove slab interval
    position = position_nudge(x = .15)
  ) +
  gghalves::geom_half_point(
    side = "l", 
    range_scale = .3, 
    alpha = .5, size = 3
  ) -> p1


g + 
  geom_violin(
    aes(fill = group, fill = after_scale(colorspace::lighten(fill, .5))),
    size = 1.2, bw = .2
  ) + 
  geom_boxplot(
    fill = "white",  size = 1.2, width = .2, 
    outlier.shape = NA, coef = 0
  ) +
  geom_point(
    position = position_jitter(width = .03, seed = 0),
    size = 5, alpha = .5
  ) +
  geom_point(
    position = position_jitter(width = .03, seed = 0),
    size = 5, stroke = .7, shape = 1, color = "black"
  ) -> p2

library(patchwork)
p1+p2