avplot lhours, msize(small) nodraw scheme(s2mono) name(avplot`y', replace) 

scatter expp1992 expp1998 sal, m(oh x) scheme(s2mono) ///
  title("School District Expenditures vs. Average Salaries")  ///
  xtitle("Average salaries") ytitle("Expenditures") ///
  legend(position(5) ring(0) label(1 "1992") label(2 "1998"))

histogram p50, normal saving(fig7_2, replace) scheme(s2mono)

	scatter lwage0 lwage1 lhours if year==`y', msize(small) ///
	        scheme(s2mono) nodraw ///
	        name(scat`y',replace) ytitle(lwage) legend(pos(5) ring(0) col(1))

tsline dDTB3, yline(0) scheme(s2mono)

foreach v of varlist lexp gnppc safewater {

label val highsal sal

scatter expp1992 expp1998 sal, m(oh x) scheme(s2mono) ///
  title("School District Expenditures vs. Average Salaries")  ///
  xtitle("Average salaries") ytitle("Expenditures") ///
  legend(position(5) ring(0) label(1 "1992") label(2 "1998"))

scatter expp1992 expp1998 sal, m(oh x) scheme(s2mono) ///
  title("School District Expenditures vs. Average Salaries")  ///
  xtitle("Average salaries") ytitle("Expenditures") ///
  legend(position(5) ring(0) label(1 "1992") label(2 "1998"))

label var p50 "Sample median"

label def satis 1 "Poor" 2 "Good" 3 "Better" 4 "Best"

label val math7sat satis

label var math7sat "Quartile of percent satisfactory, 7th grade"

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

savedresults save ivrobust e()

savedresults save ivreg2cue e()

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

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

rolling _b _se, window(90) saving(capm, replace) nodots: qreg ibm spx

qui rolling r(mean), window(3) saving(ldinvest, replace): summarize LD.linvest

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

	graph combine scat`y' avplot`y', nodraw col(2) saving(lwage`y', replace) ///
	      ti("Log wage vs. log hours for `y'") 

savedresults compare ivhomo e(), include(macros: depvar scalar: N j matrix: b V) ///
             tol(1e-7) verbose

savedresults compare ivreg2cue e(), include(macros: depvar scalar: N j matrix: b V) ///
             tol(1e-4) verbose

	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}$") 

// graph the scatter of elasticities vs. percentiles of stratio
// as well as the predictions with rline 
// and the 95% confidence bands with connected
twoway (scatter `eyex' `x', ms(Oh) yscale(range(-0.5 -2.0)) ylabel(,nogrid)) ///
   (rline `seyex1' `seyex2' `x') ///
   (connected `y' `x', yaxis(2) yscale(axis(2) range(18000 35000))), ///
   ytitle(Elasticity of price v. student/teacher ratio)

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

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

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

tsline ssq, ti("Monthly volatility from Merton method") scheme(s2mono)

graph bar p50, over(compname) scheme(s2mono) nodraw horiz name(fig9_2aa,replace)

scatter expp1992 expp1998 sal, m(oh x) scheme(s2mono) ///
  title("School District Expenditures vs. Average Salaries")  ///
  xtitle("Average salaries") ytitle("Expenditures") ///
  legend(position(5) ring(0) label(1 "1992") label(2 "1998"))

esttab one two, mtitles("OLS" "IV") stat(nu nuse, labels("Elas.(pop)" "SE"))

g lo = sum - 1.96 * se

rolling sum=r(sum) se=r(se), window(30): ///
myregress wpi L(1/4).wpi t, lagvar(wpi) nlags(4)

	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}$") 

		eststo, addscalars(eta_educ eta[1, 1] etase eta[2, 1] exper r(estimate) se r(se))

set seed 2007062926

set seed 2007062926

set seed 20071110

scalar semean = r(sd) / sqrt(r(N))

	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}$") 

list, sep(0)

list state pop smallpop largepop, sep(0)

list state medage pop if inrange(pop, 5000, 9999), sep(0)

list if beginunemp, sep(0)

list ym mnr eoq air in 1/16, sep(4) noobs

list yq airsum in 1/16, sep(4) noobs

list in 1/16, sep(4) noobs

list make price weight length in 1/10, sep(0) subvarname

list t tenn quarter tennA3 year tennA12 in 1/12, sep(3)

list t tenn quarter tennA3 tennS3 tennF3 tennL3 in 1/12, sep(3)

list set1permno set1mu set1n set1max in 1/50, noobs sep(0)

list pop mo1 so2 mv1, sep(0)

list month sumchange if sumchange < ., sep(0) noobs

list, sepby(state)

list region state pop totpop, sepby(region)

list region popsize state pop meanpop2, sepby(region popsize)

list, sepby(employer) noobs

// rename family fam
// drop nadult-oldestsibling
list fam person female age, sepby(fam) noobs

list, noobs sepby(state)

list in 1/51, sepby(qmdy)

quietly generate seqdate = cond(notmiss,sum(notmiss),.)

quietly generate month = mofd(caldate) if seqdate < .

tsset seqdate