
% calculates the win-tie-loss values for gac and gac2 trees


mre31 = mre0;    % the mre values for gac 1
mre32 = mreGac2;

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


% start calculating win tie loss values
for i = 1:1
    for j = 1:2
        
        if i ~= j   % if i and j are different
            
            for k = 1:size(mre31,1)
                eval(['[P,H] = RANKSUM(mre3',num2str(i),'(',num2str(k),',:),mre3',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(mre3',num2str(i),'(',num2str(k),',:));']);
                    eval(['medj = median(mre3',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

