gen country = substr(elecstud,1,3) gen outside = (party==0) gen gr_lead = gr * lead gen unem_lead = unem * lead gen gr_loc_md_lead = gr_loc_md * lead gen gr_loc_pc_lead = gr_loc_pc * lead gen gr_loc_tr_lead = gr_loc_tr * lead gen gr_glob_md_lead = gr_glob_md * lead gen gr_glob_pc_lead = gr_glob_pc * lead gen gr_glob_tr_lead = gr_glob_tr * lead gen unem_loc_md_lead = unem_loc_md * lead gen unem_loc_pc_lead = unem_loc_pc * lead gen unem_loc_tr_lead = unem_loc_tr * lead gen unem_glob_md_lead = unem_glob_md * lead gen unem_glob_pc_lead = unem_glob_pc * lead gen unem_glob_tr_lead = unem_glob_tr * lead * basic model clogit votemulti outside lead dist gr_loc_pc_lead gr_glob_pc_lead unem_loc_pc_lead unem_glob_pc_lead, group(id) * cluster by election study -- definitely do at least this! clogit votemulti outside lead dist gr_loc_pc_lead gr_glob_pc_lead unem_loc_pc_lead unem_glob_pc_lead, group(id) cluster(elecstud) * questionable whether we have enough clusters clogit votemulti outside lead dist gr_loc_pc_lead gr_glob_pc_lead unem_loc_pc_lead unem_glob_pc_lead, group(id) cluster(country) * random effcets hard to estimate / fixed effects don't make sense