Mercurial > repos > alenail > pieplot_macs

changeset 9:4bbea0c6efd5 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 19:24:26 -0400
parents 37b6e6c43ef8
children 9847c403d2d8
files pieplot_macs/pieplots_macs.py
diffstat 1 files changed, 2 insertions(+), 12 deletions(-) [+]
line wrap: on
line diff
--- 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)