如何从 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)
但只计算累积风险比就可以了