A variable-rate quantitative trait evolution model using penalized-likelihood

Liam J. Revell

doi:10.7717/peerj.11997

journal article Open Access Aug 17, 2021

A variable-rate quantitative trait evolution model using penalized-likelihood

Liam J. Revell

PeerJ Vol. 9 pp. e11997 · PeerJ

View at Publisher Save 10.7717/peerj.11997

Abstract

In recent years it has become increasingly popular to use phylogenetic comparative methods to investigate heterogeneity in the rate or process of quantitative trait evolution across the branches or clades of a phylogenetic tree. Here, I present a new method for modeling variability in the rate of evolution of a continuously-valued character trait on a reconstructed phylogeny. The underlying model of evolution is stochastic diffusion (Brownian motion), but in which the instantaneous diffusion rate (σ
2
)
also
evolves by Brownian motion on a logarithmic scale. Unfortunately, it’s not possible to simultaneously estimate the rates of evolution along each edge of the tree
and
the rate of evolution of σ
2
itself using Maximum Likelihood. As such, I propose a penalized-likelihood method in which the penalty term is equal to the log-transformed probability density of the rates under a Brownian model, multiplied by a ‘smoothing’ coefficient, λ, selected by the user. λ determines the magnitude of penalty that’s applied to rate variation between edges. Lower values of λ penalize rate variation relatively little; whereas larger λ values result in minimal rate variation among edges of the tree in the fitted model, eventually converging on a single value of σ
2
for all of the branches of the tree. In addition to presenting this model here, I have also implemented it as part of my
phytools
R package in the function
multirateBM
. Using different values of the penalty coefficient, λ, I fit the model to simulated data with: Brownian rate variation among edges (the model assumption); uncorrelated rate variation; rate changes that occur in discrete places on the tree; and no rate variation at all among the branches of the phylogeny. I then compare the estimated values of σ
2
to their known true values. In addition, I use the method to analyze a simple empirical dataset of body mass evolution in mammals. Finally, I discuss the relationship between the method of this article and other models from the phylogenetic comparative methods and finance literature, as well as some applications and limitations of the approach.

Topics

No keywords indexed for this article. Browse by subject →

References

63

[1]

Azzalini "The R package mnormt: the multivariate normal and t distributions. R package verison 2.0.2" (2020)

[2]

Beaulieu "Modeling stabilizing selection: expanding the Ornstein-Uhlenbeck model of adaptive evolution" Evolution (2012) 10.1111/j.1558-5646.2012.01619.x

[3]

Becker "Maps: draw geographical maps. R package version 3.3.0" (2018)

[4]

Beckers "The constant elasticity of variance model and its implications for option pricing" The Journal of Finance (1980) 10.1111/j.1540-6261.1980.tb03490.x

[5]

Blomberg "Beyond Brownian motion and the Ornstein-Uhlenbeck process: stochastic diffusion models for the evolution of quantitative characters" The American Naturalist (2020) 10.1086/706339

[6]

Butler "Phylogenetic comparative analysis: a modeling approach for adaptive evolution" The American Naturalist (2004) 10.1086/426002

[7]

Caetano "ratematrix: an R package for studying evolutionary integration among several traits on phylogenetic trees" Methods in Ecology and Evolution (2017) 10.1111/2041-210x.12826

[8]

Chasalow "Combinat: combinatorics utilities. R package version 0.0-8" (2012)

[9]

Cox "Notes on option pricing I: constant elasticity of diffusions" (1975)

[10]

Eastman "A noval comparative method for identifying shifts in the rate of character evolution on trees" Evolution (2011) 10.1111/j.1558-5646.2011.01401.x

[11]

A Leisurely Look at the Bootstrap, the Jackknife, and Cross-Validation

Bradley Efron, Gail Gong

The American Statistician 1983 10.1080/00031305.1983.10483087

[12]

Evans "Climate, niche evolution, and diversification of the bird-cage evening primroses (Oenothera, sections Anogra and Kleinia)" The American Naturalist (2009) 10.1086/595757

[13]

Fatone "The use of statistical tests to calibrate the normal SABR model" Journal of Inverse and Ill-Posed Problems (2013) 10.1515/jip-2012-0093

[14]

Fatone "The calibration of some stochastic volatility models used in mathematical finance" Open Journal of Applied Sciences (2014) 10.4236/ojapps.2014.42004

[15]

Phylogenies and the Comparative Method

Joseph Felsenstein

The American Naturalist 1985 10.1086/284325

[16]

Felsenstein (2004)

[17]

Garamszegi (2014) 10.1007/978-3-662-43550-2

[18]

Garland "Procedures for the analysis of comparative data using phylogenetically independent contrasts" Systematic Biology (1992) 10.1093/sysbio/41.1.18

[19]

Gilbert "numDeriv: Accurate Numerical Derivatives. R package version 2016.8-1.1" (2019)

[20]

Goulet "expm: Matrix Exponential, Log, ‘etc’. R package version 0.999-6" (2021)

[21]

Hagan "Managing smile risk" Wilmott Magazine (2002)

[22]

STABILIZING SELECTION AND THE COMPARATIVE ANALYSIS OF ADAPTATION

Thomas F. Hansen

