Mercurial > repos > nilesh > rseqc

diff test-data/output.clipping_profile.r @ 49:6b33e31bda10 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/output.clipping_profile.r	Thu Jul 16 17:43:43 2015 -0400
@@ -0,0 +1,5 @@
+pdf("output.clipping_profile.pdf")
+read_pos=c(0,1,2,3,4,5,6,7,8,9,44,45,46,47,48,49,50)
+count=c(16,12,11,8,6,5,1,1,1,1,1,2,2,2,3,4,4)
+plot(read_pos,1-(count/40),col="blue",main="clipping profile",xlab="Position of reads",ylab="Mappability",type="b")
+dev.off()