Tutorials › Competing Risks and Multi-State Models › State Aggregation and Hidden Heterogeneity

States, Transitions and the Markov Assumption · Tutorial 253 of 699

State Aggregation and Hidden Heterogeneity

See how state aggregation can conceal clinically important history and make a broad state’s transition behavior depend on who is currently in it.

Intermediate 12 min read

What You'll Learn

  • Explain why combining scientifically distinct situations can conceal different transition patterns.
  • Calculate an aggregate transition rate from subgroup-specific rates and changing subgroup weights.
  • Distinguish a broad-state risk from risks in the clinically distinct situations it combines.
  • Recognize data patterns that suggest hidden heterogeneity.
  • Compare broad-state summaries with stratification and history-aware descriptions.

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.

Definition: State aggregation occurs when two or more scientifically distinct situations are represented by one model state. Hidden heterogeneity is variation in transition behavior among individuals grouped into that state, often because they differ in prior history or in an unrecorded clinical characteristic.

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:

$$ \lambda_G(t)=w_A(t)\lambda_A+w_B(t)\lambda_B. $$

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:

$$ w_A(t)= \frac{w_A(0)e^{-\lambda_A t}} {w_A(0)e^{-\lambda_A t}+w_B(0)e^{-\lambda_B t}}. $$

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\):

$$ w_A(1)= \frac{0.50e^{-0.20}} {0.50e^{-0.20}+0.50e^{-0.05}} = \frac{0.4094}{0.4094+0.4756} \approx 0.4626. $$

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:

$$ \lambda_G(1) =0.46257(0.20)+0.53743(0.05) =0.1193855 \approx 0.1194\text{ per year}. $$

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:

$$ P(\text{death by 1 day}\mid\text{ward only}) = \frac{0.02}{0.22}(1-e^{-0.22}) \approx 0.01795284 \approx 0.01795. $$

The intensive-care total rate is \(0.035272+0.064728=0.100000\) per day. Its one-day death probability is:

$$ P(\text{death by 1 day}\mid\text{intensive care}) = \frac{0.064728}{0.100000}(1-e^{-0.100000}) \approx 0.0615968. $$

With equal subgroup proportions at entry, the combined one-day probability is the weighted average:

$$ P(\text{death by 1 day}\mid\text{hospitalized at entry}) = 0.50(0.01795284)+0.50(0.0615968) \approx 0.0397748 \approx 0.03977. $$

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:

$$ w_{\text{first}}(1)= \frac{0.50e^{-0.30}} {0.50e^{-0.30}+0.50e^{-0.60}} = \frac{0.740818}{0.740818+0.548812} \approx 0.5745. $$

The share with a prior relapse is therefore \(1-0.5745=0.4255\). The aggregate instantaneous relapse rate is:

$$ \lambda_{\text{relapse}}(1) =0.5745(0.30)+0.4255(0.60) \approx 0.4277\text{ per year}. $$

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.

Common mistake: “Everyone currently labeled hospitalized is in the same state, so the state has one fixed death rate.” This reasoning treats the label as if it guaranteed identical transition behavior. The correction is to ask whether the label combines subgroups with different rates and whether their proportions change. If they do, the aggregate rate can vary with composition even while subgroup rates remain constant.

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.

  1. 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?
  2. 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.
  3. What assumptions justify using \(\frac{\lambda_D}{\lambda_{\text{total}}}(1-e^{-\lambda_{\text{total}}t})\) for a subgroup’s death probability?
  4. 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.

Clinical Biostats