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