LIBNAME A "\\H123456789\Users\User1\Documents\Lesotho DHS 2014";

DATA WORK.LESOTHO;
	SET A.LSKR71FL;

	weight = (V005 / 1000000);
	psu    = V021;
	strata = V023;

/** CHILD'S AGE IN MONTHS **/
	child_age=(V008-B3);

	IF child_age < 24;
	*IF child_age LT 6;

/** CHILD LIVES WITH MOTHER (RESPONDENT) **/
	IF b9 EQ 0;

/* REMOVING SCIENTIFIC NOTATION FROM RESULTS OUTPUT */
PROC TEMPLATE;
	EDIT BASE.FREQ.CROSSTABFREQS;
	EDIT FREQUENCY;
	FORMAT=BEST12.;
	END;
	END;
RUN;

/** FINDING THE YOUNGEST CHILD USING BIDX **/

PROC SORT DATA=WORK.LESOTHO;
	BY V001 V002 V003;
RUN;

PROC SQL;
	TITLE 'MINIMUM BIDX OR YOUNGEST CHILD';
	CREATE TABLE A.bidx_query AS
		SELECT V001,V002,V003,BIDX,M4,M38, MIN(BIDX) AS MINBIDX format=COMMA16.
		FROM WORK.LESOTHO
		GROUP BY V001,V002,V003
		ORDER BY V001,V002,V003;
RUN;

/* MERGING */
PROC SORT DATA=A.bidx_query;
	BY V001 V002 V003;
RUN;

PROC SORT DATA=WORK.LESOTHO;
	BY V001 V002 V003;
RUN;

/** MERGING/COMBINING THE DATA FILES **/

DATA A.LESOTHO14_MERGED;
	MERGE WORK.LESOTHO (IN=x) A.bidx_query (IN=y);
	BY V001 V002 V003;
RUN;

DATA WORK.LESOTHO14_MERGED;
	MERGE WORK.LESOTHO (IN=x) A.bidx_query (IN=y);
	BY V001 V002 V003;
RUN;

/*********************************************************************************/

/** START HERE **/

DATA WORK.LESOTHO14_ExclBF;
	SET WORK.LESOTHO14_MERGED;

	IF (BIDX <= MINBIDX);

	FORMAT diet$29. feeding$26. exclusive_feed$6. none_breast$6. child_age_grp$5. age_at_birth$8. birth_order$3.;

/** CREATE FEEDING VARIABLES **/
	water=0;
	liquids=0;
	milk=0;
	solids=0;
	breast=0;
	bottle=0;

/** TO DETERMINE IF CHILD IS GIVEN WATER, SUGAR WATER, JUICE, TEA OR OTHER **/
	IF V409 IN(1:7)		THEN	water=1;

/** IF GIVEN OTHER LIQUIDS **/
	Array liquids2 {9} $ v409a v410 v410a v412c v413 v413a v413b v413c v413d;  

	DO a=1 to 9;
		IF 1<=liquids2{a}<=7 THEN liquids=1;
	End;

/** IF GIVEN POWDER/TINNED milk, FORMULA OR FRESH milk **/
	Array milk2 {4} $ v411 v411a v412 v414p;  

	DO b=1 to 4;
		IF 1<=milk2{b}<=7 THEN milk=1;
	End;

/** IF STILL BREASTFEEDING **/
	IF M4 IN(95)		THEN	breast=1;

/** IF WAS EVER BOTTLE FED **/
	IF M38 IN(1)		THEN	bottle=1;

/** IF GIVEN ANY SOLID FOOD **/
	Array solids2 {26} $ V414A V414B V414C V414D V414E V414F V414G V414H V414I V414J V414K V414L V414M V414N V414O V414P V414Q V414R V414S V414T V414U V414V V414W V414X V414Y V414Z;  

	DO c=1 to 26;
		IF 1<=solids2{c}<=7 THEN solids=1;
		IF v412a EQ 1 OR v412b EQ 1 THEN solids=1;
	End;

/** FIRST INDICATOR VARIABLE **/
diet=7;
	IF water=0 & liquids=0 & milk=0 & solids=0	THEN diet = "0. given only water";
	IF water=1 & liquids=0 & milk=0 & solids=0	THEN diet = "1. given only liquids";
	IF           liquids=1 & milk=0 & solids=0 	THEN diet = "2. given only milks";
	IF                       milk=1 & solids=0 	THEN diet = "3. given only solids";
	IF                       milk=0 & solids=1 	THEN diet = "4. given only milk and solids";
	IF                       milk=1 & solids=1 	THEN diet = "5. not still breastfeeding";
	IF 								breast=0	THEN diet = "6. not now being breastfed";

/** EBF **/
ebf=0;
	IF diet = "0. given only water" THEN ebf = 1;

