4
$\begingroup$

I use Cox regression (proportional hazards) to model survival for a cohort of patients. Patients are censored (alive (0), dead (1)).

I was wondering how Cox regression uses censored data intuitively. I thought when alive (0), Cox model will just ignore them, but apparently it is not so simple.

For example, I used the Cox model with all patients (both alive and dead), and the Cox model with only the patients annotated as dead. Apparently, the results were quite different (likelihood ratio test for the whole model, Wald test for the individual covariates).

In brief, intuitively, how does the Cox model uses the censored data and how the censoring affects the results?

$\endgroup$
1
  • 1
    $\begingroup$ If you're unfamiliar with how censored data are used in survival models it might be better to start with either a Kaplan Meier (or even a parametric estimate of survival) for a single group $\endgroup$
    – Glen_b
    Commented Aug 21, 2018 at 13:05

2 Answers 2

4
$\begingroup$

Inference in survival analysis is based on the survival likelihood.

The survival likelihood differs from the classical likelihood due to the presence of censored data:

  • the relevant information contained in an uncensored data ($\delta = 1$) is that the event occurred at the observed time ($y$); such data contributes to the likelihood via the density function, $f(y)$ (just like in the classical likelihood);
  • the relevant information contained in a censored data ($\delta = 0$) is that the event time exceeds the censoring time $y$; such data thus contributes to the likelihood via the survival function, $S(y)$.

Thus, the survival likelihood for a sample of size N is $$ \prod_{i=1}^N f(y_i)^{\delta_i}\, S(y_i)^{1 - \delta_i} $$


Note 1: The above likelihood is derived under working assumptions (e.g. independent and non-informative censoring).

Note 2: The Cox method (semi-parametric) uses a partial likelihood rather than the above likelihood. However, the idea was more easily explained like that.

$\endgroup$
1
0
$\begingroup$

Thanks ocram and Glen_b,

I tried the Kaplan Meier and I think there I understood what is going on. The censored observation do not seem to affect the KM curve at all at the point they occur but only at the next uncensored observation.

In Cox regression, I understand that both censored and uncensored affect the result. So, based on ocram's answer they seem to affect the survival likelihood in different ways (density or survival function). Can you explain intuitively why it has to be defined this way?

Also survival function gives the probability that a patient will survive beyond any given specified time?

Thank you

$\endgroup$
1
  • 1
    $\begingroup$ Your statement about KM is wrong. The magnitude of the steps depends on the amount of censoring. You should start by reading any introductory textbook on survival analysis. Good luck. $\endgroup$
    – ocram
    Commented Aug 22, 2018 at 5:22

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.