如何为二元逻辑模型制作校正曲线

时间:2018-08-25 02:51:51

标签: r calibration

我进行了二进制逻辑回归来构建诺模图。然后绘制校准曲线以重采样模型校准。代码如下:

library(rms)
f1 <- lrm(Recurrence_1_year~AFP_Categ_200+Tumor_Diameter_Categ+Macro_Tumor_Thrombus,
           data = Training_nom,
           x=TRUE,
           y=TRUE)     
plot(cal,
  subtitles=FALSE,                    
  conf.int=TRUE,                              
  cex.subtitles=.75,                  
  riskdist=FALSE,                     
  add=FALSE,
  scat1d.opts=list(nhistSpike=1000), 
  par.corrected=list(col="red",lty=4,lwd=4))

运行代码时,R-studio界面中存在许多警告,这些警告类似于:

1: In plot.window(...) : "conf.int" not graphic paremeters
2: In plot.window(...) : "cex.subtitles"not graphic paremeters
3: In plot.window(...) : "riskdist"not graphic paremeters
4: In plot.window(...) : "add"not graphic paremeters
5: In plot.window(...) : "par.corrected"not graphic paremeters
6: In plot.xy(xy, type, ...) : "conf.int"not graphic paremeters
7: In plot.xy(xy, type, ...) : "cex.subtitles"not graphic paremeters
8: In plot.xy(xy, type, ...) : "riskdist"not graphic paremeters
9: In plot.xy(xy, type, ...) : "add"not graphic paremeters
10: In plot.xy(xy, type, ...) : "par.corrected"not graphic paremeters
11: In axis(side = side, at = at, labels = labels, ...) : "conf.int"not graphic paremeters
12: In axis(side = side, at = at, labels = labels, ...) : "cex.subtitles"not graphic paremeters
13: In axis(side = side, at = at, labels = labels, ...) : "riskdist"not graphic paremeters
14: In axis(side = side, at = at, labels = labels, ...) : "add"not graphic paremeters
15: In axis(side = side, at = at, labels = labels, ...) :  "par.corrected"not graphic paremeters
16: In axis(side = side, at = at, labels = labels, ...) : "conf.int"not graphic paremeters
17: In axis(side = side, at = at, labels = labels, ...) : "cex.subtitles"not graphic paremeters
18: In axis(side = side, at = at, labels = labels, ...) : "riskdist"not graphic paremeters
19: In axis(side = side, at = at, labels = labels, ...) : "add"not graphic paremeters
20: In axis(side = side, at = at, labels = labels, ...) : "par.corrected"not graphic paremeters
21: In box(...) : "conf.int"not graphic paremeters
22: In box(...) : "cex.subtitles"not graphic paremeters
23: In box(...) : "riskdist"not graphic paremeters
24: In box(...) : "add"not graphic paremeters
25: In box(...) : "par.corrected"not graphic paremeters
26: In title(...) : "conf.int"not graphic paremeters
27: In title(...) : "cex.subtitles"not graphic paremeters
28: In title(...) : "riskdist"not graphic paremeters
29: In title(...) : "add"not graphic paremeters
30: In title(...) : "par.corrected"not graphic paremeters

警告中的所有参数列表也在 rms 包中的 calibrate 功能指南中列出。在调整情节方面起着重要作用。因此,我想知道在进行计算时是否有任何错误。还有一些人遇到相同的问题。您能否分享解决问题的方法。非常感谢。

0 个答案:

没有答案