是什么 ”R2R2" R 中 coxph 模型摘要中给出的值

机器算法验证 r 生存 r平方 cox模型
2022-01-22 11:40:11

是什么R2 R中的coxph模型摘要中给出的值?例如,

Rsquare= 0.186   (max possible= 0.991 )

我愚蠢地把它写成了手稿R2价值和评论者跳上它说他不知道类似物R2 来自为 Cox 模型开发的经典线性回归的统计数据,如果有,请提供参考。任何帮助都会很棒!

2个回答

使用getS3method("summary","coxph")您可以查看它是如何计算的。

相关代码行如下:

logtest <- -2 * (cox$loglik[1] - cox$loglik[2])
rval$rsq <- c(rsq = 1 - exp(-logtest/cox$n), maxrsq = 1 - 
        exp(2 * cox$loglik[1]/cox$n))

这里cox$loglik是“一个长度为 2 的向量,其中包含具有初始值和系数最终值的对数似然”(参见 参考资料?coxph.object),并且cox$n是“拟合中使用的观察次数”。

除以n总结中的观察次数coxph是错误的,应该是未经审查的事件数;参见 O'Quigley 等人。(2005 年)解释了比例风险模型中的随机性医学统计 p。479-489。