In this GitHub repository, you will find all the material that we will be using for the practical sessions focused on Bayesian inference. We will cover the following content:
- Tuesday - PART 1: Intro to MCMC: we will use R to learn more about Markov Chain Monte Carlo (MCMC) with a small dataset analysed under the
K80model. - Tuesday - PART 2: Bayesian phylogeny inference with
PhyloBayes: we will learn how to usePhyloBayes(Lartillot and Philippe, 2004) to infer a phylogeny under a Bayesian approach. - Friday - Bayesian timetree inference with
PAML: we will learn how to useMCMCtree, one of the variousPAMLprograms (Yang 2007), to infer species divergence times. When escalating to phylogenomic data, calculating the exact likelihood calculation during the MCMC may be somewhat prohibitive (computationally very demanding!). To tackle this limitation, we shall learn how to useCODEML(if protein data) orBASEML(if nucleotide data) to infer the branch lengths, the gradient, and the Hessian thatMCMCtreewill subsequently use to approximate the likelihood calculation during the timetree inference (dos Reis and Yang, 2011. Such an approximation makes it possible to analyse phylogenomic data within a reasonable amount of time!
You can click the links you see above to navigate through the content of the practical sessions. Hope you enjoy them and... Happy inference! :)
ⓒ Dr Sandra Álvarez-Carretero | @sabifo4