fpr = [];
mcc = [];
numSamples = length(scores);
conditionPositive = sum(labels(:));
conditionNegative = numSamples - conditionPositive;

for tmpThreshold = linspace(0,1,100)

    tmpScores = scores > tmpThreshold;
    predictedPositive = sum(tmpScores);
    predictedNegative = numSamples - predictedPositive;
    MCCDenominator = sqrt(conditionPositive*conditionNegative*predictedPositive*predictedNegative);

    tp = sum(tmpScores & labels);
    fn = conditionPositive - tp;
    fp = predictedPositive - tp;
    tn = conditionNegative - fp;

    tmpFpr = fp/conditionNegative;
    tmpMCC = tp*tn - fp*fn;
    if(MCCDenominator ~= 0)
        tmpMCC = tmpMCC / MCCDenominator;
    end

    fpr = [fpr; tmpFpr];
    mcc = [mcc; tmpMCC];

end