用这个:
%的代码
函数(auc,玻璃钢,tpr) = fastAUC(标签、分数、plot_flag)
%函数(auc,玻璃钢,tpr) = myauc(标签、分数、plot_flag)
%
%这个函数计算m AUC值m排名列表。
% n是排名项的数量。
% m是许多不同的排名。
%
%输入:标签是nXm二进制逻辑。
%的分数是nXm真实。对于高AUC更高的分数
% = = 1标签。
% plot_flag:二进制国旗,如果真的,那么m将绘制ROC曲线
%(默认错误)。
%
%输出:auc是真正mx₁, ROC曲线下的面积。
%玻璃钢nXm真实,假阳性。
% tpr nXm是真实的,真正的积极的利率。
如果~ (“plot_flag”,“var”)
plot_flag = 0;
结束
如果~ islogical(标签)
错误(“标签输入应该逻辑”);
结束
如果~ isequal(大小(标签),大小(分数))
错误(标签和分数应该有相同的大小);
结束
[n m] =大小(标签);
num_pos =总和(标签);
如果任何(num_pos = = 0)
错误(“没有积极的标签进入”);
结束
如果任何(num_pos = = n)
错误(“没有负面的标签进入”);
结束
[~,scores_si] =排序(分数,“下”);
清晰的分数
scores_si_reindex = scores_si +的(n - 1) * (0: m - 1) * n;
l =标签(scores_si_reindex);
清晰的scores_si标签
tp = cumsum (l = = 1,1);
fp = repmat ((1: n)”, [1 m]) tp;
num_neg = n-num_pos;
玻璃钢= bsxfun (fp, @rdivide num_neg);%假阳性率
tpr = bsxfun (tp, @rdivide num_pos);%真阳性率
%绘制ROC曲线
如果plot_flag = = 1
情节(玻璃钢,tpr);
包含(“假阳性”);
ylabel (“真阳性”);
结束
auc = (tpr。* ((diff (fp) = = 1);0(1米)])。/ num_neg;
