Tutorials › Exam ASTAM - Advanced Short-Term Actuarial Mathematics › Choosing a Restricted-Mean Horizon Under Censoring

Exam ASTAM · Kaplan-Meier and Nelson-Aalen · Tutorial 255 of 979

Choosing a Restricted-Mean Horizon Under Censoring

Choose and interpret a defensible restricted-mean horizon when follow-up is incomplete, and recognize where Kaplan-Meier estimates stop being supported by observations.

Intermediate 9 min read

What You'll Learn

  • Define the restricted mean survival time at a specified horizon.
  • Select a horizon that is supported by observed follow-up.
  • Choose a common horizon for comparing groups with different follow-up.
  • Interpret the Kaplan-Meier area without extending it beyond the data.
  • Identify why sparse follow-up near a horizon can make an estimate unstable.
  • Avoid selecting a horizon to favor a particular comparison.

Why the Horizon Matters

The restricted mean survival time summarizes survival over a specified time window. Its value depends on the window’s endpoint, so choosing that endpoint is part of the analysis—not a cosmetic detail. When follow-up is incomplete, the horizon must also be supported by the observed data. Otherwise, a calculation may look precise while relying on survival behavior that was never observed.

The Kaplan-Meier (KM) curve estimates survival using observed event times and the people still at risk at those times. Between observed events, the estimated curve is horizontal. That convention lets us calculate an area, but it does not mean we have observed what happens indefinitely after the last person leaves follow-up. In particular, a flat tail after the last follow-up time is not evidence that no later failures would occur.

Definition: For a chosen horizon \(\tau\), the restricted mean survival time (RMST) is the area under the survival function from time 0 through \(\tau\). Its Kaplan-Meier estimate is
$$ \widehat{\mu}(\tau)=\int_0^\tau \widehat{S}_{KM}(t)\,dt. $$
The result is expressed in the same time units as \(t\). It is the estimated mean survival time accumulated within the first \(\tau\) time units, not an estimate of unrestricted mean survival time.

The calculation uses the KM curve as a step function. To interpret that curve as an estimate of the event-time survival distribution, the usual non-informative censoring condition is required: censoring must not systematically hide different event risks among people who remain under observation and those who are censored. In the simplest formulation, event times and censoring times are independent. If censoring is informative, the ordinary KM estimator may not represent the target survival distribution.

Choosing a Supported Horizon

A defensible horizon is one for which there is follow-up information through the end of the interval. For one group, the largest observed follow-up time is a practical upper bound on a data-supported horizon. For a comparison of groups, a common horizon should be supported in every group; a conservative choice is no later than the shortest maximum follow-up among them. The horizon should also be chosen for the question being asked, preferably before examining which value produces the most favorable result.

Being within the observed follow-up range is necessary, but it does not guarantee a stable estimate. If only a few people remain under observation near \(\tau\), the estimated curve near the endpoint can be sensitive to individual events. Thus horizon selection should consider both the follow-up limit and the number at risk near that limit. A shorter horizon with broader support may be more interpretable than a longer one based on very sparse follow-up.

Key distinction: The KM curve is constant between observed event times because its estimate changes only at events. After the last observed follow-up time, a constant continuation is a mathematical convention, not an estimate supported by later observations. Do not use that continuation to claim an RMST through a later horizon.

In practice, state the selected horizon and its units, explain why it is supported, and report the RMST for that horizon. When comparing groups, use the same \(\tau\) for each group; otherwise, the difference combines survival differences with differences in the time window being summarized.

Why the RMST Is an Area

The area formula follows from the definition of a nonnegative event time \(T\). For each individual, the time contributed before the horizon is \(\min(T,\tau)\). At any time \(t\) between 0 and \(\tau\), this quantity includes a small interval of length \(dt\) exactly when the person survives beyond \(t\). Averaging those interval contributions gives the integral of the probability of surviving beyond \(t\).

$$ E[\min(T,\tau)] =\int_0^\tau P(T>t)\,dt =\int_0^\tau S(t)\,dt. $$

The first equality can be seen by writing \(\min(T,\tau)\) as the accumulated indicator of survival through time \(t\): it equals \(\int_0^\tau I(T>t)\,dt\), apart from endpoint conventions that do not affect the integral. Taking expectations gives the displayed result. Replacing \(S(t)\) with the estimated KM step function gives the estimated RMST.

