Mercurial > repos > nilesh > rseqc
comparison test-data/output.junctionSaturation_plot.r @ 49:6b33e31bda10 draft
Uploaded tar based on https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc 1a3c419bc0ded7c40cb2bc3e7c87bfb01ddfeba2
| author | lparsons |
|---|---|
| date | Thu, 16 Jul 2015 17:43:43 -0400 |
| parents | |
| children | 09846d5169fa |
comparison
equal
deleted
inserted
replaced
| 48:2e6190c29c54 | 49:6b33e31bda10 |
|---|---|
| 1 pdf('output.junctionSaturation_plot.pdf') | |
| 2 x=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100) | |
| 3 y=c(0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1) | |
| 4 z=c(0,0,0,0,0,0,1,1,1,1,1,1,1,2,2,2,2,2,2,3) | |
| 5 w=c(0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,2) | |
| 6 m=max(0,0,0) | |
| 7 n=min(0,0,0) | |
| 8 plot(x,z/1000,xlab='percent of total reads',ylab='Number of splicing junctions (x1000)',type='o',col='blue',ylim=c(n,m)) | |
| 9 points(x,y/1000,type='o',col='red') | |
| 10 points(x,w/1000,type='o',col='green') | |
| 11 legend(5,0, legend=c("All junctions","known junctions", "novel junctions"),col=c("blue","red","green"),lwd=1,pch=1) | |
| 12 dev.off() |