Competing Risks (竞争风险)

参考资料:

1. Competing Risks(Germ´an Rodr´ıguez)

2. Competing Risks - What, Why, When and How? (Sally R. Hinchliffe)

1. 简介与例子引入

competing risks指造成failure的原因有多种可能。

以放置节育环为例，这里的failure指的应该是节育环没能在原定期限前一直发挥好作用(即避孕作用)。导致failure的risk包括意外怀孕、节育环排出、医学原因移除、个人原因移除等；由不同原因造成的failure，视为不同类的failure。

对于这类问题，研究的兴趣点主要在于以下三项:

1) 研究变量向量$x$和特定failure的risk(指造成failure的为特定原因的风险risk)之间的关系；

2) 分析在某类failure高风险的人群是否也存在另一类failure的高风险；

3) 估计某类failure在排除其他原因的情况下的risk。

如果有不同类failure相互独立的假设，问题3)可以被回答被计算(独立的话，条件概率和自己本身的概率是一样的)，而问题2)会因此得到否定的回答。

 

2. 标记说明和基本公式

$T$: Survival Time(生存时间)

$J$: The type of failure(failure的类型，即因为何种原因造成的failure).$j\in \{1,2,...,M\}$

$x$: A vector of covariate(协变量的向量)

总体风险比率(overall hazard rate)仍像一般情形的定义:

\[\lambda(t,x)=\lim_{\Delta t \to 0}\frac{P\{t\leq T < t+\Delta t|T\geq t, x\}}{\Delta t}\]

特定原因导致failure的瞬时风险(cause-specific hazard rate)则为:

\[\lambda_j(t,x)=\lim_{\Delta t \to 0}\frac{P\{t\leq T < t+\Delta t, J=j|T\geq t, x\}}{\Delta t}\]

这里假设failure只能由其中一种原因导致。如果存在两种以上failure同时发生(同时有两种最开始定义的原因导致了failure)，把这种情况定义成新的一类failure(两种原因同时出现时，定为新的原因，和原来的两种独立造成failure的情况区分开来)，使假设依然成立。

这样子，就有

\[\lambda(t,x)=\sum_{j=1}^{m} \lambda_j(t,x)\]

根据上式及生存函数的定义得总体生存函数(overall survival function)为

\[S(t,x) = e^{-\Lambda(t,x)} = e^{-\int_0^t \lambda(u,x)du}\]

(当$x$为time-varying变量时，函数形式会有变化)

类似地，定义

\[S_j(t,x) = e^{-\Lambda_j(t,x)},\]

\[\Lambda_j(t,x)=e^{-\int_0^t \lambda_j(u,x)du}\]

$*$当failure类型数$M>1$时，一般来说，$S_j(t,x)$不能被理解成生存函数survival function，但可以推出$S_j(t,x)$和$S(t,x)$之间的关系:

\begin{eqnarray*}
S(t,x) & = & e^{-\int_0^t \lambda(u,x)du} \\
& = & e^{-\int_0^t \sum_j \lambda_j(u,x)du} \\
& = & e^{-\sum_j \int_0^t \lambda_j(u,x)du} \\
& = & \prod_j e^{-\int_0^j \lambda_j(u,x)du} \\
& = & \prod_j S_j(t,x)
\end{eqnarray*}

同样地，定义cause-specific的t时刻密度函数:

\[f_j(t,x)= \lim_{\Delta t \to 0} \frac{P\{t\leq T \leq t+\Delta t, J=j|x\}}{\Delta t}\]

\[f_j(t,x)= \lambda_j(t,x)S(t,x).\]