An implementation of Kmult (Blomberg et al. 2003; Adams 2014) in which Kmult is compared to the null distribution generated by simulation on a star phylogeny (rather than phylogenetic permutation). Additionally, the expectation of the ratio of raw to phylogenetic mean squared error (the denominator of the formula for K) is estimated via phylogenetic simulation under Brownian motion.

1 |

`model` |
An object of class |

`nsim` |
The number of iterations for phylogenetic simulation. |

`plot` |
Whether or not to plot the null and alternative distribution of Kmult. |

An object of class `compare.model`

.

Eric W. Goolsby

Golsby E.W. 2015. Likelihood-based parameter estimation for high-dimensional phylogenetic comparative models: overcoming the limitations of 'distance-based' methods. In review.

Adams D.C. 2014. A generalized K statistic for estimating phylogenetic signal from shape and other high-dimensional multivariate data. Syst. Biol. 63:685-697.

`evo.model`

`compare.models`

`physignal`

1 2 3 4 5 | ```
rand.data <- sim.traits()
null.model <- evo.model(tree = rand.data$tree,
Y = rand.data$trait_data,method = "Pairwise ML")
K.mult(model = null.model,nsim = 100)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.