filename ped414 'p414.dat';
options center linesize=78 nodate;
data;
infile ped414;
input @4 y row col;
if row=1 and col=1 then d11=1; else d11=0;
if row=1 and col=2 then d12=1; else d12=0;
if row=2 and col=1 then d21=1; else d21=0;
if row=2 and col=2 then d22=1; else d22=0;
r1=d12-d11;
r2=d11+d21;
c1=d11+d12;
c2=d21-d11;
rc1=d12+d11;
rc2=d11+d21;
o1=(d11+d12)-(d21+d22);
o2=(d11+d21)-(d12+d22);
o3=(d11-d12)-(d21-d22);
proc print;
proc reg;
  model y=d11 d12 d21 d22;
  model y=d11 d12 d21/corrb;
  model y=r1 r2;
  model y=c1 c2;
  model y=rc1 rc2;
  model y=o1 o2 o3/corrb;
  model y=o2 o3;
  model y=o1 o3;
  model y=o1 o2;
proc anova;
  class row col;
  model y=row col row*col;
run;
quit;