This function takes a modified "phylo" object as input and returns a set of so-called phylogenetic covariance matrices as a list: one for each mapped state.
multiC(tree, internal=FALSE)a phylogeny with mapped discrete state in a modified object of class "phylo" (e.g., see read.simmap.
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.