ggplot2: 添加显著性

1. ggpubr

https://blog.csdn.net/zhouhucheng00/article/details/106391872/
https://www.jianshu.com/p/77f12664540b
https://www.jianshu.com/p/18c29f89c7a1

差异统计

> head(kk)
        gene         RPKM tissue  sample
2       18-w  1.324219345      A Drone_1
3      5-HT1  0.115036222      A Drone_1
4 5-HT2alpha  0.156672970      A Drone_1
5  5-HT2beta  0.000000000      A Drone_1
6      5-ht7  0.003033164      A Drone_1
7         A4 12.485018064      A Drone_1

> ggplot()+geom_boxplot(data=kk, aes(x=tissue, y=RPKM, fill=sample), outlier.color = NA)+scale_y_log10()+theme_cowplot()+scale_color_manual(values = hue_pal()(4))  '#示例一

> Statistics_result_of_t_test <- as.data.frame(compare_means(RPKM~sample, kk, method = "t.test", paired = FALSE, group.by = "tissue"))  #进行差异统计
> Statistics_result_of_t_test
   tissue  .y.  group1  group2          p p.adj p.format p.signif method
1       A RPKM Drone_1 Drone_2 0.91564803     1    0.916       ns T-test
2       A RPKM Drone_1 Drone_3 0.40329980     1    0.403       ns T-test
3       A RPKM Drone_1 Drone_4 0.26966910     1    0.270       ns T-test
4       A RPKM Drone_2 Drone_3 0.35598501     1    0.356       ns T-test
5       A RPKM Drone_2 Drone_4 0.23113466     1    0.231       ns T-test
6       A RPKM Drone_3 Drone_4 0.82903388     1    0.829       ns T-test
7       H RPKM Drone_1 Drone_2 0.65314577     1    0.653       ns T-test
8       H RPKM Drone_1 Drone_3 0.10048293     1    0.100       ns T-test
9       H RPKM Drone_1 Drone_4 0.09835627     1    0.098       ns T-test
10      H RPKM Drone_2 Drone_3 0.12251987     1    0.123       ns T-test
11      H RPKM Drone_2 Drone_4 0.12185849     1    0.122       ns T-test
12      H RPKM Drone_3 Drone_4 0.94264035     1    0.943       ns T-test
13      T RPKM Drone_1 Drone_2 0.62284377     1    0.623       ns T-test
14      T RPKM Drone_1 Drone_3 0.24356479     1    0.244       ns T-test
15      T RPKM Drone_1 Drone_4 0.19417880     1    0.194       ns T-test
16      T RPKM Drone_2 Drone_3 0.32061625     1    0.321       ns T-test
17      T RPKM Drone_2 Drone_4 0.26185337     1    0.262       ns T-test
18      T RPKM Drone_3 Drone_4 0.93052349     1    0.931       ns T-test

2. ggsignif

https://blog.csdn.net/lalaxumelala/article/details/109389602