我有以下型号:
coxph(Surv(fulength, mortality == 1) ~ pspline(predictor))
Run Code Online (Sandbox Code Playgroud)
其中,fulength是一个随访的持续时间(包括死亡率),预测因子是死亡率的预测因子.
上面命令的输出是这样的:
coef se(coef) se2 Chisq DF p
pspline(predictor), line 0.174 0.0563 0.0562 9.52 1.00 0.002
pspline(predictor), nonl 4.74 3.09 0.200
Run Code Online (Sandbox Code Playgroud)
如何绘制这个模型,以便在y轴上获得95%置信区间和风险比的漂亮曲线?我的目标是类似于此:
