# HG changeset patch
# User francesco_lapi
# Date 1722438102 0
# Node ID 4f40b848d55a2e131d3927ad01c8a33146241b9a
# Parent  f2be75f8ef6012e4de5380464bc31f3ed2e774d4
Uploaded

diff -r f2be75f8ef60 -r 4f40b848d55a marea_2/flux_to_map.py
--- a/marea_2/flux_to_map.py	Wed Jul 31 14:57:06 2024 +0000
+++ b/marea_2/flux_to_map.py	Wed Jul 31 15:01:42 2024 +0000
@@ -738,10 +738,11 @@
         try: 
             # fallthrough is intended, regular scores need to be computed when tips aren't net but RAS datasets aren't used
             p_value, z_score = computePValue(l1, l2)
-            avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2))
+            avg1 = sum(l1) / len(l1)
+            avg2 = sum(l2) / len(l2)
+            avg = fold_change(avg1, avg2)
             if not isinstance(z_score, str) and max_z_score < abs(z_score): max_z_score = abs(z_score)
-            tmp[reactId] = [float(p_value), avg, z_score]
-        
+            tmp[reactId] = [float(p_value), avg, z_score, avg1, avg2]
         except (TypeError, ZeroDivisionError): continue
     
     return tmp, max_z_score