39、生存分析：Cox比例风险模型及R语言应用

生存分析：Cox比例风险模型及R语言应用

1. 风险比例关系

对公式进行简单变形可得：$h_m(t) = HR \times h_n(t)$ ，从这个公式可以清晰地看出，个体 $m$ 的风险与个体 $n$ 的风险成正比，且比例 $HR$ 与时间无关。

2. 一般风险比的解释

比例风险（PH）假设的有效性使得比较时间相关风险变得简单。我们可以比较想要对比的条件下的任意风险比，将这些条件称为“治疗组”和“对照组”，它们对应的风险表示如下：
- $h(t, X_{treatment}) = h_0(t)\exp(\sum_{i} \beta_iX_{treatment}^i)$
- $h(t, X_{control}) = h_0(t)\exp(\sum_{i} \beta_iX_{control}^i)$

假设PH假设成立，它们的风险比在时间上是恒定的：$HR(T \ vs \ C) = \frac{h(t, X_{treatment})}{h(t, X_{control})}$ 。这个风险比给出了治疗组和对照组在任何时间点 $t$ 的瞬时相对风险，而相对风险（RR）则是对时间累积风险的量化。风险比有以下解释：
|风险比情况|解释|
| ---- | ---- |
|$HR(T \ vs \ C) > 1$|治疗组的风险高于对照组，对照组更有利|
|$HR(T \ vs \ C) = 1$|治疗组和对照组无差异|
|$HR(T \ vs \ C) < 1$|对照组的风险高于治疗组，治疗组更有利|

例如，当 $HR(T \ vs \ C) = 1.3$ 时，治疗