是什么 R中的coxph模型摘要中给出的值?例如,
Rsquare= 0.186 (max possible= 0.991 )
我愚蠢地把它写成了手稿价值和评论者跳上它说他不知道类似物 来自为 Cox 模型开发的经典线性回归的统计数据,如果有,请提供参考。任何帮助都会很棒!
是什么 R中的coxph模型摘要中给出的值?例如,
Rsquare= 0.186 (max possible= 0.991 )
我愚蠢地把它写成了手稿价值和评论者跳上它说他不知道类似物 来自为 Cox 模型开发的经典线性回归的统计数据,如果有,请提供参考。任何帮助都会很棒!
使用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
是“拟合中使用的观察次数”。
除以总结中的观察次数coxph
是错误的,应该是未经审查的事件数;参见 O'Quigley 等人。(2005 年)解释了比例风险模型中的随机性医学统计 p。479-489。