******************************************************************************** /* Citation: Oxford Poverty and Human Development Initiative (OPHI), University of Oxford. 2019 Global Multidimensional Poverty Index - Syria PAPFAM 2009 [STATA do-file]. Available from OPHI website: http://ophi.org.uk/ For further queries, contact: ophi@qeh.ox.ac.uk */ ******************************************************************************** clear all set more off set maxvar 10000 set mem 500m *** Working Folder Path *** global path_in "../rdta/Syria PAPFAM 2009" global path_out "cdta" global path_ado "ado" ******************************************************************************** *** SYRIA PAPFAM 2009 *** ******************************************************************************** ******************************************************************************** *** Step 1: Data preparation *** Selecting variables from KR, BR, IR, & MR recode & merging with PR recode ******************************************************************************** ******************************************************************************** *** Step 1.1 CHILDREN's RECODE (0-60 months) ******************************************************************************** use "$path_in/Syria 2009 - HR.dta", clear rename _all, lower desc agemonth keep if agemonth>=0 & agemonth<=60 count /*NOTE: The data indicate a total of 15,941 children aged 0-5 years. In terms of age in months, these children are between 0-72 months old. We compute underweight and stunting measures for children up to 60 months. We compute BMI-for-age for children from 61-72 months old. This is done in the next section, that is, Step 1.1b The final sample count of children aged 0-60 months that is included in the Global MPI estimation for Syria PAPFAM 2009 is 16,631 children*/ *** Generate individual unique key variable required for data merging *** cluster=cluster number; *** hhnum=household number; *** h108b=line number of eligible child gen double ind_id = cluster*1000000 + hhnum*100 + h108b format ind_id %20.0g label var ind_id "Individual ID" duplicates report ind_id //NOTE: No duplicate observations gen child_CH=1 //Generate identification variable for observations in child recode count if h102==1 /*NOTE: In the context of Syria PAPFAM 2009, all children aged 0-60 months are permenant residents of their HH*/ *** Check the variables that WHO ado needs to calculate the z-scores: *** sex, age, weight, height, measurement, oedema & child sampling weight *** Variable: SEX *** lookfor sex tab h103,miss codebook h103,tab(30) tab h103, nol clonevar gender = h103 desc gender tab gender *** Variable: AGE *** tab agemonth, miss //Child's age is measured in months clonevar age_months = agemonth gen str6 ageunit = "months" lab var ageunit "Months" *** Variable: BODY WEIGHT (KILOGRAMS) *** tab h604, miss clonevar weight = h604 tab h607 h604 if h604==., miss //h607: Result of child measurement desc weight sum weight *** Variable: HEIGHT (CENTIMETERS) tab h605, miss clonevar height = h605 tab h607 h605 if h605==., miss desc height sum height *** Variable: MEASURED STANDING/LYING DOWN *** codebook h606, tab (10) gen measure = "l" if h606==1 //Child measured lying down replace measure = "h" if h606==2 //Child measured standing up replace measure = " " if h606==9 | h606==0 | h606==. //Replace with " " if unknown desc measure tab measure *** Variable: OEDEMA *** lookfor oedema gen oedema = "n" //It assumes no-one has oedema desc oedema tab oedema *** Variable: SAMPLING WEIGHT *** gen sw = hhweight //For household sample weight desc sw sum sw /* For this part of the do-file we use the WHO Anthro and macros. This is to calculate the z-scores of children under 5. Source of ado file: http://www.who.int/childgrowth/software/en/ */ *** Next, indicate to STATA where the igrowup_restricted.ado file is stored: adopath + "$path_ado/igrowup_stata" *** We will now proceed to create three nutritional variables: *** weight-for-age (underweight), *** weight-for-height (wasting) *** height-for-age (stunting) /* We use 'reflib' to specify the package directory where the .dta files containing the WHO Child Growth Standards are stored. Note that we use strX to specify the length of the path in string. If the path is long, you may specify str55 or more, so it will run. */ gen str100 reflib = "$path_ado/igrowup_stata" lab var reflib "Directory of reference tables" /* We use datalib to specify the working directory where the input STATA dataset containing the anthropometric measurement is stored. */ gen str100 datalib = "$path_out" lab var datalib "Directory for datafiles" /* We use datalab to specify the name that will prefix the output files that will be produced from using this ado file (datalab_z_r_rc and datalab_prev_rc)*/ gen str30 datalab = "children_nutri_syr" lab var datalab "Working file" /*We now run the command to calculate the z-scores with the adofile */ igrowup_restricted reflib datalib datalab gender age_months ageunit weight /// height measure oedema sw /*We now turn to using the dta file that was created and that contains the calculated z-scores */ use "$path_out/children_nutri_syr_z_rc.dta", clear *** Standard MPI indicator *** //Takes value 1 if the child is under 2 stdev below the median & 0 otherwise gen underweight = (_zwei < -2.0) replace underweight = . if _zwei == . | _fwei==1 lab var underweight "Child is undernourished (weight-for-age) 2sd - WHO" tab underweight, miss gen stunting = (_zlen < -2.0) replace stunting = . if _zlen == . | _flen==1 lab var stunting "Child is stunted (length/height-for-age) 2sd - WHO" tab stunting, miss gen wasting = (_zwfl < - 2.0) replace wasting = . if _zwfl == . | _fwfl == 1 lab var wasting "Child is wasted (weight-for-length/height) 2sd - WHO" tab wasting, miss count if _fwei==1 | _flen==1 /*Note: In the context of Syria PAPFAM 2009, 448 children are replaced as '.' because they have extreme z-scores which are biologically implausible */ *** Destitution indicator *** //Takes value 1 if the child is under 3 stdev below the median & 0 otherwise gen underweight_u = (_zwei < -3.0) replace underweight_u = . if _zwei == . | _fwei==1 lab var underweight_u "Child is undernourished (weight-for-age) 3sd - WHO" gen stunting_u = (_zlen < -3.0) replace stunting_u = . if _zlen == . | _flen==1 lab var stunting_u "Child is stunted (length/height-for-age) 3sd - WHO" gen wasting_u = (_zwfl < - 3.0) replace wasting_u = . if _zwfl == . | _fwfl == 1 lab var wasting_u "Child is wasted (weight-for-length/height) 3sd - WHO" //Retain relevant variables: keep ind_id child_CH underweight* stunting* wasting* order ind_id child_CH underweight* stunting* wasting* sort ind_id save "$path_out/SYR09_CH.dta", replace //Erase files from folder: erase "$path_out/children_nutri_syr_z_rc.xls" erase "$path_out/children_nutri_syr_prev_rc.xls" erase "$path_out/children_nutri_syr_z_rc.dta" ******************************************************************************** *** Step 1.1b CHILDREN's RECODE (61-72 months) ******************************************************************************** use "$path_in/Syria 2009 - HR.dta", clear rename _all, lower keep if agemonth>=61 & agemonth<=72 count /*NOTE: The final sample count of children aged 61-72 months that is included in the Global MPI estimation for Syria PAPFAM 2009 is 3,010 children*/ *** Generate individual unique key variable required for data merging *** cluster=cluster number; *** hhnum=household number; *** h108b=line number of eligible child gen double ind_id = cluster*1000000 + hhnum*100 + h108b format ind_id %20.0g label var ind_id "Individual ID" duplicates report ind_id duplicates tag ind_id, gen(duplicates) tab h108b if duplicates!=0 bys ind_id: gen line = (_n) replace ind_id = cluster*1000000 + hhnum*190 + line if duplicate!=0 //We assume consecutive hh line starting at 190 duplicates report ind_id //No duplicates at this stage gen child_CH=1 //Generate identification variable for observations in child recode count if h102==1 /*NOTE: In the context of Libya PAPFAM 2014, all children aged 61-72 months are permenant residents of their HH*/ *** Check the variables that WHO ado needs to calculate the z-scores: *** sex, age, weight, height, measurement, oedema & child sampling weight *** Variable: SEX *** lookfor sex tab h103,miss codebook h103,tab(30) tab h103, nol clonevar gender = h103 desc gender tab gender *** Variable: AGE *** tab agemonth, miss //Child's age is measured in months clonevar age_months = agemonth gen str6 ageunit = "months" lab var ageunit "Months" *** Variable: BODY WEIGHT (KILOGRAMS) *** tab h604, miss clonevar weight = h604 tab h607 h604 if h604==., miss //h607: Result of child measurement desc weight sum weight *** Variable: HEIGHT (CENTIMETERS) tab h605, miss clonevar height = h605 tab h607 h605 if h605==., miss desc height sum height *** Variable: Oedema *** lookfor oedema gen oedema = "n" //It assumes no-one has oedema desc oedema tab oedema *** Variable: Sampling weight *** gen sw = hhweight //For household sample weight desc sw sum sw /* For this part of the do-file we use the WHO AnthroPlus software. This is to calculate the z-scores for young children aged 60-72 months. Source of ado file: https://www.who.int/growthref/tools/en/ */ *** Next, indicate to STATA where the igrowup_restricted.ado file is stored: adopath + "$path_ado/who2007_stata" /* We use 'reflib' to specify the package directory where the .dta files containing the WHO Growth reference are stored. Note that we use strX to specity the length of the path in string. */ gen str100 reflib = "$path_ado/who2007_stata" lab var reflib "Directory of reference tables" /* We use datalib to specify the working directory where the input STATA data set containing the anthropometric measurement is stored. */ gen str100 datalib = "$path_out" lab var datalib "Directory for datafiles" /* We use datalab to specify the name that will prefix the output files that will be produced from using this ado file*/ gen str30 datalab = "children_nutri_syr" lab var datalab "Working file" /*We now run the command to calculate the z-scores with the adofile */ who2007 reflib datalib datalab gender age_months ageunit weight /// height oedema sw /*We now turn to using the dta file that was created and that contains the calculated z-scores to compute BMI-for-age*/ use "$path_out/children_nutri_syr_z.dta", clear gen z_bmi = _zbfa replace z_bmi = . if _fbfa==1 lab var z_bmi "z-score bmi-for-age WHO" *** Standard MPI indicator *** gen low_bmiage = (z_bmi < -2.0) /*Takes value 1 if BMI-for-age is under 2 stdev below the median & 0 otherwise */ replace low_bmiage = . if z_bmi==. lab var low_bmiage "Teenage low bmi 2sd - WHO" *** Destitution indicator *** gen low_bmiage_u = (z_bmi < -3.0) /*Takes value 1 if BMI-for-age is under 3 stdev below the median & 0 otherwise */ replace low_bmiage_u = . if z_bmi==. lab var low_bmiage_u "Teenage very low bmi 3sd - WHO" tab low_bmiage, miss tab low_bmiage_u, miss //Retain relevant variables: keep ind_id child_CH age_months low_bmiage* order ind_id child_CH age_months low_bmiage* sort ind_id save "$path_out/SYR09_CH_6Y.dta", replace //Erase files from folder: erase "$path_out/children_nutri_syr_z.xls" erase "$path_out/children_nutri_syr_prev.xls" erase "$path_out/children_nutri_syr_z.dta" ******************************************************************************** *** Step 1.2 BR - BIRTH RECODE *** (All females 15-49 years who ever gave birth) ******************************************************************************** /* Note: There is no birth history data file for Syria PAPFAM 2009. Hence this section has been deactivated */ ******************************************************************************** *** Step 1.3 IR - WOMEN's RECODE *** (All eligible females 15-49 years in the household) ******************************************************************************** use "$path_in/Syria 2009 - Wom.dta", clear rename _all, lower *** Generate individual unique key variable required for data merging *** cluster=cluster number; *** hhnum=household number; *** w_ln=respondent's line number gen double ind_id = cluster*1000000 + hhnum*100 + w_ln format ind_id %20.0g label var ind_id "Individual ID" duplicates report ind_id gen women_WM=1 //Identification variable for observations in IR recode //Retain relevant variables: keep ind_id women_WM wmweight wresult w201 w206 w207a w207b /// w103c w208 wresult w124 w125 w104 w105 w106 w107 order ind_id women_WM wmweight wresult w201 w206 w207a w207b /// w103c w208 wresult w124 w125 w104 w105 w106 w107 sort ind_id save "$path_out/SYR09_WM.dta", replace ******************************************************************************** *** Step 1.4 HH - Household's recode *** ******************************************************************************** use "$path_in/Syria 2009 - HH.dta", clear rename _all, lower *** Generate individual unique key variable required for data merging *** cluster=cluster number; *** hhnum=household number; gen double hh_id = cluster*100 + hhnum format hh_id %20.0g lab var hh_id "Household ID" save "$path_out/SYR09_HH.dta", replace ******************************************************************************** *** Step 1.5 HR - Household Member's recode **** ******************************************************************************** use "$path_in/Syria 2009 - HR.dta", clear rename _all, lower *** Generate a household unique key variable at the household level using: ***cluster=cluster number ***hhnum=household number gen double hh_id = cluster*100 + hhnum format hh_id %20.0g label var hh_id "Household ID" *** Generate individual unique key variable required for data merging using: *** cluster=cluster number; *** hhnum=household number; *** ln=respondent's line number. gen double ind_id = cluster*1000000 + hhnum*100 + ln format ind_id %20.0g label var ind_id "Individual ID" ******************************************************************************** *** Step 1.6 DATA MERGING ******************************************************************************** *** Merging WM Recode ***************************************** merge 1:1 ind_id using "$path_out/SYR09_WM.dta" tab wresult women_WM, miss col bys hh_id: egen temp=sum(women_WM) tab q101w temp, miss tab q101c temp, miss count if temp==0 & q101c >=1 //NOTE: There is 497 women not eligible but with child measures drop temp _merge erase "$path_out/SYR09_WM.dta" *** Merging HH Recode ***************************************** merge m:1 hh_id using "$path_out/SYR09_HH.dta" tab result if _m==2 drop if _merge==2 //Drop households that were not interviewed drop _merge erase "$path_out/SYR09_HH.dta" *** Merging children 0-60 months ***************************************** merge 1:1 ind_id using "$path_out/SYR09_CH.dta" drop _merge erase "$path_out/SYR09_CH.dta" *** Merging children 61-72 months ***************************************** merge 1:1 ind_id using "$path_out/SYR09_CH_6Y.dta" drop _merge erase "$path_out/SYR09_CH_6Y.dta" sort ind_id ******************************************************************************** *** Step 1.7 KEEPING ONLY DE JURE HOUSEHOLD MEMBERS *** ******************************************************************************** //Permanent (de jure) household members clonevar resident = h102 codebook resident, tab (10) label var resident "Permanent (de jure) household member" drop if resident!=1 tab resident, miss /*Note: The Global MPI is based on de jure (permanent) household members only. As such, non-usual residents will be excluded from the sample. However, in Libya PAPFAM 2009, all householders are permanent members. */ ******************************************************************************** *** Step 1.8 SUBSAMPLE VARIABLE *** ******************************************************************************** /* In the context of Syria PAPFAM 2009, height and weight measurements were collected from all children (0-5). As such there is no presence of subsample. */ gen subsample = . label var subsample "Households selected as part of nutrition subsample" tab subsample, miss ******************************************************************************** *** Step 1.9 CONTROL VARIABLES ******************************************************************************** /* Households are identified as having 'no eligible' members if there are no applicable population, that is, children 0-5 years, adult women 15-49 years or adult men. These households will not have information on relevant indicators of health. As such, these households are considered as non-deprived in those relevant indicators.*/ *** No Eligible Women ***************************************** gen fem_eligible = (women_WM==1) bys hh_id: egen hh_n_fem_eligible = sum(fem_eligible) //Number of eligible women for interview in the hh gen no_fem_eligible = (hh_n_fem_eligible==0) //Takes value 1 if the household had no eligible females for an interview tab no_fem_eligible, miss tab no_fem_eligible q101c, miss /* NOTE: There is 497 individuals living in households without eligible women but have child who was eligible for anthropometric measures.*/ lab var no_fem_eligible "Household has no eligible women" *** No Eligible Men ***************************************** //NOTE: Syria PAPFAM 2009 have no male recode file gen no_male_eligible = . lab var no_male_eligible "Household has no eligible man" tab no_male_eligible, miss *** No Eligible Children 0-60 months ***************************************** gen child_eligible = 0 replace child_eligible = 1 if q101c>=1 & (agemonth>=0 & agemonth<=60) bys hh_id: egen hh_n_children_eligible = sum(child_eligible) //Number of eligible children for anthropometrics gen no_child_eligible = (hh_n_children_eligible==0) //Takes value 1 if there were no eligible children for anthropometrics tab no_child_eligible,miss lab var no_child_eligible "Household has no children under 5 eligible" *** No Eligible Children 61-72 months ***************************************** gen child_eligible_6y = 0 replace child_eligible_6y = 1 if q101c>=1 & (agemonth>=61 & agemonth<=72) bys hh_id: egen hh_n_children_eligible_6y = sum(child_eligible_6y) //Number of eligible children for anthropometrics gen no_child_eligible_6y = (hh_n_children_eligible_6y==0) //Takes value 1 if there were no eligible children for anthropometrics tab no_child_eligible_6y, miss lab var no_child_eligible_6y "Household has no children 5-6 years eligible" *** No Eligible Women and Men *********************************************** /*Syria PAPFAM 2009 enumerated men, as household members but did not collect child mortality information from men. As such this variable is generated as an empty variable*/ gen no_adults_eligible = . lab var no_adults_eligible "Household has no eligible women or men" tab no_adults_eligible, miss *** No Eligible Children and Women *********************************************** /*NOTE: In the DHS datasets, we use this variable as a control variable for the nutrition indicator if nutrition data is present for children and women. However, in this dataset, nutrition data only covers children. So we generate an empty variable */ gen no_child_fem_eligible = . lab var no_child_fem_eligible "Household has no children or women eligible" tab no_child_fem_eligible, miss *** No Eligible Women, Men or Children *********************************************** /*NOTE: In the DHS datasets, we use this variable as a control variable for the nutrition indicator if nutrition data is present for children, women and men. There is no data for men. So we generate this variable as an empty variable. */ gen no_eligibles = . lab var no_eligibles "Household has no eligible women, men, or children" tab no_eligibles, miss *** No Eligible Subsample ***************************************** //Note that PAPFAM surveys do not collect hemoglobin data from householders gen no_hem_eligible = . lab var no_hem_eligible "Household has no eligible individuals for hemoglobin measurements" tab no_hem_eligible, miss drop fem_eligible hh_n_fem_eligible child_eligible hh_n_children_eligible sort hh_id ind_id ******************************************************************************** *** Step 1.11 RENAMING DEMOGRAPHIC VARIABLES *** ******************************************************************************** //Sample weight clonevar weight = hhweight label var weight "Sample weight" //Area: urban or rural codebook area , tab (5) replace area=0 if area==2 label define lab_area 1 "urban" 0 "rural" label values area lab_area label var area "Area: urban-rural" tab area, miss //Sex of household member codebook h103 clonevar sex = h103 label var sex "Sex of household member" //Age of household member codebook xh105a, tab (100) clonevar age = xh105a label var age "Age of household member" //Age group recode age (0/4 = 1 "0-4")(5/9 = 2 "5-9")(10/14 = 3 "10-14") /// (15/17 = 4 "15-17")(18/59 = 5 "18-59")(60/max=6 "60+"), gen(agec7) lab var agec7 "age groups (7 groups)" recode age (0/9 = 1 "0-9") (10/17 = 2 "10-17")(18/59 = 3 "18-59") /// (60/max=4 "60+"), gen(agec4) lab var agec4 "age groups (4 groups)" recode age (0/17 = 1 "0-17") (18/max = 2 "18+"), gen(agec2) lab var agec2 "age groups (2 groups)" //Total number of de jure hh members in the household gen member = 1 bysort hh_id: egen hhsize = sum(member) label var hhsize "Household size" tab hhsize, miss drop member //Subnational region /* (RG!): There is no district info for Syria PAPFAM 2009, but according to Wikipedia (!) Syria is divided in 14 Governorated (https://en.wikipedia.org/wiki/Governorates_of_Syria) and 61 districts. No var exists for districts but there is one var for Governorates: gov. Thus, I used gov and the Subnational representation for Syria PAPFAM 2009 is at Governorate level. */ codebook gov, tab (99) clonevar region = gov lab var region "Region for subnational decomposition" codebook region, tab (99) label define lab_reg /// 1 "Damascus" /// 2 "Halab" /// 3 "Rif Dimashq" /// 4 "Homs" /// 5 "Hama" /// 6 "Latakia" /// 7 "Idlib" /// 8 "Al-Hasakah" /// 9 "Deir ez-Zor" /// 10 "Tartus" /// 11 "Raqa" /// 12 "Daraa" /// 13 "As-Suwayda" /// 14 "Quneitra" label values region lab_reg ******************************************************************************** *** Step 2 Data preparation *** *** Standardization of the 10 Global MPI indicators *** Identification of non-deprived & deprived individuals ******************************************************************************** ******************************************************************************** *** Step 2.1 Years of Schooling *** ******************************************************************************** tab h110a h110b, miss //h110a: School Attendance; h110b: Highest certificate obtained tab age h110b if h111ba==1, miss /*h111a: Curretly attended (ever attended in 2009/2010); h11ba: Level attended in 2009/2010 */ rename edulevel edulevel_ori gen edulevel = h110b //Highest educational level attended replace edulevel = . if h110b==. | h110b==8 //8: DK, check that that is also the case in your data bys h110a: tab h110b, miss replace edulevel = 0 if h110a == 3 /*h110a: School attendance. Replacing as no level of education those individuals, who never attended to school */ label define lab_edulevel 0 "None" 1 "Attended school " 2 "Primary" /// 3 "Preparatory" 4 "Basic" 5 "Secondary" /// 6 "Middle Institute" 7 "University +" label values edulevel lab_edulevel label var edulevel "Highest educational level attended" gen eduhighyear = 0 //Highest grade of education completed// replace eduhighyear = 0 if edulevel == 0 replace eduhighyear = 1 if edulevel == 1 replace eduhighyear = 6 if edulevel == 2 /*Primary school 6 years; Basic education: 9 years (primary 6 + middle 3 years); secondary education: 3 years */ replace eduhighyear = 6+3 if edulevel == 3 | edulevel== 4 /*Preparatory education = middle education in the old education system in Syria */ replace eduhighyear = 6+3+3 if edulevel == 5 replace eduhighyear = 13 if (edulevel == 6 | edulevel == 7) replace eduhighyear = . if edulevel==. //These are considered missing values replace eduhighyear = 0 if h110a == 3 //Never attended school lab var eduhighyear "Highest year of education completed" gen eduhighyear2 = eduhighyear /*Using information on attendance to complete highest grade of education completed */ tab eduhighyear*, miss replace eduhighyear = h111bb if h111a == 1 & h111ba==1 & h111bb<88 tab eduhighyear*, miss replace eduhighyear = h111db if eduhighyear<6 & h111a== 2 & h111c == 1 & h111da==1 & h111db<88 //only 413 cases that did no attend the year of survey but did the previous year tab eduhighyear*, miss replace eduhighyear = h111bb+9 if h111a == 1 & (h111ba==2 | h111ba==3) & h111bb<88 //add 9 years for secondary or middle education tab eduhighyear*, miss replace eduhighyear = h111bb+9 if h111a== 2 & h111c == 1 & (h111da==2 | h111da==3) & h111db<88 //only 736 cases that did no attend the year of survey but did the previous year tab eduhighyear*, miss replace eduhighyear = h111bb+12 if h111a == 1 & h111ba==4 & h111bb<88 //add 9 years for middle education tab eduhighyear*, miss replace eduhighyear = h111db+12 if h111a== 2 & h111c == 1 & h111da==4 & h111db<88 // only 736 cases that did no attend the year of survey but did the previous year *We retrieve info for women eligible replace eduhighyear = w106 if women_WM==1 & w105==1 //Primary level attended replace eduhighyear = w106+6 if women_WM==1 & w105==2 //Preparatory level attended replace eduhighyear = w106+9 if women_WM==1 & w105==3 //Secondary level replace eduhighyear = w106+9 if women_WM==1 & w105==4 //Middle institute replace eduhighyear = w106+12 if women_WM==1 & w105==5 //University replace eduhighyear = 0 if women_WM==1 & w104==3 //Never attended school ** Cleaning inconsistencies replace eduhighyear = 0 if age < 10 /*The variable "eduyears" was replaced with a '0' given that the criteria for this indicator is household member aged 10 years or older */ ** Now we create the years of schooling gen eduyears = eduhighyear replace eduyears = . if eduyears>30 replace eduyears = . if edulevel==. ** Checking for further inconsistencies replace eduyears = . if age<=eduyears & age>0 replace eduyears = 0 if age<10 lab var eduyears "Total number of years of education accomplished" /*A control variable is created on whether there is information on years of education for at least 2/3 of the household members aged 10 years and older */ gen temp = 1 if eduyears!=. & age>=10 & age!=. bysort hh_id: egen no_missing_edu = sum(temp) /*Total household members who are 10 years and older with no missing years of education */ gen temp2 = 1 if age>=10 & age!=. bysort hh_id: egen hhs = sum(temp2) //Total number of household members who are 10 years and older replace no_missing_edu = no_missing_edu/hhs replace no_missing_edu = (no_missing_edu>=2/3) /*Identify whether there is information on years of education for at least 2/3 of the household members aged 10 years and older */ tab no_missing_edu, miss label var no_missing_edu "No missing edu for at least 2/3 of the HH members aged 10 years & older" drop temp temp2 hhs *** Standard MPI *** /*The entire household is considered deprived if no household member aged 10 years or older has completed SIX years of schooling.*/ ******************************************************************* gen years_edu6 = (eduyears>=6) /* The years of schooling indicator takes a value of "1" if at least someone in the hh has reported 6 years of education or more */ replace years_edu6 = . if eduyears==. bysort hh_id: egen hh_years_edu6_1 = max(years_edu6) gen hh_years_edu6 = (hh_years_edu6_1==1) replace hh_years_edu6 = . if hh_years_edu6_1==. replace hh_years_edu6 = . if hh_years_edu6==0 & no_missing_edu==0 lab var hh_years_edu6 "Household has at least one member with 6 years of edu" tab hh_years_edu6, miss *** Destitution MPI *** /*The entire household is considered deprived if no household member aged 10 years or older has completed at least one year of schooling.*/ ******************************************************************* gen years_edu1 = (eduyears>=1) replace years_edu1 = . if eduyears==. bysort hh_id: egen hh_years_edu_u = max(years_edu1) replace hh_years_edu_u = . if hh_years_edu_u==0 & no_missing_edu==0 lab var hh_years_edu_u "Household has at least one member with 1 year of edu" ******************************************************************************** *** Step 2.2 Child School Attendance *** ******************************************************************************** codebook h110a, tab (10) clonevar attendance = h110a //1=attending, 0=not attending recode attendance (2=0) (3=0) //2='attended in the past'; 3='never attended' replace attendance = . if attendance==9 //Missing values replaced *** Standard MPI *** /*The entire household is considered deprived if any school-aged child is not attending school up to class 8. */ ******************************************************************* gen child_schoolage = (age>=6 & age<=14) /* Note: In Syrian Arab Republic, the official school entrance age is 6 years. So, age range is 6-14 (=6+8) Source: "http://data.uis.unesco.org/?ReportId=163" Go to Education>Education>System>Official entrance age to primary education. Look at the starting age and add 8. */ /*A control variable is created on whether there is no information on school attendance for at least 2/3 of the school age children */ count if child_schoolage==1 & attendance==. //Understand how many eligible school aged children are not attending school gen temp = 1 if child_schoolage==1 & attendance!=. bysort hh_id: egen no_missing_atten = sum(temp) /*Total school age children with no missing information on school attendance */ gen temp2 = 1 if child_schoolage==1 bysort hh_id: egen hhs = sum(temp2) //Total number of household members who are of school age replace no_missing_atten = no_missing_atten/hhs replace no_missing_atten = (no_missing_atten>=2/3) /*Identify whether there is missing information on school attendance for more than 2/3 of the school age children */ tab no_missing_atten, miss label var no_missing_atten "No missing school attendance for at least 2/3 of the school aged children" drop temp temp2 hhs bysort hh_id: egen hh_children_schoolage = sum(child_schoolage) replace hh_children_schoolage = (hh_children_schoolage>0) //Control variable: //It takes value 1 if the household has children in school age lab var hh_children_schoolage "Household has children in school age" gen child_not_atten = (attendance==0) if child_schoolage==1 replace child_not_atten = . if attendance==. & child_schoolage==1 bysort hh_id: egen any_child_not_atten = max(child_not_atten) gen hh_child_atten = (any_child_not_atten==0) replace hh_child_atten = . if any_child_not_atten==. replace hh_child_atten = 1 if hh_children_schoolage==0 replace hh_child_atten = . if hh_child_atten==1 & no_missing_atten==0 /*If the household has been intially identified as non-deprived, but has missing school attendance for at least 2/3 of the school aged children, then we replace this household with a value of '.' because there is insufficient information to conclusively conclude that the household is not deprived */ lab var hh_child_atten "Household has all school age children up to class 8 in school" tab hh_child_atten, miss /*Note: The indicator takes value 1 if ALL children in school age are attending school and 0 if there is at least one child not attending. Households with no children receive a value of 1 as non-deprived. The indicator has a missing value only when there are all missing values on children attendance in households that have children in school age. */ *** Destitution MPI *** /*The entire household is considered deprived if any school-aged child is not attending school up to class 6. */ ******************************************************************* gen child_schoolage_6 = (age>=6 & age<=12) /*Note: in Syria, the official school entrance age is 6 years. So, age range for destitute measure is 6-12 (=6+6) */ /*A control variable is created on whether there is no information on school attendance for at least 2/3 of the children attending school up to class 6 */ gen temp = 1 if child_schoolage_6==1 & attendance!=. bysort hh_id: egen no_missing_atten_u = sum(temp) /*Total school age children attending up to class 6 with no missing information on school attendance */ gen temp2 = 1 if child_schoolage_6==1 bysort hh_id: egen hhs = sum(temp2) /*Total number of household members who are of school age attending up to class 6 */ replace no_missing_atten_u = no_missing_atten_u/hhs replace no_missing_atten_u = (no_missing_atten_u>=2/3) /*Identify whether there is missing information on school attendance for more than 2/3 of the school age children attending up to class 6 */ tab no_missing_atten_u, miss label var no_missing_atten_u "No missing school attendance for at least 2/3 of the school aged children" drop temp temp2 hhs bysort hh_id: egen hh_children_schoolage_6 = sum(child_schoolage_6) replace hh_children_schoolage_6 = (hh_children_schoolage_6>0) lab var hh_children_schoolage_6 "Household has children in school age (6 years of school)" gen child_atten_6 = (attendance==1) if child_schoolage_6==1 replace child_atten_6 = . if attendance==. & child_schoolage_6==1 bysort hh_id: egen any_child_atten_6 = max(child_atten_6) gen hh_child_atten_u = (any_child_atten_6==1) replace hh_child_atten_u = . if any_child_atten_6==. replace hh_child_atten_u = 1 if hh_children_schoolage_6==0 replace hh_child_atten_u = . if hh_child_atten_u==0 & no_missing_atten_u==0 /*If the household has been intially identified as deprived, but has missing school attendance for at least 2/3 of the school aged children, then we replace this household with a value of '.' because there is insufficient information to conclusively conclude that the household is deprived */ lab var hh_child_atten_u "Household has at least one school age children up to class 6 in school" tab hh_child_atten_u, miss ******************************************************************************** *** Step 2.3 Nutrition *** ******************************************************************************** /*Please note that the PAPFAM datasets do not collect nutrition data from adults. In the context of countries with PAPFAM datasets, the entire household is considered deprived if any child 0-60 months for whom there is nutritional information is malnourished or children 61-72 months who have low BMI-for-age in the household.*/ ******************************************************************************** *** Step 2.3a Child 0-60 months Nutrition *** ******************************************************************************** /* Syria PAPFAM 2009 collected nutrition data from children 0-60 months. In this section, the construction of the nutrition indicator will be on children 0-60 months. Households with no eligible children will receive a value of 1 */ *** Standard MPI: Child Underweight Indicator *** ************************************************************************ bysort hh_id: egen temp = max(underweight) gen hh_no_underweight = (temp==0) //Takes value 1 if no child in the hh is underweight replace hh_no_underweight = . if temp==. replace hh_no_underweight = 1 if no_child_eligible==1 //Households with no eligible children will receive a value of 1 lab var hh_no_underweight "Household has no child underweight - 2 stdev" drop temp *** Standard MPI: Child Stunting Indicator *** ************************************************************************ bysort hh_id: egen temp = max(stunting) gen hh_no_stunting = (temp==0) //Takes value 1 if no child in the hh is stunted replace hh_no_stunting = . if temp==. replace hh_no_stunting = 1 if no_child_eligible==1 lab var hh_no_stunting "Household has no child stunted - 2 stdev" drop temp *** Standard MPI: Child Either Stunted or Underweight Indicator *** ************************************************************************ gen uw_st = 1 if stunting==1 | underweight==1 replace uw_st = 0 if stunting==0 & underweight==0 replace uw_st = . if stunting==. & underweight==. bysort hh_id: egen temp = max(uw_st) gen hh_no_uw_st = (temp==0) //Takes value 1 if no child in the hh is underweight or stunted replace hh_no_uw_st = . if temp==. replace hh_no_uw_st = 1 if no_child_eligible==1 //Households with no eligible children will receive a value of 1 lab var hh_no_uw_st "Household has no child underweight or stunted" drop temp *** Destitute MPI: Child Underweight Indicator *** ************************************************************************ bysort hh_id: egen temp = max(underweight_u) gen hh_no_underweight_u = (temp==0) //Takes value 1 if no child in the hh is underweighted, 0 if at least one is replace hh_no_underweight_u = . if temp==. replace hh_no_underweight_u = 1 if no_child_eligible==1 //Households with no eligible children will receive a value of 1 lab var hh_no_underweight_u "Destitute: Household has no child underweight" drop temp *** Destitute MPI: MPI Child Stunting Indicator *** ************************************************************************ bysort hh_id: egen temp = max(stunting_u) gen hh_no_stunting_u = (temp==0) //Takes value 1 if no child in the hh is underweighted; 0 if at least one is replace hh_no_stunting_u = . if temp==. replace hh_no_stunting_u = 1 if no_child_eligible==1 //Households with no eligible children will receive a value of 1 lab var hh_no_stunting_u "Destitute: Household has no child stunted" drop temp *** Destitute MPI: Child Either Underweight or Stunted Indicator *** ************************************************************************ gen uw_st_u = 1 if stunting_u==1 | underweight_u==1 replace uw_st_u = 0 if stunting_u==0 & underweight_u==0 replace uw_st_u = . if stunting_u==. & underweight_u==. bysort hh_id: egen temp = max(uw_st_u) gen hh_no_uw_st_u = (temp==0) //Takes value 1 if no child in the hh is underweight or stunted replace hh_no_uw_st_u = . if temp==. replace hh_no_uw_st_u = 1 if no_child_eligible==1 //Households with no eligible children will receive a value of 1 lab var hh_no_uw_st_u "Destitute: Household has no child underweight or stunted" drop temp ******************************************************************************** *** Step 2.3b Child 5-6 years Nutrition *** ******************************************************************************** /* Syria PAPFAM 2009 collected nutrition data from children under 6. In this section, the construction of the nutrition indicator will be for children between 61-72 months. Households with no eligible children will receive a value of 1 */ *** Standard MPI: Child BMI-for-age Indicator *** ************************************************************************ bysort hh_id: egen temp = max(low_bmiage) gen hh_no_low_bmiage = (temp==0) //Takes value 1 if no child in the hh has low BMI-for-age replace hh_no_low_bmiage = . if temp==. replace hh_no_low_bmiage = 1 if no_child_eligible_6y==1 //Households with no eligible children will receive a value of 1 lab var hh_no_low_bmiage "Household has no child low BMI-for-age" drop temp *** Destitute MPI: Child BMI-for-age Indicator*** ************************************************************************ bysort hh_id: egen temp = max(low_bmiage_u) gen hh_no_low_bmiage_u = (temp==0) //Takes value 1 if no child in the hh has low BMI-for-age replace hh_no_low_bmiage_u = . if temp==. replace hh_no_low_bmiage_u = 1 if no_child_eligible_6y==1 //Households with no eligible children will receive a value of 1 lab var hh_no_low_bmiage_u "Destitute: Household has no child low BMI-for-age" drop temp ******************************************************************************** *** Step 2.3c Household Nutrition Indicator *** ******************************************************************************** *** Standard MPI *** /* Members of the household are considered deprived if the household has a child under 5 whose height-for-age or weight-for-age is under two standard deviation below the median, or has teenager with BMI-for-age that is under two standard deviation below the median, or has adults with BMI threshold that is below 18.5 kg/m2. Households that have no eligible adult AND no eligible children are considered non-deprived. The indicator takes a value of missing only if all eligible adults and eligible children have missing information in their respective nutrition variable. */ ************************************************************************ gen hh_nutrition_uw_st = 1 replace hh_nutrition_uw_st = 0 if hh_no_uw_st==0 | hh_no_low_bmiage==0 replace hh_nutrition_uw_st = . if hh_no_uw_st==. & hh_no_low_bmiage==. replace hh_nutrition_uw_st = 1 if no_child_eligible==1 & no_child_eligible_6y==1 /*We replace households that do not have the applicable population, as non-deprived in nutrition*/ lab var hh_nutrition_uw_st "Household has no child underweight or stunted" tab hh_nutrition_uw_st, miss *** Destitution MPI *** /* Members of the household are considered deprived if the household has a child under 5 whose height-for-age or weight-for-age is under three standard deviation below the median, or has teenager with BMI-for-age that is under three standard deviation below the median, or has adults with BMI threshold that is below 17.0 kg/m2. Households that have no eligible adult AND no eligible children are considered non-deprived. The indicator takes a value of missing only if all eligible adults and eligible children have missing information in their respective nutrition variable. */ ************************************************************************ gen hh_nutrition_uw_st_u = 1 replace hh_nutrition_uw_st_u = 0 if hh_no_uw_st_u==0 | hh_no_low_bmiage_u==0 replace hh_nutrition_uw_st_u = . if hh_no_uw_st_u==. & hh_no_low_bmiage_u==. replace hh_nutrition_uw_st_u = 1 if no_child_eligible==1 & no_child_eligible_6y==1 /*We replace households that do not have the applicable population, as non-deprived in nutrition*/ lab var hh_nutrition_uw_st_u "Household has no child underweight or stunted (destitute)" tab hh_nutrition_uw_st_u, miss ******************************************************************************** *** Step 2.4 Child Mortality *** ******************************************************************************** /*In the context of Syria PAPFAM 2009, information on child mortality was collected only from women */ codebook w206 w207a w207b // w206: Had children who died // w207a: number of sons who have died // w207b: number of daughters who have died egen temp_f = rowtotal(w207a w207b), missing //Total child mortality reported by eligible women replace temp_f = 0 if (w201==1 & w206!=1) | (w201==2 & w206!=1) replace temp_f = 0 if w201==. & w206==. & h106==1 & temp_f==. /*Assign a value of "0" for: - all eligible women who ever had a life birth but reported no child death - all eligible women who never had a life birth and reported no death (presumably this group are women who never ever gave birth) - all elegible but never married women were not asked birth history questions and hence we assume there is no child mortality among this group*/ bysort hh_id: egen child_mortality_f = sum(temp_f), missing lab var child_mortality_f "Occurrence of child mortality reported by women" tab child_mortality_f, miss drop temp_f egen child_mortality = rowmax(child_mortality_f) lab var child_mortality "Total child mortality within household reported by women & men" tab child_mortality, miss *** Standard MPI *** /*Deprived if any children under 18 died in the household in the last 5 years from the survey year */ ************************************************************************ /*In the case of Syria, there is no birth history data. This means, there is no information on the date of death of children who have died. As such we are not able to construct the indicator on child mortality under 18 that occurred in the last 5 years. Instead, we identify individuals as deprived if any children died in the household */ gen hh_mortality = (child_mortality==0) /*Household is replaced with a value of "1" if there is no incidence of child mortality*/ replace hh_mortality = . if child_mortality==. replace hh_mortality = 1 if no_fem_eligible==1 //Household is replaced with a value of "1" if there is no eligible women lab var hh_mortality "Household had no child mortality" tab hh_mortality, miss gen hh_mortality_u18_5y = . lab var hh_mortality_u18_5y "Household had no under 18 child mortality in the last 5 years" *** Destitution MPI *** *** (same as standard MPI) *** ************************************************************************ clonevar hh_mortality_u = hh_mortality ******************************************************************************** *** Step 2.5 Electricity *** ******************************************************************************** /* Note: Syria PAPFAM 2009 has no direct question on whether household has electricity or not. As the best alternative, the electricity indicator for Syria PAPFAM 2009 was drawn from the h517 variable: Main type of lighting. The categoreis are: Electricity; Kerosene; Gas; Oil/Candles; Other; No lighting. As such, the category 'Electricity' is recoded as 'Yes electricity' and all other categories are recoded as 'No electricity' */ *** Standard MPI *** /*Members of the household are considered deprived if the household has no electricity */ *************************************************** lookfor electricity lighting codebook h517, tab (10) //Main type of lighting in Syria PAPFAM 2009 clonevar electricity = h517 recode electricity (2/8=0) codebook electricity, tab (10) label var electricity "Household has electricity" *** Destitution MPI *** *** (same as standard MPI) *** *************************************************** gen electricity_u = electricity label var electricity_u "Household has electricity" ******************************************************************************** *** Step 2.6 Sanitation *** ******************************************************************************** /* Improved sanitation facilities include flush or pour flush toilets to sewer systems, septic tanks or pit latrines, ventilated improved pit latrines, pit latrines with a slab, and composting toilets. These facilities are only considered improved if it is private, that is, it is not shared with other households. Source: https://unstats.un.org/sdgs/metadata/files/Metadata-06-02-01.pdf Note: In cases of mismatch between the country report and the internationally agreed guideline, we followed the report. */ lookfor toilet clonevar toilet = h514 codebook toilet, tab (30) codebook h515, tab(10) recode h515(2=0)(9=.),gen(shared_toilet) replace shared_toilet=1 if h514==5 & h515==. tab h515 shared_toilet, miss *** Standard MPI *** /*Members of the household are considered deprived if the household's sanitation facility is not improved (according to the SDG guideline) or it is improved but shared with other households*/ ******************************************************************** /* NOTE: The toilet categories for Syria PAPFAM 2009 are different from the standardised version found in DHS and MICS. The categories are: 1 FT connected 2 FT not connected 3 Toilet connected 4 Toilet connected to closed pit 5 Public toilet 6 Open air 96 Other The categories of public toilet, open air & other are coded as non-improved sanitation. */ gen toilet_mdg = toilet<=4 & shared_toilet!=1 /*Household is assigned a value of '1' if it uses improved sanitation and does not share toilet with other households */ replace toilet_mdg = 0 if toilet<=4 & shared_toilet==1 /*Household is assigned a value of '0' if it uses improved sanitation but shares toilet with other households */ replace toilet_mdg = . if toilet==. //Household is assigned a value of '.' if it has missing information lab var toilet_mdg "Household has improved sanitation with MDG Standards" tab toilet toilet_mdg, miss *** Destitution MPI *** /*Members of the household are considered deprived if household practises open defecation or uses other unidentifiable sanitation practises */ ******************************************************************** gen toilet_u = . replace toilet_u = 0 if toilet==6 | toilet==96 /*Household is assigned a value of '0' if it practises open defecation or others */ replace toilet_u = 1 if toilet!=6 & toilet!=96 & toilet!=. /*Household is assigned a value of '1' if it does not practise open defecation or others */ lab var toilet_u "Household does not practise open defecation or others" tab toilet toilet_u, miss ******************************************************************************** *** Step 2.7 Drinking Water *** ******************************************************************************** /* Improved drinking water sources include the following: piped water into dwelling, yard or plot; public taps or standpipes; boreholes or tubewells; protected dug wells; protected springs; packaged water; delivered water and rainwater which is located on premises or is less than a 30-minute walk from home roundtrip. Source: https://unstats.un.org/sdgs/metadata/files/Metadata-06-01-01.pdf Note: In cases of mismatch between the country report and the internationally agreed guideline, we followed the report. */ lookfor water clonevar water = h508 codebook water, tab (30) clonevar timetowater = h510 codebook timetowater, tab (999) gen ndwater = . //No observation for non-drinking water *** Standard MPI *** /* Members of the household are considered deprived if the household does not have access to improved drinking water (according to the SDG guideline) or safe drinking water is at least a 30-minute walk from home, roundtrip */ ******************************************************************** /* In the case of Syria PAPFAM 2009, non deprived if water is 1 "piped supply", 2 "public tap", 3 "artesian well", 5 "supervised spring", 9 "rainwater" 11 "bottled water", Deprived if water is 4 "regular well", 6 "unsupervised spring", 7 "river", 8 "lake", 10 "tanker truck", 96 "other" */ gen water_mdg = 1 if (water>=1 & water <=3)| water==5 | water==9 | water==11 replace water_mdg = 0 if water== 4 | water==6 | water==7 | water==8 | /// water==10 | water==96 replace water_mdg = 0 if water_mdg==1 & timetowater >= 30 & timetowater!=. & /// timetowater!=998 //Deprived if water is at more than 30 minutes' walk (roundtrip) replace water_mdg = . if water==. | water==99 lab var water_mdg "Household has drinking water with MDG standards (considering distance)" tab water water_mdg, miss *** Destitution MPI *** /* Members of the household is identified as destitute if household does not have access to safe drinking water, or safe water is more than 45 minute walk from home, round trip.*/ ******************************************************************** gen water_u = . replace water_u = 1 if (water>=1 & water <=3)| water==5 | water==9 | water==11 replace water_u = 0 if water== 4 | water==6 | water==7 | water==8 | /// water==10 | water==96 replace water_u = 0 if water_u==1 & timetowater>45 & timetowater!=. & /// timetowater!=998 replace water_u = . if water==99 | water==. lab var water_u "Household has drinking water with MDG standards (45 minutes distance)" tab water water_u, miss ******************************************************************************** *** Step 2.8 Housing *** ******************************************************************************** /* Members of the household are considered deprived if the household has a dirt, sand or dung floor */ lookfor floor clonevar floor = h503 codebook floor, tab (10) gen floor_imp = 1 replace floor_imp = 0 if floor==1 | floor==6 replace floor_imp = . if floor==. | floor==9 replace floor_imp = 0 if floor==. & h601 >=5 lab var floor_imp "Household has floor that it is not earth/sand/dung" tab floor floor_imp, miss /* Members of the household are considered deprived if the household has wall made of natural or rudimentary materials */ lookfor wall gen wall = . gen wall_imp = . lab var wall_imp "Household has wall that it is not of low quality materials" /* Members of the household are considered deprived if the household has roof made of natural or rudimentary materials */ lookfor roof gen roof = . gen roof_imp = . lab var roof_imp "Household has roof that it is not of low quality materials" *** Standard MPI *** /* Members of the household is deprived in housing if the roof, floor OR walls are constructed from low quality materials. Since Syria PAPFAM 2009 do not have information on walls and roof, we replace the MPI indicator on housing with information on floor.*/ ************************************************************** gen housing_1 = floor_imp lab var housing_1 "Household has roof, floor & walls that it is not low quality material" tab housing_1, miss *** Destitution MPI *** /* Members of the household is deprived in housing if two out of three components (roof and walls; OR floor and walls; OR roof and floor) the are constructed from low quality materials. Since Syria PAPFAM 2009 do not have information on walls and roof, we replace the MPI indicator on housing with information on floor.*/ ************************************************************** gen housing_2 = floor_imp rename housing_2 housing_u lab var housing_u "Household has one of three aspects(either roof,floor/walls) that is not low quality material" ******************************************************************************** *** Step 2.9 Cooking Fuel *** ******************************************************************************** /* Solid fuel are solid materials burned as fuels, which includes coal as well as solid biomass fuels (wood, animal dung, crop wastes and charcoal). Source: https://apps.who.int/iris/bitstream/handle/10665/141496/9789241548885_eng.pdf */ lookfor fuel clonevar cookingfuel = h519 codebook cookingfuel, tab(9) *** Standard MPI *** /* Members of the household are considered deprived if the household uses solid fuels and solid biomass fuels for cooking. */ ***************************************************************** gen cooking_mdg = 1 replace cooking_mdg = 0 if cookingfuel==5 replace cooking_mdg = . if cookingfuel==. lab var cooking_mdg "Household has cooking fuel by MDG standards" /* Non deprived if: gas from cylendre; gas; kaz/ kerosene; other Deprived if: coal; wood */ tab cookingfuel cooking_mdg, miss *** Destitution MPI *** *** (same as standard MPI) *** **************************************** gen cooking_u = cooking_mdg lab var cooking_u "Household uses clean fuels for cooking" ******************************************************************************** *** Step 2.10 Assets ownership *** ******************************************************************************** //Check that for standard assets in living standards: "no"==0 and yes=="1" codebook h525_2 h525_1 h525_10 h525_11 h525_5 h529_2 h529_1 h525_14 h529_12 replace h525_1 = 0 if h525_1 == 2 replace h525_2 = 0 if h525_2 == 2 replace h525_10 = 0 if h525_10 == 2 replace h525_11 = 0 if h525_11 == 2 replace h525_5 = 0 if h525_5 == 2 replace h529_2 = 0 if h529_2 == 2 replace h529_1 = 0 if h529_1 == 2 replace h525_14 = 0 if h525_14 == 2 replace h529_12 = 0 if h529_12 == 2 label define yesno 0 "no" 1 "yes" label values h525_1 yesno label values h525_2 yesno label values h525_10 yesno label values h525_11 yesno label values h525_5 yesno label values h529_2 yesno label values h529_1 yesno label values h525_14 yesno label values h529_12 yesno codebook h525_2 h525_1 h525_10 h525_11 h525_5 h529_2 h529_1 h525_14 h528 h529_12 clonevar television = h525_2 gen bw_television = . clonevar radio = h525_1 clonevar telephone = h525_10 clonevar mobiletelephone = h525_11 clonevar refrigerator = h525_5 clonevar car = h529_2 clonevar bicycle = h529_1 gen motorbike = . //Syria PAPFAM 2009 has no data on ownership of motorcycle clonevar computer = h525_14 lookfor internet /*Syria PAPFAM 2009 has no data on internet ownership but has data on number of internet users within household. Number of internet users was used as a proxy for ownership of internet*/ clonevar internet = h528 replace internet = 1 if internet>0 & internet <99 replace internet = . if internet==99 replace internet = 2 if internet==0 label define internet 1"Yes" 2"No", modify label values internet internet tab internet, miss gen animal_cart = . //Syria PAPFAM 2009 has no observation for animal cart lookfor land /*Syria PAPFAM 2009 has two separate variables for land - farm and non-farming. To be consistent with the DHS and MICS surveys, the focus is only on farm land, that is land used for agricultural purposes. */ clonevar land = h529_12 //Ensure the category "no" is recoded as "0" tab land, miss //Syria PAPFAM 2009 has no data on landsize gen landsize = . gen landsizemin=. label var landsizemin "owns 3 ha or more land" foreach var in television radio telephone mobiletelephone refrigerator /// car bicycle motorbike computer internet land /// animal_cart { replace `var' = 0 if `var'==2 replace `var' = . if `var'==9 | `var'==99 | `var'==8 | `var'==98 } //Missing values replaced //Group telephone and mobiletelephone as a single variable replace telephone=1 if telephone==0 & mobiletelephone==1 replace telephone=1 if telephone==. & mobiletelephone==1 //Label indicators lab var television "Household has television" lab var radio "Household has radio" lab var telephone "Household has telephone (landline/mobilephone)" lab var refrigerator "Household has refrigerator" lab var car "Household has car" lab var bicycle "Household has bicycle" lab var motorbike "Household has motorbike" lab var computer "Household has computer" lab var animal_cart "Household has animal cart" *** Standard MPI *** /* Members of the household are considered deprived in assets if the household does not own more than one of: radio, TV, telephone, bike, motorbike, refrigerator, computer or animal cart and does not own a car or truck.*/ ***************************************************************************** egen n_small_assets2 = rowtotal(television radio telephone refrigerator bicycle motorbike computer animal_cart), missing lab var n_small_assets2 "Household Number of Small Assets Owned" gen hh_assets2 = (car==1 | n_small_assets2 > 1) replace hh_assets2 = . if car==. & n_small_assets2==. lab var hh_assets2 "Household Asset Ownership: HH has car or more than 1 small assets incl computer & animal cart" *** Destitution MPI *** /* Members of the household are considered deprived in assets if the household does not own any assets.*/ ***************************************************************************** gen hh_assets2_u = (car==1 | n_small_assets2>0) replace hh_assets2_u = . if car==. & n_small_assets2==. lab var hh_assets2_u "Household Asset Ownership: HH has car or at least 1 small assets incl computer & animal cart" ******************************************************************************** *** Step 2.11 Rename and keep variables for MPI calculation ******************************************************************************** //Retain data on sampling design: clonevar strata = gov clonevar psu = cluster codebook strata psu //Retain year, month & date of interview: clonevar year_interview = datey clonevar month_interview = datem clonevar date_interview = intc *** Rename key global MPI indicators for estimation *** /* Note: In the case of Syria PAPFAM 2009, there is no birth history file. We are not able to identify whether child mortality occured in the last 5 years preceeding the survey date. As such, for the estimation, we use the indicator 'hh_mortality' that represent all child mortality that was ever reported. */ recode hh_mortality (0=1)(1=0) , gen(d_cm) recode hh_nutrition_uw_st (0=1)(1=0) , gen(d_nutr) recode hh_child_atten (0=1)(1=0) , gen(d_satt) recode hh_years_edu6 (0=1)(1=0) , gen(d_educ) recode electricity (0=1)(1=0) , gen(d_elct) recode water_mdg (0=1)(1=0) , gen(d_wtr) recode toilet_mdg (0=1)(1=0) , gen(d_sani) recode housing_1 (0=1)(1=0) , gen(d_hsg) recode cooking_mdg (0=1)(1=0) , gen(d_ckfl) recode hh_assets2 (0=1)(1=0) , gen(d_asst) *** Rename key global MPI indicators for destitution estimation *** recode hh_mortality_u (0=1)(1=0) , gen(dst_cm) recode hh_nutrition_uw_st_u (0=1)(1=0) , gen(dst_nutr) recode hh_child_atten_u (0=1)(1=0) , gen(dst_satt) recode hh_years_edu_u (0=1)(1=0) , gen(dst_educ) recode electricity_u (0=1)(1=0) , gen(dst_elct) recode water_u (0=1)(1=0) , gen(dst_wtr) recode toilet_u (0=1)(1=0) , gen(dst_sani) recode housing_u (0=1)(1=0) , gen(dst_hsg) recode cooking_u (0=1)(1=0) , gen(dst_ckfl) recode hh_assets2_u (0=1)(1=0) , gen(dst_asst) *** Total number of missing values for each variable *** mdesc psu strata area age /// d_cm d_nutr d_satt d_educ d_elct d_wtr d_sani d_hsg d_ckfl d_asst *** Keep main variables require for MPI calculation *** keep hh_id ind_id psu strata subsample weight area region agec4 agec2 /// d_cm d_nutr d_satt d_educ d_elct d_wtr d_sani d_hsg d_ckfl d_asst order hh_id ind_id psu strata subsample weight area region agec4 agec2 /// d_cm d_nutr d_satt d_educ d_elct d_wtr d_sani d_hsg d_ckfl d_asst *** Generate coutry and survey details for estimation *** char _dta[cty] "Syria" char _dta[ccty] "SYR" char _dta[year] "2009" char _dta[survey] "PAPFAM" char _dta[ccnum] "760" char _dta[type] "micro" char _dta[class] "old_survey" *** Sort, compress and save data for estimation *** sort ind_id compress la da "Micro data for `_dta[ccty]' (`_dta[ccnum]') from `c(current_date)' (`c(current_time)')." save "$path_out/syr_papfam09.dta", replace