generate R5pc_1 = (p1 < 0.05) 

forvalues i = 1/`=r(N)' {

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

generate hiprecip = (precip > r(mean)) & !missing(precip)

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

generate byte highsal = avgsalry > r(mean) & !missing(avgsalry)

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

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

pctrange tot_day if percap > r(mean) & !missing(percap)

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

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

	local pc`p' = r(p`p')

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

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

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

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

		local p50 = r(p50)

		local r_invest = r(rho)

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

mat temp = r(N)

    matrix `v'2 = r(p)

scalar zmu = r(mean)

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

tsline invrate if inrange(company, 1, 4) & year >= 1937, by(company, ti(" ")) ///
   yline(0) ysc(range(-0.2 0.3)) ytick(-0.2(0.1)0.3) ylab(-0.2(0.1)0.3) ///
   scheme(s2mono)

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

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

label var invrate "Three-year average investment rate"

label variable `x' "Student/teacher ratio (percentiles `pct')" 

	local rc "`rc' (`i'=`c')"

egen nschool = rcount(child1-child12), cond(@ > 5 & @ < 19)

local RS real scalar

textbarplot compname beta_k, recast(dot) xline(0) ///
scheme(s2mono) nodraw name(fig9_2b,replace)

generate medagebrack = recode(medage, 29, 30, 31, 32, 33)

display "Regressors: `regressors'"

merge rep78 using repstats, uniqusing

bysort rep78: egen iqr = pctrange(price) if inrange(rep78, 3, 5)

tabstat iqr if inrange(rep78, 3, 5), by(rep78)

sort rep78

statsby mean=r(mean) sd=r(sd) n=r(N), by(agegrp sex) saving(bpstats, replace): summarize bp

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

outtable using ch9.03t, mat(allbeta) replace ///
format(%9.4f) center caption("Grunfeld company estimates")

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

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

save cb5along, replace

outsheet using census2c, replace

save ex4.1, replace

save ch45,replace

save mathpnl_state, replace

	save mathpnl`i'm, replace

save mathpnl_long2, replace

save check82, replace

save myvolat, replace

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

save f10_4_1, replace

save f10_4_2a, replace

save f10_4_3a, replace

save ch14.2,replace

save fem2,replace

quietly log using pctrange, replace

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

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

file open `fh' using matout.raw, write replace

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

duplicates report state

merge rep78 using repstats, uniqusing

  di as smcl "{txt}   package {res:`package'} from {res:`from'}"

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

matrix result = iqanx1, iqanx2 

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

svmat double retcorr