Mercurial > repos > alenail > pieplot_macs

changeset 6:d7baa700a930 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author alenail
date Thu, 14 Apr 2016 17:45:12 -0400
parents 01748da366c3
children b9b2be216c05
files pieplot_macs/pieplots_macs.py pieplot_macs/pieplots_macs.xml
diffstat 2 files changed, 15 insertions(+), 9 deletions(-) [+]
line wrap: on
line diff
--- a/pieplot_macs/pieplots_macs.py	Thu Apr 14 17:40:18 2016 -0400
+++ b/pieplot_macs/pieplots_macs.py	Thu Apr 14 17:45:12 2016 -0400
@@ -41,7 +41,16 @@
         chrom = words[2]
 
 
-        if not macsFlag:
+        if macsFlag:
+            print "MacsFlag off"
+            start = int(words[3])-1
+            dist = abs(int(words[12]))
+            maptype = words[14]
+            if maptype == 'gene':
+                maptype = words[15]
+
+        else:
+            print "macsFlag on"
             try:
                 start = int(words[3])
                 dist = abs(int(words[15]))
@@ -51,13 +60,8 @@
             except:
                 pass
 
-        else:
-            start = int(words[3])-1
-            dist = abs(int(words[12]))
-            maptype = words[14]
-            if maptype == 'gene':
-                maptype = words[15]
-
+        print "dist:"
+        print dist
 
         if chrom not in peaks:
             #new chrom
@@ -120,6 +124,8 @@
 
     args=parser.parse_args()
 
+    print args.MACS
+
     map_peaks(args.genefile, args.peakfile, args.outfile, args.MACS)
 
 
--- a/pieplot_macs/pieplots_macs.xml	Thu Apr 14 17:40:18 2016 -0400
+++ b/pieplot_macs/pieplots_macs.xml	Thu Apr 14 17:45:12 2016 -0400
@@ -18,5 +18,5 @@
   <outputs>
     <data name="out" format="pdf" hidden="false" />
   </outputs>
-  <help>help</help>
+  <help></help>
 </tool>