最近在逛R语言数据分析与可视化是发现了有个ggpubr包可以绘制箱线图的显著性标记，感觉蛮有用的，然后就在网上搜了下显著性标记相关资料，发现除了ggpubr包外，发现一个ggsignif包！查了下原来这个包是专门配套ggplot2来给图添加显著性标记的，所以收藏并使用了下 包的使用示例即说明可以参看Github链接：ggsignif
还有一个公众号也对这包的使用进行了介绍：神奇小工具丨ggplot绘图显著性添加工具---ggsignif
下面是我对这个包的初步使用
具体函数的参数使用可以看ggsignif包的pdf文档
先下载使用数据，这里就拿TCGA数据来练练手，当然不是常规的去官网下载数据啦，而是使用RTCGA包，然后利用RTCGA.mRNA包里的各个样本的表达数据，先加载会用到的包
library(RTCGA)
library(RTCGA.mRNA)
library(dplyr)
library(ggplot2)
library(ggsignif)获取乳腺癌的表达谱数据，并挑取RCA1,BRCA2,TP53三个基因在各个样本中的表达数据作为测试数据
exprSet <- BRCA.mRNA %>% as_tibble()
data <- select(exprSet, BRCA1,BRCA2,TP53) %>% stack()
colnames(data) <- c("values", "sample")绘制个箱线图瞧瞧
p <- ggplot(data, aes(x = sample, y = values)) +
geom_boxplot() +
theme_classic()
p
然后使用ggsignif包绘制有显著性标记的箱线图，先设定比较组，并选择差异分析的统计学方法，绘图
p <- ggplot(data, aes(x = sample, y = values)) +
geom_boxplot() +
theme_classic() +
geom_signif(comparisons = compared_list, test = wilcox.test, step_increase = 0.2)
p
如果像一些图上并不是显示p值，而是显示***，那么可以再添加map_signif_level参数
p <- ggplot(data, aes(x = sample, y = values)) +
geom_boxplot() +
theme_classic() +
geom_signif(comparisons = compared_list, test = wilcox.test, step_increase = 0.2, map_signif_level = T)
p
如果并不想标注p值，而是想标注认为设定的数字或者字符，那么设定annotation参数以及搭配其他一些定位参数，如下
p <- ggplot(data, aes(x = sample, y = values)) +
geom_boxplot() +
theme_classic() +
geom_signif(annotations = c("First", "Second"), y_position = c(2.5, 3.5), xmin = c(1,2), xmax = c(2,3))
p
虽然网上介绍的ggsignif包都是对boxplot进行操作的，但是作者并没有说局限于箱线图（小提琴图）等，说明还支持其他ggplot2的图，所以我也对柱状图尝试了下
先生成一个测试数据框
df <- data.frame(state = rep(c("A", "B", "C", "D", "E"), 2),
value = c(1:5,3:7),
group = c(rep("one", 5), rep("two", 5))
)当然我们需要知道的是柱状图是不能求显著性的，但是还可以用来跟上图一样的人为标注以示说明，如下
p <- ggplot(df, aes(x = state, y = value, fill = group)) +
geom_bar(stat = "identity", position = "dodge") +
theme_classic() +
geom_signif(annotations = c(0.001, "***"), y_position = c(4,5), xmin = c(0.75, 1.75), xmax = c(1.25, 2.25),
tip_length = c(c(0.2, 0.05),c(0.2, 0.05)), vjust = -1)
p
以上是对ggsignif包的简单尝试，这个ggplot2插件真心不错
本文出自于http://www.bioinfo-scrounger.com转载请注明出处