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/