Why a Broad State Can Hide Important Differences
A multi-state model begins with a scientific decision: what situations count as distinct states? The definition must be useful, but it can also be too broad. For example, “hospitalized” might include both ward-only care and intensive care, while “in remission” might include a first remission and a remission after relapse. Each label describes a current situation, yet the people grouped under it may have different prior histories and different prospects.
The consequence is not merely loss of detail in a diagram. If the combined situations have different transition patterns, the apparent transition rate out of the broad state can depend on the mix of situations currently represented there. That mix can change over time or differ between groups. A single rate attached to the broad label may therefore conceal clinically relevant variation.
Aggregation does not automatically make a model unsuitable. It may be reasonable when the distinctions are irrelevant to the scientific question, cannot be measured reliably, or would make the model impractically complex. The key is to understand what the broad label combines and what that combination means for interpretation.
How Mixtures Produce an Aggregate Transition Rate
Suppose a broad state \(G\) combines two subgroups, \(A\) and \(B\). For simplicity, consider one possible transition out of \(G\), with subgroup-specific instantaneous rates \(\lambda_A\) and \(\lambda_B\). An instantaneous rate describes the event rate at a moment among individuals who are still in the relevant subgroup and state. Let \(w_A(t)\) be the proportion of those currently in \(G\) who belong to \(A\); then \(w_B(t)=1-w_A(t)\).
At time \(t\), the aggregate rate is the weighted average of the subgroup rates:
This expression is a composition-weighted rate, not necessarily a fixed property of the broad state. If \(A\) and \(B\) have different rates, their proportions among those who remain in \(G\) can change as events occur. In that case, \(\lambda_G(t)\) changes even if each subgroup’s own rate stays constant.
The changing weights can be derived directly. Suppose entry to \(G\) occurs at time 0, everyone is initially in \(G\), the initial proportion in \(A\) is \(w_A(0)\), and each subgroup has a constant rate \(\lambda_A\) or \(\lambda_B\) of leaving \(G\). The proportions still in each subgroup at time \(t\), before normalizing, are proportional to \(w_A(0)e^{-\lambda_A t}\) and \(w_B(0)e^{-\lambda_B t}\). Thus:
This follows because \(e^{-\lambda_A t}\) and \(e^{-\lambda_B t}\) are the subgroup-specific probabilities of remaining in \(G\) to time \(t\) under the stated constant-rate model. Normalizing the two remaining subgroup masses gives the displayed conditional proportion. If \(A\) has the higher rate, members of \(A\) tend to leave sooner, so the remaining broad state becomes increasingly weighted toward \(B\).
With competing types of exit, the same principle applies to the rate of a particular transition: weight each subgroup’s transition-specific rate by its current share of the broad state. The weights then reflect survival in the state against all exits, not just the transition of interest. These calculations rely on having well-defined subgroup membership and appropriate subgroup-specific rates; they do not make a broad state homogeneous by themselves.
Worked Example: A Broad State With Changing Composition
Suppose a hypothetical “stable chronic condition” state contains two groups in equal proportions at entry. Group A has a constant rate of leaving the state of \(0.20\) per year, and group B has a constant rate of \(0.05\) per year. Assume exponential time to exit within each group and no other source of entry during this calculation. Find the composition of the broad state after one year and its aggregate instantaneous exit rate then.
At \(t=1\) year, the unnormalized remaining shares are \(0.50e^{-0.20}\) for A and \(0.50e^{-0.05}\) for B. Using \(e^{-0.20}\approx0.8187\) and \(e^{-0.05}\approx0.9512\):
Thus about 46.26% of those still in the broad state at one year are in A, and about 53.74% are in B. The aggregate instantaneous exit rate is:
At entry the aggregate rate was \(0.50(0.20)+0.50(0.05)=0.1250\) per year. It has fallen because the higher-rate group is less common among those who remain. The subgroup rates did not change; the composition did. Reporting only one rate for “stable chronic condition” would hide this pattern.
Broad Clinical Labels and Different Risks
A broad health state can combine situations with different probabilities of particular outcomes, even over a fixed follow-up interval. A useful analysis therefore asks not only “What is the risk in this state?” but also “Which clinical situations contribute to that risk, and in what proportions?” The overall risk is meaningful for a specified population and time, but it may not describe either subgroup well.
Worked Example: “Hospitalized” Does Not Specify One Mortality Risk
Consider a hypothetical cohort in which half of patients in the broad “hospitalized” state are on a ward only, and half are in intensive care. For ward-only patients, suppose the constant competing rates are \(0.20\) per day for discharge and \(0.02\) per day for death. For intensive-care patients, suppose the corresponding rates are \(0.035272\) per day for discharge and \(0.064728\) per day for death. Assume the rates are constant over the one-day interval, each patient is initially in the stated subgroup, and there are no other exits. Calculate the one-day death probability in each subgroup and in the combined population.
For constant competing exponential rates, if the death rate is \(\lambda_D\) and the total exit rate is \(\lambda_{\text{total}}\), the probability of death by time \(t\) is \(\frac{\lambda_D}{\lambda_{\text{total}}}(1-e^{-\lambda_{\text{total}}t})\). The ward-only total rate is \(0.20+0.02=0.22\) per day, so:
The intensive-care total rate is \(0.035272+0.064728=0.100000\) per day. Its one-day death probability is:
With equal subgroup proportions at entry, the combined one-day probability is the weighted average:
The combined risk is about 3.977% over one day, but it is not a good description of either subgroup: the ward-only risk is about 1.795%, while the intensive-care risk is about 6.160%. The assumptions justify the competing-rate calculations within each subgroup. The combined calculation additionally assumes the stated 50:50 composition for the target population. A hospital with a different mix of ward-only and intensive-care patients would generally have a different broad-state risk.
A similar issue arises when a label records current status but omits how that status was reached. “In remission” may combine first remission with remission after relapse and salvage therapy. “Discharged” may combine a first discharge with discharge after readmission. In each case, a reported risk for the combined label depends on the mixture of histories in the population being summarized.
Worked Example: Remission Risk Depends on the Mix of Histories
Suppose a hypothetical “in remission” state combines patients in first remission, with a constant relapse rate of \(0.30\) per year, and patients in remission after a prior relapse, with a constant relapse rate of \(0.60\) per year. Assume equal proportions enter the broad state at time 0, exponential time to relapse within each history group, and no competing exits. What is the aggregate instantaneous relapse rate among those still in remission after one year?
The first-remission share among those still in remission after one year is:
The share with a prior relapse is therefore \(1-0.5745=0.4255\). The aggregate instantaneous relapse rate is:
At entry, the corresponding rate was \(0.50(0.30)+0.50(0.60)=0.4500\) per year. By one year, more of the higher-rate group has relapsed and left the state, so the remaining remission group has a lower aggregate rate. “In remission” does not identify a single relapse rate without also specifying the mixture of histories and the time since entry.
Recognizing Hidden Heterogeneity
A broad state deserves scrutiny when its occupants have different clinical pathways, prognosis, treatment, or time since an important event. The following checks can help identify whether aggregation is masking differences:
- List what the label combines. Ask whether it includes distinct disease stages, care settings, treatment histories, or routes into the state.
- Compare transition patterns across known histories. If sufficiently observed subgroups show different risks or rates, a single broad-state summary will depend on their proportions.
- Check whether the mixture changes over follow-up. Selective exit from the state can alter subgroup shares even when each subgroup’s rate is constant.
- Compare populations with different compositions. A broad-state risk may differ between clinics or treatment groups simply because their subgroup mixes differ. That is not, by itself, evidence that the underlying subgroup-specific rates differ.
- Look for clinical meaning in the omitted distinction. A numerical difference matters especially when it changes interpretation, prediction, or a decision relevant to the research question.
These checks require care. A subgroup difference in observed rates can reflect sparse data, measurement error, or other variables that differ between groups, rather than the recorded history alone. Conversely, failure to detect a difference is not proof that histories are interchangeable, especially when subgroup samples are small. Clinical knowledge and data quality should inform the decision, not a single significance test.
Alternatives, Limitations, and a Common Mistake
The simplest approach is to keep the broad state and report its transition behavior for a clearly described target population and follow-up period. This can be appropriate when the broad average answers the scientific question. When the differences matter, alternatives include reporting subgroup-specific transition summaries, retaining a relevant history variable in the analysis, or defining more specific states. Each option requires adequate information and a clear scientific rationale; splitting states indiscriminately can create sparse data and an unnecessarily complicated model.
A broad-state model can also be useful as a deliberate simplification, provided its limitations are stated. Its average transition behavior may be population-dependent, and applying it to a population with a different mix can mislead. A single average can conceal both high- and low-risk groups, and it may be unsuitable for questions that require individual prediction or comparison of clinical pathways. The introductory material on state occupation versus event history distinguishes current-state summaries from the sequence of events that led to them; aggregation is one reason those histories can remain important.
The central practical question is not whether every possible distinction must become a separate state. It is whether the distinctions being combined are scientifically relevant to the transition questions being asked. State definitions are simplifications of clinical reality; making that simplification explicit helps prevent a broad label from being mistaken for a homogeneous clinical condition.
Check Your Understanding
Answer each question using the definitions and calculations above. Brief solutions follow the questions.
- A broad state contains two groups in equal proportions at entry. Their constant exit rates are \(0.10\) and \(0.30\) per year. Which group will make up a larger share of those still in the state after one year, and why?
- Two clinics report different one-year mortality probabilities for patients labeled “hospitalized.” Give one compositional explanation that does not require different subgroup-specific mortality probabilities.
- What assumptions justify using \(\frac{\lambda_D}{\lambda_{\text{total}}}(1-e^{-\lambda_{\text{total}}t})\) for a subgroup’s death probability?
- At entry, two remission-history groups have relapse rates of \(0.20\) and \(0.50\) per year and equal shares. Under constant rates and no competing exits, should the high-rate group’s share among those remaining in remission rise or fall over time? Explain.
Solutions. (1) The group with rate \(0.10\) per year will make up a larger share. Its remaining proportion is multiplied by \(e^{-0.10}\), which is greater than \(e^{-0.30}\), so fewer members of that group leave by one year. (2) One clinic may have a larger proportion of intensive-care patients, while the other has more ward-only patients. If those subgroups have different risks, their weighted overall risks can differ even when risks within each subgroup are the same across clinics. (3) The model assumes constant cause-specific rates during the interval, competing exponential exit times represented by those rates, and no unmodeled exits that affect the total rate. The formula computes the probability that death is the first exit by time \(t\). (4) The high-rate group’s share should fall: its members leave remission more quickly on average, so they become less represented among those who remain.