具有13个独立协变量的Cox模型的r-偏差残差诊断-有什么方法可以得到每个协变量的图吗



我是生存统计和这些论坛的新手,所以请原谅我不知道行话或正确的查询方式。我会全力以赴。

我的目标是进行一些考克斯回归,并想对其拟合优度进行一些基本的诊断。我有大约13个协变量,当我进行ggcoxdiagnostics时,我能够得到偏差残差图。然而,这个图代表了整个考克斯模型,我发现它是扭曲的。因此,我想对Cox模型中的每一个协变量进行同样的分析,以了解哪个协变量突出,并且可能需要变换(对数、三次样条曲线、多项式等)。

我根本不知道该怎么做。我应该使用什么代码?

库("surfminer")

res.cox=髋关节(Surv(TimeAxis,new_totmot)~pwvX_AUS+MAP_PWV_AUS+A_HRcarX_AUS+sex_AUS+age_scr_AUS+BMI_AUS+glukos_0_AUS+Current_smoking_AUS+BPlowering_AUS+Ripidlowing_AUS+Hypertension_AUS+Diabetes_confirmed_AUS+Had_CV_Ttreatment+FH,数据=A)

ggcoxdiagnostics(res.cox

可以

(1)试图通过使用其他类型的残差(=分数或Schoenfeld残差)为每个变量获得单独的残差。有了它们,每个变量都会返回一个值,所以这应该有效,而且使用它们非常简单:

ggcoxdiagnostics(res.cox, type ="score")
ggcoxdiagnostics(res.cox, type ="schoenfeld")

(2)坚持使用偏差残差,一次省略或修改/转换一个变量,并检查残差是否以图形方式改善了结果。

不能-为具有偏差残差的受试者获得多个残差:对于偏差残差,您只能为每个受试者返回一个残差值。文件中指出,"…对于…偏差残差,返回的对象是一个向量,每个主题有一个元素"。

最新更新