2

我正在尝试使用 Python Lifelines 包来校准和使用 Cox 比例风险模型。

所以，结果总结是：

coef  exp(coef)  se(coef)        z      p  lower 0.95  upper 0.95 
PD    -1.1446     0.3183    0.0814 -14.0563 0.0000     -1.3042     -0.9850  ***
oil   -0.1275     0.8803    0.0016 -79.2128 0.0000     -0.1306     -0.1243  ***
curr  -0.1353     0.8735    0.0020 -67.3416 0.0000     -0.1392     -0.1313  ***
matur -0.0002     0.9998    0.0000 -13.6039 0.0000     -0.0002     -0.0002  ***
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Concordance = 0.602

然后我需要计算部分危险，例如，使用我的生存数据选项卡的第一行，即：

PD  oil curr    Durat   binar   matur
0   0.135760    62.799048   59.004243   1.446575    0   179

它必须是这样的：

PD = 0.13576
oil = 62.799048
curr = 59.004243
matur = 179
np.exp(-1.1446*PD - 0.1275*oil - 0.1353*curr -0.0002*matur)

等于 9.387106981409155e-08，所以它是非常小的数字，导致所有 t 的生存概率等于 1.0。但是当我得到一个 cph.predict_partial_hazard(cox_surv) 方法时，它给了我 0.32 之类的值，我认为这是正确的数字。例如，我们有 Baseline SP = 0.7 和 (0.7^0.32) 到部分危险，我们将得到类似 0.892136633056215 的东西，这是正常的。错误是什么？我们如何以正确的方式计算部分风险？非常感谢！

4

1 回答 1

1

生命线的作者在这里。

生命线中的部分危险是通过首先对变量进行贬义来计算的，因此在生命线中，计算类似于

np.exp(-1.1446*(PD-mean_PD) - 0.1275*(oil-mean_oil) - 
          0.1353*(curr-mean_curr) -0.0002*(matur-mean_matur))

这可能会给您带来更大的部分风险（但是所有科目的相对排名保持不变）。

文档字符串没有提到这一点，但是，我会在 v0.15.0 中修复它。

于 2018-11-02T16:47:21.553 回答