Because the KM estimate is piecewise constant, each part of the integral is a rectangle. If the estimated survival is \(s_j\) over an interval of length \(L_j\), that interval contributes \(L_js_j\). For a horizon between event times, the final rectangle ends at \(\tau\), not at the next event time.

Worked Example: Choosing a Horizon Before the Follow-Up Limit

In a hypothetical study of a device-replacement program, eight units are followed for time to first failure. One unit fails at time 1, one is censored at time 2, two fail at time 3, and the remaining four are censored at time 4. Estimate the RMST through \(\tau=3.5\). Is a horizon of 5 supported?

Method and risk sets. The KM estimate changes at failure times. At time 1, eight units are at risk and one fails, so the survival multiplier is \(7/8\). At time 2, the censoring reduces the later risk set but does not cause a KM drop. At time 3, six units are at risk: the original eight less the failure at 1 and the censoring at 2. Two fail, giving a multiplier of \(4/6=2/3\). The estimated survival after time 3 is therefore \((7/8)(2/3)=7/12\).

$$ \widehat{S}_{KM}(t)= \begin{cases} 1, & 0\le t<1,\\ 7/8, & 1\le t<3,\\ 7/12, & 3\le t\le 4. \end{cases} $$

Area through 3.5. The three interval contributions are \(1(1)=1\), \(2(7/8)=7/4\), and \(0.5(7/12)=7/24\). Thus

$$ \widehat{\mu}(3.5) =1+\frac{7}{4}+\frac{7}{24} =\frac{24+42+7}{24} =\frac{73}{24} \approx 3.0417. $$

The estimate is \(3.0417\) time units, rounded to four decimal places. A horizon of 5 is not supported: the latest observed follow-up is time 4. Extending the horizontal KM segment from 4 to 5 would add \(1(7/12)=7/12\), but that area would describe unobserved follow-up. The common mistake is to treat a horizontal tail as evidence of continued survival; here it is only the estimator’s value through the observed limit.

Common Horizons for Group Comparisons

A group-specific RMST answers how much survival time, on average, is accumulated within a given window for that group. To compare groups fairly, calculate both areas through the same \(\tau\). If one group has follow-up only through time 4 and another through time 6, a comparison through time 6 is not supported for both groups. One option is a common horizon no later than time 4, subject to adequate follow-up near that endpoint.

The difference \(\widehat{\mu}_A(\tau)-\widehat{\mu}_B(\tau)\) is interpreted in time units. A negative value means the estimated accumulated survival time through \(\tau\) is smaller for group A. It does not, by itself, establish why the groups differ or describe survival after \(\tau\).

Worked Example: Comparing Two Groups at a Supported Horizon

In a hypothetical comparison of two maintenance schedules, Group A has KM survival 1 from time 0 to 1, 0.8 from 1 to 3, and 0.5 from 3 to 5. Its five remaining units are censored at time 5. Group B has survival 1 from 0 to 2 and 0.75 from 2 to 4; its six remaining units are censored at time 4. Estimate each RMST through the common horizon \(\tau=4\), then calculate A minus B.

Horizon check. The maximum observed follow-up is 5 in A and 4 in B, so time 4 is within the observed follow-up range of both groups. It is the largest common endpoint available from those limits; whether it is sufficiently stable should also be assessed using the numbers at risk near 4.

Group A. The survival levels over the interval are 1 for length 1, 0.8 for length 2, and 0.5 for length 1. Therefore,

$$ \widehat{\mu}_A(4) =1(1)+2(0.8)+1(0.5) =1+1.6+0.5 =3.1. $$

Group B. The survival levels are 1 for length 2 and 0.75 for length 2. Therefore,

$$ \widehat{\mu}_B(4) =2(1)+2(0.75) =2+1.5 =3.5. $$

Difference and interpretation. A minus B is \(3.1-3.5=-0.4\) time units. Through time 4, Group A has an estimated 0.4 fewer units of accumulated survival time than Group B. The result does not justify extending either group’s curve past time 4. In particular, a longer follow-up range for A does not make a later comparison fair if B has no observations beyond 4.

Why Extrapolation Is Not Supported

