This function conducts the 'D-test' of Huelsenbeck et al. (2003).
Dtest(t1, t2, nsim=100, ...)
set of stochastic map trees (i.e., object of class "multiSimmap"
for character 1. Note that t1
and t2
should be of the same length.
set of stochastic map trees (i.e., object of class "multiSimmap"
for character 2. Note that t1
and t2
should be of the same length.
number of simulations to use in the test.
arguments to be passed internally to make.simmap
. Note that (for now) these must be the same for both t1
and t2
(that is to say, we are not able to assume different trait evolution models for each tree).
An object of class "Dtest"
.
Note that this function has been included without much testing, and so the user should be wary.
Huelsenbeck, J. P., R. Neilsen, and J. P. Bollback (2003) Stochastic mapping of morphological characters. Systematic Biology, 52, 131-138.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.