# HG changeset patch
# User bimib
# Date 1721739259 0
# Node ID 11f9ed694ad42c69997e139651a388cc33dd09fd
# Parent  50d883e234779b5c07e6701398e7c60de40255b1
Uploaded

diff -r 50d883e23477 -r 11f9ed694ad4 marea_2/marea.py
--- a/marea_2/marea.py	Mon Jul 22 11:28:36 2024 +0000
+++ b/marea_2/marea.py	Tue Jul 23 12:54:19 2024 +0000
@@ -218,7 +218,7 @@
     elif avg2 == 0:
         return 'INF'
     else:
-        return math.log(avg1 / avg2, 2)
+        return (avg1 - avg2) / abs(avg2)
     
 def fix_style(l :str, col :Optional[str], width :str, dash :str) -> str:
     """
@@ -722,7 +722,7 @@
                 p_value = computePValue(nets1, nets2)
                 avg1 = sum(nets1)   / len(nets1)
                 avg2 = sum(nets2)   / len(nets2)
-                net = (avg1 - avg2) / abs(avg2)
+                net = fold_change(avg1, avg2)
                 
                 if math.isnan(net): continue
                 tmp[reactId[:-1] + "RV"] = [p_value, net, avg1, avg2]