R/cladesize.R

"cladesize" <-
function(tree) {
	
	if (class(tree)!='treeshape') {
		stop("invalid arguments")
	}
	
	node=sample(1:nrow(tree$merge),1)
	clade=smaller.clade.spectrum(tree)
	child.number<-clade[node,1]	
	child.number
}

Try the apTreeshape package in your browser

Any scripts or data that you put into this service are public.

apTreeshape documentation built on Jan. 8, 2021, 2:07 a.m.