Tutorials › Biostatistics › Restricted Mean Survival Time (RMST)

Survival Analysis

Restricted Mean Survival Time (RMST)

A practical guide to Restricted Mean Survival Time, including its definition, interpretation as an average survival time up to a prespecified horizon, calculation from Kaplan-Meier curves, treatment-group comparisons, confidence intervals, hypothesis testing, clinical interpretation, and implementation in R.

Advanced 18 min read

What You'll Learn

  • What Restricted Mean Survival Time means and why it is useful
  • How RMST is calculated as the area under a survival curve
  • How to compare RMST between treatment groups
  • Why RMST can be useful when proportional hazards is questionable
  • How the truncation time changes the interpretation of RMST
  • How to perform RMST analyses in R and report the results clinically

Introduction

Survival analysis is often summarized using the Kaplan-Meier survival curve, the median survival time, or a hazard ratio from a Cox proportional hazards model.

These methods are extremely useful, but each answers a somewhat different question.

The Restricted Mean Survival Time (RMST) provides another way to summarize survival that is especially attractive when the proportional hazards assumption is questionable, the median survival is not reached, or the survival curves differ in a way that is not well represented by a single hazard ratio.

The basic idea is remarkably intuitive: RMST is the average amount of survival time experienced during a prespecified follow-up period.

Key idea: If the analysis is restricted to a time horizon \( \tau \), the RMST is the area under the survival curve from time 0 through \( \tau \). It therefore has a direct interpretation in units of time.

What Is Restricted Mean Survival Time?

Let \(S(t)\) denote the survival function:

\[ S(t)=P(T>t) \]

where \(T\) is the event time.

Choose a clinically meaningful restriction time \( \tau \). The Restricted Mean Survival Time is:

\[ \operatorname{RMST}(\tau) = \int_0^\tau S(t)\,dt \]

This is the expected survival time accumulated up to \( \tau \).

If time is measured in months, RMST is measured in months. If time is measured in years, RMST is measured in years.

Interpretation: An RMST of 18 months at \( \tau=24 \) means that the average participant is estimated to survive approximately 18 months during the first 24 months of follow-up.

Why Is RMST Called "Restricted"?

The ordinary mean survival time is:

\[ E(T) = \int_0^\infty S(t)\,dt \]

In many survival studies, however, the entire survival distribution cannot be observed.

For example, a substantial proportion of patients may remain alive at the end of follow-up.

The unrestricted mean survival time may therefore be undefined or impossible to estimate reliably.

RMST avoids this problem by restricting attention to a finite time horizon:

\[ E\{\min(T,\tau)\} = \int_0^\tau S(t)\,dt \]

This quantity is finite whenever \( \tau \) is finite.

The Geometric Interpretation

The most intuitive way to understand RMST is geometrically.

Imagine the Kaplan-Meier survival curve plotted against time. The RMST up to \( \tau \) is the area underneath the curve from 0 to \( \tau \).

1
Start at time 0, where the estimated survival probability is approximately 1.
2
Follow the Kaplan-Meier survival curve through the selected follow-up period.
3
Draw a vertical boundary at the prespecified restriction time \( \tau \).
4
Calculate the area under the survival curve up to \( \tau \).
5
Interpret that area as the average survival time accumulated during the restricted period.

RMST and the Kaplan-Meier Curve

The Kaplan-Meier estimator is a step function. Suppose the estimated survival curve remains at a particular value between two event times.

The area under the curve can therefore be calculated as a sum of rectangles.

If the survival estimate is \(S(t_i)\) between consecutive event times, the contribution to RMST is approximately:

\[ S(t_i)(t_{i+1}-t_i) \]

Summing these contributions gives the estimated RMST.

Practical interpretation: RMST does not require the survival curve to have a particular parametric form. It can be calculated directly from the estimated survival curve.

Why RMST Is Useful in Clinical Trials

RMST is particularly useful when the treatment effect is not adequately summarized by a proportional hazards model.

For example, survival curves may:

  • Separate gradually over time.
  • Show delayed treatment effects.
  • Cross one another.
  • Have treatment effects that change over time.
  • Show early harm followed by later benefit.
  • Show long-term plateaus.

In such situations, a single hazard ratio may provide an incomplete description of the treatment effect.

RMST Versus the Median Survival Time

The median survival time is the time at which the survival probability reaches 0.50.

\[ S(t_{0.5})=0.5 \]

The RMST asks a different question.

\[ \operatorname{RMST}(\tau) = \int_0^\tau S(t)\,dt \]

The median identifies one point on the survival distribution. RMST summarizes the entire survival curve up to \( \tau \).

Measure Question Answered
Median survival When has approximately 50% of the population experienced the event?
RMST How much survival time is accumulated on average up to \( \tau \)?
Hazard ratio How do the instantaneous event rates compare under a proportional hazards model?

When the Median Survival Is Not Reached

Suppose a study follows patients for 24 months, but more than half of the patients in the treatment group remain alive.

The median survival time may then be reported as:

Not reached.

That statement is informative, but it does not provide a numerical summary of how much survival time was observed during the study period.

RMST can still be calculated if a suitable restriction time is chosen.

Example: If \( \tau=24 \) months and the treatment-group RMST is 20.1 months, the study can report an estimated 20.1 months of average survival during the first 24 months even if the median survival has not been reached.

Choosing the Restriction Time \( \tau \)

The choice of \( \tau \) is one of the most important decisions in an RMST analysis.

The restriction time should generally be chosen based on clinical and study considerations rather than selected after examining the treatment effect.

