Mercurial > repos > youngkim > ezbamqc

view ezBAMQC/test-data/output/data/smp1.read_distr.r @ 20:9de3bbec2479 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author youngkim
date Thu, 31 Mar 2016 10:10:37 -0400
parents dfa3745e5fd8
children
line wrap: on
line source

png("/sonas-hs/bsr/hpc/data/yjin/test_BAMqc/exp/test1/figs/smp1.read_distr.png",width=500,height=500,units="px")
M=c(24849,1533,16115,8729,205,453,7563,14335)
Mname=c("CDS","5UTR","3UTR","Intron","TSS_Up_1Kb","TES_Down_1Kb","rRNA","Others")
val = barplot(M,xlab="",space=1,ylab="Read Counts",col="blue",border="NA")
text(x=seq(val[1],val[8],by=2),y=rep(0,8),srt=60,adj=0,offset=2,pos=1,xpd=T,labels=Mname)
dev.state = dev.off()