fixing empty intervals bug

popsim-consortium · Jul 30, 2024 · 0bdce29 · 0bdce29
1 parent 5083613
commit 0bdce29
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/workflows/ts2fs.py b/workflows/ts2fs.py
@@ -43,7 +43,7 @@ def get_expected_netural_subs(ts, mutation_rate):
             if mtypes['is_neutral']:
                 annot_intervals = np.array(dfe['intervals'])
                 if annot_intervals.shape[0] == 0:
-                    np.array([[0, sequence_length]])
+                    annot_intervals = np.array([[0, sequence_length]])
                 #assumes netural is first proportion in list
                 proportion = dfe['proportions'][0]
                 ratemap = mask_to_ratemap(annot_intervals, sequence_length, mutation_rate, proportion)