Question: How do I show bootstrap values in the phylogenetic tree using R?
gravatar for lhernanj
2.5 years ago by
lhernanj0 wrote:

Hello guys, I am having a difficult time figuring out how to show the bootstrap values in my tree. I am looking for an easy answer because I am not expert in R. I just done a RaxML tree and exported my results to R. I got constructed a fine figure, but I don't know how to show bootstrap values in that. Please look at my R code:

raxml.tree <- read.tree(file.path("RAxML_bipartitionsBranchLabels.job3"))
raxml.bootstrap <- read.tree(file.path("RaxML_bootstrap.job2"))

raxml<- rep("black", length(raxml.tree$tip.label))
categorias<- c("A5190.1","Z21068.1")

for(i in 1:length(categorias)){
  raxml[grep(categorias[i], raxml.tree$tip.label)] <- colorcategorias[i]
edgecol<- rep("black", nrow(raxml.tree$edge))
edgecol[23:61]<- "blue" 
edgecol[73:75]<- "blue"

plot(unroot(raxml.tree),title("A", line=-2, adj=0.3), use.edge.length =TRUE, tip.color= raxml, type="unrooted", show.tip.label = TRUE, cex= 0.7, no.margin= TRUE, lab4ut = "axial", edge.width = 2, edge.color = edgecol, show.node.label=TRUE)

as you can see I open my bootstrap file but I am not sure about how should I manipulate it.

I would thank some of help.

phytools phylogenetic tree R • 1.5k views
ADD COMMENTlink modified 2.1 years ago • written 2.5 years ago by lhernanj0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2286 users visited in the last hour