MrAIC is an open source and command line based Perl script for calculating AIC, AICc, BIC, and Akaike weights for nucleotide substitution models. Likelihood scores under different models are estimated using PHYML (Guindon and Gascuel, 2003). Input is DNA data in PHYML format.
If the argument -modeltest is parsed, 56 models (the ones tested in Modeltest [Posada and Crandall, 1998]) are evaluated in PHYML. Default is to test the 24 models that can be specified in MrBayes v3. These are JC, F81, K2P (aka K80), HKY, SYM and GTR, each combined with Propinv (I) and/or Gamma (G).
A difference between Modeltest and MrAIC.pl is that MrAIC.pl does not evaluate all models on the same, approximate topology. Instead, PHYML is used to try to find the maximum of the likelihood function under all models. This is necessary for finding AIC, AICc, or BIC for the models.
Requirements:
· Perl
What`s New in This Release: [ read full changelog ]
· The output from PHYML (>v.3.) changed. Commands adjusted to work with the new PHYML version. Thanks to Liliana Davalos.