Skip to content

Commit 2e4f09a

Browse files
Minor tweak
1 parent 694812d commit 2e4f09a

File tree

1 file changed

+3
-3
lines changed

1 file changed

+3
-3
lines changed

sc2ts/inference.py

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -403,9 +403,9 @@ def match_samples(
403403
num_threads=None,
404404
):
405405
run_batch = samples
406-
406+
407407
mu = 0.125 ## FIXME
408-
for k in range(num_mismatches):
408+
for k in range(int(num_mismatches)):
409409
# To catch k mismatches we need a likelihood threshold of mu**k
410410
likelihood_threshold = mu**k - 1e-15
411411
logger.info(f"Running match={k} batch of {len(run_batch)} at threshold={likelihood_threshold}")
@@ -808,7 +808,7 @@ def solve_num_mismatches(k, num_sites, mu=0.125):
808808
# values of k <= 1 are not relevant for SC2 and lead to awkward corner cases
809809
assert k > 1
810810

811-
denom = (1 - mu) ** k
811+
denom = (1 - mu) ** k
812812
r = mu**k / denom
813813

814814
# Add a little bit of extra mass for recombination so that we deterministically

0 commit comments

Comments
 (0)