Following the result of Baele et al. , we used Bayes factors  calculated by marginal likelihoods derived from path sampling (PS)  and stepping-stone sampling (SS)  to compare the parametric fit of three clock models: exponential, lognormal and random. Since our sampling included all recognized species of Coptis and Xanthorhiza, a birth-death tree prior was used.