Hi Yan.

In case it is useful I just pushed an update to phytools in which I have added an as.multiPhylo method for objects of class "phylo". This will also return an object of class "multiPhylo" if that is already the object class. This means that:


obj<-as.multiPhylo(read.nexus(...))

in which ... are the arguments to read.nexus, will return an object of class "multiPhylo" with length 1 if the file contains only one tree - but will also return an object of class "multiPhylo" if the file has 2 or more trees.

To get this one needs to install phytools from GitHub using devtools:

library(devtools)
install_github("liamrevell/phytools")

All the best, Liam

Liam J. Revell, Associate Professor of Biology
University of Massachusetts Boston
web: http://faculty.umb.edu/liam.revell/
email: liam.rev...@umb.edu
blog: http://blog.phytools.org

On 12/27/2016 5:11 PM, Yan Wong wrote:

On 27 Dec 2016, at 22:08, Emmanuel Paradis <emmanuel.para...@ird.fr> wrote:

Hi Yan,

Yes, you are right: this modification can be done. In the meantime, you can fix 
the code with:

fix(read.nexus)

Find this line (#117):

   if (Ntree == 1) {

and change it to:

   if (FALSE) {

save and close. Tell me if you still have problem.

Yes, that works fine. Thanks a lot.

Yan
_______________________________________________
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/


_______________________________________________
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