1

如何从 coxph 获得累积危险比的输出?

我有我的模型

coxfit <- coxph(Surv(combi$survival, combi$realdeath) ~ combi$meanrr, combi)

并获得输出

# Call:
# coxph(formula = Surv(combi$survival, realdeath) ~ combi$meanrr, 
#     data = combi)
# 
#                   coef exp(coef)  se(coef)     z       p
# combi$meanrr -0.004140  0.995869  0.000905 -4.57 4.8e-06
# 
# Likelihood ratio test=21.9  on 1 df, p=2.94e-06
# n= 311, number of events= 70 

hr = exp(-0.004140)我知道我可以使用和手动计算危险比ci= exp(-0.004140-1.96*0.000905)exp(-0.004140+1.96*0.000905)

这将给我增加一个单位的置信区间的 HR。

是否有一个函数可以将结果作为矢量或 data.frame 给我?如果我可以定义感兴趣的单位的数量,比如

somefunction(coxfit, unit_step)

但只计算累积风险比就可以了

4

2 回答 2

2

summary(coxfit) 给出了正确答案

coxph(formula = Surv(combi$survival, realdeath) ~ meanrr, data = combi)

  n= 311, number of events= 70 

             coef  exp(coef)   se(coef)      z Pr(>|z|)    
meanrr -0.0041399  0.9958687  0.0009055 -4.572 4.83e-06 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

       exp(coef) exp(-coef) lower .95 upper .95
meanrr    0.9959      1.004    0.9941    0.9976

Concordance= 0.661  (se = 0.035 )
Rsquare= 0.068   (max possible= 0.92 )
Likelihood ratio test= 21.86  on 1 df,   p=2.938e-06
Wald test            = 20.9  on 1 df,   p=4.83e-06
Score (logrank) test = 21.24  on 1 df,   p=4.051e-06

其中,增加一个单位的风险比为 exp(coef),置信区间为下/上 .95

于 2016-08-05T07:17:29.590 回答
2

basehaz()您可以使用或提取累积风险survfit()

为您的模型键入str(basehaz(coxfit))str(survfit(coxfit))在您的情况下将揭示数据对象的结构,并提供线索basehaz(coxfit)$hazardsurvfit(coxfit)$cumhaz产生您所追求的矢量/数据框友好输出。

帮助说明?basehaz会告诉您这survfit()是首选方法,实际上后者提供了全面的输出,包括 CI。

于 2016-08-04T17:09:53.893 回答