{smcl}
{com}{sf}{ul off}{txt}
{com}. 
. **** REPLICATION OF TABLE 1 ****
. 
. 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. RUN REGRESSIONS
>         2. TABLE FOR PAPER
> */
. 
. ******************************
. ***** 1. RUN REGRESSIONS *****
. ******************************
. 
. // Improt Microcensus data
. 
. use $data\EEREV-D-21-00004_Data_Micro1936.dta, clear
{txt}
{com}. 
. keep if male==1
{txt}(7,640 observations deleted)

{com}. 
. // Controls vector
. 
. global controls "frenchnat germannat protestant other speakfrench hhold_head hhold_size"
{txt}
{com}. 
. // Full sample (eligibility)
. 
. reghdfe index eligible, cluster(birthyr) noabsorb
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 1 iterations)
{res}
{txt}HDFE Linear regression{col 51}Number of obs{col 67}= {res}     7,590
{txt}Absorbing 1 HDFE group{col 51}F({res}   1{txt},{res}     92{txt}){col 67}= {res}     55.24
{txt}Statistics robust to heteroskedasticity{col 51}Prob > F{col 67}= {res}    0.0000
{txt}{col 51}R-squared{col 67}= {res}    0.1048
{txt}{col 51}Adj R-squared{col 67}= {res}    0.1047
{txt}{col 51}Within R-sq.{col 67}= {res}    0.1048
{txt}{col 1}Number of clusters ({res}birthyr{txt}) {col 30}= {res}        93{txt}{col 51}Root MSE{col 67}= {res}    0.1993

{txt}{ralign 78:(Std. err. adjusted for {res:93} clusters in birthyr)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}       index{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 4}eligible {c |}{col 14}{res}{space 2} .1482915{col 26}{space 2} .0199525{col 37}{space 1}    7.43{col 46}{space 3}0.000{col 54}{space 4} .1086641{col 67}{space 3} .1879188
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .0015131{col 26}{space 2} .0006074{col 37}{space 1}    2.49{col 46}{space 3}0.015{col 54}{space 4} .0003068{col 67}{space 3} .0027195
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estimates store OLS
{txt}
{com}. qui sum index if e(sample)
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}.04650856
{txt}
{com}. 
. // Within municipality
. 
. reghdfe index eligible age age2 $controls, cluster(birthyr) a(com_hist)
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 1 iterations)
{res}
{txt}HDFE Linear regression{col 51}Number of obs{col 67}= {res}     7,590
{txt}Absorbing 1 HDFE group{col 51}F({res}  10{txt},{res}     92{txt}){col 67}= {res}     11.08
{txt}Statistics robust to heteroskedasticity{col 51}Prob > F{col 67}= {res}    0.0000
{txt}{col 51}R-squared{col 67}= {res}    0.1315
{txt}{col 51}Adj R-squared{col 67}= {res}    0.1148
{txt}{col 51}Within R-sq.{col 67}= {res}    0.1103
{txt}{col 1}Number of clusters ({res}birthyr{txt}) {col 30}= {res}        93{txt}{col 51}Root MSE{col 67}= {res}    0.1981

{txt}{ralign 78:(Std. err. adjusted for {res:93} clusters in birthyr)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}       index{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 4}eligible {c |}{col 14}{res}{space 2} .1292856{col 26}{space 2} .0208803{col 37}{space 1}    6.19{col 46}{space 3}0.000{col 54}{space 4} .0878155{col 67}{space 3} .1707557
{txt}{space 9}age {c |}{col 14}{res}{space 2}-.0015417{col 26}{space 2} .0010687{col 37}{space 1}   -1.44{col 46}{space 3}0.153{col 54}{space 4}-.0036642{col 67}{space 3} .0005808
{txt}{space 8}age2 {c |}{col 14}{res}{space 2}  .000014{col 26}{space 2} .0000101{col 37}{space 1}    1.38{col 46}{space 3}0.171{col 54}{space 4}-6.16e-06{col 67}{space 3} .0000341
{txt}{space 3}frenchnat {c |}{col 14}{res}{space 2} .0199409{col 26}{space 2}  .015251{col 37}{space 1}    1.31{col 46}{space 3}0.194{col 54}{space 4}-.0103489{col 67}{space 3} .0502307
{txt}{space 3}germannat {c |}{col 14}{res}{space 2} .1062829{col 26}{space 2} .0824743{col 37}{space 1}    1.29{col 46}{space 3}0.201{col 54}{space 4}-.0575183{col 67}{space 3}  .270084
{txt}{space 2}protestant {c |}{col 14}{res}{space 2}  .004189{col 26}{space 2} .0090286{col 37}{space 1}    0.46{col 46}{space 3}0.644{col 54}{space 4}-.0137426{col 67}{space 3} .0221206
{txt}{space 7}other {c |}{col 14}{res}{space 2}-.0341502{col 26}{space 2} .0144782{col 37}{space 1}   -2.36{col 46}{space 3}0.020{col 54}{space 4}-.0629052{col 67}{space 3}-.0053951
{txt}{space 1}speakfrench {c |}{col 14}{res}{space 2} .0240412{col 26}{space 2} .0064918{col 37}{space 1}    3.70{col 46}{space 3}0.000{col 54}{space 4}  .011148{col 67}{space 3} .0369343
{txt}{space 2}hhold_head {c |}{col 14}{res}{space 2} .0073702{col 26}{space 2}  .010834{col 37}{space 1}    0.68{col 46}{space 3}0.498{col 54}{space 4} -.014147{col 67}{space 3} .0288875
{txt}{space 2}hhold_size {c |}{col 14}{res}{space 2} .0017035{col 26}{space 2} .0012546{col 37}{space 1}    1.36{col 46}{space 3}0.178{col 54}{space 4}-.0007883{col 67}{space 3} .0041953
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-.0032605{col 26}{space 2} .0219381{col 37}{space 1}   -0.15{col 46}{space 3}0.882{col 54}{space 4}-.0468315{col 67}{space 3} .0403105
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{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}    com_hist{col 14}{c |}{space 1}      134{col 27}{space 1}        0{col 39}{result}{space 1}      134{col 53}{text} {col 54}{c |}
{res}{col 1}{text}{hline 13}{c BT}{hline 12}{hline 12}{hline 14}{hline 1}{c BRC}
{res}{txt}
{com}. estimates store OLS_FE
{txt}
{com}. qui sum index if e(sample)
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}.04650856
{txt}
{com}. test age=age2=frenchnat=germannat=other=protestant=speakfrench=hhold_head=hhold_size=0

