上游数据示例：

image.png

基因A

data<-read.table("D:/work_tmp/column_diagram_plot/Gene1.percent.txt",header=T,sep="\t")
data$sample=as.character(data$sample)
data$sample=factor(data$sample,levels=c("6","7","9","11","14","17","19","20","22","24","26","27","28"))
pdf("Gene1.VAF.pdf", width=12, height=6)
ggplot(data = data, mapping = aes(x = sample, y =percent,fill=Item)) + geom_bar(stat = 'identity',position="dodge")+xlab('')+ylab('VAF of Gene1')+scale_fill_discrete(name = " ")+ theme(axis.title.x =element_text(size=14))
dev.off()

结果示例：

image.png

基因B

data1<-read.table("D:/work_tmp/column_diagram_plot/Gene2.percent.txt",header=T,sep="\t")
data1$sample=as.character(data1$sample)
data1$sample=factor(data1$sample,levels=c("3","5","8","25"))
pdf("Gene2.VAF.pdf", width=12, height=6)
ggplot(data = data1, mapping = aes(x = sample, y =percent,fill=Item)) + geom_bar(stat = 'identity',position="dodge")+xlab('')+ylab('VAF of Gene2')+scale_fill_discrete(name = " ")+ theme(axis.title.x =element_text(size=14))
dev.off()

image.png

两幅图组合在一起

a<-ggplot(data = data, mapping = aes(x = sample, y =percent,fill=Item)) + geom_bar(stat = 'identity',position="dodge")+xlab('')+ylab('VAF of Gene1')+scale_fill_discrete(name = " ")+ theme(axis.title.x =element_text(size=14)) + guides(fill=FALSE)

b<-ggplot(data = data1, mapping = aes(x = sample, y =percent,fill=Item)) + geom_bar(stat = 'identity',position="dodge")+xlab('')+ylab('VAF of Gene2')+scale_fill_discrete(name = " ")+ theme(axis.title.x =element_text(size=14))

pdf("arrange.EGFR.KRAS.VAF.pdf", width=12, height=6)
ggarrange(a,b,labels = c("B","C"),ncol = 2, nrow = 1)
dev.off()

示例：

image.png

累积柱状图

image.png

data<-read.table("D:/work_tmp/column_diagram_plot/data.txt",header=T,sep="\t")
data$sample=as.character(data$sample)
data$sample=factor(data$sample,levels=c("1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24","25","26","27","28"))
 pdf("mutation_statistics.pdf", width=15, height=6)
ggplot(data = data, mapping = aes(x = sample, y =num,fill=Item)) + geom_bar(stat = 'identity')+xlab('')+ylab('Number of sample')+scale_fill_discrete(name = " ")
dev.off()

示例：

image.png