我一直在试图确定在coxph.detail(fit)$hazard
中使用的公式。它是如何从coxph.detail
DOC计算目前尚不清楚对我说:。
“的危险增加注意风险的危害和方差,始终是一些特定的未来受到这个程序使用对象$平均为未来的主题“。
我一直在试图确定在coxph.detail(fit)$hazard
中使用的公式。它是如何从coxph.detail
DOC计算目前尚不清楚对我说:。
“的危险增加注意风险的危害和方差,始终是一些特定的未来受到这个程序使用对象$平均为未来的主题“。
回答以避免冗长的评论。从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
(不粘贴整个代码以避免重复)。
更重要的是,如果你看看小插曲,你会看到两个引用是值得一读:
安德森,体育和吉尔,R.(1982)。 Cox的计数过程回归模型,一个大样本研究。统计年鉴10,1100-1120。生存数据建模:扩展考克斯模型。施普林格出版社,2000年
关系= “布瑞斯罗夫”,coxph.detail(FIT)$危险是在时间t的事件由exp的总和除以数量(适合$ linear.predictors),总和超过了风险主体。请注意,fit $ linear.predictors是CENTERED。
哪个R包? “生存”? – rbm
如果是这样,那么你总是可以通过'functionBody(survival :: coxph.detail)' – rbm
查看实际代码调用functionBody(survival :: coxph.detail)并不显示它的计算方式,因为C函数被调用来计算它。 – Bill