{p 0 7}{space 1}{text:( 1)}{space 1} {res}age - age2 = 0{p_end}
{p 0 7}{space 1}{text:( 2)}{space 1} age - frenchnat = 0{p_end}
{p 0 7}{space 1}{text:( 3)}{space 1} age - germannat = 0{p_end}
{p 0 7}{space 1}{text:( 4)}{space 1} age - other = 0{p_end}
{p 0 7}{space 1}{text:( 5)}{space 1} age - protestant = 0{p_end}
{p 0 7}{space 1}{text:( 6)}{space 1} age - speakfrench = 0{p_end}
{p 0 7}{space 1}{text:( 7)}{space 1} age - hhold_head = 0{p_end}
{p 0 7}{space 1}{text:( 8)}{space 1} age - hhold_size = 0{p_end}
{p 0 7}{space 1}{text:( 9)}{space 1} age = 0{p_end}

{txt}       F(  9,    92) ={res}    2.25
{txt}{col 13}Prob > F ={res}    0.0255
{txt}
{com}. scalar F_stat=r(p)
{txt}
{com}. estadd scalar F_stat

{txt}added scalar:
             e(F_stat) =  {res}.02552621
{txt}
{com}.                 
. // Within yob
. 
. reghdfe index $controls, cluster(birthyr) a(com_hist birthyr) // keepsing // male eligible age age2 
{res}{txt}(dropped 2 {browse "http://scorreia.com/research/singletons.pdf":singleton observations})
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 6 iterations)
{res}
{txt}HDFE Linear regression{col 51}Number of obs{col 67}= {res}     7,588
{txt}Absorbing 2 HDFE groups{col 51}F({res}   7{txt},{res}     90{txt}){col 67}= {res}      2.22
{txt}Statistics robust to heteroskedasticity{col 51}Prob > F{col 67}= {res}    0.0394
{txt}{col 51}R-squared{col 67}= {res}    0.1743
{txt}{col 51}Adj R-squared{col 67}= {res}    0.1483
{txt}{col 51}Within R-sq.{col 67}= {res}    0.0013
{txt}{col 1}Number of clusters ({res}birthyr{txt}) {col 30}= {res}        91{txt}{col 51}Root MSE{col 67}= {res}    0.1944

{txt}{ralign 78:(Std. err. adjusted for {res:91} clusters in birthyr)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}       index{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 3}frenchnat {c |}{col 14}{res}{space 2} .0180231{col 26}{space 2} .0148684{col 37}{space 1}    1.21{col 46}{space 3}0.229{col 54}{space 4}-.0115156{col 67}{space 3} .0475619
{txt}{space 3}germannat {c |}{col 14}{res}{space 2} .1090726{col 26}{space 2} .0842755{col 37}{space 1}    1.29{col 46}{space 3}0.199{col 54}{space 4}-.0583554{col 67}{space 3} .2765006
{txt}{space 2}protestant {c |}{col 14}{res}{space 2} .0043153{col 26}{space 2} .0089904{col 37}{space 1}    0.48{col 46}{space 3}0.632{col 54}{space 4}-.0135456{col 67}{space 3} .0221763
{txt}{space 7}other {c |}{col 14}{res}{space 2}-.0268937{col 26}{space 2} .0148876{col 37}{space 1}   -1.81{col 46}{space 3}0.074{col 54}{space 4}-.0564706{col 67}{space 3} .0026831
{txt}{space 1}speakfrench {c |}{col 14}{res}{space 2} .0106739{col 26}{space 2} .0045627{col 37}{space 1}    2.34{col 46}{space 3}0.022{col 54}{space 4} .0016093{col 67}{space 3} .0197386
{txt}{space 2}hhold_head {c |}{col 14}{res}{space 2} .0023552{col 26}{space 2} .0063799{col 37}{space 1}    0.37{col 46}{space 3}0.713{col 54}{space 4}-.0103195{col 67}{space 3} .0150299
{txt}{space 2}hhold_size {c |}{col 14}{res}{space 2}-.0001624{col 26}{space 2} .0010197{col 37}{space 1}   -0.16{col 46}{space 3}0.874{col 54}{space 4}-.0021883{col 67}{space 3} .0018634
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .0215201{col 26}{space 2} .0161225{col 37}{space 1}    1.33{col 46}{space 3}0.185{col 54}{space 4}-.0105101{col 67}{space 3} .0535503
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{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}    com_hist{col 14}{c |}{space 1}      134{col 27}{space 1}        0{col 39}{result}{space 1}      134{col 53}{text} {col 54}{c |}
{res}{col 1}{text}     birthyr{col 14}{c |}{space 1}       91{col 27}{space 1}       91{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}
* = FE nested within cluster; treated as redundant for DoF computation
{res}{txt}
{com}. estimates store com_ctr
{txt}
{com}. qui sum index if e(sample)
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}.04652082
{txt}
{com}. test frenchnat=germannat=other=protestant=speakfrench=hhold_head=hhold_size=0 // age=age2=

