{smcl}
{com}{sf}{ul off}{txt}
{com}. 
. **** REPLICATION OF TABLE 3 ****
. 
. clear all
{res}{txt}
{com}. set mem 2g
{txt}{bf:set memory} ignored.
{p 4 4 2}
Memory no longer
needs to be set in modern Statas;
memory adjustments are performed on the fly
automatically.
{p_end}

{com}. set more off
{txt}
{com}. 
. /*      
>         OUTLINE
>         1. REGRESSIONS
>         2. TABLE FOR PAPER
> */
. 
. **************************
. ***** 1. REGRESSIONS *****
. **************************
. 
. // Import data
. 
. use $data\EEREV-D-21-00004_Data_Municipalities.dta, clear
{txt}
{com}. 
. keep if annee==1965
{txt}(11,053 observations deleted)

{com}. 
. // Macros
. 
. global geogr_endow              "water lelevation_mean lelevation_std ldist_de" 
{txt}
{com}. global pop                              "perc_b_1928_1945 perc_b_1946_1964 perc_b_1965_1980 perc_b_1981_1996 perc_male_counter"
{txt}
{com}. global demog_chara              "perc_fore $pop"
{txt}
{com}. global socec_aspec              "perc_bacplus perc_bluecollar perc_chom lrev_median lrev_sd"
{txt}
{com}. 
. // Control vectors
. 
. global sigma                    "pop_jew p_french_name lpop_1936 perc_fore_1936 perc_dial_1936 perc_disp_1946 $geogr_endow perc_droi_1936_ins" //
{txt}
{com}. global w                                "$demog_chara $socec_aspec" //
{txt}
{com}. global dialect                  "dialect_lux dialect_mos dialect_rhe dialect_pal dialect_als dialect_hau dialect_lor dialect_fra"
{txt}
{com}. 
. // Rdd
. 
. global rdd2                     "rdd_lat rdd_lon rdd_lat_2 rdd_lon_2 rdd_lat_lon"
{txt}
{com}. 
. // First column
. 
. ivreghdfe perc_i_full perc_e_1908_1927 if annee==1965 & sample==1, cluster(ct_1936_enc)
{res}
{txt}OLS estimation
{hline 14}

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity and clustering on ct_1936_enc

Number of clusters (ct_1936_enc) = {col 33}{res}    32{txt}{col 55}Number of obs = {res}     462
{txt}{col 55}F(  1,    31) = {res}  124.73
{txt}{col 55}Prob > F      = {res}  0.0000
{txt}Total (centered) SS     = {res} 9364.269099{txt}{col 55}Centered R2   = {res}  0.1903
{txt}Total (uncentered) SS   = {res} 37028.70761{txt}{col 55}Uncentered R2 = {res}  0.7952
{txt}Residual SS             = {res} 7582.262273{txt}{col 55}Root MSE      = {res}   4.051

{txt}{hline 17}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 18}{c |}{col 30}    Robust
{col 1}     perc_i_full{col 18}{c |} Coefficient{col 30}  std. err.{col 42}      z{col 50}   P>|z|{col 58}     [95% con{col 71}f. interval]
{hline 17}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
perc_e_1908_1927 {c |}{col 18}{res}{space 2} .3724989{col 30}{space 2} .0327919{col 41}{space 1}   11.36{col 50}{space 3}0.000{col 58}{space 4}  .308228{col 71}{space 3} .4367698
{txt}{space 11}_cons {c |}{col 18}{res}{space 2} .5854486{col 30}{space 2} .7816901{col 41}{space 1}    0.75{col 50}{space 3}0.454{col 58}{space 4}-.9466359{col 71}{space 3} 2.117533
{txt}{hline 17}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
Included instruments:{col 23}perc_e_1908_1927
{hline 78}

{com}. estimates store FS1
{txt}
{com}. test perc_e_1908_1927

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 = 0{p_end}

{txt}{col 12}chi2(  1) ={res}  129.04
{txt}{col 10}Prob > chi2 =  {res}  0.0000
{txt}
{com}. estadd scalar F_fs=r(chi2) //

{txt}added scalar:
               e(F_fs) =  {res}129.03781
{txt}
{com}. sum perc_i_full  if e(sample)

{txt}    Variable {c |}        Obs        Mean    Std. dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 1}perc_i_full {c |}{res}        462      7.7382    4.506989          0    39.3405
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}7.7381999
{txt}
{com}. 
. // Second column
. 
. ivreghdfe perc_i_full perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922 perc_b_1923_1932 if annee==1965 & sample==1, cluster(ct_1936_enc)
{res}
{txt}OLS estimation
{hline 14}

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity and clustering on ct_1936_enc

Number of clusters (ct_1936_enc) = {col 33}{res}    32{txt}{col 55}Number of obs = {res}     462
{txt}{col 55}F(  4,    31) = {res}   31.84
{txt}{col 55}Prob > F      = {res}  0.0000
{txt}Total (centered) SS     = {res} 9364.269099{txt}{col 55}Centered R2   = {res}  0.2545
{txt}Total (uncentered) SS   = {res} 37028.70761{txt}{col 55}Uncentered R2 = {res}  0.8115
{txt}Residual SS             = {res} 6980.779031{txt}{col 55}Root MSE      = {res}   3.887

