Mercurial > repos > bimib > cobraxy

changeset 197:1306b3543e57 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author francesco_lapi
date Thu, 21 Nov 2024 10:43:48 +0000
parents 07539860dbca
children 1dda40a60037
files COBRAxy/flux_to_map.py
diffstat 1 files changed, 5 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/COBRAxy/flux_to_map.py	Thu Nov 21 09:42:58 2024 +0000
+++ b/COBRAxy/flux_to_map.py	Thu Nov 21 10:43:48 2024 +0000
@@ -497,8 +497,8 @@
             except ZeroDivisionError: pass
         # TODO CHECK RV
         #if not reactionId.endswith("_RV"): # RV stands for reversible reactions
-        #    Arrow(width, ArrowColor.fromFoldChangeSign(foldChange)).styleReactionElements(metabMap, reactionId)
-        #    continue
+        #   Arrow(width, ArrowColor.fromFoldChangeSign(foldChange)).styleReactionElements(metabMap, reactionId)
+        #   continue
         
         #reactionId = reactionId[:-3] # Remove "_RV"
         
@@ -726,9 +726,10 @@
             p_value, z_score = computePValue(l1, l2)
             avg1 = sum(l1) / len(l1)
             avg2 = sum(l2) / len(l2)
-            avg = fold_change(avg1, avg2)
+            f_c = 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, avg1, avg2]
+            print(reactId, float(p_value), f_c, z_score, avg1, avg2)
+            tmp[reactId] = [float(p_value), f_c, z_score, avg1, avg2]
         except (TypeError, ZeroDivisionError): continue
     
     return tmp, max_z_score