AIC, BIC, codes for stata

sort ccode year 

by ccode: gen lag4annualcoups = annualcoups[_n-4]

by ccode: gen lagannualcoups = annualcoups[_n-3]

by ccode: gen lag2annualcoups = annualcoups[_n-2]

by ccode: gen lag1annualcoups = annualcoups[_n-1]

 xtreg mobile l.mobile c.v2x_polyarchy##c.gdp_growth_percent lngdppc lagannualcoups,fe

*estat ic

est sto coup1

xtreg mobile l.mobile c.v2x_polyarchy##c.gdp_growth_percent lngdppc lag2annualcoups,fe

*estat ic

est sto coup2

xtreg mobile l.mobile c.v2x_polyarchy##c.gdp_growth_percent lngdppc lag1annualcoups,fe

*estat ic

est sto coup3

xtreg mobile l.mobile c.v2x_polyarchy##c.gdp_growth_percent lngdppc lag4annualcoups,fe

*estat ic

est sto coup4

xtreg mobile l.mobile c.v2x_polyarchy##c.gdp_growth_percent lngdppc,fe

*estat ic

est sto nocoup

est table coup1 coup2 coup3 coup4 nocoup, stats(r2_p aic bic)

esttab m1 m2 m3 using logit.rtf,  se aic obslast scalar(F) bic r2