2017-02-14 4 views
2
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) 

mit dem obigen Code, es soll beide kalibrieren Plots auf die gleiche Figur hinzufügen. Es ist jedoch berichtet Fehler:Wie kann ich mehrere Eichkurven zum gleichen Diagramm hinzufügen

Warnmeldungen:

1: In plot.window (...): "add" ist nicht grafische Parameter
2: In plot.xy (xy, Typ , ...): "add" ist kein grafischer Parameter
3: In der Achse (Seite = Seite, bei = at, Etiketten = Etiketten, ...): "add" ist nicht grafischer Parameter 4: In Achse (side = Seite, at = at, Labels = Labels, ...): "add" ist nicht graphischer Parameter
5: In Box (...): "add" ist kein grafischer Parameter eter
6: Titel (...): "add" ist nicht grafische Parameter

Oder gibt es andere Methoden für diesen Zweck? Ich möchte auch verschiedene Farben für x1 und x2 Kurven angeben. Ich kann das Argument par.corrected=list(col="blue") nicht angeben. Es funktioniert nicht.

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 

Warnmeldungen:

1: In plot.window (...): "par.corrected" ist nicht grafische Parameter 2: In plot.xy (xy, Typ. ..): "par.corrected" ist kein grafischer Parameter 3: In Achse (Seite = Seite, at = at, Etiketten = Etiketten, ...):
"par.corrected" ist kein grafischer Parameter 4: In Achse (Seite = Seite, bei = bei, Etiketten = Etiketten, ...): "par.corrected" ist nicht grafisch Parameter 5: Im Feld (...): "par.corrected" ist kein grafischer Parameter 6: Titel (...): „par.corrected“ ist nicht grafischer Parameter

Antwort

1

Die plot.calibrate Funktion im rms Paket nicht add=TRUE implementieren, zum Teil, weil Histogramme Risikoverteilungen zeigen würden kollidieren. Wenn Sie sich den Code anschauen, werden Sie feststellen, dass es nicht schwer ist, eine eigene Plotfunktion zu erstellen, um die benötigten Plots zu erstellen.

Verwandte Themen