每日R语言绘图（一）——barplotStat条形图统计

绘出后的图如下：

首先设立一个工作文件夹，里面放代码与数据输入文件，如图

#设置工作目录

setwd("C:\\Users\\Qptdada\\OneDrive\\桌面\\bioplot\\01.barplotStat条形图统计")

#读取输入文件

rt=read.table("input.txt", header=T, sep="\t",comment.char = "", check.names =FALSE)

tb=table(c(as.vector(rt[,1]),as.vector(rt[,2])))    #对两列进行统计

#tb=table(as.vector(rt[,1]))        #对一列进行统计//不用的代码加#

tb=sort(tb,decreasing =T) #频率高的排前面

#输入每个基因的邻接点节点数目

outTab=as.data.frame(tb)

colnames(outTab)=c("Gene","Count")

write.table(outTab,file="statResult.xls",sep="\t",quote=F,row.names=F)

#定义柱状图显示基因数目

showNum=30

if(nrow(tb)<showNum){

showNum=nrow(outTab)

}

n=as.matrix(tb)[1:showNum,]

#绘制柱状图

pdf(file="barplot.pdf",width=7,height=6)

par(mar=c(5,7,2,3),xpd=T)

bar=barplot(n,horiz=TRUE,col="skyblue",names=FALSE,xlim=c(0,ceiling(max(n)/5)*5),xlab="Number of adjacent nodes")

text(x=n*0.95,y=bar,n)

text(x=-0.2,y=bar,label=names(n),xpd=T,pos=2)

dev.off()

#####################

下面有演示代码和演示用的数据文件

链接：https://pan.baidu.com/s/1ku50aLOBvE67c1EKSecoNw

提取码：cnxf

--来自百度网盘超级会员V4的分享