matrix b = e(b)

local regressors: colnames e(b)

	matrix beta`i' = e(b)

nl ces @ lnoutput capital labor, parameters(b0 rho delta) ///
 initial(b0 0 rho 1 delta 0.5)

matrix b0[1,1] = ln(-1*b0[1,1])

generate byte beginspell = employer != employer[_n-1]

file write `hh' "\begin{table}[htbp]\caption{Director-level variables}\bigskip" _newline

generate byte beginunemp = missing(employer) & ((employer != employer[_n-1]) | (_n == 1))

label def tlab 0 "below mean temp" 1 "above mean temp"

cd bench/pctrange

cd bench/pctrange

label var _b_spx "quantile reg IBM beta"

label var beta_k "beta(lagged kstock)"

file write `hh' "\begin{table}[htbp]\caption{Director-level variables}\bigskip" _newline

capture erase bizday1.log.tex

capture erase bizday2.log.tex

capture erase bizday3.log.tex

capture erase bizday4.log.tex

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

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

sysuse bplong

use bpress, clear

use bpress, clear

use bpstats, clear

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

ttest bp, by(sex)

ttest bp, by(sex)

    qui ttest `v', by(insider)

egen totcand = total(last), by(state)

note popsize: variable separating states by population size

by region: egen meanpop = mean(pop)

by fam: egen nadult = total(age >= 18)

by fam: egen nboys = total(age <= 17 & !female) if age <= 17

by fam: g older = (age <= 17 & age[_n+1] <= 17) * (age < age[_n+1])

collapse (sum) airsum=air, by(qtr)

quietly by month: generate adjchange = sum(dDTB3/sqrt(dcal))

statsmat iq if age < 35, by(anxiety) stat(N mean sd) matrix(iqanx1)

separate expp, by(year) veryshortlabel

separate lwage, by(black | hisp) veryshortlabel

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

tabstat last, stat(sum) nototal by(state)

by region: summarize medage marr divr

tabstat pop medage, by(netmarr2x)

tabstat pop, stat(mean min max) by(size)

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)

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

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

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

bys cdate hiord: egen retport = mean(ret)

bysort state cand: generate byte last=(_n == _N)

bysort region popsize: egen meanpop2 = mean(pop)

bysort spellnr: egen meanwage = mean(wage)

qui bysort mon: generate ssq = sqrt(cumds2) if _n == _N

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

bysort qm: egen rollover1 = mean(roll1date)

// calculate price series as settle(near) for pre-rollover dates
bysort qm: generate futprice = settle if contract == near & qmdy < rollover

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

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

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

bysort highsal: pctrange regday specneed bilingua occupday tot_day tchratio

bysort popsize: summarize medage

bysort state cand: generate byte last=(_n == _N)

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

generate byte anyschool = nschool > 0

generate byte beginunemp = missing(employer) & ((employer != employer[_n-1]) | (_n == 1))

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

display _newline "Run `date' at `time' on Stata/`flav' version `vers', memory = `mem' bytes"