2017-02-13 1 views
0

den Code Betrachten Sie unten eine ROC-Kurve zu zeichnen:Wie ändert man Achsentitel der ROC-Kurve im Paket mada?

library(mada) 
data(AuditC) 
fit.reitsma <- reitsma(AuditC) 
plot(fit.reitsma) 

Dieser Code setzt die vertikalen und horizontalen Achsen Titel „Empfindlichkeit“ und „False Positive Rate“ bzw. standardmäßig. Ich habe versucht, die Titel zu ändern unter dem Code:

plot(fit.reitsma, xlab="Specifity") 

Aber das gibt den Fehler:

Error in plot.default(c(2, 2), ylim = ylim, xlim = xlim, xlab = "False Positive Rate", : 
formal argument "xlab" matched by multiple actual arguments 

scheint es keine Berechtigung geben die Titel zu ändern. Irgendwelche Hinweise?

Antwort

3

Mit dem Schritt-für-Schritt-Diagramm können Sie die Achsen ändern. Wenn Sie an dem Beispiel in der Readme von reitsma-class {mada}

#calculate a SROC curve 
my.sroc <- sroc(fit.reitsma) 
# calculate a confidence region but do not plot it 
my.cr <- ROCellipse(fit.reitsma) 
# plot the confidence region in ROC space as a line 
plot(my.cr$ROCellipse, type = "l", xlim = c(0,1), ylim = c(0,1), xlab="Specificity", ylab="Sensitivity") 
lines(my.sroc) 

enter image description here

+0

Vielen Dank folgen. Das funktioniert gut. Hast du jetzt eine Idee, wie man die Vorhersage-Region zum Plot hinzufügen kann? Ich weiß, dass die angepassten Werte dafür durch 'f = fit.reitsma $ passed.values' gespeichert werden können, aber keine Ahnung haben, wie man es benutzt. –

+0

Ich würde nur die Ausgänge außerhalb von Mada als Datenframe nehmen und ggplot2 verwenden. – thisisrg

Verwandte Themen