# HG changeset patch # User bimib # Date 1587676661 14400 # Node ID 19d704e249778933fa3a41f46bb989eb4bd868be # Parent aa09186cbc9fb68dd76208f86e0f0743d2df7a8f Uploaded diff -r aa09186cbc9f -r 19d704e24977 Marea/marea.py --- a/Marea/marea.py Thu Apr 23 17:03:48 2020 -0400 +++ b/Marea/marea.py Thu Apr 23 17:17:41 2020 -0400 @@ -685,11 +685,16 @@ tmp = {} count = 0 max_F_C = 0 + + primo = -1 - for l1, l2 in zip(rest, class_pat.get(single_cluster)): + for l1, l2 in zip(class_pat.get(single_cluster), rest): try: stat_D, p_value = st.ks_2samp(l1, l2) avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2)) + if primo == -1: + primo = 0 + print(avg) if not isinstance(avg, str): if max_F_C < abs(avg): max_F_C = abs(avg)