function mcdplot
global STATS

bool=STATS(:,7);
ig=find(bool);
ib=find(~bool);

a=STATS(:,1);
b=STATS(:,2);
phi00=STATS(:,3);
phi01=STATS(:,4);
phi10=STATS(:,5);
phi11=STATS(:,6);

da = abs(phi00-phi01);
db = abs(phi10-phi11);
Dab0 = abs(phi00-phi10);
Dab1 = abs(phi01-phi11);


keyboard

return