/** SKILLED BIRTH ATTENDANT **/
DATA SLDHS2013_1;
	SET WORK.SLDHS_2013;
	WHERE m3a in(0,1,9) AND m3b in(0,1,9) AND m3c in(0,1,9) AND v208 in(1,2,3,4);
	/* THE CODE ABOVE SELECTS THE SAMPLE TO MATCH THE TABLE IN THE DHS REPORT */

	/* NOW, TO DERIVE THE DEPENDENT VARIABLE (Skilled_birth)*/

		IF	m3a = 1	or	m3b = 1	or	m3c = 1	THEN	skilled = 1; else
		IF	m3a 	or m3b 		or m3c NE 1 THEN 	skilled = 0;

		FORMAT skilled_birth$32.;
			IF skilled =	0	THEN 	Skilled_birth = "0. Non-skilled provider";
			IF skilled =	1	THEN 	Skilled_birth = "1. Skilled provider";

			WEIGHT = (v005 / 1000000);
	PROC TEMPLATE;
	EDIT BASE.FREQ.CROSSTABFREQS;
	EDIT FREQUENCY;
	FORMAT=BEST12.;
	END;
	END;
RUN;

/* CHECK IF IT MATCHES THE RESULTS ON THE DHS REPORT*/

/* Percentage delivered by a skilled provider in Sierra Leone DHS 2013 report = 59.7% */

PROC FREQ DATA=WORK.SLDHS2013_1;
	TABLE Skilled_birth /NOCOL NOROW;
	WEIGHT WEIGHT;
RUN;

/* CHECK WITH OTHER VARIABLES */

PROC FREQ DATA=WORK.SLDHS2013_1;
	TABLE V025 * Skilled_birth /NOCOL NOFREQ NOPERCENT;
	TABLE V024 * Skilled_birth /NOCOL NOFREQ NOPERCENT;
	TABLE V190 * Skilled_birth /NOCOL NOFREQ NOPERCENT;
	WEIGHT WEIGHT;
RUN;