src/hg/tcga/scripts/segToBed15.py 1.2
1.2 2009/09/18 23:40:43 jsanborn
fixed bug with sparse data
Index: src/hg/tcga/scripts/segToBed15.py
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/tcga/scripts/segToBed15.py,v
retrieving revision 1.1
retrieving revision 1.2
diff -b -B -U 4 -r1.1 -r1.2
--- src/hg/tcga/scripts/segToBed15.py 6 Aug 2009 22:56:50 -0000 1.1
+++ src/hg/tcga/scripts/segToBed15.py 18 Sep 2009 23:40:43 -0000 1.2
@@ -17,9 +17,9 @@
sample = data[0]
chrom = "chr" + data[1]
start = int(data[2])
stop = int(data[3])
- val = float(data[5])
+ val = float(data[4])
if sample not in segData:
segData[sample] = {}
@@ -118,9 +118,13 @@
s += '1000' + '\t' + '+' + '\t' + start + '\t' + stop + '\t'
s += '0' + '\t' + '1' + '\t' + length + '\t' + '0,' + '\t' + numSamples + '\t'
vals = []
for sample in samples:
+ if c in mapped[sample]:
val = mapped[sample][c].vals[i]
+ else:
+ val = MISSING_VAL
+
if val == MISSING_VAL:
valStr = ''
else:
valStr = str(val)