Tutorials › Competing Risks and Multi-State Models › Bounds Under Unmeasured Informative Censoring

Practical Handling of Censoring · Tutorial 563 of 699

Bounds Under Unmeasured Informative Censoring

Learn to calculate and interpret risk bounds under unmeasured informative censoring, and see precisely which additional assumptions narrow or widen them.

Advanced 13 min read

What You'll Learn

  • Define a fixed-horizon cause-specific event outcome when follow-up is censored.
  • Derive sharp worst-case bounds without assuming censoring is independent of the event outcome.
  • Calculate how many censored outcomes would have to be events to attain an upper bound.
  • Narrow bounds using explicit restrictions on risk among censored individuals.
  • Distinguish conditional risk-ordering assumptions from marginal ones.
  • Separate identification bounds from sampling uncertainty and point estimates.

Why Bound the Risk?

A fixed-horizon event risk is often the quantity of interest: for example, the probability that a person experiences a particular event by five years. When some people are censored before five years, their event status at that horizon is unknown. If censoring may depend on unmeasured outcomes, the observed people with known status need not represent those whose status is missing. A single risk estimate may then rely on assumptions that the data cannot verify.

Bounds make this uncertainty explicit. They give a range of risks compatible with the observed outcomes and stated assumptions, rather than choosing one value for the unobserved outcomes. The range may be wide; that is informative, not a calculation failure. Its width shows how much the conclusion depends on censoring assumptions.

Definition: Fix a time horizon \(\tau\) and a cause \(k\) in a competing-risks setting. Let \(Y=1\) if cause \(k\) occurs by \(\tau\), and \(Y=0\) otherwise. Let \(R=1\) if \(Y\) is known from observed follow-up, and \(R=0\) if the person is censored before their status through \(\tau\) is known. Write \(q=P(R=1)\), \(p_{\mathrm{obs}}=P(Y=1\mid R=1)\), and \(p_{\mathrm{mis}}=P(Y=1\mid R=0)\). The target risk is \(P(Y=1)\).

For this definition, an observed competing event before \(\tau\) gives \(Y=0\): it rules out cause \(k\) as the first event by the horizon. A person known to remain event-free through \(\tau\) also has \(Y=0\). A person censored before their status through \(\tau\) is known has an unknown \(Y\). The target is the cumulative incidence of cause \(k\) by \(\tau\), not a cause-specific hazard or a hypothetical risk in a world where competing events are removed.

This setup assumes that the target population and horizon are defined, event type and timing are recorded correctly when observed, and censoring is the reason status is unknown. It does not assume that censoring is independent of the unobserved outcome. The bounds below describe what is identified at the population level under these conditions. Their sample analogues are estimates of the bound endpoints and have sampling uncertainty.

Worst-Case Bounds

The observed data identify the outcome risk among people with known status, but not the risk among those censored early. Since a binary risk must lie between zero and one, \(0\leq p_{\mathrm{mis}}\leq1\). The law of total probability gives the target risk as the weighted average of these two risks:

$$ P(Y=1)=q\,p_{\mathrm{obs}}+(1-q)\,p_{\mathrm{mis}}. $$

Setting the unknown risk to its smallest possible value gives the lower bound; setting it to its largest gives the upper bound. These are sharp bounds: without further restrictions, both endpoint assignments are compatible with the observed data.

Theorem: With binary \(Y\), known-outcome fraction \(q\), and no restriction on outcomes among those censored before status is known, the identified interval for the risk is
$$ q\,p_{\mathrm{obs}} \ \leq\ P(Y=1)\ \leq\ q\,p_{\mathrm{obs}}+(1-q). $$

The lower endpoint counts the observed cause-\(k\) events and assumes no censored person has cause \(k\) by \(\tau\). The upper endpoint adds every person with unknown status to those events and assumes all censored people have cause \(k\) by \(\tau\). These are deliberately extreme assignments, not claims that either is plausible.

Worked Example: Worst-Case Bounds With a Competing Event

In a hypothetical cohort of 200 people, status through \(\tau\) is known for 120: 24 have cause \(k\) by \(\tau\), 36 have a competing event first, and 60 remain event-free through \(\tau\). The other 80 are censored before their status through \(\tau\) is known. Calculate the worst-case bounds for the cause-\(k\) risk.

Among people with known status, \(p_{\mathrm{obs}}=24/120=0.20\), and \(q=120/200=0.60\). The competing events and known event-free outcomes count as \(Y=0\). Thus the lower endpoint is

$$ q\,p_{\mathrm{obs}}=0.60(0.20)=0.12. $$

For the upper endpoint, assign \(Y=1\) to all 80 censored people:

