royston
位于 survival
包(package)。 说明
计算 Royston 和 Sauerbrei 提出的 D 统计量以及几个伪 R 平方值。
用法
royston(fit, newdata, ties = TRUE, adjust = FALSE)
参数
fit |
考克斯适合 |
newdata |
可选的验证数据集 |
ties |
对风险评分中的联系进行修正 |
adjust |
调整可能的过度拟合 |
细节
这些值称为伪 R-squared,因为它们仅涉及线性预测变量,而不涉及结果。 R.D
是对应于 Royston 和 Sauerbrei 统计数据的值。 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.
相关用法
- R rotterdam Royston 和 Altman (2013) 使用的乳腺癌数据集
- R rats Mantel 等人的大鼠治疗数据
- R residuals.survreg 计算“survreg”对象的残差
- R reliability 可靠性数据集
- R ridge 岭回归
- R ratetableDate 将日期对象转换为费率表形式
- R ratetable 允许模型中使用ratetable()项
- R rhDNase rhDNA酶数据集
- R ratetables 预期生存和人年函数的人口普查数据集
- R retinopathy 糖尿病性视网膜病变
- R residuals.survfit 来自 survfit 对象的 IJ 残差。
- R residuals.coxph 计算“coxph”拟合的残差
- R rats2 大鼠数据来自 Gail 等人。
- R rttright 计算重新分配到正确的权重
- R hoel 小鼠癌症数据
- R survcondense 缩短 (time1, time2) 生存数据集
- R myeloid 急性粒细胞白血病
- R tobin 托宾的托比特数据
- R pseudo 生存的伪值。
- R levels.Surv 返回多状态 Surv 对象的状态
- R diabetic 糖尿病视网膜病变
- R pbc 梅奥诊所原发性胆汁性胆管炎数据
- R plot.survfit survfit 对象的绘图方法
- R kidney 肾导管数据
- R stanford2 更多斯坦福心脏移植数据
注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Compute Royston's D for a Cox model。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。