%mre21 = mre0;
mre21 = mreGac2;
mre22 = mre1;
mre23 = mre3;
mre24 = mre5;
mre25 = mre7;
mre26 = mre9;
mre27 = mrex;

win = zeros(7,1);
tie = zeros(7,1);
loss = zeros(7,1);

% start calculating win tie loss values
for i = 1:7
    
    for j = 1:7
        
        if i ~= j   % if i and j are different
            
            for k = 1:size(mre1,1)
                eval(['[P,H] = RANKSUM(mre2',num2str(i),'(',num2str(k),',:),mre2',num2str(j),'(',num2str(k),',:));']);
                if H == 0   % if they are the same
                    eval(['tie(',num2str(i),') = tie(',num2str(i),') +1;']);
                    eval(['tie(',num2str(j),') = tie(',num2str(j),') +1;']);
                else
                    % calculate median values
                    eval(['medi = median(mre2',num2str(i),'(',num2str(k),',:));']);
                    eval(['medj = median(mre2',num2str(j),'(',num2str(k),',:));']);
                    if medi < medj
                        eval(['win(',num2str(i),') = win(',num2str(i),') +1;']);
                        eval(['loss(',num2str(j),') = loss(',num2str(j),') +1;']);
                    else
                        eval(['win(',num2str(j),') = win(',num2str(j),') +1;']);
                        eval(['loss(',num2str(i),') = loss(',num2str(i),') +1;']);
                    end
                end
            end
            
        end
    end
end

win
tie
loss
% 
% % % plot the win loss tie values 
% % 
%     hold on;
%     
%     title('Win/Loss/Tie values for COC81');
%     
%     subplot(3,2,1);
%     bar([1:8],win);
%     xlabel('gac1,gac2,k=1,k=2,k=4,k=8,k=16,k=''best k'' respectively');
%     ylabel('Win Value');
%     
%     subplot(3,2,2);
%     bar([1:8],loss);
%     xlabel('gac1,gac2,k=1,k=2,k=4,k=8,k=16,k=''best k'' respectively');
%     ylabel('Loss Value');
%     
%     subplot(3,2,3);
%     bar([[1:8]],tie);
%     xlabel('gac1,gac2,k=1,k=2,k=4,k=8,k=16,k=''best k'' respectively');
%     ylabel('Tie Value');
%     
%     hold off;
% 
% 
