tabstat country  , by(country) stat(mean)

table est_size3, cont(freq mean country )

* 
by country: sum country if country < 29


***** now compute the correlations
*** for S1
tetrachoric woco_ext KCLIMATE_D if (country==5 | country==16 | country==19 | country==20)
*by country: sum country if (country==5 | country==16 | country==19 | country==20)
by country: tetrachoric woco_ext KCLIMATE_D if (country==5 | country==16 | country==19 | country==20)
*sum KCLIMATE_D KOSICK  KORETEN KOLOMOT
by country: tetrachoric woco_ext KOSICK  if (country==5 | country==16 | country==19 | country==20)
by country: tetrachoric woco_ext KORETEN if (country==5 | country==16 | country==19 | country==20)
by country: tetrachoric woco_ext KOLOMOT if (country==5 | country==16 | country==19 | country==20)


 
*** for S2
tetrachoric union_ext KCLIMATE_D if (country==13 | country==18 | country==27)
*by country: sum country if (country==13 | country==18 | country==27)
by country: tetrachoric union_ext KCLIMATE_D if (country==13 | country==18 | country==27)

by country: tetrachoric union_ext KOSICK  if (country==13 | country==18 | country==27)
by country: tetrachoric union_ext KORETEN if (country==13 | country==18 | country==27)
by country: tetrachoric union_ext KOLOMOT if (country==13 | country==18 | country==27)



***
***Table 1 for 5 country clusters
***

log using mylog_full_MM_vanderberg_A, replace
*****
***** creating subsets for new runs and replicating TABLE 1
*****
* cluster 1_germanic (CUM LU): ..... & (country==5 | country==19 | country==16| country==20)
* First, replicate table 1/CASE A:
set more off
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)

* now cluster 2_French (sans lux): BE, EL, ES, FR, IT, PT: .... & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)

* now cluster 3_Anglo-saxon: IE, UK: .... & (country==7 | country==28)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)

* now cluster 4_SCANDINAVIAN/NORDIC: SE, FI, DN: .... & (country==4 | country==26 | country==27)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)

* now cluster 5_transition (10): BG, CZ, EE, LV, LT, HU, PL, RO, SK, SI, (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)

log close

log using mylog_full_MM_vanderberg_B, replace
*****
***** creating subsets for new runs and replicating TABLE 1
*****
* cluster 1_germanic (cum LU): ..... & (country==5 | country==19 | country==20)
* First, replicate table 1/CASE B:
set more off
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)

* now cluster 2_French (sans lux): BE, EL, ES, FR, IT, PT: .... & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)

* now cluster 3_Anglo-saxon: IE, UK: .... & (country==7 | country==28)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)

* now cluster 4_SCANDINAVIAN/NORDIC: SE, FI, DN: .... & (country==4 | country==26 | country==27)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)

* now cluster 5_transition (10): BG, CZ, EE, LV, LT, HU, PL, RO, SK, SI, (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)

log close


log using mylog_full_MM_vanderberg_C, replace
*****
***** creating subsets for new runs and replicating TABLE 1
*****
* cluster 1_germanic (cum LU): ..... & (country==5 | country==19 | country==20)
* First, replicate table 1/CASE C:
set more off
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)

* now cluster 2_French (sans lux): BE, EL, ES, FR, IT, PT: .... & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)

* now cluster 3_Anglo-saxon: IE, UK: .... & (country==7 | country==28)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)

* now cluster 4_SCANDINAVIAN/NORDIC: SE, FI, DN: .... & (country==4 | country==26 | country==27)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)

* now cluster 5_transition (10): BG, CZ, EE, LV, LT, HU, PL, RO, SK, SI, (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)

log close


log using mylog_full_MM_vanderberg_D, replace
*****
***** creating subsets for new runs and replicating TABLE 1
*****
* cluster 1_germanic (cum LU): ..... & (country==5 | country==19 | country==20)
* First, replicate table 1/CASE D:
set more off
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==5 | country==16| country==19 | country==20) || country: , intmethod(mcaghermite)

