Rank substitution models, with the margin over the runners-up
select_substitution_modelCompare substitution models for aligned DNA or protein sequences, returning ranked models by AIC/AICc/BIC with the winner's margin and agreement between criteria.
Instructions
Compare substitution models and report how much the winner won by.
A single model name reads as a finding. The ranking, the delta to the next model, and whether AIC/AICc/BIC agree are what make it one.
Args:
fasta: Aligned sequences in FASTA, nucleotide or protein as declared by
sequence_type.
criterion: "AIC", "AICc" or "BIC". BIC penalises parameters more heavily.
seed: Fixes the engine's search.
top_n: How many ranked models to return.
sequence_type: "dna" (default) or "protein". Ranks within that molecule
type's model set — nucleotide and protein models are not comparable.
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| seed | No | ||
| fasta | Yes | ||
| top_n | No | ||
| criterion | No | AIC | |
| sequence_type | No | dna |
Output Schema
| Name | Required | Description | Default |
|---|---|---|---|
| seed | Yes | ||
| ranking | Yes | ||
| warnings | Yes | ||
| alignment | Yes | ||
| criterion | Yes | ||
| best_model | Yes | ||
| criteria_agree | Yes | ||
| best_by_criterion | Yes | ||
| n_models_compared | Yes | ||
| indistinguishable_from_best | Yes |