R/intnode.times.R

In sebastianduchene/NELSI: NELSI: Nucleotide EvoLution Simulator

intnode.times <-
function(phylo){
    di.phylo <- dist.nodes(phylo)
    root.phylo <- phylo$edge[, 1][!(phylo$edge[, 1] %in% phylo$edge[, 2])][1]
    phylo.depth <- max(di.phylo[as.numeric(colnames(di.phylo)) == root.phylo, ])
    node.times <- phylo.depth - di.phylo[as.numeric(rownames(di.phylo)) == root.phylo, (length(phylo$tip.label) + 1):nrow(di.phylo)]
    return(node.times)
}