R/reorder.phylo.R

   1 ## reorder.phylo.R (2010-04-02)
   2
   3 ##   Internal Reordering of Trees
   4
   5 ## Copyright 2006-2010 Emmanuel Paradis
   6
   7 ## This file is part of the R-package `ape'.
   8 ## See the file ../COPYING for licensing issues.
   9
  10 reorder.phylo <- function(x, order = "cladewise", ...)
  11 {
  12     order <- match.arg(order, c("cladewise", "pruningwise"))
  13     if (!is.null(attr(x, "order")))
  14         if (attr(x, "order") == order) return(x)
  15     nb.node <- x$Nnode
  16     if (nb.node == 1) return(x)
  17     nb.tip <- length(x$tip.label)
  18     nb.edge <- dim(x$edge)[1]
  19     neworder <- if (order == "cladewise")
  20       .C("neworder_cladewise", as.integer(nb.tip),
  21          as.integer(x$edge[, 1]), as.integer(x$edge[, 2]),
  22          as.integer(nb.edge), integer(nb.edge),
  23          PACKAGE = "ape")[[5]]
  24     else
  25       .C("neworder_pruningwise", as.integer(nb.tip),
  26          as.integer(nb.node), as.integer(x$edge[, 1]),
  27          as.integer(x$edge[, 2]), as.integer(nb.edge),
  28          integer(nb.edge), PACKAGE = "ape")[[6]]
  29     x$edge <- x$edge[neworder, ]
  30     if (!is.null(x$edge.length))
  31       x$edge.length <- x$edge.length[neworder]
  32     attr(x, "order") <- order
  33     x
  34 }