Common considerations include:

  • Planned study follow-up.
  • A clinically meaningful time horizon.
  • The point by which adequate follow-up is expected in both groups.
  • Trial-specific estimand considerations.
  • Availability of reliable risk-set information.
Important: The RMST depends on \( \tau \). Therefore, an RMST analysis should always report the restriction time explicitly.

A Simple Example

Suppose a clinical trial follows patients for 24 months. Consider a hypothetical treatment group whose survival curve has the following simplified structure:

Interval Estimated Survival
0 to 6 months 1.00
6 to 12 months 0.90
12 to 18 months 0.75
18 to 24 months 0.60

The RMST is the sum of the areas of the four rectangles:

\[ \operatorname{RMST}(24) = (6)(1.00) + (6)(0.90) + (6)(0.75) + (6)(0.60) \]

Therefore:

\[ \operatorname{RMST}(24) = 6+5.4+4.5+3.6 = 19.5 \]

The estimated RMST is therefore:

\[ \boxed{\operatorname{RMST}(24)=19.5\text{ months}} \]

The interpretation is: the average participant experiences approximately 19.5 months of survival during the first 24 months.

Comparing Two Treatment Groups

Suppose a randomized trial compares treatment A with treatment B. Let:

\[ \mu_A(\tau) = \operatorname{RMST}_A(\tau) \]

and:

\[ \mu_B(\tau) = \operatorname{RMST}_B(\tau) \]

A natural treatment effect measure is the RMST difference:

\[ \Delta_{\mathrm{RMST}}(\tau) = \mu_A(\tau)-\mu_B(\tau) \]

If treatment A has an RMST of 19.5 months and treatment B has an RMST of 16.8 months at \( \tau=24 \) months:

\[ \Delta_{\mathrm{RMST}}(24) = 19.5-16.8 = 2.7\text{ months} \]
Clinical interpretation: At 24 months, treatment A provides an estimated 2.7 additional months of average survival compared with treatment B, within the restricted 24-month period.

The RMST Difference Is Often Easy to Explain

One advantage of RMST is that its treatment effect is expressed directly in time units.

Compare:

Measure Example Interpretation
Hazard ratio 0.72 Relative instantaneous event-rate comparison under the model
RMST difference +2.7 months Average additional survival time through the selected horizon

For many clinical audiences, "2.7 additional months of average survival through 24 months" can be more immediately interpretable than a hazard ratio.

RMST Ratio

A second possible treatment effect is the RMST ratio:

\[ \operatorname{RMST\ Ratio} = \frac{\mu_A(\tau)}{\mu_B(\tau)} \]

If:

\[ \mu_A(24)=19.5 \qquad \mu_B(24)=16.8 \]

then:

\[ \frac{19.5}{16.8} \approx1.161 \]

Thus the RMST under treatment A is approximately 16.1% greater than that under treatment B over the first 24 months.

Difference versus ratio: The RMST difference is expressed in units of time and is often easier to interpret clinically. The RMST ratio provides a relative measure and may be useful when a proportional comparison is desirable.

RMST and the Area Between Two Survival Curves

The difference in RMSTs has another useful geometric interpretation.

If treatment A has survival curve \(S_A(t)\) and treatment B has survival curve \(S_B(t)\), then:

\[ \Delta_{\mathrm{RMST}}(\tau) = \int_0^\tau \left[ S_A(t)-S_B(t) \right]dt \]

Therefore, the RMST difference is the signed area between the two survival curves up to \( \tau \).

This interpretation becomes particularly important when the curves cross.

What Happens When Survival Curves Cross?

Suppose treatment A has better survival early in follow-up but treatment B has better survival later.

The survival curves may cross.

A single hazard ratio can be difficult to interpret in this setting because the relative treatment effect changes over time.

RMST remains well defined.

\[ \Delta_{\mathrm{RMST}}(\tau) = \int_0^\tau [S_A(t)-S_B(t)]dt \]

Positive and negative regions of the difference in survival probabilities are automatically incorporated into the integrated treatment effect.

Important: A statistically significant RMST difference does not imply that one treatment has a constant advantage at every point in time. It describes the integrated difference in survival through the selected restriction time.

RMST and the Proportional Hazards Assumption

The Cox proportional hazards model assumes that the hazard ratio is constant over time:

\[ \frac{h_A(t)}{h_B(t)} = \theta \]

If proportional hazards is approximately true, the hazard ratio can be a useful summary.

If proportional hazards is not plausible, however, a single hazard ratio may not adequately describe the treatment effect.

RMST does not require a constant hazard ratio.

Practical principle: RMST should not be thought of as a replacement for the Cox model in every analysis. Rather, it is an alternative estimand that can provide a more direct and robust summary when the treatment effect varies over time.

RMST Does Not Assume Proportional Hazards

The definition of RMST is:

\[ \operatorname{RMST}(\tau) = \int_0^\tau S(t)\,dt \]

No proportional hazards assumption appears in this definition.

The analysis is based directly on the survival functions.

This makes RMST particularly attractive for:

  • Delayed treatment effects.
  • Crossing survival curves.
  • Non-proportional hazards.
  • Immuno-oncology trials with delayed separation.
  • Studies with long-term survival plateaus.
  • Other settings in which a single hazard ratio is difficult to interpret.

RMST and Delayed Treatment Effects

Consider an immunotherapy trial in which the treatment effect is minimal during the first several months but becomes increasingly apparent later.

The proportional hazards assumption may be questionable.

An RMST analysis can quantify the accumulated survival benefit over a specified horizon without requiring the treatment effect to be constant.

