******************************************************************************** /* Suggested citation: Oxford Poverty and Human Development Initiative (OPHI), University of Oxford. 2019 Global Multidimensional Poverty Index - Philippines DHS 2017 [STATA do-file]. Retrieved from: https://ophi.org.uk/multidimensional-poverty-index/mpi-resources/ For further queries, please contact: ophi@qeh.ox.ac.uk */ ******************************************************************************** clear all set more off set maxvar 10000 set mem 500m *** Working Folder Path *** global path_in "T:/GMPI 2.0/rdta/Philippines DHS 2017" global path_out "G:/GMPI 2.0/cdta" global path_ado "T:/GMPI 2.0/ado" ******************************************************************************** *** PHILIPPINES DHS 2017 *** ******************************************************************************** ******************************************************************************** *** Step 1: Data preparation *** Selecting variables from KR, BR, IR, & MR recode & merging with PR recode ******************************************************************************** /*Philippines DHS 2017: Anthropometric data was not collected from chidren, women and men. The variables exist, but it has no observations */ ******************************************************************************** *** Step 1.1 PR - INDIVIDUAL RECODE *** (Children under 5 years) ******************************************************************************** //No anthropometric data for children ******************************************************************************** *** Step 1.2 BR - BIRTH RECODE *** (All females 15-49 years who ever gave birth) ******************************************************************************** /*The purpose of step 1.2 is to identify children under 18 who died in the last 5 years prior to the survey date.*/ use "$path_in/PHBR70FL.dta", clear *** Generate individual unique key variable required for data merging *** v001=cluster number; *** v002=household number; *** v003=respondent's line number gen double ind_id = v001*1000000 + v002*100 + v003 format ind_id %20.0g label var ind_id "Individual ID" desc b3 b7 gen date_death = b3 + b7 //Date of death = date of birth (b3) + age at death (b7) gen mdead_survey = v008 - date_death //Months dead from survey = Date of interview (v008) - date of death gen ydead_survey = mdead_survey/12 //Years dead from survey gen age_death = b7 label var age_death "Age at death in months" tab age_death, miss //Check whether the age is in months codebook b5, tab (10) gen child_died = 1 if b5==0 //1=child dead; 0=child alive replace child_died = 0 if b5==1 replace child_died = . if b5==. label define lab_died 1 "child has died" 0 "child is alive" label values child_died lab_died tab b5 child_died, miss /*NOTE: For each woman, sum the number of children who died and compare to the number of sons/daughters whom they reported have died */ bysort ind_id: egen tot_child_died = sum(child_died) egen tot_child_died_2 = rsum(v206 v207) //v206: sons who have died; v207: daughters who have died compare tot_child_died tot_child_died_2 //In Philipiness DHS 2017, these figures are identical //Identify child under 18 mortality in the last 5 years gen child18_died = child_died replace child18_died=0 if age_death>=216 & age_death<. label values child18_died lab_died tab child18_died, miss bysort ind_id: egen tot_child18_died_5y=sum(child18_died) if ydead_survey<=5 /*Total number of children under 18 who died in the past 5 years prior to the interview date */ replace tot_child18_died_5y=0 if tot_child18_died_5y==. & tot_child_died>=0 & tot_child_died<. /*All children who are alive or who died longer than 5 years from the interview date are replaced as '0'*/ replace tot_child18_died_5y=. if child18_died==1 & ydead_survey==. //Replace as '.' if there is no information on when the child died tab tot_child_died tot_child18_died_5y, miss bysort ind_id: egen childu18_died_per_wom_5y = max(tot_child18_died_5y) lab var childu18_died_per_wom_5y "Total child under 18 death for each women in the last 5 years (birth recode)" //Keep one observation per women bysort ind_id: gen id=1 if _n==1 keep if id==1 drop id duplicates report ind_id gen women_BR = 1 //Identification variable for observations in BR recode //Retain relevant variables keep ind_id women_BR childu18_died_per_wom_5y order ind_id women_BR childu18_died_per_wom_5y sort ind_id save "$path_out/PHL17_BR.dta", replace ******************************************************************************** *** Step 1.3 IR - WOMEN's RECODE *** (Eligible female 15-49 years in the household) ******************************************************************************** /*The purpose of step 1.3 is to identify all deaths that are reported by eligible women.*/ use "$path_in/PHIR70FL.dta", clear *** Generate individual unique key variable required for data merging *** v001=cluster number; *** v002=household number; *** v003=respondent's line number gen double ind_id = v001*1000000 + v002*100 + v003 format ind_id %20.0g label var ind_id "Individual ID" duplicates report ind_id tab v012, miss codebook v201 v206 v207,tab (999) //Women 15-49 years provided information on child mortality gen women_IR=1 //Identification variable for observations in IR recode keep ind_id women_IR v003 v005 v012 v201 v206 v207 order ind_id women_IR v003 v005 v012 v201 v206 v207 sort ind_id save "$path_out/PHL17_IR.dta", replace ******************************************************************************** *** Step 1.4 PR - INDIVIDUAL RECODE *** (Girls 15-19 years in the household) ******************************************************************************** //No anthropometric data ******************************************************************************** *** Step 1.5 MR - MEN'S RECODE ***(Eligible man 15-64 years in the household) ******************************************************************************** //No male recode ******************************************************************************** *** Step 1.6 PR - INDIVIDUAL RECODE *** (Boys 15-19 years in the household) ******************************************************************************** //No anthropometric data ******************************************************************************** *** Step 1.7 PR - HOUSEHOLD MEMBER'S RECODE ******************************************************************************** use "$path_in/PHPR70FL.dta", clear *** Generate a household unique key variable at the household level using: ***hv001=cluster number ***hv002=household number gen double hh_id = hv001*10000 + hv002 format hh_id %20.0g label var hh_id "Household ID" codebook hh_id *** Generate individual unique key variable required for data merging using: *** hv001=cluster number; *** hv002=household number; *** hvidx=respondent's line number. gen double ind_id = hv001*1000000 + hv002*100 + hvidx format ind_id %20.0g label var ind_id "Individual ID" codebook ind_id sort hh_id ind_id ******************************************************************************** *** Step 1.8 DATA MERGING ******************************************************************************** *** Merging BR Recode ***************************************** merge 1:1 ind_id using "$path_out/PHL17_BR.dta" drop _merge erase "$path_out/PHL17_BR.dta" *** Merging IR Recode ***************************************** merge 1:1 ind_id using "$path_out/PHL17_IR.dta" tab women_IR hv117, miss col tab ha65 if hv117==1 & women_IR ==., miss //Total number of eligible women not interviewed drop _merge erase "$path_out/PHL17_IR.dta" /*Check if the number of women in BR recode matches the number of those who provided birth history information in IR recode. */ count if women_BR==1 count if v201!=0 & v201!=. & women_IR==1 /*Check if the number of women in BR and IR recode who provided birth history information matches with the number of eligible women identified by hv117. */ count if hv117==1 count if women_BR==1 | v201==0 count if (women_BR==1 | v201==0) & hv117==1 tab v201 if hv117==1, miss /*Note: Some 2% eligible women did not provide information on their birth history. This will result in missing value for the child mortality indicator that we will construct later */ sort ind_id ******************************************************************************** *** Step 1.9 KEEP ONLY DE JURE HOUSEHOLD MEMBERS *** ******************************************************************************** /*The Global MPI is based on de jure (permanent) household members only. As such, non-usual residents will be excluded from the sample. */ clonevar resident = hv102 codebook resident, tab (10) label var resident "Permanent (de jure) household member" drop if resident!=1 tab resident, miss /*Philippines DHS 2017: 996 (0.83%) individuals who were non-usual residents were dropped from the sample*/ ******************************************************************************** *** Step 1.10 KEEP HOUSEHOLDS SELECTED FOR ANTHROPOMETRIC SUBSAMPLE *** *** if relevant ******************************************************************************** /*In a number of DHS surveys, only a subsample of households were selected for anthropometric measure. In such cases, the Global MPI estimation is based on this subsample. As such, we only retain the relevant subsample of households for analyses. */ /* Philippines DHS 2017: The survey did not collect anthropometric data. Hence no presence of subsample. */ gen subsample=. label var subsample "Households selected as part of nutrition subsample" ******************************************************************************** *** Step 1.11 CONTROL VARIABLES ******************************************************************************** /* Households are identified as having 'no eligible' members if there are no applicable population, that is, children 0-5 years or adult women 15-49 years. 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 15-59 years *** for adult nutrition indicator *********************************************** //Philippines DHS 2017: No anthropometric data gen no_fem_nutri_eligible = . lab var no_fem_nutri_eligible "Household has no eligible women for anthropometric" *** No eligible women 15-49 years *** for child mortality indicator ***************************************** tab hv117, miss tab hv103 hv117, miss /*Only women who slept the night before in the household were eligible for female interview (women's questionnaire). */ gen fem_eligible = (hv117==1) bysort 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 women for an interview lab var no_fem_eligible "Household has no eligible women for interview" drop fem_eligible hh_n_fem_eligible tab no_fem_eligible, miss *** No eligible men 15-54 years *** for adult nutrition indicator (if relevant) *********************************************** //Philippines DHS 2017: No anthropometric data. gen no_male_nutri_eligible = . lab var no_male_nutri_eligible "Household has no eligible men for anthropometric" *** No eligible men 15-54 years *** for child mortality indicator (if relevant) ***************************************** //Philippines DHS 2017: No male recode. gen no_male_eligible = . lab var no_male_eligible "Household has no eligible man" *** No eligible children under 5 *** for child nutrition indicator ***************************************** gen no_child_eligible = . lab var no_child_eligible "Household has no children eligible" *** No eligible women and men *** for adult nutrition indicator *********************************************** gen no_adults_eligible = . lab var no_adults_eligible "Household has no eligible women or men" *** No Eligible Children and Women *** for child and women nutrition indicator *********************************************** gen no_child_fem_eligible = . lab var no_child_fem_eligible "Household has no children or women eligible" *** No Eligible Women, Men or Children *** for nutrition indicator *********************************************** gen no_eligibles = . lab var no_eligibles "Household has no eligible women, men, or children" *** No Eligible Subsample *** for hemoglobin ***************************************** gen no_hem_eligible = . lab var no_hem_eligible "Household has no eligible individuals for hemoglobin measurements" sort hh_id ind_id ******************************************************************************** *** Step 1.12 RENAMING DEMOGRAPHIC VARIABLES *** ******************************************************************************** //Sample weight /*Note: DHS sample weight are calculated to six decimals but are presented in the standard recode files without the decimal point. As such, all DHS weights should be divided by 1,000,000 before applying the weights to calculation or analysis. */ desc hv005 clonevar weight = hv005 replace weight = weight/1000000 label var weight "Sample weight" //Area: urban or rural desc hv025 codebook hv025, tab (5) clonevar area = hv025 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" //Sex of household member codebook hv104 clonevar sex = hv104 label var sex "Sex of household member" //Age of household member codebook hv105, tab (100) clonevar age = hv105 replace age = . if age>=98 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)" //Marital status of household member clonevar marital = hv115 codebook marital, tab (10) recode marital (0=1)(1=2)(8=.) label define lab_mar 1"never married" 2"currently married" 3"widowed" /// 4"divorced" 5"not living together" label values marital lab_mar label var marital "Marital status of household member" tab hv115 marital, miss //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 /* According to Philippines DHS 2017 (pp. 31 of the PDF): "[...] The Philippines has 17 administrative regions, namely the National Capital Region (NCR), Cordillera Administrative Region (CAR), Region I (Ilocos Region), Region II (Cagayan Valley), Region III (Central Luzon), Region IV-A (CALABARZON), MIMAROPA Region, Region V (Bicol Region), Region VI (Western Visayas), Region VII (Central Visayas), Region VIII (Eastern Visayas), Region IX (Zamboanga Peninsula), Region X (Northern Mindanao), Region XI (Davao Region), Region XII (SOCCSKSARGEN), Caraga Region, and the Autonomous Region in Muslim Mindanao (ARMM) [...] The sampling scheme provides data representative of the country as a whole, for urban and rural areas separately, and for each of the country’s administrative regions." */ clonevar region = hv024 codebook region, tab (99) lab var region "Region for subnational decomposition" label define lab_reg /// 1"Ilocos" 2"Cagayan Valley" 3"Central Luzon" /// 4"Calabarzon" 5"Bicol" 6"Western Visayas" /// 7"Central Visayas" 8"Eastern Visayas" 9"Zamboanga Peninsula" /// 10"Northern Mindanao" 11"Davao" 12"Soccsksargen" /// 13"National Capital Region" 14"Cordillera Admin Region" /// 15"ARMM" 16"Caraga" 17"Mimaropa" 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 *** ******************************************************************************** codebook hv108, tab(30) clonevar eduyears = hv108 //Total number of years of education replace eduyears = . if eduyears>30 //Recode any unreasonable years of highest education as missing value replace eduyears = . if eduyears>=age & age>0 /*The variable "eduyears" was replaced with a '.' if total years of education was more than individual's age */ replace eduyears = 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 */ /*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" ******************************************************************************** *** Step 2.2 Child School Attendance *** ******************************************************************************** codebook hv121, tab (99) clonevar attendance = hv121 recode attendance (2=1) label define lab_attend 1 "currently attending" 0 "not currently attending" label values attendance lab_attend label var attendance "Attended school during current school year" tab attendance, miss *** 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) /* In Philipiness, the official school entrance age to primary school is 6 years. So, age range is 6-14 (=6+8) Source: "http://data.uis.unesco.org/?ReportId=163" Also, in Philippines DHS report 2017 (pp. 41 of the PDF), school attendance is measured for Children age 6-11 for primary school. */ /*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!=. /*Generate a variable that captures the number of eligible school aged children who are attending school */ 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) //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 ******************************************************************************** *** Step 2.3 Nutrition *** ******************************************************************************** /*Philippines DHS 2013 has no information on nutrition.*/ gen hh_nutrition_uw_st = . ******************************************************************************** *** Step 2.4 Child Mortality *** ******************************************************************************** codebook v206 v207 /* v206: number of sons who have died v207: number of daughters who have died */ egen temp_f = rowtotal(v206 v207), missing //Total child mortality reported by eligible women replace temp_f = 0 if v201==0 //This line replaces women who have never given birth 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 *** /* Members of the household are considered deprived if women in the household reported mortality among children under 18 in the last 5 years from the survey year. Members of the household is considered non-deprived if eligible women within the household reported (i) no child mortality or (ii) if any child died longer than 5 years from the survey year or (iii) if any child 18 years and older died in the last 5 years. In adddition, members of the household were identified as non-deprived if eligible men within the household reported no child mortality in the absence of information from women. Households that have no eligible women or adult are considered non-deprived. The indicator takes a missing value if there was missing information on reported death from eligible individuals. */ ************************************************************************ tab childu18_died_per_wom_5y, miss /* The 'childu18_died_per_wom_5y' variable was constructed in Step 1.2 using information from individual women who ever gave birth in the BR file. The missing values represent eligible woman who have never ever given birth and so are not present in the BR file. But these 'missing women' may be living in households where there are other women with child mortality information from the BR file. So at this stage, it is important that we aggregate the information that was obtained from the BR file at the household level. This ensures that women who were not present in the BR file is assigned with a value, following the information provided by other women in the household.*/ replace childu18_died_per_wom_5y = 0 if v201==0 /*Assign a value of "0" for: - all eligible women who never ever gave birth */ replace childu18_died_per_wom_5y = 0 if no_fem_eligible==1 /*Assign a value of "0" for: - individuals living in households that have non-eligible women */ bysort hh_id: egen childu18_mortality_5y = sum(childu18_died_per_wom_5y), missing replace childu18_mortality_5y = 0 if childu18_mortality_5y==. & child_mortality==0 /*Replace all households as 0 death if women has missing value and men reported no death in those households */ label var childu18_mortality_5y "Under 18 child mortality within household past 5 years reported by women" tab childu18_mortality_5y, miss gen hh_mortality_u18_5y = (childu18_mortality_5y==0) replace hh_mortality_u18_5y = . if childu18_mortality_5y==. lab var hh_mortality_u18_5y "Household had no under 18 child mortality in the last 5 years" tab hh_mortality_u18_5y, miss ******************************************************************************** *** Step 2.5 Electricity *** ******************************************************************************** *** Standard MPI *** /*Members of the household are considered deprived if the household has no electricity */ *************************************************** clonevar electricity = hv206 codebook electricity, tab (9) label var electricity "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. */ desc hv205 hv225 clonevar toilet = hv205 clonevar shared_toilet = hv225 codebook shared_toilet, tab(99) *** 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*/ ******************************************************************** codebook toilet, tab(99) gen toilet_mdg = ((toilet<23 | toilet==41) & 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<23 | toilet==41) & shared_toilet==1 /*Household is assigned a value of '0' if it uses improved sanitation but shares toilet with other households */ replace toilet_mdg = 0 if toilet == 14 | toilet == 15 /*Household is assigned a value of '0' if it uses non-improved sanitation: "flush to somewhere else" and "flush don't know where" */ replace toilet_mdg = . if toilet==. | toilet==99 //Household is assigned a value of '.' if it has missing information lab var toilet_mdg "Household has improved sanitation" tab toilet toilet_mdg, 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. */ clonevar water = hv201 clonevar timetowater = hv204 clonevar ndwater = hv202 tab hv202 if water==71, miss /*Households using bottled water are only considered to be using improved water when they use water from an improved source for cooking and personal hygiene. This is because the quality of bottled water is not known.*/ *** 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 */ ******************************************************************** codebook water, tab(99) gen water_mdg = 1 if water==11 | water==12 | water==13 | water==14 | /// water==21 | water==31 | water==41 | water==51 | water==71 replace water_mdg = 0 if water==32 | water==42 | water==43 | /// water==61 | water==62 | water==96 /*Deprived if it is "unprotected well", "unprotected spring", "tanker truck" "surface water (river/lake, etc)", "cart with small tank","other" */ codebook timetowater, tab(9999) replace water_mdg = 0 if water_mdg==1 & timetowater >= 30 & timetowater!=. & /// timetowater!=996 & timetowater!=998 & timetowater!=999 /*Deprived if water source is 30 minutes or more from home, roundtrip. Please check the value assigned to 'in premises'. If this is different from 996, add to the condition accordingly */ replace water_mdg = . if water==. | water==99 replace water_mdg = 0 if water==71 & /// (ndwater==32 | ndwater==42 | ndwater==43 | /// ndwater==61 | ndwater==62 | ndwater==96) /*Households using bottled water for drinking are classified as using an improved or unimproved source according to their water source for non-drinking activities */ lab var water_mdg "Household has safe drinking water on premises" tab water water_mdg, miss ******************************************************************************** *** Step 2.8 Housing *** ******************************************************************************** /* Members of the household are considered deprived if the household has a dirt, sand or dung floor */ clonevar floor = hv213 codebook floor, tab(99) gen floor_imp = 1 replace floor_imp = 0 if floor<=12 | floor==96 //Deprived if mud/earth, sand, dung, other replace floor_imp = . if floor==. | floor==99 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 walls made of natural or rudimentary materials. We followed the report's definitions of natural or rudimentary materials. */ clonevar wall = hv214 codebook wall, tab(99) gen wall_imp = 1 replace wall_imp = 0 if wall<=26 | wall==96 /*Deprived if no wall, cane/palms/trunk, mud/dirt, grass/reeds/thatch, pole/bamboo with mud, stone with mud, plywood, cardboard, carton/plastic, uncovered adobe, canvas/tent, unburnt bricks, reused wood, other */ replace wall_imp = . if wall==. | wall==99 lab var wall_imp "Household has wall that it is not of low quality materials" tab wall wall_imp, miss /* Members of the household are considered deprived if the household has roof made of natural or rudimentary materials. We followed the report's definitions of natural and rudimentary materials. */ clonevar roof = hv215 codebook roof, tab(99) gen roof_imp = 1 replace roof_imp = 0 if roof<=24 | roof==96 /*Deprived if no roof, thatch/palm leaf, mud/earth/lump of earth, sod/grass, plastic/polythene sheeting, rustic mat, cardboard, canvas/tent, wood planks/reused wood, unburnt bricks, other */ replace roof_imp = . if roof==. | roof==99 lab var roof_imp "Household has roof that it is not of low quality materials" tab roof roof_imp, miss *** Standard MPI *** /* Members of the household is deprived in housing if the roof, floor OR walls are constructed from low quality materials.*/ ************************************************************** gen housing_1 = 1 replace housing_1 = 0 if floor_imp==0 | wall_imp==0 | roof_imp==0 replace housing_1 = . if floor_imp==. & wall_imp==. & roof_imp==. lab var housing_1 "Household has roof, floor & walls that it is not low quality material" tab housing_1, miss ******************************************************************************** *** 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 cooking combustible clonevar cookingfuel = hv226 *** Standard MPI *** /* Members of the household are considered deprived if the household uses solid fuels and solid biomass fuels for cooking. */ ***************************************************************** codebook cookingfuel, tab(99) gen cooking_mdg = 1 replace cooking_mdg = 0 if cookingfuel>5 & cookingfuel<95 replace cooking_mdg = . if cookingfuel==. lab var cooking_mdg "Household cooks with clean fuels" /* Deprived if: coal/lignite, charcoal, wood, straw/shrubs/grass, agricultural crop, animal dung */ tab cookingfuel cooking_mdg, miss ******************************************************************************** *** Step 2.10 Assets ownership *** ******************************************************************************** /*Assets that are included in the global MPI: Radio, TV, telephone, bicycle, motorbike, refrigerator, car, computer and animal cart */ *** Television/LCD TV/plasma TV/color TV/black & white tv lookfor tv television plasma lcd codebook hv208 clonevar television = hv208 lab var television "Household has television" *** Radio/walkman/stereo/kindle lookfor radio walkman stereo codebook hv207 clonevar radio = hv207 lab var radio "Household has radio" *** Handphone/telephone/iphone/mobilephone/ipod lookfor telephone téléphone mobilephone ipod codebook hv221 hv243a clonevar telephone = hv221 replace telephone=1 if telephone!=1 & hv243a==1 //hv243a=mobilephone. Combine information on telephone and mobilephone. tab hv243a hv221 if telephone==1,miss lab var telephone "Household has telephone (landline/mobilephone)" *** Refrigerator/icebox/fridge lookfor refrigerator réfrigérateur codebook hv209 clonevar refrigerator = hv209 lab var refrigerator "Household has refrigerator" *** Car/van/lorry/truck lookfor car voiture truck van codebook hv212 clonevar car = hv212 lab var car "Household has car" *** Bicycle/cycle rickshaw lookfor bicycle bicyclette codebook hv210 clonevar bicycle = hv210 lab var bicycle "Household has bicycle" *** Motorbike/motorized bike/autorickshaw lookfor motorbike moto codebook hv211 clonevar motorbike = hv211 lab var motorbike "Household has motorbike" *** Computer/laptop/tablet lookfor computer ordinateur laptop ipad tablet codebook hv243e clonevar computer = hv243e lab var computer "Household has computer" *** Animal cart lookfor cart codebook hv243c clonevar animal_cart = hv243c lab var animal_cart "Household has animal cart" foreach var in television radio telephone refrigerator car /// bicycle motorbike computer animal_cart { replace `var' = . if `var'==9 | `var'==99 | `var'==8 | `var'==98 } //Missing values replaced *** 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" ******************************************************************************** *** Step 2.11 Rename and keep variables for MPI calculation ******************************************************************************** //Retain DHS wealth index: desc hv270 clonevar windex=hv270 desc hv271 clonevar windexf=hv271 //Retain data on sampling design: desc hv022 hv021 clonevar strata = hv022 clonevar psu = hv021 //Retain year, month & date of interview: desc hv007 hv006 hv008 clonevar year_interview = hv007 clonevar month_interview = hv006 clonevar date_interview = hv008 //Generate missing variables gen religion_men = . gen ethnic_men = . gen insurance_men = . *** Rename key global MPI indicators for estimation *** /* Note: Philippines DHS 2017 has no anthropometric data. The global MPI for this country dataset will exclude nutrition, hence will only include 9 out of the 10 global MPI indicators in the estimation. */ recode hh_mortality_u18_5y (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) *** Generate coutry and survey details for estimation *** char _dta[cty] "Philippines" char _dta[ccty] "PHL" char _dta[year] "2017" char _dta[survey] "DHS" char _dta[ccnum] "608" char _dta[type] "micro" *** Sort, compress and save data for estimation *** sort ind_id compress la da "Micro data for `_dta[ccty]' (`_dta[ccnum]'). Last save: `c(filedate)'." save "$path_out/phl_dhs17.dta", replace