Plot results from TreeDriftTest
Plot which labels reject drift hypothesis.
PlotTreeDriftTest(test.list, tree, ...)
test.list
: Output from TreeDriftTesttree
: phylogenetic tree...
: adition arguments to plotNo return value, called for plot side effects
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)
DriftTest TreeDriftTest
Diogo Melo
Useful links