Repository 'pieplot_macs'
hg clone https://toolshed.g2.bx.psu.edu/repos/alenail/pieplot_macs

Changeset 9:4bbea0c6efd5 (2016-04-14)
Previous changeset 8:37b6e6c43ef8 (2016-04-14) Next changeset 10:9847c403d2d8 (2016-04-14)
Commit message:
Uploaded
modified:
pieplot_macs/pieplots_macs.py
b
diff -r 37b6e6c43ef8 -r 4bbea0c6efd5 pieplot_macs/pieplots_macs.py
--- a/pieplot_macs/pieplots_macs.py Thu Apr 14 19:05:32 2016 -0400
+++ b/pieplot_macs/pieplots_macs.py Thu Apr 14 19:24:26 2016 -0400
[
@@ -23,7 +23,7 @@
 from matplotlib import pyplot as plt
 matplotlib.rcParams['pdf.fonttype']=42
 matplotlib.rcParams['font.size']=14
-# import pandas as pd
+import pandas as pd
 from argparse import ArgumentParser
 
 def map_peaks(gene,peak,outfile,macsFlag):
@@ -41,31 +41,23 @@
         chrom = words[2]
 
         if macsFlag:
-            print "macsflag on"
             start = int(words[3])-1
             dist = abs(int(words[12]))
-            print dist
             maptype = words[14]
             if maptype == 'gene':
                 maptype = words[15]
 
         else:
-            print "macsflag off"
             try:
-                print words
                 start = int(words[3])
                 dist = abs(int(words[15]))
-                print dist
                 maptype = words[16]
                 if maptype == 'gene':
                     maptype = words[17]
             except:
-                print 'trigger except'
+                print 'there probably were not enough columns in the genefile'
                 pass
 
-        print "dist:"
-        print dist
-        print "that was dist"
 
         if chrom not in peaks:
             #new chrom
@@ -128,8 +120,6 @@
 
     args=parser.parse_args()
 
-    print args.MACS
-
     map_peaks(args.genefile, args.peakfile, args.outfile, args.MACS)