coxph.detail$hazard的公式是什么？

Question

我一直试图确定在coxph.detail(fit)$hazard中使用的公式。我不清楚如何从coxph.detail文档中计算它：

“危险增量。请注意，危险和风险的变化总是针对某一特定的未来主题。这个例程使用对象$mean作为未来的主题。”

Answer 1

回答以避免冗长的评论。从functionBody(survival::coxph.detail)你会看到它正在调用Ccoxdetail

 ff <- .C(Ccoxdetail, ...

然后查看包源，可在

2.41-3.tar.gz

解压缩和解压缩以找到文件survival\src\coxdetail.c

(没有粘贴整个代码以避免重复)。

更重要的是，如果您查看这个小片段，您将看到两个值得阅读的参考资料：

安徒生，P.和吉尔，R. (1982)。Cox回归模型的计数过程，一个大样本的研究。统计年鉴10,1100-1120。 Therneau，T.，Grambsch，P.，模拟生存数据:扩展Cox模型。斯普林格-维拉格，2000年。

Answer 2

对于ties="breslow"，coxph.detail(fit)$hazard是时间t除以exp(fit$linear.predictors)之和的事件数，其中和在处于风险的对象之上。注意，fit$linear.predictors是居中的。