質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

Q&A

解決済

1回答

1357閲覧

ROCカーブの描画とAUCの計算のプログラム

t.enigma

総合スコア0

R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

0グッド

0クリップ

投稿2021/06/11 22:34

function (val, y, N=0) { # val は判別関数 F(x) の値,y はクラスラベル m = sum(y == 0) #y=0 のサンプル数
n = sum(y == 1) #y=1 のサンプル数
TPR = 1 # TPRの初期値
FPR = 1 # FPRの初期値
AUC=0 #AUC の初期値
uni_val=sort(unique(val)) #閾値 c の集合 plot(1,1,type="n",xlim=c(0,1),ylim=c(0,1),xlab="FPR",ylab="TPR") #描画 for (i in 1:length(uni_val)) {
cri=uni_val[i]
FPR[i+1]=?
TPR[i+1]=?
AUC=?
segments(FPR[i],TPR[i],FPR[i+1],TPR[i+1]) #ROC カーブを描画
}
title(paste("ROC curve","(T=",N,")")) #図にタイトルを追加 text(0.8,0.2,paste("AUC=",round(AUC,3))) # 図に AUC の値を描画 return(list(AUC=AUC))
}

このコードの?の部分に何をいれればいいのかわからないのでお時間があれば教えてください。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

1T2R3M4

2021/06/12 01:46

引用も出典も書いてないのでご自身で書かれたものなのでしょうが、 設計ではどのような計算をすることになっているのでしょうか。
t.enigma

2021/06/12 12:53

自分で書きました。
1T2R3M4

2021/06/12 12:55

設計はどうなっているのですか。
t.enigma

2021/06/12 13:51

正直あまり理解できず友人にききながらできたので設計の部分はよくわかっていません。申し訳ございません。
guest

回答1

0

自己解決

function (val, y, N=0)
{
m = sum(y == 0) #y=0 のサンプル数
n = sum(y == 1) #y=1 のサンプル数
TPR = 1 # TPR の初期値
FPR = 1 # FPR の初期値
AUC=0 #AUC の初期値
uni_val=sort(unique(val)) #閾値 c の集合
plot(1,1,type="n",xlim=c(0,1),ylim=c(0,1),xlab="FPR",ylab="TPR") #描画
for (i in 1:length(uni_val))
{
cri=uni_val[i]
FPR[i+1]=ifelse(i==length(uni_val),0,sum(y[-c(0:cri)]==0))/m
TPR[i+1]=ifelse(i==length(uni_val),0,sum(y[-c(0:cri)]==1))/n
AUC= AUC + (FPR[i]-FPR[i+1])*TPR[i]*ifelse(TPR[i]==TPR[i+1],1,1/2)
segments(FPR[i],TPR[i],FPR[i+1],TPR[i+1]) #ROC カーブを描画
}
title(paste("ROC curve","(T=",N,")")) #図にタイトルを追加
text(0.8,0.2,paste("AUC=",round(AUC,3))) # 図に AUC の値を描画
return(list(AUC=AUC))
}

このコードで解決しました。丸投げな質問になってしまい申し訳ございません。

投稿2021/06/12 12:54

t.enigma

総合スコア0

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問