For example:

\[ \Delta_{\mathrm{RMST}}(24)=3.1\text{ months} \]

can be interpreted directly as an average 3.1-month survival advantage through 24 months.

RMST at Different Restriction Times

Because RMST depends on \( \tau \), different restriction times can produce different treatment effects.

For example:

Restriction Time RMST Difference
12 months +0.4 months
18 months +1.6 months
24 months +3.1 months
36 months +4.2 months

This tells a much richer story than one hazard ratio.

The treatment effect appears to accumulate over time.

Choosing a Clinically Meaningful Horizon

Suppose a trial is expected to have meaningful follow-up for approximately 36 months. An RMST analysis might therefore specify:

\[ \tau=36\text{ months} \]

The resulting estimand is:

\[ \operatorname{RMST}(36) = \int_0^{36}S(t)\,dt \]

The treatment comparison becomes:

\[ \Delta_{\mathrm{RMST}}(36) = \operatorname{RMST}_A(36) - \operatorname{RMST}_B(36) \]

This is the estimated difference in average survival time accumulated during the first 36 months.

Why \( \tau \) Should Usually Be Prespecified

If investigators examine many possible restriction times and select the one that produces the most favorable treatment result, the nominal statistical properties of the analysis can be compromised.

The restriction time should therefore generally be defined before the confirmatory analysis.

Good practice: Document the rationale for \( \tau \) in the protocol or statistical analysis plan.

RMST and Censoring

Clinical survival data commonly contain right-censored observations. For example, a patient may leave the study before experiencing the event.

The Kaplan-Meier estimator accounts for right censoring under the usual independent-censoring framework.

RMST can then be estimated by integrating the Kaplan-Meier survival estimate.

Thus, censoring does not mean that patients must simply be excluded from the RMST analysis.

Key point: RMST is an estimand based on the survival distribution, while the estimated survival distribution is obtained using methods that account for censoring.

Administrative Censoring

Suppose the study ends at a planned database lock. Patients who have not experienced the event by that time may be administratively censored.

This is compatible with an RMST analysis provided the analysis appropriately handles the resulting right-censored observations.

The restriction time should also be chosen so that the survival curve can be estimated with adequate information through the selected horizon.

Why Very Late \( \tau \) Can Be Problematic

Suppose almost all patients have been censored by 60 months. The estimated survival curve beyond the last few events may then be based on a very small risk set.

An RMST calculated at an excessively late \( \tau \) can therefore become unstable or heavily dependent on sparse late follow-up.

This is one reason that \( \tau \) should be selected using clinical and statistical considerations.

RMST Compared With Restricted Mean Time Lost

Another useful quantity is the restricted mean time lost (RMTL).

Over the interval from 0 to \( \tau \):

\[ \operatorname{RMTL}(\tau) = \tau-\operatorname{RMST}(\tau) \]

RMTL can be interpreted as the average amount of time lost to the event during the restricted period.

For example, if:

\[ \tau=24 \]

and:

\[ \operatorname{RMST}(24)=18 \]

then:

\[ \operatorname{RMTL}(24) = 24-18 = 6\text{ months} \]

The average participant loses approximately six months of potential survival during the first 24 months.

RMTL and RMST Differences

Because:

\[ \operatorname{RMTL}(\tau) = \tau-\operatorname{RMST}(\tau) \]

the difference in RMTL between two groups is the negative of the RMST difference:

\[ \operatorname{RMTL}_A(\tau) - \operatorname{RMTL}_B(\tau) = -\Delta_{\mathrm{RMST}}(\tau) \]

Therefore, a positive RMST difference corresponds to less restricted mean time lost.

RMST as an Alternative Estimand

An important conceptual point is that RMST is not merely a different way of calculating a hazard ratio.

It defines a different treatment-effect estimand.

The hazard ratio compares hazards.

The RMST difference compares average survival time accumulated over a specified period.

Estimand Primary Quantity
Hazard ratio Relative hazard
Median survival difference Difference between survival medians
RMST difference Difference in average survival through \( \tau \)
RMST ratio Ratio of average survival through \( \tau \)
RMTL difference Difference in average time lost through \( \tau \)

Hypothesis Testing With RMST

A common hypothesis is:

\[ H_0: \Delta_{\mathrm{RMST}}(\tau)=0 \]

against:

\[ H_A: \Delta_{\mathrm{RMST}}(\tau)\ne0 \]

A one-sided alternative may also be used when the direction of benefit is prespecified.

The analysis estimates the RMST in each treatment group and the difference between them, together with an appropriate standard error and confidence interval.

Confidence Intervals

Suppose the estimated treatment effect is:

\[ \widehat{\Delta}_{\mathrm{RMST}}=2.7 \]

with a 95% confidence interval of:

\[ (0.8,\;4.6) \]

The clinical interpretation is: the estimated average survival advantage through the specified restriction time is 2.7 months, with a 95% confidence interval from 0.8 to 4.6 months.

Because the confidence interval excludes zero, the corresponding two-sided hypothesis test would reject the null hypothesis at the 5% level under the usual inferential framework.

RMST Standard Errors

Several methods can be used to estimate the variance of the RMST estimator. One common approach is based on the variance of the area under the Kaplan-Meier curve.

The details can become mathematically involved, but the practical objective is straightforward:

  • Estimate RMST.
  • Estimate its standard error.
  • Construct a confidence interval.
  • Compare treatment groups.

Software packages can perform these calculations without requiring investigators to manually derive the variance estimator.

Worked Two-Group Example

Suppose a randomized oncology trial evaluates overall survival through 24 months. The estimated RMSTs are:

