@@ -84,6 +84,7 @@ export(groupClade)

 export(groupOTU)

 export(gzoom)

 export(inset)

+export(is.ggtree)

 export(mask)

 export(merge_tree)

 export(msaplot)

@@ -1,5 +1,7 @@

 CHANGES IN VERSION 1.7.4

 ------------------------

+ o is.ggtree function to test whether object is produced by ggtree <2016-12-06, Tue>

+ o now branch.length can set to feature available in phylo4d@data and yscale is supported for phylo4d object <2016-12-06, Tue>

  o bug fixed of rm.singleton.newick, remove singleton parent instead of singleton <2016-12-01, Thu>

  o reorder phylo to postorder before ladderrize <2016-11-28, Mon>

  o allow yscale to use data stored in phylo4d object <2016-11-24, Thu>

@@ -1,6 +1,6 @@

 ##' drawing phylogenetic tree from phylo object

 ##'

-##' 

+##'

 ##' @title ggtree

 ##' @param tr phylo object

 ##' @param mapping aes mapping

@@ -43,13 +43,13 @@ ggtree <- function(tr,

                    branch.length  = "branch.length",

                    ndigits        = NULL,

                    ...) {

-    

+

     layout %<>% match.arg(c("rectangular", "slanted", "fan", "circular", "radial", "unrooted"))

     if (is(tr, "r8s") && branch.length == "branch.length") {

         branch.length = "TREE"

     }

-    

+

     if(yscale != "none") {

         ## for 2d tree

         layout <- "slanted"

@@ -76,20 +76,22 @@ ggtree <- function(tr,

     } else {

         multiPhylo <- FALSE

     }

-    

+

     p <- p + geom_tree(layout=layout, multiPhylo=multiPhylo, ...)

     p <- p + theme_tree()

-    

+

     if (layout == "circular" || layout == "radial") {

         p <- layout_circular(p)

         ## refer to: https://github.com/GuangchuangYu/ggtree/issues/6

         ## and also have some space for tree scale (legend)

-        p <- p + ylim(0, NA) 

+        p <- p + ylim(0, NA)

     } else if (layout == "fan") {

         p <- layout_fan(p, open.angle)

     }

+    class(p) <- c("ggtree", class(p))

+

     return(p)

 }

@@ -394,8 +394,11 @@ fortify.phylo4 <- function(model, data, layout="rectangular", yscale="none",

 ##' @method fortify phylo4d

 ##' @export

 fortify.phylo4d <- function(model, data, layout="rectangular", yscale="none",

-                            ladderize=TRUE, right=FALSE, mrsd=NULL, ...) {

-    res <- fortify.phylo4(model, data, layout, yscale="none", ladderize, right, mrsd, ...) # not apply yscale at this moment

+                            ladderize=TRUE, right=FALSE, branch.length="branch.length",

+                            mrsd=NULL, ...) {

+    phylo <- set_branch_length(model, branch.length)

+    res <- fortify(phylo, data, layout, branch.length=branch.length,

+                   ladderize, right, mrsd, ...)

     tdata <- model@data[match(res$node, rownames(model@data)), , drop=FALSE]

     df <- cbind(res, tdata)

     scaleY(as.phylo.phylo4(model), df, yscale, layout, ...)

@@ -685,9 +688,3 @@ fortify.phyloseq <- function(model, data, layout="rectangular",

 ##     ggplot(df) + geom_tree()

 ## }

-

-

-calculate_angle <- function(data) {

-    data$angle <- 360/(diff(range(data$y)) + 1) * data$y

-    return(data)

-}

@@ -65,7 +65,7 @@ setMethod("get.tree", signature(object="phylo"),

           function(object, ...) {

               return(object)

           })

-         

+        

 ##' @rdname get.tree-methods

 ##' @exportMethod get.tree

 setMethod("get.tree", signature(object="data.frame"),

@@ -85,4 +85,5 @@ setMethod("get.tree", signature(object="data.frame"),

               phylo$Nnode <- sum(!object[, "isTip"])

               class(phylo) <- "phylo"

               return(phylo)

-          })

\ No newline at end of file

+          })

+

@@ -15,19 +15,11 @@ setMethod("groupClade", signature(object="codeml"),

 ##' @rdname groupClade-methods

 ##' @exportMethod groupClade

-setMethod("groupClade", signature(object="gg"),

+setMethod("groupClade", signature(object="ggtree"),

           function(object, node, group_name) {

-              groupClade.ggplot(object, node, group_name)

+              groupClade.ggtree(object, node, group_name)

           })

-##' @rdname groupClade-methods

-##' @exportMethod groupClade

-setMethod("groupClade", signature(object="ggplot"),

-          function(object, node, group_name) {

-              groupClade.ggplot(object, node, group_name)

-          })

-

-

 ##' @rdname groupClade-methods

 ##' @exportMethod groupClade

 setMethod("groupClade", signature(object="jplace"),

@@ -38,7 +30,7 @@ setMethod("groupClade", signature(object="jplace"),

 ##' group selected clade

 ##'

-##' 

+##'

 ##' @rdname groupClade-methods

 ##' @exportMethod groupClade

 setMethod("groupClade", signature(object="nhx"),

@@ -73,7 +65,7 @@ groupClade.phylo <- function(object, node, group_name) {

             clade$tip.label

         })

     }

-    

+

     groupOTU.phylo(object, tips, group_name)

 }

@@ -88,7 +80,7 @@ groupClade_ <- function(object, node, group_name) {

 }

-groupClade.ggplot <- function(object, nodes, group_name) {

+groupClade.ggtree <- function(object, nodes, group_name) {

     df <- object$data

     df[, group_name] <- 0

     for (node in nodes) {

@@ -32,18 +32,12 @@ setMethod("groupOTU", signature(object="codeml_mlc"),

           }

           )

-##' @rdname groupOTU-methods

-##' @exportMethod groupOTU

-setMethod("groupOTU", signature(object="gg"),

-          function(object, focus, group_name, ...) {

-              groupOTU.ggplot(object, focus, group_name, ...)

-          })

 ##' @rdname groupOTU-methods

 ##' @exportMethod groupOTU

-setMethod("groupOTU", signature(object="ggplot"),

+setMethod("groupOTU", signature(object="ggtree"),

           function(object, focus, group_name="group", ...) {

-              groupOTU.ggplot(object, focus, group_name, ...)

+              groupOTU.ggtree(object, focus, group_name, ...)

           })

@@ -185,7 +179,7 @@ groupOTU_ <- function(object, focus, group_name, ...) {

 }

-groupOTU.ggplot <- function(object, focus, group_name, ...) {

+groupOTU.ggtree <- function(object, focus, group_name, ...) {

     df <- object$data

     df[, group_name] <- 0

     object$data <- groupOTU.df(df, focus, group_name, ...)

@@ -38,7 +38,7 @@ gzoom.phylo <- function(phy, focus, subtree=FALSE, widths=c(.3, .7)) {

     invisible(list(p1=p1, p2=p2))

 }

-gzoom.ggplot <- function(tree_view, focus, widths=c(.3, .7), xmax_adjust=0) {

+gzoom.ggtree <- function(tree_view, focus, widths=c(.3, .7), xmax_adjust=0) {

     node <- MRCA(tree_view, focus)

     cpos <- get_clade_position(tree_view, node)

     p2 <- with(cpos, tree_view+

@@ -51,9 +51,9 @@ gzoom.ggplot <- function(tree_view, focus, widths=c(.3, .7), xmax_adjust=0) {

 ##' @rdname gzoom-methods

 ##' @exportMethod gzoom

 ##' @param xmax_adjust adjust xmax (xlim[2])

-setMethod("gzoom", signature(object="gg"),

+setMethod("gzoom", signature(object="ggtree"),

           function(object, focus, widths=c(.3, .7), xmax_adjust=0) {

-              gzoom.ggplot(object, focus, widths, xmax_adjust)

+              gzoom.ggtree(object, focus, widths, xmax_adjust)

           })

similarity index 95%

rename from R/tree-utilities.R

rename to R/tidytree.R

@@ -1,6 +1,6 @@

 ##' convert tip or node label(s) to internal node number

 ##'

-##' 

+##'

 ##' @title nodeid

 ##' @param x tree object or graphic object return by ggtree

 ##' @param label tip or node label(s)

@@ -8,9 +8,9 @@

 ##' @export

 ##' @author Guangchuang Yu

 nodeid <- function(x, label) {

-    if (is(x, "gg")) 

+    if (is(x, "gg"))

         return(nodeid.gg(x, label))

-    

+

     nodeid.tree(x, label)

 }

@@ -40,14 +40,14 @@ reroot_node_mapping <- function(tree, tree2) {

         ip <- getParent(tree, k)

         if (node_map[ip, "visited"])

             next

-        

+

         cc <- getChild(tree, ip)

         node2 <- node_map[cc,2]

         if (anyNA(node2)) {

             node <- c(node, k)

             next

         }

-        

+

         to <- unique(sapply(node2, getParent, tr=tree2))

         to <- to[! to %in% node_map[,2]]

         node_map[ip, 2] <- to

@@ -75,7 +75,7 @@ layout.unrooted <- function(tree) {

     df[root, "start"] <- 0

     df[root, "end"]   <- 2

     df[root, "angle"] <- 0

-    

+

     nb.sp <- sapply(1:N, function(i) length(get.offspring.tip(tree, i)))

     nodes <- getNodes_by_postorder(tree)

@@ -83,22 +83,22 @@ layout.unrooted <- function(tree) {

     for(curNode in nodes) {

         curNtip <- nb.sp[curNode]

         children <- getChild(tree, curNode)

-        

+

         start <- df[curNode, "start"]

         end <- df[curNode, "end"]

-        

+

         if (length(children) == 0) {

             ## is a tip

             next

         }

-        

+

         for (i in seq_along(children)) {

             child <- children[i]

             ntip.child <- nb.sp[child]

-            

+

             alpha <- (end - start) * ntip.child/curNtip

             beta <- start + alpha / 2

-            

+

             length.child <- df[child, "length"]

             df[child, "x"] <- df[curNode, "x"] + cospi(beta) * length.child

             df[child, "y"] <- df[curNode, "y"] + sinpi(beta) * length.child

@@ -107,9 +107,9 @@ layout.unrooted <- function(tree) {

             df[child, "end"] <- start + alpha

             start <- start + alpha

         }

-        

+

     }

-    

+

     return(df)

 }

@@ -118,7 +118,7 @@ getParent.df <- function(df, node) {

     res <- df$parent[i]

     if (res == node) {

         ## root node

-        return(0) 

+        return(0)

     }

     return(res)

 }

@@ -165,10 +165,10 @@ get.offspring.df <- function(df, node) {

     return(sp)

 }

-    

+

 ##' extract offspring tips

 ##'

-##' 

+##'

 ##' @title get.offspring.tip

 ##' @param tr tree

 ##' @param node node

@@ -188,17 +188,17 @@ get.offspring.tip <- function(tr, node) {

 ##' calculate total number of nodes

 ##'

-##' 

+##'

 ##' @title getNodeNum

 ##' @param tr phylo object

-##' @return number 

+##' @return number

 ##' @author Guangchuang Yu

 ##' @export

 getNodeNum <- function(tr) {

     Ntip <- length(tr[["tip.label"]])

     Nnode <- tr[["Nnode"]]

     ## total nodes

-    N <- Ntip + Nnode 

+    N <- Ntip + Nnode

     return(N)

 }

@@ -233,7 +233,7 @@ getSibling <- function(tr, node) {

     if (node == root) {

         return(NA)

     }

-    

+

     parent <- getParent(tr, node)

     child <- getChild(tr, parent)

     sib <- child[child != node]

@@ -272,9 +272,9 @@ getNodeName <- function(tr) {

     return(nodeName)

 }

-##' get the root number 

+##' get the root number

+##'

 ##'

-##' 

 ##' @title getRoot

 ##' @param tr phylo object

 ##' @return root number

@@ -286,7 +286,7 @@ getRoot <- function(tr) {

     ## 2nd col is child,

     if (!is.null(attr(tr, "order")) && attr(tr, "order") == "postorder")

         return(edge[nrow(edge), 1])

-    

+

     parent <- unique(edge[,1])

     child <- unique(edge[,2])

     ## the node that has no parent should be the root

@@ -306,7 +306,7 @@ get.trunk <- function(tr) {

 ##' path from start node to end node

 ##'

-##' 

+##'

 ##' @title get.path

 ##' @param phylo phylo object

 ##' @param from start node

@@ -333,7 +333,7 @@ get.path_length <- function(phylo, from, to, weight=NULL) {

     if (is.null(weight)) {

         return(length(path)-1)

     }

-    

+

     df <- fortify(phylo)

     if ( ! (weight %in% colnames(df))) {

         stop("weight should be one of numerical attributes of the tree...")

@@ -350,7 +350,7 @@ get.path_length <- function(phylo, from, to, weight=NULL) {

         ee <- get_edge_index(df, path[i], path[i+1])

         res <- res + df[ee, weight]

     }

-    

+

     return(res)

 }

@@ -362,7 +362,7 @@ getNodes_by_postorder <- function(tree) {

 getXcoord2 <- function(x, root, parent, child, len, start=0, rev=FALSE) {

     x[root] <- start

     x[-root] <- NA  ## only root is set to start, by default 0

-        

+

     currentNode <- root

     direction <- 1

     if (rev == TRUE) {

@@ -395,12 +395,12 @@ getXcoord_no_length <- function(tr) {

     cl <- split(child, parent)

     child_list <- list()

     child_list[as.numeric(names(cl))] <- cl

-    

+

     while(anyNA(x)) {

         idx <- match(currentNode, child)

         pNode <- parent[idx]

         ## child number table

-        p1 <- table(parent[parent %in% pNode]) 

+        p1 <- table(parent[parent %in% pNode])

         p2 <- table(pNode)

         np <- names(p2)

         i <- p1[np] == p2

@@ -419,8 +419,8 @@ getXcoord_no_length <- function(tr) {

         currentNode <- unique(c(currentNode, newNode))

     }

-    x <- x - min(x) 

-    return(x)    

+    x <- x - min(x)

+    return(x)

 }

@@ -439,7 +439,7 @@ getXcoord <- function(tr) {

 }

 getXYcoord_slanted <- function(tr) {

-    

+

     edge <- tr$edge

     parent <- edge[,1]

     child <- edge[,2]

@@ -470,7 +470,7 @@ getYcoord <- function(tr, step=1) {

     cl <- split(child, parent)

     child_list <- list()

     child_list[as.numeric(names(cl))] <- cl

-    

+

     y <- numeric(N)

     tip.idx <- child[child <= Ntip]

     y[tip.idx] <- 1:Ntip * step

@@ -486,19 +486,19 @@ getYcoord <- function(tr, step=1) {

         ## idx <- sapply(pNode, function(i) all(child[parent == i] %in% currentNode))

         idx <- sapply(pNode, function(i) all(child_list[[i]] %in% currentNode))

         newNode <- pNode[idx]

-        

+

         y[newNode] <- sapply(newNode, function(i) {

             mean(y[child_list[[i]]], na.rm=TRUE)

-            ##child[parent == i] %>% y[.] %>% mean(na.rm=TRUE)           

+            ##child[parent == i] %>% y[.] %>% mean(na.rm=TRUE)

         })

-        

+

         currentNode <- c(currentNode[!currentNode %in% unlist(child_list[newNode])], newNode)

         ## currentNode <- c(currentNode[!currentNode %in% child[parent %in% newNode]], newNode)

         ## parent %in% newNode %>% child[.] %>%

         ##     `%in%`(currentNode, .) %>% `!` %>%

         ##         currentNode[.] %>% c(., newNode)

     }

-    

+

     return(y)

 }

@@ -507,7 +507,7 @@ getYcoord_scale <- function(tr, df, yscale) {

     N <- getNodeNum(tr)

     y <- numeric(N)

-    

+

     root <- getRoot(tr)

     y[root] <- 0

     y[-root] <- NA

@@ -539,7 +539,7 @@ getYcoord_scale <- function(tr, df, yscale) {

 getYcoord_scale2 <- function(tr, df, yscale) {

     root <- getRoot(tr)

-    

+

     pathLength <- sapply(1:length(tr$tip.label), function(i) {

         get.path_length(tr, i, root, yscale)

     })

@@ -566,7 +566,7 @@ getYcoord_scale2 <- function(tr, df, yscale) {

         } else {

             ordered_tip <- c(ordered_tip[1:ii],sib)

         }

-        

+

         ii <- ii + nn + 1

     }

@@ -576,20 +576,20 @@ getYcoord_scale2 <- function(tr, df, yscale) {

     N <- getNodeNum(tr)

     y <- numeric(N)

-    

+

     y[root] <- 0

     y[-root] <- NA

     ## yy <- df[, yscale]

     ## yy[is.na(yy)] <- 0

-    

+

     for (i in 2:length(long_branch)) {

         y[long_branch[i]] <- y[long_branch[i-1]] + df[long_branch[i], yscale]

     }

-    

+

     parent <- df[, "parent"]

     child <- df[, "node"]

-    

+

     currentNodes <- root

     while(anyNA(y)) {

         newNodes <- c()

@@ -618,14 +618,14 @@ getYcoord_scale2 <- function(tr, df, yscale) {

 getYcoord_scale_numeric <- function(tr, df, yscale, ...) {

     df <- .assign_parent_status(tr, df, yscale)

     df <- .assign_child_status(tr, df, yscale)

-    

+

     y <- df[, yscale]

     if (anyNA(y)) {

         warning("NA found in y scale mapping, all were setting to 0")

         y[is.na(y)] <- 0

     }

-    

+

     return(y)

 }

@@ -655,7 +655,7 @@ getYcoord_scale_numeric <- function(tr, df, yscale, ...) {

     if (!is.null(yscale_mapping)) {

         yy <- yscale_mapping[yy]

     }

-    

+

     na.idx <- which(is.na(yy))

     if (length(na.idx) > 0) {

         tree <- get.tree(tr)

@@ -680,7 +680,7 @@ getYcoord_scale_numeric <- function(tr, df, yscale, ...) {

 getYcoord_scale_category <- function(tr, df, yscale, yscale_mapping=NULL, ...) {

     if (is.null(yscale_mapping)) {

         stop("yscale is category variable, user should provide yscale_mapping,

-             which is a named vector, to convert yscale to numberical values...") 

+             which is a named vector, to convert yscale to numberical values...")

     }

     if (! is(yscale_mapping, "numeric") ||

         is.null(names(yscale_mapping))) {

@@ -694,11 +694,11 @@ getYcoord_scale_category <- function(tr, df, yscale, yscale_mapping=NULL, ...) {

             df[ii, yscale] <- df[ii, "node"]

         }

     }

-    

+

     ## assign to parent status is more prefer...

     df <- .assign_parent_status(tr, df, yscale)

     df <- .assign_child_status(tr, df, yscale, yscale_mapping)

-    

+

     y <- df[, yscale]

     if (anyNA(y)) {

@@ -733,37 +733,43 @@ calculate_branch_mid <- function(res) {

 set_branch_length <- function(tree_object, branch.length) {

-    phylo <- get.tree(tree_object)

-    

+    if (is(tree_object, "phylo4d")) {

+        phylo <- as.phylo.phylo4(tree_object)

+        d <- tree_object@data

+        tree_anno <- data.frame(node=rownames(d), d)

+    } else {

+        phylo <- get.tree(tree_object)

+    }

+

     if (branch.length %in%  c("branch.length", "none")) {

         return(phylo)

     }

-

     ## if (is(tree_object, "codeml")) {

     ##     tree_anno <- tree_object@mlc@dNdS

     ## } else

-    

+

     if (is(tree_object, "codeml_mlc")) {

         tree_anno <- tree_object@dNdS

     } else if (is(tree_object, "beast")) {

         tree_anno <- tree_object@stats

     }

+

     if (has.extraInfo(tree_object)) {

         tree_anno <- merge(tree_anno, tree_object@extraInfo, by.x="node", by.y="node")

     }

     cn <- colnames(tree_anno)

     cn <- cn[!cn %in% c('node', 'parent')]

-    

+

     length <- match.arg(branch.length, cn)

     if (all(is.na(as.numeric(tree_anno[, length])))) {

         stop("branch.length should be numerical attributes...")

     }

-    

+

     edge <- as.data.frame(phylo$edge)

     colnames(edge) <- c("parent", "node")

-    

+

     dd <- merge(edge, tree_anno,

                 by.x  = "node",

                 by.y  = "node",

@@ -772,7 +778,7 @@ set_branch_length <- function(tree_object, branch.length) {

     len <- unlist(dd[, length])

     len <- as.numeric(len)

     len[is.na(len)] <- 0

-    

+

     phylo$edge.length <- len

     return(phylo)

@@ -785,7 +791,7 @@ re_assign_ycoord_df <- function(df, currentNode) {

         idx <- sapply(pNode, function(i) with(df, all(node[parent == i & parent != node] %in% currentNode)))

         newNode <- pNode[idx]

         ## newNode <- newNode[is.na(df[match(newNode, df$node), "y"])]

-        

+

         df[match(newNode, df$node), "y"] <- sapply(newNode, function(i) {

             with(df, mean(y[parent == i], na.rm = TRUE))

         })

@@ -795,3 +801,20 @@ re_assign_ycoord_df <- function(df, currentNode) {

     return(df)

 }

+

+##' test whether input object is produced by ggtree function

+##'

+##'

+##' @title is.ggtree

+##' @param x object

+##' @return TRUE or FALSE

+##' @export

+##' @author guangchuang yu

+is.ggtree <- function(x) inherits(x, 'ggtree')

+

+

+

+calculate_angle <- function(data) {

+    data$angle <- 360/(diff(range(data$y)) + 1) * data$y

+    return(data)

+}

similarity index 74%

rename from R/NHX.R

rename to R/treeio.R

@@ -1,3 +1,13 @@

+filename <- function(file) {

+    ## textConnection(text_string) will work just like a file

+    ## in this case, just set the filename as ""

+    file_name <- ""

+    if (is.character(file)) {

+        file_name <- file

+    }

+    return(file_name)

+}

+

 ##' read nhx tree file

 ##'

 ##'

@@ -103,3 +113,47 @@ setMethod("get.fields", signature(object="nhx"),

               get.fields.tree(object)

           }

           )

+

+

+Ntip <- function(tree) {

+    phylo <- get.tree(tree)

+    length(phylo$tip.label)

+}

+

+Nnode <- function(tree, internal.only=TRUE) {

+    phylo <- get.tree(tree)

+    if (internal.only)

+        return(phylo$Nnode)

+

+    Ntip(phylo) + phylo$Nnode

+}

+

+

+has.extraInfo <- function(object) {

+    if (!is.tree(object)) {

+        return(FALSE)

+    }

+

+    if (! .hasSlot(object, "extraInfo")) {

+        return(FALSE)

+    }

+

+    extraInfo <- object@extraInfo

+

+    if (nrow(extraInfo) > 0) {

+        return(TRUE)

+    }

+

+    return(FALSE)

+}

+

+##' @importFrom methods .hasSlot is missingArg new slot slot<-

+has.slot <- function(object, slotName) {

+    if (!isS4(object)) {

+        return(FALSE)

+    }

+    .hasSlot(object, slotName)

+    ## slot <- tryCatch(slot(object, slotName), error=function(e) NULL)

+    ## ! is.null(slot)

+}

+

@@ -1,26 +1,6 @@

-Ntip <- function(tree) {

-    phylo <- get.tree(tree)

-    length(phylo$tip.label)

-}

-

-Nnode <- function(tree, internal.only=TRUE) {

-    phylo <- get.tree(tree)

-    if (internal.only)

-        return(phylo$Nnode)

-    Ntip(phylo) + phylo$Nnode

-}

-filename <- function(file) {

-    ## textConnection(text_string) will work just like a file

-    ## in this case, just set the filename as ""

-    file_name <- ""

-    if (is.character(file)) {

-        file_name <- file

-    }

-    return(file_name)

-}

 ##' @importFrom ggplot2 last_plot

@@ -32,15 +12,6 @@ get_tree_view <- function(tree_view) {

 }

-##' @importFrom methods .hasSlot is missingArg new slot slot<-

-has.slot <- function(object, slotName) {

-    if (!isS4(object)) {

-        return(FALSE)

-    }

-    .hasSlot(object, slotName)

-    ## slot <- tryCatch(slot(object, slotName), error=function(e) NULL)

-    ## ! is.null(slot)

-}

 has.field <- function(tree_object, field) {

     if ( ! field %in% get.fields(tree_object) ) {

@@ -80,24 +51,6 @@ has.field <- function(tree_object, field) {

     return(res)

 }

-has.extraInfo <- function(object) {

-    if (!is.tree(object)) {

-        return(FALSE)

-    }

-

-    if (! .hasSlot(object, "extraInfo")) {

-        return(FALSE)

-    }

-

-    extraInfo <- object@extraInfo

-

-    if (nrow(extraInfo) > 0) {

-        return(TRUE)

-    }

-

-    return(FALSE)

-}

-

 append_extraInfo <- function(df, object) {

     if (has.extraInfo(object)) {

         info <- object@extraInfo

@@ -4,7 +4,7 @@ ggtree: an R package for visualization and annotation of phylogenetic trees with

 [![releaseVersion](https://img.shields.io/badge/release%20version-1.6.4-green.svg?style=flat)](https://bioconductor.org/packages/ggtree) [![develVersion](https://img.shields.io/badge/devel%20version-1.7.4-green.svg?style=flat)](https://github.com/GuangchuangYu/ggtree) [![Bioc](http://www.bioconductor.org/shields/years-in-bioc/ggtree.svg)](https://www.bioconductor.org/packages/devel/bioc/html/ggtree.html#since) [![total](https://img.shields.io/badge/downloads-16878/total-blue.svg?style=flat)](https://bioconductor.org/packages/stats/bioc/ggtree) [![month](https://img.shields.io/badge/downloads-1621/month-blue.svg?style=flat)](https://bioconductor.org/packages/stats/bioc/ggtree) <img src="logo.png" align="right" />

-[![Project Status: Active - The project has reached a stable, usable state and is being actively developed.](http://www.repostatus.org/badges/latest/active.svg)](http://www.repostatus.org/#active) [![codecov](https://codecov.io/gh/GuangchuangYu/ggtree/branch/master/graph/badge.svg)](https://codecov.io/gh/GuangchuangYu/ggtree) [![Last-changedate](https://img.shields.io/badge/last%20change-2016--12--01-green.svg)](https://github.com/GuangchuangYu/ggtree/commits/master) [![GitHub forks](https://img.shields.io/github/forks/GuangchuangYu/ggtree.svg)](https://github.com/GuangchuangYu/ggtree/network) [![GitHub stars](https://img.shields.io/github/stars/GuangchuangYu/ggtree.svg)](https://github.com/GuangchuangYu/ggtree/stargazers) [![Awesome](https://cdn.rawgit.com/sindresorhus/awesome/d7305f38d29fed78fa85652e3a63e154dd8e8829/media/badge.svg)](https://awesome-r.com/#awesome-r-graphic-displays)

+[![Project Status: Active - The project has reached a stable, usable state and is being actively developed.](http://www.repostatus.org/badges/latest/active.svg)](http://www.repostatus.org/#active) [![codecov](https://codecov.io/gh/GuangchuangYu/ggtree/branch/master/graph/badge.svg)](https://codecov.io/gh/GuangchuangYu/ggtree) [![Last-changedate](https://img.shields.io/badge/last%20change-2016--12--06-green.svg)](https://github.com/GuangchuangYu/ggtree/commits/master) [![GitHub forks](https://img.shields.io/github/forks/GuangchuangYu/ggtree.svg)](https://github.com/GuangchuangYu/ggtree/network) [![GitHub stars](https://img.shields.io/github/stars/GuangchuangYu/ggtree.svg)](https://github.com/GuangchuangYu/ggtree/stargazers) [![Awesome](https://cdn.rawgit.com/sindresorhus/awesome/d7305f38d29fed78fa85652e3a63e154dd8e8829/media/badge.svg)](https://awesome-r.com/#awesome-r-graphic-displays)

 [![platform](http://www.bioconductor.org/shields/availability/devel/ggtree.svg)](https://www.bioconductor.org/packages/devel/bioc/html/ggtree.html#archives) [![Build Status](http://www.bioconductor.org/shields/build/devel/bioc/ggtree.svg)](https://bioconductor.org/checkResults/devel/bioc-LATEST/ggtree/) [![Linux/Mac Travis Build Status](https://img.shields.io/travis/GuangchuangYu/ggtree/master.svg?label=Mac%20OSX%20%26%20Linux)](https://travis-ci.org/GuangchuangYu/ggtree) [![AppVeyor Build Status](https://img.shields.io/appveyor/ci/Guangchuangyu/ggtree/master.svg?label=Windows)](https://ci.appveyor.com/project/GuangchuangYu/ggtree) [![install with bioconda](https://img.shields.io/badge/install%20with-bioconda-green.svg?style=flat)](http://bioconda.github.io/recipes/bioconductor-ggtree/README.html)

@@ -18,7 +18,7 @@ Please cite the following article when using `ggtree`:

 **G Yu**, DK Smith, H Zhu, Y Guan, TTY Lam<sup>\*</sup>. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. ***Methods in Ecology and Evolution***. *accepted*

-[![doi](https://img.shields.io/badge/doi-10.1111/2041--210X.12628-green.svg?style=flat)](http://dx.doi.org/10.1111/2041-210X.12628) [![citation](https://img.shields.io/badge/cited%20by-1-green.svg?style=flat)](https://scholar.google.com.hk/scholar?oi=bibs&hl=en&cites=7268358477862164627) [![Altmetric](https://img.shields.io/badge/Altmetric-274-green.svg?style=flat)](https://www.altmetric.com/details/10533079)

+[![doi](https://img.shields.io/badge/doi-10.1111/2041--210X.12628-green.svg?style=flat)](http://dx.doi.org/10.1111/2041-210X.12628) [![citation](https://img.shields.io/badge/cited%20by-1-green.svg?style=flat)](https://scholar.google.com.hk/scholar?oi=bibs&hl=en&cites=7268358477862164627) [![Altmetric](https://img.shields.io/badge/Altmetric-284-green.svg?style=flat)](https://www.altmetric.com/details/10533079)

 ------------------------------------------------------------------------

@@ -1,5 +1,5 @@

 % Generated by roxygen2: do not edit by hand

-% Please edit documentation in R/AllGenerics.R, R/NHX.R, R/RAxML.R, R/ape.R, R/beast.R, R/codeml.R, R/codeml_mlc.R, R/hyphy.R, R/jplace.R, R/paml_rst.R, R/phangorn.R, R/r8s.R

+% Please edit documentation in R/AllGenerics.R, R/RAxML.R, R/ape.R, R/beast.R, R/codeml.R, R/codeml_mlc.R, R/hyphy.R, R/jplace.R, R/paml_rst.R, R/phangorn.R, R/r8s.R, R/treeio.R

 \docType{methods}

 \name{get.fields}

 \alias{get.fields}

@@ -18,8 +18,6 @@

 \usage{

 get.fields(object, ...)

-\S4method{get.fields}{nhx}(object, ...)

-

 \S4method{get.fields}{raxml}(object, ...)

 \S4method{get.fields}{apeBootstrap}(object, ...)

@@ -39,6 +37,8 @@ get.fields(object, ...)

 \S4method{get.fields}{phangorn}(object, ...)

 \S4method{get.fields}{r8s}(object, ...)

+

+\S4method{get.fields}{nhx}(object, ...)

 }

 \arguments{

 \item{object}{one of \code{jplace}, \code{beast}, \code{hyphy}, \code{codeml}, \code{codeml_mlc}, \code{paml_rst} object}

@@ -1,5 +1,5 @@

 % Generated by roxygen2: do not edit by hand

-% Please edit documentation in R/tree-utilities.R

+% Please edit documentation in R/tidytree.R

 \name{get.offspring.tip}

 \alias{get.offspring.tip}

 \title{get.offspring.tip}

@@ -1,5 +1,5 @@

 % Generated by roxygen2: do not edit by hand

-% Please edit documentation in R/tree-utilities.R

+% Please edit documentation in R/tidytree.R

 \name{get.path}

 \alias{get.path}

 \title{get.path}

@@ -1,5 +1,5 @@

 % Generated by roxygen2: do not edit by hand

-% Please edit documentation in R/tree-utilities.R

+% Please edit documentation in R/tidytree.R

 \name{getNodeNum}

 \alias{getNodeNum}

 \title{getNodeNum}

@@ -1,5 +1,5 @@

 % Generated by roxygen2: do not edit by hand

-% Please edit documentation in R/tree-utilities.R

+% Please edit documentation in R/tidytree.R

 \name{getRoot}

 \alias{getRoot}

 \title{getRoot}

@@ -7,8 +7,7 @@

 \alias{groupClade,beast-method}

 \alias{groupClade,codeml-method}

 \alias{groupClade,codeml_mlc-method}

-\alias{groupClade,gg-method}

-\alias{groupClade,ggplot-method}

+\alias{groupClade,ggtree-method}

 \alias{groupClade,hyphy-method}

 \alias{groupClade,jplace-method}

 \alias{groupClade,nhx-method}

@@ -34,9 +33,7 @@ groupClade(object, node, group_name = "group", ...)