unzipfile "output7717034280080927434.zip", replace
use ESS1-6e01_0_F1, clear
*use cntry essround stfjbot wkhct using ESS1-6e01_0_F1, clear
// Select Round 5
keep if inlist(essround, 5)
// Fix country variable
encode cntry, gen(country)
// Satisfaction with WLB
rename stfjbot swlb
// Contracted working hours
recode wkhct ( 0/20 = 1 "Marginal part-time") ///
( 21/34.75 = 2 "Substantial part-time") ///
(35/220 = 0 "Full-time") ///
, gen(parttime)
label var parttime "Working hours"
// Professional status
generate professional = .
replace professional = (iscoco >= 1000 & iscoco <= 2999)
label var professional "Professional status"
label define professional 0 "Non-prof." 1 "Professional"
label val professional professional
// Sex
recode gndr ( 2 = 1 "Female") ///
( 1 = 0 "Male") ///
(.a = .a "No answer") ///
, gen(female) label(female)
label var female "Gender"
// Select cases: Those in paid work and with a partner
keep if partner == 1 & mnactic == 1
// Are part-time workers more satisfied with their WLB than full-time employees?
qui regress swlb i.parttime i.country, cluster(country)
estimates store m1
// Hypothesis 1b: The effect is stronger for marginal PT than substantial PT.
qui test 1.parttime == 2.parttime
local F = round(r(F), .001)
local p = round(r(p), .001)
local r = r(df)
qui margins, at(parttime=(0 1 2))
marginsplot, recast(bar) ///
title("SWLB difference between full-time and part-time workers", size(large)) ///
plotopts(fcolor(gs14) lcolor(black)) ytitle("Predicted SWLB") ylabel(6 (.5) 7.5, format(%6.1f) grid) ///
name(noint, replace) xtitle("") ysize(3) ///
note(" " ///
"Marginal and substantial part-time" ///
"differs significantly:" ///
"{it:F}(`r', `r(df_r)') = `F', {it:p} = `p'", ///
pos(11) ring(0) bmargin(small)) ///
nodraw
// Are professional part-time workers less satisfied than non-professional part-time workers?
qui regress swlb i.parttime##i.professional i.country, cluster(country)
estimates store m2
qui margins, at(parttime=(0 1 2) professional=(0 1))
marginsplot, recast(bar) xdimension(professional) ///
bydimension(parttime) byopts(row(1) noiyaxes imargin(zero) ///
title("Interaction part-time status and professional status")) ///
subtitle(, pos(6)) /// // Place label of by dimensions below plot
plotopts(fcolor(gs14) lcolor(black)) ytitle("Predicted SWLB") ylabel(6 (.5) 7.5, format(%6.1f)) ///
name(int1, replace) xtitle("") ysize(3) nodraw
// Are women working part-time more satisfied with their SWLB than part-time working men?
qui regress swlb i.parttime##i.female i.country, cluster(country)
estimates store m3
qui margins, at(parttime=(0 1 2) female=(0 1))
marginsplot, recast(bar) xdimension(female) ///
bydimension(parttime) byopts(row(1) noiyaxes imargin(zero) ///
title("Interaction part-time status and gender")) ///
subtitle(, pos(6)) /// // Place label of by dimensions below plot
plotopts(fcolor(gs14) lcolor(black)) ytitle("Predicted SWLB") ylabel(6 (.5) 7.5, format(%6.1f)) ///
name(int2, replace) xtitle("") ysize(3) nodraw
// Output table and figure
esttab m1 m2 m3 using test.tex, compress replace se label nomtitles ///
indicate(Country dummies = *country) ///
varwidth(30) interaction(" X ") ///
title(Regression table\label{tab1}) ///
booktabs
graph combine noint int1 int2, col(1) ysize(9) ///
note("95% CI's based on cluster-robust standard errors")
