* copy, paste, and run alloc macro first

data stratinfo;
input group $ popsize s2 c;
cards;
TownA 155 25  9
 
TownB 62  225 9
 
Rural 93  100 16;

* optimal;


%alloc(type = optimal, strata=group, setup=stratinfo,param=mean, npop=popsize, var = s2, cost = c, bound = 2); 


* proportional allocation ;
%alloc(type = proportional, strata=group, setup=stratinfo,
          param=mean, npop=popsize, var = s2, cost = c, bound = 2);

* Neyman allocation ;
%alloc(type = Neyman, strata=group, setup=stratinfo,param=mean, npop=popsize, var = s2, cost = c, bound = 2);

* solving for a fixed cost, calculations on page 136 ;
%alloc(type = optimal, opt = fixedc, strata=group, setup=stratinfo,
param=mean, npop=popsize, var = s2, cost = c, maxcost = 500) ;