phylobase (version 0.6.8)

printphylo4: print a phylogeny

Description

Prints a phylo4 or phylo4d object in data.frame format with user-friendly column names

Usage

printphylo4(x, edgeOrder=c("pretty","real"), printall = TRUE)

Arguments

x

a phylo4 tree or phylo4d tree+data object

edgeOrder

in the data frame returned, the option 'pretty' returns the internal nodes followed by the tips, the option 'real' returns the nodes in the order they are stored in the edge matrix.

printall

default prints entire tree. printall=FALSE returns the first 6 rows

Value

• A data.frame with a row for each node (descendant), sorted as follows: root first, then other internal nodes, and finally tips. The returned data.frame has the following columns:
• label
  Label for the taxon at the node (usually species name).
• node
  Node number, i.e. the number identifying the node in x@edge.
• ancestor
  Node number of the node's ancestor.
• branch.length
  The branch length connecting the node to its ancestor (NAs if missing).
• node.type
  "root", "internal", or "tip". (internally generated)
• data
  phenotypic data associated with the nodes, with separate columns for each variable.

Details

This is a user-friendly version of the tree representation, useful for checking that objects were read in completely and translated correctly. The phylogenetic tree is represented as a list of numbered nodes, linked in a particular way through time (or rates of evolutionary change). The topology is given by the pattern of links from each node to its ancestor. Also given are the taxon names, node type (root/internal/tip) and phenotypic data (if any) associated with the node, and the branch length from the node to its ancestor. A list of nodes (descendants) and ancestors is minimally required for a phylo4 object.

Examples

tree.phylo <- ape::read.tree(text="((a,b),c);")
tree <- as(tree.phylo, "phylo4")
##plot(tree,show.node=TRUE) ## plotting broken with empty node labels: FIXME
tip.data <- data.frame(size=c(1,2,3), row.names=c("a", "b", "c"))
treedata <- phylo4d(tree, tip.data)
plot(treedata)
print(treedata)