options nocenter ps=1000;
%let nsims=20;
%let n=3;
%let sigma=2;
%let seed=123;
data fake;
   do dataset=1 to &nsims;
   do trt=1 to 10;
      do rep=1 to &n;
	  y=10 + rannor(&seed)*σ
	  output;
	  end;
   end;
   end;
run;

proc mixed;
   by dataset;
   class trt;
   model y=trt;
   lsmeans trt/diff;
   ods output diffs=diffs;
run;
data diffs;
   set diffs;
   sig=probt<.05;
run;
proc freq;
   tables dataset*sig;
run;