Substitution Model¶
A phylogenetic probability model that assigns rates or probabilities to changes among a finite set of biological character states along evolutionary branches, converting branch duration into transition probabilities for likelihood, distance, simulation, and ancestral inference.
Core Idea¶
A substitution model is a probability model for changes among a finite set of biological character states along evolutionary lineages. The states may be nucleotides, amino acids, codons, or coded morphological character values. The model specifies how an ancestral state at one end of a phylogenetic branch becomes a descendant state at the other end as evolutionary opportunity accumulates. Combined with a phylogenetic tree, branch lengths, and a distribution at the root, it yields probabilities for the character patterns observed at the tips.
Scope of Application¶
The home scope is molecular evolution and statistical phylogenetics. Nucleotide models range from the equal-rate, equal-frequency Jukes–Cantor model through transition/transversion-sensitive and unequal-frequency models to the general time-reversible family. Amino-acid models use a twenty-state alphabet, often estimating exchangeabilities empirically from large protein collections. Codon models use sense codons and can parameterize synonymous versus nonsynonymous change, transition bias, codon frequencies, and lineage- or site-specific selective regimes.
Clarity¶
The abstraction clarifies what “the sequence evolved by this tree” means probabilistically. A topology alone says which lineages connect; it does not say how likely an \(A\to G\) change is, how repeated hits accumulate, or what a branch length measures. The substitution model supplies that missing law.
Manages Complexity¶
The Markov transition law compresses an enormous space of possible histories. Even one site may change many times on many branches, but likelihood calculation need not enumerate every continuous-time event history. Matrix exponentiation integrates over the unobserved substitutions on each branch, and the pruning algorithm sums over possible ancestral states at internal nodes.
Abstract Reasoning¶
The simplest distance-correction example shows why the abstraction is needed. Under JC69, if \(p\) is the observed fraction of differing nucleotide sites between two sequences, the estimated evolutionary distance is
Knowledge Transfer¶
Within evolutionary biology, the same abstraction transfers literally across nucleotide, amino-acid, codon, binary, and morphological alphabets. The matrix dimension, parameterization, and biological interpretation change, but the roles remain: finite states, branch duration, conditional transition probabilities, composition over a tree, and an observation likelihood.
Transfer also occurs across inference tasks. The kernel used to estimate a tree can simulate alignments for method validation, reconstruct ancestral states on a fixed tree, calculate corrected distances, and form null distributions for goodness-of-fit tests. That recurrence is stronger than merely sharing vocabulary.
Relationships to Other Abstractions¶
Current abstraction Substitution Model Domain-specific
Parents (1) — more general patterns this builds on
-
Substitution Model is a kind of Markov Process Prime
Substitution Model strictly specializes Markov Process.
Hierarchy paths (4) — routes to 4 parentless roots
- Substitution Model → Markov Process → Stochastic Process
- Substitution Model → Markov Process → State and State Transition → Phase Space
- Substitution Model → Markov Process → Probability → Measure → Set and Membership
- Substitution Model → Markov Process → Probability → Measure → Aggregation → Micro Macro Linkage
Neighborhood in Abstraction Space¶
Substitution Model sits in a sparse region of the domain-specific corpus (86th percentile for distinctiveness): few abstractions share its structure, so a faithful description tends to retrieve it precisely.
Family — Unclustered & Miscellaneous (1565 abstractions)
Nearest neighbors
- Clade — 0.81
- Polyphyly — 0.80
- Kleene–Brouwer Order — 0.80
- Phylogenetic bracketing — 0.80
- Balding–Nichols Model — 0.79
Computed from structural-signature embeddings · 2026-09-08