examples/simulations/NELSI_RCodeOnly/get.df.R

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)
}
josephwb/chronos documentation built on Jan. 27, 2023, 2:35 a.m.