sis_test | R Documentation |
Compute multiple tests based on sister group comparisons
sis_test(pairs, drop_matches = TRUE, warn = TRUE)
pairs |
Data.frame with one row per sister group comparison, with one column for number of taxa in state 0, and one column for the number of taxa in state 1. |
drop_matches |
Drop sister group comparisons with equal numbers of taxa |
warn |
Some tests will fail with warnings (too few sister groups or other reasons). Setting this to FALSE will suppress those |
A vector with the results of many tests, as well as summary data for the comparisons
data(geospiza, package="geiger")
cleaned <- sis_clean(geospiza$phy, geospiza$dat)
phy <- cleaned$phy
traits <- cleaned$traits
trait <- sis_discretize(traits[,1])
sisters <- sis_get_sisters(phy)
sisters_comparison <- sis_format_comparison(sisters, trait, phy)
pairs <- sis_format_simpified(sisters_comparison)
sis_test(pairs)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.