Hi,
A tree has n terminal nodes (aka tips) and m internal nodes (aka nodes simply).
In the edge matrix, the tips are numbered 1:n and the nodes are numbered
(n+1):(n+m) (same than n+1:m). n and m can be found with:
n <- Ntip(tree) # or length(tree$tip.label)
m <- Nnode(tree) # or tree$Nnode
For `ape::dist.nodes()`, how can one match the output matrix rows/columns with
the node IDs in the tree (eg., the tip.labels)? I cannot just use
`ape::cophenetic()` in my particular situation. The docs for
`ape::dist.nodes()` state:
```
… in the case of dist.nodes, the numbers of the tips and
