**** LHS wants a regression that is purely cross-sectional, he says let's regress by year 
**** first we'll need annual average beta 

use "/Users/natashasarin/Dropbox/Z__Natamir/Natasha_Research/Banks/Bailout_2016/data_20160808.dta", clear

by year permno, sort: egen mean_beta = mean(beta_log_fwd)
drop output_indicator 
by year permno, sort: gen output_indicator = _n 
keep if output_indicator == 1 
by year, sort: reg mean_beta bhck7204_inv

save "/Users/natashasarin/Dropbox/Z__Natamir/Natasha_Research/Banks/Bailout_2016/data_20160815.dta", replace  

statsby "regress mean_beta bhck7204_inv" _b _se rsq=e(r2) adjrsq=e(r2_a) fstat=e(F) resid=e(df_r), by (year) 

gen tval = b_bhck7204_inv/se_bhck7204_inv
gen pval = (2 * ttail(resid, abs(b_bhck7204_inv/se_bhck7204_inv)))  


outsheet using test, replace
type test.out 



**** below is scrap materials 

*** we want to loop through years 
/*egen order = group(year), label 
su order, meanonly 
forval y = 1/`r(max)' { 
reg mean_beta bhck7204_inv 
outreg2 bhck7204_inv using "`y'"
}*/

save "/Users/natashasarin/Dropbox/Z__Natamir/Natasha_Research/Banks/Bailout_2016/data_20160815.dta", replace 