# HG changeset patch
# User alenail
# Date 1460676266 14400
# Node ID 4bbea0c6efd5b0628e2379b4d49340c899332ea3
# Parent  37b6e6c43ef85d4d9b467f283057d842a001e029
Uploaded

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)