source: prextra/kannc.m @ 7

%KANNC 1-Nearest Neighbor Classifier using ANN Matlab Wrapper
%
%   W = KANNC(A,,K,OPTION)
%
% INPUT
%   A      Dataset
%   K      Number of desired neighbours
%   OPTION Options for ANNQUERY
%
% OUTPUT
%   W      NN classifier using the ANN Query package
%
% DESCRIPTION  
% This is the nearest neighbor implementation for PRTools using the ANN
% Matlab Wrapper package. It should be in the path. If needed download it
% http://webscripts.softpedia.com/scriptDownload/ANN-MATLAB-Wrapper-Download-33976.html
%
% SEE ALSO
% MAPPINGS, DATASETS, ANNQUERY, KNNC

% Copyright: R.P.W. Duin, r.p.w.duin@prtools.org
% Faculty EWI, Delft University of Technology
% P.O. Box 5031, 2600 GA Delft, The Netherlands


function out = kannc(a,knn,opt)

        prtrace(mfilename);
        
        name = 'KANNC';
        % No input data, return an untrained classifier.
        if nargin < 3, opt = []; end
        if nargin < 2, knn = 1; end
        if (nargin == 0) | (isempty(a))
                out = mapping(mfilename,'untrained',{knn,opt});
                out = setname(out,name);
        elseif isdataset(a) & ~ismapping(knn) % training
                islabtype(a,'crisp');
                isvaldfile(a,1,2); % at least 1 object per class, 2 classes
                a = testdatasize(a);
                a = testdatasize(a,'objects');
                a = seldat(a);    % get labeled objects only
                [m,k,c] = getsize(a);
                v.data = a;
                v.knn  = knn;
                v.opt  = opt;
                out = mapping(mfilename,'trained',v,getlablist(a),k,c);
                out = setname(out,name);
        elseif nargin == 2 & ismapping(knn)  % execution
                % to avoid confusion, rename trained mapping
                w = knn;
                v = +w;   % get datafield
                [k,c] = size(w);
                nlab = getnlab(v.data);
                J = annquery(+v.data',(+a)',v.knn)';
                n = size(J,1); % no of test objects
                J = nlab(J);
                if size(J,2) == 1
                        out = zeros(2,n);
                        out([1:2:n*2]+J'-1) = ones(1,n);
                        out = out';
                else
                        out = hist(J',[1:c])';
                end
                % Use Bayes estimators for posteriors
                out = (out+ones(size(out,1),c))/(v.knn+c);
                out = setdat(a,out,w);
        else
                error('Illegal input')
        end