Solved – Estimating median survival times from Kaplan-Meier plot inspection

estimationkaplan-meiermediansurvival

I've gone through the various questions relating to Kaplan-Meier plots and survival estimates, but I haven't really been able to find anything to help with this specific scenario.

Sometimes, when reading medical publications you realise that the authors haven't reported median survival times of subgroups, but have simply plotted them in Kaplan-Meier curves. Drawing a line at 0.5 on the y-axis and dropping a vertical line to the x-axis when said line intersects with the Kaplan-Meier curve will give you a rough estimate of the median survival time and help you interpret the potential "difference", if the paper otherwise isn't too enlightening.

However, how do you estimate the median survival time, if say one subgroup never reaches the 0.5 mark during the analysis period?

Attached is such a situation taken from this paper:

Epidermal Growth Factor Receptor, p53 Mutation, and Pathological Response Predict Survival in Patients with Locally Advanced Esophageal Cancer Treated with Preoperative Chemoradiotherapy, Michael K Gibson et al., Clin Cancer Res 2003;9:6461-6468

a plot

I mean regardless of if the cohort with "no mutation" reaches 0.5 or not, you could still calculate a median survival time for that cohort if you had the raw data!
I appreciate that the method of "drawing a line" and eyeballing isn't exact science is it, but sometimes its the best you can do.

Any pointers helping with this issue would be great and any ideas for a work around would be great.

Many thanks in advance,

Oliver

Best Answer

The time at which the Kaplan-Meier survival curve crosses the 50% line is the non-parametric estimate of the median survival time. If the Kaplan-Meier curve does not cross the 50% line, then the non-parametric estimate is not defined. At this stage, I can see two simple options:

  • use another quantile (e.g. 0.75) to compare the two groups;
  • approximate the survival curve by means of a parametric fit and derive the median survival time using the model.

As a side note, the typical underlying assumption in survival analysis is that any subject will eventually experience the event of interest if he is observed for a sufficient long period of time. In other words, the survival curve will reach zero (and in particular cross the 50% line) if the observation period lasts a sufficient long time. In contrast, cure modelling deals with situations where the survival curve eventually reaches a plateau. In that case, it is not correct to use a parametric model to extrapolate the survival curve below the plateau.

Related Question