N <- length(phy$tip.label)
bt <- sort(branching.times(phy))
g <- rev(c(bt[1], diff(bt))) # internode intervals are from past to present
N <- length(phy$tip.label)
bt <- sort(branching.times(phy))
g <- rev(c(bt[1], diff(bt))) # internode intervals are from past to present