ResearchTrend.AI
  • Papers
  • Communities
  • Events
  • Blog
  • Pricing
Papers
Communities
Social Events
Terms and Conditions
Pricing
Parameter LabParameter LabTwitterGitHubLinkedInBlueskyYoutube

© 2025 ResearchTrend.AI, All rights reserved.

  1. Home
  2. Papers
  3. 1905.12146
35
31

Gradients do grow on trees: a linear-time O(N){\cal O}\hspace{-0.2em}\left( N \right)O(N)-dimensional gradient for statistical phylogenetics

29 May 2019
X. Ji
Zhenyu Zhang
Andrew J Holbrook
A. Nishimura
G. Baele
A. Rambaut
P. Lemey
M. Suchard
ArXiv (abs)PDFHTML
Abstract

Calculation of the log-likelihood stands as the computational bottleneck for many statistical phylogenetic algorithms. Even worse is its gradient evaluation, often used to target regions of high probability. Order O(N){\cal O}\hspace{-0.2em}\left( N \right)O(N)-dimensional gradient calculations based on the standard pruning algorithm require O(N2){\cal O}\hspace{-0.2em}\left( N^2 \right)O(N2) operations where N is the number of sampled molecular sequences. With the advent of high-throughput sequencing, recent phylogenetic studies have analyzed hundreds to thousands of sequences, with an apparent trend towards even larger data sets as a result of advancing technology. Such large-scale analyses challenge phylogenetic reconstruction by requiring inference on larger sets of process parameters to model the increasing data heterogeneity. To make this tractable, we present a linear-time algorithm for O(N){\cal O}\hspace{-0.2em}\left( N \right)O(N)-dimensional gradient evaluation and apply it to general continuous-time Markov processes of sequence substitution on a phylogenetic tree without a need to assume either stationarity or reversibility. We apply this approach to learn the branch-specific evolutionary rates of three pathogenic viruses: West Nile virus, Dengue virus and Lassa virus. Our proposed algorithm significantly improves inference efficiency with a 126- to 234-fold increase in maximum-likelihood optimization and a 16- to 33-fold computational performance increase in a Bayesian framework.

View on arXiv
Comments on this paper