Hi,

To get a tree with a basal multichotomy as a rooted tree, you have to set its root edge:


R> tr <- read.tree(text = "(A,B,(C,D));")
R> is.rooted(tr)
[1] FALSE
R> tr$root.edge <- 0
R> is.rooted(tr)
[1] TRUE

This could be done directly in the Newick string:

R> tr <- read.tree(text = "(A,B,(C,D)):0;")
R> is.rooted(tr)
[1] TRUE

Best,

Emmanuel

Le 19/01/2018 à 23:08, Yan Wong a écrit :
Hi

I have a set of rooted trees, some of which have polytomies at the root, like

(A,B,(C,D));
((A,B),(C,D));

These are not treated as rooted by default, so I have tried to root them using

root(tree, node = Ntip(phy) + 1, resolve.root = TRUE), but this gives an error, 
explained here:

https://www.mail-archive.com/r-sig-phylo@r-project.org/msg03805.html

Note that in my case there is no particular reason to specify either A, B, or 
(C,D) as an outgroup

I should explain that what I am trying to do is to calculate distance metrics 
between these trees (using metrics that are valid with polytomies, such as the 
KC metric), such that the pairwise distances between (A,B,(C,D)) and 
(B,A,(C,D)) and ((C,D), B, A) etc are all 0, but the distance between, say, 
(A,B,(C,D)) and ((A,B),(C,D)) or perhaps (A,C,(B,D)); is not zero.

So how can I specify a polytomy at the root, with the tree format used in R?

Cheers

Yan Wong
Oxford
_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/


Pour nous remonter une erreur de filtrage, veuillez vous rendre ici : 
http://f.security-mail.net/304anPlLmqF



_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/

Reply via email to