如何將幾條校準曲線添加到同一圖上

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（...）： 「添加」 是不是圖形參數
2：在plot.xy（XY，型，...）：「add」不是圖形參數
3：在軸（side = side，at = at，labels = labels，...）中：「add」不是 圖形參數
4：In （...）：「add」不是 圖形參數
5：在框（...）中：「add」不是圖形參數ETER
6：在標題（...）：「添加」不是圖形參數

還是有什麼其他的方法用於此目的？我也想爲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，類型。 ..）：「par.corrected」不是圖形參數 3：在軸（side = side，at = at，labels = labels，...）中：
「par.corrected」不是圖形參數4：In軸（邊=邊，在= at，標籤=標籤，...）：「par.corrected」不是圖形 參數5：在框（...）中：「par.corrected」不是圖形參數 6：在標題（...）：「par.corrected」不是圖形參數

Answer 1

的plot.calibrate功能在rms包不執行add=TRUE，部分原因是直方圖，顯示風險分佈會發生碰撞。通過查看代碼，您將看到，編寫自己的繪圖功能來製作所需的繪圖並不困難。