Files: OUSimTraj.R: code to simulate trajectories of 1000 replicates under the OU process (Fig. 1) *.xml: XML inputs for BEAST2 analyses. Abbreviations: div: diversified extant sampling rnd: random extant sampling cst: constant-rate FBD iid10: skyline FBD with independent and identically distributed rates among 10 time intervals ou10: skyline FBD with Ornstein-Uhlenbeck smoothing hyperprior rates among 10 time intervals ou20: skyline FBD with Ornstein-Uhlenbeck smoothing hyperprior rates among 20 time intervals drs: d-r-s parameterization for the birth, death and sampling rates, instead of direct parameterization *.pdf: Maximum clade credibility trees Abbreviations: same as above