横向堆积柱状图,并添加标签
setwd('D:\\10_NM_MH_3\\03_diffpeak\\02_diff')
data<-read.csv('1plot_peak.txt',sep = '\t',header = 0)
library(reshape2)
library(ggpubr)
library(stringr)
colnames(data)<-c('baifen','feature','number','Sample')
unique(data$feature)
# data$feature<-str_replace(data$feature,"Distal Intergenic","Intergenic")
data$feature<-str_replace(data$feature,"Downstream ","Downstream (<=3kb)")
data$feature<-str_replace(data$feature,"Promoter","Promoter (<=1kb)")
unique(data$feature)
unique(data$Sample)
data$Sample<-str_replace(data$Sample,"BCD-bised1",'BbMHSs')
data$Sample<-str_replace(data$Sample,"BCD-bised2",'CbMHSs')
data$Sample<-str_replace(data$Sample,"BCD-common",'cMHSs')
data$feature<- factor(data$feature,
levels= c("Distal Intergenic","Downstream (<=3kb)","3' UTR","Intron ","Exon ","5' UTR",
"Promoter (<=1kb)"),ordered = FALSE)
class(data$baifen)
head(data)
data$baifen<-round(data$baifen,2)
data$ratio<-paste(data$baifen*100, "%", sep='')
> head(data)
baifen feature number Sample ratio
1 0.04 3' UTR 66 BbMHSs 4%
2 0.02 5' UTR 35 BbMHSs 2%
3 0.38 Distal Intergenic 666 BbMHSs 38%
4 0.14 Downstream (<=3kb) 244 BbMHSs 14%
5 0.03 Exon 59 BbMHSs 3%
6 0.08 Intron 144 BbMHSs 8%
ggplot(data,aes(x=Sample,fill=feature,y=number))+
geom_bar(stat = "identity",position ='fill',colour="black")+
# geom_bar(position = position_stack(), stat ="identity", width = .7) +
labs(x='',y='')+
# facet_grid(Sample~.,scales="free",space="free")+
scale_fill_brewer(type="seq",palette="Set2",)+
guides(fill = guide_legend(title = '',reverse=TRUE))+
geom_text(aes(label = ratio), position = 'fill',hjust=1.5,vjust=0.5)+
# geom_text(aes(label=number), #position=position_dodge(width=0.9),
# angle =0,colour='black',,
# fontface = 'bold')+
# scale_fill_hue('',type="seq",palette="Set1")+
coord_flip()+theme_bw()+
theme(axis.title = element_text(face = 'bold',
size = "16",color = "black"),
# legend.position = 'right',
axis.text.x = element_text(color = "black",face = 'bold',
size = 13, hjust = 0.5, vjust = 0.5),
axis.text.y = element_text(face = 'bold',size =13,color = "black"),
legend.text = element_text(face = 'bold',color = "black",size = 10),
legend.title = element_text(face = 'bold',color = "black",size = 10),
panel.grid.major.y = element_blank(),
panel.grid.minor.y = element_blank(),
panel.grid.major.x = element_blank(),
panel.grid.minor.x = element_blank(),
legend.position="right",
strip.text.x = element_text( face="bold",size = 14),
#size=8,angle=75
strip.text.y = element_text( face="bold",size = 14),
panel.border = element_rect(colour = "black", fill=NA, size=1),
strip.background = element_rect(colour = "white", fill = "grey"),
plot.title = element_text(face = "bold",color = "black",lineheight=.8, hjust=0.5, size =11),)