ECG-Kit 1.0

%FEATRANK Feature ranking on individual performance for classification
% 
% 	[I,F] = FEATRANK(A,CRIT,T)
% 	  I   = A*FEATRANK([],CRIT,T)
% 	  I   = A*FEATRANK(CRIT,T)
% 
% INPUT
%   A      input dataset
%   CRIT   string name of a method or untrained mapping, default 'NN'
%   T      validation dataset (optional)
%
% OUTPUT
%   I      vector with sorted feature indices
%   F      vector with criteria values
%
% DESCRIPTION
% Feature ranking based on the training dataset A. CRIT determines 
% the criterion used by the feature evaluation routine feateval. If 
% the dataset T is given, it is used as test set for feateval. In I 
% the features are returned in decreasing performance. In F the 
% corresponding values of feateval are given. Default: crit='NN'.
% 
% SEE ALSO (<a href="http://37steps.com/prtools">PRTools Guide</a>)
% MAPPINGS, DATASETS, FEATEVAL, FEATSELO, FEATSELB, FEATSELF,
% FEATSELP, FEATSELM

% Copyright: R.P.W. Duin, duin@ph.tn.tudelft.nl
% Faculty of Applied Sciences, Delft University of Technology
% P.O. Box 5046, 2600 GA Delft, The Netherlands

function [I,F] = featrank(varargin)

argin = shiftargin(varargin,'char');
argin = setdefaults(argin,[],'NN',[]);
if mapping_task(argin,'definition')
  I = define_mapping(argin,'fixed');
else
  [a,crit,t] = deal(argin{:});
	[m,k,c] = getsize(a);
	F = zeros(1,k);
	isvaldfile(a,1,2); % at least 1 object per class, 2 classes
	a = testdatasize(a);
	iscomdset(a,t);
	
	if isempty(t)
		for j = 1:k
			F(j) = feateval(a(:,j),crit);
		end
	else
		% run the criterion on the validation set
		for j = 1:k
			F(j) = feateval(a(:,j),crit,t(:,j));
		end
	end
	
	[F,I] = sort(-F);
	F = -F;
end
return