Mercurial > repos > nilesh > rseqc
comparison test-data/output.junctionSaturation_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 |
comparison
equal
deleted
inserted
replaced
50:f242ee103277 | 51:09846d5169fa |
---|---|
1 pdf('output.junctionSaturation_plot.pdf') | 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) | 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) | 3 y=c(0,0,0,0,0,0,1,1,1,1,1,1,1,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) | 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) | 5 w=c(0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,2) |
6 m=max(0,0,0) | 6 m=max(0,0,0) |
7 n=min(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)) | 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') | 9 points(x,y/1000,type='o',col='red') |
10 points(x,w/1000,type='o',col='green') | 10 points(x,w/1000,type='o',col='green') |