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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.