# HG changeset patch
# User francesco_lapi
# Date 1746440123 0
# Node ID 99703bcaa2901a1251231c6b598db007be2437d9
# Parent  655e33ca9225d57f3812be88f30220d10f01102e
Uploaded

diff -r 655e33ca9225 -r 99703bcaa290 COBRAxy/marea.py
--- a/COBRAxy/marea.py	Mon Mar 10 15:54:03 2025 +0000
+++ b/COBRAxy/marea.py	Mon May 05 10:15:23 2025 +0000
@@ -761,6 +761,7 @@
             # fallthrough is intended, regular scores need to be computed when tips aren't net but RAS datasets aren't used
             p_value, z_score = computePValue(l1, l2)
             avg = fold_change(sum(l1) / len(l1), sum(l2) / len(l2))
+            if math.isnan(p_value) or (isinstance(avg, float) and math.isnan(avg)): continue
             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]