参考资料:
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).\]