DM 'LOG;CLEAR;OUT;CLEAR';
OPTIONS NODATE NONUMBER NOCENTER PS=54 LS=70;

*******************************************************;
*** Systematic sample analysis for equal-sized PSUs ***;
*******************************************************;

DATA systmtc2;
  M0 = 400;    * number of secondary sampling units (SSUs) in population;
  n  = 6;      * number of primary sampling units (PSUs) sampled;
  m  = 8;      * number of SSUs in a PSU;

  wgt = M0/(n*m);
  DO psu = 1 to n;
  DO ssu = 1 to m;
     INPUT y @@; OUTPUT;
  END; END;
DATALINES;
20 19 26 31 21 31 40 40 16 26 22 28 28 31 38 37 
25 26 34 31 26 33 40 43 24 30 29 35 29 36 38 41 
18 22 28 34 29 33 46 45 26 27 32 44 37 45 46 46
;

*** TOTAL = number of PSUs in the population ***;

PROC SURVEYMEANS DATA=systmtc2 TOTAL=50 MEAN CLM SUM CLSUM;
     VAR y;
     CLUSTER psu;
     WEIGHT wgt;
TITLE 'Systematic Sample from Figure 9b';
RUN;