From 881100c4ebec9638a497b0b0045a4f8360e52a2f Mon Sep 17 00:00:00 2001 From: Yasamin Tabatabaee Date: Thu, 26 Oct 2023 20:30:48 -0500 Subject: [PATCH] changing mu for terminal branches --- castles.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/castles.py b/castles.py index bfe13e1..eb7eb89 100644 --- a/castles.py +++ b/castles.py @@ -243,7 +243,7 @@ def castles(st, gts, tns): # compute internal branch (Table S3 in paper) delta = safe_div(lm_i - ln_i, ln_i) if lm_i > ln_i else 1e-03 l_est = 1/6 * (3 * delta + np.sqrt(3 * delta * (4 + 3 * delta))) * ln_i if ln_i > 0 else 1e-06 - mu1_est = ln_i #l_est / d_est + mu1_est = l_est / d_est # compute terminal branches (Table S3 in paper, unbalanced) l_a_est = ln_a + (mu1_est * (d_est - p_est) + (lm_a - ln_a) * (