Group RMST Through 24 Months
Experimental treatment 19.4 months
Control 16.7 months

The RMST difference is:

\[ \widehat{\Delta}_{\mathrm{RMST}} = 19.4-16.7 = 2.7\text{ months} \]

Suppose the 95% confidence interval is:

\[ (0.9,\;4.5) \]

The result can be reported as:

Clinical interpretation: Through 24 months, the experimental treatment was associated with an estimated 2.7 additional months of average survival compared with control (95% CI: 0.9 to 4.5 months).

What If the Median Is Not Reached?

Suppose median overall survival is:

Group Median OS RMST at 24 Months
Experimental Not reached 19.4 months
Control 18.2 months 16.7 months

The RMST provides a numerical treatment comparison even though the experimental median is not estimable from the observed follow-up.

Important distinction: "Median not reached" does not mean that survival is infinite or that no treatment comparison is possible. It means that the estimated survival curve did not fall to 0.50 during the available follow-up.

RMST and Hazard Ratios Can Tell Different Stories

Suppose a trial reports:

\[ HR=0.78 \]

and:

\[ \Delta_{\mathrm{RMST}}(24)=2.9\text{ months} \]

These are not competing estimates of the same quantity. The hazard ratio describes a relative hazard under the Cox model, while the RMST difference describes accumulated survival time through 24 months.

Both can be reported if both estimands are clinically and statistically relevant.

When the Hazard Ratio Can Be Difficult to Interpret

Consider survival curves with delayed separation. During the first several months, the hazards may be similar. Later, the treatment group may experience substantially fewer events.

A single hazard ratio averages information in a way that may obscure the timing of the treatment effect.

RMST provides a direct time-based summary of the resulting survival curves.

RMST in Immuno-Oncology

RMST is particularly relevant to immuno-oncology because delayed treatment effects and non-proportional hazards can occur.

For example, suppose:

  • The treatment effect is modest during the first six months.
  • Survival curves begin separating after six months.
  • A long-term survival plateau develops in the treatment group.

A proportional hazards model may not provide a complete description of this pattern.

An RMST analysis at 24 or 36 months can summarize the accumulated survival advantage.

Clinical interpretation: RMST is particularly attractive when the question is "How much additional survival time is gained over a clinically meaningful period?" rather than simply "What is the hazard ratio?"

RMST and Non-Proportional Hazards

Non-proportional hazards occur when:

\[ HR(t) \]

changes meaningfully with time.

Examples include:

  • Delayed treatment effects.
  • Early treatment toxicity followed by later benefit.
  • Crossing hazards.
  • Mixture-cure-type behavior.
  • Long-term treatment plateaus.

RMST remains defined regardless of whether the hazard ratio is constant.

RMST Does Not Eliminate the Need to Inspect Survival Curves

Although RMST is useful, investigators should not analyze RMST in isolation.

A complete survival analysis should generally consider:

  • Kaplan-Meier curves.
  • Numbers at risk.
  • Event counts.
  • Censoring patterns.
  • Median survival when estimable.
  • RMST and its confidence interval.
  • Other clinically relevant estimands.

The graphical survival curves provide important context for interpreting the RMST.

Numbers at Risk Matter

Late survival estimates can become unstable when few patients remain under observation.

For this reason, an RMST analysis should not simply report a large value of \( \tau \) without examining the amount of information available near the end of the restricted interval.

Practical rule: Always examine the Kaplan-Meier curve and risk table before selecting or interpreting a late restriction time.

RMST With a Treatment Effect That Changes Sign

Suppose treatment A is worse than treatment B early but better later. The survival curves may cross.

The RMST difference integrates both periods:

\[ \Delta_{\mathrm{RMST}}(\tau) = \int_0^\tau [S_A(t)-S_B(t)]dt \]

The early disadvantage contributes negatively to the integral, while the later advantage contributes positively.

The final RMST difference therefore represents the net accumulated survival difference over the selected period.

Interpret carefully: If curves cross, the RMST difference may hide important differences in when benefit and harm occur. The underlying survival curves should therefore always be examined.

RMST Is Not a Probability

One common interpretation error is to treat RMST as if it were a survival probability.

It is not.

If:

\[ \operatorname{RMST}(24)=18.2 \]

the result means 18.2 months of average restricted survival time. It does not mean an 18.2% probability of surviving.

RMST Is Not the Mean of Observed Survival Times

Another common mistake is to calculate the ordinary arithmetic mean of observed follow-up times.

That approach is generally inappropriate in censored survival data because censored patients do not have fully observed event times.

RMST instead uses the estimated survival distribution while accounting for right censoring.

RMST and the Survival Probability at \( \tau \)

The survival probability at the restriction time is:

\[ S(\tau)=P(T>\tau) \]

RMST is:

\[ \operatorname{RMST}(\tau) = \int_0^\tau S(t)\,dt \]

These quantities answer different questions.

Measure Question
\(S(\tau)\) What proportion is expected to remain event-free beyond \( \tau \)?
RMST(\(\tau\)) How much event-free time is accumulated on average through \( \tau \)?

RMST and Restricted Mean Time Lost

The two quantities are directly related:

\[ \operatorname{RMST}(\tau) + \operatorname{RMTL}(\tau) = \tau \]

This identity can be useful when communicating results to clinical audiences.

If patients spend more time event-free, they necessarily spend less restricted time having experienced the event.

R Implementation

Several R packages can be used for RMST analysis. A commonly used approach is the survRM2 package.

