组间差异分析就要这样可视化！

欢迎关注”生信修炼手册”!

在之前的文章中，我们分享了多个基因差异分析的可视化，使用的是ggpubr这个R包，ggpubr在标记p值时，可以根据指定的差异分组自动添加组间的连线，非常方便，但是无法指定添加的p值的位置，在某些时候会缺乏灵活性，今天要介绍的是另外一个R包ggsignif，其帮助手册链接如下

https://cran.r-project.org/web/packages/ggsignif/vignettes/intro.html

首先我们用示例数据跑一跑

> library(ggplot2)
> library(ggsignif)
> head(iris)
  Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1          5.1         3.5          1.4         0.2  setosa
2          4.9         3.0          1.4         0.2  setosa
3          4.7         3.2          1.3         0.2  setosa
4          4.6         3.1          1.5         0.2  setosa
5          5.0         3.6          1.4         0.2  setosa
6          5.4         3.9          1.7         0.4  setosa
> ggplot(iris, aes(x = Species, y = Sepal.Length)) +
+   geom_boxplot() +
+   geom_signif(
+     comparisons = list(c("versicolor", "virginica"))
+ )

效果图