X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=R%2Froot.R;h=aadd6c15454454947c4e5fc1e2d3a3d0124ede26;hb=3ad385892d75db5c646c92f0f631ae9c5e3da4f6;hp=3abdf650e8ab06091ae20a75e91051b80ccc2b8b;hpb=507aa18c4e3f9312efcb07b90766df1158a39402;p=ape.git diff --git a/R/root.R b/R/root.R index 3abdf65..aadd6c1 100644 --- a/R/root.R +++ b/R/root.R @@ -1,4 +1,4 @@ -## root.R (2009-07-06) +## root.R (2009-09-09) ## Root of Phylogenetic Trees @@ -297,14 +297,19 @@ root <- function(phy, outgroup, node = NULL, resolve.root = FALSE) phy$edge[, 1] <- newNb[phy$edge[, 1]] if (!is.null(phy$node.label)) { + #browser() newNb <- newNb[-(1:n)] if (fuseRoot) { newNb <- newNb[-1] phy$node.label <- phy$node.label[-1] } phy$node.label <- phy$node.label[order(newNb)] - if (resolve.root) - phy$node.label <- c(phy$node.label[1], NA, phy$node.label[-1]) + if (resolve.root) { + phy$node.label[is.na(phy$node.label)] <- phy$node.label[1] + phy$node.label[1] <- NA + ##phy$node.label <- c(phy$node.label[1], NA, phy$node.label[-1]) + ##phy$node.label <- c("NA", phy$node.label) + } } phy }