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/chronosCI documentation built on Jan. 30, 2023, 5:34 a.m.