R语言绘柱状图

         输入txt为两列，对齐全是各个名字，或者输入一列名字

1.  
   #设置工作目录
2.  
   setwd("F:\\R\\02柱状图")
3.  
    
4.  
   #读取输入文件
5.  
   rt=read.table("input.txt", header=T, sep="\t",comment.char = "", check.names =FALSE)
6.  
   #rt=read.table("也可以输入文件绝对地址", header=T, sep="\t",comment.char = "", check.names =FALSE)
7.  
   tb=table(c(as.vector(rt[,1]),as.vector(rt[,2]))) #对两列进行统计
8.  
   #tb=table(as.vector(rt[,1])) #对一列进行统计//不用的代码加#
9.  
   tb=sort(tb,decreasing =T) #频率高的排前面
10.  
     
11.  
    #输入每个基因的邻接点节点数目
12.  
    outTab=as.data.frame(tb)
13.  
    colnames(outTab)=c("Gene","Count")
14.  
    write.table(outTab,file="statResult.xls",sep="\t",quote=F,row.names=F)
15.  
     
16.  
    #定义柱状图显示基因数目
17.  
    showNum=30
18.  
    if(nrow(tb)<showNum){
19.  
    showNum=nrow(outTab)
20.  
    }
21.  
    n=as.matrix(tb)[1:showNum,]
22.  
     
23.  
    #绘制柱状图
24.  
    pdf(file="barplot.pdf",width=7,height=6)#设置图布大小
25.  
    par(mar=c(5,7,2,3),xpd=T)
26.  
    bar=barplot(n,horiz=TRUE,col="skyblue",names=FALSE,xlim=c(0,ceiling(max(n)/5)*5),xlab="Number of adjacent nodes")#y轴刻度是否显示，图的颜色，x轴刻度，x轴标题
27.  
    text(x=n*0.95,y=bar,n)
28.  
    text(x=-0.2,y=bar,label=names(n),xpd=T,pos=2)
29.  
    dev.off()

这篇好文章是转载于：学新通技术网