* now cluster 2_French (SANS lux): BE, EL, ES, FR, IT, PT: .... & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==1 | country==8 | country==9 |country==10 |country==12 | country==22) || country: , intmethod(mcaghermite)

* now cluster 3_Anglo-saxon: IE, UK: .... & (country==7 | country==28)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==7 | country==28) || country: , intmethod(mcaghermite)

* now cluster 4_SCANDINAVIAN/NORDIC: SE, FI, DN: .... & (country==4 | country==26 | country==27)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==4 | country==26 | country==27) || country: , intmethod(mcaghermite)

* now cluster 5_transition (10): BG, CZ, EE, LV, LT, HU, PL, RO, SK, SI, (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)
meglm KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 & (country==2 | country==3 | country==6 | country==14 | country==15 | country==17 | country==21 | country==23 | country==24 | country==25) || country: , intmethod(mcaghermite)

log close


***
***Table 1, margins					
***
*para correr na feuc (com NOVEMBER_mm_dataset_2015.dta + new stats.do

*FIRST  Case A 
log using margins_A_MM, replace
set more off
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype woco_ext union_ext)  predict(mu fixedonly)
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  ud_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  centr_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.br woco_ext union_ext)  predict(mu fixedonly)

log close

log using margins_BCD_MM, replace
melogit KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype woco_ext union_ext)  predict(mu fixedonly)
melogit KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  ud_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  centr_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KOSICK i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.br woco_ext union_ext)  predict(mu fixedonly)


melogit KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype woco_ext union_ext)  predict(mu fixedonly)
melogit KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  ud_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  centr_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KORETEN i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.br woco_ext union_ext)  predict(mu fixedonly)


melogit KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype woco_ext union_ext)  predict(mu fixedonly)
melogit KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  ud_visser    woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  ud_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  centr_visser woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  centr_visser woco_ext union_ext)  predict(mu fixedonly)
melogit KOLOMOT i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.br                   woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.br woco_ext union_ext)  predict(mu fixedonly)

log close



***** ref #1/point 4/ suggestion
list woco_ext union_ext both1 union woco only_union only_woco if country==1

sum woco_ext union_ext both1 only_union only_woco if country==1
***
*** RE-RUN TABLE 1
***
* replication melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype               woco_ext union_ext q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
* for belgium toy
list  only_union only_woco both1 if country==1
* run with only_union only_woco both1 VERSUS NONE
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  only_union only_woco both1        q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country==1 || country: , intmethod(mcaghermite)


drop    both1_w
gen     both1_w=0
replace both1_w=1 if both1==1 & woco_ext==1

gen     both1_u=0
replace both1_u=1 if both1==1 & union_ext==1

list woco_ext union_ext both1 union woco only_union only_woco both1_w both1_u if country==1
sum woco_ext union_ext both1 only_union only_woco both1_w both1_u if country==1
by country: sum woco_ext union_ext both1 only_union only_woco both1_w both1_u
*fine
* run with nly_union only_woco both1 both1_w both1_u VERSUS NONE

*** FOR ALL COUNTRIES
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  only_union only_woco both1                 q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  only_union only_woco both1_w both1_u       q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)



***** PARA A FEUC

*** FOR ALL COUNTRIES BUT ONLY FOR KCLIMATE, COLUMN (1)
log using REFEREE_1_KCLIMATE, replace
set more off

melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  only_union only_woco both1                 q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  only_union only_woco both1)  predict(mu fixedonly)
melogit KCLIMATE_D i.NACE6_R1 APRIVATE i.est_size3 ASINGLE i.cbtype  only_union only_woco both1_w both1_u       q33perm_g q33wom_g q33univ_g q33pt_g E0A_A E0A_B E0A_C E0A_D E0A_E qh13 qh15 HVBPRES HVPINPER HVPGRPE HVPPRSH HVPSHOW               if country<29 || country: , intmethod(mcaghermite)
margins, dydx(i.cbtype  only_union only_woco both1_w both1_u)  predict(mu fixedonly)

log close