如何在同一个图中添加几条校准曲线答案

【问题标题】：how can I add several calibrate curve to the same plot如何在同一个图中添加几条校准曲线
【发布时间】：2017-02-14 22:15:46
【问题描述】：

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：在 plot.xy(xy, type, ...) 中：“add”不是图形参数
3：在轴（边=边，在=处，标签=标签，...）：“添加”不是图形参数
4：在轴（边=边，在=处，标签=标签，...）：“添加”不是图形参数
5：在框（...）中：“添加”不是图形参数
6：在title(...)中：“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.corrected”不是图形参数 2：在 plot.xy(xy, type, ...) 中：“par.corrected”不是图形参数 3：在轴（边 = 边，在 = 在，标签 = 标签，...）：
“par.corrected”不是图形参数 4：在轴（边 = 边，在 = 在，标签 = 标签，...）：“par.corrected”不是图形参数 5：在框（...）中：“par.corrected”不是图形参数 6：在标题（...）中：“par.corrected”不是图形参数

【问题讨论】：

标签： r logistic-regression rms

【解决方案1】：

rms 包中的plot.calibrate 函数没有实现add=TRUE，部分原因是显示风险分布的直方图会发生冲突。通过查看代码，您会发现编写自己的绘图函数来制作您需要的绘图并不难。

【讨论】：

【解决方案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"])

【讨论】：