Compare series

scomp.py

from numpy import genfromtxt, average, median, std, percentile
from sys import argv, exit

if len(argv) < 3:
    print 'Usage: {0} <series1.txt> <series2.txt> [map lambda]'.format(argv[0])
    exit(1)

lam = None
if len(argv) == 4:
    lam = eval(argv[3])

data = [genfromtxt(argv[1]), genfromtxt(argv[2])]
if callable(lam):
    data = map(lam, data)

fmt = '{0:.4} vs {1:.4} ({2:+.3}%)'

def pct(bef, aft):
    return (aft-bef)/bef*100

def summary(series):
    [af, ag] = map(average, series)
    [mf, mg] = map(median, series)
    [sf, sg] = map(std, series)
    print 'Average    : '+ fmt.format(af, ag, pct(af, ag))
    print 'Median     : '+ fmt.format(mf, mg, pct(mf, mg))
    print 'Std dev    : '+ fmt.format(sf, sg, pct(sf, sg))
    print 'Percentiles:'
    for p in [5, 10, 90, 95]:
        [pf, pg] = map(lambda x:percentile(x, p), series)
        print ('    {3}% : '+ fmt).format(pf, pg, pct(pf, pg), p)

summary(data)