qui replace ndays = r(N) in `n'

statsmat iq if age >= 35, by(anxiety) stat(N mean sd) matrix(iqanx2)

* N
makematrix t2a, from(r(Stat1) r(Stat2)): tabstat insider, by(insider) stat(N) save

tabstat medage, by(region) statistics(N mean sd min max) 

assert r(unique_value) == r(N)

matrix temp = r(N)

    matrix `v'1 = r(N)

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

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

local ncomp : word count `colist'

local ncty: word count country

	qui replace ndays = r(N) in `n'

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

// qmdy is the first date when the near contract is no longer quoted
// that minus 3 trading days is the target rollover date
bysort cmdy (qmdy): generate rolldate = qmdy[_n-3] if near[_n] != near[_n-1] 

bysort cmdy (qmdy): generate nnear = contract if near[_n] != near[_n-1]

use NEstates, clear

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

label values netmarr2x netmarr2xc

label define netmarr2xc 1 "marr > 2 divr" 2 "marr <= 2 divr"

generate byte newengland = inlist(state, "Connecticut", "Maine", ///
"Massachusetts", "New Hampshire", "Rhode Island", "Vermont")

recode cc `rc', gen(newcc)

tabulate newcc

list state medage pop if newengland, sep(0)

bysort qm: egen nextnear = mean(nnear)

// calculate return series for the rollover date
bysort qm: generate settleprev = settle if contract == nextnear & qmdy == rollover1

local nind : word count `indices'

type nlces.ado

// fixup for first obs
replace nnear = . in 1

type nneighbor.ado

matrix accum C = weight age iq, dev nocons

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

bootstrap diff=(r(mu_1)-r(mu_2)), nodots: ttest bp, by(sex)

simulate p25=r(p25) p50=r(p50) p75=r(p75), reps(1000) nodots ///
         saving(pctiles, replace): simpctile

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

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

// 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)

mat list grades, noheader

  cap noi spinst_itsp_wrk egenmore http://fmwww.bc.edu/repec/bocode/e

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

  cap noi spinst_itsp_wrk mat2txt http://fmwww.bc.edu/repec/bocode/m

  cap noi spinst_itsp_wrk moremata http://fmwww.bc.edu/repec/bocode/m

  cap noi spinst_itsp_wrk onespell http://fmwww.bc.edu/repec/bocode/o

  cap noi spinst_itsp_wrk outtable http://fmwww.bc.edu/repec/bocode/o

  cap noi spinst_itsp_wrk ranktest http://fmwww.bc.edu/repec/bocode/r

  cap noi spinst_itsp_wrk statsmat http://fmwww.bc.edu/repec/bocode/s

  cap noi spinst_itsp_wrk tscollap http://fmwww.bc.edu/repec/bocode/t

ml maximize, nolog

esttab using esttab_example.tex, ar2 label se nostar nodepvars brackets ///
nomtitles title("Models of sulphur dioxide concentration") booktabs ///
alignment(D{.}{.}{-1}) replace

local anxlevel None Mild Moderate Severe 

local anxlevel None Mild Moderate Severe 

local anxlevel None Mild Moderate Severe 

local anx None Mild Moderate Severe

list, sepby(employer) noobs

list, sep(0) noobs

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

list, noobs

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 qmdy contract futprice futret if !missing(futprice) in 1/62,  noobs sepby(qm)

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

pctrange regday specneed bilingua occupday tot_day tchratio, mat(MCAS) noprint

    matrix `v'3 = 1 - normprob(abs(r(z)))

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