{p 0 7}{space 1}{text:( 1)}{space 1} {res}frenchnat - germannat = 0{p_end}
{p 0 7}{space 1}{text:( 2)}{space 1} frenchnat - other = 0{p_end}
{p 0 7}{space 1}{text:( 3)}{space 1} frenchnat - protestant = 0{p_end}
{p 0 7}{space 1}{text:( 4)}{space 1} frenchnat - speakfrench = 0{p_end}
{p 0 7}{space 1}{text:( 5)}{space 1} frenchnat - hhold_head = 0{p_end}
{p 0 7}{space 1}{text:( 6)}{space 1} frenchnat - hhold_size = 0{p_end}
{p 0 7}{space 1}{text:( 7)}{space 1} frenchnat = 0{p_end}

{txt}       F(  7,    90) ={res}    2.22
{txt}{col 13}Prob > F ={res}    0.0394
{txt}
{com}. scalar F_stat=r(p)
{txt}
{com}. estadd scalar F_stat

{txt}added scalar:
             e(F_stat) =  {res}.03938286
{txt}
{com}. 
. // Only eligible
. 
. reghdfe index age age2 $controls if eligible==1, cluster(birthyr) a(com_hist) // keepsing
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 1 iterations)
{res}
{txt}HDFE Linear regression{col 51}Number of obs{col 67}= {res}     2,303
{txt}Absorbing 1 HDFE group{col 51}F({res}   9{txt},{res}     19{txt}){col 67}= {res}      6.88
{txt}Statistics robust to heteroskedasticity{col 51}Prob > F{col 67}= {res}    0.0002
{txt}{col 51}R-squared{col 67}= {res}    0.1021
{txt}{col 51}Adj R-squared{col 67}= {res}    0.0431
{txt}{col 51}Within R-sq.{col 67}= {res}    0.0247
{txt}{col 1}Number of clusters ({res}birthyr{txt}) {col 30}= {res}        20{txt}{col 51}Root MSE{col 67}= {res}    0.3492

{txt}{ralign 78:(Std. err. adjusted for {res:20} clusters in birthyr)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}       index{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}age {c |}{col 14}{res}{space 2} .0346422{col 26}{space 2} .0324638{col 37}{space 1}    1.07{col 46}{space 3}0.299{col 54}{space 4}-.0333053{col 67}{space 3} .1025897
{txt}{space 8}age2 {c |}{col 14}{res}{space 2}-.0008673{col 26}{space 2} .0006191{col 37}{space 1}   -1.40{col 46}{space 3}0.177{col 54}{space 4}-.0021632{col 67}{space 3} .0004286
{txt}{space 3}frenchnat {c |}{col 14}{res}{space 2} .0478196{col 26}{space 2} .0524002{col 37}{space 1}    0.91{col 46}{space 3}0.373{col 54}{space 4}-.0618551{col 67}{space 3} .1574944
{txt}{space 3}germannat {c |}{col 14}{res}{space 2} .1385904{col 26}{space 2}  .138399{col 37}{space 1}    1.00{col 46}{space 3}0.329{col 54}{space 4}-.1510821{col 67}{space 3} .4282629
{txt}{space 2}protestant {c |}{col 14}{res}{space 2} .0018264{col 26}{space 2} .0287847{col 37}{space 1}    0.06{col 46}{space 3}0.950{col 54}{space 4}-.0584208{col 67}{space 3} .0620735
{txt}{space 7}other {c |}{col 14}{res}{space 2}-.0792935{col 26}{space 2} .0303048{col 37}{space 1}   -2.62{col 46}{space 3}0.017{col 54}{space 4}-.1427222{col 67}{space 3}-.0158649
{txt}{space 1}speakfrench {c |}{col 14}{res}{space 2} .0390032{col 26}{space 2} .0184571{col 37}{space 1}    2.11{col 46}{space 3}0.048{col 54}{space 4}  .000372{col 67}{space 3} .0776344
{txt}{space 2}hhold_head {c |}{col 14}{res}{space 2} .0407289{col 26}{space 2} .0299807{col 37}{space 1}    1.36{col 46}{space 3}0.190{col 54}{space 4}-.0220214{col 67}{space 3} .1034791
{txt}{space 2}hhold_size {c |}{col 14}{res}{space 2} .0004328{col 26}{space 2} .0038473{col 37}{space 1}    0.11{col 46}{space 3}0.912{col 54}{space 4}-.0076197{col 67}{space 3} .0084852
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-.2371316{col 26}{space 2} .4336956{col 37}{space 1}   -0.55{col 46}{space 3}0.591{col 54}{space 4}-1.144867{col 67}{space 3} .6706037
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{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}    com_hist{col 14}{c |}{space 1}      134{col 27}{space 1}        0{col 39}{result}{space 1}      134{col 53}{text} {col 54}{c |}
{res}{col 1}{text}{hline 13}{c BT}{hline 12}{hline 12}{hline 14}{hline 1}{c BRC}
{res}{txt}
{com}. estimates store eligible_ctr
{txt}
{com}. sum index if e(sample)

{txt}    Variable {c |}        Obs        Mean    Std. dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}index {c |}{res}      2,303    .1498046    .3569573          0          1
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}.1498046
{txt}
{com}. margins, dydx(age) atmeans //  
{res}
{txt}Conditional marginal effects{col 58}{lalign 13:Number of obs}{col 71} = {res}{ralign 5:2,303}
{txt}{col 1}Model VCE: {res:Robust}

{txt}{p2colset 1 13 13 2}{...}
{p2col:Expression:}{res:Linear prediction, predict()}{p_end}
{p2col:dy/dx wrt:}{res:age}{p_end}
{p2colreset}{...}
{lalign 4:At: }{space 0}{lalign 11:age} = {res:{ralign 8:23.96049}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:age2} = {res:{ralign 8:612.2584}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:frenchnat} = {res:{ralign 8:.9913157}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:germannat} = {res:{ralign 8:.0034737}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:protestant} = {res:{ralign 8:.578376}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:other} = {res:{ralign 8:.0069475}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:speakfrench} = {res:{ralign 8:.8749457}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:hhold_head} = {res:{ralign 8:.0785931}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:hhold_size} = {res:{ralign 8:5.070343}} {txt:(mean)}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26} Delta-method
{col 14}{c |}      dy/dx{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}age {c |}{col 14}{res}{space 2} .0346422{col 26}{space 2} .0324638{col 37}{space 1}    1.07{col 46}{space 3}0.286{col 54}{space 4}-.0289856{col 67}{space 3} .0982701
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{txt}
{com}. test age=age2=frenchnat=germannat=other=protestant=speakfrench=hhold_head=hhold_size=0

