Prune Trees To Overlapping Taxa
I have some phylogenetic trees with partial overlapping taxa. I would like to prune these trees, say 2 of them, to their overlapping taxa in order to be able to compare them. Is there any built in function for this in R, Python, etc (or I need to invent the wheel)?
Thanks!
2 answers
A Python approach using ETE:
from ete3 import Tree
newick1 = "((A, B, (C, D)), (E, F));"
newick2 = "((A, B, (C, G)), (E, T));"
t1 = Tree(newick1)
t2 = Tree(newick2)
common_leaves = set(t1.get_leaf_names()) & set(t2.get_leaf_names())
t1.prune(common_leaves)
t2.prune(common_leaves)
t1.write()
Out[16]: '((A:1,B:1,C:1)1:1,E:1);'
t2.write()
Out[17]: '((A:1,B:1,C:1)1:1,E:1);'
In R: match vectors of names between the two trees, then use drop.tip() inape.
• 0 views
•
link
updated
by
Ram
•
written
by
Brice Sarver
•
written
by
Brice Sarver
Log in to answer this question.
More posts like this
-
Force same scale on different phylogenetic trees
written by hakon.kaspersen •
-
Drop tip from multiPhylo object in R
written by phylofun
-
prunning posterior distribution of trees in R
written by phylofun
-
Comparative measure of phylogenetic tree resolution
written by jnf3769 •
-
annotating Multiple Phylogenetic trees on the basis of color
written by Sourabh J •
-
Resolving multifurcating at root?
written by rezam •
-
Pruning phylogenetics trees to reduce complexity
written by tim_mci •
-
construct a phylogenetic tree from a set of local trees
written by sam •
-
Program for searching for a specific partition within a set of phylogenetic trees
written by rwn
-
Effective Scoring Function To Compare Multiple Sequence Alignments Or Phylogenetic Trees
written by Pappu
