OPTIONS PS=54 LS=64 NODATE NONUMBER;

DATA skinfold; INPUT y @@; CARDS;
12 33 45 21 34 55 18 23 22 27
62 31 49 50 40 33 10 26 17 25
;
******************************;
*** BOOTSTRAPPING THE MEAN ***;
******************************;
DATA boot;
     array x(20) x1-x20;
     n = 20;          * n = # of data points;
     bssamp = 100;    * bssamp= # of boostrap samples;

     DO j = 1 to bssamp;
     DO i = 1 to n;
        readit = CEIL(RANUNI(0)*n);
        SET skinfold POINT=readit NOBS=n;
        x(i) = y;
     END;
        bootmn= MEAN(OF x1-x20); output;
     END;
     STOP;
KEEP bootmn;

PROC SORT DATA=boot; BY bootmn;
PROC PRINT DATA=boot;
TITLE '100 BOOTSTRAP MEANS';

PROC MEANS DATA=boot; VAR bootmn;

RUN;