This function takes a modified "phylo" object as input and returns a set of so-called phylogenetic covariance matrices (e.g., see vcv.phylo) as a list: one for each mapped state.
multiC(tree, internal=FALSE)an object of class "simmap" consisting of a phylogeny with a mapped discrete character.
logical value indicating whether or not internal nodes should be returned.
A list of matrices.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.