The analysis begins with survival data containing:

  • A follow-up time.
  • An event indicator.
  • A treatment-group indicator.

For example:

library(survival)
library(survRM2)

time <- c(
  5, 8, 11, 14, 17,
  19, 22, 24, 24, 24,
  7, 10, 13, 16, 18,
  21, 23, 24, 24, 24
)

status <- c(
  1, 1, 0, 1, 0,
  1, 0, 1, 0, 0,
  1, 0, 1, 0, 1,
  0, 1, 0, 0, 0
)

arm <- c(
  1, 1, 1, 1, 1,
  1, 1, 1, 1, 1,
  0, 0, 0, 0, 0,
  0, 0, 0, 0, 0
)

Here:

  • time is the observed follow-up time.
  • status equals 1 for an event and 0 for censoring.
  • arm identifies the treatment group.

Choosing the Restriction Time in R

Suppose the clinically relevant horizon is 24 months:

tau <- 24

The RMST analysis should then be explicitly interpreted as an analysis through 24 months.

Kaplan-Meier Estimation in R

The survival curves can first be estimated using the survival package:

fit <- survfit(
  Surv(time, status) ~ arm
)

summary(fit)

The corresponding Kaplan-Meier curves can be plotted:

plot(
  fit,
  xlab = "Time",
  ylab = "Survival probability",
  lwd = 2
)

The graph should be examined before interpreting the RMST.

RMST Analysis With survRM2

A common implementation is:

result <- rmst2(
  time = time,
  status = status,
  arm = arm,
  tau = 24
)

result

The output provides estimates of the restricted mean survival times and comparisons between the two groups.

Interpreting RMST Output

Suppose the software produces:

Measure Estimate
RMST, treatment 19.2 months
RMST, control 16.5 months
RMST difference 2.7 months
95% CI 0.5 to 4.9 months
P-value 0.016

The appropriate interpretation is:

Example interpretation: Through 24 months, the estimated restricted mean survival time was 19.2 months in the treatment group and 16.5 months in the control group, corresponding to an RMST difference of 2.7 months (95% CI: 0.5 to 4.9; P=0.016).

R Code for a Reproducible RMST Analysis

library(survival)
library(survRM2)

# Survival object
surv_obj <- Surv(
  time = time,
  event = status
)

# Kaplan-Meier curves
km_fit <- survfit(
  surv_obj ~ arm
)

# RMST comparison
rmst_result <- rmst2(
  time = time,
  status = status,
  arm = arm,
  tau = 24
)

print(rmst_result)

RMST Using a Kaplan-Meier Step Function

The RMST can also be understood computationally by summing the areas under the Kaplan-Meier step function.

Suppose the survival curve has event times:

\[ t_1,t_2,\ldots,t_k \]

with corresponding survival estimates:

\[ \widehat S(t_1), \widehat S(t_2), \ldots, \widehat S(t_k) \]

The estimated RMST is the accumulated area up to \( \tau \).

Conceptually:

\[ \widehat{\operatorname{RMST}}(\tau) = \sum_i \widehat S(t_i) \left[ \min(t_{i+1},\tau)-t_i \right] \]

where only intervals within the restricted follow-up period contribute.

A Small Computational Example

Suppose the estimated survival function is:

Interval Survival
0–5 1.00
5–10 0.90
10–15 0.80
15–20 0.65

For \( \tau=20 \), the area is:

\[ (5)(1.00) + (5)(0.90) + (5)(0.80) + (5)(0.65) \]

Therefore:

\[ 5+4.5+4+3.25 = 16.75 \]

Thus:

\[ \boxed{\operatorname{RMST}(20)=16.75\text{ time units}} \]

RMST in a Clinical Trial Statistical Analysis Plan

If RMST is a prespecified primary or key secondary estimand, the statistical analysis plan should clearly define the analysis.

At minimum, specify:

  • The survival endpoint.
  • The event definition.
  • The censoring rules.
  • The analysis population.
  • The restriction time \( \tau \).
  • The RMST estimand.
  • The treatment contrast.
  • The confidence interval method.
  • The hypothesis test.
  • Any multiplicity adjustment.
  • Handling of intercurrent events.
  • Sensitivity analyses.

RMST as a Primary Endpoint

RMST can sometimes be used as a primary estimand rather than merely as a supplementary analysis.

If the primary question is: "How much additional survival time does treatment provide during the first 24 months?" then:

\[ \Delta_{\mathrm{RMST}}(24) \]

may be a natural primary treatment-effect measure.

This is especially attractive when a constant hazard ratio is not scientifically plausible.

RMST as a Secondary or Sensitivity Analysis

In many trials, the Cox hazard ratio remains the primary analysis while RMST is included as a complementary analysis.

For example:

Analysis Purpose
Kaplan-Meier curve Visual description of survival
Cox hazard ratio Relative hazard comparison
RMST difference Absolute survival-time comparison through \( \tau \)
Median survival Time to 50% survival when estimable

This combination can provide a much more complete picture of treatment benefit.

RMST and Estimands

Modern clinical trial analysis emphasizes explicitly defining the treatment effect of interest.

An RMST estimand should therefore be described precisely.

For example:

Example estimand: The difference in restricted mean overall survival through 24 months between patients randomized to experimental treatment and control, using the treatment policy strategy for intercurrent events and accounting for right censoring according to the prespecified survival-analysis methodology.

The exact wording should reflect the study's endpoint, intercurrent-event strategy, population, and analysis framework.

RMST and Intercurrent Events

In a clinical trial, events such as treatment discontinuation, switching, rescue therapy, or subsequent therapy can complicate interpretation.