{p 0 7}{space 1}{text:( 1)}{space 1} {res}age - age2 = 0{p_end}
{p 0 7}{space 1}{text:( 2)}{space 1} age - frenchnat = 0{p_end}
{p 0 7}{space 1}{text:( 3)}{space 1} age - germannat = 0{p_end}
{p 0 7}{space 1}{text:( 4)}{space 1} age - other = 0{p_end}
{p 0 7}{space 1}{text:( 5)}{space 1} age - protestant = 0{p_end}
{p 0 7}{space 1}{text:( 6)}{space 1} age - speakfrench = 0{p_end}
{p 0 7}{space 1}{text:( 7)}{space 1} age - hhold_head = 0{p_end}
{p 0 7}{space 1}{text:( 8)}{space 1} age - hhold_size = 0{p_end}
{p 0 7}{space 1}{text:( 9)}{space 1} age = 0{p_end}

{txt}       F(  9,    19) ={res}    6.88
{txt}{col 13}Prob > F ={res}    0.0002
{txt}
{com}. scalar F_stat=r(p)
{txt}
{com}. estadd scalar F_stat

{txt}added scalar:
             e(F_stat) =  {res}.00021559
{txt}
{com}. 
. // Within birthyr eligible
. 
. reghdfe index $controls if eligible==1, cluster(birthyr) a(com_hist birthyr) // keepsing
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 6 iterations)
{res}
{txt}HDFE Linear regression{col 51}Number of obs{col 67}= {res}     2,303
{txt}Absorbing 2 HDFE groups{col 51}F({res}   7{txt},{res}     19{txt}){col 67}= {res}      1.98
{txt}Statistics robust to heteroskedasticity{col 51}Prob > F{col 67}= {res}    0.1116
{txt}{col 51}R-squared{col 67}= {res}    0.1306
{txt}{col 51}Adj R-squared{col 67}= {res}    0.0656
{txt}{col 51}Within R-sq.{col 67}= {res}    0.0014
{txt}{col 1}Number of clusters ({res}birthyr{txt}) {col 30}= {res}        20{txt}{col 51}Root MSE{col 67}= {res}    0.3450

{txt}{ralign 78:(Std. err. adjusted for {res:20} clusters in birthyr)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}       index{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 3}frenchnat {c |}{col 14}{res}{space 2} .0575122{col 26}{space 2} .0513894{col 37}{space 1}    1.12{col 46}{space 3}0.277{col 54}{space 4}-.0500471{col 67}{space 3} .1650716
{txt}{space 3}germannat {c |}{col 14}{res}{space 2} .1314519{col 26}{space 2}  .148152{col 37}{space 1}    0.89{col 46}{space 3}0.386{col 54}{space 4}-.1786338{col 67}{space 3} .4415375
{txt}{space 2}protestant {c |}{col 14}{res}{space 2} .0061547{col 26}{space 2} .0275548{col 37}{space 1}    0.22{col 46}{space 3}0.826{col 54}{space 4}-.0515181{col 67}{space 3} .0638274
{txt}{space 7}other {c |}{col 14}{res}{space 2}-.0740841{col 26}{space 2} .0307371{col 37}{space 1}   -2.41{col 46}{space 3}0.026{col 54}{space 4}-.1384175{col 67}{space 3}-.0097506
{txt}{space 1}speakfrench {c |}{col 14}{res}{space 2} .0257872{col 26}{space 2} .0176661{col 37}{space 1}    1.46{col 46}{space 3}0.161{col 54}{space 4}-.0111884{col 67}{space 3} .0627628
{txt}{space 2}hhold_head {c |}{col 14}{res}{space 2} .0253033{col 26}{space 2} .0310484{col 37}{space 1}    0.81{col 46}{space 3}0.425{col 54}{space 4}-.0396818{col 67}{space 3} .0902885
{txt}{space 2}hhold_size {c |}{col 14}{res}{space 2} .0002348{col 26}{space 2} .0037857{col 37}{space 1}    0.06{col 46}{space 3}0.951{col 54}{space 4}-.0076887{col 67}{space 3} .0081583
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .0635485{col 26}{space 2} .0629781{col 37}{space 1}    1.01{col 46}{space 3}0.326{col 54}{space 4}-.0682662{col 67}{space 3} .1953632
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{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}    com_hist{col 14}{c |}{space 1}      134{col 27}{space 1}        0{col 39}{result}{space 1}      134{col 53}{text} {col 54}{c |}
{res}{col 1}{text}     birthyr{col 14}{c |}{space 1}       20{col 27}{space 1}       20{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}
* = FE nested within cluster; treated as redundant for DoF computation
{res}{txt}
{com}. estimates store yob_ctr
{txt}
{com}. qui sum index if e(sample)
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}.1498046
{txt}
{com}. test frenchnat=germannat=other=protestant=speakfrench=hhold_head=hhold_size=0

{p 0 7}{space 1}{text:( 1)}{space 1} {res}frenchnat - germannat = 0{p_end}
{p 0 7}{space 1}{text:( 2)}{space 1} frenchnat - other = 0{p_end}
{p 0 7}{space 1}{text:( 3)}{space 1} frenchnat - protestant = 0{p_end}
{p 0 7}{space 1}{text:( 4)}{space 1} frenchnat - speakfrench = 0{p_end}
{p 0 7}{space 1}{text:( 5)}{space 1} frenchnat - hhold_head = 0{p_end}
{p 0 7}{space 1}{text:( 6)}{space 1} frenchnat - hhold_size = 0{p_end}
{p 0 7}{space 1}{text:( 7)}{space 1} frenchnat = 0{p_end}

