R/id.subtrees.R

In iteRates: Parametric rate comparison

id.subtrees <-
function(tree){
	if((length(tree$edge.length)-length(tree$tip.label)+1!=length(tree$tip.label)+1)==FALSE)stop("\nTree not fully bifurcating")
	tree$node.label<-1:length(subtrees2.6(tree))
	return(list(tree=tree,subtree=subtrees2.6(tree)))
	}