fix IR

dlcgold · dlcgold · commit cd3243aa80f5 · 2024-08-27T08:34:53.000+02:00
diff --git a/exps/2-dm-real/workflow/scripts/parse_res.py b/exps/2-dm-real/workflow/scripts/parse_res.py
@@ -86,7 +86,7 @@ def parse_pantas(path):
         if abs(e.dpsi) < FILTER:
             continue
         if e.etype == "IR":
-            e.dpsi = 1 - e.dpsi
+            e.dpsi = -e.dpsi # Pantas alternative isoform = rMATS canonical isoform
         event_pantas[e.etype].append(e)
     return event_pantas
 
diff --git a/scripts/remap.py b/scripts/remap.py
@@ -187,7 +187,7 @@ def main(args):
                 junction2_nodes = "."
                 W1, W2 = W2, W1
                 psi_c1, psi_c2 = psi_c2, psi_c1
-                dpsi = 1.0-float(dpsi)
+                dpsi = -float(dpsi)
 
         htype = "reference"
         if annotation_type == "annotated":