{txt}       F(  7,    19) ={res}    1.98
{txt}{col 13}Prob > F ={res}    0.1116
{txt}
{com}. scalar F_stat=r(p)
{txt}
{com}. estadd scalar F_stat

{txt}added scalar:
             e(F_stat) =  {res}.11157123
{txt}
{com}. 
. // Only always-takers
. 
. reghdfe index age age2 $controls if eligible==0 & male==1, cluster(birthyr) a(com_hist) // keepsing
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 1 iterations)
{res}
{txt}HDFE Linear regression{col 51}Number of obs{col 67}= {res}     5,287
{txt}Absorbing 1 HDFE group{col 51}F({res}   9{txt},{res}     72{txt}){col 67}= {res}      1.00
{txt}Statistics robust to heteroskedasticity{col 51}Prob > F{col 67}= {res}    0.4452
{txt}{col 51}R-squared{col 67}= {res}    0.0527
{txt}{col 51}Adj R-squared{col 67}= {res}    0.0265
{txt}{col 51}Within R-sq.{col 67}= {res}    0.0284
{txt}{col 1}Number of clusters ({res}birthyr{txt}) {col 30}= {res}        73{txt}{col 51}Root MSE{col 67}= {res}    0.0384

{txt}{ralign 78:(Std. err. adjusted for {res:73} clusters in birthyr)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}       index{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}age {c |}{col 14}{res}{space 2} .0001121{col 26}{space 2} .0001717{col 37}{space 1}    0.65{col 46}{space 3}0.516{col 54}{space 4}-.0002302{col 67}{space 3} .0004544
{txt}{space 8}age2 {c |}{col 14}{res}{space 2}-1.54e-06{col 26}{space 2} 1.64e-06{col 37}{space 1}   -0.94{col 46}{space 3}0.351{col 54}{space 4}-4.82e-06{col 67}{space 3} 1.73e-06
{txt}{space 3}frenchnat {c |}{col 14}{res}{space 2} .0022823{col 26}{space 2} .0021904{col 37}{space 1}    1.04{col 46}{space 3}0.301{col 54}{space 4}-.0020841{col 67}{space 3} .0066486
{txt}{space 3}germannat {c |}{col 14}{res}{space 2}  .108076{col 26}{space 2} .0997515{col 37}{space 1}    1.08{col 46}{space 3}0.282{col 54}{space 4} -.090775{col 67}{space 3}  .306927
{txt}{space 2}protestant {c |}{col 14}{res}{space 2} .0009974{col 26}{space 2} .0016199{col 37}{space 1}    0.62{col 46}{space 3}0.540{col 54}{space 4}-.0022318{col 67}{space 3} .0042266
{txt}{space 7}other {c |}{col 14}{res}{space 2}-.0003009{col 26}{space 2}  .001898{col 37}{space 1}   -0.16{col 46}{space 3}0.874{col 54}{space 4}-.0040845{col 67}{space 3} .0034828
{txt}{space 1}speakfrench {c |}{col 14}{res}{space 2} .0022187{col 26}{space 2} .0021958{col 37}{space 1}    1.01{col 46}{space 3}0.316{col 54}{space 4}-.0021585{col 67}{space 3} .0065959
{txt}{space 2}hhold_head {c |}{col 14}{res}{space 2} -.002551{col 26}{space 2} .0026663{col 37}{space 1}   -0.96{col 46}{space 3}0.342{col 54}{space 4}-.0078661{col 67}{space 3} .0027642
{txt}{space 2}hhold_size {c |}{col 14}{res}{space 2} -.000231{col 26}{space 2} .0002824{col 37}{space 1}   -0.82{col 46}{space 3}0.416{col 54}{space 4}-.0007939{col 67}{space 3}  .000332
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-.0005212{col 26}{space 2} .0051217{col 37}{space 1}   -0.10{col 46}{space 3}0.919{col 54}{space 4}-.0107312{col 67}{space 3} .0096887
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{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}    com_hist{col 14}{c |}{space 1}      134{col 27}{space 1}        0{col 39}{result}{space 1}      134{col 53}{text} {col 54}{c |}
{res}{col 1}{text}{hline 13}{c BT}{hline 12}{hline 12}{hline 14}{hline 1}{c BRC}
{res}{txt}
{com}. estimates store always_ctr
{txt}
{com}. sum index if e(sample)

{txt}    Variable {c |}        Obs        Mean    Std. dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}index {c |}{res}      5,287    .0015131    .0388734          0          1
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}.00151315
{txt}
{com}. margins, dydx(age) atmeans //  
{res}
{txt}Conditional marginal effects{col 58}{lalign 13:Number of obs}{col 71} = {res}{ralign 5:5,287}
{txt}{col 1}Model VCE: {res:Robust}

{txt}{p2colset 1 13 13 2}{...}
{p2col:Expression:}{res:Linear prediction, predict()}{p_end}
{p2col:dy/dx wrt:}{res:age}{p_end}
{p2colreset}{...}
{lalign 4:At: }{space 0}{lalign 11:age} = {res:{ralign 8:45.96614}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:age2} = {res:{ralign 8:2630.268}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:frenchnat} = {res:{ralign 8:.9899754}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:germannat} = {res:{ralign 8:.0035937}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:protestant} = {res:{ralign 8:.578589}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:other} = {res:{ralign 8:.0073766}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:speakfrench} = {res:{ralign 8:.1706071}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:hhold_head} = {res:{ralign 8:.6413845}} {txt:(mean)}
{lalign 4:}{space 0}{lalign 11:hhold_size} = {res:{ralign 8:4.4216}} {txt:(mean)}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26} Delta-method
{col 14}{c |}      dy/dx{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}age {c |}{col 14}{res}{space 2} .0001121{col 26}{space 2} .0001717{col 37}{space 1}    0.65{col 46}{space 3}0.514{col 54}{space 4}-.0002244{col 67}{space 3} .0004486
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{txt}
{com}. test age=age2=frenchnat=germannat=other=protestant=speakfrench=hhold_head=hhold_size=0

