JMODELTEST(1)

πŸ’– NAME πŸ’–

jmodeltest β€” bestie, this is your HPC-powered glow-up for picking models of nucleotide substitution ✨


✨ SYNOPSIS ✨

jmodeltest -d sequenceFileName [-getPhylip] [-ckp checkpointFileName.ckp] [-n executionName] [-t fixed|BIONJ|ML] [-u userTreeFileName] [-o outputFileName] [-S NNI|SPR|BEST] [-AIC] [-AICc] [-BIC] [-DT] [-c confidenceInterval] [-s 3|5|7|11|203] [-f] [-i] [-g numberOfCategories] [-uLNL] [-dLRT] [-h confidenceInterval] [-hLRT] [-O {ftvwxgp}] [-a] [-z] [-p] [-v] [-w] [-tr numberOfThreads] [-machinesfile machinesFileName]

🌸 DESCRIPTION 🌸

okay so jModelTest is literally the tool you use to do statistical selection of best-fit models of nucleotide substitution. it's giving "which model ate and which one flopped" energy 🫢

it implements five different model selection strategies:

  • hierarchical likelihood ratio tests (hLRT)
  • dynamical likelihood ratio tests (dLRT)
  • Akaike information criterion (AIC)
  • Bayesian information criterion (BIC)
  • a decision theory method (DT)

it also hands you estimates of model selection uncertainty, parameter importances and model-averaged parameter estimates, including model-averaged tree topologies. iconic πŸ‘‘

jModelTest 2 brings High Performance Computing (HPC) capabilities to the party, plus extra features like new strategies for tree optimization, model-averaged phylogenetic trees (both topology and branch length), heuristic filtering and automatic logging of user activity. no cap πŸ”₯


πŸ’… OPTIONS πŸ’…

  • -a β€” estimate model-averaged phylogeny for each active criterion (e.g., -a) (default is false)
  • -AIC β€” calculate the Akaike Information Criterion (e.g., -AIC) (default is false)
  • -AICc β€” calculate the corrected Akaike Information Criterion (e.g., -AICc) (default is false)
  • -BIC β€” calculate the Bayesian Information Criterion (e.g., -BIC) (default is false)
  • -DT β€” calculate the decision theory criterion (e.g., -DT) (default is false)
  • -c confidenceInterval β€” confidence interval (e.g., -c 90) (default is 100)
  • -ckp checkpointFileName β€” loads a checkpointing file
  • -d sequenceFileName β€” input data file (e.g., -d data.phy)
  • -dLRT β€” do dynamical likelihood ratio tests (e.g., -dLRT) (default is false)
  • -f β€” include models with unequals base frecuencies (e.g., -f) (default is false)
  • -g numberOfCategories β€” include models with rate variation among sites and number of categories (e.g., -g 8) (default is false & 4 categories)
  • -G threshold β€” heuristic search. Requires a threshold > 0 (e.g., -G 0.1)
  • -getPhylip β€” converts the input file into phylip format
  • -h confidenceInterval β€” confidence level for the hLRTs (e.g., -a0.002) (default is 0.01)
  • -H informationCriterion β€” information criterion for clustering search (AIC, AICc, BIC). (default is BIC) this argument applies only for 203 substitution schemes (e.g., -s 203 -H AIC)
  • -help β€” displays this help message
  • -hLRT β€” do hierarchical likelihood ratio tests (default is false) hypothesis testing order can be specified with -O argument
  • -i β€” include models with a proportion invariable sites (e.g., -i) (default is false)
  • -machinesfile manchinesFileName β€” gets the processors per host from a machines file
  • -n executionName β€” execution name for appending to the log filenames (default: current time yyyyMMddhhmmss)
  • -o outputFileName β€” set output file (e.g., -o jmodeltest.out)
  • -O hypothesisOrder β€” hypothesis order for the hLRTs (e.g., -hLRT -O gpftv) (default is ftvwxgp)
f=freq, t=titvi, v=2ti4tv(subst=3)/2ti(subst>3), w=2tv, x=4tv,
g=gamma, p=pinv

this argument is used only if -hLRT argument is set. 'f','t','v','g','p' are mandatory in any order. 'w' is optional, and 'x' requires 'w' to be present. thus, length should be 5, 6 (including 'w') or 7 (including both 'w' and 'x') e.g., -hLRT -O gpfvwxt

  • -p β€” calculate parameter importances (e.g., -p) (default is false)
  • -r β€” backward selection for the hLRT (e.g., -r) (default is forward)
  • -s numberOfSubstitutionSchemes β€” number of substitution schemes (e.g., -s 11) (it has to be 3,5,7,11,203; default is 3)
  • --set-local-config localConfigurationFile β€” set a local configuration file in replacement of conf/jmodeltest.conf
  • --set-property propertyName=propertyValue β€” set a new value for a property contained in the configuration file (conf/jmodeltest.conf)
  • -S NNI|SPR|BEST β€” tree topology search operation option (NNI (fast), SPR (a bit slower), BEST (best of NNI and SPR)) (default is BEST)
  • -t fixed|BIONJ|ML β€” base tree for likelihood calculations (e.g., -t BIONJ)
    • fixed β€” (common BIONJ-JC topology)
    • BIONJ β€” (Neighbor-Joining topology)
    • ML β€” (Maximum Likelihood topology) (default)
  • -tr numberOfThreads β€” number of threads to execute (default is 4)
  • -u treeFileName β€” user tree for likelihood calculations (e.g., -u data.tre)
  • -uLnL β€” calculate delta AIC,AICc,BIC against unconstrained likelihood (e.g., -uLnL)

(default is false if the input alignment has gaps or ambiguous characters)

  • -v β€” do model averaging and parameter importances (e.g., -v) (default is false)
  • -w β€” write PAUP block (e.g., -w) (default is false)
  • -z β€” strict consensus type for model-averaged phylogeny (e.g., -z) (default is majority rule)

🌈 EXAMPLE 🌈

jmodeltest -d sequenceFileName -i -f -g 4 -BIC -AIC -AICc -DT -v -a -w

πŸŽ€ AUTHOR πŸŽ€

This manpage was written by Andreas Tille for the Debian distribution and can be used for any other usage of the program.


k that's the whole thing bestie πŸ’– this page is from jmodeltest 2.1.10 and it's dated August 2016. go pick some models, you got this πŸ¦‹βœ¨