Boxplot+geom_signif

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)",title = markers[i])+
  theme(panel.background = element_blank(),
         axis.line = element_line(colour="black"),
         axis.text.x=element_text(size=14,angle = 75,vjust = 0.5,colour = 'black'),
         axis.text.y=element_text(size=14,colour = 'black'), 
         axis.title.y=element_text(size = 14,colour = 'black'), 
         axis.title.x=element_text(size = 14,colour = 'black'), 
         plot.title = element_text(size=18,face="bold",hjust = 0.5,colour = 'black'),
         #panel.grid.major = element_blank(), 
         #panel.grid.minor = element_blank())
  )

pp