How does ggtree interprets branch support value (e.g bootstrap value) given as node label?
1
3
Entering edit mode
6.4 years ago
Chirag Parsania ★ 2.0k

Hi,

I came across the paper describing issues interpreting phylogenetic tree due to branch support given as node label in newick format.

Let me describe the problem very briefly, which is explained in depth in the paper given above. Many tree output programmes gives branch support values (e.g bootstrap value) as node label. However, many tree manipulation softwares interpret node labels as node labels and not branch support values. Therefore, while rerooting the tree it misplace the branch support values, which leads to different tree interpretation.

I wonder, how ggtree addresses this issue ? below, I explained the situation with example

library("ggtree")
## random tree with random bootstap values given as node label in newick format 
sample_tree <- read.tree(text = "((C,D)90,(A,(B,X)95)100,E)85;")

## reroot from E 
reroot_e <- root(phy = sample_tree , outgroup = "E") 
ggtree(reroot_e)+ geom_tiplab(size = 7, col ="blue") +geom_nodelab(size = 7, col= "red") + theme_tree2() +  xlim(0,4.5)

root_e

## reroot the tree from X 
re_root_x <- root(phy = sample_tree , outgroup = "X")
ggtree(re_root_x)+ geom_tiplab(size = 7, col ="blue") +geom_nodelab(size = 7, col= "red") + theme_tree2() + xlim(0,4.5)

root_x

I wonder,are these trees have same bootstrap interpretation or they are different ??

@Guangchuang Yu : you might be the best person to answer this. Can someone tag him on Biostar?? Any help would be highly appreciated

ggtree phylogenetic newick R • 6.2k views
ADD COMMENT
0
Entering edit mode

You have already tagged Guangchuang in this post. Unfortunately he does not appear to have logged into Biostars for a month+.

Can you edit the title to make it concise?

ADD REPLY
2
Entering edit mode
6.4 years ago
Guangchuang Yu ★ 2.6k

This is an issue of re-root, you should ask the author of ape for the ape::root function.

ADD COMMENT
0
Entering edit mode

Hi Guangchuang,

Thanks for your reply. I didn't realise that root belongs ape. Anyway, I checked the documentation of ape::root. In fact after the critical review rerooting function of ape have been updated.

ADD REPLY

Login before adding your answer.

Traffic: 2133 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6