matrix t2a = t2a, temp

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

mat t2aa = t2a \ (t2aa, temp)

matrix rownames t2aa = Observations CEO

file write `hh' "Observations (\$N\$) & " (t2aa[1, 1]) " & " (t2aa[1, 2]) " & " (t2aa[1, 3]) " \\" _newline

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

    matrix t2bb = (nullmat(t2bb) \ `v')

matrix list t2bb

file write `hh' "Departed (0/1) & "  %7.3f (t2bb[`i', 1]) " & "  %7.3f (t2bb[`i', 2]) " & " (t2bb[`i', 3])  " & " %7.3f (t2bb[`i', 4])  " & " %7.3f (t2bb[`i', 5]) " \\" _newline

matrix colnames t2bb = Outsider Insider N t_pval MW_pval

estimates table model4 model1 model3 model2, stat(r2_a rmse ll) /// 
b(%7.3g) star label title("Models of sulphur dioxide concentration") 

outtable using ch9.02t, mat(table1) replace ///
format(%9.3f) center caption("Grunfeld company statistics")

		matrix table1[`i', 1] = r(p25)

		matrix table1[`i', 3] = r(p75)

		matrix table1[`i', 5] = r(rho)

		matrix table1[`i', 7] = _se["L.kstock"]

matrix colnames table1 = p25 p50 p75 r_invest r_kap beta_k se_k

local rows = rowsof(table1)

		file write `fh' (table1[`i', `j']) " "

local cols = colsof(table1)

local rown : rownames table1

file open `hh' using table2.tex, write replace

merge tabrv using f10_4_2a, update unique

generate tabrv = upper(substr(town, 1, 3))

scatter so2 precip, msize(tiny) mlabel(tabrv) mlabsize(vsmall) ///
by(hightemp) saving(fig10_4_2, replace) scheme(s2mono)

scatter so2 precip, msize(tiny) mlabel(tabrv) mlabsize(vsmall) ///
by(hightemp) saving(fig10_4_2, replace) scheme(s2mono)

sort tabrv

sort tabrv

sort tabrv

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

egen indivcount = tag(state cand)

label def ind 1 Kappa 2 Lambda 3 Nu 4 Xi 5 Tau 6 Upsilon 7 Phi 8 Chi 9 Psi

format date %td

format end %tdMonYY

label var avgsal "Average teacher salary"

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

list t tenn quarter tennAlog3 tennSsqrt3 in 1/12, sep(3)

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

list t tenn quarter tennAlog3 tennSsqrt3 in 1/12, sep(3)

local inv6 "Board Tenure (Years)"

file write `hh' "\hline & Mean & Mean & \$N\$ & \$t\$ test & Mann--Whitney \\" _newline

do test1

do test3

type testall.do

capture erase bizday1.log.tex

capture erase bizday3.log.tex

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

    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."

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

display "There are `wds' anxiety levels:"

parrot "This is Stata..."

label var invrate "Three-year average investment rate"

estimates store three

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

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

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

mat list result, ti("IQ by anxiety score for <35, >=35 years age") 

local timestamp "Matrix created `c(current_date)' `c(current_time)'"

keep if tin(1aug2005,)

scatter so2 precip, msize(tiny) mlabel(tabrv) mlabsize(vsmall) ///
by(hightemp) saving(fig10_4_2, replace) scheme(s2mono)

scatter so2 precip, msize(tiny) mlabel(tabrv) mlabsize(vsmall) ///
by(hightemp) saving(fig10_4_2, replace) scheme(s2mono)

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

twoway rarea lo hi end, color(gs12) title("Sum of moving lag coefficients, approx. 95% CI") ///
||  tsline sum, legend(off) scheme(s2mono)

label values hightemp tlab

// compute appropriate t-statistic  for 95% confidence interval
scalar tmfx = invttail(e(df_r), 0.975)

   quietly replace `seyex2' = `eyex' - tmfx*se[1, 1] in `i'