{p 0 7}{space 1}{text:( 1)}{space 1} {res}age - age2 = 0{p_end}
{p 0 7}{space 1}{text:( 2)}{space 1} age - frenchnat = 0{p_end}
{p 0 7}{space 1}{text:( 3)}{space 1} age - germannat = 0{p_end}
{p 0 7}{space 1}{text:( 4)}{space 1} age - other = 0{p_end}
{p 0 7}{space 1}{text:( 5)}{space 1} age - protestant = 0{p_end}
{p 0 7}{space 1}{text:( 6)}{space 1} age - speakfrench = 0{p_end}
{p 0 7}{space 1}{text:( 7)}{space 1} age - hhold_head = 0{p_end}
{p 0 7}{space 1}{text:( 8)}{space 1} age - hhold_size = 0{p_end}
{p 0 7}{space 1}{text:( 9)}{space 1} age = 0{p_end}

{txt}       F(  9,    72) ={res}    1.00
{txt}{col 13}Prob > F ={res}    0.4452
{txt}
{com}. scalar F_stat=r(p)
{txt}
{com}. estadd scalar F_stat

{txt}added scalar:
             e(F_stat) =  {res}.44517674
{txt}
{com}. 
. // Within birthyr eligible always-takers
. 
. reghdfe index $controls if eligible==0, cluster(birthyr) a(com_hist birthyr) // keepsing
{res}{txt}(dropped 2 {browse "http://scorreia.com/research/singletons.pdf":singleton observations})
{res}{txt}({browse "http://scorreia.com/research/hdfe.pdf":MWFE estimator} converged in 6 iterations)
{res}
{txt}HDFE Linear regression{col 51}Number of obs{col 67}= {res}     5,285
{txt}Absorbing 2 HDFE groups{col 51}F({res}   7{txt},{res}     70{txt}){col 67}= {res}      0.87
{txt}Statistics robust to heteroskedasticity{col 51}Prob > F{col 67}= {res}    0.5312
{txt}{col 51}R-squared{col 67}= {res}    0.0618
{txt}{col 51}Adj R-squared{col 67}= {res}    0.0228
{txt}{col 51}Within R-sq.{col 67}= {res}    0.0272
{txt}{col 1}Number of clusters ({res}birthyr{txt}) {col 30}= {res}        71{txt}{col 51}Root MSE{col 67}= {res}    0.0384

{txt}{ralign 78:(Std. err. adjusted for {res:71} clusters in birthyr)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}       index{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 3}frenchnat {c |}{col 14}{res}{space 2} .0040363{col 26}{space 2} .0029149{col 37}{space 1}    1.38{col 46}{space 3}0.171{col 54}{space 4}-.0017773{col 67}{space 3} .0098499
{txt}{space 3}germannat {c |}{col 14}{res}{space 2} .1116838{col 26}{space 2}  .102059{col 37}{space 1}    1.09{col 46}{space 3}0.278{col 54}{space 4}-.0918665{col 67}{space 3}  .315234
{txt}{space 2}protestant {c |}{col 14}{res}{space 2}  .001411{col 26}{space 2} .0017142{col 37}{space 1}    0.82{col 46}{space 3}0.413{col 54}{space 4} -.002008{col 67}{space 3} .0048299
{txt}{space 7}other {c |}{col 14}{res}{space 2} .0004245{col 26}{space 2} .0017738{col 37}{space 1}    0.24{col 46}{space 3}0.812{col 54}{space 4}-.0031133{col 67}{space 3} .0039623
{txt}{space 1}speakfrench {c |}{col 14}{res}{space 2} .0006489{col 26}{space 2} .0020579{col 37}{space 1}    0.32{col 46}{space 3}0.753{col 54}{space 4}-.0034554{col 67}{space 3} .0047532
{txt}{space 2}hhold_head {c |}{col 14}{res}{space 2}-.0018277{col 26}{space 2} .0024302{col 37}{space 1}   -0.75{col 46}{space 3}0.455{col 54}{space 4}-.0066746{col 67}{space 3} .0030192
{txt}{space 2}hhold_size {c |}{col 14}{res}{space 2}-.0000578{col 26}{space 2} .0003297{col 37}{space 1}   -0.18{col 46}{space 3}0.861{col 54}{space 4}-.0007153{col 67}{space 3} .0005997
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-.0023854{col 26}{space 2} .0037759{col 37}{space 1}   -0.63{col 46}{space 3}0.530{col 54}{space 4}-.0099162{col 67}{space 3} .0051455
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{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}    com_hist{col 14}{c |}{space 1}      134{col 27}{space 1}        0{col 39}{result}{space 1}      134{col 53}{text} {col 54}{c |}
{res}{col 1}{text}     birthyr{col 14}{c |}{space 1}       71{col 27}{space 1}       71{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}
* = FE nested within cluster; treated as redundant for DoF computation
{res}{txt}
{com}. estimates store always_yob_ctr
{txt}
{com}. qui sum index if e(sample)
{txt}
{com}. scalar mean_y=r(mean)
{txt}
{com}. estadd scalar mean_y

{txt}added scalar:
             e(mean_y) =  {res}.00151372
{txt}
{com}. test frenchnat=germannat=other=other=speakfrench=hhold_head=hhold_size=0

