Mercurial > repos > xuebing > sharplabtool
diff collapseBed2.py @ 14:76e1b1b21cce default tip
Deleted selected files
author | xuebing |
---|---|
date | Tue, 13 Mar 2012 19:05:10 -0400 |
parents | 292186c14b08 |
children |
line wrap: on
line diff
--- a/collapseBed2.py Sat Mar 10 08:17:36 2012 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,36 +0,0 @@ -''' -collapse intervals -''' - -def collapseInterval_strand(filename,c_strand,c_score): - # keeping max column c - uniqintv = {} - data = {} - f = open(filename) - header = f.readline() - if 'chr' in header: - flds = header.strip().split('\t') - key = '\t'.join([flds[0],flds[1],flds[2],flds[c_strand]]) - uniqintv[key] = float(flds[c_score]) - data[key] = flds - for line in f: - flds = line.strip().split('\t') - key = '\t'.join([flds[0],flds[1],flds[2],flds[c_strand]]) - if not uniqintv.has_key(key): - uniqintv[key] = float(flds[c_score]) - data[key] = flds - elif uniqintv[key] < float(flds[c_score]): - uniqintv[key] = float(flds[c_score]) - data[key] = flds - - f.close() - for key in uniqintv.keys(): - print '\t'.join(data[key]) - -import sys - -if sys.argv[2] == '0':#ignore strand - sys.argv[2] = 1 -if sys.argv[3] == '0':# ignore score - sys.argv[3] = 2 -collapseInterval_strand(sys.argv[1],int(sys.argv[2])-1,int(sys.argv[3])-1)