我正在尝试使用rms
包绘制受限三次样条模型。但是我找不到任何方法来调整我的 cox 比例风险模型,我只能得到未调整的拟合。这是我的代码:
library(survival)
library(rms)
dd <- datadist(Cox9)
options(datadist="dd")
fit <- cph(Surv(follcox,evento) ~ rcs(G_VINO,3))
plot(Predict(fit_vino), lty=1, lwd=3, ylim=c(-0.5,1.0),xlim = c(0,50), col="white")
通过这个编码,我得到了未调整的样条模型。我想知道如何添加混杂变量来调整模型。我试过了:
fit_vino_adj <- cph(Surv(follcox,evento) ~rcs(G_VINO+edad0+actfis+energia))
plot(Predict(fit_vino_adj), lty=2, lwd=2)
但这给了我每个变量的样条模型,任何人都知道如何调整我的模型?