/** SECOND INDICATOR VARIABLE (FEEDING) **/
feeds=5;
	IF diet="6. not now being breastfed"	THEN feeds = 0;
	IF diet="0. given only water"			THEN feeds = 1;
	IF diet="1. given only liquids"			THEN feeds = 2;
	IF diet="2. given only milks"			THEN feeds = 3;
	IF diet="3. given only solids"			THEN feeds = 4;

	IF feeds=0	THEN feeding = "0. Not breastfeeding";
	IF feeds=1	THEN feeding = "1. Exclusive breastfeeding";
	IF feeds=2	THEN feeding = "2. +Water";
	IF feeds=3	THEN feeding = "3. +Liquids";
	IF feeds=4	THEN feeding = "4. +Other Milk";
	IF feeds=5	THEN feeding = "5. +Solids";
	LABEL  feeding  ="Breastfeeding status";

/** EXCLUSIVE BREASTFEEDING **/
exclusive_feed="0. No";
	IF feeds IN(1)	THEN	exclusive_feed = "1. Yes";	
	LABEL  exclusive_feed  ="Exclusive breastfeeding";

/** NOT BREASTFEEDING **/
none_breast="0. No";
	IF feeds IN(0)	THEN	none_breast = "1. Yes";	
	LABEL  none_breast  ="Not breastfeeding";

/** CREATING THE PREDOMINANTLY BREASTFEEDING VARIABLE **/
predom=0;
	IF feeds IN(1:3)	THEN	predom = 1;	

/** CHILD AGE GROUPS **/
	IF child_age in(0:1)		THEN child_age_grp = "00-01";
	IF child_age in(2:3)		THEN child_age_grp = "02-03";
	IF child_age in(4:5)		THEN child_age_grp = "04-05";
	IF child_age in(6:8)		THEN child_age_grp = "06-08";
	IF child_age in(9:11)		THEN child_age_grp = "09-11";
	IF child_age in(12:17)		THEN child_age_grp = "12-17";
	IF child_age in(18:23)		THEN child_age_grp = "18-23";

/** MATERNAL AGE AT BIRTH **/
	agebirth = int((B3-V011)/12);

	IF agebirth <20				THEN	age_at_birth ="1. <20";
	IF agebirth IN(20:34)		THEN	age_at_birth ="2. 20-34";
	IF agebirth >=35			THEN	age_at_birth ="3. 35-49";

/** BIRTH ORDER **/
	birth_order1 = bord;
	
	IF B0 IN(2) 				THEN	birth_order1 = bord - 1;
	IF B0 IN(3) 				THEN	birth_order1 = bord - 2;

	IF birth_order1 IN(1)		THEN	birth_order ="1";
	IF birth_order1 IN(2:3)		THEN	birth_order ="2-3";
	IF birth_order1 IN(4:5)		THEN	birth_order ="4-5";
	IF birth_order1 >=6			THEN	birth_order ="6+";
	
	DROP a b c feeds agebirth birth_order1;
RUN;

PROC FREQ DATA=WORK.LESOTHO14_EXCLBF;
	WEIGHT	WEIGHT;
	TITLE "Table 11.3 Breastfeeding status by age";
	TABLE child_age_grp * feeding 			/NOCOL NOFREQ NOPERCENT;
RUN;

/*
PROC SURVEYFREQ DATA=WORK.LESOTHO14_EXCLBF;
	STRATA strata;
	CLUSTER psu;
	WEIGHT WEIGHT;
	TITLE "Table 11.3 Breastfeeding status by age";
	TABLE child_age_grp * feeding 			/NOSTD ROW CL ALPHA=0.05;
RUN;
*/

PROC FREQ DATA=WORK.LESOTHO14_EXCLBF;
	WEIGHT	WEIGHT;
	TITLE "EXCLUSIVE BREASTFEEDING";
	TABLE child_age_grp * exclusive_feed	/NOCOL NOFREQ NOPERCENT;
RUN;

/********************************************************************/

PROC FREQ DATA=WORK.LESOTHO14_EXCLBF;
	WEIGHT	WEIGHT;
	TITLE "TABLE 11.3 - N";
	TABLE 	child_age_grp * feeding			/NOCOL NOROW NOPERCENT;
	TABLE age_at_birth * feeding			/NOCOL NOROW NOPERCENT;
	TABLE birth_order * feeding				/NOCOL NOROW NOPERCENT;

	TABLE 	child_age_grp * exclusive_feed	/NOCOL NOROW NOPERCENT;
	TABLE age_at_birth * exclusive_feed		/NOCOL NOROW NOPERCENT;
	TABLE birth_order * exclusive_feed		/NOCOL NOROW NOPERCENT;
	TABLE 	V190 * exclusive_feed			/NOCOL NOROW NOPERCENT;
RUN;

/********************************************************************/

PROC SURVEYLOGISTIC DATA=WORK.LESOTHO14_EXCLBF;
	STRATA strata;
	CLUSTER psu;
	WEIGHT weight;
	TITLE "REGRESSION: EXCLUSIVE BREASTFEEDING";
	CLASS B4(REF=FIRST) age_at_birth(REF=FIRST) birth_order(REF=FIRST) V190(REF=FIRST) V025(REF=LAST) V024(REF=LAST) /PARAM=REF;
	MODEL exclusive_feed (EVENT = '1. Yes') = B4 age_at_birth birth_order V190 V025 V024;
	ODS OUTPUT OddsRatios=OR2 ;
RUN;

/* GOOD LUCK */