A = [9,7,4,6,3] B = [2,3,3,8,4] import scipy.stats as st def avg(lst): return float(sum(lst)) / len(lst) print("A:", avg(A)) print("B:", avg(B)) print("Speedup: %.2f" % ((avg(A) / avg(B)))) stddev = st.nanstd print("A: %.2f" % stddev(A)) print("B: %.2f" % stddev(B)) print(avg(A)-stddev(A), "<", avg(B)+stddev(B), " (overlap!)") pval = st.ttest_ind(A,B)[1] print("p-value: %.0f%%" % (pval*100)) pval = st.ttest_ind(A*3,B*3)[1] print("p-value: %.2f%%" % (pval*100))