Modeling the evolution of rates of continuous trait evolution

Rates of phenotypic evolution vary markedly across the tree of life, from the accelerated evolution apparent in adaptive radiations to the remarkable evolutionary stasis exhibited by so-called “living fossils”. Such rate variation has important consequences for large-scale evolutionary dynamics, generating vast disparities in phenotypic diversity across space, time, and taxa. Despite this, most methods for estimating trait evolution rates assume rates vary deterministically with respect to some variable of interest or change infrequently during a clade’s history. These assumptions may cause underfitting of trait evolution models and mislead hypothesis testing. Here, we develop a new trait evolution model that allows rates to vary gradually and stochastically across a clade. Further, we extend this model to accommodate generally decreasing or increasing rates over time, allowing for flexible modeling of “early/late bursts” of trait evolution. We implement a Bayesian method, termed “evolv..., For the empirical example analyzing cetacean body length data (\"05_cet_prep_data.R\" and \"06_cet_figures.R\"): the cetacean phylogeny is from a recent study (see Lloyd et al. 2021 and associated dryad repository). We used the tree \"Cetacea_Safe_Extant_MCC.tre\" located in the compressed file \"Cetacean_Metatree_Data.zip\" in \"Cetacean Metatree Data/Time TreeInference/Safe\". Most cetacean body length data comes from the \"whales\" dataset included in the R package geiger (after installing geiger, you can load this yourself using the code \"data('whales')\"). These data represent maximum adult female lengths in meters after a natural log transformation. All other cetacean body length data was found in various primary literature sources and hard-coded into an R script. When possible, we also based length data on maximum adult female lengths, though in some cases sex information was unavailable and/or all measured female specimens were juvenile. See Table S1 in the online appendix (\"appendix.pdf\") f..., The README file contains more detailed information on scripts, folders, and file names, including the meaning of various notations and their associated inputs/outputs.
All relevant scripts/data are included in the compressed file \"evorates.zip\". Scripts are in the main folder of the compressed file and come with a numeric prefix corresponding to the order in which they should be run. Data, intermediate files (i.e., the output of a previous script that will be used as input for a later script), and figures/tables are organized into various folders generally relating to the corresponding part of the analysis (e.g., files beginning with \"cet_\" correspond to the empirical analysis; see Code/Software for more detailed information on which files each folder contains).
For convenience, the most processed representations of the simulated data and associated model fits (\"inf.res.full\", \"tru.res.full\", and \"prior.sens.res.full\"), Hamiltonian Monte Carlo samples/accompanying info for models fitted...