255
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to searchno edit summary
fo = open(fout, 'w')
fo = open(fout, 'w')
for ln in myopen(fin):
for ln in myopen(fin):
f = ln.strip().split()
l = len(f)
f = f[:6] + [ f[i*2] + ' '+f[i*2 +1] for i in xrange(3, l/2 )]
fo.write('\t'.join(f[:6]))
fo.write('\t')
if len(f[6:]) != len(flag):
die('Inconsistent length of ped and map files')
newMarker = [m for i, m in enumerate(f[6:]) if flag[i]]
fo.write('\t'.join(newMarker))
fo.write('\n')
#print marker[:10]
#die('test')
return True
return True
