Abstract:Background: Continuous traits evolution of a group of taxa that are correlated through a phylogenetic tree is commonly modelled using parametric stochastic differential equations to represent deterministic change of trait through time, while incorporating noises that represent different unobservable evolutionary pressures. Often times, a heterogeneous Gaussian process that consists of multiple parametric sub-processes is often used when the observed data come from a very diverse set of taxa. In the maximum likelihood setting, challenges can be found when exploring the involved likelihood surface and when interpreting the uncertainty around the parameters. Results: We extend the methods to tackle inference problems for mixed Gaussian phylogenetic models (MGPMs) by implementing a Bayesian scheme that can take into account biologically relevant priors. The posterior inference method is based on the Population Monte Carlo (PMC) algorithm that are easily parallelized, and using an efficient algorithm to calculate the likelihood of phylogenetically correlated observations. A model evaluation method that is based on the proximity of the posterior predictive distribution to the observed data is also implemented. Simulation study is done to test the inference and evaluation capability of the method. Finally, we test our method on a real-world dataset. Conclusion: We implement the method in the R package bgphy, available at <a class="link-external link-http" href="http://github.com/bayubeta/bgphy" rel="external noopener nofollow">this http URL</a>. Simulation study demonstrates that the method is able to infer parameters and evaluate models properly, while its implementation on the real-world dataset indicates that a carefully selected model of evolution based on naturally occurring classifications results in a better fit to the observed data.

Efficient Bayesian species tree inference under the multi-species coalescent

Bayesian Inference of Species Trees from Multilocus Data

Inferring Species Trees Directly from Biallelic Genetic Markers: Bypassing Gene Trees in a Full Coalescent Analysis

Computing the probability of gene trees concordant with the species tree in the multispecies coalescent

Coalescent-based species tree estimation: a stochastic Farris transform

Identifiability of speciation times under the multispecies coalescent

Approximate Bayesian computation for Markovian binary trees in phylogenetics

An Annealed Sequential Monte Carlo Method for Bayesian Phylogenetics

PhyloAcc-GT: A Bayesian method for inferring patterns of substitution rate shifts on targeted lineages accounting for gene tree discordance

Efficient Bayesian Inference of General Gaussian Models on Large Phylogenetic Trees

ASTRAL: genome-scale coalescent-based species tree estimation

Computational Performance and Statistical Accuracy of *BEAST and Comparisons with Other Methods

An Efficient Coalescent Epoch Model for Bayesian Phylogenetic Inference

A simulation approach for change-points on phylogenetic trees

Bayesian inference of phylogenetic distances: revisiting the eigenvalue approach

A Variational Approach to Bayesian Phylogenetic Inference

The Multilocus Multispecies Coalescent: A Flexible New Model of Gene Family Evolution

Accelerating Bayesian inference of dependency between complex biological traits

Bayesian inference of mixed Gaussian phylogenetic models

MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space

Sequential importance sampling for multi-resolution Kingman-Tajima coalescent counting