Tutorials › Competing Risks and Multi-State Models › History Dependence in Illness-Death Data

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

History Dependence in Illness-Death Data

Use illness-death examples to see when people in the same current state can have different future outcomes because their histories differ.

Intermediate 11 min read

What You'll Learn

  • Define history dependence for future transitions in illness-death data.
  • Distinguish effects of illness duration from changes over calendar time.
  • Calculate and compare short-term outcomes for groups with different histories.
  • Explain how changing history composition can change a pooled state-specific risk.
  • Recognize limitations of crude comparisons and common Markov-model pitfalls.

Why the Same State May Not Mean the Same Future

A state records a person’s current status, but it may not record how that status was reached. Two people who are both ill at a particular time might differ in how long they have been ill, how many previous episodes they have had, or what treatment they have received. If these differences predict what happens next, the current state alone does not summarize the information relevant to their future.

This tutorial focuses on that pattern in illness-death data. The central question is not whether a person’s history is associated with an outcome in general, but whether people with the same current state at the same time have different future outcomes depending on their earlier histories. Fixing the observation time matters: a transition probability can change over calendar time in a time-inhomogeneous model without depending on earlier history.

Definition: Let \(X(t)\) denote a person’s state at time \(t\), and let \(H_t\) denote the observed state and event history up to time \(t\). There is history dependence relative to the current-state model when the conditional distribution of future states, given \(X(t)\) and time \(t\), still differs according to \(H_t\). For a possible transition from state \(i\) to state \(j\), this can be expressed as a transition intensity \(\lambda_{ij}(t\mid H_t)\) that varies between histories having the same current state \(X(t)=i\) and the same time \(t\).

In the basic illness-death structure, state 0 denotes being free of the illness, state 1 denotes being ill, and state 2 denotes death, which is absorbing. The usual transitions are \(0\to1\), \(0\to2\), and \(1\to2\). Extensions can include recovery or repeated episodes. The issue here is not which transitions are allowed; it is whether the current state and time are enough to describe their future occurrence.

For a time-inhomogeneous Markov model, the future conditional on the past depends on the current state and the current time, not on the earlier path. In intensity notation, this means that, once \(X(t)=i\) and \(t\) are fixed, earlier histories do not further change \(\lambda_{ij}\). History dependence is a failure of that state-only condition. It does not, by itself, identify the cause of the difference or show that history has a causal effect.

Three Ways Histories Can Differ

The following examples use hypothetical data to illustrate comparisons within the same current illness state. In each example, suppose all 40 people in each group are in state 1 at the stated landmark, are followed for the same horizon, and have complete follow-up. The displayed proportions are descriptive probabilities over that horizon in the hypothetical group; they are not fitted transition-intensity estimates. For simplicity, assume the illness-death structure has no recovery during the follow-up window. If recovery is possible, the outcome would need to be defined to account for that competing transition.

Worked Example: Duration of the Current Illness

At a common landmark date, two hypothetical groups are both currently ill. Group A has been in the illness state for 6 months; Group B has been in it for 3 years. During the following 12 months, 6 of 40 people in Group A die, compared with 18 of 40 in Group B. Calculate and compare the observed 12-month death proportions.

For Group A,

$$ \hat{p}_A=\frac{6}{40}=0.15 $$

For Group B,

$$ \hat{p}_B=\frac{18}{40}=0.45 $$

The difference is \(0.45-0.15=0.30\), or 30 percentage points. In this hypothetical sample, the group with longer recorded illness duration had the higher observed 12-month death proportion. Because the groups share a current state and landmark time, the contrast is consistent with duration-related history dependence. It does not establish that duration itself caused the difference: the groups could also differ in severity, prior care, or other characteristics.

Worked Example: A Previous Illness Episode

Consider a hypothetical illness that can remit and later recur. At one landmark, 40 people are in an active illness episode for the first time, and 40 are in an active episode after at least one earlier episode. All are currently in state 1. Over the next 6 months, 4 people in the first-episode group and 14 people in the recurrent-episode group die while still in the illness state. Calculate the two observed proportions and their difference.

For the first-episode group, the observed proportion is \(4/40=0.10\). For the recurrent-episode group, it is \(14/40=0.35\). Thus the difference is \(0.35-0.10=0.25\), or 25 percentage points. The count of previous episodes distinguishes the histories even though both groups are currently ill.

The comparison is about an observed association, not a causal effect of recurrence. For example, people with recurrent illness may have more severe disease or longer overall illness histories. The example also defines the event specifically as death while in state 1; if people can recover during follow-up, that transition must be recorded and incorporated rather than silently treated as continued illness.

Worked Example: Earlier Treatment Exposure

In a hypothetical study, 40 people currently in state 1 received a specified treatment before the landmark, while 40 people in the same state did not receive it. During the following 90 days, 5 people in the previously treated group and 13 in the untreated group die. Suppose follow-up is complete and the event is death from the current illness state. Calculate the observed 90-day proportions and interpret the contrast.

The treated group’s proportion is \(5/40=0.125\), or 12.5%. The untreated group’s proportion is \(13/40=0.325\), or 32.5%. The observed difference is \(0.325-0.125=0.200\), or 20 percentage points, with the higher proportion in the untreated group.

This difference could reflect treatment history, but it could also reflect who received treatment, when treatment began, or how severe their illness was. In particular, treatment may have been given more often to people with worse illness, or treatment recipients may have differed in other ways. The data as stated do not justify concluding that treatment reduced mortality. They do illustrate that a state-only description can omit information associated with subsequent outcomes.

