tscollap (last) aireoq=air (sum) airsum=air (mean) airavg=air, to(q) gen(yq)

generate qtr = qofd(dofm(ym))

generate qtr = qofd(dofm(ym))

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

format qtr %tq

label var _b_spx "quantile reg IBM beta"

tsmktim quarter, start(1978q1) // create quarterly calendar var

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

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

label var math4sat "Quartile of percent satisfactory, 4th grade"

	qui regress `c' `allbut' date

// reestimate using IV
qui ivregress 2sls so2 (pop = temp precip) manuf wind

qui probit hiso2 pop manuf hiprecip

    qui ranksum `v', by(insider)

	qui replace gdptot = gdptot + `c'

     qui replace nschool = nschool + inrange(`v', 6, 18)

     qui replace anyschool = max(anyschool, inrange(`v', 6, 18))

qui generate dscorr = 100 * (p - p[_n-1])

qui generate dscorrd = dscorr / sqrt(deld)

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

qui gen ndays = .

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

qui generate highcorr = .

	qui replace model = "`model'" in `i'

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

qui summarize so2, detail

qui levelsof hiord, local(indices)

	qui correlate index1-index`nind' ret if permno == `f'

	quietly {

	quietly {

// run regression
quietly regress price nox dist rooms stratio

quietly regress so2 temp pop wind

quietly regress so2 temp manuf pop wind precip days

eststo: quietly regress so2 temp pop wind

eststo: quietly regress so2 temp manuf pop wind precip days

quietly ivreg2 lw s expr tenure rns smsa  (iq=med kww age mrt), gmm2s 

quietly ivreg2 lw s expr tenure rns smsa  (iq=med kww age mrt), cue 

quietly reshape wide

// compute elasticities at those points
   quietly mfx compute, eyex at(mean stratio = `pc`p'')  

   quietly replace `eyex' = eta[1, 1] in `i'    // and save in tempvar eyex 

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

quietly replace mvalue = . in 55

quietly replace kstock = . in 94

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

quietly generate dcal = D.caldate if seqdate < .

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

quietly pctrange tot_day

quietly mygmm2s lw, endog(iq) inexog(s expr tenure rns smsa) ///
                    exexog(med kww age mrt) robust

// summarize, detail computes percentiles of stratio
quietly summarize stratio if e(sample), detail 

quietly log close