There was a bug when picking the outgroup for summarisation of trees,…

… such tip was hard-coded as t1 and therefore the bug went silent when using simulated trees from ape. However on trees with tip labels other than tN, it failed as t1 was not found to be a tip. Raising patch number.

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: tbea
 Title: Tools for Pre- and Post-processing in Bayesian Evolutionary Analyses
-Version: 1.1.1
+Version: 1.1.2
 Authors@R: c(person("Gustavo A.", "Ballen", email = "gaballench@gmail.com", role = c("aut", "cre")), 
  person("Sandra", "Reinales", email = "spreinalesl@gmail.com", role = c("aut")))
 Description: Package for bayesian inference in phylogenetics and evolution.
@@ -16,7 +16,7 @@ Imports: ape, Rfit, Hmisc, boot, graphics, stats, utils
 License: file LICENSE
 Encoding: UTF-8
 LazyData: true
-RoxygenNote: 7.1.2
+RoxygenNote: 7.2.3
 URL: https://github.com/gaballench/tbea
 BugReports: https://github.com/gaballench/tbea/issues
 Suggests: 

R/summaryBrlen.R

@@ -36,10 +36,12 @@ summaryBrlen <- function(mphy, method) {
  if (inherits(mphy, "multiPhylo")) {
  mphy <- list(mphy)
  }
+ # pick an arbitrary branch for rooting
+ outgroup <- mphy[[1]][[1]]$tip.label[1]
  summary_trees <- vector("list", length=length(mphy))
  iter <- 1
  for (i in mphy) {
- rooted_tr <- root(i, outgroup="t1", resolve.root=TRUE)
+ rooted_tr <- root(i, outgroup=outgroup, resolve.root=TRUE)
  for (k in seq_along(rooted_tr)) {
  rooted_tr[[k]] <- makeNodeLabel(rooted_tr[[k]], method="md5sum")
  }