Practice

  • The phangorn library in R contains functions to perform maximum likelihood based analyses of sequence data
    • While not the fastest program, it is very flexible
    • Can fit a wide variety of substitution models
  • For large-scale phylogenetic inference, there are specialised (external) software packages available
    • RAxML/ExaML
    • PhyML
    • IQTREE
  • As with alignment, we can use R to format the data, and call these external programs
  • There is also an interface to RAxML in the ips package

phangorn

  • Takes data in a special format, phyDat
    • Sequences read in using ape are in DNAbin format
    • There is a convenience function called as.phyDat that can convert between them
  • Can compare between different models using modelTest
  • Can fit models (and trees) using pml

RAxML

  • RAxML is one of the fastest, most accurate phylogeny reconstruction software packages available
  • It does not implement a wide variety of models, instead, it assumes a fairly complex one, which will probably be the best model if the dataset is variable enough

PhyML

  • PhyML is slower than RAxML, but offers a wider variety of models
  • Also offers an interactive (text) interface

IQTREE

  • IQTREE is a fairly new phylogenetic software package
  • It does implement a variety of models, and includes model selection and a fast bootstrap-like procedure