pbd_bootstrap function

Bootstrap analysis under protracted birth-death model of diversification