Skip to content

Commit 10c2c60

Browse files
committed
update ratio inconsistent peaks
1 parent 3e3af60 commit 10c2c60

1 file changed

Lines changed: 5 additions & 5 deletions

File tree

  • src/metrics/ratio_inconsistent_peaks

src/metrics/ratio_inconsistent_peaks/script.py

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -202,12 +202,12 @@
202202

203203
# case 1 or 3 where we have consistent peaks in unintegrated
204204
if peaks_u_s1 == peaks_u_s2:
205-
if peaks_s1 != peaks_s2:
206-
n_case3 += 1
207-
case_details["case3"].append((donor, celltype, marker))
208-
else:
205+
if peaks_s1 == peaks_u_s1 and peaks_u_s2 == peaks_s2:
209206
n_case1 += 1
210207
case_details["case1"].append((donor, celltype, marker))
208+
else:
209+
n_case3 += 1
210+
case_details["case3"].append((donor, celltype, marker))
211211
else:
212212
print(
213213
"WARNING! Inconsistent peaks detected in unintegrated data (case 2 or 4). Skipping calculation",
@@ -232,7 +232,7 @@
232232
)
233233
metric_val = np.nan
234234
else:
235-
metric_val = n_case3 / (n_case1 + n_case3)
235+
metric_val = n_case1 / (n_case1 + n_case3)
236236

237237

238238
print("Write output AnnData to file", flush=True)

0 commit comments

Comments
 (0)