Embeddability and rate identifiability of Kimura 2-parameter matrices
Marta Casanellas1, Jesús Fernández-Sánchez2, Jordi Roca-Lacostena2
1Dpt. Matemàtiques, Universitat Politècnica de Catalunya and BGSMath, Diagonal 647, 08028, Barcelona, Spain. marta.casanellas@upc.edu.
This study characterizes embeddable nucleotide substitution matrices for the K80 model, identifying conditions for rate identifiability in phylogenetic analysis. It reveals a subset of embeddable matrices with non-identifiable rates, impacting parameter estimation.
Area of Science:
- Computational Biology
- Evolutionary Biology
- Mathematical Biology
Background:
- The embeddability of substitution matrices, determining if a Markov process has a continuous-time realization, remains an open problem for 4x4 matrices.
- Understanding embeddability and rate identifiability is crucial for accurate phylogenetic inference and evolutionary modeling.
Purpose of the Study:
- To fully characterize the set of embeddable K80 Markov matrices.
- To identify the subset of embeddable K80 matrices for which evolutionary rates are identifiable.
- To investigate the implications of non-identifiable rates in parameter estimation for phylogenetics.
Main Methods:
- Mathematical analysis of K80 model matrices to define conditions for embeddability.
- Characterization of the parameter space for embeddable and identifiable K80 matrices.
- Computation of relative volumes of these sets to quantify their prevalence.
Main Results:
- Complete characterization of embeddable K80 matrices and those with identifiable rates.
- Identification of an open subset of embeddable matrices with non-identifiable rates, including those with positive eigenvalues and diagonal largest column properties.
- Calculation of the relative volumes of embeddable K80 matrices and embeddable matrices with identifiable rates.
Conclusions:
- The study resolves the embedding problem for the K80 model and its submodels.
- The findings provide critical insights into rate identifiability issues in phylogenetic parameter estimation.
- This work concludes the embedding problem for the K81 model and related submodels.
Related Concept Videos
One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation
On...
Multi-input and Multi-variable systems
In the absence of...
Parameters Affecting Nonlinear Elimination: Zero-Order Input, First-Order Absorption and Two-Compartment Model
When a drug is administered through a constant intravenous infusion and eliminated via nonlinear pharmacokinetics, it follows zero-order input. For example, oral drugs undergo first-order absorption upon administration and are eliminated through nonlinear pharmacokinetics.
In the case of subcutaneously administered drugs,...
Routh-Hurwitz Criterion II
The first scenario occurs when a singular zero appears in the first column of the Routh table. This situation creates a division by zero issues. To resolve this, a small positive or negative number, denoted as epsilon (∈), is substituted for the zero. The stability analysis proceeds by assuming a sign for ∈. If ∈ is positive, any sign change in the first...
Pole and System Stability
Simple poles are unique roots of the denominator polynomial. Each simple pole corresponds to a distinct solution to the system's characteristic equation, typically resulting in exponential decay terms in the system's...
Inertia Tensor
The diagonal components of the inertia tensor matrix represent the moments of inertia concerning the principal axes of the object. These primary axes are defined as the axes where the object experiences the least...


