我正试图找到最佳的变量转换,以在Cox比例风险模型的对数风险或对数累积风险中产生线性。作为这些研究数据,我无法在此处发布。为此,我正在尝试使用此代码针对残差绘制变量。
cox_mod_spline = coxph(Surv(timespan_censored,status)~ risk_factor, data = df)
res = residuals(cox_mod_spline, type = "martingale")
df$risk_factor
plot(na.omit(df$risk_factor), res)
但是我收到此错误消息:Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differ
实际上是当我输入此代码时:
length(df$risk_factor)
length(res)
我知道
[1] 587
[1] 577
分别
我还检查了NA
中没有df$risk_factor
鉴于残差是从变量本身创建的,为什么残差和变量的长度不同?