
Last chance! 50% off unlimited learning
Sale ends in
branchClasses(tree, whichExtant = NULL, tol = 0.01)
#simulated example
set.seed(444)
taxa <- simFossilTaxa(p=0.1,q=0.1,nruns=1,mintaxa=20,maxtaxa=30,maxtime=1000,maxExtant=20)
tree <- taxa2phylo(taxa)
brlenRes <- branchClasses(tree)
#see frequency histograms of branch lengths
layout(1:4)
for(x in 1:length(brlenRes)){
hist(brlenRes[[x]],main="Branch Lengths",xlab=names(brlenRes)[x])
}
#see frequency histograms of branch depths
layout(1:4)
for(x in 1:length(brlenRes)){
hist(as.numeric(names(brlenRes[[x]])),main="Branch Depths",xlab=names(brlenRes)[x])
}
layout(1)
Run the code above in your browser using DataLab