급하면 코딩이 잘 되는 법이다. AUC를 구해 달라는 연락을 받고 인터넷 서핑을 통하여 해결했다. AUC 구하는 패키지는 많이 있는데, 잘 알려진 pROC 패키지를 이용했다.
dt %>% roc('group', 'variable')
plot(roc1, print.auc=T, print.thres=T)해석은 variable이 34.654일 때 sensitivity 0.850, specifity 0.897, AUC 0.927이다.
AUC의 95% 신뢰구간과 p-value는 따로 구한다.
ci(roc1) #95% CI
var(ci(roc1)) # p-value