R ape: find the branch length of a particular tip label
Hello-
I have a tree in R:
library(ape)
tree = read.tree(file)
and I would like to associate each tip label in tree$tip.label with the length of the branch immediately leading up to that tip. I know that I can pull out the branch lengths with tree$edge.length, but does anyone know how to associate these lengths with tip labels?
Thanks, Alissa
• 8,807 views
•
link
updated
by
Hend •
•
written
by
alissamariewilliams2015 •
•
written
by
alissamariewilliams2015 •
1 answer
Log in to answer this question.
More posts like this
-
Colour tip labels in a tree based on wildcard/expression using ggtree
written by mrmrwinter •
-
Plot phylogenetic tree with all ancestor nodes
written by AfinaM •
-
plotting trees in R: node labels from nexus tree files
written by jmah •
-
Adding branch lengths to a topology
written by treesandthings •
-
Combine node/branch labels from multiple trees into one
written by al-ash
-
Figtree - Make branches the same colour as their tip labels
written by wl284 •
-
The fonts and size of a tree displayed by ggtree
written by sunyeping •
-
Phylogenetics in R: obtaining internal branch lengths
written by spiral01
-
Plotting traits on phylogenetic tree in R - issue with aligning traits
written by ch_at •
-
Cut a tree based on time t
written by tutoring891 •
What do you mean by "associate these lengths with tip labels"? If it's for plotting, you should be able to do something like:
Or if it's about how to retrieve which value in tree$edge.length corresponds to the which tree leave, check tree$edge. It's a two column matrix where each row is an edge whose length is in the same position in the tree$edge.length vector. The first column gives the parent node and the second gives the child node. Nodes are numbered such that the leaves are numbered first then the internal nodes so n leaves would be numbered 1:n and the internal nodes would have id>n. So
tree$edge.length[tree$edge[,2] <= Ntip(tree)]should give you the lengths of all branches ending at a leave.Thank you for your response! I have plotted the trees with tip labels and branch lengths before, so my question here is more related to your latter answer. I see that I can get the lengths of all branches ending at a leaf, but is there a way to also pull out the name on each of those leaves? As in, can I use a built-in function to retrieve the name associated with a particular branch length? If that makes sense. I want something conceptually like:
where the tip names are in the same order as their respective branches, as outputted with
What about tree$tip.label? In this vector the labels are in the order of the corresponding node ID, i.e. tree$tip.label[2] is the label for node with ID 2.