R/get.df.R

Defines functions get.df

get.df <- 
function(tre){
	brlens <- as.numeric(scale(tre$edge.length, scale = F, center = T))
	df <- sum(brlens[which(tre$edge[,2] <= Ntip(tre))]) - sum(brlens)
	return(df)
}
sebastianduchene/NELSI documentation built on Aug. 18, 2022, 11:45 p.m.