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;
   plot='LDL_decrease      ';
   z=21.19+15.45*rannor(3132017);
   y=z;
   output;
  end;
 end;
run;
*proc contents;run;
data ldl_sim1;
 set ldl_sim;
 drop y;
run;
data ldl_sim;
 set ldl_sim;
 keep plot y;
run;
proc summary data=ldl_sim1;
 var z;
 by sample;
 output out=ldl_means mean=y;
run;
data ldl_means;
 set ldl_means;
 plot='LDL_decrease_means';
 keep plot y;
run;
*proc contents;run;
data ldl_all;
 set ldl_sim ldl_means;
run;
proc sgpanel data=ldl_all;
 panelby plot/rows=2 columns=1;
 histogram y;
run;