4
$\begingroup$

I was hoping someone could help me with this problem in the cox proportional hazards model.

I am given the following setup.

T is a non-negative random variable with continous distribution and hazard function $\lambda_T(t)$. T has density $f_T(t) = \lambda_T(t) S(t)$ and $S(t) = P(T>t)$. Also $F(t) = P(T \leq t)$ is the distribution function.

If I have $n$ observations of $T$. Note no censoring is assumed. Can anyone tell me how I arrive at an $\textbf{efficient influence function}$ for $S(t_0)$ where $t_0$ is a fixed time point.

note $\sqrt{n} ( \hat{S(t_0)} - S(t_0) ) = \sum_{i=1}^{n} \phi(T_i)$ where we have $\phi(T_i)$ is the influence function. This leads to an efficient estimator of $\hat{S}$

$\endgroup$
11
  • $\begingroup$ is there a particular text book you are working from in regards to this? $\endgroup$ Apr 14 '21 at 18:46
  • $\begingroup$ what does $O_p$ here represent? $\endgroup$ Apr 14 '21 at 18:51
  • $\begingroup$ $O_p$ is a term converging to zero. $\endgroup$
    – nalen
    Apr 14 '21 at 18:54
  • $\begingroup$ No not a particular book actually. $\endgroup$
    – nalen
    Apr 14 '21 at 18:55
  • $\begingroup$ As this is a Cox regression model, are you asking about estimating the influence of removing one event time-point on the baseline survival function? Put another way, would the individual with the event at that time point still be included at prior times (effectively censored at what was really the event time)? $\endgroup$
    – EdM
    Apr 14 '21 at 21:54
0
$\begingroup$

I'm not much of an expert on influence functions; I'll start with the working definition provided in this answer by Michael Chernick: "The influence function for a parameter...essentially measures the difference between the parameter estimate when the data point is included compared with when it is left out."

In your case you want to know how removing particular event times from the observation set (maybe more precisely, making small changes in observed event times $T_i$) affect an estimate of survival at a particular time, $\hat S(t_0)$. In your situation with a non-parametric survival function estimate,* that might be the Kaplan-Meier estimate, or the survival function derived from the Nelson-Aalen estimate of cumulative hazard. So ask yourself the following questions:

If $T_i > t_0$, is $\hat S(t_0)$ affected if you omit observation $i$ or make a (small) change in its observed time?

If $T_i = t_0$ (an event perhaps of 0 probability in principle, but maybe of some practical interest), what happens to $\hat S(t_0)$ if you omit observation $i$ or make a (small) change in its observed time?

If $T_i < t_0$, what happens to $\hat S(t_0)$ if you omit observation $i$ or make a (small) change in its observation time?

The Wikipedia entry shows the derivation of the Kaplan-Meier estimate based on maximum likelihood, which might help put the above into a more formal argument.

Although you ask in the context of no censoring, also consider what happens to $\hat S(t_0)$ if there are small changes in censoring times that aren't close to $t_0$.


*Although the question was originally posed in terms of a Cox regression, discussion in comments clarified that the question is about a non-parametric estimate of a single survival curve. A "semi-parametric" Cox regression makes no parametric assumptions about the baseline hazard, with parametric modeling of covariate effects on hazard. If the "influence function" is defined in terms of small changes in observed event times with unaltered covariate values, this type of argument can be extended to Cox models. In Cox models, however, the "influence" of interest is generally in how each of $n$ individual cases, with associated covariate values, affects estimates of each of the $p$ regression coefficients.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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