#!/net/em-stor4/Volumes/data/wangyr/local/bin/python2.7
from argparse import ArgumentParser
from sys import stderr, stdout
import myutil
if __name__=="__main__":
parser = ArgumentParser()
parser.add_argument("infile", help="")
parser.add_argument("-c", "--column", default=2, type=int, help="")
opts = parser.parse_args()
data = myutil.get_column_as_list( opts.infile, opts.column, "float" )
zscore = myutil.cal_zscores_from_list( data )
lines = open( opts.infile, "r" ).readlines()
assert len(lines) == len(data), "discrepancy in lines( %s ) and data( %s )" %(len(lines), len(data))
for idx in range( 0, len(lines) ):
print lines[idx].strip(), zscore[idx]
parser = ArgumentParser()
parser.add_argument("-f", "--files", nargs="+", required=True, help="structures you would like to superimpose to")
parser.add_argument("-c", "--column", required=True, type=int, help="column")
parser.add_argument(
"-m",
"--stats_funcnames",
nargs="+",
default=["mean", "median", "min", "max", "sum"],
help="stats you would like to calculate: mean/ median/ min/ max/ percentile/",
)
parser.add_argument("-l", "--label", action="store_true")
# parser.add_argument("-p", "--percentile", help="what percentile?")
opts = parser.parse_args()
stdout.write("#filename %s\n" % " ".join(opts.stats_funcnames))
# stdout.write( opts.tag )
for f in opts.files:
assert exists(f)
stdout.write(f)
nums = myutil.floatList(myutil.get_column_as_list(f, opts.column))
for stats_funcname in opts.stats_funcnames:
if nums:
stdout.write(" %5.1f" % cal_stats(stats_funcname, nums))
else:
stdout.write(" NA")
if opts.label:
stdout.write(" (%s)" % stats_funcname)
stdout.write("\n")