ROC receiver operator curve and derived statistics. [...] = roc(d,c); d DATA c CLASS, vector with 0 and 1 [...]=roc(d1,d0); d1 DATA of class 1 d2 DATA of class 0 [SEN, SPEC, TH, ACC, AUC,Yi,idx]=roc(...); OUTPUT: SEN sensitivity SPEC specificity TH Threshold ACC accuracy AUC area under ROC curve Yi max(SEN+SPEC-1), Youden index c TH(c) is the threshold that maximizes Yi