$$ q\,p_{\mathrm{obs}}+(1-q) =0.60(0.20)+0.40 =0.12+0.40 =0.52. $$

The worst-case identified interval is therefore 0.12 to 0.52, or 12% to 52% by \(\tau\). Its width is 40 percentage points, exactly the fraction of the cohort whose horizon outcome is unknown. The interval is not a 95% confidence interval; it represents uncertainty about the risk from missing outcomes under the stated assumptions.

How Assumptions Narrow the Bounds

A narrower interval requires information or assumptions restricting \(p_{\mathrm{mis}}\). Any such restriction should be stated in a form that can be examined and varied. For instance, suppose subject-matter knowledge supports that the risk among censored people is within \(\delta\) of the observed-group risk, where \(\delta\) is a chosen nonnegative risk difference. Then

$$ \max(0,p_{\mathrm{obs}}-\delta) \ \leq\ p_{\mathrm{mis}}\ \leq\ \min(1,p_{\mathrm{obs}}+\delta). $$

Substitution into the total-risk identity yields bounds. The restriction is a sensitivity assumption, not a consequence of the observed data. At \(\delta=1\), it imposes no effective restriction and recovers the worst-case interval. At \(\delta=0\), it assumes equal risks in the two groups and identifies the risk as \(p_{\mathrm{obs}}\).

Worked Example: A Risk-Difference Sensitivity Restriction

Use the hypothetical data from the previous example, for which \(q=0.60\) and \(p_{\mathrm{obs}}=0.20\). Suppose investigators consider it plausible that the censored-group risk differs from the observed-group risk by no more than 0.10 in either direction. Then \(\delta=0.10\), so

$$ 0.10=\max(0,0.20-0.10) \ \leq\ p_{\mathrm{mis}}\ \leq\ \min(1,0.20+0.10)=0.30. $$

The total-risk lower endpoint is \(0.60(0.20)+0.40(0.10)=0.12+0.04=0.16\). The upper endpoint is \(0.60(0.20)+0.40(0.30)=0.12+0.12=0.24\). Under this explicit restriction, the risk lies between 16% and 24% by \(\tau\). The narrower range comes from the assumed 10-percentage-point limit, not from additional observed outcomes.

A different restriction is an ordering assumption. For example, \(p_{\mathrm{mis}}\leq p_{\mathrm{obs}}\) says that, on average, censored people are no more likely to have the event by \(\tau\) than people with known status. Combining that assumption with \(p_{\mathrm{mis}}\geq0\) gives

$$ q\,p_{\mathrm{obs}} \ \leq\ P(Y=1)\ \leq\ q\,p_{\mathrm{obs}}+(1-q)p_{\mathrm{obs}} =p_{\mathrm{obs}}. $$

The reverse ordering, \(p_{\mathrm{mis}}\geq p_{\mathrm{obs}}\), instead gives the interval from \(p_{\mathrm{obs}}\) to \(q\,p_{\mathrm{obs}}+(1-q)\). Either ordering may be plausible in a particular setting, but neither follows merely from observing censoring. An assumed direction of selection should not be presented as a finding established by the observed outcomes.

Worked Example: Why Conditional Assumptions Matter

Consider a hypothetical cohort of 200 people divided into two equally sized baseline risk groups. In the lower-risk group, 80 of 100 have known outcomes, including 8 cause-\(k\) events; 20 are censored. In the higher-risk group, 40 of 100 have known outcomes, including 16 cause-\(k\) events; 60 are censored. Suppose, as a sensitivity assumption, that in each group separately the censored-group risk is no greater than that group’s observed-group risk.

The observed-group risks are \(8/80=0.10\) and \(16/40=0.40\). The overall lower bound is the 24 observed cause-\(k\) events divided by 200, or \(0.12\). For the upper bound under the group-specific assumptions, assign censored people the largest allowed risk in their group. The lower-risk group contributes at most \(20(0.10)=2\) additional events; the higher-risk group contributes at most \(60(0.40)=24\). Therefore,

$$ \text{upper bound} =\frac{24+2+24}{200} =\frac{50}{200} =0.25. $$

The resulting interval is 12% to 25%. This calculation illustrates why an assumption must specify whether it is marginal or conditional on covariates: risk ordering within each group does not imply the same ordering after the groups are combined. It also shows that stratifying alone does not make missing outcomes known. The gain here is from the stated group-specific restrictions, not from stratification by itself.

Derivation and Interpretation

The bounds follow directly from partitioning the population into those with known and unknown horizon outcomes. Let \(R\) indicate whether the outcome is known. Then \(Y=1\) occurs either among people with \(R=1\) or among people with \(R=0\), so

