Bio_informatics_science_technology_image_beyond_1/2025-09-05_分别基于R的ggplot2和tidyplots包做的一个柱状图.qmd at main · iamhealthy/Bio_informatics_science_technology_image_beyond_1 · GitHub

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
87
88
89
90
91
92
93
94
95
---
title: ggplot2 vs tidyplots（柱状图）
subtitle: 分别基于ggplot2包和tidyplots包做的一个柱状图
date: 2025-09-05
toc-depth: 4
toc-expand: true
lang: en
---

2025-08-30的推文[电泳条带强度值分析](2025-08-30_电泳条带强度值分析.qmd)中的柱状图是用R的tidyplots包[@Engler2025]做的。Tidyplots是ggplot2包[@ggplot2]的一个延伸（extension）包。

在本推文中，将分别用ggplot2和tidyplots包来做这个柱状图，并比较两者的代码和效果。

测试数据下载链接：https://pan.baidu.com/s/1-hwYc2fkvuq_fPbtj5nrZQ?pwd=gkq7

## 1. 读取数据

```{r}
library(readr) # 读取csv文件

file_name <- "raw_data/2025-08-30_bands.csv" # 数据文件名
tbl <- file_name |> read_csv(show_col_types = FALSE) # 读取数据

tbl |> dim() # 数据维度(5行9列)
tbl |> names() # 数据列名
```

## 2. 用ggplot2包做柱状图

```{r}
library(ggplot2) # 加载ggplot2包
```

ggplot2包的作者：Hadley Wickham

<center>
![Hadlay Wickham。图源：https://hadley.nz/](images/hadley.jpg){width="40%"}
</center>

```{r}
tbl |>
  ggplot(aes(x = file_name, y = Band_ratio)) +
  geom_bar(stat = "identity", aes(fill = file_name), width = 0.5) +
  geom_text(aes(label = round(Band_ratio, 2)), vjust = -0.8, size = 5) +
  geom_point() +
  geom_line(group = 1) +
  labs(x = "Increased exposure time", y = "Band ratio (upper band / lower band)") +
  theme_classic() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1), legend.position = "none",
    axis.title = element_text(size = 10), axis.text = element_text(size = 8)) +
  scale_y_continuous(limits = c(0, max(tbl$Band_ratio) + 0.5), expand = c(0, 0),
    breaks = c(seq (0, max(tbl$Band_ratio) + 0.5, 0.5))) +
  scale_x_discrete(labels = c("20240104-3" = "Exposure 1",
                              "20240104-4" = "Exposure 2",
                              "20240104-5" = "Exposure 3",
                              "20240104-6" = "Exposure 4",
                              "20240104-7" = "Exposure 5"))
```

## 3. 用tidyplots包做柱状图

```{r}
library(tidyplots) # 加载tidyplots包
```

Tidyplots包的作者：Jan Broder Engler

<center>
![Jan Broder Engler。图源：https://jbengler.de/](images/broder_avatar.jpg){width="40%"}
</center>

```{r}
tbl |>
  tidyplot(x = file_name, y = Band_ratio, color = file_name) |>
  add_mean_bar() |>
  add_mean_value(color = "black", accuracy = 0.01, fontsize = 8) |>
  add_data_points(color = "black") |>
  add_mean_line(group = 1, color = "black") |>
  adjust_x_axis_title("Increased exposure time") |>
  adjust_y_axis_title("Band ratio (upper band / lower band)") |>
  adjust_x_axis(rotate_labels = 45) |>
  adjust_y_axis(limits = c(0, max(tbl$Band_ratio) + 0.5), breaks = c(seq(0, max(tbl$Band_ratio) + 0.5, 0.5))) |>
  rename_x_axis_labels(new_names = c("20240104-3" = "Exposure 1",
                                     "20240104-4" = "Exposure 2",
                                     "20240104-5" = "Exposure 3",
                                     "20240104-6" = "Exposure 4",
                                     "20240104-7" = "Exposure 5")) |>
  remove_legend()
```

更喜欢哪个风格？

[给我买杯茶🍵](给我买杯茶.qmd)

## References