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:在box(...)中:“添加”不是图形参数
6:在title(...)中:“添加”不是图形参数


还是为此目的有其他方法吗?我也想为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:在轴上(side = side,at = at,labels = labels,...):
“ par.corrected”不是图形参数4:在轴上(侧=侧,在=
处,标签=标签,...):“ par.corrected”不是图形的
参数5:在box(...)中:“ par.corrected”不是图形参数
6:在title(...)中:“ par.corrected”不是图形参数

最佳答案

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

关于r - 如何将多个校准曲线添加到同一图,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/42237268/

10-12 17:53