Simulate stochastic character map on a tree.
make.simmap(): This function simulates a stochastic character mapping on a tree.
Input
mtrees<-make.simmap(tree,x,model="SYM",nsim=1)
"tree" is a phylogenetic tree in "phylo" format;
"x" is a vector with the states at the tips of the tree for a discretely valued character;
"model" can be any model in the "ape" function ace();
"nsim" is the number of simulations.
Output
make.simmap() outputs an R "phylo" or "multiPhylo" object with:
"$maps" a list of named vectors containing the times spent in each state on each branch, in the order in which they occur;
"$mapped.edge" a matrix containing the total time spent in each state along each edge of the tree.
Citation
Function written by Liam J. Revell 2011.
home
Last modified July 13, 2011.