PlotTreeDriftTest | R Documentation |
Plot which labels reject drift hypothesis.
PlotTreeDriftTest(test.list, tree, ...)
test.list |
Output from TreeDriftTest |
tree |
phylogenetic tree |
... |
adition arguments to plot |
No return value, called for plot side effects
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.