data ldl;
 input subject baseline follow_up LDL_decrease;
cards;
  1  160.5  168.1  -7.6
  2  195.3  181.4  13.9
  3  181.7  154.6  27.1
  4  175.1  160.3  14.8
  5  198.3  192.0   6.3
  6  215.5  173.5  42.0
  7  227.9  186.2  41.7
  8  201.7  183.2  18.5
  9  161.5  130.3  31.2
 10  189.0  165.0  24.0
;
run;
proc univariate data=ldl;
 var LDL_decrease;
run;
data ldl_sim;
 label y='LDL Decrease';
 do sample=1 to 10000;
  do obs=1 to 10;
   y=21.19+15.45*rannor(3132017);
   output;
  end;
 end;
run;
proc summary data=ldl_sim;
 var y;
 by sample;
 output out=ldl_means mean=ybar;
run;
data ldl_means;
 set ldl_means;
 label ybar='Mean LDL Decrease';
run;
proc univariate data=ldl_means;
	ods select histogram;
	var ybar;
	histogram ybar/vscale=count endpoints=0 to 44 by 2;
run;
proc univariate data=ldl_sim;
	ods select histogram;
	var y;
	histogram y/vscale=count endpoints=-48 to 92 by 2;
run;
proc univariate data=ldl_means;
	ods select histogram;
	var ybar;
	histogram ybar/vscale=count endpoints=-48 to 92 by 2;
run;
* Standardized hist;
data sldl_sim;
 label y='Standardized Mean LDL Decrease';
 do sample=1 to 10000;
  do obs=1 to 10;
   y=rannor(3132017);
   output;
  end;
 end;
run;
proc univariate data=sldl_sim;
	ods select histogram;
	var y;
	histogram y/vscale=count endpoints=-4 to 4 by .25;
run;