data_select$group_list1<-group_list1
data_select$group_list2<-group_list2
colnames(data_select1)<-c('value','group_list1','group_list2')
group_list1为小分组
group_list2为大分组
library(ggplot2)
library('preprocessCore')
library(pheatmap)
library(gplots)
library(ggsignif)
compaired <- list(c('Naive','Primed'))
pp = ggplot(data_select,aes(x =group_list2,y=value,fill=group_list1))+
stat_boxplot(geom = "errorbar",width=0.75,color="black")+
geom_boxplot(size=1,aes(fill=group_list1,color='black'))+
scale_color_manual(values=c("black","black"))+
geom_signif(comparisons = compaired,step_increase = 0.1,map_signif_level = F,test = wilcox.test)+
labs(x=c(""),y="Log2(FPKM+1)"