как сделать калибровочную кривую для бинарной логистической модели - PullRequest
Новая
       8

как сделать калибровочную кривую для бинарной логистической модели

0 голосов
/ 25 августа 2018

Я делаю бинарную логистическую регрессию, чтобы построить номограмму. А затем выполните калибровку кривой для повторной выборки калибровки модели. код приведен в списке ниже: 

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 пакетах, которые также сыграли важную роль в получить корректировку для сюжета. Поэтому мне интересно знать, есть ли ошибка, когда я делаю расчет. И некоторые сталкиваются с той же проблемой. Не могли бы вы поделиться своим методом для решения проблемы. Большое спасибо.

...