Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Topic: What exactly the ROC curve can tell us or can be inferred?
Replies: 0

 Aaronne Posts: 110 Registered: 6/2/11
What exactly the ROC curve can tell us or can be inferred?
Posted: Mar 18, 2013 8:50 AM

Hi Smart Guys,

I wrote some codes to run a linear discriminant analysis based classification:

%% Construct a LDA classifier with selected features and ground truth information
LDAClassifierObject = ClassificationDiscriminant.fit(featureSelcted, groundTruthGroup, 'DiscrimType', 'linear');
LDAClassifierResubError = resubLoss(LDAClassifierObject);

Thus, I can get

Resubstitution Error of LDA (Training Error): 1.7391e-01
Resubstitution Accuracy of LDA: 82.61%
Confusion Matrix of LDA:
14 3
1 5

Then I run a ROC analysis for the LDA classifier:

% Predict resubstitution response of LDA classifier
[LDALabel, LDAScore] = resubPredict(LDAClassifierObject);

% Fit probabilities for scores (the groundTruthGroup contains lables either 'Good' or 'Bad')
[FPR, TPR, Thr, AUC, OPTROCPT] = perfcurve(groundTruthGroup(:,1), LDAScore(:,1), 'Good');

I have got:

OPTROCPT =

0.1250 0.8667

Therefore, we can get:

Accuracy of LDA after ROC analysis: 86.91%
Confusion Matrix of LDA after ROC analysis:
13 1
2 7

My questions are:

1. After ROC analysis we obtained a better accuracy, when we report the accuracy of the classifier, which value we should use? What exactly the ROC curve can tell us or can be inferred? Can we say after ROC analysis we found a better accuracy of the LDA classifier?

2. Why the ROC can produce a better accuracy for the classifier, but the original ClassificationDiscriminant.fit can't?

3. I have also done a cross validation for the LDA classifier, like

cvLDAClassifier = crossval(LDAClassifierObject, 'leaveout', 'on');

Then how to get the ROC analysis for the cross validation? 'resubPredict' method seems only accept 'discriminant object' as input, then how can we get the scores?

4. classperf function of Matlab is very handy to gather all the information of the classifier, like

%% Get the performance of the classifier
LDAClassifierPerformace = classperf(groundTruthGroup, resubPredict(LDAClassifierObject));

However, anyone knows how to gather these information such as accuracy, FPR, etc. for the cross validation results?

Thanks very much. I am really looking forward to see the reply to above questions.

A.