Why History Composition Changes a Pooled Risk

A simple calculation shows why history dependence matters even when the goal is only to describe outcomes for everyone currently ill. Suppose the two duration groups in the first example have observed 12-month death proportions of 0.15 and 0.45. If the groups are equally represented, the pooled proportion is 0.30. But if the mix of durations changes, the pooled proportion changes too, even when the group-specific proportions stay the same.

$$ P(\text{death by 12 months}\mid X(t)=1) =\sum_h P(\text{death by 12 months}\mid X(t)=1,H_t=h)\, P(H_t=h\mid X(t)=1) $$

For an equal mixture of short- and long-duration histories, the calculation is

$$ 0.50(0.15)+0.50(0.45)=0.30. $$

If 80% of the currently ill group has the longer-duration history and 20% the shorter-duration history, then

$$ 0.20(0.15)+0.80(0.45)=0.39. $$

The pooled proportion rises from 0.30 to 0.39 because the history mix changes; the two group-specific proportions have not changed. This is a mixture identity, not evidence that one history caused another outcome. It explains why a prediction based only on current state can be inadequate when the distribution of earlier histories differs across populations or times.

The identity follows by partitioning the group currently in state 1 into history categories \(h\). The probability of the outcome for the whole group is the sum, over those categories, of the category-specific probability multiplied by the fraction of the group in that category. If the category-specific probabilities differ, changing those fractions changes the pooled probability. If the probabilities are equal across histories, this particular source of variation disappears.

How to Assess the Pattern

A useful descriptive assessment begins by defining the landmark time, current state, future outcome, follow-up horizon, and history feature before comparing groups. For example, “death within 12 months among people in state 1 on the landmark date, by duration in state 1 at that date” is more precise than “death by illness duration.” The groups must be compared at the same time scale and under a clear rule for transitions such as recovery.

A practical sequence is:

1
Fix current state and time.
Choose a landmark and include people who occupy the same state at that landmark.
2
Specify the earlier-history contrast.
Record a defined feature, such as elapsed illness duration, number of prior episodes, or prior treatment exposure.
3
Define the future outcome.
State the transition and horizon, and account for recovery or other competing transitions if they can occur.
4
Compare and interpret.
Report group-specific outcomes and their contrast, then consider whether other differences or follow-up issues could explain the pattern.

If histories with the same current state show systematically different future transition patterns, that is evidence against the adequacy of a state-only Markov description for the question and data being considered. Relevant checks include whether the apparent contrast persists at multiple landmark times and whether groups have comparable measurement and follow-up. Sparse observations, inaccurate onset dates, or different loss-to-follow-up patterns can make a contrast unstable or misleading. A single small-sample difference is a reason to investigate, not a definitive diagnosis.

Interpretation, Alternatives, and Limitations

A state-only Markov model treats people in the same state at the same time as having the same transition behavior, regardless of how they arrived there. When duration, recurrence, or treatment history is informative, that assumption may be too restrictive. One response is to describe or model the relevant history alongside current state; another is to use a time scale or model structure suited to the feature under study. These alternatives require clear definitions and adequate data. They do not make confounding, measurement error, or incomplete follow-up disappear.

The comparison is not simply between “Markov” and “wrong.” A state-only model may be a useful approximation for a particular prediction target, population, and time horizon even if some history differences exist. Conversely, a strong association between history and outcome does not alone identify the best model. The choice depends on the scientific question, the histories observed, and the predictions or summaries required. The introductory material on Markov assumptions distinguishes dependence on current state from dependence on earlier history; here the illness-death examples make that distinction concrete.

A central pitfall is to compare two groups at different calendar times and attribute their different outcomes to illness duration. Calendar-time changes in care or population composition can produce differences even without duration dependence. Another is to treat crude group proportions as causal treatment effects. A third is to assume that a transition out of illness is impossible because it was not tabulated. Define the state process and outcome explicitly, and distinguish descriptive evidence from causal claims.

Key takeaway: Two people can occupy the same illness state at the same time yet have different future outcomes because their durations, previous episodes, or treatment histories differ. Such differences challenge a state-only summary of the future, but their interpretation requires attention to time, competing transitions, other group differences, and data quality.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. At the same landmark, two groups are both in state 1. One has been ill for 1 year and the other for 4 years. What feature would make a future-outcome difference evidence of history dependence relative to a current-state model?
  2. In a hypothetical 6-month comparison, 9 of 30 people with a first episode and 18 of 60 people with a recurrent episode die. Calculate both observed proportions. Do they differ?
  3. Why can a pooled death proportion change when the mix of short- and long-duration illness histories changes, even if each group’s outcome proportion is stable?
  4. Why does a lower observed death proportion among previously treated people not, by itself, establish a protective causal effect of treatment?

Solutions. 1. At the same calendar time and current state, the future outcome distribution or transition intensity should differ systematically by illness duration; a difference between groups observed at different times would not isolate that feature. 2. The proportions are \(9/30=0.30\) and \(18/60=0.30\); they are equal in this hypothetical comparison, so it shows no difference in these observed proportions. 3. The pooled proportion is a weighted average of the history-specific proportions. Changing the weights changes the average when the history-specific proportions are unequal. 4. Treatment recipients and nonrecipients may differ in severity, eligibility, timing, or other factors that influence outcomes, so the crude contrast is an association and may be confounded.

Clinical Biostats