data data2;
seed = 34191;
phi = 0.5;
do site = 1 to 3;
do trt = 0 to 1;
do subj = 1 to 6;
u = ranuni(54469);
output;
end;
end;
end;
drop subj;
run;
proc sort data=data2;
by site u;
run;
data data2;
set data2;
by site;
retain subj;
if first.site then subj=0;
else subj=subj+1;
subjid = site*100 + subj;
format trt trtFmt. site siteFmt.;
run;
18
proc print data=data2 split='*';
pageby site;
by site;
id subjid;
var trt;
label subjid='Subject*Number*-------'
trt='Treatment*Assignment*----------';
run;
proc freq data=data2;
tables site*trt / nopercent nocol norow;
run;