Simulates the time before the splitting of the marks into two different fragments, and the size of the fragment at this time.
simulate.Tau.X(epsilon, x, alpha, beta, n, ab = FALSE, eta = 1, x.ab = 1, lambda = NULL)
Minimum size of unsampled splits (see appendix 1)
Initial fragment size
Clade age-richness index
Imbalance index
Number of marks
Boolean, should the abundances be computed as well? (default to FALSE)
Clade abundance-richness index (if ab)
Initial abundance
Optional, vector of lambda_epsilon returned by the function lambda_N. If set to NULL (the default) it will be computed inside the function
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