DR function

Calculate the tip diversification rates (DR) for a phylogenetic tree