We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 13272fb commit 23f935fCopy full SHA for 23f935f
treetime/treetime.py
@@ -816,9 +816,9 @@ def generate_subtree(self, parent):
816
# branches without mutations are ready to coalesce -- others have to mutate first
817
ready_to_coalesce = [b for b in branches_alive if mutations_per_branch.get(b.name,0)==0]
818
if hasattr(self, 'merger_model') and (self.merger_model is not None):
819
- coalescent_rate = 0.5*len(ready_to_coalesce)*dummy_coalescent_rate + mutation_rate
820
- else:
821
coalescent_rate = self.merger_model.branch_merger_rate(t) + mutation_rate
+ else:
+ coalescent_rate = 0.5*len(ready_to_coalesce)*dummy_coalescent_rate + mutation_rate
822
823
total_mutations = np.sum([mutations_per_branch.get(b.name,0) for b in branches_alive])
824
n_branches_w_mutations = len(branches_alive) - len(ready_to_coalesce)
0 commit comments