Mercurial > repos > bimib > marea

comparison Marea/marea.py @ 69:19d704e24977 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author bimib
date Thu, 23 Apr 2020 17:17:41 -0400
parents a61733753aec
children 84da46340c39
comparison
equal deleted inserted replaced
68:aa09186cbc9f 69:19d704e24977
683 rest = rest + i 683 rest = rest + i
684 684
685 tmp = {} 685 tmp = {}
686 count = 0 686 count = 0
687 max_F_C = 0 687 max_F_C = 0
688 688
689 for l1, l2 in zip(rest, class_pat.get(single_cluster)): 689 primo = -1
690
691 for l1, l2 in zip(class_pat.get(single_cluster), rest):
690 try: 692 try:
691 stat_D, p_value = st.ks_2samp(l1, l2) 693 stat_D, p_value = st.ks_2samp(l1, l2)
692 avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2)) 694 avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2))
695 if primo == -1:
696 primo = 0
697 print(avg)
693 if not isinstance(avg, str): 698 if not isinstance(avg, str):
694 if max_F_C < abs(avg): 699 if max_F_C < abs(avg):
695 max_F_C = abs(avg) 700 max_F_C = abs(avg)
696 tmp[ids[count]] = [float(p_value), avg] 701 tmp[ids[count]] = [float(p_value), avg]
697 count += 1 702 count += 1