$$ \begin{aligned} P(Y=1) &=P(Y=1,R=1)+P(Y=1,R=0)\\ &=P(R=1)P(Y=1\mid R=1) +P(R=0)P(Y=1\mid R=0)\\ &=q\,p_{\mathrm{obs}}+(1-q)p_{\mathrm{mis}}. \end{aligned} $$

Because \(Y\) is binary, the only restriction available without further assumptions is \(0\leq p_{\mathrm{mis}}\leq1\). Multiplication by the nonnegative fraction \(1-q\) preserves these inequalities, giving

$$ q\,p_{\mathrm{obs}} \leq q\,p_{\mathrm{obs}}+(1-q)p_{\mathrm{mis}} \leq q\,p_{\mathrm{obs}}+(1-q). $$

This justification also explains sharpness. One can complete the unobserved outcomes by assigning \(Y=0\) to every censored person to attain the lower endpoint, or \(Y=1\) to every censored person to attain the upper endpoint, without changing any observed outcome. Additional restrictions exclude some such completions and can narrow the interval.

Diagnostics, Alternatives, and Limitations

Before interpreting a bound, check what counts as an unknown outcome. A competing event observed before censoring is a known zero for a cause-specific cumulative incidence; it must not be relabelled as censoring. Likewise, people known to be event-free through \(\tau\) have known zero outcomes. Counting either group as missing needlessly widens the bound. Conversely, treating a person censored before \(\tau\) as known event-free understates uncertainty.

The worst-case interval is especially informative when few people are censored before the horizon: its width is \(1-q\). When censoring is extensive, the interval may be too wide for a useful conclusion. Check the number and fraction with unknown status, the event coding, and whether the chosen horizon leaves enough follow-up to observe the outcome. If bounds are calculated within covariate groups, inspect whether each group has observed outcomes; a group with no known outcomes has no observed-group risk to support an ordering restriction.

These bounds answer a different question from approaches that produce a point estimate. Inverse probability of censoring weighting and related methods use assumptions such as conditional independent censoring, together with positivity, to estimate a risk. Selection models and pattern-mixture models make different assumptions about unobserved outcomes. Such assumptions can yield narrower estimates, but the observed data alone do not verify them. Bounds can complement those approaches by making the role of assumptions visible; they do not automatically validate a point estimate or provide a sampling confidence interval.

The basic bounds also have limitations. They use only the binary outcome at a chosen horizon, so they do not identify the full event-time distribution or the timing of events after censoring. They can be wide, and assumptions used to narrow them may be difficult to justify. In a multi-state setting, the same logic applies to a binary indicator such as being in a specified state at \(\tau\): known state membership contributes observed outcomes, while censored individuals have unknown membership. The target and what constitutes a known status must be defined for that estimand.

Common Mistake

  • Wrong reasoning: “Censoring was unrelated to the measured covariates, so the observed event proportion is the true risk.” This does not rule out dependence between censoring and unmeasured outcomes.
  • Correction: Independence from measured covariates is not, by itself, evidence that \(p_{\mathrm{mis}}=p_{\mathrm{obs}}\). Without a justified censoring assumption or a restriction on \(p_{\mathrm{mis}}\), report the worst-case bounds. If an assumption is used to narrow them, state it and show how results change under alternatives.
Key takeaway: With unrestricted outcomes among people censored before the horizon, the risk is bounded below by observed events per population member and above by that quantity plus the fraction censored. A narrower interval always reflects additional information or an additional assumption.

Check Your Understanding

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

  1. If 75% of a cohort has known status through \(\tau\), and 10% of those people have cause \(k\) by \(\tau\), what are the worst-case risk bounds?
  2. In a competing-risks analysis, should an observed competing event before \(\tau\) be counted as an unknown outcome for the cause-\(k\) risk? Why?
  3. If \(p_{\mathrm{obs}}=0.30\), \(q=0.80\), and the assumption is \(p_{\mathrm{mis}}\leq p_{\mathrm{obs}}\), calculate the resulting interval.
  4. What does the width of the unrestricted bound measure in this setup?

Solutions. (1) The lower bound is \(0.75(0.10)=0.075\); the upper bound is \(0.075+(1-0.75)=0.325\), or 7.5% to 32.5%. (2) No. An observed competing event establishes that cause \(k\) was not the first event by \(\tau\), so it is a known zero for this outcome. (3) The lower endpoint is \(0.80(0.30)=0.24\), and the upper endpoint is \(0.24+0.20(0.30)=0.30\). Thus the interval is 24% to 30%. (4) Its width is \(1-q\), the fraction whose status through \(\tau\) is unknown; each such person could, without further assumptions, have either value of the binary outcome.

Clinical Biostats