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曲线呢?可以通过参数设置来改变吗?谢谢!