clear all *** Import Data *** use "/Users/Nicola/Box Sync/2_Euro_Melle/Aggregate stuff/Output/output_aggregate1.dta" xtset country_name year cd "/Users/Nicola/Box Sync/2_Euro_Melle/JEPP Submission/TABLES_and_GRAPHS" * Labeling variable as asdoc does not allow to change variable names inside gen net_support_lag = L.net_support label variable net_support_lag "Net Support (Lag)" label variable net_support "Net Support" label variable fisconso_dummy_l "Austerity (Lag)" label variable EB_dummy_l "EB Austerity (Lag)" label variable TB_dummy_l "TB Austerity (Lag)" label variable exchange_rate_l "Exchange rate (Lag)" label variable unempl_l "Unemployment Rate (Lag)" label variable inflation_l "Inflation Rate (Lag)" label variable gdp_pc_growth_l "GDP p.c. growth (Lag)" *** Table 1 (fisconso) * With Lagged dependent variable * Null effects asdoc xtreg net_support net_support_lag fisconso_dummy_l if eurozone == 1, fe vce(cluster country) nest save(Regression1) replace fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag fisconso_dummy_l _cons) title(Table 1) label add(Lagged controls , No, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support net_support_lag fisconso_dummy_l i.year if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag fisconso_dummy_l _cons) title(Table 1) label add(Lagged controls , No, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) asdoc xtreg net_support net_support_lag fisconso_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag fisconso_dummy_l _cons) title(Table 1) label add(Lagged controls , Yes, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support i.year net_support_lag fisconso_dummy_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag fisconso_dummy_l _cons) title(Table 1) label add(Lagged controls , Yes, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) *** Table 2 (EB and TB) * With lagged dependent variable. * Negative at 10% and null effects asdoc xtreg net_support net_support_lag EB_dummy_l if eurozone == 1, fe vce(cluster country) nest reset save(Regression2) replace fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support net_support_lag EB_dummy_l i.year if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) asdoc xtreg net_support net_support_lag EB_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l _cons) label add(Lagged controls , Yes, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support i.year net_support_lag EB_dummy_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l _cons) title(Table 2) label add(Lagged controls , Yes, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) * Null effects asdoc xtreg net_support net_support_lag TB_dummy_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l TB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support net_support_lag TB_dummy_l i.year if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l TB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) asdoc xtreg net_support net_support_lag TB_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l TB_dummy_l _cons) label add(Lagged controls , Yes, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support i.year net_support_lag TB_dummy_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(net_support_lag EB_dummy_l TB_dummy_l _cons) title(Table 2) label add(Lagged controls , Yes, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) ********************************************* APPENDIX *************************************************************************** *** table 1A *** Without Lagged dependent variable asdoc xtreg net_support fisconso_dummy_l if eurozone == 1, fe vce(cluster country), nest reset save(Regression_Appendix1a) replace fs(10) stat(N, r2_a) font(Times New Roman) keep(fisconso_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support fisconso_dummy_l i.year if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(fisconso_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) asdoc xtreg net_support fisconso_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(fisconso_dummy_l _cons) label add(Lagged controls , Yes, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support i.year fisconso_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country), nest fs(10) stat(N, r2_a) font(Times New Roman) keep(fisconso_dummy_l _cons) title(Table 1A) label add(Lagged controls , Yes, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) *** table 2a *** Negative at 10% and null effects asdoc xtreg net_support EB_dummy_l if eurozone == 1, fe vce(cluster country), nest reset save(Regression_Appendix2a) replace fs(10) stat(N, r2_a) font(Times New Roman) keep(EB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support EB_dummy_l i.year if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(EB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) asdoc xtreg net_support EB_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country), nest fs(10) stat(N, r2_a) font(Times New Roman) keep(EB_dummy_l _cons) label add(Lagged controls , Yes, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support i.year EB_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(EB_dummy_l _cons) title(Table 2A) label add(Lagged controls , Yes, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) *** table 3A *** Null effects asdoc xtreg net_support TB_dummy_l if eurozone == 1, fe vce(cluster country), nest reset save(Regression_Appendix3a) replace fs(10) stat(N, r2_a) font(Times New Roman) keep(TB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support TB_dummy_l i.year if eurozone == 1, fe vce(cluster country), nest fs(10) stat(N, r2_a) font(Times New Roman) keep(TB_dummy_l _cons) label add(Lagged controls , No, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) asdoc xtreg net_support TB_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country), nest fs(10) stat(N, r2_a) font(Times New Roman) keep(TB_dummy_l _cons) label add(Lagged controls , Yes, Country FE, Yes, Year FE, No) cnames(Net Support) drop(_cons) asdoc xtreg net_support i.year TB_dummy_l exchange_rate_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, fe vce(cluster country) nest fs(10) stat(N, r2_a) font(Times New Roman) keep(TB_dummy_l _cons) title(Table 3A) label add(Lagged controls , Yes, Country FE, Yes, Year FE, Yes) cnames(Net Support) drop(_cons) ************* de Chaisemartin, C and D'Haultfoeuille, X (2020a). American Economic Review, vol. 110, no. 9. Two-Way Fixed Effects Estimators with Heterogeneous Treatment Effects. ************ Null effects did_multiplegt net_support country_name year fisconso_dummy_l if eurozone == 1, cluster(country_name) breps(150) graphoptions(nodraw) did_multiplegt net_support country_name year fisconso_dummy_l if eurozone == 1, controls(exchange_rate_l unempl_l inflation_l gdp_pc_growth_l) cluster(country_name) breps(150) graphoptions(nodraw) did_multiplegt net_support country_name year EB_dummy_l if eurozone == 1, cluster(country_name) breps(150) graphoptions(nodraw) did_multiplegt net_support country_name year EB_dummy_l if eurozone == 1, controls(exchange_rate_l unempl_l inflation_l gdp_pc_growth_l) cluster(country_name) breps(150) graphoptions(nodraw) did_multiplegt net_support country_name year TB_dummy_l if eurozone == 1, cluster(country_name) breps(150) graphoptions(nodraw) did_multiplegt net_support country_name year TB_dummy_l if eurozone == 1, controls(exchange_rate_l unempl_l inflation_l gdp_pc_growth_l) cluster(country_name) breps(150) graphoptions(nodraw) **** Dynamic fixed effects models for short T. Estimator proposed in Bun, M. J. G., and M. A. Carree. 2005. Bias-corrected estimation in dynamic panel data models. Journal of Business & Economic Statistics 23: 200-210. * table 5A xtdpdbc net_support fisconso_dummy_l if eurozone == 1, vce(robust) small xtdpdbc net_support fisconso_dummy_l if eurozone == 1, vce(robust) teffects small xtdpdbc net_support fisconso_dummy_l unempl_l inflation_l gdp_pc_growth_l exchange_rate_l if eurozone == 1, vce(robust) small xtdpdbc net_support fisconso_dummy_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, vce(robust) teffects small * table 6A xtdpdbc net_support EB_dummy_l if eurozone == 1, vce(robust) small xtdpdbc net_support EB_dummy_l if eurozone == 1, vce(robust) teffects small xtdpdbc net_support EB_dummy_l unempl_l inflation_l gdp_pc_growth_l exchange_rate_l if eurozone == 1, vce(robust) small xtdpdbc net_support EB_dummy_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, vce(robust) teffects small * table 7A xtdpdbc net_support TB_dummy_l if eurozone == 1, vce(robust) small xtdpdbc net_support TB_dummy_l if eurozone == 1, vce(robust) teffects small xtdpdbc net_support TB_dummy_l unempl_l inflation_l gdp_pc_growth_l exchange_rate_l if eurozone == 1, vce(robust) small xtdpdbc net_support TB_dummy_l unempl_l inflation_l gdp_pc_growth_l if eurozone == 1, vce(robust) teffects small ************************************************************************************************************************************* ************ Import & Clean Data *** use "/Users/Nicola/Box Sync/2_Euro_Melle/data from Melle/harmonised_EB_2004-2021_v3-0-0.dta", clear save working.dta, replace cd "/Users/Nicola/Box Sync/2_Euro_Melle/JEPP Submission/TABLES_and_GRAPHS" ************ Control variables chosen to coincide with Banducci et al. 2009 & Roth et al. 2016 *** recode curr 2=0 // 1= for, 0= against single currency gen unemp=. // 1 = unemployment, 0 = employed or retired recode unemp .=1 if occup==540 recode unemp .=0 gen retired=. // 1 = retired, 0 = employed or unemployed recode retired .=1 if occup==530 recode retired .=0 recode p_28 2=0 // perception of inflation 1=very important, 0=not at all recode citizen_h 2=0 3=0 // EU identity, 1 = often feels like EU citizen gen goodmem = mem // membership a good thing 1 = good, 0 = neither or bad recode goodmem 2=0 3=0 gen badmem = mem // membership a bad thing 1 = bad, 0 = neither or good recode badmem 1=0 2=0 3=1 recode educ -99951 = . // level of education recode gender 2=0 3=. // gender male=1, female=0 gen age1 = age/10 // age in decades // lr = left-right self-placement > lower values are left, higher right-wing attitudes save working.dta, replace ************ Merge with Country-Year level Data for macroeconomic variables*** rename country country_encode gen country = "missing" replace country = "AUT" if country_encode==18 replace country = "BEL" if country_encode==2 replace country = "BG" if country_encode==29 replace country = "CYP" if country_encode==19 replace country = "CZ" if country_encode==20 replace country = "DEU" if country_encode==4 replace country = "DNK" if country_encode==7 replace country = "EE" if country_encode==21 replace country = "ESP" if country_encode==12 replace country = "FIN" if country_encode==16 replace country = "FRA" if country_encode==1 replace country = "GBR" if country_encode==9 replace country = "GRC" if country_encode==11 replace country = "HR" if country_encode==32 replace country = "HU" if country_encode==22 replace country = "IRL" if country_encode==8 replace country = "ITA" if country_encode==5 replace country = "LT" if country_encode==24 replace country = "LV" if country_encode==23 replace country = "LUX" if country_encode==6 replace country = "PL" if country_encode==26 replace country = "PRT" if country_encode==13 replace country = "RO" if country_encode==30 replace country = "SI" if country_encode==28 replace country = "SK" if country_encode==27 replace country = "SWE" if country_encode==17 ************ * Nicola's adding: Northern Ireland as GBR, netherkabds, and East Germany replace country = "GBR" if country_encode==10 replace country = "NDL" if country_encode==3 replace country = "DEU" if country_encode==14 ************ merge m:1 year country using "/Users/Nicola/Box Sync/2_Euro_Melle/Aggregate stuff/Output/output_aggregate1.dta" keep if _merge==3 drop _merge save working.dta, replace * Center ideology so that it is symmetrical around the middle value gen lr2 = lr - 5.5 if lr != . * Keep only variables of interest keep ebid intmonth curr unemp retired p_28 educ gender age1 lr2 inflation_l gdp_pc_growth_l eurozone fisconso_dummy_l EB_dummy_l TB_dummy_l year unempl_l citizen_h goodmem badmem country_encode country label variable fisconso_dummy_l "Austerity (Lag)" label variable EB_dummy_l "EB Austerity (Lag)" * Misssingness patterns by year. * 2019: 50% missingness in goodmem, badmem. Notice that all EB models wille xclude 2019 as well because there is no eurozone country that did EB fisconsolidation that year. * 2014-2018: 100% missingness in goodmem and badmem. * 2012-2013: no currency question; in the Spring one there is also no ideology question. * 2012: 100% missingness in goodmem, badmem * 2004: 100% missingness in perception of inflation * To avoid too many NAs, main analysis does no tinlcude: goodmem, badmem. **************************************** Baseline results with no interaction*********************************** * Relative to our previous submission: we include unemplyment aggregate. ************ Table 3 and Figures 1 and 2 * Model 1 eststo clear * fisconso_dummy_l: -.0764791 beta / .0135098 SE / -5.66 t stat/ 0.000 p value loc karibou = "unemp retired educ gender age1 lr2 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" eststo: melogit curr `karibou' i.fisconso_dummy_l `nat' i.year if eurozone == 1 || country_encode: , startvalues(zero) * delta probability = .6349596 - .6291376 margins i.fisconso_dummy_l, atmeans marginsplot, /// plot1(m(D) mc(gs1) mlc(gs1) lc(gs1%0)) /// ci1(lc(gs1)) /// title("") /// ylab(, angle(0)) graphr(c(white)fc(white)lc(white)) /// plotr(c(white)fc(white)lc(white)) /// scheme(s2mono) /// xlab(0 "No" 1 "Yes") xti(Lagged Fiscal Consolidation) /// yscale(range(0.625(0.002)0.650)) ylab(0.625(0.002)0.650) yti(Marginal Probability) /// saving(margins_fisconso1.gph, replace) estadd local Individual_controls "Yes" estadd local Country_controls_lag "Yes" estadd local Country_RE "Yes" estadd local Year_FE "Yes" * Model 2 * EB -.0756886 .0161415 -4.69 0.000 loc karibou = "unemp retired educ gender age1 lr2 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" eststo: melogit curr `karibou' i.EB_dummy_l `nat' i.year if eurozone == 1 & TB_dummy_l == 0 || country_encode: , startvalues(zero) margins i.EB_dummy_l, atmeans marginsplot, /// plot1(m(D) mc(gs1) mlc(gs1) lc(gs1%0)) /// ci1(lc(gs1)) /// title("") /// ylab(, angle(0)) graphr(c(white)fc(white)lc(white)) /// plotr(c(white)fc(white)lc(white)) /// scheme(s2mono) /// xlab(0 "No" 1 "Yes") xti(Lagged Expediture Based Consolidation) /// yscale(range(0.625(0.002)0.650)) ylab(0.625(0.002)0.650) yti(Marginal Probability) /// saving(margins_EB1.gph, replace) estadd local Individual_controls "Yes" estadd local Country_controls_lag "Yes" estadd local Country_RE "Yes" estadd local Year_FE "Yes" esttab using table1_mlm.rtf, se(3) stats(Individual_controls Country_controls_lag Country_RE Year_FE N N_clust ll, labels("Individual-level controls" "Country-level controls (Lag)" "Country RE" "Year FE" "N Individuals" "N Country" "Log Likelihood")) /// varwidth(28) modelwidth(17) label nobaselevels star(* 0.10 ** 0.05 *** 0.001) /// keep(1.fisconso_dummy_l 1.EB_dummy_l _cons) /// title(This is a regression table) eqlabel(none) /// mtitles("Fiscal Austerity" "EB Austerity") eststo clear /// Concern for inflation and fiscal consolidation notice constantonly rather than zero as maximization at times * .1945687 (0.000) / -.0151286 not SS / -.0347234 not SS loc karibou = "unemp retired educ gender age1 lr2" loc nat = "inflation_l gdp_pc_growth_l unempl_l" eststo: melogit curr `karibou' i.p_28##i.fisconso_dummy_l `nat' i.year if eurozone == 1 || country_encode: , startvalues(constantonly) estadd local Individual_controls "Yes" estadd local Country_controls_lag "Yes" estadd local Country_RE "Yes" estadd local Year_FE "Yes" /// Gender and fiscal consolidation (Melle) * With unemployment aggregate loc karibou = "unemp retired educ age1 lr2 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" eststo: melogit curr `karibou' i.gender##i.fisconso_dummy_l `nat' i.year if eurozone == 1 || country_encode: , startvalues(zero) estadd local Individual_controls "Yes" estadd local Country_controls_lag "Yes" estadd local Country_RE "Yes" estadd local Year_FE "Yes" /// Unemployment and fiscal consolidation * unemp -.3084749 (0.001) fisconso -.0338434 (0.015) interaction -.1127516 (0.007) loc karibou = "retired educ gender age1 lr2 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" eststo: melogit curr `karibou' i.unemp##i.fisconso_dummy_l `nat' i.year if eurozone == 1 || country_encode: , startvalues(zero) estadd local Individual_controls "Yes" estadd local Country_controls_lag "Yes" estadd local Country_RE "Yes" estadd local Year_FE "Yes" esttab esttab, se(3) stats(Individual_controls Country_controls_lag Country_RE Year_FE N N_clust ll, labels("Individual-level controls" "Country-level controls (Lag)" "Country RE" "Year FE" "N Individuals" "N Country" "Log Likelihood")) /// varwidth(28) modelwidth(17) label nobaselevels star(* 0.10 ** 0.05 *** 0.001) /// drop(inflation_l gdp_pc_growth_l unempl_l retired educ age1 lr2) /// title(This is a regression table) eqlabel(none) /// mtitles("Fiscal Austerity" "EB Austerity") eststo clear **** Only significant interaction is unemployment. Since it is the last estimation we can call margins directly margins fisconso_dummy_l, at(unemp=(0 1)) marginsplot, recast(scatter) /// plot1(m(D) mc(gs1) mlc(gs1) lc(gs1%0)) /// ci1(lc(gs1)) /// title("") /// ylab(, angle(0)) graphr(c(white)fc(white)lc(white)) /// plotr(c(white)fc(white)lc(white)) /// scheme(s2mono) /// xlab(0 "No" 1 "Yes") xti(Lagged Fiscal Consolidation) /// yscale(range(0.58(0.004)0.64)) ylab(0.58(0.004)0.64) yti(Predicted Probability) /// saving(margins_unempXfisconso1.gph, replace) *** Now same a sabove but with EB /// Concern for inflation and fiscal consolidation --> chenge lr2 --> notice constantonly rather than zero as maximization * Next inlcude unemployment aggregate .2044272 (0.001) / -.1126836 (0.001) / -.0272462 (notSS) loc karibou = "unemp retired educ gender age1 lr2" loc nat = "inflation_l gdp_pc_growth_l unempl_l" melogit curr `karibou' i.p_28##i.EB_dummy_l `nat' i.year if eurozone == 1 & TB_dummy_l == 0 || country_encode: , startvalues(zero) /// Unemployment and fiscal consolidation * * * unemp -.3268559 0.000 / EB_dummy_l -.1307627 0.000 / interaction -.0505722 (notSS) loc karibou = "retired educ gender age1 lr2 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" melogit curr `karibou' i.unemp##i.EB_dummy_l `nat' i.year if eurozone == 1 & TB_dummy_l == 0 || country_encode: , startvalues(zero) /// Gender and fiscal consolidation loc karibou = "unemp retired educ age1 lr2 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" melogit curr `karibou' i.gender##i.EB_dummy_l `nat' i.year if eurozone == 1 & TB_dummy_l == 0 || country_encode: , startvalues(zero) ************************************ Triple interaction*********************************** loc karibou = "unemp retired educ gender age1 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" melogit curr `karibou' c.lr2##c.lr2##i.fisconso_dummy_l `nat' i.year if eurozone == 1 || country_encode: , startvalues(zero) margins fisconso_dummy_l, at (lr2=(-4.5(1)4.5) lr2=(-4.5(1)4.5)) marginsplot, /// saving(fisconsoXlrXlr.gph, replace) *************************************** Now EB_dummy_l *********************************** loc karibou = "unemp retired educ gender age1 p_28" loc nat = "inflation_l gdp_pc_growth_l unempl_l" melogit curr `karibou' c.lr2##c.lr2##i.EB_dummy_l `nat' i.year if eurozone == 1 & TB_dummy_l == 0 || country_encode: , startvalues(constantonly) margins EB_dummy_l, at (lr2=(-4.5(1)4.5) lr2=(-4.5(1)4.5)) marginsplot, /// saving(EBXlrXlr.gph, replace)