Mercurial > repos > iuc > multiqc

comparison test-data/junction_saturation.txt @ 23:d01bcfd04189 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/multiqc commit 8a91236cee4d408ae2b53a3e9b6daebc332d631a
author iuc
date Sat, 10 Dec 2022 11:20:55 +0000
parents
children
comparison
equal deleted inserted replaced
22:85613e9ad431 23:d01bcfd04189
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,0,0,0,0,0,0,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,2,2,2,2,2,2,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()