Mercurial > repos > bimib > cobraxy

changeset 509:5956dcf94277 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author francesco_lapi
date Wed, 01 Oct 2025 15:34:21 +0000
parents ca98c149ec61
children
files COBRAxy/ras_generator.py
diffstat 1 files changed, 1 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/COBRAxy/ras_generator.py	Wed Oct 01 14:21:26 2025 +0000
+++ b/COBRAxy/ras_generator.py	Wed Oct 01 15:34:21 2025 +0000
@@ -191,6 +191,7 @@
     for rule, reaction_ids in dict_rule_reactions.items():
         if len(rule) != 0:
             # there is one gene at least in the formula
+            rule = rule.replace("-","_")
             rule_split = rule.split()
             rule_split_elements = list(set(rule_split))                                # genes in formula
             
@@ -220,8 +221,6 @@
                             else:
                                 # complex expression (e.g. A or (B and C))
                                 data = count_df_filtered.loc[genes_in_count_matrix]  # dataframe of genes in the GPRs
-
-                                rule = rule.replace("-", "_")
                                 tree = ast.parse(rule, mode="eval").body
                                 values_by_cell = [dict(zip(data.index, data[col].values)) for col in data.columns]
                                 for j, values in enumerate(values_by_cell):