fit_hbd_psr_parametric function

Fit parameterized pulled speciation rates of birth-death models.