Re: Constructing Wealth Incidence Curve [message #29353 is a reply to message #29351] |
Tue, 04 June 2024 08:42 |
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
|
|
|