Plot results from TreeDriftTest
Plot which labels reject drift hypothesis.
PlotTreeDriftTest(test.list, tree, ...)
test.list |
Output from TreeDriftTest |
tree |
phylogenetic tree |
... |
adition arguments to plot |
Diogo Melo
DriftTest TreeDriftTest
library(ape) data(bird.orders) tree <- bird.orders mean.list <- llply(tree$tip.label, function(x) rnorm(5)) names(mean.list) <- tree$tip.label cov.matrix.list <- RandomMatrix(5, length(tree$tip.label)) names(cov.matrix.list) <- tree$tip.label sample.sizes <- runif(length(tree$tip.label), 15, 20) test.list <- TreeDriftTest(tree, mean.list, cov.matrix.list, sample.sizes) PlotTreeDriftTest(test.list, tree)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.