{p 0 7}{space 1}{text:( 1)}{space 1} {res}frenchnat - germannat = 0{p_end}
{p 0 7}{space 1}{text:( 2)}{space 1} frenchnat - other = 0{p_end}
{p 0 7}{space 1}{text:( 3)}{space 1} frenchnat - other = 0{p_end}
{p 0 7}{space 1}{text:( 4)}{space 1} frenchnat - speakfrench = 0{p_end}
{p 0 7}{space 1}{text:( 5)}{space 1} frenchnat - hhold_head = 0{p_end}
{p 0 7}{space 1}{text:( 6)}{space 1} frenchnat - hhold_size = 0{p_end}
{p 0 7}{space 1}{text:( 7)}{space 1} frenchnat = 0{p_end}
{txt}       Constraint 3 dropped

       F(  6,    70) ={res}    0.98
{txt}{col 13}Prob > F ={res}    0.4465
{txt}
{com}. scalar F_stat=r(p)
{txt}
{com}. estadd scalar F_stat

{txt}added scalar:
             e(F_stat) =  {res}.44653555
{txt}
{com}. 
. ******************************
. ***** 2. TABLE FOR PAPER *****
. ******************************
. 
. // Fixed effects
. 
. estfe . *, labels(com_hist "Municipality (1936) FE" ct_1936_enc "\textit{c -(}Canton{c )-} FE" birthyr "Year-of-birth FE")
{res}{txt}
{com}. return list

{txt}macros:
        r(indicate_fe) : "{res} "Municipality (1936) FE=0.com_hist"  "Year-of-birth FE=0.birthyr"  "_cons=0._cons"{txt}"

{com}. 
. // Labels
. 
. la var eligible         "Born 1908-1927"
{txt}
{com}. la var male             "[0.3em] Male"
{txt}
{com}. la var hhold_head       "[0.3em] Household head"
{txt}
{com}. la var speakfrench      "[0.3em] French-speaking"
{txt}
{com}. la var hhold_size       "[0.3em] Household size"
{txt}
{com}. la var age                      "[0.3em] Age"
{txt}
{com}. la var age2                     "[0.3em] Age $ \times $ Age"
{txt}
{com}. la var germannat        "[0.3em] German national"
{txt}
{com}. la var frenchnat        "[0.3em] French national"
{txt}
{com}. la var protestant       "[0.3em] Protestant"
{txt}
{com}. la var other            "[0.3em] Other religions"
{txt}
{com}. 
. // Table
. 
. estout OLS OLS_FE com_ctr eligible_ctr yob_ctr always_ctr always_yob_ctr ///
>         using "$results\\$program.tex", replace style(tex) type label ///
>         cells(b(star fmt(%9.3f)) se(par fmt(%9.3f)))  ///
>         collabels(none) ///
>         mlabels(, span none prefix(\multicolumn{c -(}@span{c )-}{c -(}c{c )-}{c -(}) suffix({c )-})) ///
>         msign(--) starlevels(* .1 ** .05 *** .01) ///
>         title(Estimation strategy: Conscription determinants) ///
>         drop(_cons 0._cons 0.com_hist 0.birthyr) ///
>         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 )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(6){c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(7){c )-} \\" /// 
>         "\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Dep. Variable{c )-}{c )-} &\multicolumn{c -(}7{c )-}{c -(}c{c )-}{c -(}Individual is conscripted{c )-}\\" /// 
>         "\cmidrule(lr){c -(}2-8{c )-}" ///
>         "\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Sample{c )-}{c )-}&\multicolumn{c -(}3{c )-}{c -(}c{c )-}{c -(}Full sample (male){c )-}     &\multicolumn{c -(}2{c )-}{c -(}c{c )-}{c -(}Eligible (male){c )-}    &\multicolumn{c -(}2{c )-}{c -(}c{c )-}{c -(}Ineligible (male){c )-}\\" /// 
>         "\cmidrule(lr){c -(}2-4{c )-} \cmidrule(lr){c -(}5-6{c )-} \cmidrule(lr){c -(}7-8{c )-}" ///
>         "\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Specification{c )-}{c )-}&\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}No{c )-}                      &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}             &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}             &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}             &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}     &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}     &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}     \\" /// 
>         "       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}controls{c )-}   &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}municipality{c )-}       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}year-of-birth{c )-}      &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}municipality{c )-}       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}year-of-birth{c )-}      &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}municipality{c )-}       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}year-of-birth{c )-}\\") /// 
>         posthead("\hline" ///
>         "\\") ///
>         prefoot("") ///
>         postfoot("\\")
{res}
\begin{c -(}table{c )-}[t!]
\begin{c -(}center{c )-}
\begin{c -(}threeparttable{c )-}
\begin{c -(}scriptsize{c )-}
\caption{c -(}Estimation strategy: Conscription determinants\label{c -(}EEREV-D-21-00004_Table1{c )-}{c )-}
\begin{c -(}tabular*{c )-}{c -(}\columnwidth{c )-}{c -(}@{c -(}\extracolsep{c -(}\stretch{c -(}1{c )-}{c )-}{c )-}l*{c -(}7{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 )-} &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(6){c )-}  &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}(7){c )-} \\
\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Dep. Variable{c )-}{c )-} &\multicolumn{c -(}7{c )-}{c -(}c{c )-}{c -(}Individual is conscripted{c )-}\\
\cmidrule(lr){c -(}2-8{c )-}
\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Sample{c )-}{c )-}&\multicolumn{c -(}3{c )-}{c -(}c{c )-}{c -(}Full sample (male){c )-}     &\multicolumn{c -(}2{c )-}{c -(}c{c )-}{c -(}Eligible (male){c )-}    &\multicolumn{c -(}2{c )-}{c -(}c{c )-}{c -(}Ineligible (male){c )-}\\
\cmidrule(lr){c -(}2-4{c )-} \cmidrule(lr){c -(}5-6{c )-} \cmidrule(lr){c -(}7-8{c )-}
\multicolumn{c -(}1{c )-}{c -(}l{c )-}{c -(}\textit{c -(}Specification{c )-}{c )-}&\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}No{c )-}                      &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}             &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}             &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}             &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}     &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}     &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}Within{c )-}     \\
       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}controls{c )-}   &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}municipality{c )-}       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}year-of-birth{c )-}      &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}municipality{c )-}       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}year-of-birth{c )-}      &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}municipality{c )-}       &\multicolumn{c -(}1{c )-}{c -(}c{c )-}{c -(}year-of-birth{c )-}\\
