R/getDesc.R

In BAMMtools: Analysis and Visualization of Macroevolutionary Dynamics on Phylogenetic Trees

#############################################################
#
#	getDesc(....)
#
#	 returns a copy of the tree with a component 'desc_set', which
#	is a vector of all desc nodes in downpass sequences order

getDesc <- function(ephy, node)
{
	ephy$desc_set <- ephy$downseq[which(ephy$downseq == node):which(ephy$downseq == ephy$lastvisit[node])]
#	if (is.null(phy$desc_set)) {
#		phy$desc_set <- node;
#	}
#	
#	if (node > length(phy$tip.label)){
#		dset <- phy$edge[,2][phy$edge[,1] == node];
#		phy$desc_set <- c(phy$desc_set, dset[1]);
#		phy <- getDesc(phy, dset[1]);
#		phy$desc_set <- c(phy$desc_set, dset[2]);
#		phy <- getDesc(phy, dset[2]);
#		
#	}
# 
 	return(ephy);
}