Internal function needed to simulate the ranked tree shape model
lambda.epsilon(epsilon, beta, n, upper = 10000)Minimum size of unsampled splits (see appendix 1)
Imbalance index
Tip number
Upper bound from the integral (see lambda_n expression, default to 10000)
\(\lambda_\epsilon = 2 \int_(\epsilon,\infty) exp(-(\beta+n+1)x)(1-exp(-x))^\beta dx\)
Maliet O., Gascuel F., Lambert A. (2018) Ranked tree shapes, non-random extinctions and the loss of phylogenetic diversity, bioRxiv 224295, doi: https://doi.org/10.1101/224295