We used the R package diversitree69to fit the models of trait evolution with maximum likelihood and reconstruct marginal ancestral state probabilities under the ML model. The two-epoch model was fit by sweeping across possible epoch boundary times at 1 my intervals between 1 and 120 mya, fitti...