time <- sort(branching.times(phy), decreasing = TRUE)
N <- 1:(length(time) + 1)
plot.default(-c(time, 0), N, xlab = xlab, ylab = ylab,
time <- sort(branching.times(phy), decreasing = TRUE)
N <- 1:(length(time) + 1)
plot.default(-c(time, 0), N, xlab = xlab, ylab = ylab,