workdir = './matfiles/';
load([workdir 'acc_svm_rbf_q.7k.mat']);
load([workdir 'OneVsRest.7k.mat']);

a = permute(mean(tacc,1),[3 2 1]);
b = permute(var(tacc,[],1),[3 2 1]);

aa = reshape(a,numel(a),1);
bb = reshape(b,numel(b),1);

[sa IIa]=sort(aa,'descend');
[Ia Ja]= ind2sub(size(a),IIa);

[sb IIb]=sort(reshape(bb,numel(b),1),'ascend');
[Ib Jb]= ind2sub(size(b),IIb);

% [c t gamma]'
params_by_acc = [params(:,Ia);trbfgamma(Ja);Ia';Ja';sa';bb(IIa)'];
params_by_var = [params(:,Ib);trbfgamma(Jb);Ib';Jb';aa(IIb)';sb'];


fd = fopen('params_svm_rbf_q.txt','w');
	fprintf(fd,'translated RBF\n');
	fprintf(fd,'params by mean accuracy (descending)\n');
	fprintf(fd,'\tc=%f t=%f g=%f mean_acc(%d,%d)=%f var=%f\n', ...
		reshape(params_by_acc(:,1:10),numel(params_by_acc(:,1:10)),1));
		
	fprintf(fd,'params by var (ascending)\n');
	fprintf(fd,'\tc=%f t=%f g=%f mean_acc(%d,%d)=%f var=%f\n', ...
		reshape(params_by_var(:,1:10),numel(params_by_var(:,1:10)),1));


	fprintf(fd,'\n\nstandard RBF\n');
	fprintf(fd,'\tg=%f mean_acc=%f var=%f\n',reshape([rbfgamma;mean(acc,1);var(acc)],3*length(rbfgamma),1));
fclose(fd);

acc_centered = (a-mean(aa))/var(bb);

figure('vis', 'off');
%figure(1);
hold on;
showmat(acc_centered);
title('accuracy (centered)');
saveas(gcf, 'acc_svm_rbf_q', 'png');

figure('vis', 'off');
%figure(2);
showmat((b-mean(bb))/var(bb));
title('variance (centered)');
saveas(gcf, 'var_svm_rbf_q', 'png');