From Visit Records to a Likelihood
When a process is observed only at scheduled visits, the data usually give a sequence of states rather than exact transition times. For example, a participant might be in state A at baseline, state B at the next visit, and state C at a later visit. A likelihood for these data must describe the probability of that observed sequence, not assign an exact time to an unobserved change of state.
The key construction is local: for each pair of successive visits, use the transition probability from the state recorded at the first visit to the state recorded at the second. Multiplying these conditional probabilities gives the likelihood contribution for the observed history, provided the state process satisfies the stated Markov assumption. The same construction works when different intervals have different transition matrices.
The matrices may be identical across visits when the transition law is unchanged and the intervals have the same relevant length. More generally, \(P_j\) can vary with the interval, such as when visits are separated by different amounts of time or when the model allows calendar-time variation. The definition is conditional on the interval information used by the model; it does not imply that all intervals share one matrix.
Deriving the Likelihood Contribution
For one subject with observed sequence \(x_0,x_1,\ldots,x_m\), begin by conditioning on the observed initial state \(x_0\). The chain rule expresses the probability of the whole sequence as a product of successive conditional probabilities. The Markov assumption then reduces each conditional probability to the appropriate transition-matrix entry.
The first equality is the chain rule and does not require a Markov model. The second equality uses the Markov assumption, which is the substantive modeling step. The final equality follows from the definition of \(P_j\). If the entries depend on an unknown parameter vector \(\theta\), this sequence probability is the subject’s likelihood contribution \(L(\theta)\), conditional on the observed initial state and the visit schedule.
If the initial state is not conditioned on and its distribution is part of the model, its probability must also be included. If \(\pi_a=P\{X(t_0)=a\}\), the probability of the full observed sequence is \(\pi_{x_0}\prod_{j=1}^{m}(P_j)_{x_{j-1},x_j}\). Omitting this initial-state factor is appropriate when the analysis conditions on the baseline state; it is not appropriate when the model claims to describe the joint distribution including that state.
Worked Example: One Three-State History
Suppose three states are labeled A, B, and C. A hypothetical study records one subject in A at baseline, B at visit 1, and C at visit 2. The two visit intervals use the same transition matrix:
Rows and columns are ordered A, B, C. Each row sums to 1; for example, the A row gives \(0.70+0.20+0.10=1.00\). Thus the matrix is a valid transition-probability matrix. Conditional on the observed baseline state A, the probability of the recorded sequence A, B, C is
The sequence contributes 0.0500 to the likelihood, conditional on starting in A. It is not the probability of a particular transition time or of exactly two transitions: the observations specify only the states at the visits.
Combining Intervals and Subjects
For several subjects whose histories are independent under the model, the joint likelihood is the product of their individual likelihood contributions. This independence is ordinarily understood conditional on modeled covariates and any shared parameters. If all subjects share one transition matrix, observed transitions can also be summarized by counts: each matrix entry is raised to the number of observed interval transitions of that type. When matrices differ by interval, the corresponding interval-specific entries must remain distinct.
Worked Example: Two Subjects With a Common Matrix
In a hypothetical panel study, two independent subjects are observed at three equally spaced visits. Both are modeled using the matrix below. Subject 1 has history A, A, B; subject 2 has history B, C, C.
For subject 1, the A-to-A and A-to-B entries give \(L_1=(0.80)(0.15)=0.120\). For subject 2, the B-to-C and C-to-C entries give \(L_2=(0.20)(0.80)=0.160\). Independence gives the joint likelihood:
Equivalently, the four observed transitions consist of one A-to-A, one A-to-B, one B-to-C, and one C-to-C transition, so \(L=P_{AA}P_{AB}P_{BC}P_{CC}\). The factorization describes the probability of these particular labeled histories; it does not require adding a combinatorial coefficient.
A useful general form makes the role of transition counts explicit. Let \(N_{j,ab}\) be the number of subjects observed in state \(a\) at the start of interval \(j\) and state \(b\) at its end. For conditionally independent subjects with common interval-specific matrices, the likelihood conditional on their observed initial states is
This is simply the individual-history product rearranged by matrix entry. A zero count contributes a factor of one. A positive count paired with a zero transition probability makes the likelihood zero, because the model assigns probability zero to an observed transition.
Worked Example: Different Matrices Across Visits
Consider a hypothetical subject observed in B, then C, then A. The first and second intervals have different transition matrices. Suppose the relevant entries are \( (P_1)_{BC}=0.30\) and \( (P_2)_{CA}=0.20\). The matrices are assumed valid, with nonnegative entries and rows summing to one; only the entries needed for this history are shown here.
Conditional on the observed starting state B, the history likelihood is
Using \(P_1\) for both transitions would instead give the wrong expression unless the model specifically states that the transition law is the same over both intervals. The visit-specific index matters: the first factor describes B-to-C over interval 1, and the second describes C-to-A over interval 2.
Assumptions, Checks, and Alternatives
The product formula is simple, but its interpretation depends on several conditions. The states must be defined consistently at every visit, and each recorded transition must correspond to the correct pair of successive visits. The transition matrices must describe the conditional probabilities for the relevant intervals and covariates. The Markov assumption must be defensible at the observation times, and the subject-level likelihood product requires independence across subjects conditional on the model.
The visit schedule and observation process also matter. The likelihood above is conditional on the stated visit times and observed states. If visit attendance, dropout, or measurement depends on unobserved state history in a way related to the parameters of interest, simply conditioning on the observed schedule may not capture the full data-generating process. Such a setting may require modeling the observation process as well; the basic transition product alone does not make informative observation harmless.
Useful checks include verifying that every row of each proposed matrix sums to one, confirming that the row index is the state before the interval and the column index the state after it, and checking that the selected matrix corresponds to the correct interval. Compare observed transitions with fitted probabilities, especially for states or transition types with sparse observations. If the same current state is followed by systematically different next-state patterns depending on earlier observed states, that is evidence against the assumed Markov structure. A zero fitted probability for a transition that appears in the data is a direct warning that the model or parameter restrictions are incompatible with those observations.
This likelihood is not an event-time density. A method for interval-censored event times targets probabilities associated with event-time intervals; here the observed data are a sequence of states, and the factors are transition probabilities between visits. A direct model for state occupation at a particular visit can target the marginal distribution at that visit, but does not by itself represent the joint probability of a complete observed path. Continuous-time transition models can impose additional structure on how probabilities vary with interval length; the transition-matrix likelihood presented here does not require that particular construction.
Why the Construction Matters
Writing the likelihood as a product makes the data-model connection transparent. Each observed step contributes exactly one probability, while the assumptions determine which probability is appropriate. This structure supports estimation by maximizing the likelihood over the transition probabilities or their model parameters. It also makes it easier to identify where a discrepancy originates: a particular interval, a particular starting state, or a particular observed destination.
The same clarity prevents overinterpretation. A panel history shows that the state differed between two visits, but it does not reveal when within the interval the change occurred. The transition-matrix likelihood uses the information actually recorded, without silently treating visit times as transition times.
Check Your Understanding
Answer each question using the transition-matrix likelihood construction.
- A subject is observed in A, B, B across three visits. Write the likelihood conditional on the initial state, allowing separate matrices for the two intervals.
- Why does a likelihood conditioned on the initial state omit the probability of that state?
- Suppose \(P_{AB}=0.25\) and \(P_{BC}=0.40\) in a common matrix. What is the conditional likelihood of the observed history A, B, C?
- What does a positive count of an observed transition with fitted probability zero imply for the likelihood?
Solutions. (1) The contribution is \((P_1)_{AB}(P_2)_{BB}\), because each factor corresponds to one successive visit pair. (2) Conditioning defines the probability distribution among histories given the starting state, so the probability of that starting state is not part of the conditional likelihood; include it only when modeling the joint distribution of the initial state and later observations. (3) \(L=(0.25)(0.40)=0.100\), conditional on starting in A. (4) The observed data receive likelihood zero under that model, indicating incompatibility between the fitted transition restrictions and the observed transition.