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.