comparison COBRAxy/marea.py @ 277:1bf66a30ab23 draft

Uploaded
author francesco_lapi
date Mon, 12 May 2025 15:36:25 +0000
parents 2748cef7fa12
children 5e8f1ab4242f
comparison
equal deleted inserted replaced
276:2748cef7fa12 277:1bf66a30ab23
765 765
766 # fallthrough is intended, regular scores need to be computed when tips aren't net but RAS datasets aren't used 766 # fallthrough is intended, regular scores need to be computed when tips aren't net but RAS datasets aren't used
767 p_value, z_score = computePValue(l1, l2) 767 p_value, z_score = computePValue(l1, l2)
768 avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2)) 768 avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2))
769 if not isinstance(z_score, str) and max_z_score < abs(z_score): max_z_score = abs(z_score) 769 if not isinstance(z_score, str) and max_z_score < abs(z_score): max_z_score = abs(z_score)
770 tmp[reactId] = [float(p_value), avg, z_score] 770 tmp[reactId] = [float(p_value), avg, z_score, sum(l1) / len(l1), sum(l2) / len(l2)]
771 771
772 except (TypeError, ZeroDivisionError): continue 772 except (TypeError, ZeroDivisionError): continue
773 773
774 return tmp, max_z_score 774 return tmp, max_z_score
775 775