The KM estimator is built from the observed risk sets and event counts. In the usual product-limit form, each observed failure time contributes a conditional survival multiplier based on the number at risk just before that time. Once everyone has left observation, there is no later risk set and no later multiplier to estimate. Thus the calculation supplies no information about the probability of survival after the last observed follow-up.

A software plot may visually continue the final horizontal segment. That display is useful for showing the last estimated level, but it does not turn an unobserved time interval into observed evidence. A parametric survival model could produce predictions beyond follow-up by imposing a distributional form; that is a different method and relies on additional assumptions. The nonparametric KM curve alone does not provide that extrapolation.

Worked Example: Correct Area, Unsupported Extension

In a hypothetical component study, nine components are followed. Four fail at time \(13/4=3.25\), and the other five are censored at time 5. Calculate the RMST through time 5. Then assess a proposed calculation through time 8 that holds survival constant after time 5.

KM step. At time 3.25, all nine components are at risk and four fail. The survival estimate drops to \(1(5/9)=5/9\). The five remaining components are observed until they are censored at time 5, so the estimated survival level from 3.25 through 5 is \(5/9\).

RMST through 5. The curve is 1 for \(13/4\) time units and \(5/9\) for \(5-13/4=7/4\) time units. Therefore,

$$ \widehat{\mu}(5) =\frac{13}{4}+\frac{7}{4}\left(\frac{5}{9}\right) =\frac{117}{36}+\frac{35}{36} =\frac{152}{36} =\frac{38}{9} \approx 4.2222. $$

The result is \(4.2222\) time units, rounded to four decimal places. The proposed extension to time 8 would add \(3(5/9)=15/9\), giving \(38/9+15/9=53/9\). That arithmetic is correct, but the extended value is not a data-supported RMST estimate: no component was observed beyond time 5. The common mistake is to confuse a correct rectangle calculation with a justified horizon.

Practical Checks Before Reporting

Before reporting an RMST, check the horizon and the interpretation as carefully as the arithmetic. A defensible report states \(\tau\), its units, the estimated area, and how the horizon relates to observed follow-up. For a group comparison, it also states that the same horizon was used in both groups.

  • Check the follow-up limit. Do not integrate a KM curve beyond the largest observed follow-up time for the population being summarized.
  • Check support near the endpoint. A horizon near the end of follow-up may be based on very few people still under observation, even if it is technically within the observed range.
  • Use a common comparison window. Different horizons produce summaries over different periods and can make a group difference hard to interpret.
  • Choose for the question, not the preferred result. Selecting \(\tau\) after seeing which horizon makes one group look better can distort the comparison.
  • Keep the meaning restricted. The RMST is expected survival accumulated up to \(\tau\); it is not the unrestricted mean unless additional conditions support that interpretation.
Key takeaway: Choose a horizon that is both relevant to the question and supported by follow-up in every population being summarized. Integrate the KM curve only through that horizon; a flat continuation beyond observed follow-up is not evidence for extrapolation.

Check Your Understanding

For each problem, check both the area calculation and whether the proposed horizon is supported by the stated follow-up.

  1. A hypothetical cohort has estimated survival 1 from time 0 to 2 and 0.7 from time 2 to 5. Follow-up ends at time 5. Calculate the RMST through time 4. (Answer: 3.4 time units.)
  2. In a hypothetical study, one group has maximum follow-up 7 and a second group has maximum follow-up 5. What is the largest common horizon allowed by these limits, before considering numbers at risk? (Answer: 5 time units.)
  3. A hypothetical KM curve has survival 0.6 from time 3 through the last observed follow-up at time 6. A proposed RMST horizon is 8. What area would the unsupported extension add, and is it justified by the KM data? (Answer: It would add \(2(0.6)=1.2\) time units; it is not justified by the observed follow-up.)
  4. Through a common horizon of 4, Group C has estimated survival 1 from 0 to 1, 0.8 from 1 to 3, and 0.5 from 3 to 4. Group D has survival 1 from 0 to 2 and 0.75 from 2 to 4. Calculate the RMST difference, C minus D. (Answer: \(-0.4\) time units.)
  5. A hypothetical study has maximum follow-up 6, but only two participants remain under observation at time 5. Is a horizon of 5 automatically a stable choice simply because it is within the observed follow-up range? (Answer: No; it is within the range, but the sparse risk set near 5 may make the estimate unstable.)
Clinical Biostats