Why Matrix Coherence Matters
A predicted transition matrix is meant to describe a probability distribution over possible states at a future time. Its entries may each look plausible when inspected separately, yet the collection can still be impossible: a probability could be negative, exceed one, or leave too much or too little probability accounted for across destinations. Checking this internal coherence is a basic audit of the prediction itself. It does not establish that the predictions agree with observed outcomes; that is a different question.
The key is to identify precisely what the entries represent. A full transition matrix describes the state occupied at a specified future time, conditional on the state at the starting time and any specified covariates or history. It is not necessarily a table of probabilities for the first transition out of the origin. This distinction matters when a process can move through intermediate states.
The conditioning information \(z\) may include covariates and, where appropriate, observed history. Predictions in a row must refer to the same origin, time interval, covariate values, and conditioning information. Otherwise, even individually valid probabilities may not form a probability distribution over one common prediction target.
Basic Coherence Conditions
For every origin state \(i\) and destination state \(j\), a valid probability satisfies \(0\le p_{ij}(s,t\mid z)\le 1\). In a complete transition matrix, the destination states are mutually exclusive at time \(t\), and exactly one of them is occupied. Therefore, each row must also satisfy
The row sum is the appropriate total because a row conditions on one particular origin state. The entries in that row give the distribution over future destinations for that origin. Adding entries from different rows has no probability interpretation unless a separate distribution over origin states has been specified.
A model may report only a subset of destinations. In that case, the displayed entries need not sum to one, but the omitted probability must be understood and accounted for. For example, if a row includes two specified destinations but omits an “other state” category, their sum can be less than one. A sum greater than one cannot be explained by omitted destinations. For a transparent audit, it is often useful to display the residual probability explicitly:
where \(\mathcal{J}\) is the set of displayed destinations. This residual is valid only if it is nonnegative and genuinely represents the omitted outcomes, rather than unrecorded probability being treated as if it had a known destination.
Worked Example: Auditing a Complete Three-State Matrix
Suppose a hypothetical model predicts state occupancy six months after prediction origin. Its states are A, B, and C, and the full matrix is
For an origin in A, the entries are all in \([0,1]\), and the row sum is \(0.70+0.20+0.10=1.00\). For an origin in B, the sum is \(0.05+0.80+0.15=1.00\). For an origin in C, it is \(0.02+0.08+0.90=1.00\). Thus all three rows meet the basic numerical coherence checks. The first row, for example, predicts a 70% probability of being in A, a 20% probability of being in B, and a 10% probability of being in C at six months, conditional on starting in A and on the specified prediction information.
Why the Conditions Hold
The conditions follow directly from the probability interpretation, not from a particular estimation method. Conditional on the specified origin and prediction information, the events “in state 1 at time \(t\),” through “in state \(K\) at time \(t\)” are mutually exclusive. Under the exhaustive-state assumption, one of them must occur. Addition of probabilities for disjoint events therefore gives
Each summand is a probability, so it lies between zero and one. This justification also identifies the assumptions behind the check: the entries must share the same conditioning information, and the destination states must be mutually exclusive and exhaustive for the target being predicted. If a state is omitted, or if the entries come from different targets, the row-sum identity may not apply to the displayed values as written.
For a prediction of the first transition out of a state, the corresponding events are instead “no outgoing transition by time \(t\)” and “the first outgoing transition by time \(t\) is to destination \(j\),” for each allowed destination \(j\). With a complete list of competing destinations, their probabilities plus the probability of remaining in the origin sum to one. This is the appropriate accounting for competing outgoing transitions. It is not automatically the same as a row of state-occupation probabilities: after leaving the origin, a process may enter an intermediate state and then move again before time \(t\).
Worked Example: Finding an Incoherent Competing-Risk Prediction
In a hypothetical prediction from state A over one year, a model reports a 0.62 probability of remaining in A, a 0.27 probability that the first outgoing transition is to B, and a 0.19 probability that the first outgoing transition is to C. Assume B and C are the complete list of competing first destinations.
The three outcomes are mutually exclusive and exhaustive under those assumptions. Their predicted total is \(0.62+0.27+0.19=1.08\), not 1. The reported probabilities therefore cannot all describe this same first-transition distribution. Each is individually between zero and one, but that does not make the collection coherent.
One can see the conflict from the event mass alone: the predicted probability of leaving A is \(1-0.62=0.38\), whereas the two reported first-transition probabilities add to \(0.27+0.19=0.46\). The outgoing probabilities exceed the available event mass by \(0.46-0.38=0.08\). A sound response is to investigate how the predictions were constructed, including whether the same horizon and conditioning information were used and whether destination probabilities were estimated jointly. Rescaling 0.27 and 0.19 to add to 0.38 would preserve their ratio, but that is an additional modeling choice, not a general statistical correction.
Coherence Across Time Intervals
A matrix may pass every row-sum check and still conflict with the model’s assumptions about transitions over time. Under a Markov model, the distribution after two consecutive intervals must agree with the distribution obtained by composing the interval-specific transition matrices. The Markov assumption here means that, conditional on the state at an intermediate time and the relevant modeled covariates, the future distribution does not additionally depend on earlier history.
when these matrices describe the same process and use compatible covariate and time specifications. The multiplication sums over every possible intermediate state: the probability of ending in \(j\) from \(i\) is the sum, over intermediate states \(k\), of the probability of reaching \(k\) and then reaching \(j\). This is a stronger check than row sums, and it depends on the Markov specification. It should not be imposed blindly when the prediction target or model does not have that property.
Worked Example: Checking Two-Interval Consistency
Consider a hypothetical two-state Markov process, where state 2 is absorbing. For each of two consecutive intervals, the model gives
Both matrices are individually coherent: each entry lies between zero and one and each row sums to one. Their two-interval prediction is
Suppose a separately produced two-interval forecast instead reports \(\begin{pmatrix}0.55&0.45\\0&1\end{pmatrix}\). This direct forecast also has valid entries and row sums of one, but it disagrees with the composed prediction: from state 1, its probability of occupying state 2 is 0.45 rather than 0.51. Under the stated Markov model and interval assumptions, the forecasts are temporally inconsistent. The discrepancy calls for checking model inputs, interval definitions, and implementation; row sums alone would not reveal it.
Practical Audit and Failure Modes
A useful audit proceeds from the definition of the prediction target to the numerical checks. It should be performed for every origin row and relevant prediction time, not just on an overall average of the matrix.
Record the origin state, horizon, covariate values, and whether the entries describe state occupancy at the horizon or the first outgoing transition.
Verify that destination categories are mutually exclusive and exhaustive, or identify explicitly what outcomes are omitted.
Flag any value below zero or above one, calculate each row sum, and retain the residual rather than silently discarding it.
Verify structural zeros, absorbing-state rows, and any cross-time consistency required by the model.
Small discrepancies may arise from rounding displayed values. For instance, a row shown as 0.333, 0.333, and 0.333 sums to 0.999; unrounded values may sum to one. The audit should therefore use unrounded predictions where available and report the display precision. Rounding cannot plausibly explain substantial discrepancies, negative probabilities, or values above one.
Structural constraints provide further checks. If a transition is impossible by definition, its probability should be zero. If a state is absorbing under the specified model, its row should assign probability one to remaining in that state at every horizon. These checks depend on the model’s state definitions and assumptions; they are not universal properties of every state with a particular label.
Coherence is necessary, not sufficient. A matrix can satisfy every numerical constraint and still be poorly estimated or systematically disagree with observed outcomes. Conversely, estimating each destination probability with a separate model can produce values that look plausible individually but fail to form a joint distribution. Joint multinomial or competing-transition approaches can enforce shared probability accounting by construction; constructing a transition matrix from a coherent transition model is another option. Independently fitted probabilities may require a revised model or a carefully justified reconciliation step. Mechanical truncation of negative values, clipping values above one, or normalizing rows changes predictions and does not by itself establish a valid statistical model.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- A complete row for three destination states is \((0.40,0.35,0.20)\). Is it coherent? State the numerical check.
- A model reports first-transition probabilities 0.25 to B and 0.15 to C, with the remaining-in-origin probability 0.60. What is the row sum, and does it pass the basic check?
- A displayed row includes probabilities 0.50 and 0.30 for two destinations, while a third destination is omitted. What residual probability remains, and what must be true for it to be interpretable?
- Under a Markov model, two consecutive interval matrices are both \(\begin{pmatrix}0.8&0.2\\0&1\end{pmatrix}\). What is the two-interval probability of moving from state 1 to state 2?
Solutions. (1) The row is not coherent as a complete distribution: its entries are in \([0,1]\), but \(0.40+0.35+0.20=0.95\), not 1. (2) The sum is \(0.25+0.15+0.60=1.00\), so it passes the basic numerical check, assuming the three outcomes are exhaustive and share the same prediction target. (3) The residual is \(1-(0.50+0.30)=0.20\). It represents the omitted destination only if that destination accounts for all omitted outcomes and the displayed probabilities refer to the same origin, horizon, and conditioning information. (4) Matrix multiplication gives the state-1-to-state-2 probability \(0.8(0.2)+0.2(1)=0.36\).