在 R 中，有没有一种方法可以根据数据框中的因子向量为 phytools 树着色？

Question

我有一棵树，上面有 18 个提示，看起来像这样；

library(phytools)
# Create tree
sim.tree<-pbtree(n=18)
plot(sim.tree)

然后我有一个带有提示名称的数据框和一个我想用来给树着色的因子向量，比如；

df<-data.frame(tip = paste0("t", 1:18),
               vector.to.color.with = as.factor(c("<10", "10-20", "10-20", "10-20", "NA", "10-20", 
                                   "10-20", "10-20", "20-35", "<10", "10-20", "<10", 
                                   "35", "20-35", "<10", "NA", "10-20", "<10")))

我尝试使用tip.color()和plotBranchbyTrait()无济于事。 我知道这可能一直都在做，我错过了一些简单的事情。 任何帮助是极大的赞赏

Answer 1

检查?plot.phylo

library(phytools)
# Create tree
df<-data.frame(tip = paste0("t", 1:18),
               vector.to.color.with = as.factor(c("<10", "10-20", "10-20", "10-20", "NA", "10-20", 
                                                  "10-20", "10-20", "20-35", "<10", "10-20", "<10", 
                                                  "35", "20-35", "<10", "NA", "10-20", "<10")))
set.seed(123)
sim.tree<-pbtree(n=18)
cols <- as.integer(df$vector.to.color.with[match(sim.tree$tip.label,df$tip)])
plot(sim.tree, tip.color = cols)