R 如何在同一绘图中添加多条校准曲线

警告信息：

1:在plot.window（…）：“par.corrected”不是图形参数2： 在plot.xy（xy，类型，…）中：“par.corrected”不是图形参数 3：在轴中（侧=侧，at=at，标签=标签，…）：
“校正部分”不是图形参数4：轴内（侧=侧，at= 在，labels=labels，…）：“par.corrected”不是图形化的 参数5:框（…）：“par.corrected”不是图形参数 6:在标题（…）：“par.corrected”不是图形参数

---

rms

软件包中的

plot.calibrate

功能未实现

add=TRUE

，部分原因是显示风险分布的直方图会发生冲突。通过查看代码，您会发现编写自己的绘图函数来绘制所需的绘图并不困难。

以下内容允许我在f1原始绘图的顶部绘制f2的偏差校正校准曲线。通过检查由

calibrate

创建的对象的内容，可以查看有关打印其他组件的信息

  plot(calibrate(f2),par.corrected=list(col="blue",lwd=2))

 #n=100   Mean absolute error=0.061   Mean squared error=0.00541
 #0.9 Quantile of absolute error=0.109

set.seed（888）
Y
  plot(calibrate(f2),par.corrected=list(col="blue",lwd=2))

 #n=100   Mean absolute error=0.061   Mean squared error=0.00541
 #0.9 Quantile of absolute error=0.109

set.seed(888)
y<-sample(c(0,1), 100, replace = TRUE)
x1<-rnorm(100,2,3)
x2<-rnorm(100,1,4)
library(rms)
f1<-lrm(y~x1,,x=TRUE,y=TRUE)
f2<-lrm(y~x2,,x=TRUE,y=TRUE)
f1.cal<-calibrate(f1)
f2.cal<-calibrate(f2)
plot(f1.cal)
lines(f2.cal[,"predy"], cal.plot[,"calibrated.corrected"])