# HG changeset patch # User francesco_lapi # Date 1721997808 0 # Node ID 2c75a3b96bc0ebd60ffd159f26514a8d293fe4ad # Parent c9908e5f107fece02e673935cc1f0d12c252a4f2 Uploaded diff -r c9908e5f107f -r 2c75a3b96bc0 marea_2/marea.py --- a/marea_2/marea.py Fri Jul 26 10:59:40 2024 +0000 +++ b/marea_2/marea.py Fri Jul 26 12:43:28 2024 +0000 @@ -732,16 +732,11 @@ reactDir = ReactionDirection.fromReactionId(reactId) # Net score is computed only for reversible reactions when user wants it on arrow tips or when RAS datasets aren't used if (ARGS.net or not ARGS.using_RAS) and reactDir is not ReactionDirection.Unknown: - try: - position = ids.index(reactId[:-1] + ('B' if reactDir is ReactionDirection.Direct else 'F')) - nets1 = np.subtract(l1, dataset1Data[position]) - nets2 = np.subtract(l2, dataset2Data[position]) - except ValueError: - if ARGS.net: - continue # we look for the complementary id, if not found we skip - else: - nets1 = l1 - nets2 = l2 + try: position = ids.index(reactId[:-1] + ('B' if reactDir is ReactionDirection.Direct else 'F')) + except ValueError: continue # we look for the complementary id, if not found we skip + + nets1 = np.subtract(l1, dataset1Data[position]) + nets2 = np.subtract(l2, dataset2Data[position]) p_value, z_score = computePValue(nets1, nets2) avg1 = sum(nets1) / len(nets1)