gen 	married_women = 0
replace married_women = 1 if worker==2

gen 	treated3=0
replace treated3=1 if experimental==1 & married_women==1

reg earnings experimental##married_women year##experimental year##married_women treated3##ib2014.year, cluster(state)

coefplot, keep(1.treated3#*) omitted baselevels cirecast(rcap) ///
    rename(1.treated3#([0-9]+).year = \1, regex) at(_coef) ///
    yline(0, lp(solid)) yline(-5000, lp(dot)) xline(2014.5, lpattern(dash)) ///
    xlab(2010(1)2019)