use f10_4_1, clear

save f10_4_2, replace

save f10_4_2a, replace

save f10_4_3, replace

save f10_4_3a, replace

    di in smcl "{p}{err}At least one of the packages was not able to be " ///
               "installed.  Try typing {cmd:spinst_itsp} again.  If that " ///
               "fails, look at the output of the error message above to " ///
               "diagnose the problem.  Perhaps you have an out-of-date " ///
               "version of the package already installed.  In that case, " ///
               "type {cmd:adoupdate} to update it."

use fem2

use fem2, clear

do fem_setup.do

tempname fh

use fictprice, clear

tsline _b_spx, yline(1) saving(fig7_1, replace) scheme(s2mono)

qnorm p50, saving(fig7_3, replace) scheme(s2mono)

graph combine fig9_2a fig9_2b, col(2) ti("Firm statistics")

graph combine fig9_2aa fig9_2bb, col(2) ti("Firm statistics") 

graph dot beta_k, over(compname) scheme(s2mono) nodraw name(fig9_2bb,replace) yline(0)

graph combine fig9_2a fig9_2d, col(2) ti("Firm statistics")

graph use fig10_4_2

	estimates notes: from file `c(filename)' saved `c(filedate)'

	estimates notes: from file `c(filename)' saved `c(filedate)'

	estimates notes: from file `c(filename)' saved `c(filedate)'

graph combine fig9_2aa fig9_2bb, col(2) ti("Firm statistics") 

maxindcorr index1-index9, ret(ret) firm(permno) gen(set1)

foreach f of local firms {

use firmstats, clear

use firmstats2, clear

local flav = cond(`c(MP)', "MP", cond(`c(SE)', "SE", "IC"))

* label var reg "Census region"
label data "1980 Census data for NE and NC states"

note medagel: median age for large states only

	esttab _all using lwage`y'.tex, replace ti("Wage equations for `y'") ///
	       nomtitles nodepvars not se noobs ar2 booktabs ///
	       scalar(eta_educ etase exper se) ///
	       addnotes("eta_educ: semi-elasticity of lwage with respect to educ" ///
	       "etase: standard error of the semi-elasticity" ///
	       "exper: effect of exper on lwage at mean exper") ///
	       substitute("_cons" "Constant" "eta_educ" "$\eta_{educ}$" ///
	       "etase" "$\eta_{se}$") 

estimates for model1 model2 model3 model4: test temp = -1.6

varextrema price if foreign

global subsample foreign

outtable using MCAS, mat(MCAS) caption("MCAS percentile ranges") format(%9.2f) nobox replace

esttab three four, margin

foreach c in US UK DE FR  {

local country US UK DE FR

local country US UK DE FR

  cap noi spinst_itsp_wrk freduse http://fmwww.bc.edu/repec/bocode/f

histogram medagebrack, discrete frequency ///
lcolor(black) fcolor(gs15) addlabels ///
addlabopts(mlabposition(6)) xtitle(Upper limits of median age) ///
title(Northeast and North Central States: Median Age) 

* CEO
makematrix t2aa, from(r(Stat1) r(Stat2)): tabstat ceo, by(insider) stat(mean) save

    capture noisily net install `package', from(`from')

    qui makematrix t2b, from(r(Stat1) r(Stat2)): tabstat `v', by(insider) stat(mean) save

  args package from

    di as smcl "{cmd}. net install `package', from(`from')

aggreg tenn, per(3) func(first) // extract first month's value

aggreg2 tenn, per(3) func(sum) trans(sqrt) // calculate quarterly sum of sqrt(tenn)

// calculate price series as settle(nnear) for rollover date et seq.
bysort qm: replace futprice = settle if contract==nextnear & qmdy >= rollover

// drop obs no longer needed, flagged by missing settle
drop if futprice == .

bysort qm: generate double futret = log(settle) - log(sprev) if qmdy == rollover

list qmdy contract futprice futret if !missing(futprice) in 1/62,  noobs sepby(qm)