R语言 生存预测模型Calibration曲线
使用rms包中的calibrate函数对生存预测模型进行校准曲线绘制,我的代码如下:
modelforcalibration<-coxph(Surv(time,censor_embolism)~cumulative_risk_stroke,data=dataforcalibration,x=1)
cal<-calibrate(modelforcalibration, cmethod = "KM", method = "boot", u = 1638, m = 400, B = 1000)
plot(cal)
得到了一张校准曲线图,横纵坐标是预测生存概率和实际生存概率
类似这张图:
我的问题:如何能绘制横纵坐标分别是预测事件发生率和实际事件发生率的calibration曲线呢?可以通过参数设置来改变吗?谢谢!