{txt}{hline 17}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 18}{c |}{col 30}    Robust
{col 1}     perc_i_full{col 18}{c |} Coefficient{col 30}  std. err.{col 42}      z{col 50}   P>|z|{col 58}     [95% con{col 71}f. interval]
{hline 17}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
perc_e_1908_1927 {c |}{col 18}{res}{space 2} .4748617{col 30}{space 2} .0451449{col 41}{space 1}   10.52{col 50}{space 3}0.000{col 58}{space 4} .3863793{col 71}{space 3} .5633441
{txt}perc_b_1903_1912 {c |}{col 18}{res}{space 2}-.0185474{col 30}{space 2} .0638622{col 41}{space 1}   -0.29{col 50}{space 3}0.771{col 58}{space 4} -.143715{col 71}{space 3} .1066203
{txt}perc_b_1913_1922 {c |}{col 18}{res}{space 2}-.2988494{col 30}{space 2} .1459345{col 41}{space 1}   -2.05{col 50}{space 3}0.041{col 58}{space 4}-.5848757{col 71}{space 3}-.0128231
{txt}perc_b_1923_1932 {c |}{col 18}{res}{space 2}-.3100704{col 30}{space 2} .1134298{col 41}{space 1}   -2.73{col 50}{space 3}0.006{col 58}{space 4}-.5323886{col 71}{space 3}-.0877521
{txt}{space 11}_cons {c |}{col 18}{res}{space 2} 5.340521{col 30}{space 2} 1.334832{col 41}{space 1}    4.00{col 50}{space 3}0.000{col 58}{space 4} 2.724299{col 71}{space 3} 7.956743
{txt}{hline 17}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
Included instruments:{col 23}perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922
{col 23}perc_b_1923_1932
{hline 78}

{com}. lincom perc_e_1908_1927+perc_b_1903_1912 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1903_1912 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      z{col 46}   P>|z|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .4563143{col 26}{space 2}  .077455{col 37}{space 1}    5.89{col 46}{space 3}0.000{col 54}{space 4} .3045053{col 67}{space 3} .6081233
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1913_1922 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1913_1922 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      z{col 46}   P>|z|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1760123{col 26}{space 2} .1286442{col 37}{space 1}    1.37{col 46}{space 3}0.171{col 54}{space 4}-.0761256{col 67}{space 3} .4281502
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1923_1932 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1923_1932 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      z{col 46}   P>|z|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1647913{col 26}{space 2} .1097662{col 37}{space 1}    1.50{col 46}{space 3}0.133{col 54}{space 4}-.0503465{col 67}{space 3} .3799291
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estimates store FS2
{txt}
{com}. test perc_e_1908_1927

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 = 0{p_end}

{txt}{col 12}chi2(  1) ={res}  110.64
{txt}{col 10}Prob > chi2 =  {res}  0.0000
{txt}
{com}. estadd scalar F_fs=r(chi2) //

{txt}added scalar:
               e(F_fs) =  {res}110.64109
{txt}
{com}. sum perc_i_full  if e(sample)

{txt}    Variable {c |}        Obs        Mean    Std. dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 1}perc_i_full {c |}{res}        462      7.7382    4.506989          0    39.3405
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}7.7381999
{txt}
{com}. 
. // Third column
. 
. ivreghdfe perc_i_full perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922 perc_b_1923_1932 $rdd2 if annee==1965 & sample==1, cluster(ct_1936_enc) a(borderfe $dialect) 
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 6 iterations)
{res}
{txt}OLS estimation
{hline 14}

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity and clustering on ct_1936_enc

Number of clusters (ct_1936_enc) = {col 33}{res}    32{txt}{col 55}Number of obs = {res}     462
{txt}{col 55}F(  9,    31) = {res}   15.93
{txt}{col 55}Prob > F      = {res}  0.0000
{txt}Total (centered) SS     = {res} 7234.578995{txt}{col 55}Centered R2   = {res}  0.1534
{txt}Total (uncentered) SS   = {res} 7234.578995{txt}{col 55}Uncentered R2 = {res}  0.1534
{txt}Residual SS             = {res} 6124.789833{txt}{col 55}Root MSE      = {res}   3.727

