Loglikelihood for diversity-dependent diversification models with multiple decoupling (rate shift) events