# HG changeset patch # User jingchunzhu@gmail.com # Date 1438317735 25200 # Node ID f55d3b781adb5454dd24ee7fb7953634fad20c7e # Parent 14c406f8fca62be91ee2f327ccd6638cb9039db0 check file input ignore bad lines` diff -r 14c406f8fca6 -r f55d3b781adb ttest/stats.py --- a/ttest/stats.py Wed Oct 22 18:14:35 2014 -0700 +++ b/ttest/stats.py Thu Jul 30 21:42:15 2015 -0700 @@ -28,6 +28,8 @@ outfile.write("%s\tStatistic\tpValue\tMedian1\tMedian2\tDelta\n" % (headerTokens[0])) for row in genomicFp: tokens = row.rstrip().split("\t") + if len(tokens)!= len(headerTokens): + continue symbol = tokens[0] groupATokens = [tokens[ii] for ii in columnsGroupA] groupBTokens = [tokens[ii] for ii in columnsGroupB]