gikoha’s blog

個人的メモがわり

RによるROC曲線

先ほどのデータを解析

source("http://aoki2.si.gunma-u.ac.jp/R/src/ROC.R", encoding="euc-jp")
ROC0(acs,survive)

AUCが最後に表示される

...
$c.index
[1] 0.6589782

ここによるとあんまり尤もらしい曲線じゃないけどね!

Epiパッケージを使うともっとすごい

require(Epi)
Epi::ROC(d2$LHratio,d2$ACS,plot="ROC")

AUCはおろか、cut-off valueまで表示してくれます