Merge pull request #717 from macs3-project/fix/macs3/pvalue_stat

fix pvalue_stat issue that will cause nan result

Author: taoliu

MACS3/Signal/ScoreTrack.py

@@ -1,6 +1,6 @@
 # cython: language_level=3
 # cython: profile=True
-# Time-stamp: <2025-11-16 13:00:03 Tao Liu>
+# Time-stamp: <2025-11-20 10:58:03 Tao Liu>
 
 """Scoring utilities for MACS3 signal tracks and peak callers.
 
@@ -414,6 +414,7 @@ def compute_pvalue(self):
         v: cnp.ndarray
         pos: cnp.ndarray
         ln: cython.long
+        tmp_l: cython.long
         i: cython.long
         prev_pos: cython.long
         chrom: bytes
@@ -429,10 +430,11 @@ def compute_pvalue(self):
                 v[i] = get_pscore(cython.cast(cython.int,
                                               (p[i] + self.pseudocount)),
                                   c[i] + self.pseudocount)
+                tmp_l = pos[i] - prev_pos
                 try:
-                    self.pvalue_stat[v[i]] += pos[i] - prev_pos
+                    self.pvalue_stat[v[i]] += tmp_l
                 except Exception:
-                    self.pvalue_stat[v[i]] = pos[i] - prev_pos
+                    self.pvalue_stat[v[i]] = tmp_l
                 prev_pos = pos[i]
 
         self.scoring_method = ord('p')