The RMST itself does not automatically determine how these events should be handled.

The estimand must specify the strategy for such events.

Important: RMST is a statistical summary of a survival function. It does not by itself resolve the broader estimand question of what should happen after an intercurrent event.

RMST and Crossover

Suppose control-group patients are allowed to switch to the experimental treatment after progression.

The observed survival curves may then reflect both randomized treatment assignment and subsequent treatment exposure.

An RMST analysis can still be performed, but the interpretation depends on the estimand and the handling of treatment switching.

RMST and Competing Risks

Standard overall survival analysis treats death from any cause as the event. For cause-specific outcomes, competing risks may require a different analysis framework.

For example, if the endpoint is cancer-specific death while non-cancer death prevents cancer death, simply applying an ordinary survival analysis may answer a different question from the desired cumulative-incidence estimand.

RMST concepts can be extended to more complex event settings, but investigators must carefully define the underlying survival or state-occupation function.

Do not automatically apply ordinary RMST code to every competing-risk question. First define the estimand and the appropriate event process.

Common RMST Mistakes

  1. Failing to report \( \tau \). An RMST value without its restriction time is incomplete.
  2. Choosing \( \tau \) after seeing the treatment results. This can introduce selection bias and invalidate prespecified inferential properties.
  3. Calling RMST a survival probability. RMST is measured in units of time.
  4. Calculating the ordinary mean of observed survival times. That ignores the structure of right-censored observations.
  5. Ignoring late sparse risk sets. Very late restriction times may depend on limited information.
  6. Reporting RMST without showing the Kaplan-Meier curves. The curves provide important context for interpreting the integrated effect.
  7. Assuming RMST and hazard ratio answer the same question. They are different estimands.
  8. Ignoring crossing survival curves. The net RMST difference can conceal the timing of benefit and harm.
  9. Reporting only a P-value. The RMST estimate and confidence interval are essential for clinical interpretation.
  10. Using an inappropriate restriction time. The horizon should be clinically meaningful and statistically supportable.

RMST Versus Cox Regression

Feature Cox Model RMST
Primary effect measure Hazard ratio Difference or ratio in restricted mean survival
Units Unitless ratio Time for RMST difference
Requires proportional hazards? Standard interpretation does No
Depends on \( \tau \)? No Yes
Easy clinical interpretation Moderate Often high
Handles crossing curves? Potentially difficult to summarize Yes, through integration
Summarizes entire restricted survival curve? Not directly Yes

RMST Versus Log-Rank Testing

The log-rank test evaluates whether survival distributions differ between groups under its corresponding assumptions and weighting structure.

RMST instead produces an effect estimate in time units.

Thus, an RMST analysis provides more than a binary conclusion about whether a difference was detected.

Good reporting: Do not report only "RMST was statistically significant." Report the estimated RMST in each group, the RMST difference, its confidence interval, the restriction time, and the P-value.

RMST and Sample Size

Sample-size planning for RMST studies differs from simple exponential-survival or proportional-hazards calculations.

The required sample size depends on factors such as:

  • The expected survival distributions.
  • The restriction time.
  • The anticipated RMST difference.
  • Event rates.
  • Censoring.
  • Allocation ratio.
  • Type I error.
  • Desired power.

Simulation is often useful when the expected treatment effect is strongly non-proportional over time.

RMST and Power Under Non-Proportional Hazards

Suppose the treatment effect is delayed. A sample-size calculation based only on a constant hazard ratio may be poorly matched to the scientific hypothesis.

An RMST-based design can instead target a clinically meaningful difference such as:

\[ \Delta_{\mathrm{RMST}}(24)=2.5\text{ months} \]

The trial can then be powered to detect that time-based difference.

Simulation for RMST Designs

A general simulation workflow is:

1
Specify plausible control and treatment survival distributions.
2
Specify the restriction time \( \tau \).
3
Generate survival times and censoring times.
4
Calculate RMST differences for simulated trials.
5
Estimate the proportion of trials rejecting the null hypothesis.
6
Adjust the sample size until the desired power is achieved.

A Simple RMST Simulation Concept

set.seed(123)

n_per_arm <- 150

# Example exponential survival assumptions
lambda_control <- 0.060
lambda_treatment <- 0.045

time_control <- rexp(
  n_per_arm,
  rate = lambda_control
)

time_treatment <- rexp(
  n_per_arm,
  rate = lambda_treatment
)

# Add administrative censoring
censor_time <- 24

obs_control <- pmin(
  time_control,
  censor_time
)

obs_treatment <- pmin(
  time_treatment,
  censor_time
)

status_control <- as.integer(
  time_control <= censor_time
)

status_treatment <- as.integer(
  time_treatment <= censor_time
)

The simulated dataset can then be analyzed using the same RMST methodology that would be applied to the actual clinical trial.

Reporting RMST in a Clinical Study

A clear table might look like:

Measure Treatment Control Difference
RMST through 24 months 19.4 16.7 +2.7
95% CI for difference — — 0.9 to 4.5
P-value — — 0.004

The corresponding text might state:

At 24 months, the restricted mean overall survival was 19.4 months in the treatment group and 16.7 months in the control group. The estimated RMST difference was 2.7 months (95% CI: 0.9 to 4.5; P=0.004), favoring treatment.

Reporting Multiple Restriction Times

If multiple restriction times are clinically useful, investigators may report a prespecified set of analyses.

Horizon RMST Treatment RMST Control Difference
12 months 10.8 10.4 +0.4
24 months 19.4 16.7 +2.7
36 months 27.2 22.9 +4.3

