TreeDriftTest | R Documentation |
Performs a regression drift test along a phylogeny using DriftTest function.
TreeDriftTest(tree, mean.list, cov.matrix.list, sample.sizes = NULL)
tree |
phylogenetic tree |
mean.list |
list of tip node means. Names must match tip node labels. |
cov.matrix.list |
list of tip node covariance matrices. Names must match tip node labels. |
sample.sizes |
vector of tip nodes sample sizes |
A list of regression drift tests performed in nodes with over 4 descendant tips.
Diogo Melo
DriftTest PlotTreeDriftTest
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)
#Ancestral node plot:
test.list[[length(test.list)]]$plot
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.