The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Topics » Wealth Index » Constructing Wealth Incidence Curve (Confusion- Using Wealth Index to Construct Wealth Incidence Curve)
Re: Constructing Wealth Incidence Curve [message #29353 is a reply to message #29351] Tue, 04 June 2024 08:42 Go to previous messageGo to previous message
tednoel is currently offline  tednoel
Messages: 12
Registered: April 2024
Member
Ok I now understand my original error. I had believed that there were numeric values associated with the "poorest" "poor" etc, categories for the wealth index. At this point I have shifted the strategy to creating wealth indices that demonstrate the change in proportion of those who are categorized as "poorest" from 2015-2021 for Senegal. I can generate the graph but there are no indices and I suspect this is an issue of scale. I know I need to incorporate survey weights in the code but I'm not sure where to do this.

*** Disaggregate Wealth Quintiles by Year
svyset [pweight=hv005], strata(hv023) psu(hv001)
gen wealth_quintile_2015 = hv270 if hv007 == 2015
gen wealth_quintile_2017 = hv270 if hv007 == 2017
gen wealth_quintile_2018 = hv270 if hv007 == 2018
gen wealth_quintile_2019 = hv270 if hv007 == 2019
gen wealth_quintile_2020 = hv270 if hv007 == 2020
gen wealth_quintile_2021 = hv270 if hv007 == 2021

****** MY FIRST ATTEMPT
*Calculate the proportion of individuals classified as "poorest" within the variable for each year
egen poverty2015 = group(wealth_quintile_2015)
egen prop_poorest_2015 = mean(wealth_quintile_2015 == 1), by(poverty2015)

egen poverty2017 = group(wealth_quintile_2017)
egen prop_poorest_2017 = mean(wealth_quintile_2017 == 1), by(poverty2017)

egen poverty2018 = group(wealth_quintile_2018)
egen prop_poorest_2018 = mean(wealth_quintile_2018 == 1), by(poverty2018)

egen poverty2019 = group(wealth_quintile_2019)
egen prop_poorest_2019 = mean(wealth_quintile_2015 == 1), by(poverty2015)

egen poverty2020 = group(wealth_quintile_2020)
egen prop_poorest_2020 = mean(wealth_quintile_2020 == 1), by(poverty2020)

egen poverty2021 = group(wealth_quintile_2021)
egen prop_poorest_2021 = mean(wealth_quintile_2021 == 1), by(poverty2021)

*** This did not work so below is my second attempt
*******

* Calculate the proportion of individuals classified as "poorest" for each year separately
egen prop_poorest_2015_2 = mean(wealth_quintile_2015 == 1 & wealth_quintile_2015)
egen prop_poorest_2017_1 = mean(wealth_quintile_2017 == 1 & wealth_quintile_2017)
egen prop_poorest_2018_1 = mean(wealth_quintile_2018 == 1 & wealth_quintile_2018)
egen prop_poorest_2019_1 = mean(wealth_quintile_2019 == 1 & wealth_quintile_2019)
egen prop_poorest_2020_1 = mean(wealth_quintile_2020 == 1 & wealth_quintile_2020)
egen prop_poorest_2021_1 = mean(wealth_quintile_2021 == 1 & wealth_quintile_2021)


twoway (line prop_poorest_2015_1 prop_poorest_2017_1 prop_poorest_2018_1 prop_poorest_2019_1 prop_poorest_2020_1 prop_poorest_2021_1, ytitle(Proportion "Poorest") xtitle(Years) legend(off))

Neither one of these attempts produces a viable graph and I think I need to incorporate the sampling weights
 
Read Message
Read Message
Read Message
Previous Topic: Comparative wealth index for Uganda (2006 & 2011)
Next Topic: Antenatal care
Goto Forum:
  


Current Time: Sat Nov 23 13:42:50 Coordinated Universal Time 2024