countSimmap {phytools}R Documentation

Counts the number of character changes on a SIMMAP style tree or set of trees

Description

This function takes a tree or a set of trees with a mapped discrete character (SIMMAP style, e.g., see make.simmap or read.simmap), and computes the total number of character changes as well as the number of character changes between all states.

Usage

countSimmap(tree, states=NULL, message=TRUE)

Arguments

tree

a single tree or a set of trees with a mapped discrete character (e.g, see make.simmap or read.simmap.

states

optional argument with the states for the mapped character. If not provided, these will be computed from the tree. This is useful if averaging across many trees, some of which may lack certain states.

message

optional logical argument indicating whether or not to return an informative message about the function output.

Value

A list with up to three elements: N is an integer value giving the total number of character changes on the tree; Tr gives the number of of transitions between row and column states (or a matrix containing both N and the transitions between states, in rows, for an object of class "multiPhylo"); and (optionally) message contains an explanatory message about the function output.

Author(s)

Liam Revell liam.revell@umb.edu

References

Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.

Examples

tree<-pbtree(n=100,scale=1)
Q<-matrix(c(-2,1,1,1,-2,1,1,1,-2),3,3)
colnames(Q)<-rownames(Q)<-c("A","B","C")
mtree<-sim.history(tree,Q)
countSimmap(mtree,states=rownames(Q))

[Package phytools version 0.4-60 Index]