compute_thres.m

关于人脸识别的一个VC程序

function [thr_add, thr_add_old] = compute_thres(resp, labels)

[resp, ind] = sort(resp);
labels      = labels(ind);

%now look for the point which cuts the label vector into 'pure' parts
num = size(labels,2);

err=zeros(num-1,1);

for cut=1:(num-1),
  err(cut) = (cut + sum(labels(1:cut)))/2 + ((num-cut) - sum(labels((cut+1):num)))/2;
end

best_ind = find(min(err) == err);

%plot(err);
thr_add = -mean((resp(best_ind) + resp(best_ind+1))/2);

best_ind = floor(mean(best_ind));

% choose thresh to the right!

thr_add_old = -(resp(best_ind)+resp(best_ind+1))/2;