现有一文件命名为“cons_count.txt”:
score A B
0 0.170682731 0.192926045
0.0-0.1 0.008032129 0.007073955
0.1-0.2 0.010040161 0.016720257
0.2-0.3 0.032128514 0.025723473
0.3-0.4 0.03815261 0.036655949
0.4-0.5 0.102409639 0.091318328
0.5-0.6 0.156626506 0.149196141
0.6-0.7 0.174698795 0.21221865
0.7-0.8 0.164658635 0.195498392
0.8-0.9 0.106425703 0.064951768
0.9-1 0.036144578 0.007717042
用四行命令可以画出两个系列的柱状图
read.table("cons_count.txt",header=T)->z
barplot(rbind(z[,2],z[,3]),names.arg=z[,1],col=c("red","green"),width=c(0.1,0.1),space=c(0,1),beside=T,border=NA,ylim=c(0,0.25))
legend("topleft",c(colnames(z)[2],colnames(z)[3]),fill=c("red","green"),bty="n",border=NA)
box()
另外,还有一个博文http://abbypan.blogspot.hk/2009/04/r.html 介绍的更详细一点。
https://wap.sciencenet.cn/blog-824692-644694.html
下一篇:
UCSC genome browser 镜像搭建