fix: make sure joint_inverse_cdf is not computed for bad branches
@@ -482,6 +482,7 @@ class ClockTree(TreeAnc):
node.joint_pos_Lx = res
node.joint_pos_Cx = res_t
+
# construct the inverse cumulant distribution for the branch number estimates
from scipy.interpolate import interp1d
if node.date_constraint is not None and node.date_constraint.is_delta:
