我试图找到在 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
鉴于残差是从变量本身创建的,为什么残差和变量的长度不同?