This can demonstrate how the accumulated treatment effect changes over time.

However, multiple horizons should be prespecified and interpreted carefully, particularly if formal hypothesis testing is performed at multiple horizons.

RMST and Long-Term Survivors

In diseases with a possible long-term survivor population, RMST can capture differences in the survival curves during a fixed clinically meaningful period.

For example, if the treatment curve develops a plateau after 24 months, an RMST through 36 months incorporates the survival advantage accumulated during that interval.

This does not prove the existence of a cure fraction by itself.

Do not overinterpret a plateau: A long-term Kaplan-Meier plateau can be suggestive of durable survival, but specialized mixture or cure-model analyses may be required to formally model a cure fraction.

RMST in Oncology

RMST is particularly relevant to oncology trials because survival curves often exhibit patterns that make a single proportional-hazards summary difficult to interpret.

Potential applications include:

  • Overall survival.
  • Progression-free survival.
  • Disease-free survival.
  • Event-free survival.
  • Time-to-treatment-failure endpoints.

The precise interpretation depends on the endpoint definition.

RMST for Progression-Free Survival

Suppose progression-free survival is the endpoint. The RMST through 24 months is:

\[ \operatorname{RMST}_{PFS}(24) = \int_0^{24}S_{PFS}(t)\,dt \]

The treatment difference can then be interpreted as the average additional progression-free time accumulated during the first 24 months.

RMST for Disease-Free Survival

For disease-free survival:

\[ \operatorname{RMST}_{DFS}(\tau) = \int_0^\tau S_{DFS}(t)\,dt \]

The same mathematical framework applies. The clinical interpretation changes according to what constitutes an event.

RMST and Time-to-Event Endpoints Beyond Oncology

RMST is not limited to oncology. It can be used for many time-to-event outcomes, including:

  • Time to hospitalization.
  • Time to cardiovascular event.
  • Time to disease recurrence.
  • Time to treatment discontinuation.
  • Time to clinical deterioration.
  • Time to composite event.

The key requirement is that the endpoint and censoring process are appropriately defined.

RMST and Composite Endpoints

If the event is a composite, the RMST summarizes time until the first component of the composite endpoint.

For example, if a composite endpoint is:

  • Death,
  • Hospitalization, or
  • Major cardiovascular event,

then the RMST represents the average event-free time up to \( \tau \) for that composite endpoint.

Advantages of RMST

  • Direct interpretation in clinically meaningful time units.
  • Does not require proportional hazards for its definition.
  • Can be estimated even when the median survival is not reached.
  • Uses information across the survival curve.
  • Can accommodate delayed treatment effects.
  • Provides an absolute treatment-effect measure.
  • Can be interpreted geometrically as area under the survival curve.

Limitations of RMST

  • The result depends on the selected restriction time.
  • Late restriction times may be unstable when few patients remain at risk.
  • It does not replace examination of the survival curves.
  • It is a different estimand from the hazard ratio and should not be treated as a direct substitute in every setting.
  • Sample-size planning can require assumptions about the full survival curves.
  • Competing risks and intercurrent events require careful estimand definition.
  • Multiple restriction times can complicate formal inference.

A Practical RMST Workflow

1
Define the time-to-event endpoint and event definition.
2
Define the analysis population and censoring rules.
3
Choose a clinically meaningful restriction time \( \tau \).
4
Estimate Kaplan-Meier survival curves.
5
Calculate RMST within each treatment group.
6
Calculate the RMST difference or ratio.
7
Calculate the confidence interval and P-value.
8
Inspect the survival curves and numbers at risk.
9
Perform sensitivity analyses when appropriate.
10
Interpret the result in clinically meaningful time units.

What Should Be Reported?

A complete RMST result should identify:

  • The restriction time \( \tau \).
  • RMST in each treatment group.
  • The treatment contrast.
  • The confidence interval.
  • The P-value when applicable.
  • The survival endpoint.
  • The analysis population.
  • The method used to estimate RMST.

Example Statistical Analysis Plan Language

A statistical analysis plan might state:

Example: The restricted mean survival time through 24 months will be estimated separately for each randomized treatment group using the Kaplan-Meier survival estimator. The primary treatment contrast will be the difference in RMST between the experimental and control groups. The difference will be reported with a two-sided 95% confidence interval and corresponding P-value. The restriction time of 24 months was selected a priori based on the planned follow-up period and clinical relevance.

RMST and Sensitivity Analyses

Sensitivity analyses may examine whether conclusions are robust to reasonable changes in assumptions.

Potential analyses include:

  • A different clinically justified restriction time.
  • Alternative censoring rules.
  • Alternative analysis populations.
  • Adjustment for important baseline covariates.
  • Alternative methods for handling intercurrent events.

Sensitivity analyses should be distinguished clearly from the prespecified primary analysis.

Covariate Adjustment

RMST comparisons can also be adjusted for baseline covariates.

For example, investigators may wish to account for:

  • Age.
  • Disease stage.
  • Baseline performance status.
  • Biomarker status.
  • Other important prognostic factors.

The exact adjustment methodology should be prespecified and aligned with the study estimand.

RMST Does Not Require a Parametric Survival Model

One of the attractions of Kaplan-Meier-based RMST is that it is fundamentally nonparametric with respect to the survival distribution.

The analysis does not require investigators to assume:

  • Exponential survival.
  • Weibull survival.
  • Log-normal survival.
  • Constant hazards.

This can be advantageous when the shape of the survival distribution is uncertain.

Parametric RMST

RMST can also be calculated from parametric survival models.

