regress weight age age2

	qui regress `c' `allbut' date

// run regression
quietly regress price nox dist rooms stratio

regress deaths duration troop

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

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

qui probit hiso2 pop manuf hiprecip

xi: regress bp i.agegrp*i.sex

ttest bp, by(sex)

ttest bp, by(sex)

	quietly regress invest L(1/2).kstock L.kapavg if company == `c'

		regress invest L.kstock L.kapavg if company == `c'

		regress lwage educ lhours black hisp if year == `y'

	regress lwage educ exper lhours if year == `y'

    qui ranksum `v', by(insider)

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

ml display

ml model d0 mynormal_lf_d0_c1 (a:) (mu: price = weight turn) /lnsigma, ///
maximize nolog diparm(lnsigma, exp label("sigma")) 

ml model lf meanvar (mu1: price = ) (mu2: price = ) /sigma1 /sigma2

ml model lf meanvar (mu1: price = ) (mu2: price = ) /sigma1 /sigma2, ///
   constraints(1)

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