R-如何在同一图中添加几个校准曲线

Question

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)
plot(calibrate(f1))
plot(calibrate(f2),add=TRUE)

使用上述代码，应该将两个校准图添加到同一图中。但是，它报告了错误：

警告消息：

1：在plot.window（...）中：" add"不是图形参数
2：in plot.xy（xy，type，...）：" add"不是图形参数
3：在轴上（侧=侧，at = at，labels = labels，...）：" add"不是 图形参数
4：在轴上（侧=侧，at = at，labels = labels，...）：" add"不是 图形参数
5：在框（...）中：" add"不是图形参数
6：在标题（...）中：" add"不是图形参数

还是为此目的还有其他方法？我也想为X1和X2曲线指定不同的颜色。我不知道如何指定参数par.corrected=list(col="blue")。它不起作用。

  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

警告消息：

1：在plot.window（...）中：" par.concrected"不是图形参数2： 在plot.xy（xy，type，...）中：" par.concorced"不是图形参数 3：在轴上（侧=侧，at = at，labels = labels，...）：
" PAR.CORNECT"不是图形参数4：在轴（侧=侧，at = AT，标签=标签，...）：" par.concorced"不是图形 参数5：在框（...）中：" par.conconcrect"不是图形参数 6：在标题（...）中：" par.concorced"不是图形参数

Answer 1

rms软件包中的plot.calibrate功能未实现add=TRUE，部分是因为显示风险分布的直方图将碰撞。通过查看代码，您会发现编写自己的绘图功能并不难以制作所需的图。

Answer 2

以下允许我在F1的原始图上绘制F2的偏置校准校准曲线。通过检查calibrate创建的对象的内容，您可以看到有关绘制其他组件的内容。

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"])