filename ped349 'p349.dat';
options center linesize=78 pagesize=65 nodate;
data;
infile ped349;
input id y grp;
if grp=1 then d1=1; else d1=0;
if grp=2 then d2=1; else d2=0;
if grp=3 then d3=1; else d3=0;
o1=d1-d2;
o2=(d1+d2)-2*d3;
o3=(d1-d2)/2;
o4=((d1+d2)/2)-d3;
proc reg;
model y=d1 d2 d3;
model y=d1 d2/covb;
model y=d2 d3;
model y=d1 d3;
model y=d1 d2 d3/noint covb;
model y=o1 o2/covb;
model y=o3 o4/covb;
run;
quit;