For a fitted survival function \( \widehat S(t) \):

\[ \widehat{\operatorname{RMST}}(\tau) = \int_0^\tau \widehat S(t)\,dt \]

A parametric model may be useful when smooth survival estimates or extrapolation are required.

Important distinction: A Kaplan-Meier RMST is based on the observed follow-up distribution, whereas a model-based RMST can depend on the assumptions of the fitted survival model.

RMST and Extrapolation

Standard Kaplan-Meier RMST is most naturally interpreted within the observed follow-up period.

If investigators want to estimate mean survival beyond the observed follow-up, they must extrapolate the survival curve.

That requires additional assumptions.

Extrapolated restricted or unrestricted mean survival should therefore be distinguished clearly from an RMST calculated entirely within observed follow-up.

RMST in Regulatory Clinical Trial Reporting

When RMST is included in a clinical trial report, the analysis should be transparent enough that the reader understands:

  • Why RMST was selected.
  • What restriction time was used.
  • How the survival curve was estimated.
  • How the treatment contrast was calculated.
  • How uncertainty was quantified.
  • How missing and censored observations were handled.
  • Whether the analysis was prespecified.

Common Interpretation Questions

Does a larger RMST always mean better treatment?

For a survival endpoint where longer survival is clinically favorable, a larger RMST generally indicates greater average survival through the selected horizon. However, the direction should always be defined according to the endpoint.

Can RMST be calculated if the median is not reached?

Yes. This is one of its important practical advantages.

Does RMST require proportional hazards?

No.

Does RMST replace the Kaplan-Meier curve?

No. It summarizes the curve; it does not replace visual inspection of the underlying survival distribution.

Does RMST replace the hazard ratio?

Not necessarily. RMST and the hazard ratio answer different questions.

The Most Important Mathematical Identity

The central equation to remember is:

\[ \boxed{ \operatorname{RMST}(\tau) = \int_0^\tau S(t)\,dt } \]

Everything else follows from this definition.

For two treatment groups:

\[ \boxed{ \Delta_{\mathrm{RMST}}(\tau) = \int_0^\tau [S_A(t)-S_B(t)]\,dt } \]

And for restricted mean time lost:

\[ \boxed{ \operatorname{RMTL}(\tau) = \tau-\operatorname{RMST}(\tau) } \]

A Compact Worked Example

Suppose a randomized trial follows patients for 24 months. The estimated RMSTs are:

\[ \operatorname{RMST}_A(24)=19.8 \]

and:

\[ \operatorname{RMST}_B(24)=17.1 \]

The treatment difference is:

\[ \Delta_{\mathrm{RMST}}(24) = 19.8-17.1 = 2.7 \]

Therefore:

\[ \boxed{ \Delta_{\mathrm{RMST}}(24)=2.7\text{ months} } \]

If the 95% confidence interval is:

\[ (0.8,\;4.6) \]

then the treatment provides an estimated 2.7 additional months of average survival through 24 months, with the confidence interval indicating that the true difference compatible with the model and data ranges from approximately 0.8 to 4.6 months.

RMST Versus Median: A Final Example

Measure Treatment Control
Median survival Not reached 18.4 months
RMST at 24 months 19.6 months 16.9 months
RMST difference +2.7 months

The median comparison alone gives an incomplete picture because the treatment median has not been reached.

The RMST comparison gives a directly interpretable numerical treatment effect.

Practical Checklist

1
Define the survival endpoint.
2
Define the event and censoring rules.
3
Choose and prespecify \( \tau \).
4
Plot Kaplan-Meier survival curves.
5
Inspect numbers at risk and late follow-up.
6
Estimate RMST in each group.
7
Calculate the RMST difference.
8
Calculate a confidence interval and P-value.
9
Interpret the result in time units.
10
Consider sensitivity analyses and alternative survival estimands when appropriate.

Bottom Line

Bottom line: Restricted Mean Survival Time is the area under the survival curve up to a prespecified time horizon. It provides an average survival-time measure that can be particularly useful when the median survival is not reached or when the proportional hazards assumption is questionable. The primary treatment contrast is often the RMST difference, which can be interpreted directly in units of time: for example, "2.7 additional months of average survival through 24 months." RMST does not require a constant hazard ratio, but its interpretation depends on the selected restriction time. A rigorous RMST analysis should therefore prespecify \( \tau \), account appropriately for censoring, examine the underlying Kaplan-Meier curves, and report the RMST estimates, treatment contrast, confidence interval, and P-value.

References

Royston, P. & Parmar, M.K.B. (2011). The use of restricted mean survival time to estimate the treatment effect in randomized clinical trials when the proportional hazards assumption is in doubt. Statistics in Medicine, 30, 2409–2421.
Uno, H., Claggett, B., Tian, L., et al. (2014). Moving beyond the hazard ratio in quantifying the between-group difference in survival analysis. Journal of Clinical Oncology, 32, 2380–2385.
Uno, H., Wittes, J., Fu, H., et al. (2015). Moving beyond the hazard ratio in quantifying the between-group difference in survival analysis. Journal of Clinical Oncology, 33, 2380–2385.
Trinquart, L., Jacot, J., Blais, C., & Chevret, S. (2016). Time-varying treatment effects and the restricted mean survival time in randomized clinical trials. Clinical Trials.
Royston, P. (2011). Estimating the restricted mean survival time using Stata. The Stata Journal.
Uno, H., Tian, L., Claggett, B., & Wei, L.J. (2014). A versatile test for equality of restricted mean survival times based on the weighted Kaplan-Meier estimators. Statistics in Medicine.