Mercurial > repos > lparsons > rseqc

annotate test-data/output.junctionSaturation_plot.r @ 1:cd2daa987b69 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit 2400cce6e3918648fda8f2a8d8485ab251a58827
author lparsons
date Tue, 07 Jul 2015 14:14:04 -0400
parents 306127635959
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
1 pdf('output.junctionSaturation_plot.pdf')
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
2 x=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100)
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
3 y=c(0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1)
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
4 z=c(0,0,0,0,0,0,1,1,1,1,1,1,1,2,2,2,2,2,2,3)
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
5 w=c(0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,2)
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
6 m=max(0,0,0)
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
7 n=min(0,0,0)
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
8 plot(x,z/1000,xlab='percent of total reads',ylab='Number of splicing junctions (x1000)',type='o',col='blue',ylim=c(n,m))
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
9 points(x,y/1000,type='o',col='red')
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
10 points(x,w/1000,type='o',col='green')
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
11 legend(5,0, legend=c("All junctions","known junctions", "novel junctions"),col=c("blue","red","green"),lwd=1,pch=1)
306127635959 planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit aeb25d807817746dd6957f30ce2070662cc10e91
lparsons
parents:
diff changeset
12 dev.off()