如何从Cox PH模型计算预测的危险率？

11

我有以下Cox PH模型：

（时间，事件）〜X + Y + Z

我想获得预测的危险率（我说的是风险率不危险比）给出的具体值X，Y，Z。我知道muhaz R软件包可以计算观察到的危险率，但是我对预测模型感兴趣。

有没有办法在R中做到这一点？

r survival hazard cox-model

— 没有存货
source

2

看来您的“方程式”模拟了二进制结果。现在，考克斯（Cox）的PH模型处理的是生存数据，它的左手边应该是危险率……我想你应该澄清

— 一下

2

Cox模型不估计基线风险函数（通常表示为

）。要从模型估计发病率，应改用参数生存模型。

λ_{0} (t)

$\lambda_0(t)$

— boscovich 2012年

1

您的问题含糊不清。危害函数是概率，而不是比率。是您想要的危险功能，还是其他？

— Fomite 2012年

16

R生存软件包中用于获得基准风险率的函数为basehaz。

$e^{\beta}$

一个简单的示例可能会有所帮助：

library(survival) #survival analysis
library(eha) #used for data 
data(oldmort) #create the data

# Create surv data set
mort <- Surv(time=oldmort$enter,time2=oldmort$exit,event=oldmort$event)

reg_fit <- coxph(formula=mort~oldmort$sex)
summary(reg_fit)

# Now get baseline curve
baseline <- basehaz(reg_fit)

# Draw baseline hazard (that's male)
plot(baseline$time, baseline$hazard, type='l',main="Hazard rates") 

# Draw female hazard
lines(baseline$time, exp(-0.1929)*baseline$hazard, col="blue")

$\exp(-0.1929)=0.8245$ summary(reg_fit) $18\%$

— 卡尔·奈特
source

13

该函数basehaz（来自上一个答案）提供累积危害，而不是危害函数（比率）。我认为这个问题是关于危害功能的。估计危害函数将需要指定平滑类型（如密度估计）。所述Muhazř包可为一个样本数据执行此操作。我不知道可以针对Cox模型中的基准危害执行此操作的功能。我也需要这个我认为我需要自行消除累积的危害。

— 214
source