Skip to Main content Skip to Navigation
Journal articles

Computing Bayes Factors using Thermodynamilmc Integration

Nicolas Lartillot 1 Hervé Philippe 2
1 MAB - Méthodes et Algorithmes pour la Bioinformatique
LIRMM - Laboratoire d'Informatique de Robotique et de Microélectronique de Montpellier
Abstract : In the Bayesian paradigm, a common method for comparing two models is to compute the Bayes factor, defined as the ratio of their respective marginal likelihoods. In recent phylogenetic works, the numerical evaluation of marginal likelihoods has often been performed using the harmonic mean estimation procedure. In the present article, we propose to employ another method, based on an analogy with statistical physics, called thermodynamic integration. We describe the method, propose an implementation, and show on two analytical examples that this numerical method yields reliable estimates. In contrast, the harmonic mean estimator leads to a strong overestimation of the marginal likelihood, which is all the more pronounced as the model is higher dimensional. As a result, the harmonic mean estimator systematically favors more parameter-rich models, an artefact that might explain some recent puzzling observations, based on harmonic mean estimates, suggesting that Bayes factors tend to overscore complex models. Finally, we apply our method to the comparison of several alternative models of amino-acid replacement. We confirm our previous observations, indicating that modeling pattern heterogeneity across sites tends to yield better models than standard empirical matrices.
Document type :
Journal articles
Complete list of metadatas

Cited literature [61 references]  Display  Hide  Download

https://hal-lirmm.ccsd.cnrs.fr/lirmm-00102630
Contributor : Christine Carvalho de Matos <>
Submitted on : Tuesday, October 16, 2018 - 10:08:33 AM
Last modification on : Monday, July 20, 2020 - 12:34:52 PM
Document(s) archivé(s) le : Thursday, January 17, 2019 - 1:22:22 PM

File

55-2-195.pdf
Files produced by the author(s)

Identifiers

Collections

Citation

Nicolas Lartillot, Hervé Philippe. Computing Bayes Factors using Thermodynamilmc Integration. Systematic Biology, Oxford University Press (OUP), 2006, 55 (2), pp.195-207. ⟨10.1080/10635150500433722⟩. ⟨lirmm-00102630⟩

Share

Metrics

Record views

175

Files downloads

51