prunning posterior distribution of trees in R
I have a posterior distribution of trees from a Bayesian analysis in nexus format. Is there a function in R to prune these trees to match my trait dataset? Something like the treedata function in the R package geiger seems to do what I want, but it does not handle multiple trees.
1 answer
keep.tip from ape and lapply should do the job:
library(ape)
trees <- .uncompressTipLabel(trees)
trees <- lapply(trees, keep.tip, trait_names)
trees
class(trees) <- "multiPhylo"
trees <- .compressTipLabel(trees)
Here trees is a multiPhylo object (you can read them in with read.nexus) and trait_names is a vector with the tip labels you want to keep.
Log in to answer this question.
More posts like this
-
Simulation of label-free bottom-up proteomics expression dataset
written by KABILAN •
-
How to export and read into R Splitstree6 stree6 files
written by hans •
-
Importing nexus-formatted tree (e.g. BEAST tree) into R as a phylo object with node labels
written by al-ash
-
Missing ',' in line when Biopython reads a nexus tree
written by Mordovian Guy •
-
plotting trees in R: node labels from nexus tree files
written by jmah •
-
Drop tip from multiPhylo object in R
written by phylofun
-
How do I plot the posterior probability distribution of a particular node from a stochastic charact…
written by Namenlos •
-
Rooting phylogenetic trees in R using Ape's root function: Trees don't seem to root.
written by spiral01
-
Prune Trees To Overlapping Taxa
written by rezam •
-
Concatenating Nexus Files (For Mrbayes, Etc)
written by Andre Elias
