R royston 计算 Cox 模型的 Royston D

说明

计算 Royston 和 Sauerbrei 提出的 D 统计量以及几个伪 R 平方值。

用法

royston(fit, newdata, ties = TRUE, adjust = FALSE)

参数

fit	考克斯适合
newdata	可选的验证数据集
ties	对风险评分中的联系进行修正
adjust	调整可能的过度拟合

细节

这些值称为伪 R-squared，因为它们仅涉及线性预测变量，而不涉及结果。 R.D 是对应于 Royston 和 Sauerbrei D 统计数据的值。 R.KO是Kent和O'Quigley提出的值，R.N是Nagelkerke提出的值，C.GH对应Goen和Heller的一致性度量。

D 的调整基于比率 r=(事件数量)/(系数数量)。对于具有足够样本量 (r>20) 的模型，调整会很小。

Nagelkirke 值是 Cox-Snell R-squared 除以缩放常数。这两个单独的值作为 2 元素向量 rsq 出现在 summary.coxph 的结果中，并在旧版本的打印例程中列为 "Rsquare" 和 "max possible"。 (自默认打印输出中被索引取代以来。)当 newdata 存在时，不会返回 Nagelkirke 估计值。

值

包含 D 值、D 的估计标准误差以及三个或四个伪 R-squared 值的向量。

例子

# An example used in Royston and Sauerbrei
pbc2 <- na.omit(pbc)  # no missing values
cfit <- coxph(Surv(time, status==2) ~ age + log(bili) + edema + albumin +
                   stage + copper, data=pbc2, ties="breslow")
royston(cfit)

参考

M. Goen and G. Heller, Concordance probability and discriminatory power in proportional hazards regression. Biometrika 92:965-970, 2005.

N. Nagelkerke, J. Oosting, J. and A. Hart, A simple test for goodness of fit of Cox's proportional hazards model. Biometrics 40:483-486, 1984.

P. Royston and W. Sauerbrei, A new measure of prognostic separation in survival data. Statistics in Medicine 23:723-748, 2004.