Mercurial > repos > nilesh > rseqc
diff test-data/output.junctionSaturation_plot_r @ 59:dbedfc5f5a3c draft
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 247059e2527b66f1dbecf1e61496daef921040c3"
author | iuc |
---|---|
date | Sat, 18 Dec 2021 19:41:19 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.junctionSaturation_plot_r Sat Dec 18 19:41:19 2021 +0000 @@ -0,0 +1,12 @@ +pdf('output.junctionSaturation_plot.pdf') +x=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100) +y=c(0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1) +z=c(0,0,0,0,0,0,1,1,1,1,1,1,1,2,2,2,2,2,2,3) +w=c(0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,2) +m=max(0,0,0) +n=min(0,0,0) +plot(x,z/1000,xlab='percent of total reads',ylab='Number of splicing junctions (x1000)',type='o',col='blue',ylim=c(n,m)) +points(x,y/1000,type='o',col='red') +points(x,w/1000,type='o',col='green') +legend(5,0, legend=c("All junctions","known junctions", "novel junctions"),col=c("blue","red","green"),lwd=1,pch=1) +dev.off()