function mcdaux
global STATS
[M,N] = size(STATS);
II = list_all_ind([2 2 2 2 2 2]) - 1;
for i=2:N-1
  for j=i+1:N-1
    %for k=j+1:N-1
      si = STATS(:,i);
      sj = STATS(:,j);
     % sk = STATS(:,k);
      %if all(si)
      %'SUCEEEEEEEEESSSSSSSSSSS!!!!!!!!'
      %end
      
      %if all(si+sj+sk)
      if all(si+sj)
	%'success!'
	%[i j]
	%fprintf('success! i= %d j=%d  k=%d \n',i,j,k);
	%fprintf('success! i= %d j=%d  \n',i,j);
	prti1 = num2str(find(II(i,:)));
	prti0 = num2str(find(~II(i,:)));
	prtj1 = num2str(find(II(j,:)));
	prtj0 = num2str(find(~II(j,:)));
	fprintf('success! i=%d: [%s] [%s] OR j=%d: [%s] [%s]  \n',i,prti1,prti0,j,prtj1,prtj0);
      end
    %end
  end
end


return