R/treebalance.R

In bcm-uga/apTreeshape: Analyses of Phylogenetic Treeshape

"treebalance" <-
function(nodes, names) {
	
#nodes est une matrice a deux colonnes
	if (!(is.matrix(nodes)) || (ncol(nodes) != 2)) {
		stop("invalid argument")
	}
#names est nu vecteur a n elements.
	if (missing(names)) {
		names=as.character(1:(nrow(nodes)+1))
	}
	
	if (length(names)!=nrow(nodes)+1){
		stop("Wrong number of elements for vector names")
	}
	
	res <- list(merge=nodes, names=names)
	class(res)<-'treebalance'
	res


	
}