我一直在尝试确定coxph.detail(fit)$hazard
中使用的公式。我不清楚它是如何从coxph.detail
doc:
"危险增量。请注意,危险的危害和方差总是针对某些特定的未来主题。此例程使用对象$ mean作为未来主题。"
答案 0 :(得分:0)
回答避免冗长的评论。从functionBody(survival::coxph.detail)
开始,您会看到它正在呼叫Ccoxdetail
:
ff <- .C(Ccoxdetail, ...
然后查看包源,可从
获取https://cran.r-project.org/src/contrib/survival_2.41-3.tar.gz
解压缩并解压缩以查找文件survival\src\coxdetail.c
(不粘贴整个代码以避免重复)。
更重要的是,如果你看一下插图,你会看到两个值得一读的参考文献:
Andersen,P。和Gill,R。(1982)。考克斯的计数过程回归模型,大样本研究。统计年鉴10,1100-1120。
Therneau,T.,Grambsch,P.,Modeling Survival Data:扩展Cox模型。施普林格出版社, 2000。
答案 1 :(得分:0)
对于tie =“breslow”,coxph.detail(fit)$ hazard是在时间t的事件数除以exp的总和(拟合$ linear.predictors),其中总和超过风险受试者。请注意,拟合$ linear.predictors是CENTERED。