2017-02-13 11 views
0

は、ROC曲線を描くために以下のコードを考えてみましょう:パッケージmadaのROCカーブの軸のタイトルを変更するにはどうすればいいですか?

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

このコードは、デフォルトでは、それぞれ「感度」と「偽陽性率」を縦軸と横軸のタイトルを設定します。私は以下のコード使用してタイトルを変更しようとしました:

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

をしかし、これはエラーを与える:

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

は、それがタイトルを変更する権限を与えるものではありませんようです。すべての手がかりは?

答えて

3

ステップバイステッププロットを使用すると、軸を変更することができます。あなたは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

のREADMEに例に沿って続く場合はありがとうございました。これはうまくいく。今、予測領域をプロットに追加する方法はありますか?私はそれに適合した値が 'f = fit.reitsma $ fitted.values'によって保存できることを知っていますが、それをどのように使用するかはわかりません。 –

+0

私はmada以外の出力をデータフレームとして受け取り、ggplot2を使用します。 – thisisrg

関連する問題