# HG changeset patch
# User francesco_lapi
# Date 1747064185 0
# Node ID 1bf66a30ab239717b8700f51433db6ad1cf65143
# Parent  2748cef7fa1286a4fd6c58a153f5cf68ea7a7aa3
Uploaded

diff -r 2748cef7fa12 -r 1bf66a30ab23 COBRAxy/marea.py
--- a/COBRAxy/marea.py	Mon May 05 10:33:30 2025 +0000
+++ b/COBRAxy/marea.py	Mon May 12 15:36:25 2025 +0000
@@ -767,7 +767,7 @@
             p_value, z_score = computePValue(l1, l2)
             avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2))
             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, sum(l1) / len(l1), sum(l2) / len(l2)]
         
         except (TypeError, ZeroDivisionError): continue
     
@@ -911,4 +911,4 @@
     print('Execution succeeded')
 ###############################################################################
 if __name__ == "__main__":
-    main()
\ No newline at end of file
+    main()