{txt}{hline 17}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 18}{c |}{col 30}    Robust
{col 1}     perc_i_full{col 18}{c |} Coefficient{col 30}  std. err.{col 42}      t{col 50}   P>|t|{col 58}     [95% con{col 71}f. interval]
{hline 17}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
perc_e_1908_1927 {c |}{col 18}{res}{space 2} .3574214{col 30}{space 2} .0479645{col 41}{space 1}    7.45{col 50}{space 3}0.000{col 58}{space 4} .2595971{col 71}{space 3} .4552456
{txt}perc_b_1903_1912 {c |}{col 18}{res}{space 2}-.0462477{col 30}{space 2} .0575956{col 41}{space 1}   -0.80{col 50}{space 3}0.428{col 58}{space 4}-.1637148{col 71}{space 3} .0712193
{txt}perc_b_1913_1922 {c |}{col 18}{res}{space 2}-.1980167{col 30}{space 2} .1687533{col 41}{space 1}   -1.17{col 50}{space 3}0.250{col 58}{space 4}-.5421914{col 71}{space 3}  .146158
{txt}perc_b_1923_1932 {c |}{col 18}{res}{space 2}-.1882621{col 30}{space 2} .1277914{col 41}{space 1}   -1.47{col 50}{space 3}0.151{col 58}{space 4}-.4488944{col 71}{space 3} .0723702
{txt}{space 9}rdd_lat {c |}{col 18}{res}{space 2} 1699.332{col 30}{space 2} 630.9022{col 41}{space 1}    2.69{col 50}{space 3}0.011{col 58}{space 4} 412.5988{col 71}{space 3} 2986.066
{txt}{space 9}rdd_lon {c |}{col 18}{res}{space 2} 708.6403{col 30}{space 2} 434.7777{col 41}{space 1}    1.63{col 50}{space 3}0.113{col 58}{space 4}-178.0947{col 71}{space 3} 1595.375
{txt}{space 7}rdd_lat_2 {c |}{col 18}{res}{space 2} -16.5091{col 30}{space 2} 6.367524{col 41}{space 1}   -2.59{col 50}{space 3}0.014{col 58}{space 4}-29.49575{col 71}{space 3} -3.52245
{txt}{space 7}rdd_lon_2 {c |}{col 18}{res}{space 2}-7.062733{col 30}{space 2} 4.794588{col 41}{space 1}   -1.47{col 50}{space 3}0.151{col 58}{space 4}-16.84136{col 71}{space 3} 2.715894
{txt}{space 5}rdd_lat_lon {c |}{col 18}{res}{space 2}-12.36626{col 30}{space 2} 8.298235{col 41}{space 1}   -1.49{col 50}{space 3}0.146{col 58}{space 4}-29.29063{col 71}{space 3} 4.558099
{txt}{hline 17}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
Included instruments:{col 23}perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922
{col 23}perc_b_1923_1932 rdd_lat rdd_lon rdd_lat_2 rdd_lon_2
{col 23}rdd_lat_lon
{help ivreg2##partial:Partialled-out}:{col 23}_cons
{col 23}nb: total SS, model F and R2s are after partialling-out;
{col 23}    any {help ivreg2##s_small:small-sample adjustments} include partialled-out
{col 23}    variables in regressor count K
{hline 78}
{res}
{txt}Absorbed degrees of freedom:
{res}{col 1}{text}{hline 13}{c TT}{hline 12}{hline 12}{hline 14}{hline 1}{c TRC}
{col 1}{text} Absorbed FE{col 14}{c |} Categories{col 27} - Redundant{col 39}  = Num. Coefs{col 54}{c |}
{res}{col 1}{text}{hline 13}{c +}{hline 12}{hline 12}{hline 14}{hline 1}{c RT}
{col 1}{text}  borderfe25{col 14}{c |}{space 1}        8{col 27}{space 1}        0{col 39}{result}{space 1}        8{col 53}{text} {col 54}{c |}
{res}{col 1}{text} dialect_lux{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text} {col 54}{c |}
{res}{col 1}{text} dialect_mos{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_rhe{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_pal{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_als{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_hau{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_lor{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_fra{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text}{hline 13}{c BT}{hline 12}{hline 12}{hline 14}{hline 1}{c BRC}
? = number of redundant parameters may be higher
{res}{txt}
{com}. lincom perc_e_1908_1927+perc_b_1903_1912 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1903_1912 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .3111736{col 26}{space 2} .0722739{col 37}{space 1}    4.31{col 46}{space 3}0.000{col 54}{space 4} .1637699{col 67}{space 3} .4585773
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1913_1922 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1913_1922 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1594046{col 26}{space 2} .1419341{col 37}{space 1}    1.12{col 46}{space 3}0.270{col 54}{space 4}-.1300719{col 67}{space 3} .4488811
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1923_1932 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1923_1932 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1691593{col 26}{space 2} .1354902{col 37}{space 1}    1.25{col 46}{space 3}0.221{col 54}{space 4}-.1071747{col 67}{space 3} .4454933
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estimates store FS3
{txt}
{com}. test perc_e_1908_1927

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 = 0{p_end}

{txt}       F(  1,    31) ={res}   55.53
{txt}{col 13}Prob > F ={res}    0.0000
{txt}
{com}. estadd scalar F_fs=r(F)

{txt}added scalar:
               e(F_fs) =  {res}55.529147
{txt}
{com}. sum perc_i_full  if e(sample)

{txt}    Variable {c |}        Obs        Mean    Std. dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 1}perc_i_full {c |}{res}        462      7.7382    4.506989          0    39.3405
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}7.7381999
{txt}
{com}. 
. // Fourth column
. 
. ivreghdfe perc_i_full perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922 perc_b_1923_1932 $rdd2 $sigma perc_prot_church if annee==1965 & sample==1, cluster(ct_1936_enc) a(borderfe $dialect pop_jew water) partial($sigma)
{res}{txt}(warning: absorbing 11 dimensions of fixed effects; check that you really want that)
{res}{txt}note: {res}variable #12{txt} is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
{txt}note: {res}variable #18{txt} is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 8 iterations)
{res}
{txt}OLS estimation
{hline 14}

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity and clustering on ct_1936_enc

Number of clusters (ct_1936_enc) = {col 33}{res}    32{txt}{col 55}Number of obs = {res}     462
{txt}{col 55}F( 10,    31) = {res}   17.99
{txt}{col 55}Prob > F      = {res}  0.0000
{txt}Total (centered) SS     = {res} 6388.091868{txt}{col 55}Centered R2   = {res}  0.1043
{txt}Total (uncentered) SS   = {res} 6388.091868{txt}{col 55}Uncentered R2 = {res}  0.1043
{txt}Residual SS             = {res} 5721.714065{txt}{col 55}Root MSE      = {res}   3.661

{txt}{hline 17}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 18}{c |}{col 30}    Robust
{col 1}     perc_i_full{col 18}{c |} Coefficient{col 30}  std. err.{col 42}      t{col 50}   P>|t|{col 58}     [95% con{col 71}f. interval]
{hline 17}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
perc_e_1908_1927 {c |}{col 18}{res}{space 2} .2607967{col 30}{space 2} .0498459{col 41}{space 1}    5.23{col 50}{space 3}0.000{col 58}{space 4} .1591353{col 71}{space 3} .3624582
{txt}perc_b_1903_1912 {c |}{col 18}{res}{space 2}-.0113742{col 30}{space 2} .0559472{col 41}{space 1}   -0.20{col 50}{space 3}0.840{col 58}{space 4}-.1254793{col 71}{space 3} .1027309
{txt}perc_b_1913_1922 {c |}{col 18}{res}{space 2}-.1472673{col 30}{space 2} .1648041{col 41}{space 1}   -0.89{col 50}{space 3}0.378{col 58}{space 4}-.4833874{col 71}{space 3} .1888529
{txt}perc_b_1923_1932 {c |}{col 18}{res}{space 2}-.1142281{col 30}{space 2} .1258467{col 41}{space 1}   -0.91{col 50}{space 3}0.371{col 58}{space 4}-.3708942{col 71}{space 3} .1424379
{txt}{space 9}rdd_lat {c |}{col 18}{res}{space 2}-490.8896{col 30}{space 2} 958.8898{col 41}{space 1}   -0.51{col 50}{space 3}0.612{col 58}{space 4}-2446.558{col 71}{space 3} 1464.779
{txt}{space 9}rdd_lon {c |}{col 18}{res}{space 2} 531.5413{col 30}{space 2} 478.1289{col 41}{space 1}    1.11{col 50}{space 3}0.275{col 58}{space 4} -443.609{col 71}{space 3} 1506.692
{txt}{space 7}rdd_lat_2 {c |}{col 18}{res}{space 2} 5.808443{col 30}{space 2} 9.391831{col 41}{space 1}    0.62{col 50}{space 3}0.541{col 58}{space 4}-13.34632{col 71}{space 3} 24.96321
{txt}{space 7}rdd_lon_2 {c |}{col 18}{res}{space 2}-2.666936{col 30}{space 2} 4.511978{col 41}{space 1}   -0.59{col 50}{space 3}0.559{col 58}{space 4}-11.86918{col 71}{space 3} 6.535303
{txt}{space 5}rdd_lat_lon {c |}{col 18}{res}{space 2}-9.986794{col 30}{space 2} 9.684488{col 41}{space 1}   -1.03{col 50}{space 3}0.310{col 58}{space 4}-29.73844{col 71}{space 3}  9.76485
{txt}perc_prot_church {c |}{col 18}{res}{space 2} .0137781{col 30}{space 2} .0062947{col 41}{space 1}    2.19{col 50}{space 3}0.036{col 58}{space 4}   .00094{col 71}{space 3} .0266161
{txt}{hline 17}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
Included instruments:{col 23}perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922
{col 23}perc_b_1923_1932 rdd_lat rdd_lon rdd_lat_2 rdd_lon_2
{col 23}rdd_lat_lon perc_prot_church
{help ivreg2##partial:Partialled-out}:{col 23}pop_jew p_french_name lpop_1936 perc_fore_1936
{col 23}perc_dial_1936 perc_disp_1946 water lelevation_mean
{col 23}lelevation_std ldist_de perc_droi_1936_ins _cons
{col 23}nb: total SS, model F and R2s are after partialling-out;
{col 23}    any {help ivreg2##s_small:small-sample adjustments} include partialled-out
{col 23}    variables in regressor count K
{hline 78}
{res}
{txt}Absorbed degrees of freedom:
{res}{col 1}{text}{hline 13}{c TT}{hline 12}{hline 12}{hline 14}{hline 1}{c TRC}
{col 1}{text} Absorbed FE{col 14}{c |} Categories{col 27} - Redundant{col 39}  = Num. Coefs{col 54}{c |}
{res}{col 1}{text}{hline 13}{c +}{hline 12}{hline 12}{hline 14}{hline 1}{c RT}
{col 1}{text}  borderfe25{col 14}{c |}{space 1}        8{col 27}{space 1}        0{col 39}{result}{space 1}        8{col 53}{text} {col 54}{c |}
{res}{col 1}{text} dialect_lux{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text} {col 54}{c |}
{res}{col 1}{text} dialect_mos{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_rhe{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_pal{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_als{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_hau{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_lor{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_fra{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text}     pop_jew{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text}       water{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text}{hline 13}{c BT}{hline 12}{hline 12}{hline 14}{hline 1}{c BRC}
? = number of redundant parameters may be higher
{res}{txt}
{com}. lincom perc_e_1908_1927+perc_b_1903_1912 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1903_1912 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .2494226{col 26}{space 2} .0675266{col 37}{space 1}    3.69{col 46}{space 3}0.001{col 54}{space 4} .1117011{col 67}{space 3}  .387144
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1913_1922 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1913_1922 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1135295{col 26}{space 2}  .137268{col 37}{space 1}    0.83{col 46}{space 3}0.415{col 54}{space 4}-.1664304{col 67}{space 3} .3934893
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1923_1932 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1923_1932 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1465686{col 26}{space 2} .1229547{col 37}{space 1}    1.19{col 46}{space 3}0.242{col 54}{space 4}-.1041991{col 67}{space 3} .3973363
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estimates store FS4
{txt}
{com}. test perc_e_1908_1927

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 = 0{p_end}

{txt}       F(  1,    31) ={res}   27.37
{txt}{col 13}Prob > F ={res}    0.0000
{txt}
{com}. estadd scalar F_fs=r(F)

{txt}added scalar:
               e(F_fs) =  {res}27.374433
{txt}
{com}. sum perc_i_full  if e(sample)

{txt}    Variable {c |}        Obs        Mean    Std. dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 1}perc_i_full {c |}{res}        462      7.7382    4.506989          0    39.3405
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}7.7381999
{txt}
{com}. 
. // Fifth column
. 
. ivreghdfe perc_i_full perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922 perc_b_1923_1932 $rdd2 $sigma $w perc_prot_church lpop_counter if annee==1965 & sample==1, cluster(ct_1936_enc) a(borderfe $dialect pop_jew water) partial($sigma $w)
{res}{txt}(warning: absorbing 11 dimensions of fixed effects; check that you really want that)
{res}{txt}Warning - collinearities detected
Vars dropped:{col 21}perc_b_1981_1996
{res}{txt}note: {res}variable #13{txt} is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
{txt}note: {res}variable #19{txt} is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 8 iterations)
{res}
{txt}OLS estimation
{hline 14}

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity and clustering on ct_1936_enc

Number of clusters (ct_1936_enc) = {col 33}{res}    32{txt}{col 55}Number of obs = {res}     462
{txt}{col 55}F( 11,    31) = {res}   11.10
{txt}{col 55}Prob > F      = {res}  0.0000
{txt}Total (centered) SS     = {res} 6025.790907{txt}{col 55}Centered R2   = {res}  0.0757
{txt}Total (uncentered) SS   = {res} 6025.790907{txt}{col 55}Uncentered R2 = {res}  0.0757
{txt}Residual SS             = {res} 5569.532457{txt}{col 55}Root MSE      = {res}   3.659

{txt}{hline 17}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 18}{c |}{col 30}    Robust
{col 1}     perc_i_full{col 18}{c |} Coefficient{col 30}  std. err.{col 42}      t{col 50}   P>|t|{col 58}     [95% con{col 71}f. interval]
{hline 17}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
perc_e_1908_1927 {c |}{col 18}{res}{space 2} .2402204{col 30}{space 2} .0537875{col 41}{space 1}    4.47{col 50}{space 3}0.000{col 58}{space 4} .1305202{col 71}{space 3} .3499207
{txt}perc_b_1903_1912 {c |}{col 18}{res}{space 2}-.0067285{col 30}{space 2}  .055808{col 41}{space 1}   -0.12{col 50}{space 3}0.905{col 58}{space 4}-.1205497{col 71}{space 3} .1070927
{txt}perc_b_1913_1922 {c |}{col 18}{res}{space 2} -.125911{col 30}{space 2}  .159967{col 41}{space 1}   -0.79{col 50}{space 3}0.437{col 58}{space 4}-.4521658{col 71}{space 3} .2003439
{txt}perc_b_1923_1932 {c |}{col 18}{res}{space 2}-.0991447{col 30}{space 2} .1352625{col 41}{space 1}   -0.73{col 50}{space 3}0.469{col 58}{space 4}-.3750144{col 71}{space 3} .1767249
{txt}{space 9}rdd_lat {c |}{col 18}{res}{space 2}-856.8427{col 30}{space 2} 1040.226{col 41}{space 1}   -0.82{col 50}{space 3}0.416{col 58}{space 4}-2978.398{col 71}{space 3} 1264.712
{txt}{space 9}rdd_lon {c |}{col 18}{res}{space 2} 330.7684{col 30}{space 2} 536.4727{col 41}{space 1}    0.62{col 50}{space 3}0.542{col 58}{space 4} -763.375{col 71}{space 3} 1424.912
{txt}{space 7}rdd_lat_2 {c |}{col 18}{res}{space 2} 9.261741{col 30}{space 2} 10.26221{col 41}{space 1}    0.90{col 50}{space 3}0.374{col 58}{space 4}-11.66817{col 71}{space 3} 30.19165
{txt}{space 7}rdd_lon_2 {c |}{col 18}{res}{space 2} -2.13023{col 30}{space 2} 4.630607{col 41}{space 1}   -0.46{col 50}{space 3}0.649{col 58}{space 4}-11.57442{col 71}{space 3} 7.313955
{txt}{space 5}rdd_lat_lon {c |}{col 18}{res}{space 2}-6.032844{col 30}{space 2}  10.7801{col 41}{space 1}   -0.56{col 50}{space 3}0.580{col 58}{space 4}-28.01901{col 71}{space 3} 15.95332
{txt}perc_prot_church {c |}{col 18}{res}{space 2} .0102945{col 30}{space 2} .0067317{col 41}{space 1}    1.53{col 50}{space 3}0.136{col 58}{space 4} -.003435{col 71}{space 3} .0240239
{txt}{space 4}lpop_counter {c |}{col 18}{res}{space 2} .0409133{col 30}{space 2} 1.058093{col 41}{space 1}    0.04{col 50}{space 3}0.969{col 58}{space 4}-2.117081{col 71}{space 3} 2.198907
{txt}{hline 17}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
Included instruments:{col 23}perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922
{col 23}perc_b_1923_1932 rdd_lat rdd_lon rdd_lat_2 rdd_lon_2
{col 23}rdd_lat_lon perc_prot_church lpop_counter
{help ivreg2##partial:Partialled-out}:{col 23}pop_jew p_french_name lpop_1936 perc_fore_1936
{col 23}perc_dial_1936 perc_disp_1946 water lelevation_mean
{col 23}lelevation_std ldist_de perc_droi_1936_ins perc_fore
{col 23}perc_b_1928_1945 perc_b_1946_1964 perc_b_1965_1980
{col 23}perc_male_counter perc_bacplus perc_bluecollar perc_chom
{col 23}lrev_median lrev_sd _cons
{col 23}nb: total SS, model F and R2s are after partialling-out;
{col 23}    any {help ivreg2##s_small:small-sample adjustments} include partialled-out
{col 23}    variables in regressor count K
Dropped collinear:{col 23}perc_b_1981_1996
{hline 78}
{res}
{txt}Absorbed degrees of freedom:
{res}{col 1}{text}{hline 13}{c TT}{hline 12}{hline 12}{hline 14}{hline 1}{c TRC}
{col 1}{text} Absorbed FE{col 14}{c |} Categories{col 27} - Redundant{col 39}  = Num. Coefs{col 54}{c |}
{res}{col 1}{text}{hline 13}{c +}{hline 12}{hline 12}{hline 14}{hline 1}{c RT}
{col 1}{text}  borderfe25{col 14}{c |}{space 1}        8{col 27}{space 1}        0{col 39}{result}{space 1}        8{col 53}{text} {col 54}{c |}
{res}{col 1}{text} dialect_lux{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text} {col 54}{c |}
{res}{col 1}{text} dialect_mos{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_rhe{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_pal{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_als{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_hau{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_lor{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text} dialect_fra{col 14}{c |}{space 1}        1{col 27}{space 1}        1{col 39}{result}{space 1}        0{col 53}{text}?{col 54}{c |}
{res}{col 1}{text}     pop_jew{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text}       water{col 14}{c |}{space 1}        2{col 27}{space 1}        1{col 39}{result}{space 1}        1{col 53}{text}?{col 54}{c |}
{res}{col 1}{text}{hline 13}{c BT}{hline 12}{hline 12}{hline 14}{hline 1}{c BRC}
? = number of redundant parameters may be higher
{res}{txt}
{com}. lincom perc_e_1908_1927+perc_b_1903_1912 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1903_1912 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .2334919{col 26}{space 2} .0645933{col 37}{space 1}    3.61{col 46}{space 3}0.001{col 54}{space 4}  .101753{col 67}{space 3} .3652308
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1913_1922 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1913_1922 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1143095{col 26}{space 2} .1372648{col 37}{space 1}    0.83{col 46}{space 3}0.411{col 54}{space 4}-.1656439{col 67}{space 3} .3942628
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. lincom perc_e_1908_1927+perc_b_1923_1932 //

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 + perc_b_1923_1932 = 0{p_end}

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1} perc_i_full{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}(1) {c |}{col 14}{res}{space 2} .1410757{col 26}{space 2} .1278356{col 37}{space 1}    1.10{col 46}{space 3}0.278{col 54}{space 4}-.1196467{col 67}{space 3}  .401798
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estimates store FS5
{txt}
{com}. test perc_e_1908_1927

{p 0 7}{space 1}{text:( 1)}{space 1} {res}perc_e_1908_1927 = 0{p_end}

{txt}       F(  1,    31) ={res}   19.95
{txt}{col 13}Prob > F ={res}    0.0001
{txt}
{com}. estadd scalar F_fs=r(F)

{txt}added scalar:
               e(F_fs) =  {res}19.946086
{txt}
{com}. sum perc_i_full  if e(sample)

{txt}    Variable {c |}        Obs        Mean    Std. dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 1}perc_i_full {c |}{res}        462      7.7382    4.506989          0    39.3405
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}7.7381999
{txt}
{com}. 
. ******************************
. ***** 3. TABLE FOR PAPER *****
. ******************************
. 
. // Fixed effects
. 
. estfe . FS*, labels(dialect_lux "Historical dialects")
{res}{txt}
{com}. return list

{txt}macros:
        r(indicate_fe) : "{res} "Historical dialects=0.dialect_lux"  "borderfe25=0.borderfe25" "dialect_mos=0.dialect_mos" "dialect_rhe=0.dialect_rhe" "diale{txt}.."

{com}. 
. // Labels
. 
. la var perc_e_1908_1927         "\\ Eligible births (\%)"
{txt}
{com}. la var perc_b_1903_1912         "[0.3em] Born 1903-1912 (\%)"
{txt}
{com}. la var perc_b_1913_1922         "[0.3em] Born 1913-1922 (\%)"
{txt}
{com}. la var perc_b_1923_1932         "[0.3em] Born 1923-1932 (\%)"
{txt}
{com}. 
. // First-Stage
. 
. estout FS1 FS2 FS3 FS4 FS5 ///
>         using "$results\\$program.tex", replace style(tex) type label ///
>         cells(b(star fmt(%9.3f)) se(par fmt(%9.3f)))  ///
>         stats(F_fs, lab("\\ First-stage $ F $-statistic") fmt(2)) ///
>         collabels(none) ///
>         mlabels(, span none prefix(\multicolumn{c -(}@span{c )-}{c -(}c{c )-}{c -(}) suffix({c )-})) ///
>         msign(--) starlevels(* .1 ** .05 *** .01) ///
>         keep(perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922 perc_b_1923_1932) ///
>         order(perc_e_1908_1927 perc_b_1903_1912 perc_b_1913_1922 perc_b_1923_1932) ///
>         title(Estimation strategy: First-stage estimates) ///
>         prehead( ///
>         "\begin{c -(}table{c )-}[t!]" ///
>         "\begin{c -(}center{c )-}" ///
>         "\begin{c -(}threeparttable{c )-}" ///
>         "\begin{c -(}scriptsize{c )-}" ///
>         "\caption{c -(}@title\label{c -(}$program{c )-}{c )-}" ///
>         "\begin{c -(}tabular*{c )-}{c -(}\columnwidth{c )-}{c -(}@{c -(}\extracolsep{c -(}\stretch{c -(}1{c )-}{c )-}{c )-}l*{c -(}@M{c )-}{c -(}c{c )-}@{c -(}{c )-}{c )-}" ///
>         "\hline" "\hline" ///
>         "&\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(1){c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(2){c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(3){c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(4){c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(5){c )-}\\" ///
>         "\cmidrule(lr){c -(}2-6{c )-}" ///
>         "\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Specification{c )-}{c )-}&\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}No controls{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Fertility{c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Lat-lon poly{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Historical{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Historical and{c )-}\\" ///
>         " &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}controls{c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}and dialect{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}controls{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}contemporary{c )-}\\" ///
>         "\cmidrule(lr){c -(}2-6{c )-}" ///
>         "& \multicolumn{c -(}5{c )-}{c -(}c{c )-}{c -(}First-Stage Estimates. Dep. Variable: Conscription proxy (\%){c )-}\\") ///
>         posthead("\hline") //
{res}
\begin{c -(}table{c )-}[t!]
\begin{c -(}center{c )-}
\begin{c -(}threeparttable{c )-}
\begin{c -(}scriptsize{c )-}
\caption{c -(}Estimation strategy: First-stage estimates\label{c -(}EEREV-D-21-00004_Table3{c )-}{c )-}
\begin{c -(}tabular*{c )-}{c -(}\columnwidth{c )-}{c -(}@{c -(}\extracolsep{c -(}\stretch{c -(}1{c )-}{c )-}{c )-}l*{c -(}5{c )-}{c -(}c{c )-}@{c -(}{c )-}{c )-}
\hline
\hline
&\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(1){c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(2){c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(3){c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(4){c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(5){c )-}\\
\cmidrule(lr){c -(}2-6{c )-}
\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Specification{c )-}{c )-}&\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}No controls{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Fertility{c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Lat-lon poly{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Historical{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Historical and{c )-}\\
 &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}controls{c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}and dialect{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}controls{c )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}contemporary{c )-}\\
\cmidrule(lr){c -(}2-6{c )-}
& \multicolumn{c -(}5{c )-}{c -(}c{c )-}{c -(}First-Stage Estimates. Dep. Variable: Conscription proxy (\%){c )-}\\
\hline
\\ Eligible births (\%)&       0.372***&       0.475***&       0.357***&       0.261***&       0.240***\\
                    &     (0.033)   &     (0.045)   &     (0.048)   &     (0.050)   &     (0.054)   \\
[0.3em] Born 1903-1912 (\%)&               &     --0.019   &     --0.046   &     --0.011   &     --0.007   \\
                    &               &     (0.064)   &     (0.058)   &     (0.056)   &     (0.056)   \\
[0.3em] Born 1913-1922 (\%)&               &     --0.299** &     --0.198   &     --0.147   &     --0.126   \\
                    &               &     (0.146)   &     (0.169)   &     (0.165)   &     (0.160)   \\
[0.3em] Born 1923-1932 (\%)&               &     --0.310***&     --0.188   &     --0.114   &     --0.099   \\
                    &               &     (0.113)   &     (0.128)   &     (0.126)   &     (0.135)   \\
\\ First-stage $ F $-statistic&      129.04   &      110.64   &       55.53   &       27.37   &       19.95   \\
{txt}(output written to {browse  `"D:\Dropbox\SVLACHOS\Research\Malgre_nous\Replication_files\Results\EEREV-D-21-00004_Table3.tex"'})

{com}. 
. estout FS1 FS2 FS3 FS4 FS5 ///
>         using "$results\\$program.tex", append style(tex) type label ///
>         cells(b(star fmt(%9.3f)) se(par fmt(%9.3f)))  ///
>         stats(mean_y N N_clust, lab("Mean dep. variable" "Observations" "Clusters") fmt(2 0 0)) ///
>         collabels(none) ///
>         mlabels(, span none prefix(\multicolumn{c -(}@span{c )-}{c -(}c{c )-}{c -(}) suffix({c )-})) ///
>         msign(--) starlevels(* .1 ** .05 *** .01) ///
>         drop(*) ///
>         indicate( ///   
>         "Border segment FE = 0.borderfe25" ///
>         "Historical dialect FE = 0.dialect_lux" ///
>         "Historical controls = perc_prot_church" ///
>         "Contemporary controls (1965) = lpop_counter" ///
>         , labels(Yes )) ///
>         title(First-stage estimates) ///
>         posthead("\\ Lat-lon polynomial &  &  & 2nd & 2nd & 2nd \\ ") ///
>         prefoot("\hline") ///
>         postfoot("\hline" "\hline" ///
>         "\end{c -(}tabular*{c )-}" /// 
>         "\end{c -(}scriptsize{c )-}" ///
>         "\begin{c -(}tablenotes{c )-}[normal,flushleft]" ///
>         "\begin{c -(}scriptsize{c )-}" ///
>         "\item \emph{c -(}Notes{c )-}:" /// 
>         "First-stage estimates of the effect of the draft rule on Wehrmacht conscription. " ///
>         "The unit of observation is a municipality. " ///
>         "Standard errors clustered at the canton level in parentheses. " ///
>         "Columns (3) to (5) include a quadratic polynomial in latitude and longitude ($ x+y+x^2+y^2+xy $). " ///
>         "First-stage F-statistic is the Kleibergen-Paap rk Wald F-statistic. " ///
>         "* significant at 10\%; ** at 5\%; *** at 1\%. " ///    
>         "\end{c -(}scriptsize{c )-}" ///
>         "\end{c -(}tablenotes{c )-}" ///
>         "\end{c -(}threeparttable{c )-}" ///
>         "\end{c -(}center{c )-}" ///
>         "\end{c -(}table{c )-}")
{res}
\\ Lat-lon polynomial &  &  & 2nd & 2nd & 2nd \\ 
Border segment FE   &               &               &         Yes   &         Yes   &         Yes   \\
Historical dialect FE &               &               &         Yes   &         Yes   &         Yes   \\
Historical controls &               &               &               &         Yes   &         Yes   \\
Contemporary controls (1965) &               &               &               &               &         Yes   \\
\hline
Mean dep. variable  &        7.74   &        7.74   &        7.74   &        7.74   &        7.74   \\
Observations        &         462   &         462   &         462   &         462   &         462   \\
Clusters            &          32   &          32   &          32   &          32   &          32   \\
\hline
\hline
\end{c -(}tabular*{c )-}
\end{c -(}scriptsize{c )-}
\begin{c -(}tablenotes{c )-}[normal,flushleft]
\begin{c -(}scriptsize{c )-}
\item \emph{c -(}Notes{c )-}:
First-stage estimates of the effect of the draft rule on Wehrmacht conscription. 
The unit of observation is a municipality. 
Standard errors clustered at the canton level in parentheses. 
Columns (3) to (5) include a quadratic polynomial in latitude and longitude ($ x+y+x^2+y^2+xy $). 
First-stage F-statistic is the Kleibergen-Paap rk Wald F-statistic. 
* significant at 10\%; ** at 5\%; *** at 1\%. 
\end{c -(}scriptsize{c )-}
\end{c -(}tablenotes{c )-}
\end{c -(}threeparttable{c )-}
\end{c -(}center{c )-}
\end{c -(}table{c )-}
{txt}(output written to {browse  `"D:\Dropbox\SVLACHOS\Research\Malgre_nous\Replication_files\Results\EEREV-D-21-00004_Table3.tex"'})

{com}. 
. cap log close
{smcl}
{com}{sf}{ul off}