cdfdist.py

   1 import sys
   2
   3 def main(argv=None):
   4     if not argv:
   5         argv = sys.argv
   6
   7     if len(argv) < 4:
   8         print "usage: python %s bins percent infile" % sys.argv[0]
   9         sys.exit(1)
  10
  11     bins = int(argv[0])
  12     percent = int(argv[1])
  13     infilename = argv[2]
  14
  15     cdfDist(bins, percent, infilename)
  16
  17
  18 def cdfDist(bins, percent, infilename):
  19     infile = open(infilename)
  20     binsList = [0] * bins
  21
  22     for line in infile:
  23         fields = line.strip().split()
  24         index = 0
  25         for binCdf in fields[-1 * bins:]:
  26             if int(binCdf) > percent:
  27                 binsList[index] += 1
  28                 break
  29
  30             index += 1
  31
  32     infile.close()
  33     print binsList
  34
  35
  36 if __name__ == "__main__":
  37     main(sys.argv)