\hline
\\
Born 1908-1927      &       0.148***&       0.129***&               &               &               &               &               \\
                    &     (0.020)   &     (0.021)   &               &               &               &               &               \\
[0.3em] Age         &               &     --0.002   &               &       0.035   &               &       0.000   &               \\
                    &               &     (0.001)   &               &     (0.032)   &               &     (0.000)   &               \\
[0.3em] Age $ \times $ Age&               &       0.000   &               &     --0.001   &               &     --0.000   &               \\
                    &               &     (0.000)   &               &     (0.001)   &               &     (0.000)   &               \\
[0.3em] French national&               &       0.020   &       0.018   &       0.048   &       0.058   &       0.002   &       0.004   \\
                    &               &     (0.015)   &     (0.015)   &     (0.052)   &     (0.051)   &     (0.002)   &     (0.003)   \\
[0.3em] German national&               &       0.106   &       0.109   &       0.139   &       0.131   &       0.108   &       0.112   \\
                    &               &     (0.082)   &     (0.084)   &     (0.138)   &     (0.148)   &     (0.100)   &     (0.102)   \\
[0.3em] Protestant  &               &       0.004   &       0.004   &       0.002   &       0.006   &       0.001   &       0.001   \\
                    &               &     (0.009)   &     (0.009)   &     (0.029)   &     (0.028)   &     (0.002)   &     (0.002)   \\
[0.3em] Other religions&               &     --0.034** &     --0.027*  &     --0.079** &     --0.074** &     --0.000   &       0.000   \\
                    &               &     (0.014)   &     (0.015)   &     (0.030)   &     (0.031)   &     (0.002)   &     (0.002)   \\
[0.3em] French-speaking&               &       0.024***&       0.011** &       0.039** &       0.026   &       0.002   &       0.001   \\
                    &               &     (0.006)   &     (0.005)   &     (0.018)   &     (0.018)   &     (0.002)   &     (0.002)   \\
[0.3em] Household head&               &       0.007   &       0.002   &       0.041   &       0.025   &     --0.003   &     --0.002   \\
                    &               &     (0.011)   &     (0.006)   &     (0.030)   &     (0.031)   &     (0.003)   &     (0.002)   \\
[0.3em] Household size&               &       0.002   &     --0.000   &       0.000   &       0.000   &     --0.000   &     --0.000   \\
                    &               &     (0.001)   &     (0.001)   &     (0.004)   &     (0.004)   &     (0.000)   &     (0.000)   \\

\\
{txt}(output written to {browse  `"D:\Dropbox\SVLACHOS\Research\Malgre_nous\Replication_files\Results\EEREV-D-21-00004_Table1.tex"'})

{com}. 
. estout OLS OLS_FE com_ctr eligible_ctr yob_ctr always_ctr always_yob_ctr ///
>         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 r2 F_stat, lab("Mean dep. variable" "Observations" "Clusters" "$ R $-squared" "Controls $ F $-test ($ p $-value)") fmt(3 %9.0fc 0 2 2)) /// 
>         collabels(none) ///
>         mlabels(, span none prefix(\multicolumn{c -(}@span{c )-}{c -(}c{c )-}{c -(}) suffix({c )-})) ///
>         msign(--) starlevels(* .1 ** .05 *** .01) ///
>         drop(*) ///
>         indicate("Municipality (1936) FE = 0.com_hist" ///      
>         "Year-of-birth FE = 0.birthyr" ///
>         , labels(Yes )) ///
>         prehead("") ///
>         posthead("") ///
>         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 )-}:" /// 
>         "OLS estimates of the likelihood of appearing in the \textit{c -(}Index of French Nationals Compelled into German Armed Forces{c )-} \citep{c -(}ministere1945basrhin, ministere1946moselle{c )-}. " ///
>         "The unit of observation is a male individual. " ///
>         "%Columns (1)-(3): all male individuals; " ///
>         "%Columns (4) and (5): eligible males (born in 1908-1927); " ///
>         "%Columns (6) and (7): ineligible males (not born in 1908-1927). " ///
>         "%Municipality fixed effects are included in Columns (2) to (7). " ///
>         "%Year-of-birth fixed effects are included in Columns (3), (5) and (7). " ///
>         "Standard errors clustered at the year-of-birth-level in parentheses. " /// 
>         "* 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}


Municipality (1936) FE &               &         Yes   &         Yes   &         Yes   &         Yes   &         Yes   &         Yes   \\
Year-of-birth FE    &               &               &         Yes   &               &         Yes   &               &         Yes   \\

\hline
Mean dep. variable  &       0.047   &       0.047   &       0.047   &       0.150   &       0.150   &       0.002   &       0.002   \\
Observations        &       7,590   &       7,590   &       7,588   &       2,303   &       2,303   &       5,287   &       5,285   \\
Clusters            &          93   &          93   &          91   &          20   &          20   &          73   &          71   \\
$ R $-squared       &        0.10   &        0.13   &        0.17   &        0.10   &        0.13   &        0.05   &        0.06   \\
Controls $ F $-test ($ p $-value)&               &        0.03   &        0.04   &        0.00   &        0.11   &        0.45   &        0.45   \\
\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 )-}:
OLS estimates of the likelihood of appearing in the \textit{c -(}Index of French Nationals Compelled into German Armed Forces{c )-} \citep{c -(}ministere1945basrhin, ministere1946moselle{c )-}. 
The unit of observation is a male individual. 
%Columns (1)-(3): all male individuals; 
%Columns (4) and (5): eligible males (born in 1908-1927); 
%Columns (6) and (7): ineligible males (not born in 1908-1927). 
%Municipality fixed effects are included in Columns (2) to (7). 
%Year-of-birth fixed effects are included in Columns (3), (5) and (7). 
Standard errors clustered at the year-of-birth-level in parentheses. 
* 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_Table1.tex"'})

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