Mercurial > repos > nilesh > rseqc
diff test-data/output.DupRate_plot.r @ 51:09846d5169fa draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
| author | iuc |
|---|---|
| date | Tue, 14 Mar 2017 10:23:21 -0400 |
| parents | 6b33e31bda10 |
| children |
line wrap: on
line diff
--- a/test-data/output.DupRate_plot.r Tue May 03 16:36:57 2016 -0400 +++ b/test-data/output.DupRate_plot.r Tue Mar 14 10:23:21 2017 -0400 @@ -4,11 +4,11 @@ seq_uniqRead=c(40) pos_occ=c(1) pos_uniqRead=c(40) -plot(pos_occ,log10(pos_uniqRead),ylab='Number of Reads (log10)',xlab='Frequency',pch=4,cex=0.8,col='blue',xlim=c(1,500),yaxt='n') +plot(pos_occ,log10(pos_uniqRead),ylab='Number of Reads (log10)',xlab='Occurrence of read',pch=4,cex=0.8,col='blue',xlim=c(1,500),yaxt='n') points(seq_occ,log10(seq_uniqRead),pch=20,cex=0.8,col='red') ym=floor(max(log10(pos_uniqRead))) -legend(300,ym,legend=c('Sequence-base','Mapping-base'),col=c('blue','red'),pch=c(4,20)) +legend(300,ym,legend=c('Sequence-based','Mapping-based'),col=c('blue','red'),pch=c(4,20)) axis(side=2,at=0:ym,labels=0:ym) -axis(side=4,at=c(log10(pos_uniqRead[1]),log10(pos_uniqRead[2]),log10(pos_uniqRead[3]),log10(pos_uniqRead[4])), labels=c(round(pos_uniqRead[1]*100/sum(pos_uniqRead)),round(pos_uniqRead[2]*100/sum(pos_uniqRead)),round(pos_uniqRead[3]*100/sum(pos_uniqRead)),round(pos_uniqRead[4]*100/sum(pos_uniqRead)))) +axis(side=4,at=c(log10(pos_uniqRead[1]),log10(pos_uniqRead[2]),log10(pos_uniqRead[3]),log10(pos_uniqRead[4])), labels=c(round(pos_uniqRead[1]*100/sum(pos_uniqRead*pos_occ)),round(pos_uniqRead[2]*100/sum(pos_uniqRead*pos_occ)),round(pos_uniqRead[3]*100/sum(pos_uniqRead*pos_occ)),round(pos_uniqRead[4]*100/sum(pos_uniqRead*pos_occ)))) mtext(4, text = "Reads %", line = 2) dev.off()