The 1995/1996 Household Income, Expenditure - (PDF, 101 mb ...

WTCALC.SAS 1
0ptions pagesize=60;
libname out 'c:\temp';
filename aO 'c:\process\sampchk.txt';
filename a1 'c:\process\sample1.txt';
/*First get all the actual data from validated final */;
/*Get the household values for cumulative age/sex vars */;
/*Persons is the total number of persons in the k-th sample HH */;
/*in the j-th sample chunk in the i-th sample sv in the h-th stratum*/;
/*sampchk.txt came from sorting the coy rectype by gOY urbrur */;
/*km sv, then proc summary var males females by same key */i
/*The actual data has 14 wrong SV or KM codes in the second */i
/*half sample. To attach the weights, these must be corrected */;
/*first. Do this to the raw data first with the program below */i
/**************************************************************
Data sortsamp;
set in.cov;
proc sort; by gOY urbrur km sv; run;
data null;
set sortsarnp;
if gov=1 and urbrur=1 and km=27 and psu=48 then sV=020;
if gov=1 and urbrur=1 and km=34 and psu=60 then sv=030;
if gov=12 and urbrur=O and km=9 and psu=21 then sv=370i
if gov=12 and urbrur=O and km=16 and psu=11 then sv=210;
if gov=12 and urbrur=O and km=16 and psu=12 then sv=390;
if gov=14 and urbrur=O and km=3 and psu=65 then sv=110;
if gov=18 and urbrur=O and km=4 and psu=134 then sv=220;
if gov=18 and urbrur=l and sv=010 and psu=154 then km=12;
if gov=19 and urbrur=O and km=6 and psu=156 then sv=60;
if gov=19 and urbrur=1 and km=l and psu=227 then sv=20;
if gov=23 and urbrur=O and km=3 and psu=186 then sV=20;
if gov=25 and urbrur=1 and km=2 and psu=198 then sv=20;
if gov=28 and urbrur=l and km=6 and psu=216 then sv=20;
if gov=34 and urbrur=O and sv=050 and psu=280 then km=9;
proc summary; var males females; by gOY urbrur km SVi run;
output out=sampchk sum=; run;
*************************************************************/;
DATA first;
infile aO missover;
input gOY 1-8 urbrur 9-17 km 18-26 sv 27-35 psu 36-44
completd 45-53 males 54-62 females 63-71;
persons=males+females;
proc sort;
by gOY urbrur km sv psu;
run;
.,
data second;
/* Weighting of Survey data based on Megill(95) */;
/* urpop is 1986 Census governorate/urban or rural population */;
/* chunks is the number of sample chunks in the stratum */;
/* svpop is the 1986 Census shiakha or village population */;
/* svhhs is the 1986 Census frame # HHs in the shiakha or village */;
infile a1 missover;
input gOY 1-2 urbrur 3-4 km 5-7 sv 9-11 name $ 13-36 psu 37-39
urpop 40-47 chunks 48-50 svpop 51-57 svhhs 58-63;
run;
data one;
set second;
if urbrur=2 then urbrur=o;