clear
// Install rcall command
*github install haghish/rcall, stable
// Clear memory of R session
rcall clear
//#install.packages("essurvey");
rcall: library(essurvey)
// Download data
rcall: download_rounds(c(9), ///
ess_email = "REGISTERED E-MAIL ADDRESS HERE", ///
output_dir = "data", ///
format = 'stata')
// Open data set
local vars cntry gndr polintr dweight
use `vars' using "data/ESS9/ESS9e01_2.dta", clear
// Variables
generate interest = (inlist(polintr, 1, 2))
replace interest = . if missing(polintr)
generate female = (gndr == 2)
// Calculate gender gap per country
tempname foo
postfile `foo' str50 cntry men women gap gap_se using deleteme.dta, replace
levelsof cntry, local(cntry)
foreach i of local cntry {
qui regress interest i.female [pw = dweight] if cntry == `"`i'"'
qui margins i.female, post
local men = _b[0.female] * 100
local women = _b[1.female] * 100
qui regress interest i.female [pw = dweight] if cntry == `"`i'"'
qui margins r.female, post
local gap = _b[r1vs0.female] * -100
local gap_se = _se[r1vs0.female] * -100
post `foo' ("`i'") (`men') (`women') (`gap') (`gap_se')
}
postclose `foo'
// Open data set with gender gap per country
use deleteme.dta, clear
kountry cntry, from(iso2c)
egen country = rank(gap), unique
labmask country, value(NAMES_STD)
generate ub = gap + 1.96 * gap_se
generate lb = gap - 1.96 * gap_se
generate y1 = -women
generate y2 = 0
generate y3 = -women + men
generate x1 = country - 0.2
generate x2 = country + 0.2
mylabels 0(20)60, myscale(-@) local(show)
local spaces = 30 * " "
twoway (scatter y1 country, msymbol(p) yaxis(1 2)) ///
(rbar y1 y2 x1, blcolor(gs10%40) bfcolor(gs10%40) barwidth(0.4)) ///
(rbar y1 y3 x2, blcolor(gs5%40) bfcolor(gs5%40) barwidth(0.4)) ///
(rbar gap y2 x2, blcolor(gs3%40) bfcolor(gs3%40) barwidth(0.4)) ///
(rspike ub lb x2) ///
, ylabel(`show', grid ang(h)) ///
ylabel(0(10)40, grid ang(h) axis(2)) ///
ytitle("% interested in politics") ///
ytitle("`spaces' Gender gap", axis(2)) ///
xlabel(1/19, noticks ang(45) valuelabel labsize(small)) ///
legend(order(2 "Men" 3 "Women" 4 "Gender gap" 5 "Gender gap 95% CI") ///
col(4) pos(12) ring(0)) name(gendergap, replace)