Evolution 1997 10.1111/j.1558-5646.1997.tb01457.x

[23]

Harmon (2019)

[24]

Harvey (1996) 10.1093/oso/9780198549857.001.0001

[25]

Lemon "Plotrix: a package in the red light district of r" R-News (2006)

[26]

LeRoy "Plant phylogenetic history explains in-stream decomposition at a global scale" Journal of Ecology (2019) 10.1111/1365-2745.13262

[27]

Lewis "A likelihood approach to estimating phylogeny from discrete morphological character data" Systematic Biology (2001) 10.1080/106351501753462876

[28]

Ligges "Scatterplot3d—an R package for visualizing multivariate data" Journal of Statistical Software (2003) 10.18637/jss.v008.i11

[29]

Martin "The cryptic origins of evolutionary novelty: 1000-fold faster trophic diversification rates without increased ecological opportunity or hybrid swarm" Evolution (2016) 10.1111/evo.13046

[30]

Meade "BayesTraits V3.0.5" (2021)

[31]

Nunn (2011) 10.7208/chicago/9780226090009.001.0001

[32]

Evolutionary Inferences from Phylogenies: A Review of Methods

Brian C. O'Meara

Annual Review of Ecology, Evolution, and Systemati... 2012 10.1146/annurev-ecolsys-110411-160331

[33]

O’Meara "Testing for different rates of continuous trait evolution using likelihood" Evolution (2006) 10.1111/j.0014-3820.2006.tb01171.x

[34]

ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R

Emmanuel Paradis, Klaus Schliep

Bioinformatics 2019 10.1093/bioinformatics/bty633

[35]

geiger v2.0: an expanded suite of methods for fitting macroevolutionary models to phylogenetic trees

Matthew W. Pennell, Jonathan M. Eastman, Graham J. Slater et al.

Bioinformatics 2014 10.1093/bioinformatics/btu181

[36]

Philippe "Phylogenomics" Annual Review of Ecology, Evolution, and Systematics (2005) 10.1146/annurev.ecolsys.35.112202.130205

[37]

Pinheiro "nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1-152" (2021)

[38]

Plummer "Coda: convergence diagnosis and output analysis for MCMC" R News (2006)

[39]

Puttick "High rates of evolution preceded the origin of birds" Evolution (2014) 10.1111/evo.12363

[40]

Qiu "clusterGeneration: Random Cluster Generation (with Specified Degree of Separation). R package version 1.3.7" (2020)

[41]

R Development Core Team (2020)

[42]

Rabosky "Automatic detection of key innovations, rate shifts, and diversity-dependence on phylogenetic trees" PLOS ONE (2014) 10.1371/journal.pone.0089543

[43]

Rabosky "Analysis and visualization of complex macroevolutionary dynamics: an example from Australian scincid lizards" Systematic Biology (2014) 10.1093/sysbio/syu025

[44]

phytools: an R package for phylogenetic comparative biology (and other things)

Liam J. Revell

Methods in Ecology and Evolution 2012 10.1111/j.2041-210x.2011.00169.x

[45]

Revell "Two new graphical methods for mapping trait evolution on phylogenies" Methods in Ecology and Evolution (2013) 10.1111/2041-210x.12066

[46]

Revell "Phylogenetic analysis of the evolutionary correlation using likelihood" Evolution (2009) 10.1111/j.1558-5646.2009.00616.x

[47]

Revell "Testing quantitative genetic hypotheses about the evolutionary rate matrix for continuous characters" Evolutionary Ecology Research (2008)

[48]

Revell "Phylogenetic signal, evolutionary process, and rate" Systematic Biology (2008) 10.1080/10635150802302427

[49]

Revell "A new phylogenetic method for identifying exceptional phenotypic diversification" Evolution (2012) 10.1111/j.1558-5646.2011.01435.x

[50]

Rohlf "Comparative methods for the analysis of continuous variables: geometric interpretations" Evolution (2001) 10.1111/j.0014-3820.2001.tb00731.x

Showing 50 of 63 references

Cited By

27

Estimating the rate of quantitative trait evolution in the presence of gene tree discordance by calculating likelihoods across trees

Yu K Mo, Matthew W Hahn · 2026

Evolutionary Journal of the Linnean...

Modeling the Evolution of Rates of Continuous Trait Evolution

Bruce S Martin, Gideon S Bradburd · 2022

Systematic Biology

A simple hierarchical model for heterogeneity in the evolutionary correlation on a phylogenetic tree

Liam J. Revell, Ken S. Toyama · 2022

PeerJ

Metrics

27

Citations

63

References

Details

Published: Aug 17, 2021
Vol/Issue: 9
Pages: e11997
License: View

Authors

L

Liam J. Revell

Department of Biology, University of Massachusetts Boston, Boston, MA, USA; Facultad de Ciencias, Universidad Católica de la Santísima Concepción, Concepción, Chile

Funding

National Science Foundation Award: DBI‐1759940

FONDECYT, Chile Award: 1201869

Cite This Article

Liam J. Revell (2021). A variable-rate quantitative trait evolution model using penalized-likelihood. PeerJ, 9, e11997. https://doi.org/10.7717/peerj.11997

A variable-rate quantitative trait evolution model using penalized-likelihood

You May Also Like