barplot()R语言绘制状图

barplot()R语言绘制状图


数据:

OTU IDCRW
Membrane Transport0.1232550610.1314609080.136381709
Carbohydrate Metabolism0.1185871220.1152349450.101467081
Amino Acid Metabolism0.0990784610.1023070980.098133001
Replication and Repair0.0640276080.081496790.077107142
Energy Metabolism0.0630605790.0669996830.055660175
Translation0.0486662250.0525558320.048639947
Poorly Characterized0.0501483020.0474659650.049391479
Metabolism of Cofactors and Vitamins0.0463398870.0430211620.043328055
Nucleotide Metabolism0.0390221580.0371012490.037240203
Cellular Processes and Signaling0.0360890060.0262093820.039451668

绘图代码

d<-read.table("data.txt",header = T,row.names = 1,sep="\t")
d<-as.matrix(t(d))
barplot(d)
box("l")
pdf(file = "bar.pdf",h=7,w=14)
cc=c("#e41a1c","#377eb8","#4daf4a")
par(mar=c(5,20,4,2))
p=barplot(d,beside = TRUE,width = 1,xlab = "Relative abundance",bty = "n", col = cc,border = NA,horiz = T,xaxs="i",names.arg = NULL,axisnames=F)
legend("topright",legend = rownames(d), 
       fill = cc,  
       bty='n',xpd=T,inset = 0.1,cex = 1.2)
axis(side = 2, at = p[2,],labels =F)
text(y=p[2,]-0.25, x=-0.005, labels=colnames(d),cex=1,  xpd=TRUE, adj=1)
#text(0.1,0.2,labels = "sss")
box(bty="l")
dev.off()
绘图结果:
attachments-2020-03-hswbdgSI5e5fac7198953.png

  • 发表于 2020-03-04 21:26
  • 阅读 ( 2483 )
  • 分类:R

0 条评论

请先 登录 后评论
omicsgene
omicsgene

生物信息

654 篇文章

作家榜 »

  1. omicsgene 654 文章
  2. 安生水 325 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. 红橙子 78 文章
  6. CORNERSTONE 72 文章
  7. rzx 67 文章
  8. xun 66 文章