程序如下：

%macro BA_ATE_EP9 ( case =, subj =, x =, y =, sp =, rp =, beta =, alpha =, sim = );

       %if &case=1 %then %do;

       %do icc=5 %to 9 %by 1;

          %do sim=1 %to ∼

          data SML_BH (keep=icc subj rep xij yij);
             icc=icc;
             s1=sqrt(icc/10);
             s2=sqrt(1-icc/10);
              do subj=1 to &sp;
                sbb=0+s1*rannor(0);
                  do rep=1 to &rp;
                    sww=0+s2*rannor(0);
                    dij=sbb+sww;
                    xij=10+rannor(0);
                    yij=xij+dij;
                  output;
                 end;
              end;
          run;

        %end;
        %end;
        %end;

%mend  BA_ATE_EP9;


%BA_ATE_EP9 ( case =1, subj =subj, x =xij, y =yij, sp =2, rp =2, beta =0.95, alpha =0.05, sim =1 );



谢谢！