%% fcurves
% Feature curves for Bayes-Normal, on the Satellite dataset.
% Average over 25 repetitions
%
% <http://37steps.com/software/prtools/ PRTools> and
% <http://37steps.com/software/prdatasets/ PRDataSets> should be in the path.
% See http://37steps.com/prtools for more. 
% Download m-file from <http://37steps.com/exam/fcurves/fcurves.m Download the m-file from here>

%%        Initialization
delfigs 
a = satellite;
a = setprior(a,getprior(a));
w = featself(a,'maha-s');
trainsize = [20 30 50 100 500];
iter = 25;

%% Create feature curves for optimized feature ranking
x = a*w;
randreset;
e1 = cell(1,numel(trainsize));
% prwaitbar calls are just used for reporting progress and may be skipped
[n,s,count] = prwaitbarinit('Processing %i curves:',2*numel(trainsize));
for j=1:numel(trainsize)
  e1{j} = clevalf(x,remclass(2)*qdc([],[],1e-6),[1:15],trainsize(j),iter,[],testd);
  count = prwaitbarnext(n,s,count);
end
plote(e1,'nolegend')
legend('train size: 20','train size: 30','train size: 50','train size: 100','train size: 500')
title('Feature curve for Satellite dataset, optimized feature order')
fontsize(14)
set(gcf,'position',[ 680 558 808 420]);

%% Create feature curves for randomized feature ranking
figure
randreset;
x = a(:,randperm(size(a,2)));
e2 = cell(1,numel(trainsize));
for j=1:numel(trainsize)
  randreset;
  e2{j} = clevalf(x,remclass(2)*qdc([],[],1e-6),[1:15],trainsize(j),iter,[],testd);
  count = prwaitbarnext(n,s,count);
end
plote(e2,'nolegend');
legend('train size: 20','train size: 30','train size: 50','train size: 100','train size: 500')
title('Feature curve for Satellite dataset, random feature order')
fontsize(14)
set(gcf,'position',[ 680 558 808 420]);